HIDROLOGIA

DE LAS QUEBRADAS SAN JERONIMO, CONDORAY Y JACAYITA

1. ANTECEDENTES

Los Flujos de Escombros denominados comúnmente “huaycos”, son eventos naturales que
suelen presentarse con frecuencia en épocas de avenidas y/o fenómenos de EL NIÑO, que
ocurren en nuestro medio. Estos eventos geodinámicos son cíclicos por lo que cada cierto
tiempo presentan un peligro inminente para las poblaciones e infraestructuras involucradas
próximas al evento.

La cuenca Cañete no es ajena a este fenómeno natural por lo que cada cierto tiempo sus
quebradas son activadas por este flujo torrentoso que generalmente arrasa literalmente
todo a su paso. Es el caso de las quebradas San Jerónimo, Condoray y Jacayita, que
sufrieron su evento más extremo el 2017, que incluso llegaron a bloquear vías de acceso.

El proyecto: “Formulación del Plan Integral para el Control de Inundaciones y Movimientos

de Masa de la Cuenca del Río Cañete”, según los TdR, contempla el desarrollo del
“Componente B”, considera la formulación de estudios a nivel de perfil de proyectos de gran
impacto, que representen la disminución del riesgo en 60% en el ámbito de intervención;
por lo que el área de Formulación ha considerado incluir a las quebradas san Jerónimo,
Condoray y Jacayita como proyectos a ser mitigado .

En este sentido, existe la propuesta de realizar estructuras de encauzamientos en estas

quebradas, para atenuar los huaycos (flujo de escombros). Por lo tanto, es imprescindible
generar información hidrológica en las secciones de interés para poder evaluar la
concepción y magnitud de las obras a proyectarse.

Como se sabe el flujo que generan estos eventos, es un flujo del tipo hiperconcentrado; que
contiene una concentración considerable de sólidos y agua. Es por este motivo que se ha
generado hidrogramas de descargas máximas de agua para diferentes periodos de retorno
(50, 100 y 500) para las quebradas. Posteriormente también se realiza el cálculo de
caudales de flujo de escombros (caudales totales de sólido y liquido) que sirve para el
diseño de estructuras de mitigación de huaycos.
2. DESCRIPCIÓN DE LA QUEBRADAS

2.1 Ubicación

Políticamente, las subcuencas que forman las quebradas San Jerónimo, Condoray y
Jacayita forma parte del distrito de Lunahuana de la provincia de Cañete, perteneciente al
departamento de Lima.

Las quebradas se encuentran ubicada sobre la margen izquierda del río cañete, estas tres
quebradas están ubicadas una de tras de otra, pero tienen como cuencas vecinas a ddd
por el norte, wwww por el sur, qqq por este y www por oeste. En el siguiente cuadro se
muestra su ubicación en coordenadas UTM:

Cuadro HI - N° 1: Coordenadas UTM de la quebradas de interés.

Quebrada Este Norte

San Jerónimo 374484 - 389765 8562466 - 8565868
Condoray
Jacayita

En el Plano 018-2019-HID - 01 del Anexo VII, se presenta la cuenca de las quebradas de

interés.

2.2 Geomorfología de la Cuenca

Las quebradas presentan una cuenca pequeña, teniendo un área que va de 8.9 a 32.6 km2;
el mayor afluente de estas cuencas son las quebradas del mismo nombre, aunque la cuenca
Condoray tiene como afluente principal a la quebrada Escalón. En el Cuadro HI - N° 2 se
resume los parámetros fisiográficos de las cuencas.

Cuadro HI - N° 2: Parámetros Geomorfológicos de las quebradas de interés.

QUEBRADA San Jerónimo Condoray Jacayita

Mayor Aflu San Jeronimo Escalon

Área (km2) 32.683 32.4 8.922
Perímetro (km) 26.888 28.87 13.313
 QUEBRADA San Jerónimo Condoray Jacayita

Longitud de Cauce (m) 12462.2 12308 3624.63

Altura Media ( msnm) 1299 1751 1389

Coeficiente de
1.3 1.3 1.2
compacidad

Pendiente 58% 81% 109%

Tiempo de Concentración
0,93 0.84 0,25
(min)
Lado Mayor (km) 10.2 11.4 4.7
Lado Menor (km) 3.1 2.9 1.8

Altura Max (msnm) 2225 2920 2562

Altura Min (msnm) 368 496 542

3. MÁXIMAS AVENIDAS

La estimación de las descargas máximas se ha realizado a partir de la información

disponible, utilizando registros de precipitación máxima de 24 horas, mediante un análisis
de frecuencias y/o un modelo precipitación – escorrentía.

Para el Análisis de Frecuencias de Distribución Extrema se utilizó el software

HYDROGNOMOM, desarrollando las siguientes distribuciones: Distribución Normal,
Distribución Log Normal, Distribución Pearson III, Distribución Log Pearson III, Distribución
Gumbel, Distribución GEV- L Moments.

Para el desarrollo del Modelo de Precipitación – Escorrentía, se ha utilizado el Software

HEC-HMS. Este modelo se calibra mediante la utilización de los registros históricos de
caudales máximos diarios anuales de la estación Socsi y de una manera referencial la
escasa información disponible de otras estaciones de control que forman parte de la red
hidrométrica .

La generación de los caudales de avenida en las secciones de cuencas de interés, se han

obtenido para diferentes periodos de retorno: 2, 5, 10, 25, 50, 100, 200 y 500 años.

3.1 Estaciones Pluviométricas

Las estaciones pluviométricas utilizadas para el análisis de Máximas Avenidas en las

cuencas San Jerónimo, Condoray y Jacayita, son las estaciones que pertenecen a la
 cuenca Cañete. La precipitación de la cuenca es registrada en la actualidad por siete
(07) estaciones pluviométricas: Yauricocha, Tanta, Carania, Vilca, Yauyos, Huangascar,
Pacaran y Cañete. Asimismo, se cuenta con registros de 4 estaciones ya desactivadas:
Colonia, Huantan, Catahuasi, Siria y Sunca. Para el presente estudio, se ha considerado
los registros de las siete primeras estaciones, más los registros de las estaciones de
Colonia y Huantan por contar con data pluviométrica de más de 20 años.

El Cuadro HI – N° 3 muestra la relación de las estaciones consideradas en el presente

estudio con sus respectivas características de código, nombre y ubicación.

Cuadro HI - N° 3: Estaciones Pluviométricas de la cuenca Río Cañete

N° Estaciones Tipo Procedencia Operación Cuenca Dpto Prov. Longitud Latitud Altitud Periodo
1 Tanta Plu Senamhi Senamhi Cañete Lima Cañete 76°01' 12°07' 4323 1964-2018
2 Carania Plu Senamhi Senamhi Cañete Lima Cañete 75°52' 12°20' 3820 1964-2018
3 Vilca Plu Senamhi Senamhi Cañete Lima Cañete 75°49' 12°06' 3832 1964-2018
4 Colonia Plu Senamhi Senamhi Cañete Lima Cañete 75°53' 12°38' 3308 1964-1988
5 Huantan Plu Senamhi Senamhi Cañete Lima Cañete 75°49' 12°27' 3219 1964-1990
6 Yauyos Plu Senamhi Senamhi Cañete Lima Cañete 75°54' 12°29' 2285 1964-2018
7 Huangascar Plu Senamhi Senamhi Cañete Lima Cañete 75°50' 12°53' 2500 1964-2018
8 Pacaran Plu Senamhi Senamhi Cañete Lima Cañete 76°03' 12°52' 684 1964-2018
9 Yauricocha Plu Senamhi Senamhi Cañete Lima Cañete 75°42 12°18' 4560 1964-2018

3.2 Precipitaciones Máximas en 24 horas

Tal como se indicó en el Cuadro HI - N°3, se reportan 09 estaciones meteorológicas que

