UNIVERSIDAD NACIONAL DE HUANCAVELICA

(Creada por Ley Nro. 25265)

Facultad de Ciencias de Ingenierı́a

Escuela Profesional De Ingenierı́a Civil-Hvca

INFORME 01:

DISEÑO DE UNA PRESA DE concreto EN EL DISTRITO

DE CASTROVIRREYNA, PROVINCIA DE
CASTROVIRREYNA, REGIÓN HUANCAVELICA

CÁTEDRA:
PRESAS Y OBRAS DE EMBALSE
CATEDRÁTICO:
Ing. GASPAR PACO, Carlos
ESTUDIANTE:
QUISPE TAIPE, Jhon Alex (2018151039)
CICLO:
X

Huancavelica, Julio 2023

 Índice general

1 OBJETIVOS 2

2 INTRODUCCIÓN 3

3 GENERALIDADES DEL PROYECTO 4

3.1 Antecedentes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4

3.2 Ubicación . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4

3.3 Geologı́a y Fisiografı́a . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4

3.4 Accesibilidad . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5

3.5 Información Hidrométrica . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5

4 ESTUDIOS PRELIMINARES 6

4.1 Delimitación de Cuenca . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6

4.1.1 Ubicación de Cuenca . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6

4.1.2 Procedimiento de la Delimitación de Cuenca . . . . . . . . . . . . . . 6

4.1.3 Parámetros Geomorfológicos de la Cuenca . . . . . . . . . . . . . . . . 10

4.2 Cálculo del Caudal Máximo (Máxima Avenida) . . . . . . . . . . . . . . . . . 11

4.2.1 Periodo de Retorno . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11

4.2.2 Riesgo de Falla . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11

4.2.3 Distribuciones de Probabilidad . . . . . . . . . . . . . . . . . . . . . . 12

4.2.4 DETERMINACIÓN DE PRECIPITACIONES . . . . . . . . . . . . . 26

ÍNDICE GENERAL i
 4.2.5 ANÁLISIS ESTADÍSTICO . . . . . . . . . . . . . . . . . . . . . . . . 28

4.2.6 PRECIPITACIÓN MÁXIMA PROBABLE . . . . . . . . . . . . . . . 30

4.2.7 Coeficiente de escorrentı́a . . . . . . . . . . . . . . . . . . . . . . . . . 36

4.2.8 MÉTODO RACIONAL . . . . . . . . . . . . . . . . . . . . . . . . . . 37

4.2.9 MÉTODO RACIONAL MODIFICADO . . . . . . . . . . . . . . . . . 38

4.3 Cálculo de Caudales Mensuales . . . . . . . . . . . . . . . . . . . . . . . . . . 39

4.3.1 Método de Lutz Schölz . . . . . . . . . . . . . . . . . . . . . . . . . . . 39

4.3.2 Generación de los Caudales Mensuales . . . . . . . . . . . . . . . . . . 47

4.3.3 Restricciones Del Modelo . . . . . . . . . . . . . . . . . . . . . . . . . 48

4.3.4 Resultado del los Caudales Medios Mensuales . . . . . . . . . . . . . . 49

4.4 Balance Hı́drico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50

4.4.1 Oferta Hı́drica . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50

4.4.2 Demanda Hı́drica . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51

4.4.3 Balance Hı́drico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55

4.5 Alturas Tı́picas de la Presa . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55

4.5.1 NAMINO (Nivel de Aguas Mı́nimas de Operación) . . . . . . . . . . . 55

4.5.2 NAMO (Nivel de Aguas Maximas de Operación) . . . . . . . . . . . . 59

4.5.3 NAME (Nivel de Aguas Máximas Extraordinarias) . . . . . . . . . . . 60

4.5.4 Borde Libre de la Presa . . . . . . . . . . . . . . . . . . . . . . . . . . 60

4.5.5 Alltura Total de la Presa . . . . . . . . . . . . . . . . . . . . . . . . . 61

4.6 Volumen de Almacenamiento de la Presa . . . . . . . . . . . . . . . . . . . . . 61

4.6.1 Curva Altura - Volumen - Área . . . . . . . . . . . . . . . . . . . . . . 61

5 DISEÑO DE PRESA DE CONCRETO 63

6 CONCLUSIONES 64

7 BIBLIOGRAFÍA 65

ÍNDICE GENERAL ii
8 ANEXOS 66

ÍNDICE GENERAL iii

 DEDICATORIA

El presente trabajo va dedicado a mis padres quienes me brindan el apoyo para poder
lograr mis metas profesionales, a los docentes de la EPICH quienes son la mejor guı́a en mi
formación profesional y también va dedicado para la comunidad estudiantil de la Universidad
Nacional de Huancavelica, exclusivamente hacia la Escuela Profesional de Ingenierıa Civil
Huancavelica (EPICH), teniendo como finalidad la ampliación de conocimientos por parte de
los estudiantes y docentes que gusten leer este contenido. Este presente informe cuenta con
información garantizada para el Curso de Presas y Obras de Embalse.

ÍNDICE GENERAL 1
 Capı́tulo 1

OBJETIVOS

Objetivo General

Diseñar una presa de tierra en la localidad de Sinto y Esmeralda, distrito y provincia

de Castrovirreyna,en la región de Huancavelica.

Objetivos Especı́ficos

Realizar un resumen y familiarizarse con los parámetros de diseño y base teórica nece-
sarios para el diseño de una presa de tierra.

Determinar el volumen de sedimentos.

Determinar el volumen util por método pico secuente.

Determinar curva elevación- área-volumen

Determinar borde libre de la presa.

Determinar el nivel del agua final en la presa y determinar la altura final de la misma

Determinar la estabilidad de la presa (volcamiento y deslizamiento)

CAPÍTULO 1. OBJETIVOS 2
 Capı́tulo 2

INTRODUCCIÓN

El agua es un recurso natural potencialmente renovable, vulnerable, indispensable

para la vida, con un alto valor considerado vital, que forma parte de todos los seres vivos y
necesaria para todo tipo de actividades humanas, sin embargo, a pesar de su importancia el
agua es uno de los recursos más deficientemente utilizados, pues muchas veces es utilizado
como un recurso inagotable.

El presente trabajo consiste en plantear un diseño de Presa de concreto, la cual va

estar ubicada en el departamento de Huancavelica, provincia y distrito de Castrovirreyna
entre las localidades, Sinto y Esmeralda.

El objetivo es determinar las dimensiones de la estructura, en base a los caudales

medio mensuales que permitan efectuar un análisis del balance de agua con fines de riego de
347 has irrigables y proyectar una adecuada estructura hidráulica que permita evacuar las
aguas de máxima avenida de tal forma que no represente un problema a la Infraestructura y
permita drenar libremente, ante la ocurrencia de caudales máximos normales y eventualmente
extraordinarios para un determinado periodo de retorno y probabilidad de ocurrencia, acorde
a consideraciones usuales de diseño y criterios asumidos.

CAPÍTULO 2. INTRODUCCIÓN 3
 Capı́tulo 3

GENERALIDADES DEL PROYECTO

3.1. Antecedentes

3.2. Ubicación

A. Ubicación Geográfica
El área de influencia para la zona de estudio se ubica en el Datum WGS 84, Zona 18
Sur en las coordenadas geográficas y UTM siguientes:

Tabla 3.1
Ubicación geográfica del proyecto

COORDENADAS
GEOGRÁFICAS UTM
Latitud Sur: 12°21’ 07.6”S Coordenada Norte: 8634510.00 S
Longitud Oeste: 75°00’ 57.0”W Coordenada Este: 498278.00 E

B. Ubicación Polı́tica

Tabla 3.2
Ubicación polı́tica del proyecto

UBICACIÓN POLÍTICA
Departamento: Huancavelica
Provincia: Tayacaja
Distrito: Acostambo
Localidad: Acostambo
Ubigeo:

3.3. Geologı́a y Fisiografı́a

CAPÍTULO 3. GENERALIDADES DEL PROYECTO 4

3.4. Accesibilidad

La vı́a más rápida de acceso al proyecto es

3.5. Información Hidrométrica

A. Fuentes de Información

La información utilizada para la realización del análisis hidrológico del área de estudio,
ha sido obtenida de documentos correspondientes a las siguientes instituciones: Servicio
Nacional de Meteorologı́a e Hidrologı́a (SENAHMI).

B. Red Hidrometeorológica

La cuenca en estudio no tiene estación meteorólogicas, pero se puede considerar como

las más cercanas a las siguientes estaciones hidrometeorológicas:

Tabla 3.3
Estaciones Hidrometeorológicas cercanas al punto de interés

LATITUD / ALTITUD PER. DE

ESTACIÓN TIPO FUENTE
LONGITUD msnm REGISTRO
12°21’38.8”S Convenc./
Acostambo 3603 SENAMHI 1989 - 2014
75°3’26.1”W Meteorológ.

CAPÍTULO 3. GENERALIDADES DEL PROYECTO 5

 Capı́tulo 4

ESTUDIOS PRELIMINARES

4.1. Delimitación de Cuenca

4.1.1. Ubicación de Cuenca

La cuenca en estudio se ubica en el departamento de Huancavelica, provincia Tayacaja

y distrito de Acostambo.

Tabla 4.1
Ubicación de la cuenca del proyecto

COORDENADAS
GEOGRÁFICAS UTM
Latitud Sur: 12°21’ 07.6”S Coordenada Norte: 8634510.00 S
Longitud Oeste: 75°00’ 57.0”W Coordenada Este: 498278.00 E

4.1.2. Procedimiento de la Delimitación de Cuenca

4.1.2.1. Obtención del DEM

Se descargó el DEM del MINAM, para poder ubicar nuestra punto de aforo de interés,
y posteriormente delimitar la cuenca.

El mapa geográfico es el ASTER GDEM (Aster Global Digital Elevation Model),

que se encuentran en formato GeoTIFF con coordenadas geográficas lat/long, con resolución
espacial de 30 metros.

CAPÍTULO 4. ESTUDIOS PRELIMINARES 6

Figura 4.1
Mapa geográfico ASTER GDEM - DEM de todo el Perú

4.1.2.2. Importación del DEM y el punto de aforo (Qgis)

Una vez identificado el punto de interés procedemos a importar el archivo DEM del
lugar (ASTGTM S14W076 dem.Tif), luego reocortamos, y reproyectamos a coordenadas
UTM y también importamos el punto de aforo.

Figura 4.2
DEM y punto de aforo del proyecto en QGIS

4.1.2.3. Delimitación en (GRASS GIS)

Una vez obtenido el DEM del proyecto, recortado y reproyectado en coordenadas

UTM, procedemos a empezar la delimitación en GRASS GIS, con la ayuda del complemento
(r.basin), el cual nos permite realizar la delimitación de una cuenca con mayor optimiza-
ción y rápida creación de capas, además también nos brinda los parámetros geomorfológicos

CAPÍTULO 4. ESTUDIOS PRELIMINARES 7

principales de la cuenca.

Figura 4.3
Interfaz de GRASS GIS

Figura 4.4
Capas de salida de la delimitación de cuenca - GRASS GIS

CAPÍTULO 4. ESTUDIOS PRELIMINARES 8

Figura 4.5
Mapa de delimitación de la micro cuenca Acostambo

Figura 4.6
Mapa de orden de rı́os de la micro cuenca Acostambo

CAPÍTULO 4. ESTUDIOS PRELIMINARES 9

4.1.3. Parámetros Geomorfológicos de la Cuenca

La morfologı́a de una cuenca queda definida por los parámetros de forma, relieve y
drenaje, para lo cual se han establecido una serie de parámetros, que a través de ecuaciones
matemáticas, sirven de referencia para la clasificación y comparación de las cuencas.

