Hidrologia e Hidraulica Del Proyecto
Hidrologia e Hidraulica Del Proyecto
Hidrologia e Hidraulica Del Proyecto
INDICE
Relación de Tablas
Relación de Cuadros
Relación de Gráficos
Para realizar el análisis hidrológico se utilizó como base una carta de modelo digital d de elevación
(DEM), proporcionado por la Instalación Satelital de Alaska (ASF), de la Administración Nacional
de Aeronáutica y el Espacio de los Estados Unidos de Norte América (NASA). Los productos ALOS
PALSAR han sido descargados EN formato Geo TIFF compatible con GIS, que pertenecen al
proyecto de Terreno Radiométricamente Corregido (RTC). Los productos RTC utilizados para el
proyecto tienen una resolución de 12.5 y 30 m para las áreas de mapeo indicadas.
También se han utilizado cartas nacionales en formato vectorial de cobertura vegetal, tipo de
suelo y número de curva, proporcionados por las instituciones gubernamentales en sus portales
de servicio.
2.1.2. DESCRIPCIÓN
Actualmente en la zona del proyecto existen estructuras de protección en la ribera del río como
muros de contención de gaviones y tramos de muros de enrocado, estas estructuras no satisfacen
con la protección integral por lo que se debe realizar obras de ampliación y un estudio completo
de la dinámica del tramo del río Torres en la parte baja de la sub cuenca del mismo nombre.
El área del proyecto se encuentra en la zona de la vertiente del Amazonas, en la Intercuenca Alto
Marañon V, más precisamente en la sub cuenca del Río Torres.
Las condiciones climatológicas en la zona de estudio son influenciadas por la altitud, presión
atmosférica y las corrientes de viento que viajan desde el atlántico, esta última condición trae en
su viaje nubes desde el atlántico sur, hacia la zona montañosa de la cordillera de los andes, donde
existe la mayor ocurrencia de precipitación a lo largo del año. Para años donde las zonas de
convergencia intertropicales (ZCIT) y la del atlántico sur (ZCAS), desarrollan episodios
convectivos sobre el continente, como en el caso del ZCAS, la ocurrencia de precipitación es
mayor sobre la Intercuenca alto Marañon V, ocurriendo precipitaciones orográficas en la parte
alta de la cuenca Río Torres.
El río torres que es el río principal de la sub cuenca del mismo nombre nace a una altitud de 4829
msnm, y es tributario del río Vizcarra en la zona urbana del distrito de Huallanca, esto ocurre a
3534 msnm en el sector de Santa Rosa.
La zona centro del Perú ha sufrido en múltiples oportunidades eventos extremos como los
conocidos fenómenos del niño, por lo que se han realizado diversos estudios con el fin de
contrarrestar los efectos destructivos producidos por la creciente de los ríos y quebradas.
• Los estudios previos no consideran el fenómeno del niño (FEN), como evento de diseño.
2.3. HIDROLOGÍA
Los datos de precipitaciones máximas en 24 horas fueron obtenidos de los registros de la estación
del SENAMHI, Estación Huallanca, ubicada en el distrito de Huallanca, provincia de Bolognesi, en
la región Ancash. Se escogió esta estación por encontrarse dentro de la sub cuenca en estudió
(Sub Cuenca del río Torres), la cual incluyen registros del año 1964 al 1976.
Precipitación
HUALLANCA Máxima en 24 286173.64 8906898.69 3815.00 1964-1976
horas
Max
Año Ene Feb Mar Abr May Jun Jul Ago Sep Oct Nov Dic
Anual
1964 21.0 19.0 24.0 11.0 42.0 1.0 13.0 13.0 15.0 22.0 23.9 15.0 42.0
1965 13.0 18.0 18.0 18.0 17.0 5.0 11.0 7.0 13.0 10.0 13.0 17.0 18.0
1966 17.0 35.0 28.0 17.7 18.5 9.0 0.0 8.5 16.5 17.6 11.0 24.0 35.0
1967 20.0 20.5 21.0 20.0 13.5 4.0 0.5 15.0 17.5 24.0 10.5 12.5 24.0
1968 19.0 21.0 25.0 29.0 13.0 7.5 10.0 16.0 15.0 22.0 14.5 20.0 29.0
1969 16.0 35.8 10.0 25.8 9.4 10.0 4.4 8.4 15.0 29.0 32.2 23.8 35.8
1970 21.8 24.2 23.0 18.2 14.0 13.8 12.2 11.4 15.0 25.4 15.0 15.3 25.4
1971 20.2 48.2 21.0 21.0 15.5 7.0 5.5 28.4 9.0 8.0 14.2 33.0 48.2
1972 21.2 20.4 21.2 18.0 12.5 0.0 4.0 14.0 11.0 16.0 28.0 18.0 28.0
1973 15.5 25.0 16.0 22.0 17.0 2.9 14.3 7.0 15.0 19.0 19.3 15.5 25.0
1974 13.0 21.0 19.0 34.0 9.7 8.0 32.0 12.8 9.0 25.0 18.0 17.0 34.0
1975 31.5 18.8 20.5 14.3 11.0 10.0 5.0 12.0 9.0 14.0 28.0 20.0 31.5
1976 26.5 19.1 21.0 25.0 12.8 9.1 3.0 9.0 9.2 23.4 23.4 13.0 26.5
Nº Datos 13 13 13 13 13 13 13 13 13 13 13 13 156.0
Promedio 19.7 25.1 20.6 21.1 15.8 6.7 8.8 12.5 13.0 19.6 19.3 18.8 16.8
Desv.Est. 5.2 9.1 4.4 6.2 8.3 4.0 8.4 5.6 3.1 6.3 7.1 5.6 1.9
Máxima 31.5 48.2 28.0 34.0 42.0 13.8 32.0 28.4 17.5 29.0 32.2 33.0 48.2
Mínima 13.0 18.0 10.0 11.0 9.4 0.0 0.0 7.0 9.0 8.0 10.5 12.5 0.0
Fuente: ANA y SENAMHI.
50.0
40.0
30.0
20.0
10.0
0.0
1962 1964 1966 1968 1970 1972 1974 1976 1978
AÑOS
Los datos dudosos (outliers) son puntos de la información que se alejan significativamente de la
tendencia de la información restante. La retención o eliminación de estos datos puede afectar
significativamente la magnitud de los parámetros estadísticos calculados para la información,
especialmente en muestras pequeñas. Los procedimientos para tratar los datos dudosos
requieren un criterio que involucra consideraciones matemáticas e hidrológicas. De acuerdo con
el Water Resources Council (1981), si la simetría de estación es mayor que +0.4, se consideran
primero las pruebas para detectar datos dudosos altos; si la asimetría de estación es menor que
-0.4, primero se consideran pruebas para detectar datos dudosos bajos. Cuando la asimetría de
la estación esta entre +-0.4, deben aplicarse pruebas para detectar datos dudosos altos y bajos
antes de eliminar cualquier dato dudoso del conjunto de datos.
Las siguientes ecuaciones de frecuencia han sido utilizadas para detectar datos dudosos altos y
bajos respectivamente.
𝑦𝐻 = 𝑦̅ + 𝐾𝑛 𝑠𝑦
𝑦𝐿 = 𝑦̅ + 𝐾𝑛 𝑠𝑦
PP (mm) PP (mm)
PARAMETROS
Observado Max Aceptable
MAX. PP24 48.20 52.70
MIN. PP24 18.00 17.09
Los sistemas hidrológicos son afectados algunas veces por eventos extremos, tales como
tormentas severas, crecientes y sequias. La magnitud de un evento extremo esta inversamente
relacionado con su frecuencia de ocurrencia, es decir, eventos muy severos ocurren con menor
frecuencia que eventos más moderados. El objetivo del análisis de frecuencia de información
hidrológica es relacionar la magnitud de los eventos extremos con su frecuencia de ocurrencia
mediante el uso de distribuciones de probabilidad.
• Distribución Normal
• Distribución Log Normal
• Distribución Gumbell
• Distribución Pearson Tipo III
• Distribución Log Pearson Tipo III
60.00 61.51
57.35
Precipitacion en mm
53.05
50.00
47.03 y = 6.1373ln(x) + 32.868
42.01 R² = 0.9985
40.00
30.00
20.00
10.00
0.00
1 10 100 1000
Se determina que la frecuencia de mejor ajuste es la que presenta la distribución Log Normal y
se calculan las precipitaciones máximas para distintos tiempos de retorno.
Los datos obtenidos, se han multiplicado por un factor de corrección como lo sugiere la
Organización Meteorológica Mundial (OMM).
Las estaciones de lluvia ubicadas en la zona, no cuentan con registros pluviográficos que permitan
obtener las intensidades máximas. Sin embargo, estas pueden calcularse a partir de las
precipitaciones máximas. Para el estudio se utilizó la metodología de Dick Peschke, la cual se
presenta la siguiente ecuación:
𝑑 0.25
𝑃𝑑 = 𝑃24ℎ ( )
1440
Donde:
𝑑 = 𝐷𝑢𝑟𝑎𝑐𝑖ó𝑛 (𝑚𝑖𝑛𝑢𝑡𝑜𝑠)
CURVA INTENSIDAD-DURACION-FRECUENCIA
ESTACIÓN HUALLANCA
140.0
5
120.0
10
Intensidad de precipitación (mm/h)
100.0 25
50
80.0 100
140
60.0 500
40.0
20.0
0.0
0.0 50.0 100.0 150.0 200.0 250.0 300.0 350.0 400.0
Duración (min)
Para calcular el tiempo de concentración critico se utilizó la fórmula de kirpich, Kirpich Modificado,
Giandiott, Temez, US. Corps of Engineers, y US. Hathaway.
Formula de Kirpich:
𝑡𝑐 = 0.000325𝐿0.77 𝑆 −0.385
0.385
𝐿3
𝑡𝑐 = 0.0195 ∗ ( )
𝐻
0.030𝐿0.76
𝑡𝑐 =
𝑆 0.19
Donde:
𝑚
𝑆 = 𝑃𝑒𝑛𝑑𝑖𝑒𝑛𝑡𝑒 𝑎 𝑙𝑜 𝑙𝑎𝑟𝑔𝑜 𝑑𝑒𝑙 𝑐𝑎𝑢𝑐𝑒 ( )
𝑚
𝑛 = 𝐹𝑎𝑐𝑡𝑜𝑟𝑒 𝑑𝑒 𝑟𝑢𝑔𝑜𝑠𝑖𝑑𝑎𝑑
Debido a que existe una relación inversa entre la duración de una tormenta y su intensidad, se
asume que la duración crítica es igual al tiempo de concentración.
A continuación, se presentan los valores de tiempo de concentración para las diferentes fórmulas
propuestas.
TIEMPO DE CONCENTRACIÓN
FORMULA KIRPICH KIRPICH (Modificado) CALIFORNIA TEMEZ
Tc 1.547 1.546 1.547 1.538
min 92.799 92.763 92.799 92.301
Tc Final 92.79 min
Se utilizo el mapa de número de curva generado por la Autoridad Nacional del Agua (ANA). La
metodología utilizada se basa en la utilización de herramientas de los sistemas de información
geográfica (SIG), los cuales permiten trabajar con datos diversos como mapas de cobertura
vegetal, suelos y modelos digitales de elevación (MDE).
• Mapa temático de cobertura vegetal y uso, fuente Ministerio del Ambiente, escala
1/100000.
• Mapa temático de suelo, fuente: FAO y UNESCO a una escala de 1:5 000 000
• Modelo Digital de Elevación (MDE), fuente SRTM resolución espacial 90 m.
A continuación, se presenta un mapa del número de curva par la sub cuenca del Río Torres.
𝑅 = 1 − (1 − 1/T)𝑛
Donde:
𝑇 ∶ 𝑃𝑒𝑟𝑖𝑜𝑑𝑜 𝑑𝑒 𝑟𝑒𝑡𝑜𝑟𝑛𝑜
En la siguiente tabla se presenta el valor de T para varios riesgos permisibles R y para la vida útil
N de la obra.
RIESGO
VIDA ÚTIL DE LAS OBRAS (n años)
ADMINSIBLE
R 1 2 3 5 10 20 25 50 100 200
0.01 100 199 299 498 995 1990 2488 4975 9950 19900
A partir de las consideraciones que presenta la tabla con respecto al riesgo admisible y el tiempo
útil de la obra se ha considerado 25 % de riesgo admisible y 40 años de vida útil, por ser la obra
un sistema de defensa ribereña, teniendo de esta manera un valor de tiempo de retorno T=140
años.
Mientras que para cálculos de socavación se considera el tiempo de retorno, T=500 años.
Para calcular el caudal de diseño se utilizó el método de precipitación - escorrentía, puesto que
en la zona de proyecto no existen datos de caudales medidos. Los métodos utilizados utilizados
se describen a continuación:
El método racional según la fórmula propuesta por Témez, adaptada para las condiciones
climáticas de España. Y permite estimar la forma sencilla caudales punta en cuencas de drenaje
naturales con áreas menores a 770 Km2 y con tiempos de concentración Tc de entre 0.25 y 24
horas, la fórmula es la siguiente:
𝑄 = 0.278𝐶𝐼𝐴𝐾
Donde:
𝑚3
𝑄: 𝑑𝑒𝑠𝑐𝑎𝑟𝑔𝑎 𝑚á𝑥𝑖𝑚𝑎 𝑑𝑒 𝑑𝑖𝑠𝑒ñ𝑜 ( )
𝑠𝑒𝑔
𝑚𝑚
𝐼: 𝐼𝑛𝑡𝑒𝑛𝑠𝑖𝑑𝑎𝑑 𝑑𝑒 𝑝𝑟𝑒𝑐𝑖𝑝𝑖𝑡𝑎𝑐𝑖ó𝑛 𝑚á𝑥𝑖𝑚𝑎 ℎ𝑜𝑟𝑎𝑟𝑖𝑎 ( )
ℎ𝑟
𝐾: 𝐶𝑜𝑒𝑓𝑖𝑐𝑖𝑒𝑛𝑡𝑒 𝑑𝑒 𝑈𝑛𝑖𝑓𝑜𝑟𝑚𝑖𝑑𝑎𝑑.
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.
Para el modelo de perdida se utilizó el método de perdida de SCS Numero de Curva, mientras
que para el método de transformación de lluvia - escorrentía se utilizó los métodos que se
muestras a continuación:
Modelo meteorológico
Especificaciones de control
ESTACIÓN
TIEMPO DE
CAUDAL DE DISEÑO (m3/seg)
CONSENTRACIÓN
HU SCS HU SNYDER HU CLARK RACIONAL MODIFICADO
25
20
15
10
5
0
-500:00 02:24 04:48 07:12 09:36 12:00 14:24 16:48 19:12
TIEMPO (min)
En el anexo, se encuentran los cálculos realizados y las salidas del programa HEC-HMS.
Considerando las condiciones geomorfológicas del río Torres, se ha calculado el caudal sólido de
acarreo a partir de las siguientes ecuaciones, que involucran únicamente el caudal líquido y
pendiente longitudinal del lecho
𝑄𝑠 = 𝑘𝑆 𝑚 𝑄
Donde:
𝑄 = 𝐶𝑎𝑢𝑑𝑎𝑙 𝑙𝑖𝑞𝑢𝑖𝑑𝑜
𝑘 𝑦 𝑚 = 𝑐𝑜𝑒𝑓𝑖𝑐𝑖𝑒𝑛𝑡𝑒𝑠
𝐶 = 𝑘𝑆 𝑚
Donde:
𝑄𝑆
𝐶= ⁄𝑄 , 𝑒𝑠 𝑙𝑎 𝑐𝑜𝑛𝑐𝑒𝑡𝑟𝑎𝑐𝑖ó𝑛 𝑑𝑒 𝑠𝑒𝑑𝑖𝑚𝑒𝑛𝑡𝑜𝑠 𝑝𝑜𝑟 𝑢𝑛𝑖𝑑𝑎𝑑 𝑑𝑒 𝑐𝑎𝑢𝑑𝑎𝑙 𝑙í𝑞𝑢𝑖𝑑𝑜.
𝐶 = 2.5𝑆 1.6
Formulade Rickenman (1991), integrando datos propios Rickenmmann, 1990 y Smart y Jaeggi
(1983) y de Meyer – Peter y Müller (1948), ascendiendo en total a 252, recomienda la ecuación
para pendientes entre el 0.1% y 20%.
𝐶 = 1.5𝑆 1.5
𝐶 = 7.35𝑆 2
En el caso de la zona del proyecto se tiene cercano a las obras proyectadas, los aportes de dos
quebradas con evidencia de material de escombros.
Respecto a la descripción de la geodinámica externa se indica que durante los momentos máxima
intensidad de lluvia se producen huaycos por las quebradas SC-01 Y SC-02, ubicadas
respectivamente, en la margen derecha, cercana al inicio de las obras de protección.
Las dos quebradas son dinámicamente activas, con arrastre de materiales que se extienden por
la menor pendiente y energía, al depositar el torrente de escombros en su cauce, a la entrega
del río Torres.
Siendo que los modelos actuales de flujos de escombros (debris Flow) se basan en modelos
reológicos y parámetros empíricos, se considera que el régimen del flujo en las cercanías del pico
es de régimen laminar.
2.4. HIDRAÚLICA
El estudio de la hidráulica fluvial del río Torres parte del entendimiento de las características
geomorfológicas, topográficas, e hidráulicas en general del cauce y las zonas de interacción, por
lo que se ha desarrollado el estudio de la hidráulica fluvial para el tramo en estudio; de esta
manera se ha utilizado el caudal de diseño estimado en el capítulo de hidrología, para determinar
el transito distribuido de crecidas, a través del río; para tal fin se ha aplicado la dinámica de
Fluidos Computacional (CFD), de esta manera la simulación numérica operada en el software Iber
3.1 ha servido para realizar el modelado 2D de flujo.
Del modelo CFD se han determinado las velocidades medias y altura hidráulica del agua, este
último parámetro determina la altura de inundación, mientras que la velocidad servirá para
analizar condiciones de socavación. Estos parámetros permitieron realizar los diseños de las obras
de infraestructura en el medio fluvial que puedan reducir los riesgos por crecida del río.
La estimación del NAME (nivel de aguas máximas extraordinarias) se ha realizado sobre la base
de:
• Caudal máximo
• Secciones transversales del cause
• Perfil del cause
Para realizar el modelo hidráulico es necesario contar con un modelo topográfico del cauce del
río, llanuras de inundación, y demás zonas aledañas, por lo que se realizó el levantamiento
batimétrico y topográfico y a partir de ello se realizó un modelo de malla de 1 m en formato ACCII
para poder generar una superficie en el software Iber 3.1.
Mientras que la rugosidad del cauce y las zonas aledañas al sistema fluvial han sido determinadas
a partir de las consideraciones de la Guía para la selección del coeficiente de rugosidad de
Manning para Canales Naturales y Llanuras de Inundación del Servicio Geológico de los Estados
Unidos, (USGS). El USGS sugiere que los valores del coeficiente de rugosidad n, pueden ser
asignados para las condiciones que existen en el momento de un determinado evento de flujo,
para condiciones promedio en un rango de etapa, o para las condiciones previstas en el omento
de un evento futuro.
Los fatores más importantes que afectan a la selección de los valores de n son el tamaño de los
materiales que componen el lecho y las márgenes del cauce, como también la forma del canal
natural, es así que para el estudio se empleara el método de Cowan desarrollado en 1956, para
a cuál utilizaremos la siguiente ecuación.
𝑛 = (𝑛𝑏 + 𝑛1 + 𝑛2 + 𝑛3 + 𝑛4 )𝑚
𝑛𝑏 = 𝑐𝑜𝑒𝑓𝑖𝑐𝑖𝑒𝑛𝑡𝑒 𝑏𝑎𝑠𝑒.
𝑛2 = 𝑣𝑎𝑙𝑜𝑟 𝑝𝑎𝑟𝑎 𝑙𝑎𝑠 𝑣𝑎𝑟𝑖𝑎𝑐𝑖𝑜𝑛𝑒𝑠 𝑒𝑛 𝑓𝑜𝑟𝑚𝑎 𝑦 𝑡𝑎𝑚𝑎ñ𝑜 𝑑𝑒 𝑙𝑎 𝑠𝑒𝑐𝑐𝑖ó𝑛 𝑡𝑟𝑎𝑛𝑠𝑣𝑒𝑟𝑠𝑎𝑙 𝑑𝑒𝑙 𝑐𝑎𝑛𝑎𝑙.
El factor de corrección para la existencia de meandros ha sido considera con un valor de 1, por
la inexistencia del mismo, mientras que para determinar los valores de 𝑛0 se ha empleado la tabla
tomada de la Guía para la selección del coeficiente de rugosidad de Manning para Canales
Naturales y Llanuras de Inundación del Servicio Geológico de los Estados Unidos. En el cuadro
Nº se muestran los datos y los resultados obtenidos.
RIBERAS Y LLANURA
COONDICIÓN DEL CANAL VALORES CAUCE
DE INUNDACIÓN
n Básico n0 0.05 0.05
Grado de irregularidad n1 menor 0.005 0.005
Variación de la sección
n2 gradual 0 0
transversal
Efecto relativo de las
n3 insignificante 0 0
obstrucciones
Vegetación n4 baja 0.005 0.01
m - 0.06 0.065
Los valores considerados se han compatibilizado con la tabla de Ven Te Chow, que establece
rangos de valores entre 0.050 a 0.080.
2.4.2. CAUDALES
Los caudales considerados para el modelo son los obtenidos en el capítulo de hidrología de un
tiempo de concentración T=140 año y T=500 años. El caudal de 140 años es un caudal que
responde a las consideraciones de riesgo admisible y años de vida útil de la obra, mientras que
el caudal de 500 años responde a un criterio de estimación de socavación en áreas donde las
obras de protección tendrán la cimentación en contacto con el cauce; de esta manera se ha
utilizado los caudales para dos escenarios, donde el hidrograma es un dato de ingreso para el
modelo hidráulico, el cual describe un régimen no permanente.
Este análisis se realizó sobre las condiciones para un flujo crítico y subcrítico, de régimen no
permanente variado, que en condiciones de estrangulamiento del cauce en zonas donde existen
emplazamientos de la sub estructura de puentes se ha producido el flujo no permanente
rápidamente variado. La condición de entrada del caudal de diseño, como un hidrograma fue
analizado en un tiempo de 2500 segundo, tránsito a través del área de estudio que ha sido
utilizado con el criterio de enrutamiento hidrológico. Este caudal corresponde a un evento
extremo estimado, por lo que los parámetros del resultado solo se presentaran hipotéticamente
en dichas condiciones, y es la misma que ayudara a diseñar el sistema de protección ribereña.
Las características del modelo hidrodinámico, esta basado en las ecuaciones 2D de Saint Venant,
así como esquemas explícitos de volumen finito en mallas no estructuradas que corresponden a
la topografía del cause y llanuras de inundación del proyecto, que fueron obtenidas mediante un
estudio de topografía y batimetría.
2.4.4. RESULTADOS
A. VELOCIDADES MAXIMAS
Los valores de velocidad media estimados varían en un rango de 0 a 7.37 m/s, estando el
valor mas alto presente en regiones de estrangulamiento del cauce en el puente los
Quenuales, y puente Santa rosa, también existe una región de estrangulamiento del cauce a
100 metros de la Urbanización Chinllillin, este estrangulamiento se da de causa natural por
la presencia de una roca de tamaño considerable, donde la sección del cauce alcanza apenas
una solera de 2 metros de ancho, condición que permite la presencia de velocidades altas del
flujo. En la zona adyacente a la urbanización Los Quenuales, donde existe estructuras de
gaviones que encausan el flujo, presenta velocidades por encima del promedio total, en un
rango de 4.08 a 5.60 m/s. Otra región donde el régimen de velocidad media se encuentra
sobre el promedio total se encuentra a 200 metros del puente Santa Rosa, zona que a sido
encausada con muros de enrocado, en ambas márgenes.
VELOCIDADES MINIMAS.
Las velocidades mínimas están presentes a lo largo del cauce cuando el transito del
hidrograma se encuentra en sus niveles más bajos. Regiones como la urbanización los
Quenuales, Santa Rosa y en la desembocadura del río Torres, hacia el río Vizcarra tienen
velocidades del flujo que se encuentran en un rango de 1.05 y 2.09 m/s, mientras que, en el
inicio del área del proyecto, en la zona de la Urbanización Chinllillin, las velocidades mínimas
son permanentes, incluso en el punto más alto del hidrograma. La explicación para los datos
estimados se sustenta en la geometría del cauce, donde se aprecia que las características
geométricas son mayores con respecto a las otras zonas del área de estudio, también la
planicie de inundación que se encuentra a ambas márgenes, permitiendo el traslado del
caudal por los canales de inundación que se forman hacia aguas abajo, inundando toda el
área adyacente a la zona de la urbanización Chinllillin.
CALADO.
El modelo a mostrado valores de calado que se encuentran en un rango de 0 a 3.79 m, donde
el valor máximo se encuentra en la zona del puente Santa Rosa, lugar donde se produce el
estrangulamiento del cauce, las condiciones geométricas dan como resultado a un régimen
no permanente rápidamente variado, el mismo que permite el incremento del calado del
agua.
En la zona próxima aguas abajo a la urbanización Chinllillin donde se ubica la roca antes
descrita en el ítem de velocidades, el calado se incrementa hasta una altura de 3.79, mientas
que para zonas donde la capacidad hidráulica de la geometría del cauce es mayor, el calado
se reduce a niveles propios del caudal de formación.
B. Socavación general.
La socavación general se ha estimado aplicando el método de Lischtvan Lebediev.
• Lischtvan Lebediev
La máxima erosión del lecho ocasionada por el paso de las descargas máximas, e la
sección de interés, fue estimada haciendo uso de la expresión propuesta por Lischtvan –
Lebediev cuya expresión matemática es la siguiente.
𝑉𝑒 = 0.68𝛾𝑑 0.28 𝐻𝑆 𝑥+1
1
+1
𝑎𝐻0 5/3 𝑋
𝐻𝑆 = ( )
0.6𝛾𝑑 1.18
𝑄𝑑
𝑎=
𝐵𝑒 𝐻𝑚 5/3 𝜇
La profundidad de socavación será:
𝑑𝑆 = 𝐻𝑆 − 𝐻0
𝑄𝑑 = 𝑐𝑎𝑢𝑑𝑎𝑙 𝑑𝑒 𝑑𝑖𝑠𝑒ñ𝑜 𝑚2
𝑚
𝐴 = á𝑟𝑒𝑎 𝑑𝑒𝑙 𝑚𝑜𝑑𝑒𝑙𝑎𝑚𝑖𝑒𝑛𝑡𝑜 ( )
𝑠
𝑚
𝑉𝑒 = 𝑣𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑚𝑒𝑑𝑖𝑎 𝑑𝑒 𝑙𝑎 𝑒𝑟𝑜𝑠𝑖ó𝑛 𝑑𝑒 𝑓𝑜𝑛𝑑𝑜 ( )
𝑠
Coeficiente “β”.
El coeficiente (β) depende del periodo de retorno del evento hidrológico extremo, según la
siguiente expresión.
𝛽 = 0.8416 + 0.0334 ln 𝑇𝑟
Coeficiente “µ”.
0.387
𝜇 = 1.00 − ∗𝑈
𝐿
Siendo U la velocidad media de aproximación del flujo y L, La longitud entre apoyos (para casos
de interés).
En la siguiente tabla se muestran los resultados de socavación para las secciones de interés.
PROFUNDIDAD DE
SECCIÓN
SOCAVACIÓN
CRITICA
(m)
TIEMPO
DE
RETORNO TR=140 TR=500
(años)
S-01 1.43 3.01
S-02 1.13 2.38
S-03 1.54 3.24
S-04 1.02 2.15
S-05 1.67 3.52
S-06 1.48 3.12
ANEXO A
HIDROLOGIA
ANEXO A.1
HIDROLOGIA ESTADISTICA
ESTUDIO HIDROLOGICO
A. DATOS
Max
Año Ene Feb Mar Abr May Jun Jul Ago Sep Oct Nov Dic
Anual
1964 21.0 19.0 24.0 11.0 42.0 1.0 13.0 13.0 15.0 22.0 23.9 15.0 42.0
1965 13.0 18.0 18.0 18.0 17.0 5.0 11.0 7.0 13.0 10.0 13.0 17.0 18.0
1966 17.0 35.0 28.0 17.7 18.5 9.0 0.0 8.5 16.5 17.6 11.0 24.0 35.0
1967 20.0 20.5 21.0 20.0 13.5 4.0 0.5 15.0 17.5 24.0 10.5 12.5 24.0
1968 19.0 21.0 25.0 29.0 13.0 7.5 10.0 16.0 15.0 22.0 14.5 20.0 29.0
1969 16.0 35.8 10.0 25.8 9.4 10.0 4.4 8.4 15.0 29.0 32.2 23.8 35.8
1970 21.8 24.2 23.0 18.2 14.0 13.8 12.2 11.4 15.0 25.4 15.0 15.3 25.4
1971 20.2 48.2 21.0 21.0 15.5 7.0 5.5 28.4 9.0 8.0 14.2 33.0 48.2
1972 21.2 20.4 21.2 18.0 12.5 0.0 4.0 14.0 11.0 16.0 28.0 18.0 28.0
1973 15.5 25.0 16.0 22.0 17.0 2.9 14.3 7.0 15.0 19.0 19.3 15.5 25.0
1974 13.0 21.0 19.0 34.0 9.7 8.0 32.0 12.8 9.0 25.0 18.0 17.0 34.0
1975 31.5 18.8 20.5 14.3 11.0 10.0 5.0 12.0 9.0 14.0 28.0 20.0 31.5
1976 26.5 19.1 21.0 25.0 12.8 9.1 3.0 9.0 9.2 23.4 23.4 13.0 26.5
Nº Datos 13 13 13 13 13 13 13 13 13 13 13 13 156.0
Promedio 19.7 25.1 20.6 21.1 15.8 6.7 8.8 12.5 13.0 19.6 19.3 18.8 16.8
Desv.Est. 5.2 9.1 4.4 6.2 8.3 4.0 8.4 5.6 3.1 6.3 7.1 5.6 1.9
Máxima 31.5 48.2 28.0 34.0 42.0 13.8 32.0 28.4 17.5 29.0 32.2 33.0 48.2
Mínima 13.0 18.0 10.0 11.0 9.4 0.0 0.0 7.0 9.0 8.0 10.5 12.5 0.0
Fuente: ANA y SENAMHI.
50.0
40.0
30.0
20.0
10.0
0.0
1962 1964 1966 1968 1970 1972 1974 1976 1978
AÑOS
ASIMETRIA DE LA ESTACION
𝑛
− 𝐿 𝑔( )
𝑉 𝐶𝑠 = 10 10
Donde:
𝐴 = −0.33 + 0.08 𝐶𝑠 𝑠𝑖 𝐶𝑠 0.90
𝐴 = −0.52 + 0.30 𝐶𝑠 𝑠𝑖 𝐶𝑠 >0.90
B= 0.94 − 0.26 𝐶𝑠 𝑠𝑖 𝐶𝑠 1.50
B= 0.55 𝑠𝑖 𝐶𝑠 1.50
𝐶𝑠 =
0.01
A= -0.33
B= 0.94
V(Cs)= 0.04
Kn= 2.175
YH= 1.72
pp 52.7 mm
YL= 1.23
pp= 17.09 mm
Hydrognomon HIDROESTA
Excel
4 2
FUNCION
DMax DMax DMax
DISTRIB.
0.096 0.096 0.137
NORMAL
DISTRIB.
0.051 0.053 0.125
LOGNORMAL
DISTRIB.
0.056 0.065 0.108
GUMBELL
DISTRIB.
0.053 0.056 0.088
PEARSON III
DISTRIB.
NO SE
LOG 0.064 0.098
AJUSTA
PEARSON III
WEBULL-NORMAL
WEIBULL-GUMBEL
WEIBULL-PERARSON
WEIBULL-LOGPERSON
DISTRIB.
DISTRIB. PMAX
T DISTRIB. DISTRIB. DISTRIB. LOG PMAX
PEARSON (corregida)
(años) NORMAL LOGNORMAL GUMBELL PEARSON mm
III mm
III
5 37.7553 37.1753 36.7703 37.35 37.319 37.18 42.01
70.00 70.87
63.49
61.51
60.00
Precipitacion en mm
57.35
53.05
50.00
47.03
y = 6.2457ln(x) + 32.556
42.01
40.00 R² = 0.9984
30.00
20.00
10.00
0.00
1 10 100 1000
ANEXO A.2
PARAMETROS DE LA CUENCA
Cuenca: Río Torres
DESCRIPCION UND VALOR
De la superficie
Área km2 115.32
Perimetro de la cuenca km 62.07
Ancho de la cuenca km 5.68
Longitud de la cuenca km 20.30
Alturas cuenca menor
Altura máxima Z mnsm 5222.00
Altura mínima Z mnsm 3534.00
Altura media Z mnsm 3886.81
Centroide
X Centroide m
Y Centroide m
Pendiente Longitudinal del río
Longitud del río m 16650.00
Cota mayor m 4829
Cota menor m 3534
Pendiente m/m 0.08
Parametros Generados
Tiempo de concentración
KIRPICH hr 1.55
KIRPICH (Modificado) hr 1.55
Giandiott hr 2.36
TEMEZ hr 1.54
Promedio hr 1.75
Promedio min 104.86
TIEMPO DE CONCENTRACIÓN
FORMULA KIRPICH KIRPICH (Modificado) CALIFORNIA TEMEZ
hr 1.547 1.546 1.547 1.538
min 92.799 92.763 92.799 92.301
Tc Promedio 92.666
Tc Menor 92.301
Tc mayor 92.799
AREA
RANGO PROMEDIO NC*A NCp
(KM2)
33 44 38.5 0.72 27.72
54 64 59 0.65 38.16
ANEXO A.3
CALCULO DE CAUDALES
o también: Y = a0 + a1 * x1 + a2 * x2
𝟏𝟎𝑲 ∗ 𝑻𝒎
𝑰=
𝒕𝒏K 2.5
= 3
m 0.1
= 2
n 0.7
= 5
CURVA INTENSIDAD-DURACION-FRECUENCIA
ESTACIÓN HUALLANCA
140.0
5
120.0
10
Intensidad de precipitación (mm/h)
100.0 25
50
80.0 100
140
60.0 500
40.0
20.0
0.0
0.0 50.0 100.0 150.0 200.0 250.0 300.0 350.0 400.0
Duración (min)
Periodo de
TR = 140
retorno (años)
Curva Intensidad -
Duración - Frecuencia (IDF)
I: Intensidad
2.5347 0.1223 máxima (mm/h)
I T: Período de
10 T
= retorno en años
t: Duración de
la precipitación
0.75 (min)
t
12
10
8
6
4
2
0
0-20
5-10
10-15
15-20
20-25
25-30
30-35
35-40
40-45
45-50
50-55
55-60
60-65
65-70
70-75
75-80
80-85
Tiempo (min)
14
12
10
8
6
4
2
0
0-20
5-10
10-15
15-20
20-25
25-30
30-35
35-40
40-45
45-50
50-55
55-60
60-65
65-70
70-75
75-80
80-85
Tiempo (min)
B. MODELO ANALIZADO
HEC-HMS
TR = 140 AÑOS
SCS
ENTRADA
SALIDA
SNYDER
ENTRADA
SALIDA
CLARCK
ENTRADA
SALIDA
ANEXO B: HIIDRAULICA
SIMULACIÓN
SECCIONES CRITICAS
MAPA DE CALADO
MAPA DE VELOCIDADES