Valerio PC

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

Universidad Nacional Mayor de San Marcos

Universidad del Perú. Decana de América


Facultad de Ciencias Físicas
Escuela Profesional de Ingeniería Mecánica de Fluidos

Estudio de hidráulica fluvial para la protección contra


inundaciones entre el tramo desembocadura –
localidad Paullo del río Cañete, provincia de Cañete –
departamento de Lima

INFORME PROFESIONAL

Para optar el Título Profesional de Ingeniero Mecánico de Fluidos

AUTOR
Carlos Ivan VALERIO PONCE

Lima, Perú

2022
Reconocimiento - No Comercial - Compartir Igual - Sin restricciones adicionales

https://creativecommons.org/licenses/by-nc-sa/4.0/
Usted puede distribuir, remezclar, retocar, y crear a partir del documento original de modo no
comercial, siempre y cuando se dé crédito al autor del documento y se licencien las nuevas
creaciones bajo las mismas condiciones. No se permite aplicar términos legales o medidas
tecnológicas que restrinjan legalmente a otros a hacer cualquier cosa que permita esta licencia.
Referencia bibliográfica

Valerio, C. (2022). Estudio de hidráulica fluvial para la protección contra


inundaciones entre el tramo desembocadura – localidad Paullo del río Cañete,
provincia de Cañete – departamento de Lima. [Informe profesional de pregrado,
Universidad Nacional Mayor de San Marcos, Facultad de Ciencias Físicas, Escuela
Profesional de Ingeniería Mecánica de Fluidos]. Repositorio institucional Cybertesis
UNMSM.
Metadatos complementarios

Datos de autor

Nombres y apellidos Carlos Ivan Valerio Ponce.


DNI
Tipo de documento de identidad

Número de documento de identidad 46446392

URL de ORCID https://orcid.org/0000-0003-1204-8477

Datos de asesor

Nombres y apellidos

Tipo de documento de identidad

Número de documento de identidad

URL de ORCID

Datos del jurado

Presidente del jurado

Nombres y apellidos Arístides Ramón Serruto Colque

Tipo de documento DNI

Número de documento de identidad 01284486

Miembro del jurado 1

Nombres y apellidos Guido Américo Rozas Olivera.

Tipo de documento DNI

Número de documento de identidad 07209368

Miembro del jurado 2

Nombres y apellidos Miguel Ernesto Arámbulo Manrique.

Tipo de documento DNI

Número de documento de identidad 07418845

Miembro del jurado 3

Nombres y apellidos

Tipo de documento DNI


Número de documento de identidad

Datos de investigación

Línea de investigación No aplica.

Grupo de investigación No aplica.

Agencia de financiamiento Sin financiamiento

Universidad Nacional Mayor de San Marcos


Ubicación geográfica de la
Latitud: -12.05777
investigación
Longitud: -77.08232

Año o rango de años en que se


2019 - 2020
realizó la investigación

URL de disciplinas OCDE Agricultura


https://purl.org/pe-repo/ocde/ford#4.01.01
Universidad Nacional Mayor de San Marcos
(Universidad del PERÚ, Decana de América)
FACULTAD DE CIENCIAS FÍSICAS
ESCUELA PROFESIONAL DE INGENIERÍA MECÁNICA DE FLUIDOS

ACTA DE SUSTENTACIÓN VIRTUAL DEL INFORME PROFESIONAL PARA OPTAR EL


TITULO PROFESIONAL DE INGENIERO MECANICO DE FLUIDOS, MODALIDAD DE
TITULACIÓN M-4 – POR EXPERIENCIA PROFESIONAL
Siendo las 14:00 p.m. del día martes 8 de noviembre de 2022, en la Sala de Sesión Virtual de la Facultad
de Ciencias Físicas, bajo la presidencia del Dr. Ing. ARÍSTIDES RAMÓN SERRUTO COLQUE y con la
asistencia de los miembros Ing. GUIDO AMÉRICO ROZAS OLIVERA y el Ing. MIGUEL ERNESTO
ARÁMBULO MANRIQUE, se dio inicio a la Sesión Pública Virtual de Sustentación del Informe
Profesional para optar el Título Profesional de Ingeniero Mecánico de Fluidos, mediante la Modalidad de
Titulación M-4 Por Experiencia Profesional, del Bachiller:

CARLOS IVAN VALERIO PONCE

El Presidente del Jurado Examinador dio lectura del Resumen del Expediente, e invitó al Bachiller
CARLOS IVAN VALERIO PONCE, a realizar la Exposición del Informe Profesional titulado “ESTUDIO
DE HIDRÁULICA FLUVIAL PARA LA PROTECCIÓN CONTRA INUNDACIONES
ENTRE EL TRAMO DESEMBOCADURA – LOCALIDAD PAULLO DEL RÍO CAÑETE,
PROVINCIA DE CAÑETE – DEPARTAMENTO DE LIMA”.
Concluida la exposición del candidato y luego de las preguntas de rigor por parte del Jurado Examinador,
el Presidente invitó al Bachiller abandonar momentáneamente la sala de sesión para dar paso a la
deliberación y calificación por parte del Jurado.
Al término de la deliberación del Jurado, se invitó al candidato a regresar a la sala de sesión para dar lectura
a la calificación final obtenida por el Bachiller, la misma que fue:
DIECISIETE (17)
El Presidente del Jurado Dr. Ing. ARÍSTIDES RAMÓN SERRUTO COLQUE, a nombre de la Nación y
de la Universidad Nacional Mayor de San Marcos, declaró al Bachiller CARLOS IVAN VALERIO
PONCE Ingeniero Mecánico de Fluidos.
Siendo las 14:58 p.m. del mismo día, se levantó la sesión.

Dr. Ing. ARÍSTIDES RAMÓN SERRUTO COLQUE Ing. GUIDO AMÉRICO ROZAS OLIVERA
PRESIDENTE DEL JURADO MIEMBRO DEL JURADO

Ing. MIGUEL ERNESTO ARÁMBULO MANRIQUE


MIEMBRO DEL JURADO

Firmado digitalmente por LOZANO


BARTRA Whualkuer Enrique FAU
20148092282 soft
Motivo: Soy el autor del documento
Fecha: 01.12.2022 18:55:17 -05:00

DR. WHUALKUER ENRIQUE LOZANO BARTRA


VICEDECANO ACADÉMICO FCF

Datos de la plataforma virtual institucional del acto de sustentación:


https://us06web.zoom.us/j/83936498351?pwd=WDdpY2pGc3piczl3M0VSWXRlZjNXZz09
ID: ID de reunión: 839 3649 8351
Código de acceso: 785878
Grabación archivada en el siguiente enlace:
https://drive.google.com/file/d/1DX2yjVdbkCwFLlIxbwBUUDYkPMQAEknH/view?usp=sharing

Ciudad Universitaria, Pabellón de Ciencias Físicas, Ca. Germán Amézaga N° 375. Lima, 1
Teléfono: 619-7000 anexo 3806 E-mail: [email protected]
Anexo 1
Informe de evaluación de Originalidad

1. Facultad de Ciencias Físicas


2. Escuela/Unidad de Posgrado
Escuela Profesional de Ingeniería Mecánica de Fluidos
3. Autoridad Académica que emite el informe de originalidad
El Director de la EPIMF
4. Apellidos y nombres de la autoridad académica
Sarango Julca Douglas Donal
5. Operador del programa Informático de similitudes
Dr. Ing. Douglas Sarango Julca
6. Documento evaluado (1)
“ESTUDIO DE HIDRÁULICA FLUVIAL PARA LA PROTECCIÓN CONTRA INUNDACIONES
ENTRE EL TRAMO DESEMBOCADURA – LOCALIDAD PAULLO DEL RÍO CAÑETE,
PROVINCIA DE CAÑETE – DEPARTAMENTO DE LIMA”
7. Autor del documento (2)
CARLOS IVAN VALERIO PONCE
8. Fecha de recepción del documento 30/09/22
9. Fecha de aplicación del programa informático de similitudes 6/10/22
10. Software utilizado
• Turnitin (X)
• Ithenticate
• Otro (especificar)
11. Configuración del programa detector de similitudes
• Excluye textos entrecomillados
• Incluye Bibliografía
• Excluye cadenas menores a 40 palabras
12. Porcentaje de similitudes según programa detector de similitudes3
7%
13. Fuentes originales de las similitudes encontradas4
14. Observaciones
Ninguna
15. Calificación de originalidad
• Documento cumple criterios de originalidad, sin observaciones (X)
• Documento cumple criterios de originalidad, con observaciones
• Documento no cumple criterios de originalidad
16. Fecha del Informe 15 /11 /22

Firmado digitalmente por SARANGO


JULCA Douglas Donal FAU
20148092282 soft
Motivo: Soy el autor del documento
Fecha: 15.11.2022 22:51:39 -05:00

DR. ING. Douglas Sarango Julca


DIRECTOR ( e ) DE LA EPIMF

1 Otro (especificar) Informe Profesional Modalidad de Titulación M-4


2 Apellidos y Nombres completos
3 En letras y números
4 Indicarlas en orden decreciente y su respectivo porcentaje
DEDICATORIA

Este trabajo lo dedico con todo mi amor y respeto a las personas que a pesar de las dificultades
siempre me acompañaron, alentaron y motivaron mi camino.

Mis padres Agustín y Maruja, que son mi inspiración; mis hermanos Judith, Gino y José, mis
sobrinas Naya y killary que son mi fortaleza y son la base de mi crecimiento profesional.

2
AGRADECIMIENTO

Agradecer a mi familia por su incondicional aliento y apoyo, a mi querida escuela profesional


Ingeniería Mecánica de Fluidos por brindarme las herramientas para mi vida profesional y a
los docentes de la facultad que se esmeraron en brindarme los conocimientos teóricos y
prácticos para desempeñarme en la vida profesional.

A las empresas Sigt S.A Ingenieros Consultores y M +K Ingenieros Asociados SAC por haberme
permitido laborar en sus instalaciones para poner en práctica lo aprendido en las aulas
universitarias y adquirir nuevos conocimientos .Debo Agradecer también al Ingeniero Tomas
Sánchez e Ingeniero Walter Chuan , por sus valiosas aportaciones y sugerencias para la
culminación de este informe profesional.

3
TABLA DE CONTENIDO

DEDICATORIA .......................................................................................................................... 2
AGRADECIMIENTO .................................................................................................................. 3
ÍNDICE DE FIGURAS .............................................................................................................. 7
RESUMEN .............................................................................................................................. 11
1.1 Introducción ............................................................................................................. 13
1.2 Contexto donde se desarrolló la experiencia. ....................................................... 14
1.2.1 Ubicación de la Zona de Estudio ................................................................ 14
1.2.2 Periodo de duración de la actividad ........................................................... 15
1.2.3 Datos de la empresa .................................................................................. 15
1.2.4 Mercado ..................................................................................................... 15
1.2.5 Misión ........................................................................................................ 16
1.2.6 Visión ......................................................................................................... 16
1.2.7 Organigrama de la Empresa ...................................................................... 16
1.3 Importancia del tema ............................................................................................... 16
1.4 Identificación del problema ..................................................................................... 17
1.5 Propuesta de Solución ............................................................................................ 17
1.6 Motivación para tratar el tema ................................................................................ 18
1.7 Objetivos del Trabajo .............................................................................................. 18
1.7.1 Objetivo Principal ....................................................................................... 18
1.7.2 Objetivos específicos ................................................................................. 18
1.8 Delimitación del informe ......................................................................................... 18
1.8.1 Delimitación temporal................................................................................. 18
1.8.2 Delimitación espacial ................................................................................. 18
1.8.3 Limitaciones ............................................................................................... 19
2. MARCO TEÓRICO ........................................................................................................... 20
2.1 Antecedentes: .......................................................................................................... 20
2.2 Marco conceptual. ................................................................................................... 21
2.2.1 Inundaciones: ............................................................................................ 21
2.2.2 Clasificación morfológica de Ríos .............................................................. 21
2.2.3 Erosión Fluvial ........................................................................................... 22
2.3 Bases teóricas Modelamiento Hidráulico............................................................... 22
2.3.1 Radio hidráulico en un canal muy ancho .................................................... 22
2.3.2 Tirante de diseño (B > 30m) ....................................................................... 23
2.3.3 Velocidad media ........................................................................................ 23
2.3.4 Ancho estable o amplitud de cauce............................................................ 24
2.3.4.1 Método de Blench – Altunin ................................................................. 24
2.3.4.2 Método de Julien y Wargadalam ......................................................... 24
2.3.4.3 Método de Altunin................................................................................ 25
2.3.4.4 Método de Pettits. ................................................................................ 26
2.3.4.5 Método de Manning Strickler. .............................................................. 26

4
2.3.4.6 Método de Simons y Henderson .......................................................... 27
2.3.5 Erosión General Transitoria ....................................................................... 27
2.3.6 Descripción del modelo numérico. ............................................................. 29
2.4 Bases teóricas Diseño Hidráulico .......................................................................... 31
2.4.1 Diques de enrocado ................................................................................... 31
2.4.2 Profundidad de la uña (P uña) ................................................................... 31
2.4.3 Borde libre (BL) .......................................................................................... 32
2.4.4 Altura del enrocado (H) .............................................................................. 33
2.4.5 Espesor mínimo del enrocado (T) .............................................................. 33
2.4.6 Ancho de la base (A).................................................................................. 33
2.5 Diámetro medio de la roca (D50) ............................................................................ 33
2.5.1 Método del U.S. Department of Transportation (U.S.D.T) .......................... 33
2.5.2 Método de Maynord ................................................................................... 35
2.5.3 Método de Isbash ...................................................................................... 35
3. INFORMACIóN BÁSICA - ESTUDIOS COMPLEMENTARIOS ........................................ 36
3.1 Información topográfica .......................................................................................... 36
3.2 Hidrología. ................................................................................................................ 37
4. METODOLOGÍA ............................................................................................................... 38
4.1 Modelamiento Hidráulico ........................................................................................ 38
4.1.1 Ubicación de los tramos a Intervenir .......................................................... 38
4.1.2 Topografía y cubierta del tramo de estudio ................................................ 39
4.1.3 Coeficiente del número de Manning (n) ..................................................... 41
4.1.4 Hidrograma de diseño ................................................................................ 45
4.1.5 Criterios empleados en el trazo de eje de estructuras de defensa ............. 45
4.1.6 Modelo Digital de Elevación (DEM) ............................................................ 54
4.1.7 Modelamiento IBER 2D.............................................................................. 55
4.2 Diseño Hidráulico .................................................................................................... 61
4.2.1 Criterios de diseño de los diques de enrocado. .......................................... 61
4.2.1.1 Trazado en perfil ................................................................................. 61
4.2.1.2 Profundidad de socavación.................................................................. 62
4.2.1.3 Altura del enrocado.............................................................................. 62
4.2.1.4 Diámetro medio de la roca (D50) . ....................................................... 62
4.2.1.5 Dimensiones transversales del dique de enrocado. ............................. 62
5. ANÁLISIS Y DISCUSIÓN DE RESULTADOS. ................................................................. 64
5.1 Modelamiento por inundación (sin estructuras de defensa) ................................ 64
5.1.1 Análisis para el Tr de 100 años .................................................................. 64
5.2 Modelamiento con estructuras de defensa ............................................................ 83
5.2.1 Simulación IBER 2D con diques para el Tr de 100 años ............................ 83
5.3 Trazos definitivos de eje de diques. ....................................................................... 85
5.4 Diseño de la Defensa Ribereña. .............................................................................. 89
5.4.1 Calculo Hidráulico del Tramo Dique 6 – (5+500 hasta 6+700) km .............. 90
5.4.1.1 Ubicación............................................................................................. 90

5
5.4.2 Cálculo de la velocidad media y número de froude. ................................... 91
5.4.2.1 Velocidad media .................................................................................. 91
5.4.2.2 Número de froude. ............................................................................... 91
5.4.3 Cálculo de la Profundidad de Socavación. ................................................. 92
5.4.3.1 Probabilidad anual en que se presenta el gasto de diseño. ................. 93
5.4.3.2 Diámetro de los granos de fondo (mm). ............................................... 94
5.4.3.3 Profundidad efectiva de socavación general del cauce principal. ........ 95
5.4.4 Cálculo de cota y Profundidad de Uña(P Uña). .......................................... 96
5.4.5 Cálculo de la altura del enrocado ............................................................... 97
5.4.5.1 Borde libre ........................................................................................... 97
5.4.5.2 Altura del enrocado (H)........................................................................ 99
5.4.6 Cálculo del diámetro medio del enrocado .................................................. 99
5.4.6.1 Maynord. ............................................................................................. 99
5.4.6.2 Método de Isbash .............................................................................. 100
5.4.6.3 Método del U.S. Department of Transportation (U.S.D.T) .................. 101
5.4.7 Diseño definitivo del Dique 6 de enrocado ............................................... 103
5.4.7.1 Espesor mínimo del enrocado (T) ...................................................... 103
5.4.7.2 Ancho de la base (A) ......................................................................... 104
6. CONCLUSIONES ........................................................................................................... 106
7. RECOMENDACIONES ................................................................................................... 107
8. REFERENCIAS BIBLIOGRAFICAS ............................................................................... 108
9. ANEXOS ........................................................................................................................ 109

6
ÍNDICE DE FIGURAS

Figura 1. Ubicación de la Cuenca del Rio Cañete. .................................................................................... 15