La cuenca en estudio (Acostambo) posee una extensión de 6.749 km2 , por lo que es
una Subcuenca, debido a que su extensión se encuentra por debajo de los 50 km2

Tabla 4.2
Parámetros geomorfolóficos de la micro cuenca (Acostambo)

PARÁMETROS MORFOMÉTRICOS DE LA CUENCA

Centroide este de la cuenca 499794.740
Centroide norte de la cuenca 8634957.050
Área de la cuenca [km^2] 6.749
Perímetro de cuenca [km] 12.950
Elevación máxima [m s.l.m.] 4350
Elevación mínima [m slm] 3894
Desnivel [m] 456
Elevación media [m s.l.m.] 4183.289
Pendiente media 15.260
Longitud del vector de dirección [km] 1.586
Orientación predominante [grados desde el norte] 0.252
Coeficiente de compacidad 4.418
Relación de circularidad 0.506
Diámetro topológico 8.000
Relación de elongación 0.867
Factor de forma 1.996
Tiempo de concentración (Giandotti) - (hr) 0.905
Longitud del canal principal [km] 3.381
Pendiente media del canal principal [porcentaje] 9.455
Longitud media de ladera [m] 22.116
Magnitud 3.000
Orden de ríos máximo (strahler) 2
Número de flujos 5
Longitud total de la corriente [km] 6.072
Frecuencia de flujo de primer orden 0.445
Densidad de drenaje [km/km^2] 0.900
Relación de bifurcación (Horton) 4.000
Relación de longitud (Horton) 32064.000
Relación de área (Horton) 6.969
Relación de pendiente (Horton) 1.487

CAPÍTULO 4. ESTUDIOS PRELIMINARES 10

4.2. Cálculo del Caudal Máximo (Máxima Avenida)

En este capı́tulo se estimarán las máximas crecidas que se puedan presentar en las
quebradas en estudio, para lo cual se ha desarrollado un modelo precipitación-escorrentı́a.

Para obtener las precipitaciones máximas de 24 horas, calculamos la precipitación

máxima ocurrida en uno de los 365 dı́as durante cada año, desde el 1989 hasta el 2014.

4.2.1. Periodo de Retorno

Se define el periodo de retorno, Tr, de un evento de cierta magnitud como el tiempo

promedio que transcurre entre la ocurrencia de ese evento y la próxima ocurrencia de ese
evento con la misma magnitud.

Se define también como el tiempo que transcurre para que un evento sea excedido
o igualado, al menos una vez en promedio. Si P es la probabilidad de excedencia, se puede
demostrar matemáticamente que:

1
Tr = (4.1)
P

Es decir que si decimos que una descarga de 300 m3/s en cierto rı́o tiene un periodo
de retorno de 5 años, significa que en promedio, cada cinco años una descarga de tal magnitud
será observada en dicho rı́o.

4.2.2. Riesgo de Falla

En el diseño de obras hidráulicas expuestas a grandes avenidas, es necesario considerar

el riesgo asociado con el valor seleccionado para el diseño. Por lo común, el ingeniero diseña
una obra para resistir una avenida de cierta magnitud.

Se define el riesgo R de un diseño como la probabilidad de que la avenida para la cual

se diseña la obra sea excedida. Se entiende que ésta es una situación de riesgo, pues la obra
se diseña para soportar cierta avenida máxima, y crecientes mayores le podrı́an hacer daño o
incluso destruirla. El riesgo R puede entonces escribirse como:

CAPÍTULO 4. ESTUDIOS PRELIMINARES 11

  n
1
R=1− 1− (4.2)
Tr

Donde:

R: Riesgo

Tr: Periodo de retorno

n: Vida útil

La confiabilidad se define como el complemento del riesgo.

Conf iabilidad = 1 − R (4.3)

Se quiere que la obra tenga un riesgo pequeño de dañarse o, lo que es lo mismo, una
alta confiabilidad. No existen pautas definidas para determinar niveles aceptables de riesgo
de falla para la mayorı́a de tipos de estructuras.

Sin embargo, las probabilidades de falla no deberı́an exceder unos cuantos puntos
de porcentaje dependiendo de las consecuencias de falla para los habitantes, aguas abajo,
el uso de tierra, y la instalación misma. Podrı́an ser también significativas las consecuencias
ambientales de falla.

4.2.3. Distribuciones de Probabilidad

La distribución de probabilidad es una herramienta fundamental para la prospectiva,

puesto que con ella es posible diseñar un escenario de acontecimientos futuros considerando
las tendencias actuales de diversos fenómenos.

4.2.3.1. Distribución Normal:

Esta distribución es conocida como la campana de Gauss, por la similitud de la

distribución de probabilidad a una campana.

CAPÍTULO 4. ESTUDIOS PRELIMINARES 12

Figura 4.7
Resultados - Distribución Normal

Función de densidad de probabilidad f(x):

 2
1 x−µ
1 −
f (x) = √ e 2 σ (4.4)
2πσ

Donde:

f (x) = función de densidad normal de la variable x

f (x) = función de densidad normal de la variable x
µ = parámetro de localización.
σ = parámetro de escala.

Estandarizando con:

CAPÍTULO 4. ESTUDIOS PRELIMINARES 13

 x−µ
z= (4.5)
σ

z2
1 −
f (z) = √ e 2 (4.6)
2π

Función de distribución acumulada F(x):

 2
1 x−µ
Z x
1 −
F (x) = √ e 2 σ dx (4.7)
−∞ 2πσ

z
z2
1 −
Z
F (z) = √ e 2 dz (4.8)
−∞ 2π

Solución aproximada según Abramowitz y Stegun (1965):

z
z2
1 −
Z
F (z) = √ e 2 dz (4.9)
−∞ 2π

Parámetros:

La distribución normal es la única función que produce los mismos resultados de paráme-
tros estimados por método de momentos y máxima verosimilitud:

N
1 X
µ= xi (4.10)
N
i=1

1
N
" #
1 X 2
σ= (xi − µ)2 (4.11)
N −1
i=1

CAPÍTULO 4. ESTUDIOS PRELIMINARES 14

4.2.3.2. Distribución Log Normal de 2 Parámetros:

Figura 4.8
Resultados - Distribución Log Normal de 2 Parámetros

Función de densidad:

 2
1 ln(x) − µy
Z x
1 −
F (x) = √ e 2 σy dx (4.12)
−∞ 2πσy

Función de distribución en términos de y=lnx:

y − µy ln(x) − µy
z= =
σy σy

Función de la distribución acumulada:

CAPÍTULO 4. ESTUDIOS PRELIMINARES 15

  2
1 ln(x) − µy
Z x
1 −
F (x) = √ e 2 σy dx (4.13)
−∞ 2πσy

 2
1 y − µy
Z y
1 −
F (y) = √ e 2 σy dy (4.14)
−∞ 2πσy

y − µy ln(x) − µy
z= =
σy σy

Se obtiene, la distribución normal estándar:

Z Z
1 1 z2 2
F (Z) = √ e− 2 ( 2
)
dz (4.15)
S 2π −∝

4.2.3.3. Distribución Log Normal de 3 Parámetros:

Figura 4.9
Resultados - Distribución Log Normal de 3 Parámetros

CAPÍTULO 4. ESTUDIOS PRELIMINARES 16

 Función de densidad:

La función de probabilidad, de la distribución es:

 2
1 ln(x−x0 )−µy
1 −2 σy
f (x) = √ e (4.16)
(x − x0 )σy 2π

Donde:

x0 = Parámetros de posición en el dominio x.

µy = Parámetros de escala en el dominio x.
σ 2 y = Parámetros de forma en el dominio x.

Estimación de parámetros, método de máxima verosimilitud:

Los parámetros de la distribución son:

n
1X
µy = ln(xi − x0 ) (4.17)
n
i=1

" n
#1
2
1X
σy2 = (ln(xi − x0 ) − µy )2 (4.18)
n
i=1

Donde:

x0 = Parámetro de posición
µy = Parámetro de escala, igual al promedio de los ln(x − x0 )
σ 2 y = Parámetro de forma ln(x − x0 )

El parámetro de posición se calcula, mediante un proceso iterativo como:

n X ln(xi − x0 )
X 1
(σy2 − µy ) + n =0 (4.19)
xi − x0 xi − x0
i=1 i=1

Función de distribución acumulada:

Se calcula haciendo la transformación, a una distribución normal estándar:

CAPÍTULO 4. ESTUDIOS PRELIMINARES 17

  2
1 ln(x − x0 ) − µy
Z z
1 −
F (z) = √ e 2 σy dz (4.20)
−∞ (x − x0 ) 2πσy

Donde:

ln(x − x0 ) − µy
z=
σy

z
z2
1 −
Z
F (z) = √ e 2 dz (4.21)
−∞ 2π

4.2.3.4. Distribución Gamma de 2 Parámetros:

Figura 4.10
Resultados - Distribución Gamma de 2 Parámetros

Función de densidad:

La función de probabilidad, de la distribución es:

CAPÍTULO 4. ESTUDIOS PRELIMINARES 18

 x
−
xγ−1 e β
F (x) = dx (4.22)
β γ Γ(γ)

Función de distribución acumulada:

Los parámetros de la distribución es:

x
−
Z x
xγ−1 e β
F (x) = dx (4.23)
0 β γ Γ(γ)

Estimación de parámetros, método de máxima verosimilitud::

En (1960) Greenwood Durand, presentan las siguientes relaciones aproximadas:

Para: 0 ≤ y ≤ 0.5772

0.5000876 + 0.1648852y − 0.0544274y 2

γ= (4.24)
y

Y para: 0.5772 < y ≤ 17

8.898919 + 9.05995y + 0.9775373y 2

γ= (4.25)
y(17.79728 + 11.968477y + y 2 )

Donde:

¯
y = ln(x̄) − ln(x) (4.26)

CAPÍTULO 4. ESTUDIOS PRELIMINARES 19

4.2.3.5. Distribución Pearson Tipo III:

Figura 4.11
Resultados - Distribución Pearson Tipo III

Función de densidad:

Se dice que una variable aleatoria x, tiene una distribución gamma de 3 parámetros si
su función de densidad de probabilidad es:

(x − x0 )
−
(x − x0 )γ−1 e β
f (x) = (4.27)
β γ Γ(γ)

Para:
x0 ≤ x < ∞
−∞ ≤ x0 < ∞
0<β<∞
0<γ<∞

Función de la distribución acumulada:

CAPÍTULO 4. ESTUDIOS PRELIMINARES 20

  2
1 ln(x) − µy
Z x
1 −
F (x) = √ e 2 σy dxΓ(γ) (4.28)
−∞ 2πσy

Donde:

x = Variable aleatoria gamma de 3 parametros o Pearson tipo III

x0 = Origen de la variable x, parámetro de posición
β = Parámetro de escala
γ =parámetro de forma
Γ(γ) =función gamma completa

La variable reducida y gamma de 3 parámetros o Pearson tipo III, es:

x − x0
y=
β

por lo que la función acumulada Pearson tipo III es :

y
y γ−1 e−y
Z
G(y) = dy (4.29)
0 Γ(γ)

CAPÍTULO 4. ESTUDIOS PRELIMINARES 21

4.2.3.6. Distribución Log Pearson Tipo III:

Figura 4.12
Resultados - Distribución Log Pearson Tipo III

Función de Densidad:

(LN (x) − x0 )
−
(LN (x) − x0 )γ−1 e β
f (x) = (4.30)
β γ Γ(γ)

Para:
x0 ≤ x < ∞
−∞ ≤ x0 < ∞
0<β<∞
0<γ<∞

Función de distribución acumulada