cuentan con la información de Precipitación máxima en 24 horas y presentan periodos de
registro variables. En el Anexo I se presentan los registros históricos de Precipitaciones
máximas en 24 horas adquiridas del SENAMHI; la mayoría posee desde el año 1964, con
excepción de la estación de Cañete que posee data desde el año 1926, existiendo en varias
estaciones algunos registros discontinuos, con datos faltantes.

3.3 Análisis de la Precipitación Máxima en 24 Horas

Para el análisis de la precipitación máxima en 24 horas, previamente se hizo el análisis de

periodos con datos inconsistentes, tal como se realizó en la precipitación total mensual,
encontrando los mismos periodos inconsistentes que no fueron tomando en cuenta para el
análisis.
 Para ello, en primer lugar, el análisis se realizó visualmente, donde se detectaron datos
atípicos, que son definidos como registros que se alejan significativamente de la tendencia
observada en la muestra (suponiendo que se encuentra normalmente distribuidos). La
utilización o eliminación de estos datos puede afectar significativamente la magnitud de los
parámetros estadísticos, especialmente en muestras pequeñas. La gran mayoría de estos
datos coinciden con los datos que fueron observados en la consistencia de la información
pluviométrica histórica, lo que podría confirmar que estos datos si son diferentes al resto.

Contando con los registros oficiales obtenidos del SENAMHI, se pudo apreciar que tenemos
como mínimo 30 años de periodo de registro en la estación Yauricocha.

Así mismo, se procedió a corregir la data por medio de un factor de corrección Bell (1969),
con el fin de calcular la precipitación máxima instantánea, así como se muestra en la
siguiente expresión:

Pmáx. Ins. = Pmáx. diaria x FC

Se utiliza un factor de corrección (FC) de 1.13 cuando los registros de precipitaciones

máximas diarias se realicen una vez al día, pero si se registran dos lecturas diarias, como
es el caso de las estaciones en la cuenca, el factor que se usaría es de 1.04. A continuación,
se muestra en el Cuadro HI - N°4, las Precipitaciones máximas instantáneas calculadas
para las estaciones utilizadas.

Cuadro HI - N° 4: Precipitaciones Máximas Instantáneas para la Cuenca del Río

Cañete (mm).

SAN PEDRO
YAUYOS PACARAN TANTA CARANIA HUANGASCAR VILCA YAURICOCHA
DE PILAS
AÑO
2285 4323 3820
680 msnm 2707 msnm 2500 msnm 3832 msnm 4560 msnm
msnm msnm msnm
1964 20.3 - - 26.4 29.5 - -
-
1965 32.7 2.2 - 35.9 46.1 15.6 - -
1966 24.2 2.6 - 27.7 26.0 26.1 - -
1967 24.5 9.2 - 29.1 19.3 36.7 - -
1968 - - - 24.6 - 13.4 - -
1969 18.1 - - 34.3 30.5 22.2 - -
1970 27.9 21.1 - 39.4 17.3 29.1 - -
1971 34.3 6.6 - 25.5 18.7 20.4 - -
1972 - 5.0 - 27.1 20.9 73.3 - -
1973 29.3 6.2 - 18.9 23.5 28.3 - -
1974 22.4 2.5 - 20.1 17.6 13.2 - -
1975 19.8 3.4 - 15.7 16.6 36.0 - -
 SAN PEDRO
YAUYOS PACARAN TANTA CARANIA HUANGASCAR VILCA YAURICOCHA
DE PILAS
AÑO
2285 4323 3820
680 msnm 2707 msnm 2500 msnm 3832 msnm 4560 msnm
msnm msnm msnm
1976 20.8 0.4 - 18.2 20.1 - - -
1977 15.4 0.2 - 17.1 18.1 30.6 - -
1978 20.9 0.2 - 17.0 16.7 51.8 - -
1979 17.6 - - 12.2 15.7 18.8 - -
1980 16.1 - - 15.0 17.8 8.8 - -
1981 23.7 - - 13.6 18.2 21.8 - -
1982 - - 13.8 16.2 17.9 - -
1983 - - 17.3 10.1 - -
1984 10.4 - - 11.8 14.8 15.5 - -
1985 - - 12.9 13.4 14.4 - -
1986 3.6 18.7 20.8 19.8
16.6 - -
1987 5.0 22.8 17.5 21.7 13.6 - 39.1
1988 3.4 18.4 14.4 34.4 21.2 - 30.0
1989 6.2 16.1 14.5 25.4 20.8 - 27.1
1990 1.2 16.4 29.1 20.8
12.9 - 32.0
1991 1.6 12.0 12.9 19.8
19.2 - 25.0
1992 6.6 1.2 16.6 15.7 5.2
10.9 - 22.4
1993 18.0 3.1 45.8 43.3 16.6 20.8 - 42.1
1994 32.8 9.4 19.8 27.5 14.7 25.0 - 22.7
1995 12.7 6.4 24.3 28.1 14.0 31.2 - 21.0
1996 25.6 13.3 33.0 16.7 23.9 18.2 17.3
1997
19.6 3.7 15.8 28.5 15.2 26.3 35.4 29.3
1998
15.3 5.7 31.3 43.5 14.7 35.2 21.1 28.7
1999 20.7 11.6 25.7 25.5 16.2 25.3 18.9 25.4
2000 13.4 4.0 11.6 30.1 28.1 31.8 20.9 60.9
2001 13.8 5.8 15.9 23.6 15.5 13.3 16.8 21.4
2002 12.1 6.1 24.4 29.3 18.4 25.8 21.6 26.8
2003
15.0 4.6 13.6 29.1 19.7 15.6 0.0 62.8
2004
14.8 - 9.0 34.2 22.3 18.4 32.4 43.0
2005
14.1 - 10.1 22.9 21.3 13.5 17.7 31.6
2006 21.4 3.6 33.4 30.7 31.3 25.4 26.5 27.2
2007 20.6 2.4 16.3 34.9 24.3 15.2 31.2 30.2
2008 20.7 2.7 28.8 0.0 22.8 25.0 23.6 16.0
2009
15.7 6.2 27.1 34.3 21.3 15.4 41.6 28.0
2010
18.7 3.6 43.9 33.8 21.1 23.5 24.3 28.7
2011
24.1 4.1 10.4 24.2 25.8 15.8 27.9 14.8
2012 34.3 14.0 42.1 29.8 23.4 32.9 27.6 36.7
2013 22.6 2.1 24.4 40.1 18.8 21.8 22.4 16.7
2014 29.1 5.9 23.2 26.3 21.1 19.2 24.8 14.7
2015 13.7 10.3 24.9 34.5 19.4 16.8 30.6 37.5
2016
18.4 1.7 15.0 32.2 17.5 17.3 33.3 33.8
2017
26.0 12.6 26.4 40.1 40.0 64.6 32.0 33.8
 SAN PEDRO
YAUYOS PACARAN TANTA CARANIA HUANGASCAR VILCA YAURICOCHA
DE PILAS
AÑO
2285 4323 3820
680 msnm 2707 msnm 2500 msnm 3832 msnm 4560 msnm
msnm msnm msnm
2018 13.5 2.7 13.7 18.0 24.8 45.0 29.7 24.8

Cada serie de precipitaciones máximas de las nueve (09) estaciones pluviométricas

seleccionadas ha sido evaluada con distintas funciones de distribución, por lo que a
continuación se describe las funciones de distribución más comunes que se aplica a los
estudios hidrológicos de eventos extremos.

3.4 Análisis Estadístico y Prueba de Bondad de Ajuste.

El análisis de frecuencia consiste en aplicar las Funciones de Distribuciones de Frecuencia

a la serie de caudales máximas determinadas.

Se han seleccionado las siguientes distribuciones, por ser las más usadas para determinar
los caudales máximos.

 Distribución Normal.

 Distribución Log Normal.

 Distribución Pearson III.

 Distribución Log Pearson III.

 Distribución Gumbel.

 Distribución GEV- L Moments.