Figura 2. Organigrama de la empresa ....................................................................................................... 16
Figura 3. Clasificación Morfológica de Ríos ............................................................................................... 21
Figura 4. Erosión Fluvial ............................................................................................................................ 22
Figura 5. Radio hidráulico en un canal muy ancho. ................................................................................... 22
Figura 6. Erosión General .......................................................................................................................... 28
Fuente: Elaboración Propia........................................................................................................................ 28
Figura 7. Muestra una red de elementos triangulares a los que se les ha dado la forma de triángulos
rectángulos ................................................................................................................................................. 30
Figura 8. Selección típica de enrocado ...................................................................................................... 31
Figura 9. Resaltada cobertura LIDAR; Desembocadura rio Cañete – Localidad PAULLO. ...................... 36
Figura 10. Batimetría en Puente Clarita ..................................................................................................... 36
Figura 11. Se observa las curvas de nivel correspondiente del tramo Desembocadura rio Cañete –
Localidad Paullo ......................................................................................................................................... 37
Figura 12. Vista satelital del tramo en estudio “Desembocadura–Bocatoma Palo Herbay”. ..................... 38
Figura 13. Vista satelital del tramo en estudio “Bocatoma Palo Herbay-Localidad Paullo”. ..................... 39
Figura 14. Vista satelital de la cobertura sobre la superficie - Bocatoma La Fortaleza. ........................... 40
Figura 15. Bocatoma La Fortaleza. ............................................................................................................ 40
Figura 16. vista satelital de la cobertura - Bocatoma la Fortaleza. ............................................................ 41
Figura 17. Granulometría del cauce y Bosque inundable de caña brava,parte baja Bocatoma La
Fortaleza. ................................................................................................................................................... 42
Figura 18. Bocatoma Imperial. ................................................................................................................... 42
Figura 19. Bosques inundables y granulometría del cauce. Ambas fotos se tomaron parte baja de la
Bocatoma Imperial ..................................................................................................................................... 43
Figura 20. Sectorización de la superficie del cauce y áreas inundables . ................................................. 44
Figura 21. Sectorización de la superficie del cauce y áreas inundables . ................................................. 44
Figura 22. hidrograma de flujo de la estación Socsi . ................................................................................ 45
Figura 23. Trazo inicial propuesto tomando como referencia la Fluviomorfología del año 2003. ............. 47
Figura 24. Vista de trazos de diques colapsados en avenidas del año 2019. ........................................... 49
Figura 25. Trazos definitivos de ejes de diques para el año 2020-Bocatoma La Fortaleza. ..................... 49
Figura 26. Trazos en la Fluviomorfología del año 2009 (figura superior) y 2013 (figura inferior). ............ 50
Figura 27. Vista de diques colapsados y transición natural del ancho teórico del cauce. ........................ 52
Figura 28. Unión de trazos definitivos de ejes de diques de ambos proyectos en la zona de transición
natural del cauce (Bocatoma Imperial). ..................................................................................................... 53
Figura 29. Modelo digital de elevación en formato tiff de 0.5 m x 0.5 m de pixel. ..................................... 54
Figura 30. Modelo digital de elevación en formato tiff de 0.4 m x 0.4 m de pixel. ..................................... 55
Figura 31. Generación de la malla estructurada de 5 m ............................................................................ 56

7
Figura 32. Generación de la malla estructurada de 4 m de espaciamiento lateral. .................................. 57
Figura 33. Vista de condición de entrada y salida en el MDT. .................................................................. 58
Figura 34. Vista de condición de entrada y salida en el MDT. .................................................................. 59
Figura 35. Seccion Tipica ........................................................................................................................... 63
Figura 36. Mapa de inundaciones en IBER 2D del tramo “Desembocadura río Cañete - Bocatoma Palo
Herbay” para el Tr de 100 años. ................................................................................................................ 65
Figura 37. Vistas satelitales de zonas en riesgo -Bocatoma La Pinta. ...................................................... 66
Figura 38. Vistas satelitales de zonas en riesgo -Bocatoma La Pinta. ...................................................... 67
Figura 39. Puente Clarita. .......................................................................................................................... 68
Figura 40. Viviendas cerca de la ribera derecha aguas arriba del Puente Clarita. ................................... 68
Figura 41. Viviendas cerca en la ribera Izquierda aguas abajo del Puente Clarita. .................................. 69
Figura 42. Viviendas cerca de la ribera Izquierda aguas abaao del Puente Clarita. ................................. 70
Figura 43. Vistas satelitales y desde el terreno mismo del colapso continuo de los diques rústicos, aguas
abajo del Puente Clarita. ............................................................................................................................ 73
Figura 44. Mapa de inundaciones en IBER 2D del tramo “Puente Paullo - Bocatoma Palo Herbay” para el
Tr de 100 años. .......................................................................................................................................... 75
Figura 45. A, B y C vistas satelitales correspondientes a los años 2015, 2017 y 2019 respectivamente de
la zona en riesgo sector Lúcumo. D fotografía de la situación actual de las riberas en riesgo en el sector
Lúcumo....................................................................................................................................................... 78
Figura 46. Vista de la inundación según la simulación IBER 2D en el sector Lúcumo. ............................ 79
Figura 47. Vista de la inundación según la simulación IBER 2D aguas abajo y arriba del Puente Socsi . 81
Figura 48. Espigones existentes que protege el banco elevado, aguas abajo del Puente Socsi. ............ 82
Figura 49. Dique enrocado junto a un dique rústico que protege la llanura inundable por donde pasa la
carretera Cañete - Yauyos, aguas abajo del Puente Socsi ....................................................................... 82
Figura 50. Vista del modelamiento hidráulico IBER 2D definitivo en presencia de diques tanto
proyectados como existentes que necesitan elevación de corona............................................................ 84
Figura 51. Vista del modelamiento hidráulico IBER 2D definitivo en presencia de diques tanto
proyectados como existentes que necesitan elevación de corona............................................................ 85
Figura 52. Trazos definitivos de ejes de estructuras de defensa ribereña (líneas de color magenta), valles
vulnerables y diques existentes que necesitan elevación de corona. ....................................................... 87
Figura 53. Trazos definitivos de ejes de estructuras de defensa ribereña (líneas de color magenta), valles
vulnerables y diques existentes que necesitan elevación de corona. ....................................................... 89
Figura 54. Ubicación del dique 6. .............................................................................................................. 90
Figura 56. Gasto de diseño ........................................................................................................................ 93
Figura 57. Esquema de la socavación del cauce de rio. ........................................................................... 95
Figura 58. Calculo de la profundidad de socavacion dique 6 ................................................................... 96
Figura 59. Cotas y profundidad de la uña dique 6 ..................................................................................... 97
Figura 60. Altura del enrocado del dique 6 ............................................................................................... 99
Figura 61. Diametro medio del enrocado del dique 6. ............................................................................ 103
Figura 62. Dimensiones del dique de enrocado del dique 6 .................................................................... 104

8
Figura 63. Representación de las dimensiones del enrocado ................................................................. 105
Figura 64. Seccion tipica del enrocado .................................................................................................... 105
Figura 65. Leyenda del diseño ................................................................................................................. 105

9
ÍNDICE DE TABLAS

Tabla 1. Inversa de coeficientes de rugosidad según material de lecho de rio ................................... 23


Tabla 2. Factores del cauce según material ........................................................................................ 24
Tabla 3. coeficiente de K según el material del cauce ......................................................................... 25
Tabla 4. Valores de K para el método de Manning-Strickler ............................................................... 26
Tabla 5. Valores de m para el método de Manning-Strickler ............................................................... 27
Tabla 6. Valores de K1(Condiciones de fondo del cauce). .................................................................. 27
Tabla 7. Tamaño de malla en una serie de estudios de ingeniería y de investigación. En los tres
primeros se empleó el programa Iber 2D, en el resto el programa River 2D. ........................................... 30
Tabla 8. Coeficiente ø en función del caudal. ...................................................................................... 32
Tabla 9. Recomendación de borde libre .............................................................................................. 32
Tabla 10. Consideraciones para el factor de seguridad. ....................................................................... 34
Tabla 11. Caudales Máximos Río Cañete - Estación Socsi (m3/s). ...................................................... 37
Tabla 12. Coeficientes de Manning para cursos de agua naturales y llanuras de inundación. ............. 43
Tabla 13. Asignación del coeficiente de Manning (n) al cauce y áreas adyacentes. ............................ 44
Tabla 14. Asignación del coeficiente de Manning (n) al cauce y áreas adyacentes. ............................ 44
Tabla 15. Parámetros hidráulicos del Tramo en estudio. ...................................................................... 48
Tabla 16. Métodos de cálculo para el ancho estable (Be) del tramo en estudio ................................... 48
Tabla 17. Parámetros hidráulicos del tramo en estudio. ........................................................................ 51
Tabla 18. Métodos de cálculo para el ancho estable (Be) del tramo en estudio ................................... 51
Tabla 19. Ventana de Datos de procesamiento en Iber 2D ................................................................... 60
Tabla 20. Ventana de Datos de procesamiento en Iber 2D ................................................................... 61
Tabla 21. Longitud total y parcial de trazos en cada margen del tramo en estudio. ............................. 86
Tabla 22. Longitud total y parcial de trazos en cada margen del tramo en estudio. ............................. 88
Tabla 23. Longitud total y parcial de trazos en cada margen del tramo en estudio. ............................. 89
Tabla 24. Valores del coeficiente β ........................................................................................................ 93
1
Tabla 25. Suelos cohesivos y no cohesivos( x,y, )...................................................................... 94
1+ 
Tabla 26. Valores del coeficiente ø en función del caudal ..................................................................... 98
Tabla 27. Recomendación de borde libre .............................................................................................. 98
Tabla 28. Consideraciones para el factor de seguridad. ..................................................................... 102

10
RESUMEN

El presente trabajo de suficiencia profesional tiene por título “ESTUDIO DE HIDRÁULICA FLUVIAL
PARA LA PROTECCIÓN CONTRA INUNDACIONES ENTRE EL TRAMO DESEMBOCADURA -
LOCALIDAD PAULLO DEL RÍO CAÑETE, PROVINCIA DE CAÑETE- DEPARTAMENTO DE LIMA” el
mismo que se realizó con la finalidad de proponer defensas ribereñas, para mitigar las
inundaciones que sufren los distritos de San Vicente, Nuevo Imperial y Lunahuaná de la provincia
de Cañete.

El objetivo del presente trabajo fue realizar un estudio de Hidráulica Fluvial para Diseñar
estructuras de defensa para el control de inundaciones en el tramo “Desembocadura río Cañete
– Puente Paullo”, para un evento de máximo de 100 años. Asimismo, se identificó mediante un
modelamiento bidimensional usando el software IBER, las zonas vulnerables de inundación para
caudales 10,100 y 500 años lo que permitió determinar los tramos donde se ubicarán las
defensas ribereñas. Los resultados obtenidos en el modelamiento se corroboraron mediante
visitas técnicas de campo realizadas, se identificaron las zonas de riesgos de las áreas
expuestas a inundación y los bancos frágiles frente a la erosión que ocasionan los flujos en
máximas avenidas, así como las variables hidráulicas como los coeficientes de resistencia al
flujo del cauce y las planicies de inundación.

Se ha estimado el ancho estable del cauce en base a las formulaciones de Simons y Henderson,
Blench Lacey, Pettits, Altunin y Julien – Wargadalam, así como observaciones de campo y de
vistas satelitales. Se ha determinado que el ancho del rio en la zona de estudio es variable para
eventos que ocurren con un periodo de retorno de 100 años, con esta información se han
determinado anchos para trazar las defensas ribereñas que protegerán las poblaciones y los
campos de cultivo contra las inundaciones

El caudal de diseño del proyecto es de1064 m3/s para un periodo de retorno de 100 años. La
pendiente fue variable para cada tramo de dique que se encuentran en el orden del 1%, se
calcularon los parámetros hidráulicos tales como tirantes Hidráulico fue un promedio de 1.40 m
y espejo de agua varia 126 m- 450 m aproximadamente, profundidades de socavación, y
dimensiones de la estructura por cada tramo del proyecto, se muestran detalladamente en los
cuadros de cálculos en los anexos.

Finalmente se identificó el riesgo de áreas expuestas a inundaciones, mediante modelamientos


IBER 2D, para caudales de 625, 1064 y 1428 m3/s correspondientes a 10, 100 y 500 años de
periodo de retorno respectivamente. También se definieron trazos de eje, para estructuras de
defensa ribereña, cuya longitud total de la margen izquierda es de 14.513 km y de la derecha de
17.406 km.

11
ABSTRACT

The present work of professional sufficiency is entitled "STUDY OF FLUVIAL


HYDRAULICS FOR FLOOD PROTECTION BETWEEN THE STREAM "DESEMBOCADURA -
PAULLO LOCATION OF THE CAÑETE RIVER, PROVINCE OF CAÑETE- DEPARTMENT OF
LIMA" which was carried out with the purpose of proposing river defenses to mitigate the floods
suffered by the districts of San Vicente, Nuevo Imperial and Lunahuaná in the province of Cañete.

The objective of this work was to carry out a river hydraulics study to design defense
structures for flood control in the "Cañete River mouth - Paullo Bridge" section, for a return period
of 100 years. Likewise, a two-dimensional modeling using IBER software was used to identify the
areas vulnerable to flooding for 10, 100 and 500-year flows, which made it possible to determine
the sections where the river defenses will be located. The results obtained in the modeling were
corroborated by technical field visits, identifying the risk zones of the areas exposed to flooding
and the banks that are fragile to erosion caused by maximum flood flows, as well as the hydraulic
variables such as the flow resistance coefficients of the channel and the flood plains.

The stable channel width was estimated based on the formulations of Simons and
Henderson, Blench Lacey, Pettits, Altunin and Julien - Wargadalam, as well as field observations
and satellite views. It has been determined that the width of the river in the study area is variable
for events that occur with a return period of 100 years, with this information widths have been
determined to draw the riparian defenses that will protect populations and crop fields against
floods.

The design flow for the project is 1064 m3/s for a return period of 100 years. The slope
was variable for each levee section, which is in the order of 1%. Hydraulic parameters were
calculated, such as an average of 1.40 m and a water mirror of approximately 126 m to 450 m,
scour depths, and dimensions of the structure for each section of the project, which are shown in
detail in the tables of calculations in the annexes.

Finally, the risk of areas exposed to flooding was identified using IBER 2D modeling for
flows of 625, 1064 and 1428 m3/s corresponding to 10, 100 and 500 year return periods,
respectively. Axis lines were also defined for riparian defense structures, with a total length of
14.513 km on the left bank and 17.406 km on the right bank.

12
GENERALIDADES

1.1 Introducción

En el Perú son usuales los problemas de inundación, los cuales se deben a consecuencia
de falta de recubrimiento de vegetación, protección de los bancos de riberas, y falta de
descolmatación de los cauces de los ríos, eliminación de desmonte por parte de la población en
el lecho del rio, pero principalmente por eventos extraordinarios desencadenados por el
Fenómeno del Niño Costero, los cuales generan caudales de gran magnitud como los ocurridos
en los años 1972,1983,1998 y 2017.

El Fenómeno del Niño Costero se produce debido a cambios bruscos en la presión


atmosférica, cambios de velocidad del viento, elevación anormal de la temperatura de la
atmósfera y de las aguas del océano pacifico.

El reciente Fenómeno del Niño del 2017, ocasiono desborde del río Cañete afectando no
solo a una gran cantidad de pobladores y propiedades en los distritos de San Vicente de Cañete,
Nuevo Imperial y Lunahuaná, sino también a las estructuras hidráulicas, zonas agrícolas
infraestructura vial y edificaciones.

El río Cañete inundó zonas agrícolas tanto en la parte media como en la parte baja y
ocasionó la pérdida de parte de la plataforma de la vía Cañete – Yauyos en varios sectores. La
erosión excesiva de las márgenes del río Cañete ocasionó la caída de muros perimétricos y puso
en peligro las edificaciones construidas cerca de la margen derecha en el tramo entre Pacarán
y Lunahuaná. Además, ocasionó la caída de un puente colgante en Catapalla.

Ante lo mencionado el consorcio CISPDR-SIGT, realizo por encargo del Programa


Subsectorial de Irrigaciones (PSI) los estudios de preinversion a Nivel de Perfil “Formulación del
plan integral para el control de inundaciones y movimientos de masa de la cuenca del río Cañete.”
y de esta manera formular un plan de acción integral para el control de estos fenómenos
presentes en los cauces y quebradas de la cuenca del río Cañete.

Con el objetivo de disminuir el riesgo de inundación, pérdida de zonas agrícolas,


estructuras hidráulicas de riego, Infraestructuras e incluso humanas, se elaboró el proyecto
“ESTUDIO DE HIDRÁULICA FLUVIAL PARA LA PROTECCIÓN CONTRA INUNDACIONES ENTRE EL TRAMO
DESEMBOCADURA - LOCALIDAD PAULLO DEL RÍO CAÑETE, PROVINCIA DE CAÑETE- DEPARTAMENTO DE
LIMA “,se plantea el diseño de defensas ribereñas “Diques con Enrocados “, con la finalidad de
controlar las inundaciones entre el tramo “Desembocadura– Puente Paullo”, para un evento
máximo de 100 años .

Así mismo se usaron información de otras especialidades directamente relacionado al


proyecto como Topografía, Hidrología, Modelamiento Hidráulico y Geotecnia.

13
La metodología empleada para el diseño de estructuras de defensa ribereña se realizó a
partir de: identificación de sectores expuestos a riesgos por inundaciones, mediante
simulaciones IBER 2D; así como, identificación de riberas sometidas a fuerte y constante erosión
lateral, mediante observaciones de la migración del río en imágenes satelitales registradas
durante la mayor cantidad de años posibles (2003 – 2020). Luego, para una mejor decisión de
estos riesgos fue corroborado mediante visitas in situ, realizadas a las zonas vulnerables
identificadas a lo largo de las riberas (áreas inundables y bancos expuestos erosión). Las
inspecciones se realizaron durante el año 2020.

El criterio del diseño de las defensas ribereñas consistió en diques de enrocados, fueron
por la disponibilidad del material del lecho del río para la conformación del dique, cercanía de
canteras para el enrocado; además, por ser estructuras que no disminuyen el área hidráulica,
cuyo factor es elemental sobre todo para los tramos estrechos.

1.2 Contexto donde se desarrolló la experiencia.

1.2.1 Ubicación de la Zona de Estudio

El proyecto se ubica en la Cuenca del río Cañete, en el tramo Desembocadura del río
Cañete- Localidad de Paullo progresiva 0 + 000 km al 31 + 100 km, abarcando distritos de San
Vicente de Cañete, Nuevo Imperial y Lunahuaná, provincia de Cañete, departamento de Lima.

Los límites de la cuenca del río Cañete son:

• Norte : Cuenca del rio Mantaro.

• Sur :Cuenca de la Quebrada Topara.

• Este :Cuenca del rio Mantaro.

• Oeste : Cuenca del rio Mala -Océano Pacifico

El flujo de aguas vierte en el Océano Pacífico. Esta cuenca tiene una superficie de 6
049.47 km2, un perímetro de 523.78 km, siendo la longitud del río Cañete es de 234.32 km.

14
Figura 1. Ubicación de la Cuenca del Rio Cañete.
Fuente: Elaboración Propia.

1.2.2 Periodo de duración de la actividad

El periodo de duración del proyecto duro aproximado de 18 meses, agosto 2019 hasta
diciembre de 2020.

1.2.3 Datos de la empresa

Razón Social : Sigt S.A Ingenieros Consultores.