La función de distribución acumulada,de la distribución log-Pearson tipo III es la inte-

gral de la ecuación, de la siguiente manera:

CAPÍTULO 4. ESTUDIOS PRELIMINARES 22

 (ln(x) − x0 )
−
Z x
(ln(x) − x0 )γ−1 e β
F (x) = dx (4.31)
x0 β γ Γ(γ)

Donde:

x = Variable aleatoria log-Pearson tipo III

x0 = Parámetro de posición
β = Parámetro de escala
γ =parámetro de forma

La variable reducida de la distribución log-Pearsón tipo III,es:

LN (x) − x0
y= (4.32)
β

Por lo que la función acumulada Pearson tipo III es :

y
y γ−1 e−y
Z
G(y) = dy (4.33)
0 Γ(γ)

CAPÍTULO 4. ESTUDIOS PRELIMINARES 23

4.2.3.7. Distribución Gumbel:

Figura 4.13
Resultados - Distribución Gumbel

Función acumulada:

La función acumulada de la distribución de Gumbel es:

x−µ
−

F (x) = e−e α (4.34)

Para: − ∝< x <∝ Donde:

0 < γ <∝: Parámetro de escala

− ∝< µ <∝: Parámetro de posición

Función de densidad

Derivando la función de distribución acumulada con respecto a x tenemos:

1 − x−µ −e x−µ
α
f (x) = e α (4.35)
α

CAPÍTULO 4. ESTUDIOS PRELIMINARES 24

 La variable reducida Gumbel:

x−µ
y= (4.36)
α

4.2.3.8. Distribución Log Gumbel:

Función de distribución acumulada:

La función de distribución acumulada, se define de la siguiente manera:

ln(x)−µ
−
F (x) = ee α
(4.37)

Para:
− ∝< x <∝

Donde: 0 < α <∝: Parámetro de escala

− ∝< µ <∝: Parámetro de posición

Variable reducida:

lnx − µ
y= (4.38)
α

Proceso de cálculo:

Para el cálculo de los parámetros de la serie de datos: x1, x2, ..xn, se convierte a sus
logaritmos, luego se calcula la media y desviación estándar, con las siguientes ecuaciones:

ˆ LA MEDIA:
P
lnx
Xlnx =
N

ˆ LA DESVIACION ESTÁNDAR

rP
(lnx − Xlnx)2
Slnx =
N −1

CAPÍTULO 4. ESTUDIOS PRELIMINARES 25

4.2.4. DETERMINACIÓN DE PRECIPITACIONES

Se procede a obtener precipitaciones diarias, mensuales, anuales extraidas de la página

de estaciones hidrometereológicas del ANA SENAMHI (1989-2014).

CAPÍTULO 4. ESTUDIOS PRELIMINARES 26

 AÑOS ENE FEB MAR ABR MAY JUN JUL AGO SET OCT NOV DIC Total PROMEDIO MAX
1989 125.43 159.31 215.3 97.7 20.3 20.8 25.5 41.7 67.1 90 45 76.8 959.44 82.078 215.3
1990 149.2 55.7 60.8 47.2 18.6 56.8 16.1 40.6 40.6 96.7 117.2 121.5 821 68.417 149.2
1991 46.2 81.1 111.9 31.1 20.2 29.9 12.4 23.85 79 27.7 65.5 50.1 555.1 48.246 111.9
1992 75.9 77.7 100.4 18.2 27.6 14.8 14.6 18 28.9 75.7 39.4 54.6 518.2 45.483 100.4
1993 181.1 136 96.1 50.2 15.7 3.2 24.4 32.6 73.1 62 162.5 115.3 952.2 79.350 181.1
1994 170.9 205.6 153 76.6 7.2 3.7 1 8.6 41.4 75.9 31.3 105.1 880.3 73.358 205.6
1995 80 151.7 158.6 38.8 6.8 12.5 11.8 12.4 41.4 31.5 92 90.3 727.8 60.650 158.6
1996 104.8 136.6 87.3 52.5 4.3 29.4 25.5 7.5 24.3 52.2 64.5 102.7 636.7 57.633 136.6
1997 137.9 115 81.3 49.7 1.9 2 4.4 28.3 76.2 44.3 77.6 96.9 715.5 59.625 137.9
1998 132.6 145.2 55.2 7.6 6.6 5 25.5 14 11.6 108.3 88.1 90.9 665.1 57.550 145.2
1999 107.7 162.4 100.2 69.2 19 25.4 10 2.2 62 48.6 61.5 73.7 741.9 61.825 162.4

CAPÍTULO 4. ESTUDIOS PRELIMINARES

2000 125.3 146.3 107.9 39.6 53.3 8.7 21.5 20.1 22.9 102.7 13.9 88.9 751.1 62.592 146.3
2001 189.5 143.2 109.6 23.6 18.7 16.9 49 16.3 34 80.4 91.1 81.2 853.5 71.125 189.5
2002 52.7 202.5 133 32.7 20.3 2.5 50 29.8 41.8 73.5 64.1 117.5 820.4 68.367 202.5
2003 69.5 119.2 134.4 63.5 10.9 29.4 25.5 45.5 51.1 17.8 90.1 87.3 689.3 62.017 134.4
2004 60.2 71.7 96.4 33.3 17.7 48 11.8 17 40.7 70.6 119.3 116.6 703.3 58.608 119.3
2005 68.5 100.8 85.7 41.2 19.4 29.4 6 4.5 24.5 93.2 56 87.7 587.5 51.408 100.8
2006 147 119.5 104.5 49.7 6.5 12.3 25.5 27.9 27.4 78.8 111.7 96 781.3 67.233 147
2007 85.5 68.7 132.8 38.2 25.5 29.4 20 7.5 34.3 67.5 85.9 87.4 653.3 56.892 132.8
2008 131.8 120.5 62.3 35.9 16.6 18.1 4 5 53 101.2 32.3 85.3 666 55.500 131.8
2009 105.3 144.6 96.3 40.3 42.5 29.4 19.8 28.9 26.9 66.2 119.5 161.6 851.9 73.442 161.6
2010 209.2 101 116.8 16.6 11.5 11.3 25.5 29.8 9.6 76.2 43.3 125.9 751.2 64.725 209.2
2011 159.5 215.3 163.7 62.1 38.6 2 17.5 6.3 57.3 75.2 69.6 126.9 994 82.833 215.3
2012 77.2 252.9 55.8 122.8 37.1 9.5 10.3 6.8 32.7 48.1 77.9 201.4 932.5 77.708 252.9
2013 143.2 150.1 120.6 52.6 50.2 6.8 8.3 29.5 30.5 42.8 53.4 131.4 819.4 68.283 150.1
2014 141.2 106.6 118.6 85 15 29.4 25.5 23.85 44.3 63.05 88.2 125.75 466.4 72.204 141.2
Minimo 46.2 55.7 55.2 7.6 1.9 2 1 2.2 9.6 17.8 13.9 50.1 466.4 45.5
promedio 118.4 134.200385 109.942308 49.0730769 20.4615385 18.7153846 18.9 20.3269231 41.4 68.1 75.4 103.798077 749.8 64.9
maximo 209.2 252.9 215.3 122.8 53.3 56.8 50 45.5 79 108.3 162.5 201.4 994 82.8

27
4.2.5. ANÁLISIS ESTADÍSTICO

En hidrologı́a, para representar las variables hidrológicas se requiere de modelos pro-

babilı́sticos. Para ello, es necesaria la selección del modelo más adecuado, lo que consiste en
elegir el tipo de modelo y estimar sus parámetros. Los modelos ası́ obtenidos nos permiten
la estimación de variables hidrológicas asociadas a diferentes probabilidades.

Se analiza las precipitaciones máximas mediante la prueba de Smirnov Kolgomorov,

obteniendose que la mejor ajuste de distribución se da en una distribución Gumbel.

CAPÍTULO 4. ESTUDIOS PRELIMINARES 28

 PRUEBA DE BONDAD DE AJUSTE KOLMOGOROV - SMIRNOV
|Fo(xm)-F(xm)| F(xm) |Fo(xm)-F(xm)| F(xm) |Fo(xm)-F(xm)| F(xm) |Fo(xm)-F(xm)| F(xm) |Fo(xm)-F(xm)| F(xm) |Fo(xm)-F(xm)|
m P=X (m3/s) P (X) =m/n+1 Z = X-Xm/S F(xm) Normal Y=LnX (x-X0) Ln(x)-X0
Normal LogNormal LogNormal Pearson III Pearson III LogPearson III LogPearson III Gumbel Gumbel LogGumbel III LogGumbel
1 100.40 0.0370 -1.513 0.0652 0.028 4.6092 0.0354 0.0017 63.636 0.0402 0.0032 3.038 0.031 0.0059 0.0201 0.01695 0.0033 0.0337
2 100.80 0.0741 -1.502 0.0665 0.008 4.6131 0.0367 0.0374 64.036 0.0417 0.0324 3.042 0.032 0.0417 0.0211 0.05293 0.0038 0.0703
3 111.90 0.1111 -1.217 0.1118 0.001 4.7176 0.0878 0.0233 75.136 0.0947 0.0164 3.146 0.084 0.0267 0.0690 0.04211 0.0411 0.0700
4 119.30 0.1481 -1.026 0.1524 0.004 4.7816 0.1384 0.0097 82.536 0.1454 0.0027 3.210 0.137 0.0109 0.1232 0.02500 0.1038 0.0443
5 131.80 0.1852 -0.705 0.2405 0.055 4.8813 0.2509 0.0657 95.036 0.2555 0.0703 3.310 0.255 0.0697 0.2500 0.06479 0.2648 0.0796
6 132.80 0.2222 -0.679 0.2486 0.026 4.8888 0.2610 0.0388 96.036 0.2653 0.0431 3.317 0.265 0.0432 0.2615 0.03926 0.2791 0.0569
7 134.40 0.2593 -0.638 0.2618 0.003 4.9008 0.2775 0.0183 97.636 0.2813 0.0221 3.329 0.283 0.0233 0.2802 0.02090 0.3021 0.0428
8 136.60 0.2963 -0.581 0.2805 0.016 4.9171 0.3007 0.0044 99.836 0.3038 0.0075 3.346 0.307 0.0102 0.3063 0.00997 0.3337 0.0374
9 137.90 0.3333 -0.548 0.2919 0.041 4.9265 0.3146 0.0188 101.136 0.3172 0.0161 3.355 0.321 0.0125 0.3219 0.01146 0.3523 0.0190
10 141.20 0.3704 -0.463 0.3217 0.049 4.9502 0.3504 0.0200 104.436 0.3519 0.0185 3.379 0.358 0.0127 0.3618 0.00856 0.3989 0.0285
11 145.20 0.4074 -0.360 0.3594 0.048 4.9781 0.3945 0.0129 108.436 0.3945 0.0129 3.407 0.403 0.0047 0.4103 0.00287 0.4532 0.0457
12 146.30 0.4444 -0.332 0.3701 0.074 4.9857 0.4066 0.0378 109.536 0.4063 0.0381 3.414 0.415 0.0294 0.4235 0.02092 0.4676 0.0231
13 147.00 0.4815 -0.314 0.3769 0.105 4.9904 0.4144 0.0671 110.236 0.4138 0.0677 3.419 0.423 0.0585 0.4319 0.04957 0.4766 0.0049
14 149.20 0.5185 -0.257 0.3986 0.120 5.0053 0.4387 0.0798 112.436 0.4374 0.0811 3.434 0.448 0.0710 0.4581 0.06045 0.5044 0.0141
15 150.10 0.5556 -0.234 0.4075 0.148 5.0113 0.4486 0.1070 113.336 0.4470 0.1085 3.440 0.458 0.0980 0.4687 0.08690 0.5155 0.0401
16 158.60 0.5926 -0.015 0.4940 0.099 5.0664 0.5401 0.0525 121.836 0.5362 0.0564 3.495 0.549 0.0435 0.5641 0.02848 0.6105 0.0179
17 161.60 0.6296 0.062 0.5247 0.105 5.0851 0.5710 0.0587 124.836 0.5665 0.0631 3.514 0.580 0.0500 0.5954 0.03424 0.6399 0.0103