Distribución Normal: La distribución normal es una distribución simétrica en forma de

campana, también conocida como Campana de Gauss. Esta distribución se define:

1 1 𝑥−𝜇 2
− ( )
𝑓(𝑥) = 𝑒 2 𝜎
√2𝜋𝜎

Dónde:
 f(x): Densidad de probabilidad.
 x: Variable

Distribución Log Normal: Esta distribución es una extensión de la distribución

Normal, en la cual los valores logarítmicos de una secuencia son considerados como
normalmente distribuidos.
 1 1 𝑦−𝜇𝑦 2
− ( )
𝑓(𝑥) = 𝑒 2𝜎 𝑦
𝑥𝜎𝑦 √2𝜋

Dónde:
 y = log(x)
 f(x): Densidad de probabilidad.
 x: Variable
 µy: Media de los logaritmos naturales de x.
 σy: Desviación estándar de los logaritmos naturales de x.

Distribución Pearson III: Esta distribución posee una gran flexibilidad y diversidad de
forma, dependiendo de los valores de sus parámetros, asimilando su utilización para
precipitaciones o caudales máximos anuales. La función de densidad de probabilidad
Pearson III se define como:

1 𝑥 − 𝛿 𝛽−1 −𝑥−𝛿
𝑓(𝑥) = { } 𝑒 𝛼
∝ Γ(𝛽) 𝛼

Dónde:
 α, β, δ son parámetros de la función
 f(x): Densidad de probabilidad.
 x: Variable

Distribución Log Pearson III: Esta distribución ha sido una de las más utilizadas en
hidrología. Dado que la mayoría de las variables hidrológicas son sesgadas.

(ln 𝑥−𝑥0 )
−
(ln 𝑥 − 𝑥0 )𝛾−1 𝑒 𝛽
𝑓(𝑥) =
𝑥𝛽 𝛾 𝛤(𝛾)

Válido para:
 x0 ≤ x < ∞ ; -∞ < x0 < ∞ ; 0 < β < ∞ ; 0 < γ < ∞
Dónde:
 x0: Parámetro de posición.
 γ: Parámetro de forma.
 β: Parámetro de escala.
 Distribución Gumbel: Esta distribución se deriva de otras distribuciones de tipo
exponencial que convergen hacia una función exponencial a medida que la variable
aleatoria (x) crece. La función de densidad de probabilidad es:

−𝛼(𝑥−𝛽) ]
𝑓(𝑥) =∝ 𝑒 [−∝(𝑥−𝛽)−𝑒

Donde a y b son los parámetros de la función. Los parámetros α y β, se estiman para

muestras muy grandes:

Distribución GEV- L Moments:

Esta distribución fue introducida por Hosking (1990) y constituye un sistema alternativo
para describir las formas de las funciones de distribución de probabilidades y se derivan
de los momentos ponderados probabilísticamente (MPP) de Greenwood et al. (1979).
Los MPP de la variable aleatoria X con una función de distribución F(X) se definen por:

Se eligió la distribución de Pearson III, por ser la que muestra mejores ajustes de acuerdo
a las pruebas de Kolmogorov – Smirnov y a las tendencias de la información de la zona de
estudio.

3.5 Pruebas de Bondad de Ajuste

Las pruebas de bondad de ajuste, consiste en comprobar gráfica y estadísticamente, si la

frecuencia de la serie analizada se ajusta a una determinada función de probabilidad teórica
seleccionada a priori en base a los parámetros estimados.

Las pruebas de ajuste estadísticos, han sido ejecutadas por el método Kolomogorov -
Smirnov a través del Software Hydrognomon. En el Anexo II, se presenta las pruebas
estadísticas de Kolomogorov – Smirnov para las estaciones y la gráfica de ajuste de los
valores calculados. En el Cuadro HI - N° 5, se presenta valores para cada estación
pluviométricas y para cada periodo de retorno.
 Cuadro HI - N° 5: Precipitación Máxima en 24 Horas estimada para diferentes
periodos de retorno – Cuenca del Río Cañete.

PERIODO DE RETORNO (años)

ESTACIÓN
2 5 10 25 50 100 200 500 1000
YAUYOS 18.6 24.1 27.4 31.1 33.7 36.2 37.8 40.7 42.8
PACARAN 3.8 7.4 10.1 13.7 16.3 18.9 22.4 26 28.7
SAN PEDRO DE PILAS 20.2 29.3 35.1 42.1 47.1 51.9 56.6 62.6 67
TANTA 26.5 34.7 39.5 44.8 48.4 51.8 53.9 57.8 60.6
CARANIA 19.7 25.9 30.4 36.1 40.4 44.6 47.9 53.3 57.3
HUANGASCAR 20.4 32.1 41 52.8 61.8 70.7 78.2 89.8 98.6
VILCA 26.3 32.2 35.5 39.4 42 44.4 46.7 49.6 51.7
YAURICOCHA 27.8 38.5 45.7 54.7 61.2 67.6 72.5 80.6 86.6
CAÑETE 1.7 4 6.2 9.4 12.0 14.7 17.5 20.9 23.7

Con la información mostrada en el cuadro anterior y con la ayuda del software de

información geográfica Qgis, se ha generado la distribución espacial de la precipitación para
cada periodo de retorno mediante la construcción de Isoyetas.

A partir de los mapas de Isoyetas para cada periodo de retorno, se ha estimado las
precipitaciones máximas para cada una de las subcuencas de interés. Los valores
calculados se muestran en el Cuadro HI - N°6.

Cuadro HI - N° 6: Precipitaciones Máximas Diarias para la Subcuenca y

Periodos de Retorno (mm)
PERIODOS DE RETORNO
N° Nombre
2 Años 5 Años 10 Años 25 Años 50 Años 100 Años 200 Años 500 Años 1000 Años

33 San Jerónimo 15.4 24.2 30.8 39.4 45.9 52.3 58.2 66.5 72.8
46a Condoray 3.1 6.7 9.6 13.7 16.7 19.7 22.8 27.6 30.8

Jacayita

3.6 Curvas de Intensidad - Duración - Frecuencia (IDF)

Es una curva que involucra una relación entre la intensidad de la lluvia, la duración y las
frecuencias o periodos de retorno apropiados. La duración se ubica en el eje horizontal y la
intensidad en el eje vertical, mostrando una seria de curvas, para cada uno de los periodos
de retorno de diseño.
A continuación, se muestran las Curvas IDF para las subcuencas de interés.
 Gráfico HI - N° 1: Curvas IDF de la subcuenca San Jerónimo.

35.00
1000 años 500 años 200 años
100 años 50 años
30.00
Intensidad de Precipitación (mm/hr)

0.17
25.00 I= 23.03 T
.0.527
Donde: t
20.00 T= años
t= minutos
15.00

10.00

5.00

0.00
0 100 200 300 400 500
Duración de Lluvia (min)

Gráfico HI - N° 2: Curvas IDF de la subcuenca Condoray.

45.00
1000 años 500 años 200 años
40.00 100 años 50 años
Intensidad de Precipitación (mm/hr)

35.00
0.17
I= 28.49 T
30.00 .0.526
Donde: t
25.00 T= años
t= minutos
20.00

15.00

10.00

5.00

0.00
0 100 200 300 400 500
Duración de Lluvia (min)

Gráfico HI - N° 3: Curvas IDF de la subcuenca Jacayita.

En el Anexo III se muestra el cuadro de lluvias e intensidades máximas para las
subcuencas de interés.

3.7 Análisis de Tormentas de Diseño

Con la información de precipitación obtenida para cada uno de los componentes del modelo
hidrológico (subcuencas), se procedió a realizar la generación de los hietogramas de
diseño; es decir, la distribución temporal de las precipitaciones máximas en 24 horas.

Para este fin, se utiliza como referencias el patrón de comportamiento de la lluvia

determinado por el Soil Conservation Service (SCS) del U.S. Departament of Agriculture
(1986). Dicha institución desarrolló hietogramas sintéticos adimensionales para utilizar en
diferentes zonas con duraciones de tormentas de 6 a 24 horas. Existen cuatro tipos de
tormentas de 24 horas de duración, llamadas Tipo I, IA, II y III, representadas como factores
de distribución de tormenta, a partir de las cuales se elaboraron los hietogramas
correspondientes.

Los tipos I y IA corresponden principalmente al clima de la vertiente del Pacifico con

inviernos húmedos y veranos secos.

Por lo antes indicado, el tipo de tormenta típica seleccionada para la distribución de la

precipitación en 24 horas de la zona de estudio, es de Tipo I. En el Cuadro HI - N° 7 se
muestra los hietogramas adimensionales de tormentas típicas 24 horas, que fue obtenida
en nuestro caso para cada hora.

La distribución de la lluvia o hietograma de diseño, necesaria para el modelo hidrológico,

se obtendrán de la multiplicación de la precipitación de diseño obtenida anteriormente por
los respectivos factores adimensionales de la distribución adimensional SCS Tipo I. Este
procedimiento se realizó para las subcuencas y los resultados se muestran en el Anexo IV.

Cuadro HI - N° 7: Hietogramas adimensionales de Tormenta Típicas 24 horas.

FACTORES ADIMENSIONALES PT/P24
Tiempo (minutos)
Tipo I Tipo IA Tipo II Tipo III
0 0 0 0 0
60 0.017 0.022 0.011 0.010
120 0.035 0.051 0.023 0.020
180 0.055 0.083 0.035 0.032
240 0.076 0.116 0.048 0.043
300 0.099 0.156 0.064 0.057
 FACTORES ADIMENSIONALES PT/P24
Tiempo (minutos)
Tipo I Tipo IA Tipo II Tipo III
360 0.125 0.204 0.080 0.072
420 0.156 0.268 0.100 0.089
480 0.194 0.425 0.120 0.115
540 0.254 0.520 0.147 0.148
600 0.515 0.577 0.181 0.189
660 0.624 0.623 0.236 0.250
720 0.682 0.664 0.663 0.600
780 0.728 0.701 0.776 0.751
840 0.766 0.736 0.825 0.811
900 0.799 0.769 0.856 0.848
960 0.83 0.800 0.881 0.886
1020 0.857 0.830 0.903 0.904
1080 0.882 0.858 0.922 0.922
1140 0.905 0.884 0.938 0.939
1200 0.926 0.908 0.953 0.957
1260 0.946 0.946 0.965 0.968
1320 0.965 0.965 0.977 0.979
1380 0.983 0.983 0.989 0.989

3.8 Parámetros del software HEC-HMS 4.3

El modelamiento hidrológico del HEC HMS 4.3, determina la escorrentía originada por las
lluvias de diseño. Este programa fue creado por el Cuerpo de Ingenieros del Ejército de los
Estados Unidos de Norte América, para generar hidrogramas de avenidas causadas por la
precipitación.

Al programa se requiere tres datos de entrada:

 Modelo de cuenca o base: modelo de pérdida de agua, modelo de transformación,

modelo de flujo de base y modelo de tránsito de avenida.

 Modelo meteorológico: donde se establecen los parámetros meteorológicos a

utilizarse por las cuencas, en este caso las tormentas de diseño para cada
subcuenca.

 Especificaciones de control

El programa HEC HMS reporta los hidrogramas de salida en forma tabular y gráfica, lo que
permite visualizar el desarrollo de los eventos; los hidrogramas de salida se prolongan por
un periodo superior a las tres horas, debido a la extensión de las cuencas.
3.8.1 Modelo de Pérdida de Agua

Para el Modelo de pérdida de agua se utilizará el modelo del Número de Curva SCS, en el
cual se debe conocer por lo menos los siguientes tipos de datos3 tipos de datos:

• El número de curva CN

• La abstracción inicial (mm)

• Porcentaje de suelo impermeable en la cuenca (%)

El número de curva CN, varía entre 0 y 100 y se aplica para condiciones antecedentes de
humedad; normales (AMC II), para condiciones seca (AMC I) o condiciones húmedas (AMC
III). Los números de curva equivalentes pueden calcularse por:

4.2𝐶𝑁(𝐼𝐼) 23𝐶𝑁(𝐼𝐼)
𝐶𝑁(𝐼) = ; 𝐶𝑁(𝐼𝐼𝐼) =
10 − 0.058𝐶𝑁(𝐼𝐼) 10 + 0.13𝐶𝑁(𝐼𝐼)

El grupo de Antecedente de Humedad (AMC) se determina según la siguiente tabla:

Cuadro HI - N° 8: Valores de Número

LLUVIA ANTECEDENTE TOTAL DE 5 DIAS (PULG)

GRUPO AMC ESTACION
ESTACION ACTIVA
INACTIVA
I Menor que 0.5 Menor que 1.4
II 0.5 a 1.1 1.4 a 2.1
III Mayor que 1.1 Mayor que 2.1

Los Números de Curva han sido tabulados por el SCS con base en el tipo de suelo y el uso
de la tierra. Se definen cuatro grupos de suelos:

 Grupo A: Arena profunda, suelos profundos depositados por el viento, limos

agregados.
 Grupo B: Suelos pocos profundos depositados por el viento, marga arenosa.
 Grupo C: Margas arcillosas, margas arenosas poco profundas, suelos con
bajo contenido orgánico y suelos con altos contenidos de arcilla.
 Grupo D: Suelos que se expanden significativamente cuando se mojan,
arcillas altamente pláticas y ciertos suelos salinos.
Los valores de CN para varias condiciones de uso de la tierra en estos tipos de suelos se
muestran en el siguiente cuadro:

Cuadro HI - N° 9: Valores de Numero de Curva (CN)

GRUPO HIDROLÓGICO DEL

DESCRIPCIÓN DEL USO DE LA TIERRA
SUELO
A B C D
Tierra cultivada:
Sin tratamiento de conservación 72 81 88 91
Con tratamiento de conservación 62 71 78 81
Pastizales:
Condiciones pobres 68 79 86 89
Condiciones óptimas 39 61 74 80
Vegas de ríos: condiciones óptimas 30 58 71 78
Bosques: 45 66 77 83
Troncos delgados, cubierta pobre, sin hierbas, cubierta buena. 25 55 70 77
Área abiertas, césped, parques, campos de golf, cementerios, etc.
Óptimas condiciones: cubierta de pasto en el 75% o más condiciones 39 61 74 80
aceptables cubierta de pasto en el 50 al 75% 49 69 79 84
Áreas comerciales de negocios (85% impermeables) 89 92 94 95
Distritos Industriales (72% impermeables) 81 88 91 93

Residencial:
Tamaño promedio del lote Porcentaje promedio impermeable
1/8 acre o menos 65
77 85 90 92
1/4 acre 38
61 75 83 87
1/3 acre 30
57 72 81 86
1/2 acre 25
54 70 80 85
1 acre 20
51 68 79 84
Parqueadores pavimentados, techos, accesos, etc. 98 98 98 98
Calles y carreteras:
Pavimentados con cunetas y alcantarillados. 98 98 98 98
Grava 76 85 89 91
Tierra 72 82 87 89

Respecto al CN se eligió el tipo II del grupo de antecedente de humedad (AMC). El CN (II)

se obtuvo en base al ANA, en su publicación de Información Geocientífica 2015-05-31, ha
establecido un visor oficial del número de curva en Perú donde se tiene un GRID con
tamaño de celdas de 90x90m. Asimismo, en las cuencas con eventos de “Huaycos”
registrados se estableció condiciones de humedad tipo III.

Según esta información y con ayuda del programa de información geográfica Qgis se
obtuvo CN para las sub cuencas que se presenta a continuación:

Cuadro HI - N° 10: Datos de CN para las sub cuencas.

Subcuenca CN
San Jerónimo 69
Condoray 85
Jacayita 85

3.8.2 Modelo de Transformación

Una vez construido el modelo hidrológico, es de necesidad la utilización de hidrogramas de

tormenta.

El hidrograma se obtiene de forma directa a través de las estadísticas de data hidrométrica

y en el caso de carecer de esta información y contar sólo con estaciones de medición de
precipitaciones, se utilizan las características morfométricas de las cuencas.

La aplicación de los hidrogramas unitarios sintéticos (HUS) en zonas donde se carece de

información hidrométrica es esencial, ya que realiza una estimación de los flujos de la
cuenca y utiliza variables de fácil medición, transformándose en una herramienta simple y
potente la cuenca; siendo de ésta, información importante, la pendiente, cotas, áreas, la
longitud del curso principal, entre otras, con fines de construir el hidrograma indirectamente.

El modelo hidrológico es de tipo agregado y está dividido en subcuencas no mayores de

617 km2; por lo que, para determinar las descargas de los ríos y quebradas se ha utilizado
el método: Hidrograma Triangular Unitario Sintético de U.S. SOIL CONSERVATION
SERVICE (SCS), ya que los datos de entrada son de disposición directa y fácil estimación
y se opta por hacer uso del hidrograma unitario triangular en el modelo HEC-HMS.

Es conveniente anotar, que otro de los métodos conocidos es el HUS de Snyder, el cual
genera dudas en el caso de las tormentas de alta intensidad, principalmente por los
mayores volúmenes de crecida obtenidos, por lo que el uso de este hidrograma no se
recomienda en intensidades mayores a 100 mm/h. Finalmente, los HUS y la metodología
del número de curva se transforman en una herramienta confiable y de fácil aplicación en
lugares donde sólo se maneja información pluvial (Referencia: Pizarro – Tapia et al. Análisis
Comparativo de tres hidrogramas unitarios sintéticos en dos cuencas mediterráneas de
Chile, región del Maule. Revista Chapingo).

3.8.3 Módulo de Tránsito de Avenidas

El módulo de tránsito de avenidas es un procedimiento matemático para predecir el cambio

en magnitud, velocidad y forma de una onda de flujo, en función del tiempo (Hidrograma de
Avenida), en uno o más puntos a lo largo de un curso de agua. Para el presente estudio se
utilizará el método Muskingum, en el cual se deben conocer 3 tipos de datos:

 K: Medida del tiempo de transporte de una onda de un punto a otro, expresándose

en unidades de tiempo (horas), esta puede ser calculada como el tiempo de retardo
del tramo evaluado. Por tanto, se definirá según la siguiente ecuación:

∆𝑥
𝐾 = 0.18 ∗ ( ) ˆ(0.76)
𝑖ˆ(0.25)
Dónde:

 ∆x : Es la distancia Máxima en Km
 i : Es la pendiente máxima de la subcuenca.

Fórmula Alternativa

𝐾 = 0.6 ∗ 𝑇𝑐

 X € [ 0.1, …..0.5 ]: Depende de la pendiente media del tramo, a más pendiente más
alta será la X . En general, se usan valores entre 0.2 para las pendientes más bajas
y 0.35 para las pendientes más altas.

En los Anexos V se muestran los valores de parámetros de tránsito de caudales

utilizados en el modelo.

3.8.4 Tiempo de Concentración y Retardo

También denominado tiempo de respuesta o de equilibrio. Llamas (1993) lo define como el

tiempo requerido para que, durante un aguacero uniforme, se alcance el estado
estacionario; es decir, el tiempo necesario para que todo el sistema (toda la cuenca)
contribuya eficazmente a la generación de flujo en el desagüe. Se atribuye muy
comúnmente el tiempo de concentración al tiempo que tarda una partícula de agua caída
en el punto de la cuenca más alejado (según el recorrido de drenaje) del desagüe en llegar
a éste.

Por tener el concepto de tiempo de concentración una cierta base física, han sido
numerosos los autores que han obtenido formulaciones del mismo, a partir de
características morfológicas geométricas de la cuenca. A continuación, se muestran las
fórmulas utilizadas en el presente estudio:

a. Fórmula de Kirpich

L0.77
Tc = 0.000325
S 0.385

Dónde:
Tc : Tiempo de concentración (h)
L : Longitud del cauce (m)
S : Pendiente de la cuenca (m/m)

b. Formula de U.S. Corps of Engineers

L0.76
Tc = 0.3 ( )
S 0.19

Dónde:
Tc : Tiempo de concentración (h)
L : Longitud del cauce (km)
S : Pendiente de la Cuenca (m/m)

c. Fórmula de Temes
0.75
𝐿
Tc = 0.3 ( 0.25 )
S

Tc : Tiempo de concentración en horas

L : Longitud del cauce principal en kilómetros
S : Diferencia de cotas sobre L en porcentaje.

d. Fórmula de Chow
 𝐿 0.64
Tc = 0.00505 ( 0.5 )
S

Tc : Tiempo de concentración en horas

L : Longitud del cauce principal en kilómetros
S : Diferencia de cotas sobre L en porcentaje.

e. RR. Rowe
0.385
𝐿3
Tc = 0.86 ( )
∆ℎ

Tc : Tiempo de concentración en horas

L : Longitud del cauce principal en kilómetros
∆h : Desnivel de L en metros.

El tiempo de concentración para el modelo precipitación – escorrentía, será determinado

por el promedio aritmético de los principales formulas mostradas (Kirpich, RR.
Rowe,Temes, Chow y Soil Conservation Service ), lo cual ha sido tomado como criterio al
existir diferencias menores del 10 %, considerando que nos da una mejor representatividad
de este parámetro para nuestra área de estudio. . En el Anexo V.1 se observa los valores
promedio del tiempo de concentración para las subcuencas.

El tiempo de retardo que es el valor que ingresa al modelo del programa HEC – HMS se
calcula por la siguiente formula:

Tr = 0.6 ∗ (Tc)

Dónde:
Tr : Tiempo de retardo (min)
Tc : Tiempo de Concentración (min)

3.8.5 Modelo Meteorológico

Respecto al modelo meteorológico, se emplearán los hietogramas obtenidos de las

tormentas SCS, considerándose una tormenta tipo I y las precipitaciones máximas para
diferentes periodos de retorno de las subcuencas, como se mencionó en el en el Capítulo
3.7. La distribución de lluvia de las subcuencas se muestra en el Anexo IV, a fin de
determinar las avenidas máximas que se producen para un periodo de retorno de 50, 100,
200,500 años.

3.8.6 Especificaciones de Control

En las especificaciones de control, se tendrá en cuenta que el tiempo de modelamiento será

mayor a la duración del Hietograma de diseño.

3.9 Resultados

Los resultados de caudales picos en 24 horas de las subcuencas de interés se muestran

en el Anexo VI. En el Anexo VII, se muestran las fotografías de la visita de campo en las
zonas de interés.
 Quebrada San Jerónimo

Para esta quebrada, los periodos de retorno de 50,100,200 y 500 años, se tiene como
resultado el siguiente hidrograma y el cuadro resumen de los caudales máximos que se
muestran a continuación:

Figura HI - N° 1: Hidrograma de la subcuenca San Jeronimo para los

periodos de retorno de 50,100,200 y 500 años

6
50 Años 100 Años 200 Años 500 Años

4
Q (M3/S)

0
12:00
15:00
18:00
21:00

12:00
15:00
18:00
21:00

12:00
15:00
18:00
21:00
0:00
3:00
6:00
9:00

0:00
3:00
6:00
9:00

0:00
3:00
6:00
9:00

0:00

TIEMPO (HR)
 Cuadro HI - N° 11: Resumen de Caudales Máximos para los periodos de
retorno de 50,100,200 y 500 años de la subcuenca San Jerónimo.

Periodo de Retorno (años) QPico(m3/s)

50 0,7
100 1,4
200 2,2
500 5,6

 Quebrada Condoray

Para esta quebrada, los periodos de retorno de 50,100,200 y 500 años, se tiene como
resultado el siguiente hidrograma y el cuadro resumen de los caudales máximos que se
muestran a continuación:
Figura HI - N° 2: Hidrograma de la subcuenca Condoray para los
periodos de retorno de 50,100,200 y 500 años

7
50 Años 100 Años 200 Años 500 Años
6

4
Q (M3/S)

0
0:00
3:00
6:00
9:00
12:00
15:00
18:00
21:00
0:00
3:00
6:00
9:00
12:00
15:00
18:00
21:00
0:00
3:00
6:00
9:00
12:00
15:00
18:00
21:00
0:00

Tiempo (hr)

Cuadro HI - N° 12: Resumen de Caudales Máximos para los periodos de retorno de

50,100,200 y 500 años de la subcuenca Condoray

Periodo de Retorno (años) QPico(m3/s)

50 0,90
100 1,70
200 3,10
500 6,30
  Quebrada Jacayita

Para esta quebrada, los periodos de retorno de 50,100,200 y 500 años, se tiene como
resultado el siguiente hidrograma y el cuadro resumen de los caudales máximos que se
muestran a continuación:
Figura HI - N° 3: Hidrograma de la subcuenca Jacayita para los periodos de retorno
de 50,100,200 y 500 años

1.2
50 Años 100 Años 200 Años 500 Años

0.8
Q (M3/S)

0.6

0.4

0.2

0
10:00
12:00
14:00
16:00
18:00
20:00
22:00

10:00
12:00
14:00
16:00
18:00
20:00
22:00
0:00
2:00
4:00
6:00
8:00

0:00
2:00
4:00
6:00
8:00

0:00
TIEMPO (HR)

Cuadro HI - N° 13: Resumen de Caudales Máximos para los periodos de retorno de

50,100,200 y 500 años de la subcuenca Jacayita

Periodo de Retorno (años) QPico(m3/s)

50 0,1
100 0,2
200 0,4
500 1,0

Es necesario indicar que los caudales obtenidos para estas subcuencas reflejan el caudal
líquido y no los caudales reales que se presentan como “flujo de escombros”, en donde se
debe incorporar las concentraciones de flujo hiperconcentrado.

Así mismo, se puede observar en las figuras anteriores, que cada una de ellas muestra una
distribución espacial diferente; lo que refleja a su vez las distintas características
geomorfológicas de cada subcuenca.
4 INFLUENCIA DEL EVENTO DE EL NIÑO

El pronóstico del Fenómeno de El Niño se basa en el análisis del estado del clima en el
Océano Pacífico Tropical. Para realizar este pronóstico, existe actualmente el Comité
Multisectorial encargado del Estudio Nacional del Fenómeno El Niño (ENFEN), que forma
parte de un esfuerzo regional en el marco del Programa "Estudio Regional del Fenómeno
El Niño" de la Comisión Permanente del Pacífico Sur CPPS compuesta por Perú, Chile,
Ecuador y Colombia.

El ENFEN sistemáticamente monitorea y estima el arribo de ondas oceánicas Kelvin cálidas

y frías a nuestras costas; esto, con el fin de prever las variaciones de la temperatura de la
superficie del mar, su impacto en las temperaturas del aire y, de ser el caso, su relación con
las precipitaciones cuando coincide con el período lluvioso.

Como indica el ENFEN, las ondas Kelvin ecuatoriales son las posibles precursoras de un
evento El Niño, ya sea en el Pacífico ecuatorial o en la costa del Perú.

Para su detección el ENFEN realiza el monitoreo de las ondas Kelvin ecuatoriales mediante
el uso de dos tipos de observaciones oceanográficas: 1) el nivel del mar mediante altimetría
satelital como JASON-2, y 2) la temperatura sub superficial para determinar la profundidad
de la termoclina (capa oceánica que separa, en forma relativamente abrupta, el agua fría
profunda y el agua más cálida cercana a la superficie).