Domicilio Legal: Av Macchu Picchu 376, Urb.San Juan Bautista-Chorrillos- Lima

Teléfono : 51-1-2549518

Ruc : 20294921312

1.2.4 Mercado

En la actualidad Sigt Ingenieros Consultores es una empresa que brinda servicios de


consultoría en ingeniería realiza estudios de Geomática (geodésicos y topográficos), minería,
obras viales e Hidráulica e Irrigaciones.

15
1.2.5 Misión

Elaborar proyectos de ingeniería con eficiencia y manteniendo los más altos estándares
de calidad, garantizando la excelencia de nuestros servicios, los cuales son respaldados por la
experiencia y capacitación constante de nuestro equipo de profesionales. Contar con los
recursos tecnológicos de última generación para mantenerse a la vanguardia y brindar servicios
competentes para alcanzar la satisfacción de nuestro cliente

1.2.6 Visión

Ser una empresa líder en brindar servicios de Consultoría en Ingeniería.

1.2.7 Organigrama de la Empresa

Gerente Dpto. de
financiero contabilidad
Gerente de Dpto. de
Dpto. tecnico
administracion logistica
Gerente
Director Contro de Dpto. de control de
general Gerente tecnico calidad
calidad
Departamento de Dpto. de costos y
mantenimiento presupuestos
Gerente de
operaciones Sub-gerencia de
equipos

Figura 2. Organigrama de la empresa

1.3 Importancia del tema

Debido a los periodos de máximas avenidas (enero-diciembre) asociado al Fenómeno


del Niño, se produce inundaciones entre el tramo “Desembocadura del rio Cañete – Localidad
de Paullo”, causando daños a las zonas agrícolas, estructuras hidráulicas, infraestructura y
incluso pérdidas humanas.

16
1.4 Identificación del problema

El tramo de estudio se encuentra en la parte baja y media de la Cuenca Cañete,


abarcando los distritos de San Vicente, Nuevo Imperial y Lunahuaná afectando:

• Centros Poblados: Santa Teresa, Urbanito Cruz, Pedro de la cruz, Boca del Rio,
ubicados en la margen derecha y La margen izquierda es vulnerable a inundaciones
entre la bocatoma Palo Herbay y el centro poblado Isla Alta están expuestas a
desbordes y inundaciones como consecuencia de intensas precipitaciones en la
cabecera de la cuenca, especialmente asociados al Fenómeno del Niño y puesto que
los diques actuales no pueden cumplir su función de evitar desbordes para un caudal
de 1064 m3/s, para un tiempo de retorno de 100 años.

• Puentes: El puente antiguo de Socsi y el puente Clarita (antigua Panamericana Sur)


son los casos más destacados. El primero por ser una estructura baja, que puede ser
colmatada en una gran avenida del río y el otro, que colapsó parcialmente en el FEN
Costero 2017 ha sido rehabilitado en un tramo, provisionalmente con un puente
metálico Bailey.

• Bocatomas: En ambas márgenes, existen estas estructuras para fines de riego del
sector agrícola como son: Nuevo Imperial, La Pinta, Fortaleza, Ungará (margen
derecha) y Palo Herbay, Ramadilla (margen izquierda) y otras rústicas. Tanto la
bocatoma Nuevo Imperial como La Pinta tienen un nivel bajo de vulnerabilidad. Por
el contrario, las demás bocatomas tienen vulnerabilidad media, debido a la posibilidad
de ser colmatadas, lo que significaría la pérdida de cultivos.

1.5 Propuesta de Solución

Se identificará las áreas de protección, usando simulaciones de inundación con el


programa IBER para tiempos de retorno de 10,50 100 y 500 años, de igual modo se corroborará,
mediante visitas in situ, los riesgos de las áreas expuestas a inundación y de bancos frágiles
frente a la erosión que ocasionan los flujos en máximas avenidas. Se definirá los trazos de las
defensas ribereñas y las variables hidráulicas como los coeficientes de resistencia al flujo del
cauce y las planicies de inundación

Se propondrá una medida estructural de protección, con fines de proteger las zonas
agrícolas, estructuras hidráulicas, infraestructuras y evitar pérdidas humanas.

17
1.6 Motivación para tratar el tema

El informe que se presenta, es parte de la experiencia profesional del autor que participó
en la elaboración del estudio de Hidráulica Fluvial, fue una oportunidad para aplicar lo aprendido
en las aulas universitarias y los conocimientos adquiridos durante los 3.5 años de experiencia
laboral.

1.7 Objetivos del Trabajo

Dentro de los objetivos que se presentaran en el presente estudio son:

1.7.1 Objetivo Principal

• Diseñar estructuras de defensa para el control de inundaciones en el tramo


“Desembocadura río Cañete – Puente Paullo”, para un periodo de retorno de 100
años.

1.7.2 Objetivos específicos

• Identificar, mediante modelamiento bidimensional IBER, las zonas vulnerables a


inundación para caudales correspondientes a 10,100 y 500 años.

• Corroborar, mediante visitas in situ, los riesgos de las áreas expuestas a inundación
y de bancos frágiles frente a la erosión que ocasionan los flujos en máximas avenidas.

• Definir los trazos de las defensas ribereñas y las variables hidráulicas como los
coeficientes de resistencia al flujo del cauce y las planicies de inundación.

1.8 Delimitación del informe

1.8.1 Delimitación temporal

El presente informe de la experiencia profesional se realizó en el periodo de marzo 2019


hasta diciembre 2020

1.8.2 Delimitación espacial

El presente informe solo abarca el componente A de la “Formulación del plan integral


para el control de inundaciones y movimientos de masa de la cuenca del río Cañete. Que
consiste en la “Creación del Servicio de Protección contra Inundaciones del río Cañete entre el
Tramo Desembocadura del rio Cañete – Localidad de Paullo en los Distritos de San Vicente de
Cañete, Nuevo Imperial y Lunahuaná de la Provincia de Cañete”.

18
1.8.3 Limitaciones

• El CONSORCIO SALZGITTER-SISA a través del “Programa de Protección de Valles y


Poblaciones Rurales Vulnerables ante Inundaciones” estaba elaborando en el año 2019-
2020 el expediente técnico definitivo “Instalación - implementación de medidas de
prevención para control de desbordes e inundaciones del rio Cañete - Provincia de
Cañete - Departamento de Lima” , como el estudio estaba avanzado habían definido
sus ejes de diques , en algunos tramos de nuestro área de influencia hubo la necesidad
de acoplarse a los trazos de este proyecto, por disposición del cliente contratante
Programa Subsectorial de Irrigaciones (PSI).

• El presente informe profesional, solo comprende los trazos y diseños del dique de
protección de los tramos críticos de las zonas de inundación, más no se considera
estudios relacionados como: paisajismo, reforestación, ambiental, ni costos y
presupuestos.

19
2. MARCO TEÓRICO

2.1 Antecedentes:

Agencia de cooperación Internacional del Japón “JICA” (2013). Proyecto de


Inversión a Nivel de Factibilidad: Instalación-Implementación de medidas de prevención para el
control de desbordes e inundaciones del rio Cañete-Provincia Cañete-Departamento de Lima
(Código SNIP N° 184598),tuvo como objetivo la protección del valle cañete ante el riesgo de
inundación por desbordes del rio cañete de los distritos de San Vicente ,Nuevo Imperial y
Lunahuaná, asimismo se enfocó el estudio en la parte baja de la cuenca, se realizó simulaciones
hidráulicas utilizando el programa HEC-RAS 2D para periodos de retorno de 10 ,25,50 años,
se determinó trabajar con un Caudal de diseño de un tiempo de retorno de 50 años.

A través de aforo y análisis de precipitaciones se obtuvo caudales de diseño 753 m3/s,y


2175 m3/s respectivamente ,se determinó trabajar con el caudal de diseño 2175 m3/s,se
determinó el ancho estable es 232 m ,se planteó las defensas ribereñas en los puntos críticos
diques con enrocados teniendo una longitud de 8km.

Tito Quispe, Y. (2017). Modelamiento Hidráulico del Rio Cañete Sector Puente Socsi-altura
Puente colgante (9 km), con fines de diseño de Defensas Ribereñas. (Tesis de Pregrado),
Universidad Nacional Agraria La Molina, Carrera de Ingeniería Agrícola, Lima-Perú.

Extracto: La finalidad de esta tesis es identificar áreas vulnerables entre los tramos
Puente Socsi-Puente colgante, tenido una longitud de 9km, de modo que realizo el modelamiento
hidráulico para eventos de máximas avenidas utilizando el programa HEC-RAS ,uso datos de
la estación hidrometría de Socsi ,calculando que el caudal de diseño para un Tiempo de retorno
de 50 años es 757.53 m3/s ,determinó que el ancho estable en este tramo es 80 m ,obtenido
con el método de Altunin, planteo obras de protección en la riberas diques con enrocados y
espigones de gaviones .

El CONSORCIO SALZGITTER-SISA (2019) a través del “Programa de Protección de


Valles y Poblaciones Rurales Vulnerables ante Inundaciones expediente técnico “Instalación -
implementación de medidas de prevención para control de desbordes e inundaciones del rio
Cañete - Provincia de Cañete - Departamento de Lima” .Determinó su caudal máximo para un
tiempo de retorno de 50 años 880 m3/s y propone para el control de inundaciones Dique con
enrocado, considera un borde libre 1.0 m, aplicando las normas aplicadas en Japón, menciona
que el nivel de fondo de la uña antisocavantes es posesionado, considerando la socavación
calculada medida desde el fondo de la parte baja de la sección del río.

20
2.2 Marco conceptual.

2.2.1 Inundaciones:

Son un fenómeno frecuente que sucede en el Perú, que consiste en el incremento de lluvias de
mayor duración e intensidad, produciendo que el caudal sobre pase la caja hidráulica del río, dando como
resultado desbordamientos. Según Rocha (1998):

Una inundación, en cambio es el desbordamiento de un rio por la incapacidad del cauce para
contener el caudal que presenta. La inundación es, pues un fenómeno de tipo hidráulico, prueba
de ellos es que pueden ocurrir inundaciones sin que haya crecidas o un evento hidrometeorológico
extraordinario. (pp.10-11).

En conclusión, no es necesario que se presente un evento extraordinario, con que tenga una caja
hidráulica deficiente el río, se producirá inundaciones.

2.2.2 Clasificación morfológica de Ríos

Existen tres tipos de ríos con diferentes alineación y geometría:

• Ríos rectos: Inusualmente se presenta en tramos cortos. “Prácticamente no existen ríos


rectos en la naturaleza, A veces sucede que existe un sistema de encausamiento recto,
constituido con diques paralelos, pero dentro de él ,para caudales menores que el de
diseño, el rio desarrolla su propia sinuosidad”(Rocha,1998,p.187).

• Ríos trenzados: son ríos que presentan cauces secundarios, en el cauce principal.
“Son cauces Inestables en el sentido de que una crecida puede cambiarlos
considerablemente, De este modo se llaman divagantes porque un brazo principal
puede encontrase tan pronto en lugar como otro”(Juan P,2001,p.27).

• Ríos meandros: Son ríos que presentan curvas continuas en su cauce, se generan por
la erosión y sedimentacion de sus margenes.

Figura 3. Clasificación Morfológica de Ríos

21
2.2.3 Erosión Fluvial

Es un fenómeno destructivo de los materiales de la corteza terrestre debido, a que implica


fracturamiento, fisuramiento, alteración física y química, los agentes erosivos son el agua, viento.
La erosión fluvial es un problema permanente, que realiza los flujos de agua sobre el suelo.

Figura 4. Erosión Fluvial

2.3 Bases teóricas Modelamiento Hidráulico

2.3.1 Radio hidráulico en un canal muy ancho

“Cuando el ancho (b) de un canal o río es mucho mayor que el tirante, se dice que es un
canal muy ancho. Esto permite hacer un cálculo más rápido y fácil del radio
hidráulico”(Rocha,2007,p.10).

Figura 5. Radio hidráulico en un canal muy ancho.

A = by
P = b + 2y
 by  y
R= =
 b + 2 y  1+ 2 y
b

En un canal muy ancho la relación b/y es muy pequeña, por lo que, se puede considerar
que:

R=Y

22
2.3.2 Tirante de diseño (B > 30m)

Se procede a calcular el tirante de diseño en función al ancho estable mayor a 30 metros,


para lo cual se ha empleado la fórmula de Manning - Strickler:

3
 Q 5
y= 0.5 
 ksBS 
Dónde:

B: Ancho efectivo del cauce (m)


Q: Caudal de diseño (m3/s)
Ks: Inversa del coeficiente de rugosidad
S: Pendiente media del tramo

Tabla 1. Inversa de coeficientes de rugosidad según material de lecho de rio

Inversa del coeficiente de Rugosidad (Ks)


Valores para Ks par Cauces Naturales Ks
Cauce con fondo solido sin irregularidades 40
Cauces de rio con acarreo irregular 33-35
cauces de ríos con vegetación 30-35
Cauces naturales con derrubio e irregularidades 30
Cauces de rio con fuerte transporte de acarreo 28
Torrentes con piedras de tamaño de una cabeza 25-28
Torrentes con derrubio grueso y acarreo móvil 19-22

2.3.3 Velocidad media

Obteniendo el tirante de diseño, se calcula la velocidad promedio, según la fórmula de


Manning-Strickler. La fórmula es una evolución de la fórmula de Chézy para el cálculo de la
velocidad del agua en canales abiertos, propuesta por el ingeniero irlandés Robert Manning en
1889:

V =K R 2/3S1/2
m s h

Dónde:

Vm: Velocidad media (m/s)


Q: Caudal de diseño (m3/s)
R: Radio hidráulico (m)

23
Ks: Inversa del coeficiente de rugosidad
S: Pendiente media del tramo

2.3.4 Ancho estable o amplitud de cauce

Es cuándo el cauce natural presenta las condiciones de equilibrio, cuando su lecho no se


erosiona, ni se sedimenta. La morfología del río es originada por el desnivel de las cotas de
terreno y la amplitud es debido a la cantidad de material en arrastre, produciendo diversas formas
del cauce natural y anchos variables. A continuación, se detallan los métodos:

2.3.4.1 Método de Blench – Altunin

Está basado en dos parámetros (Fb: factor de fondo y Fs: factor de orilla) y consideran
la aglomeración de los sedimentos trasladado en suspensión, el espesor de las partículas de
fondo y la resistencia de las orillas a ser erosionadas.

QxFb
B = 1.81
Fs

Tabla 2. Factores del cauce según material

FACTOR DE FONDO (FB) FACTOR DE ORILLA (FS)

0.80-MAT.FINOS (DM<0.50 mm) 0.10-MAT.SUELTOS


1.20-MAT.GRUESOS(DM<0.50 mm) 0.20-MAT.LIGERAMENTE COHESIVOS
0.30-MAT.COHESIVOS

2.3.4.2 Método de Julien y Wargadalam

En el año 2002 Julien y Wargadalam mejoraron los estudios de anchos estables (B = W).

Plantea lo siguiente:

2+ 4 m −4 m −1− 2 m
W = 1.33Q 5+ 6 m
ds 5+ 6 m
S 5+ 6 m

1
m=
 12.2h 
ln  
 ds 

24
Dónde:

W: Ancho efectivo del cauce asociado al caudal (m)

Q: Descarga de diseño (m3/s)

ds: Diámetro medio del material del lecho de río (m)

h: Calado medio del cauce (m)

2.3.4.3 Método de Altunin

Este método se recomienda para cauces con material no cohesivos (grueso, bolonería).
La ecuación de este método es la siguiente:

A Q
B=
S 0.2

En la ecuación propuesta para evaluar el coeficiente A, “n” es el valor de rugosidad y K


valor que varía de acuerdo a la resistencia de las orillas.

( 3/ ( 3+5 m ) )
A = ( n  K 5/3 )

Los valores “K” a adoptarse dependen de las características del cauce, el cual se
muestran en la Tabla 3 :

Tabla 3. coeficiente de K según el material del cauce

Valores del
Material de la Orilla del Rio
Coeficiente
3a4 “Material de cauce muy resistente”.
16 a 20 “Material Fácilmente erosionable”.
8 a 12 “Material aluvial”.
10 “ Valor practico”.

Fuente: Ven Te Chow

Los valores de “m” son de acuerdo a:

0.5 Ríos de Montaña


m 0.7 Cauces arenosos
1 Cauces aluviales

25
2.3.4.4 Método de Pettits.

El método de Pettit se diferencia de otros métodos, incluye el perímetro mojado y el radio


hidráulico en lugar del ancho y calado.

B = 4.44  Q0.5
Dónde:

B : Ancho efectivo del cauce asociado al caudal (m)

Q : Descarga de diseño (m3/s)

2.3.4.5 Método de Manning Strickler.

El método considera el tipo de río basándose en un criterio geográfico geomorfológico,


así que establece tres tipos de ríos: zona montañosa, zona intermedia y de planicie.

3
B = ( Q1/2  S 1/5 )( n  K 5/3 ) 3+5m
Donde:

Q : Caudal de diseño (m3/s)


B : amplitud del cauce (m)
S : Pendiente media del tramo
n : Valor de rugosidad de Manning
K : Valor de material del cauce (tabla N° 4)
m : Coeficiente de tipo de río (tabla N° 5)

Tabla 4. Valores de K para el método de Manning-Strickler

Coeficiente de material de cauce K


Material de cauce muy resistente 3a4
Material fácilmente erosionable 16 a 20
Material aluvial 8 a 12
Valor práctico 10

Fuente: Ven Te Chow

26
Tabla 5. Valores de m para el método de Manning-Strickler

COEFICIENTE DE TIPO DE RIO m


Para ríos de montaña 0.50
Para cauces arenosos 0.70
Para caudales aluviales 1.00

Fuente: Ven Te Chow

2.3.4.6 Método de Simons y Henderson

Este método es la aplicación de los ríos (India y Estados Unidos) considera un fondo del
lecho de arena, y las orillas de material cohesivos (arcilla, limo) o no cohesivos(gravas).
1
B = k1  Q 2

Donde:

Q : Caudal de diseño (m3/s)


B : amplitud estable del cauce (m).
K1 : Coeficiente de material del cauce (tabla N° 6)

Tabla 6. Valores de K1(Condiciones de fondo del cauce).

Condiciones de Fondo de Rio K1


Fondo y orillas de arena 5.7
Fondo de arena y orillas de material cohesivo 4.2
Fondo y orillas de material cohesivo 3.6
Fondo y orilla de grava 2.9
Fondo de arena y orillas de material no cohesivo 2.8

