Science">
Elqui Yeye Pari Condori
Elqui Yeye Pari Condori
Elqui Yeye Pari Condori
AUTOR:
ASESOR:
LÍNEA DE INVESTIGACIÓN:
LIMA – PERÚ
2022
i
Dedicatoria
ii
Agradecimiento
iii
Índice de Contenidos
Carátula ............................................................................................................ i
Dedicatoria ........................................................................................................ ii
Resumen .......................................................................................................... ix
Abstract ............................................................................................................ x
I. INTRODUCCIÓN ......................................................................................... 11
3.6 Procedimientos....................................................................................... 39
V. DISCUSIÓN ............................................................................................... 65
VI. CONCLUSIONES...................................................................................... 76
REFERENCIAS .............................................................................................. 78
ANEXOS......................................................................................................... 79
iv
Índice de tablas
Figura N° 4: Mapa político del Perú, de la Región de Puno y la cuenca del rio Ilave
............................................................................................................................. 40
Figura N° 23: Menú de captura de imagen de la cuenca del río Ilave en GeoHMS
............................................................................................................................. 61
Figura N° 26: Menú de reconocimiento mini cuencas de la cuenca del río Ilave en
GeoHMS ............................................................................................................... 63
Figura N° 28: Menú de reportes de caudales máximos de la cuenca del río Ilave en
HecHMS ............................................................................................................... 64
vii
Gumbel ................................................................................................................. 68
viii
Resumen
ix
Abstract
The hydrological study is important, for purposes of design of hydraulic works and
in issues of prevention of risk management, the maximum flow rate that is recorded
during the unusual increase by extreme events, exceeding the normal average
values. The main objective is to determine the maximum design flows using
conventional probabilistic methods in the Ilave basin, Puno - 2022. The methods
used were the Gumbel distribution; Normal and Gamma (Pearson Type III). It has
been considered the maximum precipitation data to 24 hours to year from Ilave
station from El Collao province, Puno region – Perú country, the estimation of the
maximum design flows were processed by the soil conservation service (SCS)
method and the unit hydrograph. Whose results were achieved trough the Gumbel
distribution an approximate maximum design flow estimate of 141.6 m3/s, using de
normal distribution, an approximate maximum design flow estimate of 175.5 m3/s
and by means of the Gamma distribution (Pearson Tipo III) an estimation of
approximate maximum flow of 168.2 m3/s. Concluding that the normal distribution is
the one with the best fit to the maximum precipitation data for 24 hours a year, by
the maximum likelihood method and by showing less variability in the estimation by
return times (TR).
x
I. INTRODUCCIÓN
11
periodos de retorno, para que esto pueda determinarse mediante aplicaciones
informáticas y modelación hidráulica. (p.10).
Los estudios hidrológicos son importantes para estimar las características de los
caudales máximos y mínimos según se desee, en determinados casos el objetivo
es estimar los caudales máximos para la planificación de obras de hidrología como
protección de ríos. proporcionar medidas de protección para evitar futuras
inundaciones en los lugares medio y bajo de la cuenca. Para crear escenarios
actuales y futuros, se deben conocer los flujos máximos por tiempo de retorno de
acuerdo con las estaciones meteorológicas disponibles en su área. Actualmente
existe poca información de la calidad de los datos meteorológicos e hidrométricos,
en algunos casos se estiman caudales medios, los cuales se incluyen en este
trabajo con cálculos, los caudales máximos son con diferentes técnicas
desarrolladas y diferentes tiempos de ejecución, lo que nos da una idea de cómo se
dará el comportamiento bajo los escenarios de implementación de medidas para los
habitantes de lugares medio y bajo de la subcuenca de Ayaviri (p.14).
En este trabajo se plantea como problema general: ¿Cuánto son los caudales
máximos de diseño usando métodos convencionales probabilísticos en la
cuenca Ilave, Puno - 2022? Como también se plantean los problemas
específicos: ¿Cuánto es los Caudales máximos de diseño mediante la
distribución Gumbel en la cuenca Ilave, Puno - 2022?, ¿Cuánto es los Caudales
máximos de diseño mediante la distribución normal en la cuenca Ilave, Puno -
2022? Y ¿Cuánto es los Caudales máximos de diseño mediante Gamma
(Pearson Tipo III) en la cuenca Ilave, Puno - 2022?
13
II. MARCO TEÓRICO
14
corresponden a la zona de Ambato . catastro cantonal, el muestreo es no
probabilístico según juicio o consideración de expertos, cuyos resultados consideran
que no existe volumen de vaciado para un vehículo durante todo el período de
lavado, de acuerdo a los datos obtenidos por el método de volumen existen valores
que son aproximados al volumen de vaciado, en cuyo caso se decide hacer el
cálculo con una regla de tres simple, tomando como referencia el tiempo total de
lavado de cada vehículo, lo que da como volumen de referencia, los caudales por
el método C es 6.98 litros/min, para el método 1 es de 1.56 lt/min, para el método 2
es de 2.32 lt/min, el caudal promedio diario con un valor de 2.32 lt/min, se concluye
que el caudal de diseño calculado es de 2.32 lt/min, pero el caudal de diseño en la
planta de tratamiento se utiliza un caudal de 3 l/min, que sirve como zona de
seguridad para el limpiador, se explicaron los procesos que requiere el agua que
ingresa a la red de alcantarillado.
15
los cuales se relacionaron con la aparición de precipitación en el CRN en los meses
anteriores, entre los principales resultados se reporta que el factor clave es la
cantidad de precipitación en la estación fría de mayo a septiembre, debido a que en
el trimestre JJA es la más importante, en este sentido pueden ocurrir caudales
extremos debido a dos procesos de gran escala que interactúan entre sí, además
de otros factores aún por explicar, y que provocan perturbaciones en las
precipitaciones invernales..
16
estudio es el río Pomabamba, se encuentra dentro de las coordenadas: -8.81318 S,
-77. 692 W y -8.8260 S, -77. 5713 W, zona urbana de Pomabamba. el resultado es
una franja con un ancho de 100 a 130 metros y una longitud de 1.967 metros, que
es la extensión actual del casco urbano de la ciudad de Pomabamba. “Los máximos
caudales pico del río Pomabamba tienen efectos proporcionales en las planicies de
inundación que pueden ser modelados hidráulicamente y evaluados con
herramientas informáticas”, se obtuv0o una serie histórica de 41 años de
precipitación máxima entre 1980-2020, durante el trabajo se determinó que la
distribución de probabilidad lognormal 2P es la que mejor se ajusta de acuerdo a
las técnicas utilizadas; con el diagrama y el método Smirnov - Kolmogorov, donde
el valor máximo ⧍ = 0,13329; Así, las mayores intensidades de diseño se obtuvieron
mediante la metodología de Dick Peschke, la cual se utilizó para determinar la curva
IDF para obtener valores de intensidad en función del tiempo, concluyendo en el
trabajo que TR en 20, 50 y 100 años se obtuvieron los caudales de 135,22, 156, 8
y 17 ,98 m3/s, que con ayuda de la modelación hidráulica se pudo calcular las áreas
de inundación en 14 .683.38 m2, 21,869.71 m2 y 21,869.71 m2..
Como objetivo estimar un modelo para análisis redes de agua potable utilizando
demandas con un enfoque probabilístico/difuso cuya finalidad será estimar la
demanda máxima, registra como metodología el uso de un modelo basado en la
simulación difusa de Monte Carlo para análisis de redes de distribución del agua
potable, como diseño de la investigación se aplica el modelo a un proyecto del
ámbito nacional, de la siguiente manera: se revisarán las características de la
población considerada en el proyecto, y la demanda considerada en el diseño se
caracterizará el requerimiento del agua potable de la población en base a sus
características, utilizando un enfoque probabilístico, se definirá la correlación entre
demandas nodales como una variable difusa, se realizará la simulación de la red
utilizando el método de análisis propuesto, se determinarán las variaciones de los
principales parámetros hidráulicos de la red (velocidades, presiones y caudales), se
realizará un análisis comparativo entre los resultados del análisis determinista y el
análisis probabilístico para evaluar la variación entre los caudales de diseño, cuya
unidad de estudio: proyectos de mejoras y mayor cobertura del sistema de
distribución de agua potable de la localidad de Pucusana, la población y muestra se
estima que la población de diseño oscila entre 10,000 y 15,000 habitantes, para la
recolección de los datos se obtendrá del expediente técnico, el cual se encuentra
publicado en la web de SEDAPAL, arribando a resultados como según el análisis
determinista, cerca del 95% de las tuberías conducen un caudal de menos de 5
litros por segundo. Asimismo, se muestra que el límite superior es ligeramente
mayor que 51 l/s. Precisamente, el caudal de la tubería principal, P-4, es de 51.64
l/s, en lo referente a la velocidad el valor estimado al intervalo de confianza máximo
está por debajo de 0.94 m/s, la velocidad máxima con el enfoque determinista es de
1.01 m/s (línea vertical verde). Puede notarse la coherencia entre ambos resultados,
se generaron también las curvas de presiones para los nodos en análisis, los
resultados para el nodo J-236, que es el que soporta mayores presiones, es igual a
50.12 m, y se encuentra por debajo del rango de presiones máximas esperadas
18
según el análisis probabilístico, puede notarse que las presiones máximas
esperadas son mayores, lo que concuerda con los menores caudales requeridos en
general, llegando a concluir que entre las ventajas más importantes que presenta el
modelo es que permite representar la demanda considerando su naturaleza
aleatoria, utilizando factores de demanda validados para la zona de estudio, se
pueden optimizar los caudales resultantes del análisis hidráulico, y como
consecuencia, mejorar el diseño de redes de abastecimiento de agua, asimismo, el
método propuesto permite establecer niveles de confianza para el diseño, al
basarse en un análisis probabilístico y posibilístico (difuso).
Como objetivo de este trabajo es utilizar modelos probabilísticos para estimar los
caudales máximos de crecida de diferentes periodos de retorno en la cuenca de
interés, considerando los métodos de recolección de información cartográfica en
formato Shapefile de MINEDU del Perú, cartas sobre el estudio. El área está en una
escala de 1:100000 m, se recolectaron los datos hidrometeorológicos necesarios
para el desarrollo de este estudio para estimar los datos máximos de precipitación
y caudales, y se realizó análisis de la frecuencia máxima de eventos hidrológicos.
Para el caudal máximo diario en la estación Puente Magdalena y las lluvias máximas
en 24 hrs en las estaciones de la cuenca en estudio se utilizó el siguiente proceder
de: elegir del máximo anual: pluviometría máxima precipitación diaria. estaciones y
el caudal máximo de la estación Puente Magdalena, determinando la idoneidad de
las distribuciones de frecuencias más frecuentes (distribución normal, protocolo de
distribución) por el método de Kolmogorov-Smirnov. , Pearson III Distribution, Log
Pearson III Distribution y Gumbel Distribution) para obtener la distribución que bien
se ajusta a los datos históricos, proporcionando datos únicos de lluvia máxima de
24 horas para Lachaqui, Canta, Huamantanga, Huaros y Carampoma, sometidos al
análisis de datos sospechosos, el análisis estadístico se realizó utilizando
distribuciones normales, distribuciones Gumbel, Log Normal, Pearson III y
LogPearson Tipo III, se realizó la prueba de bondad de ajuste de Kolgomorov
Smirnov y se determinaron las distribuciones de buen ajuste, la distribución log-
Normal más cercana, que se determinó que la superficie de la cuenca hidrográfica
es de 1339, 9 km2 y la circunferencia de 282,83 km con un factor de forma de F =
0,16; Por lo tanto, se llega a la conclusión de que es una cuenca inmensa y alargada,
y se encontró que las lluvias de las estaciones pluviométricas de Canta, Lachaqui,
Huaros y Huamantanga coinciden con la curva de distribución LogNormal,
LogPearson III, LogNormal y Pearson III.
20
Para el análisis de precipitaciones máximas diarias y caudales máximos
Carmona (2020) formula en el trabajo de investigación presentado:
Como bases teóricas se tiene a Mejía (2012) que escribe sobre las distribuciones
de probabilidad manifestando que:
−𝛼𝛼(𝑥𝑥−𝛽𝛽)
𝐹𝐹(𝑥𝑥) = 𝑒𝑒 −𝑒𝑒 (1)
Para:
dF ( x)
f ( x) = (2)
dx
1 −α ( x − β )
f ( x) = e −α ( x − β ) − e (3)
α
y = α (x − β ) (4)
22
De la cual, la función de densidad minimizada Gumbel es:
−y
g ( y) = e− y −e (5)
−y
G( y) = e−e (6)
F ( x) = G ( y ) (7)
y la función:
y = α (x − β ) (8)
y
x=β+ (9)
α
Media:
c
E ( x) = X = β + (10)
α
1 1 1
c = lim n → ∞ 1 + + + ........ + − ln n
2 3 n (11)
c = 0.5772156649
0.5772156649
X =β+ (12)
α
Variancia:
23
π2
E [( x − E ( x)) ] = S = 2 2 2
(13)
α *6
De la cual se reporta:
1.2825
α=
S (14)
0.57721
β=X− = X − 0.45S
S (15)
1 n
β = ln n (16)
α ( −αx i )
∑
e
i =1
∑ [e ] ∑x
n n
( −αx i )
* xi i
1
+ i =1
− i =1
=0 (17)
α n
n
∑e
i =1
( −αx i )
Donde:
α y β son parámetros
n = muestra de datos.
xi = ésimo dato
2
1 X-µ
1 −
σ (18)
f (X) = e 2
σ 2π
24
1 n
µ= ∑ Xi
n i =1
(19)
1/ 2
1 n
σ = ∑ (X i - µ )2 (20)
n - 1 i =1
Donde:
( x −δ )
1 x − δ β1 −1 −
f ( x) = ( ) e α
(21)
αΓ( β 1 ) α
Para:
δ ≤x<∞
−∞ <δ < ∞
0 <α < ∞
0 < β1 < ∞
x ( x −δ ) β1 −1
1 − x −δ
αΓ( β 1 ) ∫0
F ( x) = e α
* dx (22)
α
Donde se denota:
25
x = Variable aleatoria de Pearson tipo III.
α = Parámetro de escala
β1 = Parámetro de forma
x −δ
y= (23)
α
y
y β1 −1e − y
F ( y) = ∫ dy (24)
0
Γ( β 1 )
Como parámetro β1, dicha variable aleatoria tienen origen en y=0 o x=δ.
Promedio:
X = δ + α * β1 (25)
Variancia:
S 2 = α 2 * β1 (26)
Sesgo:
2
Cs = g = (27)
β1
26
4
β1 = 2 (28)
Cs
α = Cs S / 2 (29)
2S
δ =X− (30)
Cs
N 2M 3
Cs = g = (31)
( N − 1)( N − 2) S 3
Dónde:
M3 =
∑ (x i − X )3
(32)
N
S =
∑ (x i − X )2
(33)
N −1
X =
∑x i
(34)
N
27
1 (35)
𝐹𝐹 = 1 −
𝑇𝑇
Para tener en cuenta el rango de retorno se usa la ecuación (36) para visualizar las
recurrencias, se muestra del siguiente modo:
𝑛𝑛 + 1
𝑇𝑇 = (36)
𝑚𝑚
Donde:
M = valor de orden
Respecto al tiempo medio, Tucci (2012) menciona que: “es el tiempo medio en
los años entre hechos o hechos iguales o superiores a un determinado rango, a
este tiempo medio se le conoce como tiempo de retorno o periodo de retorno”
(p.34).
Por su parte para los estudios de los distintos fenómenos hidrológicos Villón
(2004) afirma que:
29
𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉 𝑑𝑑𝑑𝑑 𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒í𝑎𝑎 𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠 𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡
𝐶𝐶 = (37)
𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉𝑉 𝑑𝑑𝑑𝑑 𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝ó𝑛𝑛 𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡
En estudios de los procesos hidrológicos Jain y Shain (2019) confirman que son:
30
hidrológico. De manera similar, la circulación del agua en la atmósfera o en la
superficie del océano que tiene un papel importante en la determinación de la
corteza terrestre, que crea patrones de los vegetales y controla los eventos
atmosféricos. El ciclo hidrológico puede analizarse como una máquina natural
(Figura 1), un sistema de operación de forma continua, destilado y bombeado. Un
grupo de ciclos naturales en los que una gran cantidad de materia se mueve
constantemente a través de la Tierra, incluido el ciclo hidrológico, el ciclo del
carbono, el ciclo del nitrógeno y varios ciclos biogeoquímicos. El ciclo del agua, es
la circulación del agua, que contiene cambios en su estado físico del agua entre los
estados líquido, sólido y gaseoso. En este ciclo se verifica el intercambio constante
de agua entre la atmósfera, la tierra (se incluye aguas superficiales y dentro de la
tierra), los océanos, los mares y los organismos de la naturaleza (p.54).
Toda investigación sobre el agua está relacionada con alguna parte del ciclo
hidrológico. En la disposición y calidad de agua están ligadas en bastantes lugares
al crecimiento demográfico de las poblaciones humanas, así como la actividad
humana cambia y modifica las características del ciclo hidrológico. Los estudios
hidrológicos sirven para muchos propósitos importantes; calcular el volumen de
captación, simular la respuesta de la captación a las cuencas hidrográficas,
confirmar los datos observados y predecir inundaciones. Los modelos son
adecuados para tareas diferentes. El estudio ayuda en tomar una decisión,
31
especialmente cuando el entendimiento del sistema es incompleta, falta información
o hay muchas posibilidades y no es posible entablar con un sistema de vida real
para entender las consecuencias de una decisión. Las plantillas son importantes en
estos casos (p.33).
33
III. METODOLOGÍA
Por su parte Arias (2012) indica en su trabajo con respecto al método científico
como:
34
Tiene como objetivo recopilar información teniendo en cuenta los conocimientos
teóricos y prácticos aprendidos durante la carrera; emprender acciones directas en
resolver los problemas que se investigan. La metodología cuantitativa consiste en
un conjunto de técnicas que pueden utilizarse para acercarse a una población y
estimar sus características, discusiones, referencias o demandas (p.67).
La exploración de la cuenca está limitada por las cabeceras y los puntos más altos
conocidos como cuencas hidrográficas. Entre sus características geomorfológicas
se analizan el área, longitud y circunferencia, pendiente de la cuenca, entre otras;
La precipitación total de 24 horas registrada para el evento es registrada por una
estación meteorológica, se calculan los valores máximos mensuales y la serie
histórica de montos máximos mensuales. El tiempo de retorno (TR) es una
estimación de la probabilidad de que ocurra un determinado evento durante un
tiempo determinado, también es el tiempo en el que la probabilidad del evento se
distribuye equitativamente entre los períodos que componen ese tiempo. Esto
también se llama el período de recurrencia. Los caudales máximos del proyecto son
operados por las dimensiones de las características geomorfológicas, que se
convierten en métodos de cálculo; estos a su vez se dividen en indicadores como
área, pendiente y frecuencia de drenaje. Además, los caudales de diseño más altos
35
están impulsados por la lluvia de 2 horas, que son métodos de cálculo; estos, a la
vez, se dividen en indicadores como alto, medio y bajo. Finalmente, los caudales
máximos de diseño se implementan a través de la dimensión del tiempo de retorno,
que son técnicas de cálculo; estos a su vez se dividen en indicadores como 2, 10,
50, 100 y 200 años de periodo de retorno (p.67).
Con lo anticipado, con el grupo de los datos de caudales de las estaciones que
se encuentran en la cuenca de Ilave indica Arias (2012) en su trabajo de
investigación que:
37
3.4 Técnicas e instrumentos de recolección de datos
En la Tabla 1 se consideran parámetros de validez para nuestro estudio.
Rango Magnitud
0.53 a menos Validez nula
0.54 a 0.65 Validez baja
0.60 a 0.65 Válida
0.66 a 0.71 Muy válida
0.72 a 0.99 Excelente válidez
1.00 Validez perfecta
Fuente: Oseda et al. (2015)
Título Apellidos y
DNI CIP Dictamen Evaluación
N° profesional nombres
1 Ingeniero Civil Deza Ramos, Darwin 42477401 128272 0.954
2 Ingeniero Civil Coyla Idme, Leonel 02146851 46961 0.927 0.933
Quispe Quea, Juan
3 Ingeniero Civil Pablo 01345604 87217 0.919
Para lo cual, el instrumento guiado por los parámetros de validez que considera
la Tabla 1. Para la confiabilidad en el estudio realizado es confiable y según
Hernández et al (2014):
38
resultados se utilicen en múltiples aplicaciones e interpretaciones. varias formas de
medir la confiabilidad con puntajes entre 0 y 1: el coeficiente alfa de Cronbach, la
medida de estabilidad, el método dividido y el coeficiente KR-20 de Kuder, que es
similar al coeficiente alfa pero se puede usar cuando la variable es dicotómica (p.87).
39
IV. RESULTADOS
Nombre de la tesis
Determinación de caudales máximos de diseño mediante métodos
convencionales probabilísticos en la cuenca Ilave, Puno – 2022
Ubicación política
La presente investigación se realizó en el distrito de Ilave, provincia de El Collao
en la Región de Puno - Perú.
Figura N° 4: Mapa político del Perú, de la Región de Puno y la cuenca del rio Ilave
Limites
Norte : Con la provincia de Puno.
Sur : Con la provincia de Chucuito.
Este : Con el lago Titicaca.
Oeste : Con la región de Moquegua.
40
Características morfológicas de la cuenca del río ilave
Ubicación geográfica
En el sistema de coordenadas UTM WGS84, Zona 19 sur, se localiza entre las
coordenadas: Este: 352,353 – 452,052 y Norte: 8’104,770 – 8’248,751.
En la Figura 6 se muestra una vista panorámica del rio Ilave cercano a la ciudad
de Ilave de la provincia de El Collao de la región de Puno – Perú. La estación
pluviométrica de Ilave se encuentra a mitad de la imagen mostrada. El momento
de tomar la imagen es en la estación de otoño (mes de agosto 2022) donde
existe escasez de precipitaciones pluviales.
años 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78
precipMax
17.81 20.5 20.6 21.16 21.5 23 23.2 24.5 25.8 26.1 26.2 27 27.4 27.6 27.91
(mm/dia/anual)
años 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93
precipMax
28.2 28.6 28.7 28.8 29 29.1 29.1 29.1 29.3 30.1 30.4 30.5 30.7 30.72 30.8
(mm/dia/anual)
años 94 95 96 97 98 99 00 01 02 03 04 05 06 07 08
precipMax
32.2 32.52 32.7 33 33.2 33.3 33.4 33.6 34 34 34.2 35.9 36 36.1 36.2
(mm/dia/anual)
años 09 10 11 12 13 14 15 16 17 18 19 20 21
precipMax
36.8 36.9 37 38.1 39.5 39.7 40.2 40.5 41 43.7 44 45 46
(mm/dia/anual)
Fuente: SENAMHI – Puno, elaboración del investigador
42
2021, un total de 58 años. De estos años registrados, en 1964 se tuvo la más
baja precipitación, como es el de 17.81mm, mientras que en el año de 2021 se
tuvo la más alta precipitación de 46.00mm. Cabe indicar que entre los años de
1977 a 1980 no se registraron las precipitaciones a nivel de la región de Puno,
por lo que se han estimado los valores para estos años a través de modelamiento
matemático lineal simple de entre los años anteriores.
25.00
20.00
15.00
10.00
5.00
0.00
1964
1966
1968
1970
1972
1974
1976
1978
1980
1982
1984
1986
1988
1990
1992
1994
1996
1998
2000
2002
2004
2006
2008
2010
2012
2014
2016
2018
2020
2022
años
43
Histórico de precipitaciones TOTALES anuales
(mm/año)
250.00
totales en mm/año
200.00
150.00
100.00
50.00
0.00
1964
1966
1968
1970
1972
1974
1976
1978
1980
1982
1984
1986
1988
1990
1992
1994
1996
1998
2000
2002
2004
2006
2008
2010
2012
2014
2016
2018
2020
2022
años
44
Hietograma de precipitación promedio
30.0
25.4
precipitación en mm/mes/anual
25.0 22.3
19.7 19.1
20.0
0.0
ene feb mar abr may jun jul ago set oct nov dic
meses (durante 58 años)
45
Hietograma de precipitación MAXIMAS
50.0 46.0
43.7
45.0 41.0 39.7
39.5
40.0
precipitación en mm
46
Figura N° 11: Menú de presentación del software HyFranPlus
47
Figura N° 12: Menú de ingreso y edición de datos en HyFranPlus
48
En la Figura 13, se muestra el menú de reportes de los resultados de las
estadísticas básicas en el software HyFranPlus, dentro de los que destacan el
número de datos ingresados, el mínimo, el máximo, el promedio, la desviación
estándar, la mediana, el coeficiente de variación, el coeficiente de sesgo y el
coeficiente de agrupamiento.
Objetivo específico 1:
Determinar los caudales máximos de diseño mediante la distribución Gumbel
U 28.57104
Alpha 6.086155
Ajuste con momentos ordinarios: como el delta teórico 0.0744, es menor que el
delta tabular 0.1786. Los datos se ajustan a la distribución Gumbel, con un nivel
49
de significación del 5%, cuyos parámetros de la distribución Gumbel con
momentos ordinarios:
Parámetro de posición (µ)= 28.8568
Parámetro de escala (alfa)= 5.0307
50
Figura N° 15: Menú de la prueba de bondad a la distribución Gumbel en HyFranPlus
caudal max
TR probabilidad Precipitación DesvStandar interval conf(95%) Diferencia m3/s
200 0.995 60.8 3.69 53.6 68 14.4 3305.7
100 0.99 56.6 3.26 50.2 63 12.8 3002.8
75 0.9867 54.8 3.09 48.8 60.9 12.1 2679.5
50 0.98 52.3 2.84 46.8 57.9 11.1 2255.7
20 0.95 46.6 2.28 42.2 51.1 8.9 1406.5
10 0.9 42.3 1.86 38.6 45.9 7.3 911.8
5 0.8 37.7 1.44 34.9 40.5 5.6
3 0.6667 34.1 1.15 31.8 36.3 4.5
2 0.5 30.8 0.936 29 32.6 3.6 141.6
Mu 31.760345
Sigma 6.451328
53
Media lineal (Xl)= 31.7605
Desviación estándar lineal (Sl)= 6.514
54
Figura N° 18: Menú de la prueba de bondad a la distribución Normal en HyFranPlus
caudal max
TR Probabilidad Precipitación DesvStandar interval conf(95%) Diferencia m3/s
200 0.995 48.4 1.77 44.9 51.9 7 1540.5
100 0.99 46.8 1.64 43.6 50 6.4 1433.4
75 0.9867 46.1 1.58 43 49.2 6.2 1340.2
50 0.98 45 1.5 42.1 48 5.9 1203.8
20 0.95 42.4 1.31 39.8 44.9 5.1 922.0
10 0.9 40 1.15 37.8 42.3 4.5 690.1
5 0.8 37.2 0.988 35.3 39.1 3.8
3 0.6667 34.5 0.886 32.8 36.3 3.5
2 0.5 31.8 0.847 30.1 33.4 3.3 175.5
55
En la Tabla 8 se detallan los resultados de las probabilidades en tiempos de
retorno de las precipitaciones máximas de 24 horas anuales con sus respectivos
intervalos de confianza, procesados a través del software HyFran Plus; además
de registrar en la última columna las estimaciones de los caudales máximos en
m3/s obtenidos a partir de las precipitaciones estimadas con la distribución
Normal, luego de usar el software Hec HMS 4.10.
56
encuentran dentro de los límites inferior y superior (color azul)
Objetivo específico 3:
Estimar los caudales máximos de diseño mediante la distribución Gamma
(Pearson Tipo III)
Alpha 1.605173
Lambda 105.553571
M -33.998018
57
parámetros, con un nivel de significación del 5%, los 2 parámetros de la
distribución Gamma con momentos ordinarios:
Parámetro de forma (gamma)= 23.9355
Parámetro de escala (beta)= 1.3269
Figura N° 20: Ajuste de precipitaciones máximas por la distribución Gamma (Pearson Tipo III)
58
Figura N° 21: Menú de la prueba de bondad a la distribución Pearson tipo III en HyFranPlus
Tabla N° 10: Probabilidades de las precipitaciones para distintos tiempos de retorno, intervalos
de confianza y caudales máximos estimados con la distribución Gamma (Pearson Tipo III)
59
de las precipitaciones máximas de 24 horas anuales con sus respectivos
intervalos de confianza, procesados a través del software HyFran Plus; además
de registrar en la última columna las estimaciones de los caudales máximos en
m3/s obtenidos a partir de las precipitaciones estimadas con la distribución
Gamma (Pearson Tipo III), luego de usar el software Hec HMS 4.10.
Las precipitaciones máximas de 24 horas estimadas por el modelo de
distribución de probabilidades Gamma (Pearson Tipo III) para tiempos de retorno
de 2, 10, 50 y 100 años se estiman en: 31.6, 40.1, 45.6 y 47.6 mm máximos de
24 horas, los cuales a su vez estiman 168.2, 699.2, 1275.3 y 1533.4 m3/s de
caudales máximos de diseño para la cuenca Ilave cerca al punto más crítico que
es el puente internacional de Ilave.
Figura N° 22: Tendencia de precipitaciones máximas por la distribución Gamma (Pearson Tipo
III)
Figura N° 24: Menú de reconocimiento de imagen de la cuenca del río Ilave en GeoHMS
61
En la Figura 24, se muestra el menú de reconocimiento de bloques en el software
GeoHMS de la cuenca del río Ilave a través de las cotas de nivel reconocidas en
el programa, siendo el segundo paso para determinar las parámetros mínimos
para el cálculo de caudales máximos, como son las longitudes, los bloques las
direcciones y las áreas de la cuenca Ilave.
Figura N° 25: Menú de reconocimiento de bloques de la cuenca del río Ilave en GeoHMS
62
Figura N° 26: Menú de reconocimiento mini cuencas de la cuenca del río Ilave en GeoHMS
Figura N° 27: Menú de reconocimiento de imagen de la cuenca del río Ilave en HecHMS
63
En la Figura 27, se muestra el menú de reconocimiento de la imagen en el
software HecHMS de la cuenca del río Ilave a través del software GeoHMS con
todas las características de reconocimiento de mini cuencas que se encontraron
como son las catorce mini cuencas en la cuenca del río Ilave. El software se
encarga de estimar a través de métodos matemáticos las aproximaciones de
caudales máximos a partir de las precipitaciones pluviales en esta cuenca. Para
ello se reportaron caudales máximos para periodos de retorno de 2, 10, 50 y 100
años, considerando distintas distribuciones de probabilidad.
Figura N° 28: Menú de reportes de caudales máximos de la cuenca del río Ilave en HecHMS
64
V. DISCUSIÓN
Discusión 1
65
Estimación de precipitaciones máximas con
la distribución Gumbel
55 52.865
50
45
50 100 200
Periodos de retorno en años
315.893
320
300 293.135
280 270.305
260
240
50 100 200
Periodos de retorno en años
66
el modelo de distribución de probabilidades de Gumbel para tiempos de retorno
de 50, 100 y 200 años se estiman 52.3, 56.6 y 60.8 mm máximos de 24 horas,
los cuales a su vez estiman 2255.7, 3002.8 y 3305.7 m3/s de caudales máximos
de diseño para la cuenca Ilave cerca al punto más crítico que es el puente
internacional de Ilave.
60
58 56.6
56
54 52.3
52
50
48
50 100 200
Periodos de retorno en años
67
En la figura 31 se muestra un crecimiento moderado de las estimaciones de
precipitaciones máximas en mm (52.30, 56.6 y 60.80) para los periodos de
retorno de 50, 100 y 200 años con la distribución de probabilidad de Gumbel.
3000 2255.7
2000
1000
0
50 100 200
Periodos de retorno en años
Por lo que existe similitud para el tiempo de retorno de 50 años para la estimación
de las precipitaciones máximas de 24 horas estimadas por el modelo de
distribución de probabilidades de Gumbel (52.86mm frente 52.30mm), mientras
que para el tiempo de retorno de 100 años se muestra también similitud en la
estimación de las precipitaciones máximas de 24 horas (57.33mm frente a
56,60mm). Pero para la estimación de caudales máximo de diseño a partir de la
distribución Gumbel existe fuerte discrepancia, tal es asi que para el tiempo de
retorno de 50 años en el trabajo de Mamani es de 270.31 m3/s, mientras que el
nuestro para el mismo tiempo de retorno es de 2255.7 m3/s. De igual modo para
el tiempo de retorno de 100 años en su trabajo muestra un caudal de retorno de
293.14 m3/s y en el nuestro para el mismo tiempo de retorno es de 3002.8 m3/s.
68
Discusión 2
Tabla N° 13: Estimaciones de precipitaciones máximas de 24 horas y caudales máximos según periodos
de retorno por la distribución Normal
42.88
43
42 41.66
41 40.32
40
39
50 100 200
Periodos de retorno en años
69
Estimación de caudales máximos con la
distribución Normal
225
219.25
210 206.161
205
200
195
50 100 200
Periodos de retorno en años
70
Estimación de precipitaciones máximas
con la distribución Normal
49 48.4
1540.5
1433.4
1500
1203.8
1000
500
0
50 100 200
Periodos de retorno en años
71
Se muestra cierta similitud para el tiempo de retorno de 50 años para la
estimación de las precipitaciones máximas de 24 horas estimadas por el modelo
de distribución de probabilidades de Normal (40.32mm frente 45.00mm),
mientras que para el tiempo de retorno de 100 años se muestra también similitud
en la estimación de las precipitaciones máximas de 24 horas (41.66mm frente a
46.8mm). Pero para la estimación de caudales máximo de diseño a partir de la
distribución Normal existe fuerte discrepancia, tal es asi que para el tiempo de
retorno de 50 años en el trabajo de Mamani es de 206.16 m3/s, mientras que el
nuestro trabajo para el mismo tiempo de retorno es de 1203.8 m3/s. De igual
modo para el tiempo de retorno de 100 años en su trabajo muestra un caudal de
retorno de 213.012 m3/s y en el nuestro para el mismo tiempo de retorno es de
1433.4 m3/s.
Discusión 3
72
Estimación de precipitaciones máximas
con la distribución Pearson Tipo III
64 61.781
310
300 293.135
290
280 270.305
270
260
250
240
50 100 200
Periodos de retorno en años
49
48 47.6
47
46 45.6
45
44
43
50 100 200
Periodos de retorno en años
74
de 50, 100 y 200 años con la distribución de probabilidad de Pearson Tipo III.
1000
500
0
50 100 200
Periodos de retorno en años
Conclusión 1
Conclusión 2
Conclusión 3
76
VII. RECOMENDACIONES
Recomendación 1
Recomendación 2
Recomendación 3
77
REFERENCIAS
Álvarez, AJ; Orduña, LM; Montañés, CG. 2014. Análisis y selección de modelos
estadísticos para el ajuste de la ley de frecuencia de caudales máximos
anuales en España. Revista Digital del Cedex (174): 5 - 31 p.
Chen, L; Wang, L. 2018. Recent advance in earth observation big data for
hydrology. Big Earth Data 2(1):86 - 107. Disponible en
78
https://doi.org/10.1080/20964471.2018.1435072
Chow, VT; Maidment, DR; Larry, WM. 1994. Hidrología aplicada. S.A., MI.
Bogotá, CO, Editorial McGrawHill Interamericana S.A. 584 p.
Fattorelli, S; Fernández, PC. 2011. Diseño Hidrológico. Paris, FR, World Water
Assessment Programme (WWAP). 531 p. (2)
Haan, CT; Barfield, BJ; Hayes, JC. 1994a. Design hydrology and sedimentology
for small catchments. Elsevier. p.
Khan, MA. 2015. Flood Scour for Bridges and Highways: Prevention and Control
79
of Soil Erosion. New York, US, McGraw - Hill Education. 384 p. Disponible
en
https://www.accessengineeringlibrary.com/content/book/9780071825078
Laura, MG. 2014. Modelo regional aplicado a las maximas avenidas en las
cuencas de los rios: Ilave, Callacame y Zapatilla. Pregrado. Puno, PE,
Universidad Nacional del Altiplano (UNA). 113 p.
Mackenzie, LD. 2019. Water and Wastewater Engineering: Design Principles and
Practice, Second Edition. New York, US, McGraw - Hill Education. 1344
p. Disponible en https://books.google.com.br/books?id=XypnwwEACAAJ
80
598 p. (42)
Ruiz, R. 2007. El método científico y sus etapas. Mexico DF, MX, Biblioteca
Lascasas. v. 3, 79 p. (3)
Singh, VP. 2017. Handbook of Applied Hydrology, Second Edition. 2nd edition.
ed. New York, US, McGraw - Hill Education. 1440 p. Disponible en
https://www.accessengineeringlibrary.com/content/book/9780071835091
Singh, VP. 2018. Entropy theory in hydrologic science and engineering. Texas,
US, McGraw - Hill Professional. 500 p.
Srinivasan, R; Ramanarayanan, TS; Arnold, JG; Bednarz, ST. 1998. Large area
hydrologic modeling and assessment part II: model application. American
Water Resources Association 34(1):91 - 101. Disponible en
https://onlinelibrary.wiley.com/doi/abs/10.1111/j.1752-
1688.1998.tb05962.x
Uitto, JI; Shaw, R. 2016. Sustainable development and disaster risk reduction.
Tokyo, JP, Springer. 287 p.
82
83
ANEXOS
Variables dependientes
Instrume
Variable Definición Conceptual Definición Operacional Dimensiones Indicadores Escala
nto
Los métodos convencionales, se operacionaliza
Según Singh (2017), los métodos mediante sus dimensiones: Distribución D1: Distribución I1: Parámetros de
convencionales son usados para estimar las Gumbel; Log Normal y Gamma (Pearson Tipo Gumbel forma
Y1
precipitaciones, así como los caudales, III), que vienen hacer los métodos de cálculos; a D2: Normal I2: Parámetros de Registro
Métodos Nominal
convencionales
también denominado volumen de agua, su vez cada una de las dimensiones, se D3:Gamma posición histórico
también llamado cantidad de agua, que se subdividen en indicadores, como parámetros de (Pearson Tipo I3: Parámetros de
mide en una unidad de tiempo. forma, parámetros de posición y parámetros de III) escala
escala
Variables independientes
Instrume
Variable Definición Conceptual Definición Operacional Dimensiones Indicadores Escala
nto
Para Chow et al. (1994) y Villón (2004), el Los Caudales máximos de diseño, se D1: I1: Área Razón
estudio de la cuenca está delimitado por el operacionaliza mediante su dimensión de Características I2: Pendiente Nominal
curso principal, y los puntos más altos geomorfológica I3: Densidad de
características geomorfológicas, que vienen
denominado la divisoria de agua. Las
hacer los métodos de cálculos; a su vez se s drenaje
características geomorfológicas se analiza el Razón
área, longitud de la cuenca y su perímetro, subdividen en indicadores, como área,
pendiente y densidad de drenaje. Ficha de
X1 pendiente, entre otros parámetros.
recopilaci
Caudales Para Chow et al. (1994) y Villón (2004), el total Los Caudales máximos de diseño, se
máximos de de lluvia en 24 horas que se registra en un operacionaliza mediante su dimensión de I1: Alta < 600 mm/año ón de
D2:
diseño evento, es registrado por las estaciones precipitación de 24h, que vienen hacer los I2: Media 601 - 799 informació
Precipitación de n Nominal
meteorológicas, se calculan las máximas
métodos de cálculos; a su vez se subdividen en mm/año
mensuales y de la serie histórica de máximas 24h
los indicadores, como alto, medio y bajo. I3: Baja > 800 mm/año
mensuales.
Los Caudales máximos de diseño, se D3: I1: Tiempo de retorno
Para Chow et al. (1994) y Villón (2004), el
Razón
período de retorno (T) es un estimativo de la
probabilidad de ocurrencia de un evento operacionaliza mediante su dimensión de Periodo de 50 Intervalo
determinado en un periodo determinado, tiempo de retorno, que vienen hacer los métodos retorno I2: Tiempo de retorno
también es la cantidad de tiempo para la cual de cálculos; a su vez se subdividen en 100
la probabilidad de ocurrencia se distribuye indicadores, como Tiempo de retorno de 50, 100 I3: Tiempo de retorno
uniformemente en los periodos que componen
y 200 años. 200
dicha cantidad de tiempo. También llamado
período de recurrencia, es considerado un
concepto estadístico.
Anexo 2. Matriz de Consistencia
TITULO: “Determinación de caudales máximos de diseño mediante métodos convencionales probabilísticos en la cuenca Ilave, Puno – 2022”.
AUTOR: Br. Paredes Quispe, Juan Reynaldo.
Problema Objetivo Hipótesis Variable Dimensiones Indicadores Metodología
5
Anexo 4: Instrumento de recolección de datos
6
7
Anexo 5: Certificado de validación del instrumento recolección de datos
8
9
Anexo 6: NORMAS Y LEYES DEL RECURSO HÍDRICO
10
Anexo 7: NORMAS Y LEYES DEL RECURSO HÍDRICO
11
Anexo 7: NORMAS Y LEYES DEL RECURSO HÍDRICO
12
Anexo 8: MAPAS DE UBICACIÓN DE LA CUENCA DEL RÍO ILAVE
15
Anexo 9: LUGARES DE ORIGEN DE LA CUENCA DEL RÍO ILAVE
16
Anexo 9: LUGARES DE ORIGEN DE LA CUENCA DEL RÍO ILAVE
17
Anexo 10: EL RÍO ILAVE PASANDO POR LA CIUDAD DE ILAVE
Imagen 1. Mirador de Ilave, se precia el puente internacional de Ilave en el cruce del río
Imagen 2. Mirador de Ilave, se aprecia el cauce del río Ilave hacia el lago Titicaca
Imagen 3. Mirador de Ilave, se aprecia el cauce del río Ilave hacia el oeste
18
RESULTADOS VERIFICADOS Y APROBADOS POR UN EXPERTO
19
20
21
22
23
24
25
26
27
FACULTAD DE INGENIERÍA Y ARQUITECTURA
ESCUELA PROFESIONAL DE INGENIERÍA CIVIL
He revisado dicho reporte y concluyo que cada una de las coincidencias detectadas no
constituyen plagio. A mi leal saber y entender la Tesis Completa cumple con todas las
normas para el uso de citas y referencias establecidas por la Universidad César Vallejo.