“INSTALACION DE DEFENSA RIBEREÑA EN EL RIO SANTA TRAMO
ASENTAMIENTO HUMANO LAS FLORES - QUEBRADA BRUJA ARMANAN,
DISTRITO DE CARHUAZ, PROVINCIA DE CARHUAZ – ANCASH”
ESTUDIO HIDROLÓGICO
I. INTRODUCCION
Dentro del Estudio de Defensa Ribereña, con la finalidad de calcular la altura de la
obra de protección, es necesario realizar el Estudio Hidrológico, de la cuenca de
aporte, que nos permitirá conocer el caudal máximo que discurre en donde se
proyecta construir dicha obra.
El Proyecto “Instalación de Defensa Ribereña en el Río Santa Tramo Asentamiento
Humano Las Flores – Quebrada Bruja Armanan, Distrito de Carhuaz, Provincia de
Carhuaz – Ancash” está ubicado en el distrito de Carhuaz, de la provincia de
Carhuaz, del departamento de Ancash.
En el estudio ha sido evaluada la descarga máxima, parámetro necesario para el
diseño de la obra de protección.
El estudio hidrológico se ha realizado en dos etapas: trabajo de campo y el
procesamiento de datos.
Al final del estudio se presentan las conclusiones y recomendaciones relacionadas
a las variables evaluadas.
II. UBICACIÓN
2.1 Ubicación Política
Departamento : Ancash
Provincia : Carhuaz
Distrito : Carhuaz
2.2 Ubicación Hidrográfica
Cuenca : Río Santa
1
2.3 Hidrografía de la Sub Cuenca en Estudio
La cuenca del río Santa está ubicada casi en su totalidad en el departamento
de Ancash entre los meridianos 78°40’ y 77°10’ y los paralelos 8°00’ y 10°15’
N y tiene un área de drenaje de 1,700 Km2. Su curso principal recorre de sus
a norte el valle formado por la Cordillera Negra, por la margen izquierda y la
Cordillera Blanca por la margen derecha, cuyas cumbres sirven al mismo
tiempo de límites de la cuenca. El río recorre por este valle un tramo de unos
150 Km., bajando en elevación desde más de 4,000 m.s.n.m. a unos 20 Km.
al norte de Cátac, hasta los 2,500 m.s.n.m. próximo a Caraz. A partir de este
lugar el río, con fuerte pendiente, corre a través del Cañón del Pato y
después de unos 120 Km. al norte de Chimbote.
El 85% del área de la cuenca del río Santa corresponde a la cuenca húmeda,
zona en la cual la precipitación pluvial representa un aporte efectivo al
escurrimiento superficial. Los deshielos de los nevados del flanco occidental
de la Cordillera Blanca inciden en la escorrentía, contribuyendo a mantener
una considerable descarga de estiaje.
Esta Cordillera, con una longitud de 170 Km., empieza en el nevado
Rajutuna, al sur, en las vecindades de la laguna Conococha y termina en el
nevado de Champará al norte. El pico más importante de esta cadena de
nevados es el Huascarán con una altitud de 6,768 m.s.n.m. Otros nevados
dignos de mencionar, como el Caullaraju, Huarapasca, Pongos, Yanamarey,
Cashan, Huantsan, Palcaraju, Hualcán, Huandoy, Yanapaccha, Pucahirca,
Alpamayo y Santa Cruz, cubren un área de aproximadamente 550 Km2.
También destacan en esta cuenca las numerosas lagunas originadas por los
deshielos de estos nevados. Las más importantes son: Conococha (donde
nace el río Santa), Querococha, Angrishcocha, Llanganuco, Parón,
Cullicocha, Safuna y Pucacocha, emplazadas sobre los 4,000 m.s.n.m.
El río Santa recibe el aporte de numerosos afluentes, siendo los más
importantes los de la margen derecha (provenientes de los nevados y
lagunas de la Cordillera Blanca). Entre ellos Pachacoto, Yanayacu, Olleros,
Quilcay, Marcará, Llanganuco, Parón, Colcas, Cedros, Quitaracsa, Manta y
Tablachaca.
2
A lo largo del valle, conocido como Callejón de Huaylas, se encuentran varias
pequeñas ciudades y poblaciones, siendo la ciudad de Huaraz la capital del
Departamento de Ancash y la más poblada. Todo el valle, desde Cátac hasta
Caraz está extensamente cultivado.
Características del Régimen Fluvial de la Cuenca
La cuenca del río Santa presenta un período bien definido de aguas altas,
entre los meses de enero y abril, un período seco entre junio y octubre, y
período de transición en los meses restantes. De la distribución anual de
caudales, se observa una fluctuación menor de los caudales en las
estaciones localizadas en los cursos de agua afluentes al río Santa por la
margen derecha, cuyas nacientes están en la Cordillera Blanca, debido al
efecto regulador de las lagunas existentes en esos ríos.
En la lámina Nº 01, se muestra la Cuenca de Aporte del río Santa, hasta la
zona donde se proyecta la defensa ribereña.
III. INFORMACION BASICA:
Las fuentes principales de información para el estudio hidrológico del Proyecto,
son el tipo cartográfico y meteorológico.
3.1 Cartografía:
- Plano catastral (1/100,000), recopilado del Instituto Geográfico Nacional
(IGN), hoja 21i Chiquián, 20i – Recuay, 20h – Huaraz, 19i – Huari, 19h -
Carhuaz.
IV. CARACTERISTICAS GEOMORFOLOGICAS:
Las principales características geomorfológicas de la Cuenca de Aporte, hasta
donde se proyecta la estructura de protección, son:
4.1 Area: El área total de drenaje de la Cuenca de Aporte hasta la zona donde se
proyecta la obra de protección es 3229.43 km2.
4.2 Perímetro: El perímetro de la Cuenca de Aporte, es 237.31 km.
3
4.3 Coeficiente de Compacidad: La Cuenca de Aporte, tiene un coeficiente de
compacidad de 1.58, el cual indica que la unidad hidrográfica es ligeramente
alargada, donde la respuesta hidrológica es rápida con respecto a una
cuenca de forma alargada.
4.4 Elevaciones de la Cuenca: La elevación mínima de la Cuenca en estudio es
2602.50 m.s.n.m. y su elevación máxima 5680 m.s.n.m.
4.5 Calculo de la Pendiente del Río Santa
La velocidad del escurrimiento de las corrientes de agua depende de la
pendiente de sus canales. A mayor pendiente mayor velocidad.
El perfil longitudinal del Río Santa va desde la cota mínima (2602.50
m.s.n.m.) a la cota máxima (4200 m.s.n.m). Recorriendo una longitud total de
110.572 km. La pendiente media del río es de 1.48%.
V. ESTUDIO DE MAXIMAS AVENIDAS
Como no se tiene estación de aforo en el lugar donde se ubicará la obra de
defensa ribereña; el cálculo de la descarga de diseño se realizó mediante el
método de la Regionalización de las Descargas Máximas Instantáneas Anuales.
Para la Regionalización de las Descargas Máximas, se consideraron 9
Subcuencas del río Santa que cuentan con estación de aforo y la Estación “La
Balsa” que controla las descargas de la cuenca del río Santa en el Callejón de
Huaylas. La ubicación de las subcuencas y de la estación “La Balsa” se indica en
el Cuadro Nº 01.
Previamente a la regionalización, se evaluó el comportamiento temporal de las
descargas de los datos de las nueve subcuencas y de la Estación “La Balsa”,
mediante modelos probabilísticos, donde la secuencia en el tiempo de las
variables no interesa y estas siguen la ley de una determinada distribución de
frecuencias o la ley de un determinado modelo probabilístico.
4
Cuadro Nº 01: Ubicación de las Estaciones de Aforo de las Sub-Cuencas
del río Santa
Altitud Latitud Longitud Area
Estación Río msnm Sur Oeste Km2
Recreta Santa 4018 10° 02’ 77° 19.2’ 290.0
Pachacoto Pachacoto 3745 9° 51’ 77° 24’ 198
Querococha Yanayacu 4037 9° 43.5’ 77° 19.8’ 63
Chancos Chancos 2872 9° 19.2’ 77° 34.7’ 210
Llanganuco Llanganuco 3916 9° 4.2’ 77° 39’ 86.9
Parón Parón 4100 9° 00’ 77° 41.1’ 48
Colcas Colcas 2048 8° 55.2’ 77° 50.5’ 226
Cedros Los Cedros 1878 8° 52.2’ 77° 49.3’ 112
Quitaracsa Quitaracsa 1480 8° 47.8’ 77° 51’ 383
La Balsa Santa 1861 8° 52.5’ 77° 49.5’ 5124
Fuente: Glaciares y Recursos Hídricos en la Cuenca del Río Santa (2003)
IRD-INRENA – UGRH – SENAMHI – EGENOR
5.1 Modelos Probabilísticos
Los modelos probabilísticos se caracterizan por ser independientes en el
tiempo, es decir son modelos donde la secuencia en el tiempo de las
variables no interesa y estas siguen la ley de una determinada distribución
de frecuencias o la ley de un determinado modelo probabilístico.
Se consideró emplear los siguientes modelos probabilísticos: Gumbel, Log
Normal de 2 parámetros, Log Normal de 3 parámetros, Log Gumbel,
Gamma de 2 parámetros y Gamma de 3 parámetros.
5.2 Selección del Modelo Probabilístico Adecuado
Para la selección del modelo probabilístico adecuado se emplean los
métodos gráfico y estadístico.
Para el presente trabajo se utilizó el método estadístico de Kolmogorov-
Smirnov.
5
La prueba de ajuste se realizó para los datos de caudales máximos de las 9
Subcuencas del río Santa que cuentan con estación de aforo y la Estación
“La Balsa” (Cuadro Nº 01), mediante la prueba estadística Smirnov -
Kolmogorov, aplicando el programa HIDROESTA, para las distribuciones
probabilísticas: Logarítmico Normal de 2 parámetros, Logarítmico Normal de
3 parámetros, Gamma de 2 parámetros, Gamma de 3 parámetros, Gumbel
y Log Gumbel.
Cuadro Nº 03: Resumen de la Prueba Kolomogorov-Smirnov
¿SE AJUSTA EL MODELO?
ESTACION Log Normal 2 Log Normal 3 Gamma 2 Gamma 3 Log
Gumbel
parámetros parámetros parámetros parámetros Gumbel
CEDROS SI SI SI SI SI SI
CHANCOS SI SI SI SI SI SI
COLCAS SI SI SI SI SI SI
LLANGANUCO SI SI SI SI SI SI
PACHACOTO SI SI SI NO SI SI
PARON SI SI SI SI SI SI
QUEROCOCHA SI SI SI NO SI SI
QUITARACSA SI NO SI NO SI SI
RECRETA SI SI SI SI SI SI
LA BALSA SI SI SI SI SI SI
Del Cuadro Nº 03, se tiene que los datos de caudales máximos de las
estaciones consideradas en el estudio, se ajustan a las distribuciones
probabilísticas analizadas, a excepción de la estaciones Pachacoto,
Querococha y Quitaracsa que no se ajustan a la distribución Gamma 3,
asimismo de la estación Quitaracsa que no se ajusta a la distribución Log
Normal de 3 parámetros.
Como se observa en el Cuadro Nº 03, según la prueba de Kolmogorov-
Smirnov, se aceptarían las distribuciones probabilísticas Log Normal de 2
parámetros, Gamma de 2 parámetros, Gumbel y Log Gumbel, por cuanto el
delta calculado es menor al delta tabular, todas consideradas dentro de un
nivel de significancia 0.05 .
La distribución probabilística más conocida y usada es la de Gumbel,
principalmente para eventos hidrológicos extremos, como en el presente caso
se trata de caudales máximos; por lo que para el presente estudio, se asume
que los datos de caudales máximos de las 9 Estaciones de las Sub Cuencas
en Estudio y de la Estación de La Balsa, se distribuyen a la distribución
probabilística Gumbel.
6
VI. Regionalización de las Descargas Máximas Instantáneas
Tomando como modelo probabilístico adecuado el modelo Gumbel, se
calcularon las descargas máximas para diferentes períodos de retorno teórico,
para todas las subcuencas analizadas y la Estación La Balsa. Los resultados
se muestran en el Cuadro Nº 04.
Cuadro Nº 04: Descargas Máximas para diferentes
Períodos de Retorno (m3/s)
PERIODO DE RETORNO (AÑOS)
ESTACION
2 5 10 20 25 50 100 250 500
CEDROS 11.08 14.85 17.35 19.75 20.51 22.85 25.18 28.24 30.55
CHANCOS 32.47 39.66 44.42 48.99 50.44 54.90 59.33 65.16 69.56
COLCAS 20.32 26.43 30.48 34.37 35.60 39.40 43.17 48.13 51.17
LLANGANUCO 6.24 7.57 8.45 9.29 9.56 10.38 11.20 12.28 13.09
PACHACOTO 26.03 34.88 40.75 46.37 48.16 53.65 59.11 66.29 71.72
PARON 2.85 3.37 3.71 4.04 4.14 4.46 4.78 5.20 5.52
QUEROCOCHA 7.27 9.01 10.17 11.27 11.62 12.70 13.78 15.19 16.25
QUITARACSA 54.70 68.45 77.56 86.29 89.06 97.60 106.07 117.23 125.65
RECRETA 29.19 41.04 48.89 56.41 58.80 66.16 73.46 83.07 90.33
LA BALSA 547.17 802.73 971.93 1134.23 1185.71 1344.31 1501.74 1709.01 1865.52
La regionalización de las descargas máximas instantáneas anuales se
realizó en función al Area de drenaje (Cuadro Nº 05) para diferentes
modelos matemáticos, obteniéndose como modelo matemático
regionalizado, el modelo matemático lineal, que resulta con el más alto
coeficiente de correlación.
Cuadro Nº 05: Descargas Máximas Anuales y Area de Drenaje
AREA Q
ESTACION 2
Km m3/s
RECRETA 290 31.17
PACHACOTO 198 28.08
QUEROCOCHA 63 7.60
CHANCOS 210 35.09
LLANGANUCO 86.9 6.49
PARON 48 2.92
COLCAS 226 21.45
CEDROS 112 11.78
LA BALSA 5124 594.68
QUITARACSA 383 57.26
7
Ecuación Potencial: Qmáxp 0.0661* A1.094
r = 0.981
Ecuación Lineal: Qmáxp 1.533 0.1159 A (*)
r = 0.9994
Donde:
Qmáx : Promedio de las descargas máximas anuales (m3/sg)
A : Area de la cuenca colectora en Km2
r : Coeficiente de correlación
Luego, se dividieron las descargas máximas ajustadas al modelo entre
el promedio de las descargas máximas anuales, obteniéndose valores
adimensionales, que se indican en el Cuadro Nº 06.
Cuadro Nº 06: Valores Adimensionales de Descargas Máximas
PERIODO DE RETORNO (AÑOS)
ESTACION
2 5 10 20 25 50 100 250 500
CEDROS 0.94 1.26 1.47 1.68 1.74 1.94 2.14 2.40 2.59
CHANCOS 0.93 1.13 1.27 1.40 1.44 1.56 1.69 1.86 1.98
COLCAS 0.95 1.23 1.42 1.60 1.66 1.84 2.01 2.24 2.39
LLANGANUCO 0.96 1.17 1.30 1.43 1.47 1.60 1.73 1.89 2.02
PACHACOTO 0.93 1.24 1.45 1.65 1.72 1.91 2.11 2.36 2.55
PARON 0.98 1.15 1.27 1.38 1.42 1.53 1.64 1.78 1.89
QUEROCOCHA 0.96 1.19 1.34 1.48 1.53 1.67 1.81 2.00 2.14
QUITARACSA 0.96 1.20 1.35 1.51 1.56 1.70 1.85 2.05 2.19
RECRETA 0.94 1.32 1.57 1.81 1.89 2.12 2.36 2.66 2.90
LA BALSA 0.92 1.35 1.63 1.91 1.99 2.26 2.53 2.87 3.14
MEDIA 0.947 1.224 1.407 1.585 1.642 1.813 1.987 2.211 2.379
DESV 0.018 0.071 0.124 0.178 0.194 0.246 0.298 0.362 0.414
Para luego realizar el análisis de regresión entre los valores
adimensionales y los períodos de retorno transformados en logaritmo
natural, obteniéndose un coeficiente de correlación altamente
significativo (r=0.9993). Dicho modelo matemático, que representa el
modelo regional para estimar las descargas máximas anuales en la
Cuenca del río Santa, es el siguiente:
8
Qmáx, T Qmáxp(0.8037 0.2562Ln(T )) (**)
r = 0.9993
Donde:
Qmáx, T : Descarga máxima para el período de retorno T (m3/sg)
Qmáxp : Promedio de la descarga máxima instantánea anual
(m3/sg)
Ln : Logaritmo natural
T : Período de retorno en años
VII. Calculo de Caudales Máximos para diferentes Períodos de Retorno
Con las ecuaciones (*) y (**) se calcularon los caudales máximos anuales
instantáneos para la cuenca de aporte hasta donde se proyecta la oba de
defensa ribereña, tomando en cuenta los períodos de retorno y el área de
drenaje de la cuenca de aporte A = 3229.43 km2.
En el Cuadro Nº 07, se muestra los resultados de caudales de diseño para
diferentes períodos de retorno y probabilidades de excedencia y no
excedencia.
Cuadro Nº 07: Caudales Máximos para diferentes Periodos de
Retorno - Cuenca de aporte al proyecto
Período de Probabilidad de Probabilidad de Caudal
Retorno (años) No Excedencia (%) Excedencia (%) Diseño (m3/s)
5 80 20 456.71
10 90 10 523.53
20 95 5 590.35
25 96 4 611.86
50 98 2 678.68
100 99 1 745.49
200 99.5 0.5 812.31
500 99.8 0.2 900.64
El grado de asociación entre el caudal de diseño y el período de retorno, se
expresa mediante la siguiente relación matemática:
9
Qd 68.628Ln(T ) 214.69 R2 = 1.00
Qd : Caudal de diseño, m3/s
T : Tiempo de Retorno, años
R2 : Coeficiente de Determinación
Mediante esta ecuación matemática, se calculan los caudales de diseño para
diferentes períodos de retorno; que es importante para dimensionar diferentes
obras hidráulicas, dentro de la Cuenca de aporte.
En el Gráfico Nº 01, se muestra la relación de Caudales de Diseño para
diferentes períodos de retorno, correspondiente a la Cuenca de aporte.
Gráfico Nº 01: Caudales de Diseño para diferentes Períodos de Retorno
Cuenca de Aporte Río Santa
VIII. Período de retorno y Riesgo de Falla
El diseño para el control del agua normalmente se basa en eventos extremos
de corta duración. Los factores determinantes en la selección del evento de
diseño son el costo y el grado de seguridad que se desea dar a la estructura
hidráulica. La magnitud óptima para el diseño es aquella que equilibra las
consideraciones de costo y seguridad.
10
Una estructura para el control de agua puede fallar, si la magnitud
correspondiente al período de retorno de diseño Tr se excede durante la vida
útil de la estructura. Este Riesgo hidrológico natural de falla se calcula
mediante:
R 1 1 P( X T )
n
1
T
1 (1 R)1 / n
Donde :
R : Riesgo permisible que el suceso ocurra en cualquier período
N : Vida útil de la estructura
P : Probabilidad de que un evento, ocurra por lo menos una vez en años
Tr : Período de Retorno
Considerando una vida útil de la estructura de 25 años, un riesgo de falla de
0.22, el tiempo de retorno de la avenida máxima resulta igual a Tr = 100 años.
Se toma un riesgo de falla del 22%, debido a que adyacente a la margen
derecha del río Santa, se ubican viviendas pertenecientes a las Urbanizaciones
de San Gabino y Las Flores, por lo que con fines de protección se ha tomado
dicho porcentaje de riesgo de falla.
Por tanto, el caudal de diseño para un período de retorno de 100 años es igual
a 745.49 m3/s.
IX. CONCLUSIONES Y RECOMENDACIONES
9.1 CONCLUSIONES:
9.1.1 El área de la Cuenca que aportará con caudal máximo, es de 3229.43
Km2.
9.1.2 El caudal máximo para el lugar donde se ubica el proyecto, con un
período de retorno igual a 100 años es igual a 745.49 m3/sg.
9.1.3 El comportamiento hidrológico de la Cuenca en estudio, se caracteriza
por presentar caudales máximos en época de lluvias, en época de
ausencia de lluvias, la Cuenca aporta con mínimo caudal.
11
9.1.4 Se ha tomado un caudal máximo para 100 años de período de retorno,
porque generalmente el cauce de los ríos se forma con el efecto de
avenidas; además, con la finalidad de proteger las viviendas ubicadas
en la margen derecha del tramo en estudio, se debe tomar en cuenta
un riesgo de falla, que en el presente estudio se toma igual al 22%.
9.2 RECOMENDACIONES:
9.2.1 Se recomienda realizar el diseño hidráulico, para el caudal máximo igual
a 745.49 m3/sg.
9.2.2 Realizar trabajos de protección, tanto aguas arriba y aguas debajo de la
estructura.
9.2.3 La construcción de la obra de defensa debe realizarse en época de
estiaje.
12