2.3.5 Erosión General Transitoria

Es un descenso temporal del nivel del lecho producida durante la ocurrencia de una
avenida. Existen varias fórmulas para la estimación de la erosión general , para el presente
estudio se utilizó el método Lischtvan – Lebediev. “El método Lischtvan – Lebediev es válido
para estimar la erosión general en suelo granular y en suelo no cohesivo” (Campaña,2002).

27
Figura 6. Erosión General
Fuente: Elaboración Propia.

1
 5
1+ x
  3
Hs =  
H o
 0.68    d 0.28 
 m 
 

Qd
= 5
Be  Y 3
m

Hs : Tirante después de la erosión (m)


Ho : Tirante máximo antes de la erosión (m)
Qd : Descarga de diseño (m3/s)
Be : Ancho efectivo del cauce asociado a Qd (m)
Ym : Tirante medio antes de la erosión
x : Coeficiente que depende del "dm" de granos
dm : Diámetro medio de las partículas del material del lecho (mm)
β : Coeficiente de frecuencia
ds : Profundidad socavada (m).

28
2.3.6 Descripción del modelo numérico.

El modelo elegido para la elaboración de las simulaciones es el modelo IBER versión


2.5.2, fue desarrollado Grupo de Ingeniería del Agua y del Medio Ambiente, GEAMA (Universidad
de A Corufia, UDC) y el Instituto FLUMEN (Universidad Politecnica de Catalunya, UPC, y Centro
Intemacional de M6todos Numéricos en IngenierIa, CIMNE). Bladé et al. (2014) describen, de
manera sucinta, las consideraciones teóricas y los métodos empleados para resolver las
ecuaciones de movimiento.

El Iber tiene un módulo hidrodinámico que determina las ecuaciones de St. Venant
bidimensionales, incorporando los efectos de la turbulencia y rozamiento superficial por viento.

Por otro lado, la resistencia al flujo se calcula con la ecuación de Manning, como se
presenta debajo:

n 2U x / U 2 /
 b, x =  gh 4
3
h
n 2Uy / U 2 /
 b, y =  gh 4
h3

Donde: h es el tirante, Ux y Uy son las velocidades horizontales promediadas en


profundidad, g es la aceleración de la gravedad, r es la densidad del agua, Zb es la cota del
fondo, ts es la fricción en la superficie libre debida al rozamiento producido por el viento, t b es la
fricción debida al rozamiento del fondo y nt es la viscosidad turbulenta.

“Las ecuaciones de aguas someras y las del modelo k-ε se resuelven mediante el método
de volúmenes finitos para mallas bidimensionales no estructuradas”. Los esquemas numéricos
utilizados en Iber son especialmente apropiados para la modelización de cambios de régimen y
de frentes seco-mojado (frentes de inundación).

La discretización del dominio espacial se realiza con volúmenes finitos en mallas no


estructuradas, pueden ser mallas no estructuradas( triangulares y cuadrangulares) y
estructuradas. La Figura 7, que se presenta debajo, muestra una malla de cálculo formada por
una Red Irregular de Triángulos Rectángulos (RTIN).

29
Figura 7. Muestra una red de elementos triangulares a los que se les ha dado la forma de triángulos rectángulos
Fuente: Blade et al. 2014.

Tamaño de los elementos de la malla para modelación hidráulica 2D

La incorporación de modelos bidimensionales en estudios fluviales es relativamente


reciente. Esto se debe, en gran parte, a que el hardware necesario para poder ejecutar estos
programas no estaba disponible, en principio, y cuando estuvo disponible, el precio fue elevado.
Posteriormente, el costo de los equipos fue disminuyendo. En tiempos recientes la disponibilidad
de computadoras personales con un hardware apropiado (por ejemplo, una computadora i9 de
décima generación, con procesadores gráficos adecuados) permite ejecutar una simulación de
un tramo de río en un plazo razonable (entre horas y dos días). El tiempo de ejecución aumenta
con la densidad de las mallas para un mismo tramo. Una simulación adecuada se logra con un
balance entre la precisión de las salidas y el tiempo de ejecución. A manera de ejemplo, se
incluyen dimensiones de la malla que se han empleado en estudios recientes en la siguiente
tabla.

Tabla 7. Tamaño de malla en una serie de estudios de ingeniería y de investigación. En


los tres primeros se empleó el programa Iber 2D, en el resto el programa River 2D.

Dimensiones
Ancho Tipo de Nivel de
Nº Río País de la malla
medio (m) estudio estudio
(m x m)
1 Arroyo de la Almucera España 250 25 x 25 Inundación Preliminar
2 Ranchería Colombia 500 23.4 x 23.4 Inundación Preliminar
1 x1 a 10 x
3
Ctalamuchita Argentina 35 10 Puentes Ingeniería
4 Ica Perú 60 2.5 a 5 Inundación Ingeniería

30
Se puede emplear el programa Iber para las siguientes aplicaciones: Simulación del flujo
en superficie libre en cauces naturales, estimar zonas inundables, cálculo hidráulico de
encauzamientos, cálculo hidráulico de redes de canales en superficie libre, cálculo de corrientes
de marea en estuarios, estabilidad de los sedimentos del lecho.

En este estudio el programa Iber se empleará fundamentalmente para definir la extensión


de las áreas de inundación y, además, estimar los tirantes y velocidades de una manera más
detallada que en un modelo unidimensional.

El transporte de sedimentos, tanto en fondo como en suspensión, se calcula mediante


ciertas ecuaciones predefinidas. El transporte en fondo se puede calcular tomando en cuenta
los siguiente métodos o información: umbral de movimiento de Shields, Meyer- Peter- Müller,
Van Rijn, etc. El transporte en suspensión considera la formula : Van Rijn, Smith McLean,
Ariathurai.

2.4 Bases teóricas Diseño Hidráulico

2.4.1 Diques de enrocado

El dique se conforma con material propio del rio o de préstamo, las rocas se acomodan
en cara húmeda del dique y así mismo se conforma la uña en el fondo del cauce, ubicados en la
parte baja del talud.

Figura 8. Selección típica de enrocado


Fuente: Elaboración Propia.

2.4.2 Profundidad de la uña (P uña)

La profundidad de socavación establece el nivel de cimentaciones de las estructuras.


Esta profundidad es de acuerdo al material del lecho del cauce del rio, el cual se considera en el
cálculo. “Para el caso de estructuras longitudinales de defensa ribereñas, el pie de éstas en su
cara húmeda está sujeta a socavación. Siendo necesario construir estructuras antisocavantes

31
que den estabilidad a las obras de gravedad, caso de diques enrocados, muros de gaviones,
etc.” (Rubén,1998, p.85). Por tanto, la altura socavada, equivalente a la profundidad de uña del
enrocado.

P uña = Cota de fondo del lecho - Cota de socavación

2.4.3 Borde libre (BL)

Para el borde libre se consideró las sugerencias de investigadores peruanos Campaña


(2004) . “El borde libre permite controlar la variación instantánea del caudal por disminución de
la velocidad y elevación del tirante”(Teran,1998,p.36).

Vm2
BL =  
2 g

H: Altura del enrocado (m)


BL: Borde libre (m)
Vm: velocidad media del flujo (m/s)
g: Aceleración de la gravedad (m/s2)
: Coeficiente en función de la máxima descarga.
El coeficiente (ø) en función del caudal .

Tabla 8. Coeficiente ø en función del caudal.

Caudal Máximo (m3/s) Coeficiente (ø )

3000-4000 2.00
2000-3000 1.70
1000-2000 1.40
500-1000 1.20
100-500 1.10

Fuente: “Diseño de defensas Ribereñas-Rubén Terán”.

Tabla 9. Recomendación de borde libre

Descarga de diseño (m3/s) Borde Libre (m)


Menos de 200 0.6
200-500 0.8
500-2000 1.08

Fuente: Campaña(2004).

32
2.4.4 Altura del enrocado (H)

Es la diferencia entre la cota del nivel de flujo y socavación, se suma el borde libre. El
valor de la cota del flujo y la velocidad se obtienen del simulamiento IBER 2D.

H = Cota de flujo - Cota de socavación + BL

2.4.5 Espesor mínimo del enrocado (T)

Las consideraciones para tener en cuenta son:

• No debe ser menor que el tamaño máximo de la roca a usarse en el enrocado.

• No debe ser menor de 1.5*D50.

• Por razones constructivas no deberá ser menor de 30 cm.

2.4.6 Ancho de la base (A)

El ancho base está en función de la profundidad de socavación (P uña), por lo que, su


valor se determina:

A= 1.5 P uña

A : Ancho del enrocado en la base (m).

P uña : Profundidad de socavación o de la uña del enrocado (m).

2.5 Diámetro medio de la roca (D50)

Hay varias metodologías que estiman el diámetro medio de la roca, de los cuales se ha
considerado tomar en cuenta los siguientes:

2.5.1 Método del U.S. Department of Transportation (U.S.D.T)

D50 = C  d 50t

 0.00594  Vm3 
D150 = c   
 ym  k1
0.5 1.5

Dónde:

D50: Diámetro medio del enrocado de protección

33
K1: Coeficiente de los taludes y ángulos de reposo del enrocado de protección
Ym: Tirante medio del flujo al pie del enrocado
V: Velocidad media del flujo al pie del enrocado
C: Coeficiente de corrección.
0.5
  sin 2   
K1 = 1 −  2  
  sin   

Dónde:

ɸ: Ángulo de reposo del material del enrocado.


Θ: Ángulo del enrocado con la horizontal.
2.12
Csg =
( ys − 1)
1.5

1.5
FS
C =(
sf )
1.2
C = Csg  Csf

Csg : Coeficiente de corrección por el peso específico


Csf : Coeficiente de corrección por factor de seguridad
ɣ𝑠 : Gravedad específica del material del enrocado
FS : Factor de seguridad, según las siguientes consideraciones:

Tabla 10. Consideraciones para el factor de seguridad.

Condición Rango del Fs

Flujo uniforme, tramos rectos 0 medianamente curvos (radio


de la curvalancho de! cauce > 30). Mínima influencia de 1.0-1.2
impacto de sedimentos y material flotante.

Flujo gradualmente variado, curvatura moderada (10 < radio


) de la curva/ancho del cauce < 30). Moderada de impacto de 1.3-1.6
sedimentos y material flotante.

Flujo rápidamente variado, curvas cefradas (radio de la


curvalancho del cauce < 10), flujos de alta turbulencia, flujo
1.6-2.0
de turbulencia mixta en estribos de puentes. Efecto
significativo de impacto de sedimentos y material flotante.

34
2.5.2 Método de Maynord

“Define el cálculo del diámetro de la roca en función al diámetro medio, usando


constantes según el talud y según la forma en planta del río”(Cconisilla,2017,p.19). Utiliza los
siguientes fórmulas y parámetros.

D50 = t  C1  F 3

C2  v
F=
g y

D50 : Diámetro medio de las rocas

y : Profundidad media del flujo

V : Velocidad media del flujo.

F : Número de Froude

Los valores recomendados de C1 y C2 se muestran a continuación:

0.28 Fondo Plano


C1 0.28 Talud 1V:3H
0.32 Talud 1V:2H
1.5 Tramos en curva
C2 1.25 Tramos rectos
2 Extremos de espigones

2.5.3 Método de Isbash

Este método se basa en la siguiente expresión:

V
D50 = 0.59 
W g

kr − ka
w=
Ka

D50 : Diámetro medio de la roca

V : Velocidad media del flujo

kr : peso específico de la roca

ka : peso específico del agua

g : Aceleración de la gravedad

35
3. INFORMACIÓN BÁSICA - ESTUDIOS COMPLEMENTARIOS

3.1 Información topográfica

El área de Topografía de la empresa usó Tecnología Lidar Scanner para el levantamiento


topográfico que comprendió : cauce de ríos, estructuras hidráulicas (bocatomas, puentes) y las
planicies de inundación, también se realizó levantamiento batimétrico del cauce.

Figura 9. Resaltada cobertura LIDAR; Desembocadura rio Cañete – Localidad PAULLO.


Fuente: Área Topografía

Figura 10. Batimetría en Puente Clarita


Fuente :Área de topografía

Como producto de dicho levantamiento se generaron ortofotos, modelo digital de


elevación (DEM) y modelo digital del terreno(MDT).El levantamiento topográfico entregado se
generó las curvas de nivel y las redes irregulares de triángulos (TIN), para representar la
morfología de la superficie. Los modelos TIN de terreno cuentan con datos geográficos digitales
contienen vectores y se construyen mediante la triangulación de vértices(puntos), que luego es
llevado al ArcGIS, con la finalidad de convertirlo en formato ASCII.

36
Figura 11. Se observa las curvas de nivel correspondiente del tramo Desembocadura rio Cañete – Localidad Paullo
Fuente: Elaboración Propia.

3.2 Hidrología.

El área de hidrología determino caudales de máximas avenidas utilizando dos


metodologías:

La primera utilizando información propia de la estación de aforos existente en la cuenca;


la segunda a partir de las precipitaciones máximas en 24 horas de las estaciones comprendidas
en el área de estudio; para lo cual se ha aplico también un modelo de generación de descargas
máximas, a partir dichos, con el Software HEC HMS.

Se presenta un resultado comparativo de ambas metodologías de los valores de


Caudales de Máximas Avenidas obtenidas para el río Cañete, en la estación Socsi, para
diferentes períodos de retorno (10, 25, 50, 100, 200, 500 y 1000 años).

Tabla 11. Caudales Máximos Río Cañete - Estación Socsi (m3/s).

Modelo HEC- Análisis de


Periodo de Retorno (años)
HMS Frecuencias Máx.
10 384.2 625.3
25 613.8 790.6
50 827.7 923.2
100 1079.4 1064.1
200 1202.9 1214
500 1464.2 1427.7
1000 1790.5 1602.1

Fuente: Área Hidrología

37
4. METODOLOGÍA

4.1 Modelamiento Hidráulico

4.1.1 Ubicación de los tramos a Intervenir

Políticamente el tramo “Desembocadura– Puente Paullo” del río Cañete, que forma parte
del proyecto, se ubica en el departamento de Lima, provincias de Cañete. La longitud del tramo
es de alrededor de 31 kilómetros.

Ubicación Política:

• Departamento : Lima

• Provincia : Cañete

• Distrito : Nuevo Imperial, San Vicente de Cañete, y Lunahuaná.

• Localidades : Isla Baja, Herbay Alto, Pampa Castilla, Herbay Bajo, Paullo,
Lúcumo, Socsi, La Toma, Caltopa, El Paraíso y Caltopilla.

PPuente Aramayo

Bocatoma Palo Herbay

Figura 12. Vista satelital del tramo en estudio “Desembocadura–Bocatoma Palo Herbay”.
Fuente: Elaboración Propia.

38
Figura 13. Vista satelital del tramo en estudio “Bocatoma Palo Herbay-Localidad Paullo”.
Fuente: Elaboración Propia.

4.1.2 Topografía y cubierta del tramo de estudio

En la parte baja la topografía del Río Cañete en este tramo es de un valle amplio, el cual
presenta una llanura frágil a inundaciones en lado derecho aguas arriba de la Bocatoma La
Fortaleza. La longitud de este tramo del río Cañete es de aproximadamente 17 km, el cual
presenta una pendiente promedio de 1.1 %. El terreno está cubierto en su mayoría por cultivos
de surcos rectos; también, en las áreas adyacentes a las riberas se hallan bosques de bastante
maleza y cañaverales.

El río según la forma en planta es clasificado como trenzado. Además, este tramo de río
presenta caudal discontinuo, presentando máximas avenidas de ocurrencia en las estaciones de
lluvia en la serranía.

39
Bocatoma Fortaleza

Bosques

Figura 14. Vista satelital de la cobertura sobre la superficie - Bocatoma La Fortaleza.

Fuente: Elaboración Propia.

Figura 15. Bocatoma La Fortaleza.


Fuente: Elaboración Propia.

La topografía del Río Cañete en este tramo es de un valle amplio desde La Bocatoma
Nuevo Imperial hacia aguas abajo y de un valle estrecho hacia aguas arriba del mismo. La parte
amplia presenta una llanura frágil a inundaciones en lado derecho. Por otro lado, la parte
estrecha presenta menos riesgo por inundación, sin embargo, sufre fuerte erosión lateral de
bancos con cultivos sobre su superficie. La longitud de todo el tramo en estudio es de
aproximadamente 14.5 km, el cual presenta una pendiente promedio de 1.2 %. Las áreas
adyacentes a las riberas están cubiertas en su mayoría por cultivos de surcos rectos; también,
hay presencia de bosques de maleza densa en la parte amplia del valle (parte baja). De otro
lado, el río según la forma en planta es trenzado parte baja de la Bocatoma Imperial. En cambio,
de la parte de arriba ,de esta bocatoma el cauce deja de ser trenzado, siendo su cauce

40
encajonado en el lecho de un valle muy cerrado. Además, este tramo de río presenta un caudal
muy variable. En la temporada de lluvia (Diciembre – Abril) los caudales son significativamente
mayores debido a lluvias en la cuenca media y alta.

Bocatoma Imperial

Bosques Ripirianos

Figura 16. vista satelital de la cobertura - Bocatoma la Fortaleza.

4.1.3 Coeficiente del número de Manning (n)

Para la rugosidad se estableció un valor teniendo en cuenta que el material del lecho es
en su mayoría canto rodado. Por otro lado, las llanuras de inundación están cubiertas de cultivos
y bosques de densa maleza y cañaverales.

41
Figura 17. Granulometría del cauce y Bosque inundable de caña brava,parte baja Bocatoma La Fortaleza.

Figura 18. Bocatoma Imperial.

42
Figura 19. Bosques inundables y granulometría del cauce. Ambas fotos se tomaron parte baja de la Bocatoma
Imperial

La cobertura de la superficie fue sectorizada en tres partes tales como cultivos, bosques
y cauce propio del río. De manera que, a cada tipo de cobertura se le asignó valores “n” de
Manning de acuerdo a las tablas desarrolladas por Chow (1959).

Tabla 12. Coeficientes de Manning para cursos de agua naturales y llanuras de


inundación.