CAPÍTULO 4. ESTUDIOS PRELIMINARES

18 162.40 0.6667 0.083 0.5329 0.134 5.0901 0.5790 0.0876 125.636 0.5745 0.0922 3.519 0.588 0.0790 0.6035 0.06318 0.6474 0.0193
19 181.10 0.7037 0.564 0.7136 0.010 5.1990 0.7435 0.0398 144.336 0.7382 0.0345 3.628 0.748 0.0441 0.7615 0.05780 0.7845 0.0808
20 189.50 0.7407 0.780 0.7823 0.042 5.2444 0.8005 0.0597 152.736 0.7960 0.0553 3.673 0.802 0.0617 0.8134 0.07269 0.8266 0.0859
21 202.50 0.7778 1.114 0.8675 0.090 5.3107 0.8687 0.0909 165.736 0.8661 0.0883 3.739 0.868 0.0899 0.8742 0.09641 0.8750 0.0973
22 205.60 0.8148 1.194 0.8838 0.069 5.3259 0.8817 0.0669 168.836 0.8796 0.0648 3.754 0.880 0.0653 0.8857 0.07088 0.8842 0.0694
23 209.20 0.8519 1.287 0.9009 0.049 5.3433 0.8954 0.0436 172.436 0.8938 0.0420 3.772 0.893 0.0414 0.8978 0.04597 0.8939 0.0421
24 215.30 0.8889 1.444 0.9256 0.037 5.3720 0.9156 0.0267 178.536 0.9148 0.0259 3.801 0.913 0.0238 0.9156 0.02675 0.9083 0.0194
25 215.30 0.9259 1.444 0.9256 0.000 5.3720 0.9156 0.0104 178.536 0.9148 0.0112 3.801 0.913 0.0132 0.9156 0.01029 0.9083 0.0176
26 252.90 0.9630 2.411 0.9921 0.029 5.5330 0.9797 0.0167 216.136 0.9812 0.0182 3.962 0.976 0.0132 0.9748 0.01188 0.9602 0.0028

X 159.1885 X LOG 5.0423 α 9.92 α 209.72 α 30.3015 α 0.1869

S 38.8631 S LOG 0.2397 β 12.3370 β 0.0165 µ 141.7000 µ 4.9344
CS 0.6349 CS LOG 0.1381 X0 36.7643 X0 1.5715
n 26 Dmáx 0.14801 0.1070 0.10853 0.09799 0.0964 0.09727
a 0.05 Dcritico > Dmax Si se ajusta Si se ajusta Si se ajusta Si se ajusta Si se ajusta Si se ajusta
Δ crítico 0.267 Mejor Ajuste 5 3 4 1 1 2

29
4.2.6. PRECIPITACIÓN MÁXIMA PROBABLE

Se determina la precipitación máxima probable para los diferentes peridodos de re-

torno, usando el factor de frecuencia para una distribución Gumbel:

P
x1
x=
n