Como se indica, se presentan dos tipos del Fenómeno El Niño y los cuales son observados
en distintas zonas del Océano Pacifico; los valores de tres meses anómalos de temperatura
de la superficie del mar en la Región “Niño 1+2” se denomina Niño Costero y los tres meses
de anomalía de la temperatura superficial del mar en la región “Niño 3.4” se denomina Niño
Global. Ambos eventos se registran mediante los índices ONI (Índice de Niño Oceánico) e
ICEN (Índice Costero de El Niño).

Bajo ese contexto, para la evaluación del impacto en la intensidad de los eventos de
precipitación en el área de estudio, se plantea realizar:

1) Análisis estadístico entre los eventos hidrometeorológicos observados en la cuenca

con registros de episodios de El Niño ONI (Oceanic Niño Index) disponibles por la National
Oceanic and Atmospheric Administration (NOAA), de esta manera correlacionar
intensidades del fenómeno del Niño pronosticadas por el ENFEN a magnitudes de caudal
para un periodo de retorno asociado.
2) Identificar eventos históricos de precipitación extrema en la cuenca del rio Cañete
para observar su distribución espacial, de esta manera se evaluará el comportamiento en
la parte baja, media y alta de la cuenca asociado a eventos de El Niño. Para esto se
utilizarán registros históricos de precipitación grillada (espacialmente distribuida) global
como GPM, CHIRPS o MSWEP.