CORRIENTES NATURALES
D-1.Corrientes menores( ancho superficial en nivel creciente<100
pies) Mínimo Normal Máximo
a. Corrientes en planicies
1.Limpias,rectas,maximo nivel, sin montículos ni pozos profundos 0.025 0.030 0.033
2.Igual al anterior, pero con más piedras y malezas. 0.03 0.035 0.04
3.Limpio, serpenteante, algunos pozos y bancos de arena. 0.033 0.04 0.045
4. igual al anterior, pero con algunos matorrales y piedras 0.035 0.045 0.05
5.Igual al anterior, niveles bajo, pendientes y secciones más ineficientes 0.04 0.048 0.055
6.igual al 4, pero con más piedras. 0.045 0.05 0.06
D-2.Planicies de inundación
a. Pastizales, sin matorrales
1. Pasto corto 0.025 0.03 0.035
2. Pasto alto 0.03 0.035 0.05
b. Áreas cultivadas
1.Sin cultivos 0.02 0.03 0.04
2.Cultivos en línea maduros 0.025 0.035 0.045
3.Campos de Cultivos maduros 0.03 0.04 0.05
c. Matorrales
1.Matorrales dispersos, mucha maleza. 0.035 0.05 0.07
2.Pocos matorrales y árboles, en invierno. 0.035 0.05 0.06
3. Pocos matorrales y árboles , en verano. 0.04 0.06 0.05
4.Matorrales medios a densos , en invierno 0.045 0.07 0.11

Fuente: Ven Te Chow (1959).

43
Los diferentes tipos de cobertura fueron materializados de color, tal como se puede
observar en la figura 20(Tramo Desembocadura río Cañete - Bocatoma Palo Herbay) y figura
21(Tramo Bocatoma Palo Herbay – Localidad Paullo). Los valores de Manning asignados a cada
sector se los indica en las siguientes tablas:

Tabla 13. Asignación del coeficiente de Manning (n) al cauce y áreas adyacentes.

Color Descripción Manning


río 0.033
cultivo 0.04
bosque 0.07

Figura 20. Sectorización de la superficie del cauce y áreas inundables .

Tabla 14. Asignación del coeficiente de Manning (n) al cauce y áreas adyacentes.

Color Descripción Manning


cultivo 0.04
cauce 0.033
bosque 0.07

Figura 21. Sectorización de la superficie del cauce y áreas inundables .


Fuente: Elaboración Propia.

44
4.1.4 Hidrograma de diseño

El hidrograma empleado para las simulaciones corresponde a períodos de retorno (Tr) de


10, 100 y 500 años. De donde, el caudal de diseño (1064 m3/s) utilizado para las estructuras de
defensa pertenece a un Tr de 100 años, tal como lo indica los TDRs.

Hidrograma de Flujo Río Cañete


1600
1400
1200
Caudales (m3/s)

1000
800 Tr 10 años

600 Tr 100 Años

400 Tr 500 Años


200
0
0 20 40 60 80 100 120
Tiempo (hrs)

Figura 22. hidrograma de flujo de la estación Socsi .


Fuente: Elaboración Propia.

4.1.5 Criterios empleados en el trazo de eje de estructuras de defensa

En tiempo de avenidas se ocasiona desbordamiento del río, inundando sobre todo entre
la Bocatoma la Fortaleza- Bocatoma Imperial por el lado derecho (según vestigios y
declaraciones de los pobladores). En el lugar se ubican extensas áreas de cultivos, el cual son
afectadas por las inundaciones, aguas arriba de la bocatoma Imperial se manifiesta erosión en
la base de los bancos; debido que, cultivos y algunas viviendas se hallan sobre la superficie de
estos bancos. Este problema ocasiona pérdidas económicas a los agricultores, así como también
pone en riesgo la vida humana. Motivo por la cual, para controlar la vulnerabilidad de las zonas
en peligro, se definieron trazos de eje de defensas ribereñas como medida de protección.

El trazo de eje de las estructuras protección se realizó de acuerdo con criterios y


experiencias del proyectista. El procedimiento se detalla cronológicamente a continuación:

Paso N°1

• Modelamiento del río en condiciones actuales de terreno con el software IBER


2D, a fin de identificar áreas inundables y posteriormente realizar trazos de diques
en las áreas de inundación del río. El modelamiento numérico de inundaciones se
explicará detalladamente más adelante.

45
Paso N°2

• Análisis de la erosión lateral en las riberas por medio de una revisión de registros
históricos en imágenes satelitales del Google Earth Pro. Es decir, identificación
de tramos donde se ha desencadenado grandes migraciones lateras a causa de
las fuertes velocidades erosivas y del material frágil de los bancos. Cabe
mencionar que, para este tramo se dispuso de imágenes satelitales registradas a
partir del año 2002 a la actualidad.

Por ejemplo:

• Aguas arriba de la Bocatoma La Fortaleza en el año 2003 se observó una enorme


migración lateral en ambas riberas (figura 23). Incluso, las áreas del lado de la
ribera derecha, donde hoy se encuentran algunas áreas cultivadas, pertenecían
para ese tiempo a zonas activas por donde pasaba el flujo. Tales zonas
actualmente están cubiertas en su mayoría por maleza y cañaverales, y en menor
grado por cultivos. De otro lado, la existencia de extensas áreas cultivadas del
lado izquierdo de la Bocatoma en ese mismo año (2003), han desaparecido en la
actualidad (2020) a causa de desbordamientos y de la fuerte erosión lateral.

• Aguas abajo de la Bocatoma Imperial, a la altura del Centro Poblado El Paraíso,


en el año 2013 se observó una enorme migración lateral en la ribera derecha.
Tiempo atrás para el año 2009, Las áreas adyacentes a esta ribera, eran cubiertas
de considerables áreas cultivadas y bosque frondosos (figura superior 24). Sin
embargo, estas áreas vegetativas correspondían a zonas riparianas, es decir,
pertenecían a zonas activas por donde pasa el flujo en grandes eventos. Puesto
que, estas áreas cultivadas desaparecieron totalmente, a causa de
desbordamientos y de la fuerte erosión lateral, para el año 2013 (figura inferior
24) y no han sido invadidas por agricultores a la actualidad (2020).

Paso N°3

• Corroboración in situ de los sectores ribereños que se encuentran en peligro, tanto


por inundación como por erosión lateral. Para ello se recolectó información directa
por parte de los pobladores de la zona, rastros de obras hidráulicas colapsadas,
áreas inundadas y de marcas de agua que han dejado las últimas inundaciones.

Paso N°4

• A partir del análisis detallado en los pasos anteriores se trazaron los ejes de las
defensas, tomando como referencia los tramos críticos corroborados. Además,

46
para el trazado se tuvo en cuenta que las obras civiles y cultivos existentes no se
vean afectados por la obra nueva. Por este motivo se propuso un trazo
conservador que no seguía de manera exacta las metodologías encontradas en
los métodos propuestos en la bibliografía revisada. Puesto que, estas teorías
elaboradas en otros países, en laboratorio o para ríos de características distintas
no cumplen en su totalidad con la dinámica del flujo que desarrollan los ríos
peruanos, no obstante, sirven de soporte para la elaboración del presente diseño
de defensas. Los ríos de montaña no han sido tan estudiados como los ríos de
baja pendiente, en los que se ha recolectado una gran cantidad de datos.

A. Tramo: Desembocadura río cañete – Bocatoma Palo Herbay

Figura 23. Trazo inicial propuesto tomando como referencia la Fluviomorfología del año 2003.
Fuente: Elaboración Propia.

Tomando como ejemplo el mismo tramo mostrado en la figura 23 el ancho entre trazos
debe ser de prioridad mayor al ancho estable calculado teóricamente. La razón es porque se
hallan actualmente vestigios de diques colapsados en casi todo el tramo en estudio; pues, el
colapso muy probablemente se debió a causa de que el ancho entre diques existentes
corresponde a un rango teórico. Esto indica que, como proyectista se debe aprender de los
errores cometidos por nuestros antecesores y no tanto afianzarse a cálculos teóricos. Por lo que,
de los diversos anchos estables teóricos se elegiría el mayor, con el fin aportar más seguridad
al diseño de estas nuevas obras proyectadas.

Aclarado lo anterior, se procedió a estimar anchos teóricos que se muestran en la tabla


17, a fin de elegir el ancho que se ajuste más a la del tramo de río en estudio. De acuerdo con
esta tabla se decidió elegir el ancho teórico de 205 m, estimado por el método de Blench y
Altunin, ya que se acomoda mejor al lecho del tramo en estudio. Cabe indicar que, al elegir este
ancho de más magnitud no quiere decir que se cumplió estrictamente con este cálculo teórico,

47
sino que en lo posible se realizaron trazos tomando como referencia el mínimo de 205 m.

Tabla 15. Parámetros hidráulicos del Tramo en estudio.

Parámetros hidráulicos
Tirante medio: h= 1.32 m
Diámetro 50: dm = 73.3 mm
Caudal diseño: Q= 1064 m3/s
Pendiente: S= 0.011
Manning: n= 0.033

Tabla 16. Métodos de cálculo para el ancho estable (Be) del tramo en estudio

METODOS Be (m) MATERIAL

Julien y Wargadalam 107 Granular

Altunin 106 Grava a más

Pettis 145 ----

Simons y Henderson 95 Grava

Blench y Altunin 205 Gruesos

Manning y Strickler 120 Aluvial grueso

Es importante mencionar que, los diques colapsados fueron construidos con un ancho
teórico muy estrecho de alrededor de 131 m. Según parece, estos diques existentes no han sido
capaz de resistir los grandes flujos desencadenados en los últimos años, sobre todo a los flujos
del Niño Costero del 2017. Tal como se observa en la figura 24, aguas arriba de La Bocatoma
La Fortaleza se puede apreciar desborde en varios tramos cortos y consecutivos del dique del
lado derecho.

Considerando los criterios técnicos anteriormente detallados se procedieron a realizar


trazos propios preliminares en los sectores en riesgo. Parte de estos trazos se pueden observar
mediante líneas de color magenta en la figura 23 y 24 correspondientes a los años 2003 y 2019,
respectivamente.

48
Figura 24. Vista de trazos de diques colapsados en avenidas del año 2019.
Fuente: Elaboración Propia.

Luego, por motivos que en este tramo de estudio se contaba con un proyecto avanzado
de protección parcial de riberas (Valles Vulnerables), hubo la necesidad de acoplarse a los trazos
de este proyecto. Adicionalmente, en la figura 25 se muestra un corto tramo aguas arriba de la
Bocatoma La Fortaleza, lugar donde se estrechó el trazo propio, con el propósito de unir y formar
ambos trazos de diques en uno solo. Es así que, bajo estas consideraciones se elaboraron trazos
definitivos en todos los sectores expuestos a graves riesgos en el tramo de estudio
“Desembocadura río Cañete - Bocatoma Palo Herbay”.

Figura 25. Trazos definitivos de ejes de diques para el año 2020-Bocatoma La Fortaleza.

B. Tramo: Bocatoma Palo Herbay-Puente Paullo

Tomando como ejemplo el mismo tramo mostrado en la figura 26: el ancho entre trazos
fue de prioridad mayor al ancho estable calculado teóricamente (Be). La razón fue porque se
hallan actualmente diques colapsados; pues, el colapso muy probablemente se debió a causa

49
de que el ancho entre diques existentes corresponde a un rango teórico. Esto indica que, como
proyectista se debe aprender de los errores cometidos en la zona del estudio y no seguir
exactamente cálculos teóricos que se derivan de pruebas de laboratorio, así como de ríos que
tienen otros parámetros. Por lo que, de los diversos anchos estables teóricos se eligió el mayor,
con el fin aportar más seguridad al diseño de estas nuevas obras proyectadas.

Sin embargo, por motivos que a partir de la Bocatoma Imperial hacia aguas arriba ancho
del cauce es encajonado, se tomó un ancho mínimo teórico calculado de Be = 95 m debido a
que en esta parte del tramo en estudio el río deja de ser un rio trenzado de llanura, y tiende a
tomar un comportamiento de rio de montaña de cauce estrecho.

Figura 26. Trazos en la Fluviomorfología del año 2009 (figura superior) y 2013 (figura inferior).
Fuente: Elaboración Propia.

50
Aclarado lo anterior, se procedió a estimar anchos teóricos que se muestran en la tabla
19, a fin de elegir los anchos que se ajusten más al amplio cauce aguas abajo de la Bocatoma
Imperial y al estrecho cauce aguar arriba de la misma. De acuerdo con esta tabla se decidió
elegir el ancho teórico máximo de Be = 205 m, estimado por el método de Blench y Altunin, para
la parte del amplio cauce trenzado; así mismo, para la parte del estrecho cauce se decidió elegir
el ancho teórico mínimo de Be = 95 m, estimado por el método de Simons y Henderson. Cabe
indicar que, al elegir los anchos de mayor y menor magnitud no quiere decir que se cumplió
estrictamente con este cálculo teórico, sino que en lo posible se realizaron trazos tomando como
referencia estos dos anchos teóricos.

Tabla 17. Parámetros hidráulicos del tramo en estudio.

Parámetros hidráulicos
Tirante medio: h= 1.44 m
Diámetro D50: dm = 73.3 mm
Caudal diseño: Q= 1064 m3/s
Pendiente: S= 0.012
Manning: n= 0.033

Tabla 18. Métodos de cálculo para el ancho estable (Be) del tramo en estudio

METODOS B (m) MATERIAL SUBTRAMO

Julien y Wargadalam 112 Granular ----

Altunin 118 Grava a más ----

Pettis 145 ---- ----

Aguas arriba
Simons y Henderson 95 Grava
Imperial
Aguas Abajo
Blench y Altunin 205 Gruesos
Imperial
Aluvial
Manning y Strickler 118 ----
grueso

En la figura 27 se observan una serie de espigones han sido construidos sobre la ribera
izquierda. Durante la visita de campo se corroboró que los espigones muestran serias
deficiencias que han conducido a su colapso. Se pudo notar que, en un proyecto anterior, se
tomó la decisión de solucionar problemas previos colocando diques delante de las puntas de los
espigones, lo que no ha funcionado en la práctica. El actual proyecto recomienda que el
diseñador debe adaptarse a la naturaleza del flujo y no al contrario. Es decir, en esta parte del

51
tramo en estudio el río presenta naturalmente un cauce muy ancho, advirtiendo que no necesita
ser forzado a ser estrechado. Dado que, los diques colapsados fueron construidos con un ancho
teórico muy estrecho de alrededor de 102 m, los restos de diques y espigones indican que no
han sido capaz de resistir la dinámica de los fuertes flujos en eventos máximos, como los del
niño 2017. Según rastros que han dejado desbordes de flujo, aguas abajo de La Bocatoma
Imperial se nota que el flujo ha destruido ciertos tramos cortos y consecutivos de los diques de
ambos lados.

Figura 27. Vista de diques colapsados y transición natural del ancho teórico del cauce.
Fuente: Elaboración Propia.

Considerando los criterios técnicos anteriormente detallados se procedió a realizar trazos


propios preliminares en los sectores en riesgo. Posteriormente, debido a que en este tramo de
estudio se contaba con un proyecto avanzado de protección parcial de riberas (Valles
Vulnerables), hubo la necesidad de reducir el ancho del cauce y acoplar los trazos propios a los
de este proyecto (ver líneas de color verde, figura 28). Otro motivo por el cual el ancho del cauce

52
se redujo fue porque aguas arriba de la Bocatoma Imperial el cauce empieza a estrecharse
naturalmente hasta menos de la mitad del cauce aguas abajo. En las imágenes 27 y 28 se
muestra un corto tramo aguas abajo de la Bocatoma Imperial, lugar donde ocurre la fuerte
transición del ancho del cauce; motivo por el cual, el ancho teórico también fue cambiado
radicalmente desde 205 m a 95 m. Razón por la cual, el trazo propio es estrechado
considerablemente con el propósito de adaptarse a la transición natural del cauce además de
unir y formar ambos trazos de diques en uno solo.

Figura 28. Unión de trazos definitivos de ejes de diques de ambos proyectos en la zona de transición natural del
cauce (Bocatoma Imperial).
Fuente: Elaboración Propia.

53
Paso N°5

Modelación en IBER 2D del río con los trazos de eje definitivos de diques y con el caudal de
diseño de 1064 m3/s, correspondiente a un Tr de 100 años. Para ello, los diques preliminares
ensamblados en la superficie de terreno fueron dimensionados provisionalmente en base a
la experiencia del proyectista. Pues, estos trazos y dimensiones sufrirían luego pequeñas
modificaciones hasta lograr cumplir con el objetivo del proyecto, es decir la protección.

Cabe indicar que, la Simulación IBER 2D fue de manera consecutiva e iterativa hasta lograr
que el flujo no genere desbordes a las áreas de protección, así mismo se recolecte una
confiable data hidráulica. Dado que de esta data (calados, áreas hidráulicas y espejos de
agua) depende un buen diseño definitivo del enrocado de los diques.

4.1.6 Modelo Digital de Elevación (DEM)

El DEM utilizado para el tramo: Desembocadura río Cañete – Bocatoma Palo Herbay fue
en formato tiff, cuyos pixeles fueron de 0.5m x 0.5m. Sin antes indicar que, la resolución del vuelo
Lidar fue 0.3 m x 0.3 m, no obstante, se disminuyó a una dimensión prudente de 0.5 m x 0.5 m
por motivos de capacidad de rendimiento de procesador.

Figura 29. Modelo digital de elevación en formato tiff de 0.5 m x 0.5 m de pixel.
Fuente: Elaboración Propia.

54
El DEM utilizado para el tramo: Bocatoma Palo Herbay – Puente Paullo fue en formato
tiff, cuyos pixeles fueron de 0.4 m x 0.4m. Sin antes indicar que, la resolución del vuelo Lidar fue
0.3 m x 0.3 m, no obstante, se disminuyó a una dimensión prudente de 0.4 m x 0.4 m por motivos
de capacidad de rendimiento de procesador.

Figura 30. Modelo digital de elevación en formato tiff de 0.4 m x 0.4 m de pixel.
Fuente: Elaboración Propia.

4.1.7 Modelamiento IBER 2D

A. Geometría del modelo 2D – editar Malla

Para empezar a realizar los cálculos IBER se necesitó delimitar un área de trabajo por
medio de un formato dxf. Luego, al área delimitada se lo asigno un determinado tamaño de malla
no estructurada. A continuación, a esta área se le asigno una elevación mediante una geometría
en formato ASCII, importada desde el programa ArcMap 10.5.

Es importantes mencionar que, durante el proceso continuo de simulaciones IBER 2D se


fueron probando diferentes tamaños de espaciamiento lateral de la malla triangular. Este
procedimiento fue con el fin de lograr que los calados del modelamiento se asemejen a las
marcas de agua halladas en las visitas en campo. De igual modo, el modelamiento IBER 2D tuvo
que ser consecutivo hasta lograr que sus resultados de desborde e inundación concuerde con
la información recibida por parte de los pobladores, rastros hallados en campo e imágenes
satelitales.