sP
n
− xm )2
i=1 (xi
S=
n−1
√
6
α= S
π
  
T
Yt = −ln ln
T −1

µ = x − 0.45S

P = µ + Yt α

T[años] Yt P[mm] P[ corregido]

2 0.367 152.806 172.671
5 1.500 187.150 211.480
10 2.250 209.889 237.175
25 3.199 238.620 269.641
50 3.902 259.935 293.726
100 4.600 281.091 317.633
200 5.296 302.171 341.453
500 6.214 329.981 372.879
1000 6.907 351.000 396.630

Las intensidades máximas pueden calculadas mediante la metodologı́a de Dick Peschke

(Guevara, 1991) que relaciona la duración de la tormenta con la precipitación máxima en 24
horas. La expresión es la siguiente:

d 0.25
Pd = P24h ( ) (4.39)
1440

Donde:
Pd = precipitación total (mm)
d = duración en minutos
P24h = precipitación máxima en 24 horas (mm)

CAPÍTULO 4. ESTUDIOS PRELIMINARES 30

 Las curvas de intensidad-duración-frecuencia, se han calculado indirectamente, me-
diante la siguiente relación:
KT m
I= (4.40)
tn

Donde:
I = Intensidad máxima (mm/h)
K, m, n = factores caracterı́sticos de la zona de estudio
T = perı́odo de retorno en años t = duración de la precipitación equivalente al tiempo de
concentración (min)

CAPÍTULO 4. ESTUDIOS PRELIMINARES 31

 INTENSIDAD DE PRECIPITACIÓN

Periodo de Retorno (Años)

Duracion Duracion
2.00 5.00 10.00 25.00 50.00 100.00 200.00 500.00 1000.00
(Horas) (min)
Intensidad (mm/hr)
0.083 5 502.98 616.03 690.88 785.45 855.61 925.25 994.64 1086.18 1155.36
0.167 10 299.07 366.29 410.80 467.03 508.75 550.16 591.41 645.85 686.98
0.250 15 220.65 270.25 303.08 344.57 375.35 405.90 436.34 476.50 506.85
0.333 20 177.83 217.80 244.26 277.70 302.50 327.13 351.66 384.02 408.48
0.417 25 150.43 184.24 206.62 234.90 255.89 276.71 297.47 324.84 345.53
0.500 30 131.20 160.69 180.21 204.88 223.18 241.35 259.45 283.33 301.37
0.583 35 116.88 143.15 160.54 182.51 198.82 215.00 231.12 252.39 268.47
0.667 40 105.74 129.50 145.24 165.12 179.87 194.51 209.10 228.34 242.89
0.750 45 96.80 118.56 132.96 151.16 164.66 178.06 191.42 209.04 222.35
0.833 50 89.44 109.55 122.86 139.68 152.15 164.54 176.87 193.15 205.46
0.917 55 83.27 101.99 114.38 130.04 141.65 153.18 164.67 179.83 191.28
1.000 60 78.01 95.55 107.16 121.82 132.71 143.51 154.27 168.47 179.20
2.000 120 46.39 56.81 63.72 72.44 78.91 85.33 91.73 100.17 106.55
3.000 180 34.22 41.92 47.01 53.44 58.22 62.96 67.68 73.91 78.61
4.000 240 27.58 33.78 37.89 43.07 46.92 50.74 54.54 59.56 63.36
5.000 300 23.33 28.58 32.05 36.43 39.69 42.92 46.14 50.38 53.59
6.000 360 20.35 24.92 27.95 31.78 34.62 37.43 40.24 43.94 46.74
7.000 420 18.13 22.20 24.90 28.31 30.84 33.35 35.85 39.15 41.64
8.000 480 16.40 20.09 22.53 25.61 27.90 30.17 32.43 35.42 37.67
9.000 540 15.01 18.39 20.62 23.45 25.54 27.62 29.69 32.42 34.49
10.000 600 13.87 16.99 19.06 21.66 23.60 25.52 27.43 29.96 31.87
11.000 660 12.92 15.82 17.74 20.17 21.97 23.76 25.54 27.89 29.67
12.000 720 12.10 14.82 16.62 18.90 20.58 22.26 23.93 26.13 27.79
13.000 780 11.39 13.96 15.65 17.79 19.38 20.96 22.53 24.61 26.17
14.000 840 10.78 13.20 14.81 16.83 18.34 19.83 21.31 23.28 24.76
15.000 900 10.24 12.54 14.06 15.98 17.41 18.83 20.24 22.10 23.51
16.000 960 9.75 11.94 13.39 15.23 16.59 17.94 19.28 21.06 22.40
17.000 1020 9.32 11.41 12.80 14.55 15.85 17.14 18.43 20.12 21.40
18.000 1080 8.93 10.93 12.26 13.94 15.19 16.42 17.65 19.28 20.51
19.000 1140 8.57 10.50 11.77 13.39 14.58 15.77 16.95 18.51 19.69
20.000 1200 8.25 10.10 11.33 12.88 14.03 15.17 16.31 17.81 18.95
21.000 1260 7.95 9.74 10.92 12.42 13.53 14.63 15.73 17.17 18.27
22.000 1320 7.68 9.41 10.55 11.99 13.06 14.13 15.19 16.58 17.64
23.000 1380 7.43 9.10 10.20 11.60 12.64 13.66 14.69 16.04 17.06
24.000 1440 7.19 8.81 9.88 11.24 12.24 13.23 14.23 15.54 16.53

CURVA IDF
1400.0

1200.0
INTENSIDAD DE PRECIPITACION (MM/H)

1000.0

800.0

600.0

400.0

200.0

0.0
0.0 200.0 400.0 600.0 800.0 1000.0 1200.0 1400.0 1600.0
TR 2 años TR 5 años TR 10 años TR 25 años TR 50 años TR 100 años TR 200 años TR 500 años TR 1000 años
DURACION (MIN)

Para el cálculo de la intensidad se tiene la siguiente ecuación potencial y se determi-

CAPÍTULO 4. ESTUDIOS PRELIMINARES 32

nara mediante una regresión lineal multiple.

KT m
I= (4.41)
tn

Una vez determinada la ecuación de la intensidad se puede generar el hietograma de

diseño mediante el método del bloque alterno, para el perı́odo de retorno de 500 años.

CAPÍTULO 4. ESTUDIOS PRELIMINARES 33

CAPÍTULO 4. ESTUDIOS PRELIMINARES 34
 PRECIPITACION DE DISEÑO

K= 1669.662 INTENSIDAD
T= 500 352.617×𝑇 0.0968
I=
m= 0.1287
𝐷0.751
n= 0.750

MÉTODO DEL BLOQUE ALTERNO

PROFUNDIDA PROFUNDIDAD
DURACION INTENSIDAD TIEMPO PRECIPITACION
ACUMULADA INCREMENTAL

min mm/hr mm mm min mm

60 172.107 172.107 172.107 0-60 4.022
120 102.311 204.622 32.514 60-120 4.299
180 75.473 226.419 21.797 120-180 4.626
240 60.820 243.279 16.860 180-240 5.018
300 51.443 257.216 13.937 240-300 5.496
360 44.866 269.194 11.978 300-360 6.097
420 39.965 279.756 10.562 360-420 6.877
480 36.155 289.239 9.483 420-480 7.939
540 33.097 297.870 8.631 480-540 9.483
600 30.581 305.809 7.939 540-600 11.978
660 28.470 313.173 7.364 600-660 16.860
720 26.671 320.050 6.877 660-720 32.514
780 25.116 326.510 6.460 720-780 172.107
840 23.758 332.607 6.097 780-840 21.797
900 22.559 338.385 5.778 840-900 13.937
960 21.493 343.881 5.496 900-960 10.562
1020 20.537 349.126 5.244 960-1020 8.631
1080 19.675 354.143 5.018 1020-1080 7.364
1140 18.892 358.956 4.813 1080-1140 6.460
1200 18.179 363.582 4.626 1140-1200 5.778
1260 17.526 368.038 4.456 1200-1260 5.244
1320 16.924 372.337 4.299 1260-1320 4.813
1380 16.369 376.492 4.155 1320-1380 4.456
1440 15.855 380.513 4.022 1380-1440 4.155

CURVA INTENSIDAD-DURACION-FRECUENCIA
I: Intensidad máxima (mm/h)
T: Período de retorno en años = 500.00 años
t: Duración de la precipitación (min)

HIETOGRAMA DE PRECIPITACIÓN DE DISEÑO

200.0
180.0
PRECIPITACIÓN (MM)

160.0
140.0
120.0
100.0
80.0
60.0
40.0
20.0
-
0-60

60-120

120-180

180-240

240-300

300-360

360-420

420-480

480-540

540-600

600-660

660-720

720-780

780-840

840-900

900-960

960-1020

1020-1080

1080-1140

1140-1200

1200-1260

1260-1320

1320-1380

1380-1440

TIEMPO (MIN)

CAPÍTULO 4. ESTUDIOS PRELIMINARES 35

4.2.7. Coeficiente de escorrentı́a

Usamos la siguiente formula para hallar el coefieciente de escorrentia.

(Pd − Po )(Pd + 23Po )

C= (4.42)
(Pd + 11Po )2

Donde:

Pd = Precitación máxima diaria (mm) = 25.468

5000


Po = Umbral de escorrentı́a = CN − 50 = 16.667

CN = Número de curva. Sacamos los datos de las siguientes tablas:

Figura 4.14
Grupo hidrológico según descripción del suelo.

CAPÍTULO 4. ESTUDIOS PRELIMINARES 36

Figura 4.15
Números de curva para distintas condiciones hidrológicas

Consideramos un CN = 79 obtenemos:

C= 0.5

4.2.8. MÉTODO RACIONAL

Teniendo de la ecuación de la intensidad se obtiene para un periodo de retorno de 500

años una intensidad de 185.493[mm/hr]. Reemplazamos los datos en la formula y obtenemos:

Q = 0.2778 ∗ C ∗ I ∗ A

Q = 0.2778 ∗ 185.493 ∗ 6.749

Q= 174[m3/s]

CAPÍTULO 4. ESTUDIOS PRELIMINARES 37

4.2.9. MÉTODO RACIONAL MODIFICADO

El caudal según el método racional modificado se calcula con la siguiente formula:

Q = 0.278CIAK (4.43)

Teniendo los siguientes datos:

Área: 6.749[km2]
Longitud de cauce mayor: 7.4[km]
Pendiente promedio del cauce: 0.095[m/m]
N° de curva: 79, extraido de Figura 4.14 Y Figura 4.15
Precipitacion maxima 24H(mm):252.9
Tiempo de concentración(horas): calculado por el método de Giandottise tiene 0.9 horas
Coeficiente de uniformidad: Resulta 1.059 con:

Tc1.25
K =1+ (4.44)
Tc1.25 + 14

Intensidad máxima(mm/h): Resulta 185.493 con:

KT m
I= (4.45)
tn

Umbral de escorrentia: Resulta 13.291 con:

5000
Po = − 50 (4.46)
CN

Coeficiente de escorrentia: Resulta 0.5 con:

(Pd − Po )(Pd + 23Po )

C= (4.47)
(Pd + 11Po )2

Con todos estos valores ya calculados se reemplaza en la ecuación de caudal y se obtiene

como caudal máximo de diseño:

Q= 184.3[m3/s]

CAPÍTULO 4. ESTUDIOS PRELIMINARES 38

4.3. Cálculo de Caudales Mensuales

A partir de este registro de precipitaciones medias mensuales areales en la Subcuenca

Presa Sinto-Esmeralda se procede a aplicar el Modelo de Lutz Scholz para la Generación de
Caudales Mensualizados.

4.3.1. Método de Lutz Schölz

El modelo hidrológico de Lutz Schölz, ha sido estudiado y calibrado en 19 cuencas

de la sierra peruana, ubicadas en las regiones de Cuzco, Huancavelica, Junı́n y Cajamarca,
y es aplicable generalmente para pequeñas y medianas unidades hidrográficas con escasa
información hidrométrica. Para el presente estudio ha sido seleccionado por que la zona de
estudio tiene las caracterı́sticas y condiciones para las cuales ha sido establecido el modelo.
El modelo de precipitación – escorrentı́a de Lutz Schölz, fue propuesta por la Misión Técnica
Alemana en 1980 para el Ex-Programa Nacional de Pequeñas y Medianas Irrigaciones – Plan
Meris II.

El modelo hidrológico se basa fundamentalmente en el balance hı́drico y en paráme-

tros parciales de tipo determinı́stico. Este método permite combinar los factores que producen
e influyen en los caudales, tales como, la precipitación, evaporación, el almacenamiento y la
función de agotamiento natural de la cuenca, para el cálculo de las descargas en forma de
un modelo matemático. El método es completo en su determinación final, vale decir analiza
los acuı́feros, los pantanos, los nevados, las lagunas, rı́os, manantes, es decir, todo el recurso
hı́drico que interviene en la cuenca, además de manera sustancial de la precipitación efectiva
que cae en la cuenca. Este modelo hidrológico, es combinado por que cuenta con una estruc-
tura determinı́stica para el cálculo de los caudales mensuales para el año promedio (Balance
Hı́drico – Modelo determinı́stico); y una estructura estocástica para la generación de series
extendidas de caudal (Proceso Markoviano - Modelo Estocástico).

Determinado el hecho de la ausencia de registros de caudal en la sierra peruana, el

modelo se desarrolló tomando en consideración parámetros fı́sicos y meteorológicos de las
cuencas, que puedan ser obtenidos a través de mediciones cartográficas y de campo. Los
parámetros más importantes del modelo son los coeficientes para la determinación de la
precipitación efectiva, déficit de escurrimiento, retención y agotamiento de las cuencas. Los
procedimientos que se han seguido en la implementación del modelo son:

Cálculo de los parámetros necesarios para la descripción de los fenómenos de escorrentı́a

CAPÍTULO 4. ESTUDIOS PRELIMINARES 39

 promedio.

Establecimiento de un conjunto de modelos parciales de los parámetros para el cálculo

de caudales en cuencas sin información hidrométrica. En base a lo anterior se realiza el
cálculo de los caudales necesarios.

Calibración del modelo y generación de caudales extendidos por un proceso Markoviano

combinado de precipitación efectiva del mes con el caudal del mes anterior.

Este modelo fue implementado con fines de pronosticar caudales a escala mensual, te-
niendo una utilización inicial en estudios de proyectos de riego y posteriormente exten-
diéndose el uso del mismo a estudios hidrológicos con prácticamente cualquier finalidad
(abastecimiento de agua, hidroelectricidad etc.). Los resultados de la aplicación del mo-
delo a las cuencas de la sierra peruana, han producido una correspondencia satisfactoria
respecto a los valores medidos.

4.3.1.1. Ecuación del Balance Hı́drico

La ecuación fundamental que describe el balance hı́drico mensual en mm/mes es la

siguiente:

CMi = Pi − Di + Gi − Ai (4.48)

Donde:

CMi :Caudal mensual (mm/mes).

Pi :Precipitación mensual sobre la cuenca (mm/mes).

Di :Déficit de escurrimiento (mm/mes).

Gi :Gasto de la retención de la cuenca (mm/mes).

Ai :Abastecimiento de la retención (mm/mes).

Asumiendo:

Que para perı́odos largos (en este caso 1 año) el Gasto y Abastecimiento de la retención
tienen el mismo valor es decir Gi = Ai.

Q = C0 ∗ C ∗ P ∗ A ∗ R (4.49)

CAPÍTULO 4. ESTUDIOS PRELIMINARES 40

Que es la expresión básica del método racional.

Dónde:

Q : Caudal (m3/s).

c’ : Coeficiente de conversión del tiempo (mes/s).

C : Coeficiente de escurrimiento.

P : Precipitación total mensual (mm/mes).

AR : Área de la cuenca (m2).

4.3.1.2. Coeficiente de Escurrimiento

Dado que no se ha podido obtener una ecuación general del coeficiente de escorrentı́a
para toda la sierra, se ha desarrollado la fórmula siguiente, que es válida para la región sur:

C = 3.16x1012 xP −0.571 xET P −3.686 (4.50)

D = −1380 + 0.87P + 1.032ET P (4.51)

Dónde:

C : Coeficiente de escurrimiento.

D : Déficit de escurrimiento (mm/año).

P : Precipitación total anual (mm/año).

ETP : Evapotranspiración anual según Hargreaves (mm/año).

4.3.1.3. Precipitación Efectiva

Para el cálculo de la precipitación efectiva, se supone que los caudales promedio

observados en la cuenca pertenecen a un estado de equilibrio entre gasto y abastecimiento de
la retención.

La precipitación efectiva se calculó para el coeficiente de escurrimiento promedio, de

tal forma que la relación entre precipitación efectiva y precipitación total resulte igual al
coeficiente de escorrentı́a.

CAPÍTULO 4. ESTUDIOS PRELIMINARES 41

 Para fines hidrológicos se toma como precipitación efectiva la parte de la precipitación
total mensual, que corresponde al déficit según el método del USBR (precipitación efectiva
hidrológica es la antı́tesis de la precipitación efectiva para los cultivos). A fin de facilitar el
cálculo de la precipitación efectiva se ha determinado el polinomio de quinto grado:

P E = a0 + a1 P + a2 P 2 + a3 P 3 + a4 P 4 + a5 P 5 (4.52)

Dónde:

PE : Precipitación efectiva (mm/mes).

P : Precipitación total mensual (mm/mes).

ai : Coeficiente del polinomio.

Las siguientes muestran los lı́mites de la precipitación efectiva y los tres juegos de
coeficientes, ai, que permiten alcanzar por interpolación valores de C, comprendidos entre
0.15 y 0.50.

Tabla 4.3
Lı́mite superior para la precipitación efectiva

Rango
N° Curva Ecuación
mm/mes

Curva I PE = P-120.6 P > 177.8

Curva II PE = P - 86.4 P > 152.4

Curva III PE = P - 59.7 P > 127.0

CAPÍTULO 4. ESTUDIOS PRELIMINARES 42

Tabla 4.4
Coeficientes para el cálculo de la precipitación efectiva

Coeficiente Curva I Curva II Curva III

ao -0.047 -0.1065 -0.4177
a1 0.0094 0.1477 0.3795
a2 -0.0005 -0.0029 -0.0101
a3 0.00002 0.00005 0.0002
a4 -0.00000005 -0.0000002 -0.0000009
a5 2.00E-10 2.00E-10 1E-09
C 0.15 0.3 0.5

El rango de aplicación de los coeficientes de la ecuación

Polinómica de la PE está comprendida para 0 < P < 250 mm

Fuente: Generación de Caudales Mensuales en la Sierra Peruana – Lutz Schölz Programa Nacional
de Pequeñas y Medianas Irrigaciones PLAN MERIS II.

De esta forma es posible llegar a la relación entre la precipitación efectiva y precipi-

tación total:

12
Q X P Ei
C= = (4.53)
P P
i=1

Dónde:

C : Coeficiente de escurrimiento.

Q : Caudal anual (m3 /s).

P : Precipitación total anual (mm).

4.3.1.4. Retención de la cuenca

Bajo la suposición de que exista un equilibrio entre el gasto y el abastecimiento de la

reserva de la cuenca y además que el caudal total sea igual a la precipitación efectiva anual,
la contribución de la reserva hı́drica al caudal se puede calcular según las fórmulas:

Ri = CMi − Pi (4.54)

CAPÍTULO 4. ESTUDIOS PRELIMINARES 43

 CMi = P Ei + Gi − Ai (4.55)

Dónde:

CMi : Caudal mensual (mm/mes).

P Ei : Precipitación Efectiva Mensual (mm/mes).

Ri : Retención de la cuenca (mm/mes).

Gi : Gasto de la retención (mm/mes).

Ai : Abastecimiento de la retención (mm/mes).

Ri : Gi para valores mayores que cero (mm/mes).

Ri : Ai para valores menores que cero (mm/mes).

Sumando los valores de G o A respectivamente, se halla la retención total de la cuenca

para el año promedio, que para el caso de las cuencas de la sierra varı́a de 43 a 188 (mm/año).

4.3.1.5. Relación entre Descargas y Retención

Durante la estación seca, el gasto de la retención alimenta los rı́os, constituyendo el

caudal o descarga básica. La reserva o retención de la cuenca se agota al final de la estación
seca; durante esta estación la descarga se puede calcular en base a la ecuación:

Qt = Q0 e−a(t) (4.56)

Dónde:

Qt : Descarga en el tiempo t (m3 /s).

Q0 : Descarga inicial (m3 /s).

a : Coeficiente de agotamiento.

t : tiempo (s).

Al principio de la estación lluviosa, el proceso de agotamiento de la reserva termina,

comenzando a su vez el abastecimiento de los almacenes hı́dricos. Este proceso está descrito

CAPÍTULO 4. ESTUDIOS PRELIMINARES 44

por un déficit entre la precipitación efectiva y el caudal real. En base a los hidrogramas
se ha determinado que el abastecimiento es más fuerte al principio de la estación lluviosa
continuando de forma progresiva pero menos pronunciada, hasta el final de dicha estación.

4.3.1.6. Ceficiente de Agotamiento

Mediante la ecuación anterior se puede calcular el coeficiente de agotamiento .a”, en

base a datos hidrométricos. Este coeficiente no es constante durante toda la estación seca, ya
que va disminuyendo gradualmente.

Con fines prácticos se puede despreciar la variación del coeficiente .a”durante la esta-
ción seca empleando un valor promedio. El coeficiente de agotamiento de la cuenca tiene una
dependencia logarı́tmica del área de la cuenca.

Los análisis de las observaciones disponibles muestran, además cierta influencia del
clima, la geologı́a y la cobertura vegetal. Se ha desarrollado una ecuación empı́rica para la
sierra peruana.

a = 3.1249x1067 xA−0.1144 xET P −19.336 xT −3.369 xR−1.429 (4.57)

En principio, es posible determinar el coeficiente de agotamiento real mediante aforos suce-

sivos en el rı́o durante la estación seca; sin embargo cuando no sea posible ello, se puede
recurrir a las ecuaciones desarrolladas para la determinación del coeficiente .a”para cuatro
clases de cuencas:

Cuencas con agotamiento muy rápido, debido a temperaturas elevadas (mayores a 10◦
C) y retención que va de reducida (50 mm/año) a mediana (80 mm/año):

a = −0.000252lnA + 0.034 (4.58)

Cuencas con agotamiento rápido, cuya retención varı́a entre 50 y 80 mm/año y vegeta-
ción poco desarrollada (puna):

a = −0.00252lnA + 0.030 (4.59)

Cuencas con agotamiento mediano, cuya retención es alrededor de 80 mm/año y vege-

CAPÍTULO 4. ESTUDIOS PRELIMINARES 45

 tación mezclada (pastos, bosques y terrenos cultivados):

a = −0.00252lnA + 0.026 (4.60)

Cuencas con agotamiento reducido por alta retención, cuya retención es arriba de 100
mm/año y vegetación mezclada:

a = −0.00252lnA + 0.023 (4.61)

Dónde:

ˆ a : Coeficiente de agotamiento por dı́a.

ˆ A : Área de la cuenca (km2 ).

ˆ ETP : Evapotranspiración potencial anual (mm/año).

ˆ T : Duración de la temporada seca (dı́as).

ˆ R : Retención total de la cuenca (mm/año).

Hay tres tipos de almacenes hı́dricos naturales que inciden en la retención de la cuenca.

Tabla 4.5
Lı́mite superior para la precipitación efectiva

Lámina acumulada
(mm/año)
TIPO Pendiente de la cuenca

2% 8% 15%

Napa Freática 300 300 300

Lagunas - Pantanos 500

Nevados 500

Fuente: Generación de Caudales Mensuales en la Sierra Peruana – Lutz Schölz Programa Nacional
de Pequeñas y Medianas Irrigaciones PLAN MERIS II.

Las respectivas extensiones o áreas fueron determinadas de los mapas y aerofoto-

grafı́as. Los almacenamientos de corto plazo no son considerados para este caso, estando los
mismos incluidos en las ecuaciones de la precipitación efectiva.

La lámina de agua Ai que entra en la reserva de la cuenca se muestra en forma de

déficit mensual de la Precipitación Efectiva PEi. Se calcula mediante la ecuación:

CAPÍTULO 4. ESTUDIOS PRELIMINARES 46

  
R
Ai = ai (4.62)
100

Dónde:

Ai : Abastecimiento mensual déficit de la precipitación efectiva (mm/mes).

ai : Coeficiente de abastecimiento ( %).

R : Retención total de la cuenca (mm/año).

4.3.1.7. Abastecimiento de la Retención

El abastecimiento durante la estación lluviosa es uniforme para cuencas ubicadas en

la misma región climática. En las regiones del Cuzco y Apurı́mac, el abastecimiento comienza
en el mes de noviembre con 5 %, alcanzando hasta enero el valor del 80 % del volumen final.
Las precipitaciones altas del mes de febrero completan el 20 % restante, y las precipitaciones
efectivas del mes de marzo escurren directamente sin contribuir a la retención.
Tabla 4.6
Almacenamiento hı́drico durante la época de lluvias

REGIÓN OCT NOV DIC ENE FEB MAR ABR

Cuzco 0 5 35 40 20 0 100

Huancavelica 10 0 35 30 20 5 100

Junín 10 0 35 30 30 5 100

Cajamarca 25 -5 0 20 25 35 100

Fuente: Generación de Caudales Mensuales en la Sierra Peruana – Lutz Schölz Programa Nacional
de Pequeñas y Medianas Irrigaciones PLAN MERIS II.

4.3.2. Generación de los Caudales Mensuales

A fin de generar una serie sintética de caudales para perı́odos extendidos, se ha imple-
mentado un modelo estocástico que consiste en una combinación de un proceso Markoviano
de primer orden, según las ecuaciones anteriores con una variable de impulso, que en este
caso es la precipitación efectiva. Con la finalidad de aumentar el rango de valores generados

CAPÍTULO 4. ESTUDIOS PRELIMINARES 47

y obtener una óptima aproximación a la realidad, se utiliza además una variable aleatoria.
La ecuación integral para la generación de caudales mensuales es:

p
Qt = B1 + B2 (Qt−1 ) + B3 (P Et ) + z(S) 1 + r2 (4.63)

Dónde:

Qt : Caudal del mes t (m3/s).

Qt−1 : Caudal del mes anterior (m3/s).

P Et : Retención total de la cuenca (mm/año).

B1 : Factor constante o caudal básico (m3/s).

Se calcula los parámetros B1, B2, B3, r y S sobre la base de los resultados del modelo para
el año promedio por un cálculo de regresión con Qt como valor dependiente y Qt-1 y PEt,
como valores independientes. El proceso de generación requiere de un valor inicial, el cual
puede ser obtenido en una de las siguientes formas:

Empezar el cálculo en el mes para el cual se dispone de un aforo.

Tomar como valor inicial el caudal promedio de cualquier mes.

Empezar con un caudal cero, calcular un año y tomar el último valor como valor Qo
sin considerar estos valores en el cálculo de los parámetros estadı́sticos del perı́odo
generado.

4.3.3. Restricciones Del Modelo

El modelo presenta ciertas restricciones de uso o aplicación tales como:

El uso de los modelos parciales, únicamente dentro del rango de calibración establecido.

Su uso es únicamente para el cálculo de caudales mensuales promedio.

Los registros generados en el perı́odo de secas presentan una mayor confiabilidad que
los valores generados para la época lluviosa.

La aplicación del modelo se restringe a las cuencas en las que se ha calibrado sus
parámetros (sierra peruana: Cusco, Huancavelica, Junı́n, Cajamarca)

CAPÍTULO 4. ESTUDIOS PRELIMINARES 48

Es importante tener en cuenta las mencionadas restricciones a fin de garantizar una buena
optimización del modelo.

4.3.3.1. Calibración

Es recomendable relizar trabajos de campo, donde se realicen mediciones de caudales

en la subcuenca estudiada, con la finalidad de calibrar los resultados del modelo. En la
calibración se modifican los valores del coeficiente de escorrentı́a, el valor del almacenamiento
y el coeficiente de agotamiento, de tal manera de encontrar la mejor combinación de estos
parámetros.

4.3.4. Resultado del los Caudales Medios Mensuales

Se procedió a realizar los cálculos necesarios para el desarrollo de la metodologı́a plan-

teada anteriormente a fin de generar las descargas medias mensuales en la zona del proyecto.
A continuación, los parámetros básicos utilizados para la micro cuenca Presa Acostambo.

Tabla 4.7
CAUDALES MEDIOS MENSUALES (Presa - Acostambo)

Fuente: Elaboración Propia

CAPÍTULO 4. ESTUDIOS PRELIMINARES 49

Tabla 4.8
CAUDALES MEDIOS MENSUALES (Presa - Acostambo)

DESCARGAS MEDIAS MENSUALES GENERADAS (m3/s) - CUENCA INTEGRAL

PERIODO: 1989 - 2014

Area de la cuenca: 6.75 Km2

AÑO Ene Feb Mar Abr May Jun Jul Ago Sep Oct Nov Dic Total
31 28 31 30 31 30 31 31 30 31 30 31 (m3/s)
1989 0.28 0.42 0.51 0.33 0.13 0.08 0.08 0.11 0.17 0.22 0.16 0.19 2.7
1990 0.34 0.22 0.17 0.15 0.08 0.14 0.07 0.10 0.11 0.22 0.30 0.31 2.2
1991 0.18 0.23 0.27 0.14 0.08 0.09 0.05 0.07 0.18 0.11 0.16 0.14 1.7
1992 0.19 0.23 0.25 0.11 0.09 0.06 0.05 0.06 0.08 0.18 0.13 0.15 1.6
1993 0.39 0.40 0.28 0.18 0.08 0.04 0.06 0.09 0.18 0.17 0.37 0.32 2.6
1994 0.41 0.56 0.43 0.27 0.09 0.04 0.02 0.03 0.10 0.18 0.12 0.24 2.5
1995 0.22 0.39 0.40 0.19 0.06 0.05 0.04 0.04 0.10 0.09 0.22 0.23 2.0
1996 0.27 0.37 0.26 0.18 0.06 0.08 0.08 0.04 0.07 0.12 0.17 0.25 1.9
1997 0.33 0.35 0.24 0.17 0.05 0.03 0.02 0.07 0.18 0.14 0.20 0.24 2.0
1998 0.32 0.41 0.21 0.08 0.04 0.03 0.07 0.05 0.05 0.23 0.24 0.24 2.0
1999 0.27 0.43 0.30 0.22 0.10 0.08 0.05 0.03 0.14 0.14 0.16 0.19 2.1
2000 0.30 0.40 0.31 0.17 0.15 0.06 0.07 0.06 0.07 0.22 0.10 0.20 2.1
2001 0.42 0.43 0.31 0.14 0.08 0.06 0.12 0.07 0.09 0.18 0.24 0.22 2.4
2002 0.16 0.48 0.37 0.17 0.09 0.04 0.11 0.09 0.12 0.18 0.18 0.28 2.3
2003 0.21 0.32 0.34 0.22 0.08 0.09 0.08 0.12 0.14 0.08 0.21 0.22 2.1
2004 0.18 0.21 0.24 0.14 0.08 0.12 0.06 0.06 0.11 0.17 0.29 0.30 2.0
2005 0.22 0.28 0.24 0.15 0.08 0.09 0.04 0.03 0.07 0.20 0.17 0.22 1.8
2006 0.34 0.35 0.29 0.18 0.06 0.05 0.07 0.08 0.09 0.18 0.28 0.26 2.2
2007 0.24 0.22 0.31 0.16 0.10 0.09 0.07 0.04 0.09 0.16 0.22 0.23 1.9
2008 0.32 0.35 0.21 0.13 0.07 0.06 0.03 0.02 0.12 0.23 0.13 0.21 1.9
2009 0.26 0.39 0.28 0.16 0.13 0.10 0.07 0.08 0.08 0.15 0.28 0.39 2.4
2010 0.50 0.36 0.31 0.12 0.06 0.05 0.07 0.08 0.05 0.17 0.14 0.28 2.2
2011 0.38 0.57 0.45 0.25 0.14 0.05 0.06 0.04 0.13 0.18 0.19 0.30 2.7
2012 0.23 0.61 0.25 0.32 0.16 0.07 0.04 0.04 0.08 0.12 0.20 0.44 2.6
2013 0.39 0.44 0.34 0.20 0.15 0.06 0.04 0.08 0.09 0.11 0.14 0.30 2.3
2014 0.35 0.33 0.31 0.26 0.10 0.09 0.08 0.07 0.12 0.16 0.22 0.30 2.4
MAX. 0.50 0.61 0.51 0.33 0.16 0.14 0.12 0.12 0.18 0.23 0.37 0.44 2.74
MED. 0.30 0.37 0.30 0.18 0.09 0.07 0.06 0.06 0.11 0.17 0.20 0.26 2.18
MIN. 0.16 0.21 0.17 0.08 0.04 0.03 0.02 0.02 0.05 0.08 0.10 0.14 1.57
D.EST 0.09 0.11 0.08 0.06 0.03 0.03 0.02 0.03 0.04 0.04 0.06 0.07 0.30

Fuente: Elaboración Propia

4.4. Balance Hı́drico

El Balance Hı́drico consiste en la aplicación del principio de la conservación de masa al

conjunto de una cuenca o a una cierta parte de ella definida por unas determinadas condiciones
de contorno durante un determinado perı́odo de tiempo en el que se realiza el balance.

OF ERT AHidrica − DEM AN DAT otal = BALAN CE.HIDRICO (4.64)

4.4.1. Oferta Hı́drica

La oferta hı́drica se refiere al caudal disponible en la cuenca, es decir los caudales

medios mensuales.

CAPÍTULO 4. ESTUDIOS PRELIMINARES 50

4.4.1.1. Según el Caudal Ecológico

Se entenderá como caudal ecológico al volumen de agua que se debe mantener en las
fuentes naturales de agua para la protección o conservación de los ecosistemas involucrados,
la estética del paisaje u otros aspectos de interés cientı́fico o cultural. Este caudal es una
demanda, pero por estar relacionada directamente con el caudal de oferta se analiza en esta
parte del cálculo, para evitar confusiones.

Los caudales ecológicos, son establecidas por el ANA (Autoridad Nacional del Agua).
Siguiendo los criterios establecidos, los caudales ecológicos corresponden al 15 % del caudal
medio mensual.

Finalmente los caudales de Oferta Hı́drica, serán los siguientes:

4.4.2. Demanda Hı́drica

Esta parte esta dedicada al cálculo de las demandas que deben ser atendidas con los
recursos provenientes del embalse, es decir con recursos hı́dricos regulados.

4.4.2.1. Demnada Agrı́cola

Para el presente proyecto las demandas serán para uso agrı́cola, de las cuales se
presentan los resultados de áreas que realmente pueden ser atendidos con la garantı́a hı́drica
mı́nima recomendable para un proyecto de riego, por lo tanto su formulación en términos
reales fue paralela a los cálculos de la simulación del balance hı́drico.

A. Patrón de Cultivo

Los cultivos más comunes sembrados en las localidades de Sinto y Esmeralda, obedecen
a los alimentos tradicionalmente consumidos por las poblaciones de los Andes y a las
condiciones medioambientales de la zona, dentro de ello tenemos: la cebada, haba, trigo,
papa, alafalfa, maiz choclo, y otros.

B. Coeficientes de Cultivo(Kc)

Para tener en cuenta los efectos de las caracterı́sticas del cultivo sobre sus necesidades
de agua, se presenta unos coeficientes de cultivo (Kc). Los valores de Kc en los que se
tienen en cuenta las caracterı́sticas de cultivo, momento de siembra, fases de desarrollo
vegetativo y las condiciones climáticas se aprecian a continuación:

CAPÍTULO 4. ESTUDIOS PRELIMINARES 51

Tabla 4.9
Kc obtenido en base a dos campañas

Fuente: Manual de uso consuntivo del agua para principales cultivos de los Andes Centrales
Peruanos.

Figura 4.16
Variación del coeficiente de cultivo en el tiempo

Fuente: Estudio FAO Riego y Drenaje.

C. Eficiencia de Riego

Para determinar la eficiencia de riego total se ha considerado uan eficiencia por asper-
sión, para lo cual se sabe que es de un 70 % para los cálculos de demanda del presente

CAPÍTULO 4. ESTUDIOS PRELIMINARES 52

 proyecto.

D. Jornada de Riego

Según las costumbres de los agricultores y de la zona se recomienda 12 h de riego por

presencia de heladas en las zonas altas y medias de la zona de cultivos.

E. Parcelamiento de las Áreas de Irrigación

Para el proyecto se cuenta con 249 has de riego para una primera etapa, y 98.75 has
en la segunda etapa, haciendo un toral de 347.5 has, comprendidos en la siembra de 11
cultivos los cuales son: la cebada, haba, trigo, papa, alafalfa, maiz choclo.

Tabla 4.10
Parcelamiento de las áreas de irrigación

Fuente: Elaboración Propia.

Tabla 4.11
Parcelamiento de las áreas de irrigación

Fuente: Elaboración Propia.

F. Cédula de Cultivos

Para el proyecto se está planteando dos campañas anuales de cultivo, por lo que los
benificiarios del proyecto podrán exportar 2 veces al año los alimentos de primera
necesidad a los principales mercados del Perú y asi mejorar su economı́a y calidad de
vida, esto hace eficiente y beneficioso al proyecto.

Tabla 4.12
Cédula de cultivos

Fuente: Elaboración Propia.

G. Cálculo de Evapotranspiración

No contamos con registros historicos de ETO, por lo que se estimó la evapotranspiración

por del método de HARGREAVES - SAMANI, los datos de entrada son: temperatura
máxima y minima, precipitación promedio mensual de la cuenca de estudio y su latitud.

Tabla 4.13
Evapotranspiración generada (mm) - Método HARGREAVES - SAMANI.

Fuente: Elaboración Propia.

CAPÍTULO 4. ESTUDIOS PRELIMINARES 53

Figura 4.17
Evapotranspiración generada (mm) - Método HARGREAVES - SAMANI.

Fuente: Elaboración Propia.

H. Cálculo de la Demanda Hı́drica (RIEGO)

Tabla 4.14
Cálculo de la Demanda Hı́drica

Fuente: Elaboración Propia.

Figura 4.18
Cálculo de la Demanda Hı́drica

Fuente: Elaboración Propia.

4.4.2.2. Pérdidas en el Embalse

Las pérdidas que se producirán en el embalse serán por evaporación e infiltración.

A. Caudal de Infiltración del Embalse

Para la determinar el caudal de infiltración se emplea el principio de Darcy que describe,

con base en experimentos de laboratorio, las caracterı́sticas del movimiento del agua a
través de un medio poroso. La expresión matemática de la Ley de Darcy es la siguiente:

Q = KiA (4.65)

Donde:

Q: Gasto, descarga o caudal en (m3/s).

K: Coeficiente de permeabilidad de Darcy, en función del material de la muestra

(m/s).

i: Gradiente hidráulica

A: Área de la sección transversal de la muestra, en m2.

Tabla 4.15
Cálculo de pérdidas por evaporación e infiltración en el embalse

Fuente: Elaboración Propia.

CAPÍTULO 4. ESTUDIOS PRELIMINARES 54

 Finalmente, la demanda hı́drica total que exige el proyecto se ha calculado sumando la
demanda por cultivos, la demanda por caudal ecológico y las pérdidas (evaporación e
infiltración) de embalse.

Tabla 4.16
Demanda Hı́drica Total del proyecto - Presa Sinto-Esmeralda

Fuente: Elaboración Propia.

4.4.3. Balance Hı́drico

Una vez determinados los caudales de demanda y caudales de oferta, se procede hacer
el balance hı́drico para ver si hay superávit o déficit de agua; esto nos ayudará a validar la
aprobación y ejecución del proyecto. Se puede observar que hay un deficit entre los meses de
mayo a septiembre, por lo que el proyecto de embalse es viable.

Figura 4.19
Balance Hı́drico del proyecto - Presa Sinto-Esmeralda

Fuente: Elaboración Propia.

4.5. Alturas Tı́picas de la Presa

4.5.1. NAMINO (Nivel de Aguas Mı́nimas de Operación)

La altura mı́nima del NAMINO (nivel de aguas mı́nimas de operación) lo calculamos

con un porcentaje del volumen util, la cual coincidirá con el nivel al que se va a ubicar la
toma que se proyectará. Esta altura asimismo dependerá del volumen mı́nimo requerido por
nuestro cultivo.
Figura 4.20
Distribución de agua en un embalse

Villamizar C. A.(1989)

CAPÍTULO 4. ESTUDIOS PRELIMINARES 55

 Consideramos un 12 % del volumen útil calculado.

12
VN AM IN O = ∗ 0.94M M C = 0.1128M M C
100

Entonces al ubicar en la altura de volumenes ubicación y la altura nos saldrı́a 2 metros

aproximandamente la cual consideramos esa altura.

4.5.1.1. Volumen Muerto de la Presa

Para el calculo determinaremos la altura que alcanza los sedimentos producto de

la erosión del suelo en la cuenca dentro del vaso de la presa, para determinar usaremos la
metodologia USLE.De acuerdo a la USLE, la tasa de pérdidas de suelo por erosión hı́drica,
T, es una función de: el poder erosivo de la lluvia, R, la erodabilidad de los suelos, K, la
cobertura vegetal, C, la práctica conservacionista, P y el factor combinado de la pendiente y
la longitud de la misma, LS.

T = RKCP LS (4.66)

El Factor R

El factor R depende de la intensidad de lluvia y de la lluvia total. El factor se estima

para cada una de las tormentas significativas ocurridas durante un registro de una
longitud aceptable, por lo general mayor de 10 años.

Figura 4.21
Mapa para valores del Factor R

Fuente: SENAMHI-DHI(2017)

Estimación del factor K

CAPÍTULO 4. ESTUDIOS PRELIMINARES 56

 Como se mencionó anteriormente, K es el factor de erodabilidad y se calcula en función
de: porcentaje de limo + arena muy fina (0.1-0.002 mm), porcentaje de arena, porcen-
taje de materia orgánica, ı́ndice de permeabilidad e ı́ndice de textura.Para el cálculo
usaremos un MAPA DE EROSABILIDAD DE PERÚ
Figura 4.22
Mapa para valores del Factor K

Fuente: SENAMHI-DHI(2017)

Factor topográfico LS

La USLE combina el efecto de la pendiente, S, con la longitud de flujo, L, en un factor

topográfico, LS. Para el cálculo usaremos el MAPA VALORES LS DEL PERÚ.

Figura 4.23
Mapa para valores del Factor topográfico LS

Fuente: SENAMHI-DHI(2017)

Factor de cultivo C

El factor de cultivo C es un ı́ndice del grado de cobertura y protección que presenta la

vegetación. Ese factor es constante para cultivos permanentes y variable para cultivos

CAPÍTULO 4. ESTUDIOS PRELIMINARES 57

 anuales.Para el cálculo usaremos el MAPA DE FACTOR DE CULTIVO PORPOR-
CIONADO POR LA SENAMHI-DHI(2017).

Figura 4.24
Mapa para valores del Factor de cultivo C

Fuente: SENAMHI-DHI(2017)

Factor de práctica conservacionista P

El factor de práctica conservacionista, P, refleja el efecto de las prácticas de conservación

de suelos que tienden a modificar la topografı́a, tales como: terrazas, cultivos en franjas
y curvas de nivel.Para el cálculo usaremos el MAPA DE FACTOR DE PRACTICA
CONSERVACIONISTA POR LA SENAMHI-DHI(2017).

Figura 4.25
Mapa para valores del Factor de práctica conservacionista P

Fuente: SENAMHI-DHI(2017)

CAPÍTULO 4. ESTUDIOS PRELIMINARES 58

Figura 4.26
Cálculo del volumen muerto del embalse - Método USLE

Fuente: Elaboración propia.

Al calcular el volumen muerto nos dio como resultado 0.10298 MMC la cual ubicaremos.
Entonces al ubicar en la altura de volumenes ubicación y la altura nos saldrá 3.5 metros
aproximandamente la cual consideramos esa altura.

4.5.2. NAMO (Nivel de Aguas Maximas de Operación)

El NAMO es el máximo nivel con que puede operar la presa para satisfacer las de-
mandas; volumen que se almacena entre el NAMO y el NAMIN o NAMINO se llama volumen
o capacidad útil y es con el que se satisfacen las demandas de agua.

La demanda hı́drica ya se calculó en capitulos anteriores, en cual se considera la de-

manda por riego, infiltración, evaporación y considerando el caudal ecológico. A continuación
se muestra el resumen de la demanda total del proyecto.

Tabla 4.17
Demanda Hı́drica Total del proyecto - Presa Sinto-Esmeralda

Fuente: Elaboración Propia.

Tambı́én ya se hizo el balance hı́drico, para el cuál se consideró todas las pérdidas y
el caudal ecológico en las demandas, los resultados se muestran a continuación.

Figura 4.27
Balance Hı́drico del proyecto - Presa Sinto-Esmeralda

Fuente: Elaboración Propia.

4.5.2.1. Volumen Útil de la Presa

Para estimar el volumen útil de la presa se usó la metodologı́a de los picos secuentes,
cuya teorı́a de esta metodologı́a se encuentra en el capitulo del marco teórico. Los resultados
se muestran a continuación.

Al realizar el calculo del volumen útil por la metodologia por PICO SECUENTE nos
sa lió un vol de 0.94 MMC. Entonces al ubicar en la altura de volumenes ubicación y
la altura nos saldrı́a 9 metros de altura.

CAPÍTULO 4. ESTUDIOS PRELIMINARES 59

Figura 4.28
Volumen Útil de Embalse (Met. Picos Secuentes) - Presa Sinto-Esmeralda

Fuente: Elaboración Propia.

Fuente: Elaboración Propia.

Fuente: Elaboración Propia.

Fuente: Elaboración Propia.

4.5.3. NAME (Nivel de Aguas Máximas Extraordinarias)

El NAME (nivel de aguas máximas extraordinarias) es el nivel maás alto que debe
alcanzar el agua en el vaso bajo cualquier condición. El volumen que queda entre este nivel y
el NAMO, llamado superalmacenamiento, sirve para controlar las avenidas que se presentan
cuando el nivel en el vaso está cercano al NAMO.

Al realizar el calculo del la altura del NAME por el método de piscina niveladora, nos
resultó 3 metros de altura, lo cual se considerará esa altura para el NAME.

A continucaicón se muestran los resultados del Transito de Caudal a través de un

embalse de detención utilizando el método de piscina niveladora.

Figura 4.29
Cálculo de El NAME - Método de piscina niveladora

Fuente: Elaboración propia.

Figura 4.30
Cálculo de El NAME - Método de piscina niveladora

Fuente: Elaboración propia.

4.5.4. Borde Libre de la Presa

Figura 4.31
Cálculo de el Borde Libre

Fuente: Elaboración propia.

CAPÍTULO 4. ESTUDIOS PRELIMINARES 60

 Al realizar el calculo del borde libre, considerando efectos del oleaje nos resultó 2
metros de altura para el borde libre.

4.5.5. Alltura Total de la Presa

Figura 4.32
ALTURA TOTAL de la Presa SINTO-ESMERALDA

Fuente: Elaboración propia.

4.6. Volumen de Almacenamiento de la Presa

4.6.1. Curva Altura - Volumen - Área

Estas curvas fueron elaboradas gracias al DEM obtenido de MINAM, para luego
obtener las curvas de nivel haciendo uso del software Qgis, para su posterior trattamiento en
el CIVIL3D, en el cual se dibujo con curvas de nivel cada medio metro (con el objetivo de
obtener mayor precisión) desde la cota 3960 m.s.n.m. hasta la cota 3984.5 m.s.n.m.

Figura 4.33
Curvas de nivel del embalse- CIVIL3D

Fuente: Elaboración propia.

Figura 4.34
Curvas de nivel del embalse- CIVIL3D

Fuente: Elaboración propia.

Para la curva Altura -Área, se midió el área que encerraba cada curva de nivel cortada
por el eje de la presa. Para la elaboración de la curva Altura - Volumen, se utilizaron las áreas
medidas anteriormente para la curva Altura - Área, asumiendo una pared vertical en el eje
de la presa. Finalmente con ayuda del software Excel se determinaron las lÍneas de tendencia
de ambas curvas ası́ como sus respectivas ecuaciones.

Figura 4.35
CURVA ALTURA - VOLUMEN - ÁREA de la Presa SINTO-ESMERALDA

Fuente: Elaboración propia.

CAPÍTULO 4. ESTUDIOS PRELIMINARES 61

Figura 4.36
CURVA ALTURA - VOLUMEN - ÁREA de la Presa SINTO-ESMERALDA

Fuente: Elaboración propia.

CAPÍTULO 4. ESTUDIOS PRELIMINARES 62

 Capı́tulo 5

DISEÑO DE PRESA DE CONCRETO

De acuerdo con los datos obtenidos en los estudios anteriores realizados, se proce-
derá a efectuar el diseño de la presa, se realizará el diseño por estabilidad (volcamiento y
deslizamiento). Se analizó en condición llena de la presa.

Para ello se cumplió con las siguientes condiciones:

A. Se debe asegurar que el nivel del agua en el almacenamiento no rebase al N.A.M.E.

(Nivel de Aguas Máximas Extraordinarias).

B. No es admisible que el concreto trabaje a esfuerzos de tensión. Para cumplir esta condi-
ción, la resultante de todas las fuerzas, arriba de un plano de corte horizontal cualquiera,
debe cortar a dicho plano dentro del tercio medio, a presa llena y a presa vacı́a.

C. La cortina debe ser segura al deslizamiento. El factor de seguridad contra el desliza-

miento (F.S.D.) en un plano de corte cualquiera, debe ser mayor de 1.5.

D. El coeficiente de seguridad contra el vuelco (F.S.V.) oscila entre 2 y 3.

Factor de seguridad al volteo (FSV)

El factor de seguridad al volteo, nos resultó FSV = 2.084, este valor se encuentra entre
2 y 3 entonces esta dentro del rango permisible, por lo que las dimensiones propuestas
son correctas.
Factor de seguridad al deslizamiento (FSD)
El factor de seguridad al deslizamiento, nos resultó FSD = 1.858, este valor es mayor
que 1.5 entonces se encuentra dentro del rango permisible, por lo que las dimensiones
propuestas son correctas.

A continuación se muestra el procedimeinto detallado del çálculo de los factores de

seguridad de la presa:

CAPÍTULO 5. DISEÑO DE PRESA DE CONCRETO 63

 Capı́tulo 6

CONCLUSIONES

En la realización del trabajo concluimos que para el diseño de un presa tenemos que
tener en consideración los diversos factores como el proceso constructivo, el costo del
proyecto, el impacto ambiental que produce ası́ mismo el tipo de presa a diseñar.

Luego de realizar los parámetros geomorfológicos de la cuenca de Sinto-Esmeralda, con

un DEM descargado de la plataforma del MINAM, se obtuvo un área de 56.85 Km2
ası́ mismo se obtuvo Relación de circularidad: 0.448 la cual nos muetra que la cuenca
es de una forma rectangular, con un potencial de creciente moderado y una producción
de cuadal alto.

Se ha logrado realizar el dimensionamiento hidráulico de los componentes principales

de la presa, tomando como demanda hı́drica principal la demanda agrı́cola, además de
considerar las pérdidas de caudal por infiltración, evaporación y caudal ecológico.

En la determinación de altura hidráulica se hizo el calculó de NAMIN considerando

el alamcenamiento de sedimentos con la metodologı́a de USLE y los mapas dados por
SENAMHI, NAMINO teniendo en cuenta que es para un proyecto de irrigación se
consideró el 12 % de del volumen util, NAMO hacinado el uso de la metodoglı́a de
PICOS SECUENTES Y NAME haciendo uso de la metodologı́a de transito de avenidas.
Todas estas haciendo una sumatoria total de 19.5 metros de altura.

En la determinación de la estabilidad de la presa, los valores para los factores de se-

guridad por volteo (FSV=2.084) y deslizamiento (FSD = 1.858) están dentro de lo
permitido por las normativas y condiciones (ambos son mayores que 1.5), las cuales
hacen a nuestra presa estable y segura para su construcción.

CAPÍTULO 6. CONCLUSIONES 64
 Capı́tulo 7

BIBLIOGRAFÍA

Gómez Hernández, A. y Antonio Mújica Garcı́a, L. (2009). CÁLCULO DE LAS FUER-

ZAS ACTUANTES EN LA CORTINA Y OBRA DE EXCEDENCIAS DE LA PRESA
MAXIMILIANO R. LÓPEZ..

Gomez, K (2022). Apuntes de clase del curso de Presas y Obras de Embalse Semestre
2022 -2. Escuela Profesional Ingenierı́a Civil. Universidad Nacional de Huancavelica.

Sabino Rojas, E., Felipe-Obando, O., y Lavado-Casimiro, W. (2017). Atlas de erosión de

suelos por regiones hidrológicas del Perú. Nota Técnica N◦ 002 SENAMHI-DHI-2017.

Eduardo Martı́nez Martı́n (20007), Diseño de Pequeñas Presa?, BUREAU OF RECLA-

MATION, traducción de la 3ra Edición Americana

Villón Béjar M. (2005), Diseño de Estructura Hidráulicas, Instituto Tecnológico de

Costa Rica, Escuela de ingeniera Agrı́cola, Lima-Perú

FAO (Fundación de la Naciones Unidas para la Alimentación y la Agricultura), Evapo-

transpiración del cultivo: guı́as para la determinación de los requerimientos de agua de
los cultivos. 322 p. Consultado el 19 de mayo del 2015.

ANA (2013). Plan Nacional de Recursos Hı́dricos Del Perú.

Vallarino, E. 2006. Tratado básico de Presas. 6 ed. Madrid, ES. Imprenta Visagrafic. v
1, tomo 1. 802 p.

CAPÍTULO 7. BIBLIOGRAFÍA 65
 Capı́tulo 8

ANEXOS

CAPÍTULO 8. ANEXOS 66