De los gráficos de Indice Costero de El Niño y superponiendo los caudales máximos diarios
reportados por el SENAMHI, se observa que existe relación entre los eventos de El Niño
Global de los años 2015, 1998 y 1982, con Índice de Niño Oceánico de 2.2, con caudales
máximos suscitados, los mismos que se aproximan a caudales de 50 años de periodo de
retorno.

En el año 2002 (según manifiesta la población en Cañete, porque no hay registro de

caudales) y 2017, se presentaron “Huaycos” en la zona media de la Cuenca y coincide con
los eventos de El Niño Costero, con Índices Costero de El Niño de 1.16 y 1.71 y en el río
Cañete se produjeron caudales máximos que se aproximan a periodos de retorno de 5
años. Con las coincidencias encontradas podemos indicar que al suscitarse alertar de El
Niño Costero de nivel Moderado y cerca al límite de nivel Fuerte, se suscitan grandes
precipitaciones en la zona media y baja de la Cuenca y se activan las cuencas con riesgo
a generar “huaycos”, ya que estas subcuencas presentan características geológicas
aparentes para que estos fenómenos ocurran; sin embargo, esto no influye en grandes
magnitudes al caudal del río Cañete.
Figura HI - N° 4: Índice Costero El Niño (ICEN)

Fuente: ICEN, 2019.

Figura HI - N° 5: Índice Niño Oceánico (ONI)

Fuente: ONI, 2019.

 Figura HI - N° 6 Índice Niño Oceánico vs Caudal Máximo Diario Socsi (m3/s):

Fuente: ONI, Huang et al. 2017 – SENAMHI - SIGT, 2019.

5. FLUJO DE ESCOMBROS – HUAYCOS

5.1 Aspectos Generales

Cuando el transporte de sólidos es tan grande y la densidad de la mezcla se eleva por

encima de la densidad del agua, se introduce los conceptos de fluidos “No-Newtonianos”.
Desde el punto de vista hidráulico, a esta rama se conoce como “Hidráulica Torrencial” o
“Flujo de Escombros”.

Según Takahashi (1980) “Los flujos de lodo y escombros, son flujos antiguos conformados
por una mezcla viscosa y sedimentos de todos los tamaños con acumulaciones de cantos
rodados que se vuelcan en el frente de la onda y forman lóbulos, detrás del cual siguen los
granos más finos”.