55
Tramo: Desembocadura río Cañete – Bocatoma Palo Herbay

El proceso de mallado definitivo fue de mallas estructuradas de 5 m de espaciamiento


lateral para tota la superficie. Previamente, se probó con mallas de 25, 10, 5 y 3. De los cuales,
la malla de 0.5 m fue suficiente para permitir que el modelo se adecue apropiadamente a la
superficie batimétrica e inundaciones, incluyendo todas las obstrucciones y demás detalles
presentes.

Figura 31. Generación de la malla estructurada de 5 m


Fuente: Elaboración Propia.

56
Tramo: Bocatoma Palo Herbay – Puente Paullo

El proceso de mallado definitivo fue de mallas estructuradas de 4 m de espaciamiento


lateral para tota la superficie; sin antes indicar que, previamente, se probó con mallas de 10 y 5.
De los cuales, la malla de 0.4 m fue suficiente para permitir que el modelo se adecue
apropiadamente a la superficie batimétrica e inundaciones, incluyendo todas las obstrucciones y
demás detalles presentes.

Figura 32. Generación de la malla estructurada de 4 m de espaciamiento lateral.


Fuente: Elaboración Propia.

B. Condiciones de contorneo

Tramo: Desembocadura río Cañete – Bocatoma Palo Herbay

La condición inicial del modelo es el hidrograma de entrada para cada uno de los
escenarios, en donde los caudales llegan a un máximo de 625 𝑚3/𝑠, para un Tr de 10 años; 1064
𝑚3/𝑠, Tr de 100 años; y 1428 𝑚3/𝑠, Tr de 500 años. A fin de atenuar el ingreso brusco del flujo,
el punto de inicio de descarga se ubicó a 2.1 k𝑚 aproximadamente aguas guas arriba de la
Bocatoma Palo Herbay. El punto de condición de salida se ubicó lo más alejado posible, es decir
hasta el litoral marino.

57
Para el caso del modelamiento de inundaciones se emplearon los caudales del
hidrograma correspondiente a períodos de retorno (Tr) de 10, 100 y 500 años. Para el caso del
modelamiento con estructuras de defensa (diques de enrocado) se emplearon solo los
correspondientes al Tr de 100 años, tal como lo indica los TDRs.

Condición
entrada

Condición
salida

Figura 33. Vista de condición de entrada y salida en el MDT.


Fuente: Elaboración Propia.

Tramo: Bocatoma Palo Herbay – Puente Paullo

La condición inicial del modelo es el hidrograma de entrada para cada uno de los
escenarios, en donde los caudales llegan a un máximo de 625 m3/s, para un Tr de 10 años;
1064 m3/s, Tr de 100 años; y 1428 m3/s, Tr de 500 años. A fin de atenuar el ingreso brusco del
flujo, el punto de inicio de descarga se ubicó a 0.57 km aproximadamente aguas arriba del Puente
Paullo. El punto de condición de salida se ubicó aproximadamente a unos 3.5 km de la parte final
del tramo (Bocatoma Palo Herbay).

Para el caso del modelamiento de inundaciones se emplearon los caudales del


hidrograma correspondiente a períodos de retorno (Tr) de 10, 100 y 500 años. Para el caso del
modelamiento con estructuras de defensa (diques de enrocado) se emplearon solo los caudales
correspondientes al Tr de 100 años, tal como lo indica los TDRs.

58
Figura 34. Vista de condición de entrada y salida en el MDT.
Fuente: Elaboración Propia.

C. Número de Courant e Incremento de Tiempo

Es el valor mínimo que necesita el programa para iterar , el incremento de tiempo máximo
(Δt) es relacionado al Número de Courant (C o CFL), obtenido de la relación entre espacio (Δx),
velocidad y tiempo. Cabe precisar que dicha relación debe ser menor igual a la unidad.

C : Número de Courant
V : Velocidad
Δx : Tamaño de malla o celda
Δt : Incremento de tiempo máximo.

Tramo: Desembocadura río Cañete – Bocatoma Palo Herbay

Luego de varias simulaciones se optó por un tamaño definitivo de malla de 5 m, de donde


se obtuvo una velocidad media de 3.5 m/s, por lo que reemplazando estos valores se obtuvo un
incremento de tiempo máximo (Δt) para el cálculo computacional como:

59
C = 0.8
V = 3.5 m/s
∆X = 5 m
∆T = 1.14 s

Tabla 19. Ventana de Datos de procesamiento en Iber 2D

Tramo: Bocatoma Palo Herbay – Puente Paullo

Luego de varias simulaciones se optó por un tamaño definitivo de malla de 4 m, de donde


se obtuvo una velocidad media de 4.2 m/s, por lo que reemplazando estos valores se obtuvo un
incremento de tiempo máximo (Δt) para el cálculo computacional como:

C = 0.8
V = 4.2 m/s
∆X = 4 m
∆T = 0.76 s

60
Tabla 20. Ventana de Datos de procesamiento en Iber 2D

4.2 Diseño Hidráulico

4.2.1 Criterios de diseño de los diques de enrocado.

El diseño de las defensas ribereñas consistió en diques de enrocados. Las


consideraciones específicas para adoptar esta alternativa fueron por la disponibilidad del
material del lecho del río para la conformación del dique, cercanía de canteras para el enrocado;
además, por ser estructuras que no disminuyen el área hidráulica, cuyo factor es elemental sobre
todo para los tramos estrechos.

Para su diseño se consideraron los siguientes criterios:

4.2.1.1 Trazado en perfil .

El nivel de la corona de los diques será establecido de acuerdo al borde Libre calculado
para un tiempo de retorno de 100 de un caudal de 1064 m3/s

Borde libre =1.20 m

Los planos de planta y perfil longitudinal presentan los diferentes niveles de la corona de
dique a lo largo de toda su longitud. Además, el nivel de cimentación de la uña anti socavante
tal como lo denomina los TDR a la protección del pie del dique.

61
4.2.1.2 Profundidad de socavación.

Para determinar la profundidad de socavación en el fondo del lecho, se usará la fórmula de


Lischtvan – Lebediev, se necesitó de parámetros hidráulicos obtenidos de campo y de
modelación definitiva IBER 2D. Esta altura de socavación representa la altura de la uña o
profundidad de entierro del enrocado.

Como es normal la socavación natural en un río varía a lo largo de su recorrido. Sin embargo,
para el diseño de un tramo de enrocado se decidió tomar las socavaciones más profundas, con
el propósito de que la uña aporte mayor seguridad al dique. Habría que mencionar también que,
los largos trazos de eje de dique tuvieron que ser parcializados en tramos menores, cuyo
recorrido presente características hidráulicas similares, principalmente, de calado y espejo de
agua. Puesto que, el diseño del enrocado depende principalmente de estos parámetros.
Posteriormente, el nivel de cimentación del enrocado se obtuvo restando la socavación a la cota
más inferior de cada sección.

4.2.1.3 Altura del enrocado.

El cálculo de la altura de la estructura se ha considerado la uña, calado y borde libre. Mencionar


que, el calado con el que se diseñó fue tomado de la sección que presenta mayor profundidad
en los tramos a intervenir , según la simulación IBER 2D. El criterio fue con el objeto de darle
mayor seguridad a la estructura hidráulica.

4.2.1.4 Diámetro medio de la roca (D50) .

Se consideró un diámetro, cuyo valor está dentro de un rango de 0.6 m a 1.2 m. Los métodos
que se emplearon en la estimación del diámetro fueron las fórmulas de Maynord e Isbash. Cabe
mencionar que, en algunos tramos muy anchos se estimaron D50 de roca mínimos de hasta 0.4
m. Así mismo, en algunos tramos muy estrechos entre diques se estimaron D50 de roca máximos
de hasta 1.8 m. Sin embargo, estos diámetros excesivos y mínimos, por experiencia del
proyectista, fueron ajustados dentro del rango mencionado (0.6-1.2).

4.2.1.5 Dimensiones transversales del dique de enrocado.

La dimensión está en función de los parámetros hidráulicos y del diámetro medio de la roca. El
ancho mínimo del enrocado fue el doble del diámetro medio de la roca (T=2 D50). La corona del
dique en sí, el cual no incluye el ancho de enrocado, lleva un ancho de 4 m. El talud de la cara
húmeda del enrocado fue de 1:2 y ambas caras del dique fue de 1:1.5. La sección típica se
muestra en la figura 35.

62
Figura 35. Seccion Tipica
Fuente: Elaboración Propia.

Luego, una vez ya definido los trazos y diseño de diques ensamblados sobre el terreno del tramo
en estudio se lleva a cabo la última modelación hidráulica IBER 2D.

63
5. ANÁLISIS Y DISCUSIÓN DE RESULTADOS.

5.1 Modelamiento por inundación (sin estructuras de defensa)

Los modelamientos hidráulicos por inundación se elaboraron para periodos de retorno


(Tr) de 10, 100 y 500 años, el cual se muestran gráficamente en el anexo mapas de
inundaciones. de estos, los resultados de las inundaciones ocurridas para el Tr de 100 años, se
analizó y discutió al detalle por corresponder al Tr de diseño de las estructuras de defensa. cabe
recalcar que, el caudal pico correspondiente al Tr de diseñó es de 1064 m3/s.

5.1.1 Análisis para el Tr de 100 años

Tramo: Desembocadura río Cañete – Bocatoma Palo Herbay

Primero antes de empezar, para un mejor entendimiento de los resultados se muestra el


mapa de inundaciones completo del tramo “Desembocadura río Cañete - Bocatoma Palo Herbay”
en la figura 36. En la figura se puede observar que ocurre bastante inundación aguas arriba de
la Bocatoma La Fortaleza, cuya ubicación es casi al centro de este tramo. De otro lado, el
desborde aguas abajo de esta bocatoma parece ser poco o casi nada según la figura, pero si se
observa minuciosamente se nota que ocurrieron leves inundaciones. Ahora bien, el motivo
principal por el que no ocurrieron inundaciones fue porque se han construido diques rústicos en
años recientes, según registros de imágenes satelitales. Estos diques han sido levantados con
material propio del río sin ningún tipo de planificación ni diseño, en consecuencia, serían a la
larga averiados con facilidad en ocurrencia de máximas avenidas. Por lo que, las áreas de cultivo
y centros poblados adyacentes a las riberas serían gravemente afectados.

64
Figura 36. Mapa de inundaciones en IBER 2D del tramo “Desembocadura río Cañete - Bocatoma Palo Herbay” para el Tr de 100 años.

Fuente: Elaboración Propia.

65
Ahora veamos, aguas abajo de la Bocatoma Palo Herbay se aprecia el desborde de flujo
a través de la rotura en un trecho de un dique existente (figura 37). El flujo que desborda afecta
áreas de cultivos adyacentes del lado de la ribera izquierda; por supuesto que, estos daños
fueron corroborados por agricultores.

Figura 37. Vistas satelitales de zonas en riesgo -Bocatoma La Pinta.

66
Unos metros más abajo de la rotura se observa unos menores desbordes que igualmente
afectan los cultivos. Sin embargo, estos desbordes de menor magnitud transitan sobre la cresta
del dique existente, el cual cuenta con una longitud aproximada de 2.5 km a partir de la rotura
indicada. Cabe aclarar que, este dique ha cumplido su función de defensa hasta la actualidad;
por lo que, se decidió mantenerlo esta estructura sin antes recomendar su mantenimiento y
elevación de cresta de 1 m.

A continuación, al final del dique existente se inicia una estructura tipo dique rústico hasta
la Bocatoma La Fortaleza. Al parecer, tal dique rustico anteriormente se hallaba reemplazado
por un dique diseñado, pero que al haber sido colapsado fue remplazado por este dique rústico
que aún permanece a la actualidad. Si bien los diques rústicos aun cumplen con proteger áreas
adyacentes, esto no sería por mucho tiempo. La razón es porque: los diques han sido construidos
rústicamente con piedra amontonada del cauce en la lateral derecha. Lo peor del caso, en la
visita de campo, antes de llegar a la Bocatoma Fortaleza, se observó que gran longitud de estos
diques han colapsado en su totalidad. Lo mencionado se puede corroborar en la figura 38:

Figura 38. Vistas satelitales de zonas en riesgo -Bocatoma La Pinta.

En el otro lado (orilla derecha, figura 38) se observa que ocurre bastante inundación a
partir de la bocatoma la fortaleza hacia arriba. La longitud de ribera por donde desborda el flujo
es de aproximadamente 3.5 km. El área inundada adyacente a esta ribera pertenece a una zona
activa, según se muestra en imágenes satelitales hasta el año 2003. En dicha zona activa o
mejor conocida como zona ripariana (riparian zone, en inglés) ha ocurrido una migración lateral
de hasta casi 300 m; sin embargo, dentro de esta zona en peligro se hallan considerables áreas
de cultivo expuestas a inundación.

En lo que respecta al tramo restante (Bocatoma La Fortaleza – Desembocadura río


Cañete), el desborde e inundación bajó considerablemente. El motivo principal por el que no
ocurrieron inundaciones fue porque se han construido grandes longitudes de diques rústicos,

67
que han sido también capturados al momento de realizar el levantamiento topográfico con vuelo
LIDAR. Sin embargo, al revisar registros históricos de imágenes satelitales se notó que los
diques han sido destruidos con facilidad después de cada avenida extraordinaria.
Posteriormente, estos diques rústicos con el afán de protección han sido reconstruidos
consecutivamente a lo largo de los años, pero no han funcionado como se esperaba.

Para ilustrar mejor, se enfoca y hace un análisis en la simulación por inundación IBER
2D, en un tramo corto aguas abajo del Puente Clarita: las viviendas ubicadas en las
inmediaciones del puente están sujetas a graves peligros por erosión lateral más que por
inundación; la razón es porque, las viviendas están asentadas sobre un área elevada, cuyo
material del banco ribereño es muy frágil para hacer frente a las fuertes velocidades del flujo en
crecidas (figura 42).

Figura 39. Puente Clarita.

Figura 40. Viviendas cerca de la ribera derecha aguas arriba del Puente Clarita.

68
Figura 41. Viviendas cerca en la ribera Izquierda aguas abajo del Puente Clarita.

A continuación de las viviendas se halla una inmensa área o llanura de inundación, pero
que no ha sido inundado por las simulaciones IBER 2D; debido que, el dique rustido existente
ha obstaculizado el paso del flujo a la llanura. Del mismo modo, la llanura de inundación ubicada
del otro lado (ribera izquierda) no fue afectada de acuerdo con las simulaciones, a causa también
de la existencia de otra obra rústica. No obstante, ocurre un desborde mínimo por la entrada de
un dren colector, que vierte al río las aguas provenientes de los cultivos. Lo mencionado se
puede observar en la figura siguiente:

69
Figura 42. Viviendas cerca de la ribera Izquierda aguas abaao del Puente Clarita.
Fuente: Elaboración Propia.

70
Con respecto a, las áreas que no fueron inundadas en la simulación IBER 2D se debió
corregir mediante la desaparición de estos muros, que representan diques rústicos. Trabajo que
debió realizarse en la superficie topográfica levantada por el vuelo LIDAR. Ahora bien, tal trabajo
no fue necesario ya que el desborde y recubrimiento por parte del agua en las llanuras de
inundación fueron observadas en imágenes satelitales de una serie de años anteriores.
Específicamente, las áreas en peligro por la presencia de inundaciones fueron detectados
gracias a imágenes satelitales disponibles desde el año 2003. De igual modo, la migración lateral
de bancos ribereños de material frágil fue identificada gracias a las imágenes satelitales.

Dicho lo anterior, para una mejor comprensión se recurrió a un análisis grafico del mismo
tramo de la simulación IBER 2D en la figura 43; esto es, el mismo tramo es aguas abajo del
Puente Clarita.

En la figura 43A correspondiente al año 2009 se aprecia un trazo de un dique de la ribera


derecha, que ha sido en la parte media totalmente destruida por el flujo. Aguas arriba y detrás
de este rastro de dique se nota la llanura de inundación totalmente descubierta, probablemente
hubo presencia de una avenida de algún año reciente. Razón por la cual, autoridades u otra
entidad se vio en la necesidad de construir el dique rustico que ha sido captado por la figura
satelital de ese año. Además, en el centro poblado cerca al Puente Clarita se observa que se
encuentra totalmente desprotegida, expuesto a riesgo por erosión. De otro lado, en la ribera
izquierda se halla la llanura de inundación ocupada por un área ripariana cubierta en su totalidad
por maleza.

Luego, en la figura 43B correspondiente al año 2017 ya no se puede visualizar ningún


tipo de estructura de defensa (dique). Los diques que fueron construidos artesanalmente en años
anteriores han desaparecido en su totalidad, ocasionando grandes áreas inundadas .Así mismo,
el área ripariana ubicada en la llanura inundable izquierda ha desaparecido en su totalidad
producto del desborde. Este fenómeno hidráulico se explica y guarda coherencia con lo visto en
la figura; puesto que, en ese año (2017) ocurrió el fenómeno del niño que acarreo graves daños
a lo largo de la costa peruana.

Ya para el año 2019, según la figura 43C, se observa que toda la ribera derecha ha sido
nuevamente protegida por diques. Lo mismo ocurrió con la llanura inundada de la ribera
izquierda, que ha sido protegida a partir de un punto donde la ribera no ha sufrido migración
lateral hacia el puente Aramayo.

El análisis explicado hasta este punto ha sido empleado para todo el tramo de estudio;
motivo por el cual, a los diques existentes no se los limpió de la superficie topográfica porque la
magnitud de inundaciones fue registrada en las imágenes satelitales.

71
72
Figura 43. Vistas satelitales y desde el terreno mismo del colapso continuo de los diques rústicos, aguas abajo del Puente Clarita.
Fuente: Elaboración Propia.

73
Tramo: Bocatoma Palo Herbay – Puente Paullo

Los resultados se muestra el mapa de inundaciones IBER 2D completo del tramo


“Bocatoma Palo Herbay - Puente Paullo” en la figura 44,se puede observar que ocurre bastante
inundación, sobre todo del lado derecho, aguas abajo del Puente Socsi. Justo en este lado
derecho se encuentran la mayor cantidad de áreas de cultivos, la Bocatoma Imperial y la
Carretera Cañete – Yauyos. Lo peor del caso, la carretera llega hacer afectada por inundación a
la altura donde se ubica la bocatoma. Por otra parte, el desborde e inundación parece ser de
poca magnitud aguas arriba del Puente Socsi, pero si se observa minuciosamente se nota que
si hay ocurrencia de menores inundaciones. El motivo principal por el cual las inundaciones son
bajas es porque el cauce empieza a ser más encajonado naturalmente, tendiendo a comportarse
como un río montañoso.

