Estudio Hidrologico-Mayunmarca
Estudio Hidrologico-Mayunmarca
Estudio Hidrologico-Mayunmarca
VEREDAS DE LAS VÍAS PRINCIPALES Y SECUNDARIAS DEL CENTRO POBLADO RURAL DE SAN
SUB
MARTÍN DE MAYUNMARCA, DISTRITO DE ANDABAMBA, PROVINCIA DE ACOBAMBA, GERENCIA DE
DEPARTAMENTO DE HUANCAVELICA” ESTUDIOS
CONTENIDO
1. MEMORIA DESCRIPTIVA...........................................................................................2
2. UBICACIÓN Y DESCRIPCION DEL AREA DE ESTUDIO......................................2
3. ANÁLISIS DE CONSISTENCIA DE LA INFORMACIÓN PLUVIOMÉTRICA.....4
3.1. Análisis visual de hidrogramas...................................................................................4
3.2. Análisis de doble masa................................................................................................6
3.3 ANÁLISIS ESTADÍSTICO.........................................................................................8
3.4. ANÁLISIS DE TENDENCIAS.................................................................................12
4. REGIONALIZACION DE LAS PRECIPITACIONES MENSUALES......................13
5. ANALISIS ESTADISTICO DE DATOS HIDROLOGICOS.......................................18
5.1.1. Distribución Normal..............................................................................................19
5.1.2. Distribución de Frecuencia Gumbel.....................................................................19
5.1.3. Distribución de Frecuencia Tipo log Normal.......................................................19
5.1.4. Distribución de Frecuencia tipo log Pearson tipo III...........................................20
5.1.5. Distribución Pearson Tipo III (Gamma 3 Parametros).......................................20
6. ANALISIS HIDROLOGICO..........................................................................................25
6.1 INTENSIDADES DE LLUVIAS............................................................................25
6.2. CAUDALES MAXIMOS.........................................................................................26
6.3. DETERMINACIÓN DEL CAUDAL DE AGUAS PLUVIALES...........................27
6.4. CALCULO DE CAUDALES:..................................................................................29
7. CONCLUSIONES...........................................................................................................31
8.ANEXOS ..................................................................................................................32
DATOS HISTORICOS REGISTRADOS EN LA ESTACION ACOBAMBA...........32
ESTUDIO HIDROLOGICO
1. MEMORIA DESCRIPTIVA
a) Pluviometría
El área en la que está enmarcado el presente estudios, abarca esencialmente toda la zona
rural y urbana del Centro poblado de San Martin De Mayunmarca y las zonas cercanas,
distrito de Andabamba:
Provincia: Acobamba
Distrito: Andabamba
Latitud Sur: 12° 41' 28" / Longitud Oeste: 74° 37' 17" / Altitud: 3460 msnm
Ubigeo: 090202
2.1 ACCESOS
Mientras que las vías de acceso al distrito de Andabamba son con las ciudades de
Acobamba, Huancayo (capital de la Región Junín) y Distrito de Paucará; realizándose
principalmente a través de las siguientes rutas:
TIPO DE DISTANCIA
DE – A : TIPO DE VIA
CARRETERA Aprox. (Km.)
HUANCAVELICA - HUINTOPATA V. NACIONAL ASFALTADA 5.46 KM
HUINTOPATA – YAULI - CCACCAMARCA V. DEPARTAMENTAL AFIRMADA 17.21 KM
AFIRMADA EN MALAS
CONDICIONES
CCACCAMARCA – AMBATO - PUCAPAMPA V. VECINAL 6.55 KM
PUCAPAMPA – CHECCO CRUZ V. DEPARTAMENTAL AFIRMADA 30.27 KM
CHECCO CRUZ - PALMIRA V. NACIONAL ASFALTADA 8.65 KM
PALMIRA - MAYUNMARCA V. VECINAL AFIRMADA 11.39 KM
TOTAL 79.53 KM
Fuente: ELABORACIÓN PROPIA
1000
800
600
400
200
0
1960 1970 1980 1990 2000 2010 2020
Tiempo (años)
1200
1000
800
600
400
200
0
1960 1970 1980 1990 2000 2010 2020
Tiempo (años)
1500
1000
500
0
1960 1970 1980 1990 2000 2010 2020
Tiempo (años)
El análisis doble masa se realizó para el periodo común de información. La estación con más datos
faltantes es la estación “Choclococha”, y el periodo en donde se registró información con regularidad
corresponde a 1993 al 2006; por lo tanto en dicho periodo se realizara el análisis de doble masa.
Del cuadro N° 02, se elaboró los gráficos de las figuras N° 03 y 04 del análisis de doble masa, con el
promedio acumulado, para determinar la estación con mayor regularidad, y luego para identificar el
periodo con posible inconsistencia.
Fig. N° 09
ANALISIS DOBLE MASA ENTRE LAS ESTACIONES
CON EL PROMEDIO ACUMULADO
PROMEDIO ACUMULADO DE LAS TRES ESTACIONES (mm)
12,000
10,000
8,000
BA
M
BA
Y
CA
O
AC
LIR
6,000
A
O CH
OC
O CL
CH
4,000
2,000
0
2,000 4,000 6,000 8,000 10,000 12,000 14,000 16,000 18,000
PRECIPITACION ACUMULADA DE CADA ESTACION (mm)
Del gráfico doble masa presentado, se observa que las tres estaciones presentan regularidad durante
el periodo analizado.
40000
35000
ACUMULADO ESTACION LIRCAY (mm)
30000
25000
20000
15000
10000
5000
LIRCAY
0
0 5000 10000 15000 20000 25000 30000 35000
ACUMULADO ESTACION INDICE (ACOBAMBA) (mm)
Figura N° 04. Análisis doble masas entre acumulado estacione índice Acobamba vs Acumulado
estación Lircay.
EVALUACIÓN
Con los datos históricos de precipitación disponibles, el tratamiento a los mismos consiste en realizar:
Formulación de la Hipótesis
La hipótesis planteada HP, respecto a la homogeneidad de medias es: “Las medias muéstrales de los
datos correspondientes a los periodos dudoso y fiable son equivalentes estadísticamente un nivel de
confianza de 100 – α ”.
En consecuencia, la hipótesis alterna HA queda definida así: “Las medias muéstrales de los datos
correspondientes a los periodos dudoso y confiable no son equivalentes estadísticamente al nivel de
confianza de 100 - α ”.
Estadísticamente, las hipótesis anteriores se expresan como sigue:
HP : X́ 1= X́ 2
HA : X́ 1 ≠ X́ 2
Se lee, “las medias de los periodos dudoso y confiable de la Estación en un estudio son
estadísticamente iguales, o son estadísticamente diferentes”.
En ingeniería hidrológica, la exigencia acerca del nivel de significación es α =0.05 % .
( X́ 1− X́ 2 )
T c= (03)
Sd
Además:
1
1 1
Sd =S p
[ +
n1 n2 ] 2
(04)
1
( n1−1 ) S 21 + ( n2 −1 ) S 22
S p= [ n1 +n 2−2 ] 2
(05)
n1
21 2
S= 1 ∑ ( x 1− X́ 1 ) (06)
n1−1 i=1
n2
2 1 2
S=
2 ∑ ( x 1− X́ 2) (07)
n2 −1 i=1
Dónde:
Si:
|T c|≤ T t ( α =5 % ) Aceptar HP
|T c|>T t ( α=5 % ) Aceptar HA
En el primer caso, se dice que la prueba T es no significativa y; en el segundo, la prueba puede ser
significativa o altamente significativa, según que la discrepancia entre los estadísticos sea ligera o
alta, respectivamente. En este último caso, implica corregir la información, eliminando el salto del
periodo dudoso (eliminación del salto) y sin tener que recurrir a la prueba F.
i) Formulación de la hipótesis.
Similarmente a como se procedió en la prueba de T, las hipótesis de esta prueba se plantean
estadísticamente de la siguiente manera:
HP : S21 =S 22 ( α =0.05 )
HA :S 21 ≠ S22 ( α =0.05 )
Se lee, “las variancias de los periodos dudoso y confiable se la Estación en estudio son
estadísticamente iguales o son estadísticamente diferentes”.
2 2
Donde S1 y S2 son las variancias muéstrales correspondientes a los periodos dudoso y confiable de
la Estación en estudio.
S 22 ( x )
F C= 2 , si S22 ( x ) > S21 ( x ) (09)
S1( x )
G . L . D.=n2−1 (10.b)
α =0.05( N c =95 %)
Donde:
G . L . N = Grados de libertad del numerador.
G . L . D= Grados de libertad del denominador.
N c = Nivel de confianza.
Con los grados de libertad del numerador y el denominador y el nivel de significación α =0.05 , se
obtiene de las tablas el valor teórico correspondiente o directamente del Software con los mismos
datos de entrada.
Si:
F c ≤ F t (95 %) Aceptar HP (muestra consistente)
F c > F t ( 95 %) Aceptar HA (muestra inconsistente)
estadísticos sea leve o exista una gran diferencia, respectivamente. En este último caso, debe
eliminarse el salto corrigiendo los datos del periodo dudoso partiendo del confiable.
Una de las metodologías empleadas para la corrección de saltos es mediante la utilización de las
siguientes expresiones:
Donde:
x 1= Dato a corregir
S= Variancia
X́ = Media
Para comprobar si la información corregida está dentro de los límites de aceptación con el 95% de
probabilidades, se analiza estadísticamente tanto la media y la desviación Standard y, comparándolo
con el periodo confiable, repitiendo el procedimiento ya descrito, vale decir desarrollando las pruebas
T y F, respectivamente, cuyos resultados deben ser confiables.
Cabe mencionar que la información corregida no será una información ideal o naturalizada, puesto
que mantendrá un nivel de incertidumbre en relación a los valores individuales es decir la magnitud
de años secos o húmedos; pero para propósito de cálculo ó prognosis quedará libre de
inconsistencia. Podemos decir entonces que la información perdida no podrá recuperarse jamás, de
allí la enorme importancia de mantener una red bien manejada y controlada.
Normalmente se asume que el periodo reciente es el más confiable, en los países desarrollados de
Europa ó América. En nuestro país (Perú) sin embargo ocurre a la inversa, puesto que la situación
económica de los últimos años y la informalidad existente hacen afirmar que los periodos recientes en
general son los menos confiables.
Es importante también tener presente que frecuentemente los periodos de análisis de media y
varianza son coincidentes esto es n1 y n2.
La serie que debería considerarse será la de valores anuales, por cuanto el criterio de aplicación de
1 - 2 y 12/22, corresponden a series independiente semanales o mensuales aunque se
incremente con ello la longitud del registro “n”.
En el periodo comprendido entre los años 1989 al 1994, se observaron precipitaciones muy
inferiores en la estación “Acobamba” con respecto a su propio comportamiento en el resto de los
años y con respecto al comportamiento de la estación “Lircay” y “Choclococha”.
Durante el periodo 1997 al 2002, se observa que las precipitaciones en “Choclococha” se
incrementaron en una proporción mayor a lo que se produjo en las estaciones de “Acobamba” y
“Lircay”.
IDENTIFICACION
De acuerdo a la visita realizada a las estaciones no se identificó justificación física para que se el
registro histórico este afectado por tendencia alguna.
EVALUACION
Dónde:
t : Tiempo cronológico
R : Coeficiente de correlación
R( N −2)1/2
Tc=
(1−R2 )1/ 2
Al comparar el Tc con el Tt, se utiliza el siguiente criterio para decidir si la tendencia en la media es
significativa o no:
Los resultados de las pruebas estadísticas realizadas que determinan si es significativo o no alguna
tendencia son las siguientes:
Cuadro N° 03.
Se generaron las precipitaciones mensuales para cada cuenca, estas se obtuvieron mediante la
regresión lineal de las tres estaciones con las que se contaba y las altitudes de estas.
Para generar los datos de precipitación mensual para la cuenca deseada, se trabaja con la altitud a
la que se encuentra dicha cuenca este valor se remplaza en la ecuación obtenida de la regresión
lineal y se obtiene la precipitación generada.
7 8 0 0 7 8 9 3
CHOCLOCOCH 171.7 87.6 26.9 13.4 20.0 35.7 61.0 69.2
A 4583 5 169.78 157.47 7 3 8 9.29 2 4 8 6 119.59 942.05
Cuadro N° 04. Media de las precipitaciones mensuales de estaciones aledaña a la zona de estudio
(mm).
VARIACION DE LA PRECIPITACION
200.00
180.00
160.00
Precipitacion (mm)
140.00
120.00
ACOBAMBA
100.00 CHOCLOCOCHA
LIRCAY
80.00
60.00
40.00
20.00
0.00
1 2 3 4 5 6 7 8 9 10 11 12
MES
950.00
850.00 f(x) = 0.14 x + 275.22
750.00 R² = 0.82
650.00
550.00
450.00
350.00
3300 3400 3500 3600 3700 3800 3900 4000 4100 4200 4300 4400 4500 4600
ALTITUD (msnm)
Sabiendo:
Se multiplica al factor de relación obtenido para cada cuenca con la precipitación tomada como base,
en nuestro caso se tomó a la estación Acobamba, debido a que se encuentra a una altitud similar
que las cuencas a que se va a trabajar:
Pf =Po∗Fr .
DONDE:
El tiempo promedio, en años, en que el valor del caudal pico de una creciente determinada
es igualado o superado una vez cada “T” años, se le denomina Período de Retorno
“T”. Si se supone que los eventos anuales son independientes, es posible calcular la
probabilidad de falla para una vida útil de n años.
El criterio de riesgo es la fijación, a priori, del riesgo que se desea asumir por el caso de que
la obra llegase a fallar dentro de su tiempo de vida útil, lo cual implica que no ocurra un
evento de
magnitud superior a la utilizada en el diseño durante el primer año, durante el segundo, y así
sucesivamente para cada uno de los años de vida de la obra.
El riesgo de falla admisible en función del período de retorno y vida útil de la obra está dado
por:
R=1−(1−1/T )n
Si la obra tiene una vida útil de n años, la fórmula anterior permite calcular el
período de retorno T, fijando el riesgo de falla admisible R, el cual es la
probabilidad de ocurrencia del pico de la creciente estudiada, durante la vida útil de
la obra. (Ver Figura Nº 11).
Existen varias formas de determinar el tiempo de concentración, ya sea haciendo uso de las
características hidráulicas de la cuenca, estimando velocidades, o haciendo uso de fórmulas empíricas
propuestas por varios autores a fin de ahorrar tiempo en su cálculo, de las cuales la Fórmula de Kirpich,
es una de las más utilizadas.
t c =0.01947∗L0.77∗S−0.385
L:= 281.89 m
S = 8.61%
Fuente. La ecuación para tiempo de concentración se tomó del Manual de Hidrología, Hidráulica y
Drenaje del MTC.
Con el fin de ajustar la serie anual de precipitación máxima diaria de la estación Acobamba en
función de distribución probabilística teórica o también llamada el uso de números aleatorios.
Por lo que para el análisis de frecuencia se aplicó el análisis de funciones de la distribución
teórica:
Normal, Gumbel, Log-normal, Log-Pearson III y Pearson III, por ser las más usadas en
Hidrología para caso de eventos máximos.
i. Distribución Normal
Donde:
f(x) = densidad de probabilidad
x = variable aleatoria
5.1.2. Distribución de Frecuencia Gumbel
Esta distribución Gumbel 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:
x e x
f ( x ) e
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 ( y y ) 2
f ( x) .exp ; x 0
2 y
2
x. y 2
Donde:
y = log(x)
f(x) = densidad de probabilidad
x = variable aleatoria
µy = media de los logaritmos naturales de x
σy = desviación estándar de los logarítmicos naturales de x
e = base de los logaritmos naturales.
Donde:
a = parámetro de escala
B = parámetro de forma
y = límite inferior de la distribución
5.1.5. Distribución Pearson Tipo III (Gamma 3 Parametros)
La función de densidad es:
−( x− x0 )
γ −1 β
( x−x 0 ) e
f ( x )= γ
β Γ (γ )
Válido para:
x0 ≤ x < ∞
-∞ < x0 < ∞
0<β<∞
0<γ<∞
Donde:
x0: origen de la variable x, parámetro de posición
γ: parámetro de forma
β: parámetro de escala.
b) Prueba Kolmogorov – Smirnov
Con un valor crítico d que depende del número de datos y el nivel de significancia
seleccionado (Tabla Nº 03). Si D<d, se acepta la hipótesis nula. Esta prueba tiene la ventaja
sobre la prueba de X2 de que compara los datos con el modelo estadístico sin necesidad
de agruparlos. La función de distribución de probabilidad observada se calcula como:
Fo(xm) = 1- m / (n+1)
TAMAÑO DE
LA
α = 0.10 α = 0.05 α = 0.01
MUESTRA
5 0.51 0.56 0.67
10 0.37 0.41 0.49
15 0.30 0.34 0.40
20 0.26 0.29 0.35
25 0.24 0.26 0.32
30 0.22 0.24 0.29
35 0.20 0.22 0.27
40 0.19 0.21 0.25
Uno de los primeros pasos en muchos proyectos de diseño es la determinación del evento de
lluvia a usar.
Las curvas intensidad – duración – frecuencia son un elemento de diseño que relacionan la
intensidad de la lluvia, la duración de la misma y la frecuencia con la que se puede
presentar, es decir su probabilidad de ocurrencia o el periodo de retorno.
Para determinar estas curvas IDF se necesita contar con registros pluviográficos de lluvia en
el lugar de interés y seleccionar la lluvia más intensa de diferentes duraciones en cada año,
con el fin de realizar un estudio de frecuencia con cada una de las series así formadas. Es
decir, se deben examinar los hietogramas de cada una de las tormentas ocurridas en un
año y de estos hietogramas elegir la lluvia correspondiente a la hora más lluviosa,
a las dos horas más lluviosas, a las tres horas y así sucesivamente. Con los valores
seleccionados se forman series anuales para cada una de las duraciones elegidas. Estas
series anuales están formadas eligiendo, en cada año del registro, el mayor valor
observado correspondiente a cada duración, obteniéndose un valor para cada año y cada
duración.
Cabe indicar que formar las series anuales es un proceso largo y laborioso, que
involucra el examen cuidadoso de los rollos pluviográficos, la lectura de los valores, la
digitación de la información, la contrastación y verificación de los valores leídos con los
registros pluviométricos cercanos y el análisis de las tormentas registradas para encontrar
los máximos valores registrados para cada una de las duraciones seleccionadas.
CURVA INTENSIDAD-DURACION-FRECUENCIA
60.0 ESTACIÓN GENERADA - MAYUNMARCA
Intensidad de precipitacion (mm/h)
50.0 2
5
40.0 Las 1
curva 0
30.0 2 s de
0
intens
20.0 2
5 idad
10.0 5 –
0 durac
0.0 ión –
0.0 5.0 10.0 15.0 20.0 25.0 30.0 35.0 40.0 45.0 50.0 55.0
Duracion (min) frecu
encia
tambi
én pueden expresarse como ecuaciones con el fin de evitar la lectura de la intensidad de
lluvia de diseño en un una gráfica. Un modelo general es el siguiente.
a
I=
( D+b)m
DURACIÓN DE LA COEFICIENTE
1 0.25
2 0.31
3 0.38
4 0.44
5 0.50
6 0.56
8 0.64
10 0.73
12 0.79
14 0.83
16 0.87
18 0.90
20 0.93
22 0.97
24 1.00
48 1.32
Fuente: Manual para el Diseño de Carreteras
6. ANALISIS HIDROLOGICO
( K T m)
I=
tn
Dónde:
ESTACIÓN ACOBAMBA
DURACIÓN PERÍODO DE RETORNO (años)
Hr min 2 5 10 20 25 50 100
0.08 5.00 34.88 43.49 48.00 51.71 52.80 55.90 58.69
0.17 10.00 29.73 37.07 40.91 44.07 45.00 47.64 50.02
0.25 15.00 26.18 32.65 36.03 38.82 39.63 41.96 44.05
0.33 20.00 23.56 29.38 32.43 34.94 35.67 37.77 39.65
0.67 40.00 21.54 26.85 29.63 31.93 32.60 34.51 36.24
0.83 50.00 19.91 24.83 27.39 29.52 30.14 31.91 33.50
1.00 60.00 16.47 20.54 22.67 24.43 24.94 26.40 27.72
Dónde:
Fuente: Tabla 1.a de la Norma Técnica O.S. 060 Drenaje Pluvial Urbano
𝑄 = 𝐶𝐼𝐴/3.6
Dónde:
COEFICIENT
SUB INTENSIDA
AREA E CAUDAL Q CAUDAL
CUENCA DE AREA (m2) D MAX.
(km2) ESCORREN. (m3/seg) (Lt/seg)
DRENAJE (mm/h)
"C"
1 19212.00 0.019 0.42 35.00 0.078 78.45
2 14769.00 0.014 0.42 35.00 0.060 60.31
3 32319.00 0.032 0.42 35.00 0.132 131.97
4 33658.00 0.034 0.42 35.00 0.137 137.44
5 44350.00 0.044 0.42 35.00 0.181 181.10
6 29210.00 0.029 0.42 35.00 0.119 119.27
7 43313.00 0.043 0.42 35.00 0.177 176.86
8 14190.00 0.014 0.42 35.00 0.058 57.94
9 38652.00 0.049 0.42 35.00 0.158 157.83
10 27729.00 0.028 0.42 35.00 0.113 113.23
11 40972.00 0.040 0.42 35.00 0.167 167.30
12 14134.00 0.014 0.42 35.00 0.058 57.71
13 9158.00 0.009 0.42 35.00 0.037 37.40
14 25546.00 0.026 0.42 35.00 0.104 104.31
15 18621.00 0.019 0.42 35.00 0.076 76.04
16 12835.00 0.013 0.42 35.00 0.052 52.41
7. CONCLUSIONES
BIBLIOGRAFIA USADA
8. ANEXOS