La composición de los flujos de escombros varía debido a que éste, es un flujo de múltiples
fases de granos desordenados con aire y agua mezclados entre sí. Entre los principales
componentes tenemos: el agua, la matriz fina (partículas de arcilla) y los granos gruesos.
El tamaño de partículas varía desde materiales finos hasta cantos rodados suspendidos en
una matriz viscosa.

Son muchos los factores y condiciones que generan este evento; por ejemplo, las fuertes
precipitaciones y la disponibilidad de material, que el material esté saturado y una
concentración de partículas tal que se comporte como un fluido No Newtoniano (relación
de esfuerzo-deformación no es lineal), tendiendo a un comportamiento reológico de un
fluido plástico-viscoso. Este material conformado por sedimentos proviene de la gran
cubierta cuaternaria y deslizamientos de suelo en laderas y material depositado en las
quebradas. Los tipos de flujos que se presentan a continuación:

- Flujo hiperconcentrado: Concentración de sólidos entre 20% y 65%

- Flujo de detritos: Concentración de sólidos sobre 65%, finos menor del 10% de
concentración de sólidos.
- Flujo de lodos: Concentración de sólidos mayor a 50%

Se debe tomar en cuenta lo complejo que es establecer los parámetros reológicos de cada
fluido; por lo tanto, se hace más confiable la adopción de la información de campo como
son las huellas máximas observadas de los últimos eventos de “huaycos”.
Se debe tomar en consideración que las secciones de flujo tipo Mud Flow, alcanzan
mayores secciones hidráulicas, pero poca capacidad erosiva. Por el contrario, los flujos tipo
Debris Flow, presentan mayores velocidades y también mayor capacidad erosiva, lo que
deberá tomarse en cuenta para los diseños de las respectivas obras de conducción y/o
protección.

5.2 Ocurrencia de Flujo de Escombros en Cuenca Cañete

Según reportes de instituciones e información de la población asentada en la cuenca del

río Cañete, específicamente en la cuenca media, se han presentado estos eventos con
daños a las viviendas en el año 2002 y en mayor medida el año 2017; sin embargo, se
presentan con una frecuencia de casi dos años en menor medida. Estos eventos coinciden
con el llamado “Niño Costero” de los años antes mencionados y que se reportan con los
Índice Costero El Niño que se presenta en el Capítulo 4.

En el evento más reciente del año 2017, se puede mencionar que se presentaron distintos
flujos, respecto al material de arrastre y la concentración volumétrica de sólidos. Estos son
evaluados tomando en cuenta los restos de material y que se observan en las secciones
mostradas en el presente capítulo.

La ocurrencia de eventos y este tipo de flujo no solo dependen de la intensidad de las

precipitaciones o caudales, sino que se deben considerar otras condiciones como:
antecedentes de precipitación, características de la cuenca y del cauce, pendiente,
existencia de material para ser removido, condiciones climáticas locales y en estos casos
eventos de El Niño (ligado a la ocurrencia de lluvias e intensidades).

5.3 Metodologías para la Estimación de Caudal del Flujo de Escombros

Para la estimación de los caudales máximos de los flujos de escombros existen 2

metodologías más usadas, la de JICA (Japan International Cooperation Agency) y la de
O’Brien (2000).
a. Metodología JICA – Takahashi (1981)

La metodología usada por JICA se basa en los desarrollos de Takahashi. El cual representa
el volumen de lodo y escombros es igual a:

𝑪𝒅
𝑽𝒔 = 𝑽𝒂 [ ]
𝑪∗ − 𝑪𝒅

Donde:
Vs: Volumen de detritos en m3
Va: Volumen de agua en m3
Cd: Densidad volumétrica del flujo de lodo y piedras
C*: Densidad volumétrica de sedimentos del lecho
Takahashi (1978, 1981) determina el caudal máximo de la corriente de detritos mediante la
siguiente expresión:

Según (Montserrat, 2005) el caudal máximo de flujos detríticos propuesta por Takahashi
(1978, 1981), la cual liga el caudal máximo del flujo detrítico con el caudal máximo líquido
de la descarga. El Gobierno de Japón (Ministerio de la Construcción) propone este método
para flujos detríticos donde el material proviene principalmente de la erosión del lecho
(Ayala, 1996).

𝑪∗
𝑸𝒅 = 𝑸𝟎 𝑷 = 𝑸𝟎 (𝑪∗
−𝑪𝒅 )

Donde:
Qd : Representa el caudal máximo del flujo detrítico
Q0 : El caudal máximo líquido
P : Corresponde a un coeficiente de amplificación,
Cd : Concentración volumétrica de sólidos y
C* : La concentración máxima de los sólidos (o concentración de
empaquetamiento). Generalmente C* toma valores del orden del 65%. Cd se estima según
la siguiente ecuación:

𝝆 𝒕𝒂𝒏𝜶
𝑪𝒅 =
𝝆𝒔 − 𝝆 (𝒕𝒂𝒏𝝓 − 𝒕𝒂𝒏𝜶)

Donde:
α corresponde al ángulo de inclinación del lecho, y ϕ al ángulo de fricción interna del material
depositado (entre 30º y 40º). Cd es denominada por Takahashi (1981) como la
concentración de equilibrio del flujo, la cual se alcanza en condiciones estacionarias.
Resultados obtenidos mediante esta ecuación muestran una buena correlación con datos
experimentales (Takahashi, 1981).

b. Metodología O’Brien (2000)

Otra metodología empleada es la propuesta por O’Brien (2000) donde el hidrograma de

flujo de detritos se estima a partir del hidrograma de flujo de agua multiplicado por un factor
que es una función de la concentración volumétrica de detritos. La ecuación utilizada por
O’Brien es la siguiente:

Qt = BF*Ql

𝟏
𝑩𝑭 =
𝟏 − 𝑪𝒗

Donde
Qt: Caudal total de flujo de detritos (m3/s)
Ql: Caudal líquido (m3/s)
BF: Factor (Bulking Factor)
C v: Concentración volumétrica de detritus.
Según la referencia de O’Brien (2000), para flujos de detritos el valor de Cv varía entre 0.20
(para flujos bajos) y 0.45 (para flujos altos).

c. VanDine (1996)

Se puede relacionar la velocidad del flujo con la pendiente, geometría del cauce, viscosidad
dinámica, peso específico y altura de escurrimiento mediante la fórmula de Poiseuille.

𝛾 sin 𝛼 ℎ2
𝑈=
𝑙𝑣

Donde;
U : velocidad de flujo
tan 𝛼 : pendiente del cauce
h : altura de escurrimiento
v : viscosidad dinámica
l : constante concerniente a la forma de la sección transversal del cauce (l=3 para
cauces amplios; l=8 para canales semicirculares)

5.4 Caudales Máximos Totales de Flujo de Escombros

Como se sabe los eventos denominados “huaycos” corresponden a un flujo no newtoniano

lo cual es muy complejo establecer a la fecha fórmula, ecuaciones y/o modelos que puedan
caracterizar adecuadamente los mismos. Inclusive en cada quebrada se pueden presentar
varios tipos de flujos dependiendo de condiciones como: frecuencia y antecedentes de
precipitación, características de la cuenca y el cauce, existencia de material para ser
removido, condiciones climáticas locales y en estos casos eventos de El Niño (ligado a la
ocurrencia de lluvias e intensidades).

La hidráulica de los flujos de escombros en la cual el transporte de sólidos es tan grande

que los conceptos de hidráulica fluvial no son aplicables. La fase sólida influye en el flujo,
es decir, no puede separarse el flujo de agua, por un lado, con sus ecuaciones de
movimiento, y el transporte sólido por otro.

En el siguiente cuadro se muestra los valores estimados para el cálculo de caudales