Adicionalmente, se debe tomar en cuenta que, desde la Bocatoma Imperial hacia aguas
arriba se encuentran numerosos espigones de gaviones y diques de enrocados diseñados como
defensas de protección. Los diques se encuentran en su mayoría cumpliendo correctamente su
función en contra la socavación, sin embargo, para que pueda cumplir su función contra
inundaciones las cotas de sus coronas deben ser incrementadas. Los bordes libres de estos
diques existentes no son suficientes para evitar inundaciones ocasionadas por un caudal de
1064 m3/s (Tr=100 años). Además, los espigones hechos de gaviones han sido inundados
totalmente por este mismo caudal. Además, la mayoría de estos espigones han sido afectados
estructuralmente en forma grave, sobre todo en las bases de las puntas, según se pudo apreciar
en las visitas de campo. Por lo indicado anteriormente, el colapso de estos espigones es
inminente.

74
Figura 44. Mapa de inundaciones en IBER 2D del tramo “Puente Paullo - Bocatoma Palo Herbay” para el Tr de 100
años.
Fuente: Elaboración Propia.

Aparte de las simulaciones IBER 2D, se recurrieron a registros históricos de imágenes


satelitales, a fin de identificar y analizar el riesgo de los sectores vulnerables. Estos sectores
fueron identificados con facilidad debido a los graves daños que han ocasionado las avenidas
máximas, por medio de inundaciones y fuertes velocidades erosivas en las bases de bancos de
material frágil.

Para enfatizar en el análisis se detalla uno de los puntos más críticos, el cual se ubica a
la altura del Centro Poblado Lúcumo, aguas abajo del Puente Paullo. En este sector, no solo las
áreas de cultivos y estructuras hidráulicas están expuestas a riesgos en máximas avenidas, sino
también se halla en riesgo la vida humana (figura 45).

En la figura 45A, correspondiente el año 2015, se observa que, en el afán de ganar


terreno, algunos pobladores han construido espigones sin ningún tipo de diseño ni planificación.
Se aprecia que el ancho del cauce del río ha sido reducido desde ambos lados de una manera

75
muy exagerada, llegando a medir 44 metros sin incluir la longitud del espigón. En todo caso, si
se restase la longitud de esta estructura el ancho libre para el paso flujo llega exageradamente
a medir alrededor de 25 metros. Esto conlleva a que en este punto se tenga un área hidráulica
exageradamente reducida para el tránsito de avenidas de gran magnitud y no podría soportar el
caudal de diseño de 1064 m3/s.

A continuación, según la figura 45B: el enorme flujo para el año 2017, fecha en que
ocurrió el ultimo niño costero, devastó con todo lo que irrumpió su paso en ambas riberas. Es
decir, la enorme avenida ocasionó la desaparición de los espigones, una vivienda, un canal de
regadío y áreas de cultivo, llegando el cauce a alcanzar un ancho de hasta 124 metros. La
inversión en espigones hechos probablemente de forma rústica no cumplió su cometido, error
del cual debió servir como base para no repetir este tipo de diseño. Sin embargo, a pesar de
haber experimentado la capacidad de destrucción que puede causar una máxima avenida se
volvió a cometer la misma irresponsabilidad, incluso llegando a poner nuevamente en riesgo la
vida humana. Se ha construido una vivienda poco después de la ocurrencia del Niño Costero y
que actualmente mantiene la ocupación.

En el año 2019, como se observa en la figura 45C, se reconstruyó el canal de regadío


que fue destruido en el fenómeno del niño del 2017. Para proteger este canal y cultivos que se
ubican en el lado cóncavo de la curvatura del río, se han construido espigones de gaviones a lo
largo de la ribera. Esta medida ingenieril fue necesaria y correcta, ya que además se protege las
enormes cárcavas que contienen cultivos y viviendas en su cima. Como se sabe, en las bases
de los bancos del lado cóncavo en una curvatura es donde se desarrollan fuertes erosiones
generadas por las líneas de corriente del flujo helicoidal principal. Así que, el sistema de
espigones diseñado de gaviones ha cumplido su función correctamente a actualidad, no
obstante, está presentando deficiencias sobre todo en las bases de las puntas. A su vez, en el
lado convexo (izquierdo o interior) de la curvatura, sin ningún tipo de diseño ni planificación, se
ha construido un dique compacto de roca para proteger la vivienda ubicada dentro del cauce.
Con este dique el ancho del cauce ha sido reducido nuevamente de manera considerable.

76
77
Figura 45. A, B y C vistas satelitales correspondientes a los años 2015, 2017 y 2019 respectivamente de la zona en
riesgo sector Lúcumo. D fotografía de la situación actual de las riberas en riesgo en el sector Lúcumo.
Fuente: Elaboración Propia.

78
La vista panorámica en la figura 45D confirma lo antes detallado. Así que, por el lado
cóncavo (derecho) es necesario reforzar la protección con diques; y por el lado convexo
(izquierdo) se debe retirar el dique existente y ensanchar los más posible el cauce antes de
proyectar una nueva defensa.

El grado de peligrosidad del sector mostrado en la figura anterior queda reconfirmado


mediante la simulación, tal como lo muestra la figura 46, se muestra que para un caudal de
diseño de 1064 m3/s: el flujo rebosa sobre el dique existente, llegando a alcanzar gran parte de
la vivienda y cultivos. Así mismo, el nivel del flujo en el otro lado cubre totalmente a los espigones,
alcanzando y afectando el canal de regadío y cultivos adyacentes.

Figura 46. Vista de la inundación según la simulación IBER 2D en el sector Lúcumo.

En las inmediaciones del Puente Socsi ,la sección del cauce no soporta el paso de un caudal de
diseño de 1064 m3/s. O sea, el excesivo caudal causa desbordes por ambas riberas a lo largo
de un considerable recorrido, tal como se observa en la figura 47. De acuerdo con esta figura
con simulación IBER 2D, el ancho del cauce en el eje del puente contiene tan solo unos 54
metros. Además, a pesar de haberse construido diques enrocados en ambas laderas, no son
capaces de evitar desbordes para el presente caudal de diseño. Según la figura, ambas riberas
quedarían expuestos a inundaciones ambas riberas aguas arriba del puente, llevándose la peor

79
parte los moradores que habitan en las viviendas.

Aguas abajo del Puente Socsi, por la ribera izquierda, se da inicio a bancos elevados de
material frágil (cárcavas) y susceptible frente a fuertes velocidades erosivas. Las bases de estos
bancos en la actualidad se hallan protegidos por una serie continua y consecutiva de espigones
hasta la altura del CP Caltopa. Sin embargo, en ciertos tramos se notó que los espigones están
empezando a presentar deficiencias estructurales a causa de la fuerte socavación. Por este
motivo, en los tramos en donde aún se presencia riesgos se consideró reforzar su protección
con trazos de dique a través de las zonas de empotramiento de los espigones.

Por la ribera derecha se cuenta con llanuras inundables cubiertas por cultivos y varias
viviendas. Como agravante, en estas zonas inundables se encuentra un largo tramo de la
carretera Cañete – Yauyos. Las inundaciones, según las simulaciones IBER 2D (figura 22),
suceden a pesar de que las riberas están protegidas en ciertos tramos con diques de enrocado.
Estas estructuras puedan que aun cumplan su función de proteger contra la erosión, sin
embargo, su borde libre de diseño sería insuficiente para el caudal de 1064 m3/s.

Cabe indicar que, gran parte de la carretera en peligro e inmediaciones de la Bocatoma


Nuevo Imperial pertenecen al tramo de estudio del proyecto “Valles Vulnerables”. Por este
motivo, en el tramo de “Valles Vulnerables, en este estudio se complementan tramos que faltaría
por cubrir el desborde de flujo. Dee otro lado, se emiten algunas recomendaciones de elevación
de cotas de corona para los diques existentes.

80
Figura 47. Vista de la inundación según la simulación IBER 2D aguas abajo y arriba del Puente Socsi

81
Se muestra la situación actual de los espigones y diques de enrocado, así como
montículos de material propio del cauce que actúan como diques rústicos en el tramo mostrado
en la figura 48.

Figura 48. Espigones existentes que protege el banco elevado, aguas abajo del Puente Socsi.

Figura 49. Dique enrocado junto a un dique rústico que protege la llanura inundable por donde pasa la carretera
Cañete - Yauyos, aguas abajo del Puente Socsi
Fuente: Elaboración Propia.

82
A partir de la Bocatoma hacia aguas abajo se omite ahondar en detalles; puesto que, el
río se ensancha de manera natural, comportándose como un río de llanura de material aluvial.
Dicho de otro modo, el flujo del río empieza a comportarse hidráulicamente de manera similar al
tramo Desembocadura río Cañete – Bocatoma Palo Herbay ya explicado.

5.2 Modelamiento con estructuras de defensa

El modelamiento hidráulico IBER 2D con defensas consistió en tres pasos


fundamentales: primero el trazo de ejes, segundo el diseño de la sección del dique y enrocado,
por último, la simulación numérica con el caudal de diseño (1064 m3/s) para el Tr 100 años según
los términos de referencia. La secuencia de los resultados finales se muestra a continuación:

5.2.1 Simulación IBER 2D con diques para el Tr de 100 años

Es importante recalcar que: después de una serie de modelamientos IBER 2D con el fin
de lograr resultados confiables que guarden concordancia con la información recolectada de las
imágenes satelitales, testimonios de pobladores de la zona y marcas de agua encontradas en
campo, se procedió a ensamblar los diques diseñados a la superficie topográfica. Una vez
ensamblados los diques se ejecutaron algunas simulaciones más hasta lograr que el desborde
del flujo sea controlado. Paralelamente, los datos hidráulicos se fueron actualizando en los
cuadros de cálculos hasta llegar al último modelamiento IBER 2D, cuyos resultados definitivos y
confiables se muestra en la figura 50, 51. Los detalles de profundidades y velocidades del flujo
encausado se muestran en el anexo mapas de modelamiento IBER 2D .

83
Figura 50. Vista del modelamiento hidráulico IBER 2D definitivo en presencia de diques tanto proyectados como existentes que necesitan elevación de corona.
Fuente: Elaboración Propia

84
Figura 51. Vista del modelamiento hidráulico IBER 2D definitivo en presencia de diques tanto proyectados como
existentes que necesitan elevación de corona.
Fuente: Elaboración

5.3 Trazos definitivos de eje de diques.

Tomando los criterios detallados de la metodología del modelamiento y luego de una


serie de simulaciones IBER 2D se definieron los trazos de diques márgenes derecha e izquierda
que se muestran en la figura 52 y 53 , sus progresivas fueron determinados proyectando sus
ubicaciones respecto al eje de Rio cuya progresiva comprende desde la desembocadura 0+000
km hasta Localidad Paullo 32 +000 km .Los detalles se muestran en los planos de diseño en el
anexo. Las progresivas, longitudes y respectivas márgenes de ejes de los diques se indican en
la tabla 22 y 23:

85
Tramo: Desembocadura río Cañete – Bocatoma Palo Herbay

Tabla 21. Longitud total y parcial de trazos en cada margen del tramo en estudio.

PROGRESIVA PROGRESIVA LONGITUD (m)


DIQUE (i)
INICIAL FINAL MARGEN IZQUIERDA MARGEN DERECHA
D0 0+00 0+814 814
D1 0+00 0+900 900
1+100 2+450 1350
D2
2+450 4+190 1740
1+150 2+450 1300
D3
2+450 3+195 745
D4 4+870 5+445 575
D5 5+590 6+645 1055
D6 6+080 6+694 614
D7 8+630 9+816 1186
D8 8+880 10+162 1282
D9 12+270 12+930 660
12+300 14+050 1750
D10 14+050 15+350 1300
15+350 15+605 255
D11 15+340 16+554 1214
TOTAL : 7108 9632

Adicionalmente, en el tramo de estudio se hallan dos diques existentes de 2.5 y 0.14 km


de longitud trabajando correctamente, sin embargo, se recomienda el mantenimiento y elevación
de 1.2 metros de corona para ambo

86
Figura 52. Trazos definitivos de ejes de estructuras de defensa ribereña (líneas de color magenta), valles vulnerables y diques existentes que necesitan elevación de corona.
Fuente: Elaboración Propia

87
Tramo: Bocatoma Palo Herbay – Puente Paullo

Tabla 22. Longitud total y parcial de trazos en cada margen del tramo en estudio.

LONGITUD (m)
DIQUE (i) PROGRESIVA INICIAL PROGRESIVA FINAL
MARGEN IZQUIERDA MARGEN DERECHA
16+460 18+000 1540
D13
18+000 19+840 1840
D14 16+560 17+530 970
18+800 20+700 1900
D15
20+700 21+833 1133
D16 21+180 21+899 719
D17 23+620 24+235 615
D18 24+045 24+141 96
25+580 26+600 1020
D19
26+600 27+127 527
D20 26+560 27+068 508
27+500 28+000 500
D21
28+000 28+268 268
D22 28+500 29+823 1323
D23 28+920 30+643 1723
D24 28+180 30+651 2471
30+700 30+950 250
D25
30+950 31+021 71
TOTAL: 9761 7713

Adicionalmente, en el tramo de estudio se hallan diques enrocados existentes trabajando


correctamente, sin embargo, se recomienda su elevación de 1.2 metros de corona para todos.
Debido que, las condiciones actuales de estos diques no pueden cumplir su función de evitar
desbordes para un periodo de retorno de 100 años de un caudal de 1064 m3/s. Los trazos
existentes y márgenes de estos diques también se muestran en la figura 53. A su vez, longitudes
y respectivas márgenes de los diques existente se indican en la tabla 23 .

88
Tabla 23. Longitud total y parcial de trazos en cada margen del tramo en estudio.

LONGITUD (m)
DE
MARGEN IZQUIERDA MARGEN DERECHA
DE-3 761
DE-4 274
DE-5 429
DE-6 1428
DE-7 254

Figura 53. Trazos definitivos de ejes de estructuras de defensa ribereña (líneas de color magenta), valles
vulnerables y diques existentes que necesitan elevación de corona.
Fuente: Elaboración

5.4 Diseño de la Defensa Ribereña.

La elección de diques de enrocado como estructuras de defensa ribereña es por, indicado


en el acápite de metodología del diseño hidráulico, ahora se procede a calcular las dimensiones
de los diques. Para ello se emplearon parámetros hidráulicos según las características del cauce
y áreas inundables del tramo de interés.

89
n = 0.033 manning de cauce
g= 9.81 m2/s
Q = 1064 m3/s
S = Variable
dm = 73.30 mm
Ks = 30 Manning invertido

5.4.1 Calculo Hidráulico del Tramo Dique 6 – (5+500 hasta 6+700) km

Se explicará los pasos de caculos hidráulicos del Dique 6 de la progresiva el 5+500 km


al 5+550 km

5.4.1.1 Ubicación

E tramo del dique 6 se encuentra entre el Puente Aramayo y Bocatoma la Fortaleza, en


la margen derecha del rio cañete ,teniendo una longitud de 1.1 km aproximadamente.

Fin Dique
Inicio Dique
N: 8548929
N: 8548712
E: 354437
E: 353343

Puente
Aramayo

Figura 54. Ubicación del dique 6.


Fuente: Elaboración propia

90
5.4.2 Cálculo de la velocidad media y número de froude.

5.4.2.1 Velocidad media

Se usa la fórmula de Manning-Strickler:

V =K R 2/3S1/2
m s h ………………(1)

Donde

Vm : Velocidad media (m/s)

Q : Caudal de diseño (m3/s) =1064 m3/s

R : Radio hidráulico (m) = A/P = 164.43/148.74

N : rugosidad =0.033
1
Ks : Inversa del coeficiente de rugosidad =( )= 30.30
0.033

S : Pendiente media del tramo = 0.0111

Remplazando los valores en la ecuación (1)

2
1
 164.43  3
Vm = 30.30     0.0111
2

 148.74 

Obtenemos la Vm= 3.41 m/s

5.4.2.2 Número de froude.

Se usa la fórmula para determinar el régimen:

V
FR =
g  DH
………(2)

Donde:

V : Velocidad media (m/s)=3.41 m/s

g : Aceleración de la gravedad = 9.81 m/s2.

Á𝑟𝑒𝑎 𝐻𝑖𝑑𝑟á𝑢𝑙𝑖𝑐𝑎( 𝐴ℎ) 164.43


Dh :Profundidad Hidráulica =( )= =1.12 m
𝐴𝑛𝑐ℎ𝑜 𝑑𝑒 𝑙𝑎 𝑠𝑢𝑝𝑒𝑟𝑓𝑖𝑐𝑖𝑒 𝐿𝑖𝑏𝑟𝑒 ( 𝑇) 147.32

91
Remplazando valores en la formula (2)

3.41
FR =
9.811.12

Fr=1.03 >1 por lo tanto es un flujo de régimen supercrítico.

Figura 55. Velocidad media y número de froude del dique

5.4.3 Cálculo de la Profundidad de Socavación.

Par estimar las profundidades de socavación se usó la fórmula de Lischtvan-Levedlev,


se presta para la morfología del rio cañete que son suelos no cohesivos -Granulares.

1
 5
 1+ x
 H o 3

Hs = ………1
 0.68    d 0.28 
 m 
 
α: Coeficiente a determinar

Qd
= 5 ………2
Be  Y 3
m

Ho : Tirante antes de la erosión(máxima, en el eje).


β: Coeficiente dependiente del periodo de retorno del caudal de diseño (T) =1/T
μ: Coeficiente de contracción.
Ym: Tirante medio (B>30)
Qd: Caudal de diseño.
dm: Diámetro medio de los granos de fondo
Be: Ancho efectivo o espejo de agua.
x: Coeficiente dependiente del diámetro promedio (Dm) de sólidos en el fondo del
cauce .

92
5.4.3.1 Probabilidad anual en que se presenta el gasto de diseño.

Obra n(años) R(%) Trmin(años) Tr diseño(años) Probabilidad anual(%)

Puente 40 25 140 150 0.67


Defensa Ribereña 40 25 140 150 0.67
Pontón 25 30 71 75 1.33
Alcantarilla Paso Qda. 25 30 71 75 1.33
Baden 25 30 71 75 1.33
Alcantarilla Alivio 15 35 35 35 2.86
Cuneta - Bordillo 15 40 30 30 3.33
Subdren 15 40 30 30 3.33
Recibido (Tr) 100 1.00

Figura 56. Gasto de diseño


Fuente: Elaboración Propia.

n: vida útil
R: riesgo admisible.
Tr: periodo de retorno
Usando la tabla 24 de valores del coeficiente β

Tabla 24. Valores del coeficiente β

Probabilidad anual (%) de que se


β
presente al gasto de diseño

100 0.77
50 0.82
20 0.86
10 0.9
5 0.94
2 0.97
1 1
0.3 1.03
0.2 1.05
0.1 1.07
Fuente: “Socavación en cauces naturales” de J.A Maza.

Donde el coeficiente β =1

93
5.4.3.2 Diámetro de los granos de fondo (mm).

Según análisis granulométricos, determinados por el área de geotecnia, realizaron


calicatas distribuidas a lo largo del cauce de los tramos de interés determinaron el diámetro
medio (D50) promedio de 73.3 mm, que sirve para estimar las erosiones y socavaciones.

Tenemos como dato dm =73.3 mm, utilizando la tabla 26 de valores de los suelos no
cohesivos (granulares)

1
Tabla 25. Suelos cohesivos y no cohesivos( x,y, ).
1+ 

Suelos Cohesivos Suelos no Cohesivos


Ys x 1/1+x Ys x 1/1+x Dm (mm) x 1/1+x Dm (mm) x 1/1+x
0.80 0.52 0.66 1.20 0.39 0.72 0.05 0.43 0.70 40.00 0.30 0.77
0.83 0.51 0.66 1.24 0.38 0.72 0.15 0.42 0.70 60.00 0.29 0.78
0.86 0.50 0.67 1.28 0.37 0.73 0.50 0.41 0.71 90.00 0.28 0.78
0.88 0.49 0.67 1.34 0.36 0.74 1.00 0.40 0.71 140.00 0.27 0.79
0.90 0.48 0.67 1.40 0.35 0.74 1.50 0.39 0.72 190.00 0.26 0.79
0.93 0.47 0.68 1.46 0.34 0.75 2.50 0.38 0.72 250.00 0.25 0.80
0.96 0.46 0.68 1.52 0.33 0.75 4.00 0.37 0.73 310.00 0.24 0.81
0.98 0.45 0.69 1.58 0.32 0.76 6.00 0.36 0.74 370.00 0.23 0.81
1.00 0.44 0.69 1.64 0.31 0.76 8.00 0.35 0.74 450.00 0.22 0.82
1.04 0.43 0.70 1.71 0.30 0.77 10.00 0.34 0.75 570.00 0.21 0.83
1.08 0.42 0.71 1.80 0.29 0.78 15.00 0.33 0.75 750.00 0.20 0.83
1.12 0.41 0.71 1.82 0.28 0.78 20.00 0.32 0.76 1000.00 0.19 0.84
1.16 0.40 0.71 2.00 0.27 0.79 25.00 0.31 0.76
Fuente: “Socavación en cauces naturales” de J.A. Maza(1968).

Observamos que nuestro dm =73.3 mm , se encuentra en el intervalo 60 <dm<90,


procederemos realizar una interpolación para obtener nuestro coeficiente x.

60 ---------> 0.29
73.3 --------> x
90 --------->0.28

 73.3 − 60 
x=   ( 0.28 − 0.29 ) + 0.29
 90 − 60 

x= 0.286---------> tomando el valor x=0.29

1/(1+x)=0.778

94
Teniendo los valores Be =147.32 m, Qd=1064 m3/s, Ym=1.64 m

Calcularemos el Coeficiente α:

Remplazando los valores en la ecuación (2):

1064
= 5
147.32 1.64 3

α= 3.16 m/s

obtenido los valores x=0.29, dm=73.3 mm α= 3.16 m/s,Ho=2.23 m, β =1 procederemos


remplazar en la ecuación (1).
1
 5
1+ 0.29
 3.16  2.23 3 
Hs =
 0.68  1 0.07 0.28 
 
 

Hs=3.68 m

5.4.3.3 Profundidad efectiva de socavación general del cauce principal.

Ds= Hs-Ho……3

Donde :

Ds: Profundidad socavado.

Ho: Tirante del agua para el cauce no socavado.

Hs: Tirante de agua par cauce socavado.

Figura 57. Esquema de la socavación del cauce de rio.


Fuente: Elaboración Propia.

95
Teniendo los valores Ho=2.23 m, HS=3.68 m, reemplazando en la ecuación (3)

Ds=3.68 m -2.23 m

Ds=1.45 m

Nota: En este caso por tratarse de un rio trenzado de espejo de agua muy amplio (B>30m)
se tomó el calado estimado a partir del método Manning-Strickler. En caso de quebradas
se tomaría el tirante correspondiente al Ah/Be.

Figura 58. Calculo de la profundidad de socavacion dique 6

5.4.4 Cálculo de cota y Profundidad de Uña(P Uña).

Cota de socavación = Cota eje de rio - P uña…1

Cota name=57.98 msnm

Cota eje rio= 56.15 msnm

Ds=Puña=1.45 m

Remplazando los valores en la ecuación (1)

Cota de socavación=56.15 m - 1.45 m

Cota de socavación = 54.70 m

96
Figura 59. Cotas y profundidad de la uña dique 6
Fuente: Elaboración Propia.

5.4.5 Cálculo de la altura del enrocado

5.4.5.1 Borde libre

Para el borde libre se consideró las sugerencias de investigadores peruanos Campaña


(2004) . “El borde libre permite controlar la variación instantánea del caudal por disminución de
la velocidad y elevación del tirante”(Teran,1998,p.36).

Vm2
BL =   ……1
2 g

H: Altura del enrocado (m)


BL: Borde libre (m)
Vm: velocidad media del flujo (m/s)
g: Aceleración de la gravedad (m/s2)
: Coeficiente en función de la máxima descarga.

El coeficiente (ø) en función del caudal se usará la tabla Rubén Terán(1998) nuestro
caudal de diseño 1064 m3/s.

97
Tabla 26. Valores del coeficiente ø en función del caudal

Caudal Máximo
Coeficiente (ø )
(m3/s)

3000-4000 2.00
2000-3000 1.70
1000-2000 1.40
500-1000 1.20
100-500 1.10

Fuente: “Diseño de defensas Ribereñas-Rubén Terán”.

Tabla 27. Recomendación de borde libre

Borde Libre
Descarga de diseño (m3/s)
(m)
Menos de 200 0.6
200-500 0.8
500-2000 1.08

Fuente: Campaña(2004).

 = 1.40

Valores obtenidos remplazamos en la ecuación (1)


3.412
BL = 1.40 
2  9.81

Obteniendo el BL calculado=0.83 m se redondeó a 1 m

Luego se procedió calcular borde libre máximo

BL máximo =1 +0.20, se consideró 0.20 m por la sedimentación.

BL máximo =1.20 para todo el tramo del interés de estudio.

98
5.4.5.2 Altura del enrocado (H)

Es la diferencia entre la cota del nivel de flujo y socavación, se suma el borde libre. El
valor de la cota del flujo y la velocidad se obtienen del simulamiento IBER 2D.

H = Cota de flujo - Cota de socavación + BL

H= 57.98 m - 54.70m +1.20 m

Finalmente, H calculado = 4.48 m

En el tramo del Dique 6,se consideró la máxima altura promedio H elegida =5.6 m

Figura 60. Altura del enrocado del dique 6


Fuente: Elaboración Propia.

5.4.6 Cálculo del diámetro medio del enrocado

Hay varias fórmulas para estimar el diámetro medio de la roca, de los cuales se ha
considerado tomar en cuenta los siguientes:

5.4.6.1 Maynord.

Establece que la roca para su espesor depende de los parámetros (talud y la forma de
planta ) del rio cuya formula es :

99
D50 = t  C1  F 3 …...1
C2  vm
F= …...2
g  ym

D50 : Diámetro medio de las rocas (m).


F : Número de Froude
Ym : Profundidad media del flujo =1.64 m m/s
Vm : Velocidad media del flujo. =3.41 m/s
Los valores recomendados de C1=0.32 y C2 =1.25 se muestran en la figura:

0.28 Fondo Plano


C1 0.28 Talud 1V:3H
0.32 Talud 1V:2H
1.5 Tramos en curva
C2 1.25 Tramos rectos
2 Extremos de espigones

Remplazando los valores en la ecuación (2):

1.25  3.41
F=
9.811.64

Número de Froude= 1.06

Para obtener el D50 remplazamos los valores en la ecuación (1):

D50 = 1.64  0.32 1.063

D50 = 0.63 m

5.4.6.2 Método de Isbash

Este método se basa en la siguiente expresión:

V2
D50 = 0.59  …...1
W g

kr − ka
w= ……2
Ka

D50 : Diámetro medio de la roca

100
V : Velocidad media del flujo=3.41 m/s
kr : peso específico de la roca =2.65
ka : peso específico del agua=1.0
g : Aceleración de la gravedad=9.81 m/s2

Los valores remplazando en la ecuación (2):


2.65 − 1
w=
1
W=1.65

Remplazando los valores en la ecuación (1):


3.412
D50 = 0.59 
1.65  9.81

D50 =0.42 m

5.4.6.3 Método del U.S. Department of Transportation (U.S.D.T)

D50 = C  d 50t
……… 1

 0.00594  Vm3 
D150 = c    ………2
 ym  k1
0.5 1.5

D50 : Diámetro medio del enrocado de protección


K1 : Coeficiente de los taludes y ángulos de reposo del enrocado de protección
Ym : Tirante medio del flujo al pie del enrocado
V : Velocidad media del flujo al pie del enrocado .
C : Coeficiente de corrección.

0.5
  sin 2   
K1 = 1 −  2   ……3
  sin   

2.12
Csg = ……4
( ys − 1)
1.5

1.5

=(
FS ……5
C sf )
1.2

101
C = Csg  Csf …...6

Dónde:

ɸ: Ángulo de reposo del material del enrocado =42.0 °


Θ: Ángulo de inclinación del talud protegido con enrocado=26.6°
Csg : Coeficiente de corrección por el peso específico
Csf : Coeficiente de corrección por factor de seguridad
ɣ𝑠 : Gravedad específica del material del enrocado =2.64 ton/m3
FS : Factor de seguridad, según la siguiente tabla 30. Fs=1.5.

Tabla 28. Consideraciones para el factor de seguridad.

Condición Rango del Fs

Flujo uniforme, tramos rectos 0 medianamente curvos (radio


de la curvalancho de! cauce > 30). Mínima influencia de 1.0-1.2
impacto de sedimentos y material flotante.

Flujo gradualmente variado, curvatura moderada (10 < radio


) de la curva/ancho del cauce < 30). Moderada de impacto de 1.3-1.6
sedimentos y material flotante.

Flujo rápidamente variado, curvas cefradas (radio de la


curvalancho del cauce < 10), flujos de alta turbulencia, flujo
1.6-2.0
de turbulencia mixta en estribos de puentes. Efecto
significativo de impacto de sedimentos y material flotante.

Fuente: Manual de Hidrología, Hidráulica y Drenaje (MTC, sf)

Obteniendo los valores K1 =0.74, Csg=1.101, Csf=1.40 y C=1.41, Ys=2.64Ton/m3


reemplazando en la ecuación (2):

 0.00594  3.413 
D150 = 1.41  1.5 
 1.64  0.74 
0.5

D´50 =0.29 m

Teniendo el valor D´50 =0.29 m ,C =1.41 remplazamos en la ecuación (1):

D50 = 1.41 0.29

102
Finalmente tenemos el D50=0.41 m

Figura 61. Diametro medio del enrocado del dique 6.


Fuente: Elaboración Propia.

5.4.7 Diseño definitivo del Dique 6 de enrocado

5.4.7.1 Espesor mínimo del enrocado (T)

Las consideraciones para tener en cuenta son:

• No debe ser menor que el tamaño máximo de la roca a usarse en el enrocado.

• No debe ser menor de 1.5 x D50.

• Por razones constructivas no deberá ser menor de 30 cm.

T = 2  D50

T = 2  0.6

Finalmente, el valor T =1.3 m

103
5.4.7.2 Ancho de la base (A)

El ancho base de la uña está en función de la profundidad de socavación (P uña), por lo


que, su fórmula es:

A = 1.5 P uña……1

A : Ancho del enrocado en la base (m) y

P uña : Profundidad de socavación o de la uña del enrocado (m).=2.1 m

Los valores remplazamos en la ecuación (1)

A = 1.5  2.1

Obteniendo que el ancho de la base del dique 6--------> A=3.10 m

Figura 62. Dimensiones del dique de enrocado del dique 6


Fuente: Elaboración Propia.

La pendiente fue variable porque cada tramo de dique, obteniendo tirantes, áreas
hidráulicas, demás parámetros, así como los resultados y cálculos de los 25 tramos se muestran
detalladamente en el anexo de cuadros de cálculos .

104
Figura 63. Representación de las dimensiones del enrocado

Figura 64. Seccion tipica del enrocado


Fuente: Elaboración Propia.

Figura 65. Leyenda del diseño


Fuente: Elaboración Propia.

105
6. CONCLUSIONES

• Se identificó el riesgo de áreas expuestas a inundaciones, mediante


modelamientos IBER, para caudales 625m3/s,1064m3/s y 1064 m3/s
correspondientes a 10,100 y 500 para un periodo de retorno respectivamente .

• Se identificó el riesgo de las riberas expuestas a erosión, basándose en la


migración lateral, observadas en el registro histórico de imágenes satelitales
disponibles desde el año 2002.

• Se definieron trazos de eje, para estructuras de defensa ribereña, cuya longitud


total de la margen izquierda es de 14.513 km y de la derecha de 17.406 km.

• Se diseñaron diques de enrocado como medidas estructurales de defensa para


el tratamiento y protección de sectores vulnerables del tramo en estudio (ver
cuadros de cálculos y planos en anexos).

• Se decidió mantener los diques existentes ubicados aguas abajo de la Bocatoma


Palo Herbay y La Pinta, sin antes recomendar su mantenimiento y elevación de
cresta de 1.2 m para evitar desbordes en ambos.

• Se necesita un volumen de material de 21840 m3, para la elevación de la corona


del dique existente de 2.5 km y un volumen de 3605 m3 para un corto tramo de
0.41 km de longitud del otro dique.

106
7. RECOMENDACIONES

• Se deben ejecutar las defensas ribereñas en los tramos vulnerables para proteger
de inundaciones a las poblaciones afectadas.

• Para impedir la degradación de las márgenes es también necesario colocar una


capa de geotextil a manera de filtro para impedir la migración de los finos.

• Sensibilizar a la población afectada por inundaciones, para que eviten la


utilización de los terrenos contiguos a la ribera del rio .

• Se debe retirar del cauce un volumen total 26820 m3 de material perteneciente a


los diques colapsados que obstruyen el libre paso del flujo en tramos estrechos.

107
8. REFERENCIAS BIBLIOGRAFICAS

Agencia de Cooperación Internacional del Japón (2013).Estudio de Factibilidad Instalación -


implementación de medidas de prevención para el control de desbordes e inundaciones
del Rio Cañete- Departamento de Lima .

Campaña, R.(2004). Introducción a los Proyectos de Defensa Ribereñas. Perú: Ministerio de


Agricultura Instituto Nacional de Recursos Naturales.

Cconislla Quispe, J. A. (2017). Protección contra inundaciones en el sector del centro poblado
de Gorgor, distrito de Gorgor, Provincia de Cajatambo, Lima. Obtenido de
http://repositorio.lamolina.edu.pe/handle/20.500.12996/3376

Consorcio SALZGITTER-SISA (2019). Instalación - implementación de medidas de prevención


de desbordes e Inundaciones del río Cañete - Provincia Cañete- Departamento de Lima
– Perú.

Gobierno Regional de Lima - Dirección regional de Agricultura (2019). Expediente definitivo:


Creación e implementación de medidas de prevención para el control de desbordes e
inundaciones del rio Cañete: tramo puente Socsi – puente san Jerónimo, en los distritos
de Lunahuaná, Pacarán y Zúñiga de la provincia de Cañete - departamento de lima.

IBER(2010). “Modelización bidimensional del flujo en lámina libre en aguas poco profundas”.

Martin, J. (2001). Ingeniería de Ríos. España: Alfa Omega Grupo Editor

Maza J.(1968). Socavación en cauces naturales. México: secretaria de obras públicas.

Rocha, F. (1998) .Introducción a la hidráulica fluvial. Perú: Facultad de Ingeniería Civil


Universidad Nacional de Ingeniería.

Rocha, F. (2007) .Hidráulica de Tuberías y Canales. Perú: Facultad de Ingeniería Civil


Universidad Nacional de Ingeniería.

Terán, R. (1998). Diseño y construcción de defensas ribereñas. Perú: Escuela Superior de


Administración de Aguas CHARLES SUTTON.

Tito Quispe, Y. A. (2017). Modelamiento hidráulico del río Cañete sector puente Socsi - altura
puente colgante (9 km), con fines de diseño de defensas ribereñas. Obtenido de
http://repositorio.lamolina.edu.pe/handle/20.500.12996/3053

108
9. ANEXOS

109
MAPAS DE INUNDACIONES EN IBER 2D

10 ,100 y 500 años

110
Tramo desembocadura-Bocatoma Palo Herbay Mapa de inundación sin estructuras 10 años

111
Tramo desembocadura-Bocatoma Palo Herbay Mapa de inundación sin estructuras 100 años

112
Tramo desembocadura-Bocatoma Palo Herbay Mapa de inundación sin estructuras 500 años

113
Tramo desembocadura-Bocatoma Palo Herbay Mapa de inundación con estructuras 100 años

114
Tramo bocatoma Palo Herbay- Localidad Paullo Mapa de inundación sin estructuras 10 años

115
Tramo bocatoma Palo Herbay- Localidad Paullo Mapa de inundación sin estructuras 100 años

116
Tramo bocatoma Palo Herbay- Localidad Paullo Mapa de inundación sin estructuras 500 años

117
Tramo bocatoma Palo Herbay- Localidad Paullo Mapa de inundación con estructuras 100 años

118
CUADROS DE CALCULOS

Diseño de diques enrocados, Tr 100 años

119
120
121
122
123
124
125
126
127
128
129
130
PLANOS DE DISEÑO

Tr 100 años

131
132
133
134
135
136
137
138
139
140
141

También podría gustarte