máximos estimados con las metodologías de Takahashi y O`Brien, que fueron corroborados
con marcas máximas levantadas topográficamente en campo, además de registros en video
publicados en las redes virtuales de algunos eventos.

Cuadro HI - N° 14: Caudales Máximos de Flujo de Escombros (m3/s) en las

Quebradas de Interés.

Periodo
de Qda San Jerónimo Qda Condoray Qda Jacayita
Retorno
Qt Qt Qt
Q Qt (m3/s) Q Qt (m3/s) Q Qt (m3/s)
3 (m3/s) (m 3
/s) (m 3
/s)
Años (m /s) (según (m3/s) (según (m3/s) (según
(según (según (según
Líquido Takahashi) Líquido Takahashi) Líquido Takahashi)
O´Brien) O´Brien) O´Brien)
50 0.7 18.2 1.8 0.9 19.5 2.3 0.1 6.5 0.3
100 1.4 20.7 3.5 1.7 22.0 4.3 0.2 7.2 0.6
200 2.2 23.0 5.5 3.1 24.5 7.8 0.4 8.1 1.1
500 5.6 26.4 14.0 6.3 27.6 15.8 1.0 9.3 2.8
Como se puede observar en el cuadro anterior, ambas metodologías arrojan valores
disímiles. Los caudales obtenidos con la metodología Takahashi, arroja valores más
cercanos a los verificados hidráulicamente mediante marcas máximas de lodos registrados
topográficamente en cada quebrada de interés. Según la metodología O’Brien, se tiende a
subestimar los caudales de escombros que se presentarían.

QUEBRADA SAN JERÓNIMO (LUNAHUANÁ)

La quebrada de San Jerónimo, que cruza el distrito de Lunahuaná, ha tenido un

comportamiento mixto, primero se comportó como un flujo de lodos; ya que eran flujos
lentos y con alta composición de arcillas y arenas finas; para luego transformarse en flujo
de escombros con gran capacidad erosiva. En campo se pueden observar aún las huellas
de lodo que son superiores a las huellas del flujo erosivo, lo que se podrá apreciar en las
vistas siguientes.

Figura HI - N° 7: Vista de Aguas arriba y aguas debajo de la quebrada San Jeronimo

Vista aguas arriba desde carretera Cañete - Yauyos. Véase Vista lateral aguas abajo desde carretera Cañete -
las marcas de lodos de último evento en marzo 2017. Yauyos. Véase marca de lodos del último evento marzo
2017.
 Cuadro HI - N° 15: Caudal Máximo Flujo de Lodos (Liquido + Sólido)

Periodo Caudal
Qt Total Volumen
de Líquido
3
(m /s) Total Area (m2) Velocidad
Retorno (m3/s)
(MMC) (m/s)
50** 0.7 18.2 0.43 25.8 0.7
100 1.4 20.7 0.48 27.3 0.8
200 2.2 23.0 0.54 28.5 0.8
500 5.6 26.4 0.62 30.3 0.9
** El evento de referencia de la evaluación de campo corresponde para Tr = 50 años.

El cálculo de caudal total de escombro corresponde a la metodología Takahashi.

Asimismo, estos caudales consideran un fluido con concentraciones altas en sólidos y finos,
donde predominan las arcillas y limos; se ha estimado un flujo de bajas velocidades (Mud
Flow), que corresponde un área hidráulica de 25.8 m2. Seguidamente se producen los flujos
con mayor capacidad erosiva y altas velocidades, en este caso alrededor de 4.7 m/s; el
mismo que debe considerarse como parámetro para el diseño de estructuras de
encauzamiento.
QUEBRADA CONDORAY

La quebrada de Condoray que cruza el distrito de Lunahuaná, han tenido un

comportamiento mixto en el evento de marzo del 2017, primero se comportó como un flujo
de lodos; ya que eran flujos lentos y con alta composición de arcillas y arenas finas; para
luego transformarse en flujo de escombros con gran capacidad erosiva. En campo se
pueden observar aún las huellas de lodo que son superiores a las huellas del flujo erosivo,
lo que se podrá apreciar en las vistas siguientes.

Figura HI - N° 8: Vista de Aguas arriba y aguas debajo de la quebrada Condoray

Vista aguas arriba desde carretera Cañete - Vista aguas abajo carretera Cañete - Yauyos.
Yauyos. Véase marcas de lodos de último evento Véase marcas de lodos de último evento en
en marzo 2017. marzo 2017.
 Cuadro HI - N° 16: Caudal Máximo Flujo de Lodos (Liquido + Sólido)

Periodo Caudal
Qt Total Volumen
de Líquido
(m3/s) Total Area Velocidad
Retorno (m3/s)
(MMC) (m2) (m/s)
50** 0.9 19.5 0.46 24.7 0.79
100 1.7 22.0 0.52 26.0 0.85
200 3.1 24.5 0.58 27.3 0.90
500 6.3 27.6 0.67 28.9 0.96

** El evento de referencia de la evaluación de campo corresponde para Tr = 50 años.

El cálculo de caudal total de escombro corresponde a la metodología Takahashi.

Asimismo, se debe remarcar que estos caudales consideran un primer tipo de flujo con
concentraciones altas en sólidos y finos, donde predominan las arcillas y limos. En este
caso, según lo observado de los rastros del evento de marzo del 2017, concentración en
volumen se ha asumido en 60% y es el flujo que alcanza los mayores tirantes, pero no así
mayores caudales debido a sus bajas velocidades (Mud Flow); alcanzando un área en la
sección de control de 24.7 m2.

Seguidamente se produjeron los flujos con mayor capacidad erosiva y altas velocidades,
por sobre los 4.0 m/s, (Debris Flow); por lo que se recomienda tomar estos parámetros para
los diseños de estructuras de encauzamiento y otros.
QUEBRADA JACAYITA

La quebrada de Jacayita, que cruza el distrito de Lunahuaná, ha tenido un comportamiento

mixto, primero se comportó como un flujo de lodos; ya que eran flujos lentos y con alta
composición de arcillas y arenas finas; para luego transformarse en flujo de escombros con
gran capacidad erosiva. En campo se pueden observar aún las huellas de lodo, luego del
cruce de la carretera hasta aguas abajo y son superiores a las huellas del flujo erosivo (los
mismos que se observan en tramos superiores.

Figura HI - N° 9: Vista de Aguas debajo de la quebrada Jacayita

Vista aguas bajo carretera Cañete - Yauyos. Se puede

observar un canal encauzado de tierra en dirección al río
Cañete.
 Cuadro HI - N° 17: Caudal Máximo Flujo de Lodos (Liquido + Sólido)

Periodo Caudal
Qt Total Volumen
de Líquido
3
(m /s) Total Area (m2) Velocidad
Retorno (m3/s)
(MMC) (m/s)
50** 0.1 6.5 0.14 9.8 0.7
100 0.2 7.2 0.15 10.2 0.7
200 0.4 8.1 0.16 10.6 0.8
500 1.0 9.3 0.18 11.2 0.8
** El evento de referencia de la evaluación de campo corresponde para Tr = 50 años.

El cálculo de caudal total de escombro corresponde a la metodología Takahashi.

Asimismo, estos caudales consideran un fluido con concentraciones altas en sólidos y finos,
donde predominan las arcillas y limos; se ha estimado un flujo de bajas velocidades (Mud
Flow), que corresponde un área hidráulica de 9.8 m2. Seguidamente se producen los flujos
con mayor capacidad erosiva y altas velocidades, en este caso alrededor de 3.0 m/s; el
mismo que debe considerarse como parámetro para el diseño de estructuras de
encauzamiento.
 ANEXOS

ANEXO I: Registro histórico de

las Precipitaciones Máximas
en 24 horas (Fuente: SENMAHI)
ANEXO II: Pruebas de
Kolmogorov – Smirnov para
las estaciones