Milner Elissbam Tesis Bachiller 2017 Part.1 PDF
Milner Elissbam Tesis Bachiller 2017 Part.1 PDF
Milner Elissbam Tesis Bachiller 2017 Part.1 PDF
TESIS:
ALTERACIÓN DE LOS CAUDALES DEL RÍO ASMAYACU POR
LA PRESENCIA DE LOS FENÓMENOS DEL NIÑO Y DE LA NIÑA
EN LOS ULTIMOS 30 AÑOS PARA LA DETERMINACIÓN DE
LAS PROBABLES ZONAS DE INUNDACIÓN EN LA CUENCA
ASMAYACU-CURAHUASI
Presentado por:
Asesor:
CUSCO – PERÚ
2017
DEDICATORIA
A mi madre Florencia.
Por haberme apoyado en todo momento, por creer en mí, por todos sus
consejos, sus valores, por la motivación constante que me ha permitido ser una
persona de bien.
A mi padre Marcelino.
A familiares y amigos
i
DEDICATORIA
ii
AGRADECIMIENTOS
Queremos agradecer a Dios en primer lugar por habernos permitido terminar una
etapa más en nuestras vidas como es la vida universitaria. A nuestros Padres
por estar siempre apoyándonos en los buenos y malos momentos.
iii
RESUMEN
iv
ABSTRACT
v
INTRODUCCIÓN
vi
ÍNDICE GENERAL
DEDICATORIA..................................................................................................... i
RESUMEN ......................................................................................................... iv
ABSTRACT ......................................................................................................... v
INTRODUCCIÓN ............................................................................................... vi
vii
2.2.1 FENÓMENOS NATURALES..................................................... 13
2.2.1.1 FENÓMENO NATURAL ........................................................ 13
2.2.1.2 FENÓMENO EL NIÑO ........................................................... 14
El Niño en el Perú: registro histórico................................. 15
2.2.1.3 Fenómeno la niña .................................................................. 18
Características del fenómeno de la Niña .......................... 18
Registro histórico del fenómeno de la niña ....................... 19
2.2.2 EL CICLO HIDROLÓGICO ....................................................... 20
2.2.3 LA CUENCA HIDROLÓGICA .................................................... 22
2.2.4 CARACTERISTICAS DEL SUELO ........................................... 22
2.2.4.1 Geología general ................................................................... 22
Litologia ............................................................................ 22
Geomorfología local.......................................................... 24
Estructura y geodinámica ................................................. 24
2.2.4.2 Geotecnia............................................................................... 25
Geotecnia de la Zona ....................................................... 25
2.2.5 CARACTERISTICAS FISIOGRAFICAS DE LA CUENCA......... 25
2.2.5.1 Superficie de la cuenca .......................................................... 25
2.2.5.2 Perímetro de la cuenca .......................................................... 26
2.2.5.3 Curva hipsométrica ................................................................ 26
2.2.5.4 Polígono de frecuencia de altitudes ....................................... 26
Altitud media ..................................................................... 27
Altitud más frecuente ........................................................ 27
2.2.6 ÍNDICES REPRESENTATIVOS ................................................ 27
2.2.6.1 Índice o factor de forma de una cuenca (kf) ........................... 27
2.2.6.2 Índice de compacidad o índice de gravelius .......................... 28
2.2.6.3 Número de orden ................................................................... 29
2.2.6.4 Rectángulo equivalente ......................................................... 29
2.2.6.5 Pendiente de la cuenca (s%): ................................................ 30
2.2.6.6 Pendiente del cauce principal ................................................ 31
Pendiente uniforme........................................................... 31
2.2.6.7 Precipitación .......................................................................... 32
Medición de la precipitación ............................................. 32
Cálculo de la precipitación media sobre una zona ........... 33
Polígonos de Thiessen ..................................................... 33
viii
2.2.6.8 Estudio de una tormenta ........................................................ 35
Intensidad ......................................................................... 35
Duración ........................................................................... 35
Frecuencia ........................................................................ 36
2.2.7 ANÁLISIS DE CONSISTENCIA DE LA INFORMACIÓN........... 36
2.2.7.1 Análisis visual grafico ............................................................. 36
2.2.7.2 Análisis estadístico ................................................................ 37
Análisis de Saltos ............................................................. 38
Análisis de Tendencias ..................................................... 39
2.2.8 COMPLETACIÓN DE DATOS HIDROMETEOROLÓGICOS ... 42
2.2.8.1 Completación de datos .......................................................... 42
2.2.9 MODELO HIDROLÓGICO ........................................................ 43
2.2.10 GENERACIÓN DE CAUDALES POR EL MÉTODO LUTZ
SCHOLZ 43
2.2.10.1 Balance hídrico ..................................................................... 43
2.2.10.2 Caudal medio mensual.......................................................... 44
2.2.10.3 Escorrentía ............................................................................ 44
Escorrentía superficial .................................................... 44
Escorrentía subterránea ................................................. 44
Coeficiente de escorrentía .............................................. 44
2.2.11 COEFICIENTE DE AGOTAMIENTO ......................................... 45
2.2.11.1 Agotamiento .......................................................................... 45
2.2.11.2 definición de coeficiente de agotamiento .............................. 45
2.2.11.3 Retención de la cuenca ......................................................... 46
2.2.11.4 Precipitación efectiva ............................................................ 46
2.2.11.5 Déficit de escurrimiento ......................................................... 46
2.2.11.6 Evapotranspiración................................................................ 47
2.2.11.7 Evapotranspiración potencial ................................................ 47
2.2.11.8 Almacenamiento ................................................................... 48
2.2.12 MODELO HIDROLÓGICO DE GENERACIÓN DE CAUDALES
MEDIOS MENSUALES DE LUTZ SCHOLZ ............................................. 48
2.2.12.1 Parámetros hidrológicos para generación de caudales......... 48
Ecuación de balance hídrico ........................................... 48
Coeficiente de escurrimiento .......................................... 49
Precipitación efectiva para el modelo Scholz ................. 52
Retención en la cuenca .................................................. 55
ix
Relación entre descarga y retención .............................. 55
Gasto de la retención...................................................... 56
Coeficiente de agotamiento ............................................ 57
Almacenamiento hídrico ................................................. 59
Abastecimiento de la retención ....................................... 59
2.2.13 CAUDALES MÁXIMOS ............................................................. 60
2.2.13.1 Periodo de retorno de una avenida ....................................... 60
2.2.13.2 Riesgo de falla de una estructura .......................................... 61
2.2.14 MODELO HIDROLÓGICO HEC-HMS....................................... 62
2.2.14.1 Parámetros empleados en el hec - hms ................................ 63
Pérdidas iniciales ............................................................ 63
Número de curva ............................................................ 64
Tiempo de retraso (Lag time) ......................................... 65
2.2.15 MÉTODO RACIONAL MODIFICADO ....................................... 65
2.2.16 SISTEMA DE MODELAMIENTO HIDROLOGICO HEC-RAS
VERSION 5.0.3 ......................................................................................... 67
2.2.16.1 Funciones que realiza el Hec Ras ......................................... 67
2.2.16.2 Ventajas y limitaciones del modelo hec-ras v5.0.3 ................ 68
2.2.16.3 Datos necesarios para el hec-ras .......................................... 70
2.2.17 INUNDACIONES ....................................................................... 72
2.2.18 ESTACIÓN METEOROLÓGICA CURAHUASI ......................... 73
2.2.18.1 Descripción y estado de la estación meteorológica Curahuasi
73
3 CAPÍTULO III: METODOLOGÍA ................................................................ 75
x
3.3.2.1 Descripción de la muestra...................................................... 78
3.3.2.2 Cuantificación de la muestra .................................................. 79
3.3.2.3 Método de muestreo .............................................................. 79
3.3.2.4 Criterios de evaluación de muestra ........................................ 79
3.3.3 CRITERIOS DE INCLUSIÓN .................................................... 79
3.4 INSTRUMENTOS .............................................................................. 80
3.4.1 INSTRUMENTOS METODOLÓGICOS O INSTRUMENTOS DE
RECOLECCIÓN DE DATOS .................................................................... 80
3.4.1.1 Formato de precipitación diaria .............................................. 80
3.4.1.2 Precipitación mensual ............................................................ 81
3.4.2 INSTRUMENTOS DE INGENIERIA .......................................... 83
3.5 PROCEDIMIENTOS DE RECOLECCIÓN DE DATOS ...................... 83
3.5.1 PRECIPITACIÓN ...................................................................... 83
3.6 PROCEDIMIENTOS DE ANÁLISIS DE DATOS ................................ 86
3.6.1 UBICACIÓN .............................................................................. 86
3.6.2 PROCESAMIENTO DE DATOS DE PRECIPITACION ............. 87
3.6.3 DETERMINACIÓN DE LAS CARACTERÍSTICAS
FISIOGRÁFICAS DE LA CUENCA ASMAYACU ...................................... 87
3.6.3.1 Delimitación de la cuenca ...................................................... 87
3.6.3.2 Cálculo del área de la cuenca ................................................ 89
3.6.3.3 Cálculo del perímetro de la cuenca ........................................ 90
3.6.3.4 Cálculo de la curva hipsometrica ........................................... 90
3.6.3.5 Cálculo de la cota media de la cuenca .................................. 93
3.6.3.6 Forma de la cuenca ............................................................... 93
Cálculo del Índice o factor de forma (𝑲𝒇) ......................... 93
Cálculo del Índice de compacidad o índice de Gravelius (𝑲𝒄)
94
Cálculo del Rectángulo Equivalente ................................. 95
3.6.3.7 Cálculo de la pendiente de la cuenca .................................... 97
3.6.3.8 Pendiente longitudinal del cauce ........................................... 98
3.6.3.9 Número de orden ................................................................... 99
3.6.4 PLUVIOMETRÍA ..................................................................... 101
3.6.4.1 Red de estaciones meteorológicas ...................................... 101
3.6.4.2 Influencia de las estaciones meteorológicas ........................ 102
Polígonos de Thiessen para la influencia de las estaciones
sobre la cuenca Asmayacu .............................................................. 102
xi
3.6.5 ANALISIS DE CONSISTENCIA DE DATOS ........................... 104
3.6.5.1 Análisis visual gráfico ........................................................... 104
3.6.5.2 Análisis estadístico de saltos y tendencias .......................... 107
Análisis de saltos ............................................................ 109
Análisis de tendencias .................................................... 115
3.6.6 COMPLETACIÓN DE DATOS DE PRECIPITACIÓN MENSUAL
CON EASY FIT ....................................................................................... 127
3.6.6.1 Histórico de datos de precipitación de la estación meteorológica
de Curahuasi ....................................................................................... 127
3.6.6.2 Completación del histórico de precipitaciones de la estación de
Curahuasi ............................................................................................ 128
3.6.7 GENERACIÓN DE LOS CAUDALES MENSUALES POR EL
MODELO LUTZ SCHOLZ ....................................................................... 144
3.6.7.1 Parámetros morfometricos de la cuenca ............................. 144
3.6.7.2 Datos básicos de la cuenca de estudio ................................ 145
3.6.7.3 Cálculo de parámetros para la obtención de caudales
mensuales ........................................................................................... 148
Cálculo de la evapotranspiración potencial anual ........... 148
Determinación del coeficiente de escurrimiento “C” y el
déficit de escurrimiento “D ............................................................... 149
Determinación del Coeficiente de agotamiento .............. 150
Cálculo del coeficiente de agotamiento mediante aforos 152
Retención de la cuenca (R) ............................................ 153
Precipitación efectiva ...................................................... 156
Gasto de la retención...................................................... 159
Abastecimiento de la retención ....................................... 161
Cálculo de los caudales mensuales................................ 163
3.6.8 CAUDALES MAXIMOS ........................................................... 169
3.6.8.1 Determinación del caudal máximo por el método racional
modificado ........................................................................................... 169
3.6.9 CÁLCULO DE CAUDALES MÁXIMOS CON EL MODELO HEC-
HMS 171
3.6.10 DETERMINACIÓN DE LAS PROBABLES ZONAS DE
INUNDACIÓN EN LA CUENCA DE ASMAYACU ................................... 193
3.6.10.1 Proceso para la determinación de las zonas de inundación con
Hec-Ras 194
Creación del proyecto ................................................... 194
Importación de secciones del Autocad civil 3D ............. 194
xii
Introducción de las secciones del río ............................ 195
Longitud entre sección y sección .................................. 202
Introducción del coeficiente de Manning ....................... 202
Introducción de los caudales para los diferentes periodos
de retorno 204
Condiciones de borde ................................................... 204
Ejecución del programa ................................................ 205
Resultados del proceso ................................................ 206
Probables zonas de inundación .................................. 209
Secciones que presentan desbordes.......................... 209
Perfil general de velocidades de los caudales ............ 212
4 CAPITULO IV: RESULTADOS ................................................................ 213
xiii
ÍNDICE DE TABLAS
Tabla 1 Cuadro de operacionalización de variables ........................................... 7
Tabla 2 Recurrencia del fenómeno el Niño por años de presentación ............. 16
Tabla 3 Resumen cronológico del fenómeno el niño según su intensidad 1847-
1997-19998 ...................................................................................................... 16
Tabla 4 Cronología de ocurrencia de los fenómenos del Niño en el Perú ....... 17
Tabla 5 Fenómenos del Niño ocurridos en los últimos 30 años ....................... 17
Tabla 6 Histórico de fenómeno de la Niña 1950-2016 ..................................... 20
Tabla 7 Numero de ocurrencias de el fenomeno la niña en los ultimos 30 años
......................................................................................................................... 20
Tabla 8 Ecuaciones del coeficiente de escurrimiento y su coeficiente de
correlación........................................................................................................ 51
Tabla 9 Límite superior para la precipitación efectiva ...................................... 53
Tabla 10 Coeficientes para el cálculo de la precipitación efectiva ................... 53
Tabla 11 Lámina de agua acumulada en los 3 tipos de almacén hídrico ......... 55
Tabla 12 Valores recomendados de riesgo admisible de obras de drenaje ..... 62
Tabla 13 VALORES DEL NUMERO DE CURVA ............................................. 64
Tabla 14 Formato de precipitación total diario ................................................. 80
Tabla 15 Formato de registro de precipitación mensual .................................. 82
Tabla 16 Registro de precipitación diaria (ejemplo: año 1985) ........................ 84
Tabla 17 Registro de precipitación mensual .................................................... 85
Tabla 18 Cálculos curva hipsometrica ............................................................. 92
Tabla 19 Cálculo de los valores li del rectángulo equivalente .......................... 96
Tabla 20 Longitudes de las franjas de las curvas de nivel de la cuenca según su
altitud ............................................................................................................... 98
Tabla 21 Precipitación total mensual (mm) - Estación Curahuasi ................. 105
Tabla 22 Precipitación mensual (mm) – Estación Curahuasi: Serie 1 y 2 ...... 108
Tabla 23 Tabla de distribución “t – Student” .................................................. 111
Tabla 24 Valores para la distribución F de Fischer (a = 0.05) ........................ 114
Tabla 25 Valores registrables entre la precipitación mensual (mm) Vs. el tiempo
(1er bloque) .................................................................................................... 116
Tabla 26 Valores registrables entre la precipitación mensual (mm) Vs. el tiempo
(2do bloque) ................................................................................................... 117
Tabla 27 Desviaciones estándar de cada período ......................................... 122
xiv
Tabla 28 Cálculo del producto t*Sp ................................................................ 123
Tabla 29 Datos de precipitación mensual de la estación Curahuasi .............. 129
Tabla 30 Resultados de la prueba de bondad de ajuste (enero) ................... 131
Tabla 31 Resultados de la prueba de bondad de ajuste (enero) ................... 133
Tabla 32 Precipitación mensual completada en mm (enero) ......................... 137
Tabla 33 Histórico de datos de precipitación mensual de la estación de Curahuasi
completado..................................................................................................... 143
Tabla 34 Parámetros morfometricos de la cuenca Asmayacu ....................... 144
Tabla 35 Coordenadas del centroide de la cuenca Asmayacu ...................... 145
Tabla 36 Cálculo de la radiación extraterrestre de la cuenca ........................ 146
Tabla 37 Cuadro resumen: datos básicos de la cuenca ................................ 147
Tabla 38 Formulas para la obtencion del coeficiente de escurrimiento .......... 149
Tabla 39 Aforos realizados al río Asmayacu .................................................. 152
Tabla 40 Cálculo de coeficiente de agotamiento por medio de aforos ........... 152
Tabla 41 Lámina de agua según la pendiente del suelo ................................ 153
Tabla 42 Áreas de los lagos y lagunas .......................................................... 155
Tabla 43 Cálculo de la retención de la cuenca Asmayacu ............................. 156
Tabla 44 Cálculo de las precipitaciones efectivas mensuales para el año 1985
....................................................................................................................... 159
Tabla 45 Valores de "i" para la temporada seca ............................................ 160
Tabla 46 Cálculo del gasto mensual (1985) ................................................... 161
Tabla 47 Coeficientes de abastecimiento propuesto por Scholz para la región
cusco .............................................................................................................. 161
Tabla 48 Coeficientes de abastecimiento propuesto para esta investigación 162
Tabla 49 Calculo del abastecimiento mensual (1985) .................................... 162
Tabla 50 Calculo de los caudales mensuales por el modelo hidrologico Lutz
Scholz ............................................................................................................ 164
Tabla 51 Generación de caudales mensuales con el modelo Lutz Scholz desde
1985 hasta el 2016......................................................................................... 166
Tabla 52 Cálculo del periodo de retorno ........................................................ 176
Tabla 53 Resultados de las pruebas de bondad de ajuste ............................ 178
Tabla 54 Precipitaciones máximas en 24 hrs según el periodo de retorno .... 179
Tabla 55 Patrón de tormenta.......................................................................... 181
Tabla 56 Resultados de los caudales máximos con el hec-hms .................... 193
xv
Tabla 57 Valores del coeficiente de Manning ................................................ 203
Tabla 58 Variación de las velocidades del flujo del río ................................... 212
Tabla 59 Caudales mensuales (Lutz Scholz) ................................................. 213
Tabla 60 Caudal máximo (método racional modificado) ................................ 214
Tabla 61 Resultados con el modelo HEC-HMS ............................................. 214
Tabla 62 Incremento de los caudales con fenómeno de la Niña .................... 224
xvi
ÍNDICE DE ILUSTRACIONES
Ilustración 1 Cuenca de Asmayacu-Curahuasi vista desde el puente................ 1
Ilustración 2 Calentamiento de las aguas del océano Pacifico ........................ 15
Ilustración 3 Variacion de la temperatura durante el fenómeno de la Niña ...... 18
Ilustración 4 Registro histórico de fenómenos de la Niña ................................ 19
Ilustración 5 Ciclo hidrológico........................................................................... 21
Ilustración 6 Curva hipsometrica ...................................................................... 26
Ilustración 7 Polígono de frecuencia de altitudes (ejemplo) ............................. 27
Ilustración 8 Representación del índice de compacidad o de grávelius ........... 28
Ilustración 9 Número de orden ......................................................................... 29
Ilustración 10 Rectángulo equivalente ............................................................. 30
Ilustración 11 Longitud del cauce principal ...................................................... 32
Ilustración 12 Polígonos de Thiessen .............................................................. 34
Ilustración 13 Seria histórica del río Bebedero ................................................. 37
Ilustración 16 Almacenamiento hídrico durante época de lluvias (coeficientes ai
%)..................................................................................................................... 60
Ilustración 17 Valores de periodo de retorno.................................................... 61
Ilustración 18 Flujograma del diseño de ingeniería .......................................... 77
Ilustración 19 Mapa de ubicación de la cuenca ............................................... 86
Ilustración 20 Delimitación de la cuenca .......................................................... 87
Ilustración 21 Delimitación de la cuenca asmayacu ......................................... 88
Ilustración 22 Calculo de áreas y perímetro en Autocad .................................. 89
Ilustración 23 Áreas según sus altitudes .......................................................... 91
Ilustración 24 Curva hipsometrica y frecuencia de altitudes ............................ 92
Ilustración 25 Calculo del factor de forma ........................................................ 94
Ilustración 26 Rectángulo equivalente ............................................................. 96
Ilustración 27 Longitud de franjas .................................................................... 97
Ilustración 28 Número de orden de la cuenca ................................................ 100
Ilustración 29 Mapa de estaciones ................................................................. 101
Ilustración 30 Cuenca Asmayacu y las estaciones meteorológicas más cercanas
....................................................................................................................... 102
Ilustración 31 Polígonos de Thiessen de las 4 estaciones adyacentes a la cuenca
....................................................................................................................... 103
Ilustración 32 Histograma Estación Curahuasi............................................... 106
xvii
Ilustración 33 Histograma con línea de tendencia en la media – Estación
Curahuasi ....................................................................................................... 120
Ilustración 34 Histograma con línea de tendencia en la Des. Estándar – Estación
Curahuasi ....................................................................................................... 126
Ilustración 35 Precipitación en (mm) enero (Easyfit) ...................................... 130
Ilustración 36 Precipitación en (mm) enero (Easyfit) ...................................... 132
Ilustración 37 Histograma de distribución log Pearson 3 (enero) ................... 134
Ilustración 38 Histograma de distribución normal (enero) .............................. 134
Ilustración 39 Histograma de distribución log normal (enero) ........................ 135
Ilustración 40 Generación de números aleatorios log Pearson 3 (enero) ...... 135
Ilustración 41 Números aleatorios generados (enero) ................................... 136
Ilustración 42 Histograma de distribución log normal (febrero) ...................... 138
Ilustración 43 Histograma de distribución log Pearson (febrero) .................... 138
Ilustración 44 Histograma de distribución normal (febrero) ............................ 139
Ilustración 45 Generación de números aleatorios log normal por Easyfit (febrero)
....................................................................................................................... 139
Ilustración 46 Números aleatorios generados (febrero) ................................. 140
Ilustración 47 Precipitación mensual completada en mm (febrero) ................ 141
Ilustración 48 Conversión de coordenadas UTM a grados minutos y segundos
....................................................................................................................... 145
Ilustración 49 Pendientes de la cuenca.......................................................... 154
Ilustración 50 Lagos y lagunas de la cuenca asmayacu ................................ 155
Ilustración 51 Precipitaciones efectivas (1985) .............................................. 159
Ilustración 52 Variación del caudal (1985) ..................................................... 165
Ilustración 53 Variación de los caudales mensuales segun el modelo Lutz Scholz
en los últimos 30 años ................................................................................... 167
Ilustración 54 Sequia severa con fenómeno del niño ..................................... 167
Ilustración 55 Incremento notable de caudales con fenómeno del Niño ........ 168
Ilustración 56 Máximos caudales sin fenómenos naturales ........................... 168
Ilustración 57 Creación del proyecto .............................................................. 172
Ilustración 58 Sistema de drenaje embalse y unión ....................................... 172
Ilustración 59 Elementos de base, pérdida, transformación, base flujo ......... 173
Ilustración 60 Mapa de numero de curva ...................................................... 174
Ilustración 61 Área de lagunas ....................................................................... 175
xviii
Ilustración 62 Máximas precipitaciones .......................................................... 177
Ilustración 63 Histograma distribución lognormal 3p ...................................... 178
Ilustración 64 Cuadro de probabilidades y precipitaciones máximas ............. 179
Ilustración 65 Modelo meteorologico .............................................................. 180
Ilustración 66 Control de especificaciones ..................................................... 181
Ilustración 67 Datos de la serie de tiempo ..................................................... 186
Ilustración 68 Simulación de los caudales máximo ........................................ 187
Ilustración 69 Pico de descarga del caudal .................................................... 188
Ilustración 70 Caudal y altura de precipitación ............................................... 189
Ilustración 71 Resultados de la descarga para un periodo de retorno=50 años
....................................................................................................................... 190
Ilustración 72 Caudales en función a los intervalos de tiempo de 6 min ........ 191
Ilustración 73 Caudal máximo para tiempo de retorno 100 años ................... 192
Ilustración 74 Caudal máximo para tiempo de retorno 348 años ................... 192
Ilustración 75 Caudal máximo para tiempo de retorno 500 años ................... 193
Ilustración 76 Creación y archivo del proyecto ............................................... 194
Ilustración 77 Importación de secciones ........................................................ 195
Ilustración 78 Seccion 1 (0+00) ...................................................................... 195
Ilustración 79 Sección 02 (0+50) .................................................................... 196
Ilustración 80 Sección 03 (0+100) .................................................................. 196
Ilustración 81 Sección 04 (0+150) .................................................................. 196
Ilustración 82 Sección 05 (0+200) .................................................................. 197
Ilustración 83 Sección 06 (0+250) .................................................................. 197
Ilustración 84 Sección 07 (0+300) .................................................................. 197
Ilustración 85 Sección 08 (0+350) .................................................................. 198
Ilustración 86 Sección 09 (0+400) .................................................................. 198
Ilustración 87 Sección 10 (0+450) .................................................................. 198
Ilustración 88 Sección 11 (0+500) .................................................................. 199
Ilustración 89 Sección 12 (0+550) .................................................................. 199
Ilustración 90 Sección 13 (0+600) .................................................................. 199
Ilustración 91 Sección 14 (0+650) .................................................................. 200
Ilustración 92 Sección 15 (0+700) .................................................................. 200
Ilustración 93 Sección 16 (0+750) .................................................................. 200
Ilustración 94 Sección 17 (0+800) .................................................................. 201
xix
Ilustración 95 Sección 18 (0+850) .................................................................. 201
Ilustración 96 Sección 19 (0+900) .................................................................. 201
Ilustración 97 Longitudes entre secciones del rio ........................................... 202
Ilustración 98 Coeficiente de Manning para cada sección del río .................. 203
Ilustración 99 Caudales máximos para diferentes periodos de retorno .......... 204
Ilustración 100 Condiciones de borde ............................................................ 205
Ilustración 101 Flujo subcritico ....................................................................... 205
Ilustración 102 Compilando Hec Ras ............................................................ 206
Ilustración 103 Modelamiento del rio en la zona llana de la cuenca .............. 209
Ilustración 104 Sección Critica 04 (0+150)..................................................... 210
Ilustración 105 Sección Critica 05 (0+200)..................................................... 210
Ilustración 106 Sección Critica 06 (0+250)..................................................... 211
Ilustración 107 Sección Critica 17 (0+800)..................................................... 211
Ilustración 108 Perfil longitudinal del río ......................................................... 212
Ilustración 109 Variación de los caudales mensuales según el modelo Lutz
Scholz en los últimos 30 años ........................................................................ 214
Ilustración 110 Sección Critica 04 (0+150)..................................................... 215
Ilustración 111 Sección Critica 05 (0+200)..................................................... 216
Ilustración 112 Zona de desborde identificada (0+200) ................................ 216
Ilustración 113 Sección Critica 06 (0+250)..................................................... 217
Ilustración 114 Sección Critica 17 (0+800)..................................................... 217
Ilustración 115 Zona de inundacion identificada (0+800) ............................... 218
xx
1 CAPÍTULO I: PLANTEAMIENTO DEL PROBLEMA
Fuente: PROPIA
1
1.1.2 FORMULACIÓN INTERROGATIVA DEL PROBLEMA
¿En qué medida se alteran los caudales del río Asmayacu por la presencia de
los fenómenos del Niño y de la Niña en los últimos 30 años y cuáles serán las
probables zonas de inundación en la cuenca de Asmayacu – Curahuasi?
a) ¿De qué manera se alteran los caudales por la presencia del fenómeno
de la Niña en la cuenca?
b) ¿De qué manera se alteran los caudales por la presencia del fenómeno
del Niño en la cuenca?
c) ¿Cuál es el caudal máximo generado por las precipitaciones durante los
últimos 30 años?
d) ¿Cuáles serán las zonas de probable inundación en la cuenca del río
Asmayacu por efectos del Fenómeno de la Niña?
2
Se dispone el Manual de Hidrología, Hidráulica y Drenaje, también el estudio
Definitivo Asmayacu – Lucmos de Plan Meris, entre otros textos que sirven como
guía en esta investigación.
3
La zona en estudio se limita a la cuenca del río Asmayacu perteneciente
al distrito de Curahuasi, provincia de Abancay, departamento de
Apurímac.
Se limita al periodo de estudio de 30 años contados desde el año 1985.
El cálculo de los caudales mensuales se limita al uso del modelo
hidrológico de Lutz Scholz
El cálculo de los caudales máximos se limita al uso del método racional
modificado y el modelo HEC-HMS
El cálculo de las zonas de inundación se limita al uso del software HEC
RAS
4
1.5 HIPÓTESIS
El fenómeno del Niño genera cambios climáticos lo cual altera el régimen de los
caudales produciendo sequias en el 50% de ocurrencias de estos fenómenos
naturales, de igual manera el fenómeno de la Niña también genera variaciones
climáticas y produce incrementos en los caudales en un 50% de ocurrencia de
estos fenómenos naturales y que estos a su vez generan desbordes en la zona
del poblado de Pisonaypata.
Indicador. - mm/h
5
PRECIPITACIÓN EN PRESENCIA DEL FENÓMENO EL NIÑO
Indicador. - mm/hr.
ÁREA DE LA CUENCA.
Indicador. - km2
PENDIENTE DE LA CUENCA
Indicador. - % (porcentaje)
CAUDAL.-
Indicador. - m3/s
6
1.6.3 OPERACIONALIZACIÓN DE VARIABLES
VARIABLES DEPENDIENTES
VARIABLE DESCRIPCION DE VARIABLE NIVEL INDICADORES
Alto
Es la cantidad de agua, medido
Caudal en volumen, que se mueve en una Medio m3/s
unidad de tiempo.
Bajo
VARIABLES INDEPENDIENTES
Se denomina pendiente a la
pendiente de la inclinación natural de la superficie
Altura % (porcentaje)
Cuenca que comprende una cuenca
respecto a la horizontal.
Fuente propia
7
2 CAPÍTULO II: MARCO TEÓRICO
LUGAR: Cusco-Perú
AÑO: 2015
RESUMEN
CONCLUSIONES:
8
Entre los caudales determinados por el HEC-HMS y los caudales
registrados por la estación hidrométrica La Angostura, en la
“cuenca alta del río Apurímac hasta la estación La Angostura”, no
existe diferencia significativa en época de estiaje mientras que en
época de lluvias existe diferencia significativa.
Entre los caudales determinados con el “Modelo Lutz Scholz” y los
caudales obtenidos con el HEC HMS, en la “cuenca alta del río
Apurímac hasta la estación La Angostura”, existe diferencia
significativa en época de estiaje como en época de lluvias.
RESUMEN
9
general de la atmósfera HadCM3. Para llevar a cabo el escalamiento de
las variables meteorológicas, se utilizó el programa SDSM 4.2, mediante
el cual se generaron 6 escenarios de cambio climático, con horizontes de
30 años (2005-2035) que fueron evaluados en el Modelo Integrado de la
cuenca. En cuanto a los resultados obtenidos para los escenarios de
cambio climático, para la temperatura en la estación Rengo se verificaron
aumentos de los valores mensuales, estacionales y anuales promedio,
mientras que para la estación Convento Viejo se verificaron disminuciones
con respecto al período observado (1980-2005); para la precipitación, en
las estaciones estudiadas se observó una tendencia al aumento en los
valores mensuales, estacionales y anuales promedio con respecto al
período observado (1971-2005). De acuerdo a esto, los caudales medios
mensuales de salida de la cuenca del río Cachapoal, se verían afectados
por el cambio climático, ya que con respecto al Período Base de
calibración (1999-2005), se obtuvieron aumentos de los valores promedio
en el período de deshielo y disminuciones en el período invernal. Para los
recursos hídricos subterráneos, no se pudo concluir que existirían efectos
asociados al cambio climático.
CONCLUSIONES:
10
TITULO: LOS FENÓMENOS DE EL NIÑO Y LA NIÑA
RESUMEN
11
CONCLUSIONES:
AÑO: 2011
RESUMEN
12
hipsométrica, la pendiente media de la cuenca por su valor de 56,4%,
indica que es una cuenca con alta pendiente.
CONCLUSIONES
13
Sin embargo, la formación de una gota de lluvia es un fenómeno natural de la
misma manera que un huracán. Esta expresión también se refiere, en general, a
los peligrosos fenómenos naturales también llamados "desastres naturales". La
lluvia, por ejemplo, no es en sí un "desastre", pero puede ser así dependiendo
de la perspectiva humana, si ciertas condiciones se reúnen. La mala planificación
urbana, con la construcción de estructuras en lugares vulnerables a
inundaciones u otras personas puede causar efectos desastrosos para los seres
humanos (SENAMHI, s.f.)
14
Ilustración 2 Calentamiento de las aguas del océano Pacifico
Fuente: SENAMHI
En los últimos cinco siglos ha habido por lo menos 120 episodios El Niño,
según las investigaciones históricas recopiladas por Quinn W., Neal V., y
Antúnez de Mayolo S. A diferencia de los países del Hemisferio Norte, en el
Perú la instrumentación adecuada para la observación del clima se inicia
recién en la segunda década del Siglo XX y de manera puntual en algunas
regiones del país. Es a partir del año 1965, cuando la mayoría de las
estaciones meteorológicas e hidrológicas son instaladas en nuestro territorio,
que se inicia la observación climática de manera sostenida. (SENAMHI, s.f.)
El Niño como tal, es conocido desde que civilizaciones pre-incas como los
Moche, los Lima y los Nazca se asentaron en las costas del Perú antiguo. La
geomorfología, los estudios de sedimentos y la paleontología señalan que el
Fenómeno El Niño ocurre por lo menos desde hace 40,000 años. Además,
estas investigaciones arqueológicas demuestran que cambios drásticos del
clima afectaron la costa central del Perú (SENAMHI, s.f.).
15
y hasta antes de El Niño extraordinario de 1997/98, ocurrieron unos 25
episodios El Niño de diferente intensidad; las referencias bibliográficas indican
que los eventos El Niño de 1891 y 1925, fueron eventos de intensidad
comparable a los de 1982/83 y 1997/98. En lo que va del siglo XXI, de acuerdo
al índice ONI (Oceanic Niño Index) de la NOAA, se han presentado cuatro
episodios El Niño en el Pacífico central; dos de intensidad débil (Años 2004/05
y 2006/07) y dos de intensidad moderada (años 2002/03 y 2009/10).
(SENAMHI, s.f.)
Fuente: CENEPRED
Fuente: CENEPRED
16
Tabla 4 Cronología de ocurrencia de los fenómenos del Niño en el Perú
Fuente: CENEPRED
Fuente: CENEPRED
17
2.2.1.3 Fenómeno la niña
18
Los vientos alisios anormalmente intensos, ejercen un mayor efecto de
arrastre sobre la superficie del océano.
Como resultado de la aparición de aguas relativamente frías a lo largo del
Ecuador, la temperatura superficial del mar disminuye por debajo del valor
medio climatológico.
Durante los eventos de La Niña las aguas calientes en el Pacífico
ecuatorial, se concentran en la región junto a Oceanía y es sobre esta
región, donde se desarrolla la nubosidad y la precipitación más intensa.
"La Niña" comenzó en 1903, y siguió en 1906, 1909, 1916, 1924, 1928, 1938,
1950, 1954, 1964, 1970, 1973, 1975, 1988, y en 1995.Siendo el más intenso
el de 1988/1989.4 (ElClima.com.mx, s.f.)
Fuente: NOAA
19
Tabla 6 Histórico de fenómeno de la Niña 1950-2016
Fenómeno La Niña
Duración
Comienzo Fin
(meses)
mar-50 feb-51 12
jun-54 mar-56 22
may-56 nov-56 7
may-64 ene-65 9
jul-70 ene-72 19
jun-73 jun-74 13
sep-74 abr-76 20
sep-84 jun-85 10
may-88 jun-89 14
sep-95 mar-96 7
jul-98 jun-00 23
dic-00 may-01 5
jul-07 feb-08 8
jul-10 abr-12 22
20
forma de precipitación. Durante su trayecto hacia la superficie de la tierra, el agua
precipitada puede volver a evaporarse o ser interceptada por las plantas o las
construcciones, luego fluye por la superficie hasta las corrientes o se infiltra. El
agua interceptada y una parte de la infiltrada y de la que corre por la superficie
se evapora nuevamente. De la precipitación que llega a las corrientes, una parte
se infiltra y otra llega hasta los océanos y otros grandes cuerpos de agua, como
presas y lagos. Del agua infiltrada, una parte es absorbida por las plantas y
posteriormente es transpirada, casi en su totalidad, hacia la atmósfera y otra
parte fluye bajo la superficie de la tierra hacia las corrientes, el mar u otros
cuerpos de agua, o bien hacia zonas profundas del suelo (percolación) para ser
almacenada como agua subterránea y después aflorar en manantiales, ríos o el
mar.
Por lo tanto, el ciclo hidrológico es un proceso donde el agua por efecto del calor
se evapora, las plantas la transpiran, se infiltra, percola, almacena por el suelo,
afloran a la superficie mediante manantes y transcurren mediante la escorrentía.
21
2.2.3 LA CUENCA HIDROLÓGICA
Litologia
22
fuertemente silicificadas ; afloran en el flaco derecho de Galantin y en
ambos flancos de la quebrada de Asmayacu .
Se halla fuertemente fracturada y cortada por numerosos diques de rocas
intrusivas ( microdioritas, tonalitas , etc) con rumbos N 110° a N 140° e
inclinaciones de 50° a 70° al NE.
Cuaternario Reciente(Qr)
23
Geomorfología local
Altas Mesetas
Zona de relieve suave ubicada en una altitud que varía de 4200 – 4500
msnm.
Valle Colgado
Se caracteriza por ser poco profundo bordeado bordeado por flancos de
pendiente regular a moderada por flancos de pendientes regular a
moderada con relieve ondulado.
Estructura y geodinámica
Estructura
Geodinámica
Los fenómenos geodinámicos externos detectados se desarrollan
mayormente en ladera de pendientes fuertes a muy fuerte. Se ha tipificado
los fenómenos siguientes:
24
Deslizamientos de escombros, deslizamientos antiguos, deslizamientos
de tierras, derrumbes, asentamientos; desprendimientos de roca; tramos
propensos a derrumbes, flujos de lodo y cárcavas activas.
2.2.4.2 Geotecnia
Geotecnia de la Zona
Rocas y Suelos
Consideraciones geotécnicas
La capacidad de carga de el deposito fluvio - Aluvial es superior a
1.2kg/cm2. Por lo tanto se sugiere cimentar a una profundidad de no
menor de de 1.5m en lecho del rio.
Estabilidad de Ladera
Las laderas morfológicamente se ubican a diferentes niveles, en torno al
riachuelo Asmayacu caracterizándose por presentar perfiles
longitudinales con pendientes que varían del 8 al 16% a lo largo del rio
Asmayacu con flujo rápido moderadamente turbulento bordeado por
flancos de pendiente que determinan en general un riesgo geodinámico
mínimo.
25
2.2.5.2 Perímetro de la cuenca
26
Ilustración 7 Polígono de frecuencia de altitudes (ejemplo)
Fuente: www.plusformacion.com
Altitud media
𝑎𝑟𝑒𝑎 𝐴
𝐾𝑓 = = 2
𝑝𝑒𝑟𝑖𝑚𝑒𝑡𝑟𝑜 ∗ 𝑝𝑒𝑟𝑖𝑚𝑒𝑡𝑟𝑜 𝐿
Si una cuenca tiene un factor F mayor que otra, existe mayor probabilidad de
tener una tormenta intensa simultanea sobre toda la extensión de la cuenca. Por
27
el contrario, si la cuenca tiene un F menor, tiene menos tendencia a concentrar
las intensidades de lluvias que una cuenca de igual área pero con un factor F
mayor (Villon Bejar, HIDROLOGIA, 2002).
𝑃
Kc = 0.28
√A
Kc = Coeficiente de compacidad
- Clase Kc2: Rango entre 1.25 y 1.5. Corresponde a forma oval redonda
a oval oblonga.
- Clase Kc3: Rango entre 1.5 y 1.75. Corresponde a forma oval oblonga a
rectangular oblonga. (Gobierno Regional Lambayeque, 2013)
28
2.2.6.3 Número de orden
29
Ilustración 10 Rectángulo equivalente
Donde:
L=Lado mayor del rectángulo en km.
l= Lado menor del rectángulo en km.
A= superficie de la cuenca en km2.
P= Perimetro de la cuenca Km.
30
Criterio de ALVORD
Donde:
S = Pendiente de la cuenca, %
D= desnivel constante entre curvas de nivel, Km
L= Longitud de las curvas de nivel dentro de la cuenca, Km
A= Área de la cuenca, Km2 (Villon Bejar, HIDROLOGIA, 2002)
Pendiente uniforme
Según Villon Bejar (2002), este método considera la pendiente del cauce, como
la relación entre el desnivel que hay entre los extremos del cauce y la proyección
horizontal de su longitud (pag. 55).
31
S = Pendiente
H = Diferencia de cotas entre los extremos del cauce, en km
L = Longitud del cauce, en km
Ilustración 11 Longitud del cauce principal
2.2.6.7 Precipitación
Medición de la precipitación
32
de medición se clasifican de acuerdo con el registro de las precipitaciones en
pluviómetros y pluviografos (Villon Bejar, HIDROLOGIA, 2002).
Polígonos de Thiessen
33
delimitar la zona de influencia de cada estación dentro del conjunto de
estaciones. (pág. 140)
34
2.2.6.8 Estudio de una tormenta
Según Villon Bejar (pág. 86) Se entiende por tormenta al conjunto de lluvias que
obedecen a una misma perturbación meteorológica y de características bien
definidas.
Intensidad
I max = P/T
Donde:
T = Tiempo (hr)
Duración
35
Frecuencia
36
Ilustración 13 Seria histórica del río Bebedero
37
saltos, tanto en la media como en la desviación estándar (Villon Bejar, Hidrologia
Estadistica, 2011).
Análisis de Saltos
Consistencia de la Media
𝑛1 1/2
𝑛2
1 1
𝑥̅2 = ∑ 𝑥𝑗 𝑆2 (𝑥) = ⌊ ∑(𝑥𝑗 − 𝑥̅2 )2 ⌋
𝑛2 𝑛2 − 1
𝑗=1 𝑖=1
Donde:
𝑥̅1 , ̅̅̅
𝑥2 = media de los periodos
n = tamaño de la muestra
38
El análisis estadístico consiste en probar, mediante la prueba F, si los valores de
las desviaciones estándar de las submuestras son estadísticamente iguales o
diferentes, con un 95% de probabilidad o con un 5% de nivel de significación.
se obtiene de las tablas F para una probabilidad del 95%, es decir, con un nivel
de significación = 0.05 y grados de libertad:
Análisis de Tendencias
Tendencia en la Media
𝑇𝑚 = 𝐴𝑚 + 𝐵𝑚 𝑡 + 𝐶𝑚 𝑡 2 + 𝐷𝑚 𝑡 3 + ⋯.
𝑇𝑚 = 𝐴𝑚 + 𝐵𝑚 𝑡
39
b. Evaluación de la tendencia Tm
𝑅√𝑛 − 2
𝑡𝑐 =
√1 − 𝑅 2
donde:
R = coeficiente de correlación
Cálculo de tt
Corrección de la información:
Reemplazando Tm
40
La tendencia en la desviación estándar Ts, se expresa en forma general por la
ecuación polinomial.
𝑇𝑠 = 𝐴𝑠 + 𝐵𝑠 𝑡 + 𝐶𝑠 𝑡 2 + 𝐷𝑠 𝑡 3 + ⋯
𝑇𝑠 = 𝐴𝑠 + 𝐵𝑠 𝑡
1/2
∑𝑛𝑖=1(𝑡𝑖 − 𝑡̅)2
𝑆𝑡 = [ ]
𝑛−1
𝑥(𝑡)
́ − 𝑇𝑚
𝑍𝑡 =
𝑇𝑆
𝑥(𝑡)
́ − 𝑇𝑚
𝑍𝑡 = . 𝑇̅𝑠 + 𝑇
̅̅̅̅
𝑚
𝑇𝑆
41
La serie Z es una serie homogénea y consistente al 95% de probabilidad.
Donde:
: Promedio de datos
42
Para la completación de datos, mediante generación de números aleatorios, se
debe probar por los test de Chi-Cuadrado o Kolmogorov Smirnov si la serie
hidrológica se ajusta a la distribución seleccionada.
- Para una mayor comprensión del rol de los componentes del ciclo hidrológico
en una determinada cuenca.
Para el modelo hidrológico de Lutz Scholz se debe tener en cuenta los siguientes
conceptos:
43
Se entiende por balance hídrico, al equilibrio entre todos los recursos
hídricos que ingresan al sistema y los que salen del mismo, en un intervalo de
tiempo determinado.
Para el balance hídrico se considera tanto las entradas como salidas del sistema
físico, almacenamiento y transporte interno.
2.2.10.3 Escorrentía
Escorrentía superficial
Escorrentía subterránea
Coeficiente de escorrentía
44
en una cuenca a lo largo de periodo de tiempo, e dividido por el volumen total
precipitado Vp.
Esta definición es general y se usa para saber la producción media anual de una
cuenca. A nivel de episodio también se pueden definir estas cantidades. Si se
limita al área de la cuenca A, el volumen precipitado por unidad de tiempo es
exactamente el producto de la intensidad de lluvia i por el área de la cuenca. El
volumen escurrido medio en la cuenca por unidad de tiempo es el caudal medio
medido en la cuenca. Así:
Donde:
2.2.11.1 Agotamiento
45
de decrecimiento exponencial del caudal del manantial que representa la
descarga del sistema, como consecuencia de la disminución de reservas que
sigue a la interrupción de la alimentación.”
46
evapotranspiración real también se conoce como déficit de escurrimiento” (
Universidad Agraria La Molina, s.f.).
2.2.11.6 Evapotranspiración
Chow (2000), indica que “Los mismos factores que dominan la evaporación
desde una superficie de agua abierta también dominan la evapotranspiración,
los cuales son: el suministro de energía y el transporte de vapor. Además, el
suministro de humedad a la superficie de evaporación es un tercer factor que se
debe tener en cuenta. A medida que el suelo se seca, la tasa de
evapotranspiración cae por debajo del nivel que generalmente mantiene en un
suelo bien humedecido”. (p.93)
47
2.2.11.8 Almacenamiento
Según Scholz, (1980), “entre los almacenes naturales que producen el efecto de
la retención en la cuenca, se puede distinguir tres tipos con mayor importancia:
acuíferos, nevados, lagunas y pantanos”. (p.14)
Este modelo hidrológico fue elaborado por el Ingeniero Lutz Scholz, miembro de
la Cooperación Alemana en 1979 – 1980, en el marco de Cooperación Técnica
de la República de Alemania a través del PLAN MERISS. 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 pueden 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. Este modelo fue implementado con fines de
pronosticar caudales a escala mensual, siendo utilizado en estudios de
proyectos de riego.
Dónde:
48
CMi = Caudal mensual (mm/mes)
Dónde:
Q = descarga (m3/s)
Coeficiente de escurrimiento
49
entre el caso de lluvias torrenciales para el cálculo de avenidas y el caso del
caudal anual determinado a partir de la precipitación total anual.
Dónde:
L = Coeficiente de Temperatura
En el sur, en la sierra el método de Turc falla por el clima muy especial, clima
montañoso con temperaturas bajas alrededor de 7c° ubicada en la zona tropical.
Por esta razón se ha buscado un método más apropiado para la región del
proyecto. Se ha probado varias ecuaciones de regresión lineal, cuadrada y
logarítmica entre el déficit D o el coeficiente C dependiente de la precipitación
y/o evaporación y/o temperatura, empleando la relación (d) entre el déficit y el
coeficiente de escurrimiento.
50
Las fórmulas para el cálculo de coeficientes de escurrimiento, para la sierra del
Perú, se muestran a continuación:
Nuestra cuenca está ubicada en la sierra sur del Perú y por ende se utiliza las
ecuaciones siguientes desarrolladas sobre la base de observaciones en la región
de Cusco y Huancavelica propuestas por Scholz.
Donde:
C: Coeficiente de escurrimiento
51
Donde:
TF = Componente de la temperatura
52
Tabla 9 Límite superior para la precipitación efectiva
Dónde:
53
Donde:
C: Coeficiente de escurrimiento
Q: Caudal anual
P: Precipitación Total
Debemos ajustar la suma de la PEi por interpolación entre dos series de valores
PE determinados por los coeficientes de la TABLA N° 5: Coeficientes para el
cálculo de la precipitación efectiva, según curva I y II o curva II y III.
54
Retención en la cuenca
Dónde:
55
Qo: Descarga inicial
a: Coeficiente de agotamiento
t: tiempo
Gasto de la retención
Según Lutz Scholz, el Gasto de la retención considera que los caudales durante
la estación seca están relacionados a los del mes anterior por la relación
siguiente:
Donde:
b0: relación entre la descarga del mes actual y del mes anterior
a: coeficiente de agotamiento
t: tiempo
56
Donde:
R: Retención de la cuenca.
Teniendo en cuenta que para el cálculo de para el mes de abril (i=1), mayo
0 (i=2) así hasta octubre (i=7).
Entonces:
Coeficiente de agotamiento
57
Agotamiento muy rápido, por temperatura elevada > 10° C y retención reducida
(50 mm/año) hasta retención mediana (80mm/año):
Donde:
58
Almacenamiento hídrico
Donde:
LL = 500 (mm/año)
Abastecimiento de la retención
59
Ilustración 14 Almacenamiento hídrico durante época de lluvias (coeficientes ai %)
Según Villon, 2002 (pág. 241) el caudal maximo sirve para diseñar las siguientes
estructuras:
Dimensiones de un cauce
Sistemas de drenaje (agrícola, aeropuerto, ciudad y carretera)
Muros de encauzamiento para proteger ciudades y plantaciones
Alcantarillas
Vertederos de demasías
Luz en puentes
Se debe calcular o estimar el caudal de diseño para esos casos. Son los
caudales máximos.
La magnitud del caudal de diseño es función directa del periodo de retorno que
se le asigne, el que a su vez depende de la importancia de la obra y de la vida
útil de esta.
60
Donde:
T: periodo de retorno.
Donde:
N: vida útil de la estructura (en años).
R: riesgo o falla de la estructura.
Ilustración 15 Valores de periodo de retorno
61
Tabla 12 Valores recomendados de riesgo admisible de obras de drenaje
El resultado del proceso del modelaje es el cálculo de los hidrógrafos del flujo en
sitios elegidos de la cuenca del río.
62
El HEC-HMS cuenta con un interfaz que proporciona al usuario:
Disponibilidad de agua
Drenaje urbano
Predicción del flujo
Reducción de los daños de las avenidas.
Operación de sistemas
Pérdidas iniciales
Villon Bejar, (2002), indica que el cálculo de las pérdidas iniciales, se determinan
a través de la siguiente expresión:
Donde:
63
S: infiltración potencial máxima (mm)
Número de curva
Fuente: (Ministerio de Transportes y Comunicaciones, 2011)
64
Tiempo de retraso (Lag time)
Según Villon Bejar, (2002), el tiempo de retraso se estima a partir del tiempo de
concentración:
Donde:
Donde:
K: Coeficiente de Uniformidad
65
Tiempo de Concentración (Tc)
Donde:
Coeficiente de Uniformidad
Donde:
Donde:
Donde:
66
Intensidad de Precipitación ( I )
Donde:
Coeficiente de Escorrentía ( C )
Donde:
Análisis hidráulico
Régimen Permanente
67
Régimen Impermanente
Transporte de Sedimentos
Reporte de resultados
Gráficos
Tablas
VENTAJAS:
68
a un flujo unidimensional es relativamente aceptable para grandes
escalas (ríos y barrancos) con precisiones poco exigentes.
Gran libertad geométrica: Permite el análisis con secciones naturales no
regulares (secciones fluviales: cauce principal y llanuras de inundación).
Es una gran ventaja sobre otro tipo de modelos hidráulicos existentes (y
mucho más rígidos).
Facilidad de creación, modificación y edición de geometrías (entorno
visual muy cómodo y rápido) e introducción de datos de rugosidad y
estructuras transversales (puentes, obras de paso, aliviaderos). Gran
comodidad de visualización de resultados y edición de figuras.
Gran capacidad de importación y exportación de datos en entorno
Windows (comunicación con Excel, Word, AutoCAD) para el post-proceso
de resultados y presentación.
Las nuevas tecnologías SIG (Sistemas de Información Geográfica), tipo
Arc-View, permiten el proceso de grandes cartografías para generar la
geometría del cauce con gran precisión, en formatos importables Hec-Ras
(.geo). Asimismo, existen extensiones para el Post-proceso de láminas de
inundación y mallas (“grids”) de inundación y velocidad.
Uso extendido en todo el mundo y gran experiencia de uso. Hec-Ras es
un modelo bien contrastado, herencia directa (y mejorada) del antiguo
HEC-2 (1984) en MSDOS. • ES GRATIS !!!!!!!(sin licencia).
INCONVENIENTES Y LIMITACIONES:
69
en cuenta la transferencia lateral de momentum del flujo. La distribución
lateral de velocidades no es correcta. En consecuencia, la solución de
flujo en curvas (método de los Flowpaths) es tremendamente aproximada.
Sólo se pueden modelar ríos y barrancos con pendientes menores de 100
(α< 100 , So<0.18) ya que no se tiene en cuenta la componente vertical
del peso de la columna de agua (cos α) en las ecuaciones.
Ni Saltos ni los obstáculos son reproducidos correctamente (Hec-Ras no
realiza balance de fuerzas). Existen métodos para simular o reproducir
dichos efectos, pero son altamente arbitrarios y dependen mucho del tipo
de flujo.
El resultado viene altamente condicionado por las consideraciones
geométricas adoptadas (Trazado de secciones, Areas Inefectivas,
pérdidas de estrechamiento y expansión, etc…). Por tanto, el resultado
sigue siendo bastante “manual” (bajo criterio del calculista)
Problemas en la elección del calado crítico. Dificultad para hallar el calado
crítico en secciones naturales complejas que contienen varios calados
críticos (terrazas fluviales). Por defecto el programa escoge el de calado
más bajo.
Siempre ofrece por defecto una solución, es decir, Hec-ras no se cuelga
nunca. En consecuencia, se debe ser crítico con el resultado numérico.
Limitación en la convergencia de la energía a 40 iteraciones en cada
sección de cálculo. Por defecto adopta el valor de calado de menor error.
Datos Geométricos
Datos de Flujo Permanente
70
Datos geométricos
Esquema de río
Morfología de la red
Secciones transversales
Singularidades
Condiciones de borde
Datos de flujo
71
En casos donde el régimen de flujo cambia de subcrítico a supercrítico, o
viceversa, el usuario debe ejecutar el programa en modo de régimen de
flujo mixto.
Condiciones de borde
2.2.17 INUNDACIONES
72
Inundaciones lentas, que son las que ocurren en las zonas planas de
los ríos y con valles aluviales extensos, los incrementos de nivel diario son
de apenas del orden de centímetros, reporta afectaciones de grandes
extensiones, pero usualmente pocas pérdidas de vidas humanas, el
tiempo de afectación puede fácilmente llegar a ser del orden de meses.
Inundaciones súbitas, que aunque las áreas de afectación son
menores, el poder destructivo es potencialmente mayor y cobra el mayor
número de vidas cuando se presentan, responden rápidamente a la
ocurrencia de fuertes precipitaciones en las partes altas de las
cuencas, los incrementos de nivel son del orden de metros en pocas
horas, y el tiempo de permanencia de estas inundaciones en las zonas
afectadas son igualmente de horas o pocos días, estas se presentan en
todas las cuencas de alta pendiente de la región Andina principalmente
(IDEAM, 2014).
Código: 000677
Tipo: Convencional-Meteorológica
Latitud: 13 33' 8''
Longitud: 72 44' 5''
Estado: Funcionando
73
El monitoreo de esta estación se encuentra a cargo del Servicio Nacional de
Meteorología e Hidrología del Perú – SENAMHI
74
3 CAPÍTULO III: METODOLOGÍA
75
3.2 DISEÑO DE LA INVESTIGACIÓN
Obtención de resultados.
Se modela el caudal
Conclusiones y recomendaciones.
76
3.2.2 DISEÑO DE INGENIERÍA
Fuente: Propia
77
3.3 POBLACIÓN Y MUESTRA
3.3.1 POBLACIÓN
3.3.2 MUESTRA
78
3.3.2.2 Cuantificación de la muestra
Para esta investigación se eligió el periodo comprendido desde el año 1985 hasta
el 2016.
Se considera como datos a los valores numéricos que están dentro del
rango promedio, los datos que estén fuera de contexto serán descartados
y reemplazados por valores hallados estadísticamente.
79
El histórico de datos diarios de precipitación que estén dentro del periodo
1985-2016.
El histórico de datos mensuales de precipitación que estén dentro del
periodo 1985-2016.
3.4 INSTRUMENTOS
(mm)
Estacion
: Latitud Departamento
Tipo : Longitud Provincia
Altitud Distrito
Año Dia ENE FEB MAR ABR MAY JUN JUL AGO SET OCT NOV DIC
1985 1
1985 2
1985 3
1985 4
1985 5
1985 6
1985 7
1985 8
1985 9
1985 10
1985 11
1985 12
1985 13
1985 14
1985 15
1985 16
1985 17
1985 18
1985 19
1985 20
1985 21
1985 22
1985 23
1985 24
1985 25
1985 26
1985 27
1985 28
1985 29
1985 30
1985 31
Fuente: PROPIA
80
3.4.1.2 Precipitación mensual
Ficha de datos de la precipitación elaborada para la obtención de datos del
SENAMHI
81
Tabla 15 Formato de registro de precipitación mensual
(mm)
82
3.4.2 INSTRUMENTOS DE INGENIERIA
Estación meteorológica.
Estación total
3.5.1 PRECIPITACIÓN
83
Tabla 16 Registro de precipitación diaria (ejemplo: año 1985)
84
Tabla 17 Registro de precipitación mensual
N°
AÑO ENE FEB MAR ABR MAY JUN JUL AGO SET OCT NOV DIC TOTAL
REG
1 1,964 42.0 90.0 162.0 33.6 13.0 0.0 5.2 0.0 60.0 22.1 34.3 42.6 504.8
2 1,965 85.2 158.7 113.0 48.0 8.0 0.0 5.0 0.0 10.2 0.0 37.9 135.8 601.8
3 1,966 83.9 125.9 86.0 16.5 31.4 0.0 0.0 4.0 18.8 136.4 140.6 53.5 697.0
4 1,967 87.7 115.5 252.1 66.6 15.6 0.2 14.7 6.3 15.3 68.0 35.5 154.5 832.0
5 1,968 140.2 112.1 139.4 14.5 0.5 0.0 6.7 5.8 4.7 42.1 131.5 61.1 658.6
6 1,969 100.7 115.4 226.9 73.6 0.8 3.3 3.5 2.5 9.0 93.9 90.2 155.4 875.2
7 1,970 171.4 74.2 73.4 44.6 10.8 1.2 0.0 0.0 28.6 54.2 39.6 130.5 628.5
8 1,971 76.0 2.0 0.0 0.0 0.0 2.0 38.0 84.8 58.2
9 1,972 117.3 121.7 98.5 48.5 0.4 0.0 12.9 57.9 12.6 20.9 23.2 90.5 604.4
10 1,973 87.7 74.9 128.9 43.7 2.8 0.0 3.9 3.0 20.9 23.2 90.5 35.8 515.3
11 1,974 80.3 108.2 104.4 31.6 8.0 7.9 5.0 14.3 0.5 11.8 17.6 30.3 419.9
12 1,975 26.6 42.4 60.6 30.6 4.7 2.0 1.3 4.8 9.5 20.6 4.3 197.0 404.4
13 1,976 159.0 119.2 141.8 49.2 0.5 7.9 2.0 7.3 30.4 27.8 25.2 65.1 635.4
14 1,977 75.0 121.1 107.1 11.6 5.1 0.0 3.8 0.0 10.6 81.6 133.2 113.5 662.6
15 1,978 124.0 170.4 93.2 25.1 1.1 0.0 0.3 16.6 20.0 145.0 77.7
16 1,979 42.3 104.3 116.8 35.9 4.5 0.0 1.2 4.8 11.7 13.4
17 1,980 46.5 99.4 83.1 9.6 4.1 0.0 5.6 0.0 0.0 26.1 63.1 137.3 474.8
18 1,981 114.7 122.7 63.7 46.8 0.0 6.0 0.9 18.3 24.5 53.0 145.7 173.1 769.4
19 1,982 206.4 177.2 113.1 30.4 0.0 3.4 0.0 13.6 4.7 8.7 125.1 55.2 737.8
20 1,983 138.5 87.5 67.8 45.8 4.5 0.0 0.0 0.0 0.0 53.7 93.8
21 1,984 182.1 157.9 0.0 0.0 0.0 0.0
22 1,985 93.3 35.9 10.9 2.5 0.3 9.5 33.4 78.3 99.8
23 1,986 98.3 96.9 157.5 34.0 6.8 0.0 0.0 9.2 29.9 26.8 34.4 54.4 548.2
24 1,987 150.4 81.4 27.1 54.1 11.6 0.0 5.0 0.0 41.3 43.1 63.9
25 1,988 178.6 113.6 124.6 71.6 9.0 0.0 0.0 0.0 0.0 32.6
26 1,989
27 1,990 37.1 48.5 27.1 24.1 20.8 0.0 0.0 0.0 0.0 21.1
28 1,991 24.7 81.1 25.0 12.6 0.0 0.0 0.0 0.0 0.0 0.0 61.3
29 1,992 40.0 61.3 29.0 37.9 0.0 0.0 0.0 24.7 10.8 42.1 54.0 85.0 384.8
30 1,993 164.5 113.2 119.4 71.0 6.4 0.0 5.1 38.8 0.0 40.1 124.0 164.6 847.1
31 1,994 157.6 81.3 59.7 19.0 0.0 0.0 0.0 0.0 0.0 229.3 95.6 155.5 798.0
32 1,995 129.1 89.7 139.5 34.5 6.8 0.0 0.0 0.0 5.9 4.4 105.8 119.2 634.9
33 1,996 120.9 165.8 77.1 48.7 0.8 0.0 0.0 21.4 17.9 64.1 62.2 63.7 642.6
34 1,997 171.3 136.6 115.1 21.2 18.2 0.0 0.0 19.2 4.7 14.2 91.3 129.3 721.1
35 1,998 176.0 106.6 85.7 12.7 2.2 3.9 0.0 1.7 2.2 35.1 35.8 99.3 561.2
36 1,999 182.3 171.1 119.8 60.0 8.1 3.8 2.0 0.0 27.7 64.9 65.1 107.0 811.8
37 2,000 202.8 180.5 85.8 22.2 7.1 11.2 16.0 14.1 15.5 51.5 40.5 124.9 772.1
38 2,001 197.2 68.5 127.9 23.3 13.5 0.0 12.6 23.1 4.0 48.5 81.1 69.3 669.0
39 2,002 114.3 158.8 107.2 87.4 15.0 7.2 32.2 12.9 27.4 102.4 108.0 106.0 878.8
40 2,003 168.2 239.4 198.1 43.3 6.3 2.0 0.4 14.4 16.6 33.0 112.0 145.2 978.9
41 2,004 123.8 134.1 67.7 47.9 17.9 4.0 14.5 12.4 43.9 46.9 75.8 171.7 760.6
42 2,005 77.7 97.9 102.0 28.7 2.7 0.0 11.2 2.1 8.7 31.5 66.9 141.4 570.8
43 2,006 183.0 116.9 103.9 77.0 0.0 6.0 0.0 1.1 3.1 51.4 87.8 120.0 750.2
44 2,007 75.6 109.1 151.6 19.1 7.0 0.0 7.8 1.5 1.0 47.4 84.1 174.5 678.7
45 2,008 107.5 127.3 109.4 18.7 8.4 6.3 0.0 11.1 5.5 48.2 93.3 152.1 687.8
46 2,009 84.2 185.3 86.5 47.1 13.0 0.0 4.6 5.8 13.5 36.5 119.4 95.0 690.9
47 2,010 175.5 91.2 114.6 35.5 7.1 0.0 1.8 4.6 9.9 43.0 58.7 113.0 654.9
48 2,011 145.9 185.7 152.4 52.2 5.4 1.5 6.8 18.6 33.4 37.2 90.4 102.8 832.3
49 2,012 135.9 240.1 209.5 69.2 0.0 5.6 7.0 1.4 16.8 57.0 36.8 233.6 1012.9
50 2,013 113.2 186.9 64.6 22.1 3.7. 1.0 2.3 27.0 7.9 51.7 124.2 146.8 747.7
51 2,014 182.0 121.4 134.9 60.4 5.0 0.0 2.1 1.6 21.4 43.1. 52.7 169.1 750.6
N° Datos 48 48 47 49 49 48 50 50 48 48 46 45 40
Media 121.4 122.7 109.7 41.4 7.2 2.0 4.1 8.2 13.1 43.7 76.2 112.6 662.3
Desv. Estandar 51.27 44.36 49.09 21.42 7.79 3.11 6.03 11.50 12.66 38.15 37.72 46.85 145.06
Coef. Variacion 42.22 36.16 44.76 51.69 108.21 156.81 146.54 140.23 96.83 87.31 49.52 41.60 21.90
Prec. Max. 206.4 240.1 252.1 93.3 35.9 11.2 32.2 57.9 60.0 229.3 145.7 233.6 252.1
Prec. Min. 24.7 42.4 25.0 9.6 0.0 0.0 0.0 0.0 0.0 0.0 4.3 30.3 0.0
Fuente: SENAMHI
85
3.6 PROCEDIMIENTOS DE ANÁLISIS DE DATOS
3.6.1 UBICACIÓN
El rio Asmayacu vierte sus aguas al rio Apurímac. Políticamente la cuenca de
Asmayacu, se encuentra ubicada en el Departamento de Apurímac, Provincia de
Abancay y Distrito de Curahuasi
86
3.6.2 PROCESAMIENTO DE DATOS DE PRECIPITACION
Fuente: PROPIA
87
Ilustración 19 Delimitación de la cuenca asmayacu
88
3.6.3.2 Cálculo del área de la cuenca
Fuente: propia
89
El área de la cuenca Asmayacu tiene un valor de:
Con los datos obtenidos en el ArcGIS como son áreas secciones según sus
altitudes se realizan los cálculos para la obtención de la curva HIPSOMETRICA
90
Ilustración 21 Áreas según sus altitudes
Fuente propia
91
Una vez obtenidos las áreas con altitudes diferentes se procede a hallar los
porcentajes que estas áreas representan del total de la cuenca, para ser
procesados en el cuadro siguiente:
Fuente propia
Con los porcentajes de las áreas según su altitud se grafica la curva hipsométrica
y la frecuencia de altitudes, de la cuenca en estudio.
Fuente propia
El mayor porcentaje de área esta entre las altitudes de 4100 y 4300 msnm.
92
3.6.3.5 Cálculo de la cota media de la cuenca
E: Elevación media
𝟏𝟏𝟐𝟗𝟓𝟏. 𝟎𝟒
𝑬=
𝟐𝟗. 𝟕𝟕
E = 3794.48 msnm
𝑨
𝑲𝒇 =
𝑳𝟐
Se conoce que:
Área: A = 29.758 km2
Longitud al punto más lejano: L = 9.793 km
93
Ilustración 23 Calculo del factor de forma
Fuente: PROPIA
Donde:
𝟐𝟗. 𝟕𝟓𝟖
𝑲𝒇 =
𝟗. 𝟕𝟗𝟑𝟐
Kf = 0.310
Según Villon esta cuenca tiene menos tendencia a concentrar las intensidades
en forma simultánea en toda la cuenca por tener un bajo valor de Kf.
𝑷
𝑲𝒄 = 𝟎. 𝟐𝟖
√𝑨
Se conoce que:
Área: A = 29.758144 km2
Perímetro: P = 26.952185 km2
Donde:
𝟐𝟔. 𝟗𝟓
𝑲𝒄 = 𝟎. 𝟐𝟖
√𝟐𝟗. 𝟕𝟓𝟖
Kc=1.527428718
94
El valor indica que la cuenca tiene una forma alargada porque el valor de Kc de
esta cuenca es mayor a 1, si fuese igual a 1 entonces tendría una forma circular.
l = 2.381 km
El producto de los lados del rectángulo debe dar por respuesta el área total de
la cuenca, de la misma forma la suma de los lados de este debe ser igual al
perímetro de la cuenca.
95
Tabla 19 Cálculo de los valores li del rectángulo equivalente
Fuente propia
Fuente propia
Los colores del rectángulo representan las áreas que componen la cuenca
diferenciadas por sus altitudes y los valores expresados dentro representan la
longitud parcial del lado mayor del rectángulo. Al sumar estos valores
obtendremos el lado mayor del rectángulo el cual es 12.498 km.
96
3.6.3.7 Cálculo de la pendiente de la cuenca
Donde:
Sc = Pendiente de la cuenca (%) :
D = Desnivel entre curvas de nivel (km)
L = Longitud total de las curvas de nivel dentro de la cuenca (km)
A = Área de la cuenca (km2)
Se conoce que:
D = 0.2 km
A = 29.758144 km2
Ilustración 25 Longitud de franjas
Fuente: PROPIA
97
Tabla 20 Longitudes de las franjas de las curvas de nivel de la cuenca según su altitud
Fuente propia
𝐒𝐜 = 𝟑𝟐. 𝟕𝟒𝟗%
𝑯
𝑺=
𝑳
Donde:
Se tiene que:
98
H=1450.295 msnm.
L= 11.997 km.
1450.295
𝑆= ∗ 100%
11.997
S=12.089 %
99
Ilustración 26 Número de orden de la cuenca
Fuente: PROPIA-ARCGIS
100
3.6.4 PLUVIOMETRÍA
Fuente: SENAMHI
101
Ilustración 28 Cuenca Asmayacu y las estaciones meteorológicas más cercanas
Fuente: PROPIA
102
Ilustración 29 Polígonos de Thiessen de las 4 estaciones adyacentes a la cuenca
Fuente: PROPIA
103
Usando el método de los polígonos de Thiessen para determinar la influencia de
las estaciones se tiene que solo la estación de Curahuasi es la que tiene
influencia sobre la cuenca de Asmayacu, debido a esto las estaciones de San
Antonio, Cunyac y Curpahuasi son desestimadas.
104
Tabla 21 Precipitación total mensual (mm) - Estación Curahuasi
N°
AÑO ENE FEB MAR ABR MAY JUN JUL AGO SET OCT NOV DIC TOTAL
REG
1 1,964 42.0 90.0 162.0 33.6 13.0 0.0 5.2 0.0 60.0 22.1 34.3 42.6 504.8
2 1,965 85.2 158.7 113.0 48.0 8.0 0.0 5.0 0.0 10.2 0.0 37.9 135.8 601.8
3 1,966 83.9 125.9 86.0 16.5 31.4 0.0 0.0 4.0 18.8 136.4 140.6 53.5 697.0
4 1,967 87.7 115.5 252.1 66.6 15.6 0.2 14.7 6.3 15.3 68.0 35.5 154.5 832.0
5 1,968 140.2 112.1 139.4 14.5 0.5 0.0 6.7 5.8 4.7 42.1 131.5 61.1 658.6
6 1,969 100.7 115.4 226.9 73.6 0.8 3.3 3.5 2.5 9.0 93.9 90.2 155.4 875.2
7 1,970 171.4 74.2 73.4 44.6 10.8 1.2 0.0 0.0 28.6 54.2 39.6 130.5 628.5
8 1,971 76.0 2.0 0.0 0.0 0.0 2.0 38.0 84.8 58.2 261.0
9 1,972 117.3 121.7 98.5 48.5 0.4 0.0 12.9 57.9 12.6 20.9 23.2 90.5 604.4
10 1,973 87.7 74.9 128.9 43.7 2.8 0.0 3.9 3.0 20.9 23.2 90.5 35.8 515.3
11 1,974 80.3 108.2 104.4 31.6 8.0 7.9 5.0 14.3 0.5 11.8 17.6 30.3 419.9
12 1,975 26.6 42.4 60.6 30.6 4.7 2.0 1.3 4.8 9.5 20.6 4.3 197.0 404.4
13 1,976 159.0 119.2 141.8 49.2 0.5 7.9 2.0 7.3 30.4 27.8 25.2 65.1 635.4
14 1,977 75.0 121.1 107.1 11.6 5.1 0.0 3.8 0.0 10.6 81.6 133.2 113.5 662.6
15 1,978 124.0 170.4 93.2 25.1 1.1 0.0 0.3 16.6 20.0 145.0 77.7 673.4
16 1,979 42.3 104.3 116.8 35.9 4.5 0.0 1.2 4.8 11.7 13.4 334.9
17 1,980 46.5 99.4 83.1 9.6 4.1 0.0 5.6 0.0 0.0 26.1 63.1 137.3 474.8
18 1,981 114.7 122.7 63.7 46.8 0.0 6.0 0.9 18.3 24.5 53.0 145.7 173.1 769.4
19 1,982 206.4 177.2 113.1 30.4 0.0 3.4 0.0 13.6 4.7 8.7 125.1 55.2 737.8
20 1,983 138.5 87.5 67.8 45.8 4.5 0.0 0.0 0.0 0.0 53.7 93.8 491.6
21 1,984 182.1 157.9 0.0 0.0 0.0 0.0 340.0
22 1,985 93.3 35.9 10.9 2.5 0.3 9.5 33.4 78.3 99.8 363.9
23 1,986 98.3 96.9 157.5 34.0 6.8 0.0 0.0 9.2 29.9 26.8 34.4 54.4 548.2
24 1,987 150.4 81.4 27.1 54.1 11.6 0.0 5.0 0.0 41.3 43.1 63.9 477.9
25 1,988 178.6 113.6 124.6 71.6 9.0 0.0 0.0 0.0 0.0 32.6 530.0
26 1,989
27 1,990 37.1 48.5 27.1 24.1 20.8 0.0 0.0 0.0 0.0 21.1 178.7
28 1,991 24.7 81.1 25.0 12.6 0.0 0.0 0.0 0.0 0.0 0.0 61.3 204.7
29 1,992 40.0 61.3 29.0 37.9 0.0 0.0 0.0 24.7 10.8 42.1 54.0 85.0 384.8
30 1,993 164.5 113.2 119.4 71.0 6.4 0.0 5.1 38.8 0.0 40.1 124.0 164.6 847.1
31 1,994 157.6 81.3 59.7 19.0 0.0 0.0 0.0 0.0 0.0 229.3 95.6 155.5 798.0
32 1,995 129.1 89.7 139.5 34.5 6.8 0.0 0.0 0.0 5.9 4.4 105.8 119.2 634.9
33 1,996 120.9 165.8 77.1 48.7 0.8 0.0 0.0 21.4 17.9 64.1 62.2 63.7 642.6
34 1,997 171.3 136.6 115.1 21.2 18.2 0.0 0.0 19.2 4.7 14.2 91.3 129.3 721.1
35 1,998 176.0 106.6 85.7 12.7 2.2 3.9 0.0 1.7 2.2 35.1 35.8 99.3 561.2
36 1,999 182.3 171.1 119.8 60.0 8.1 3.8 2.0 0.0 27.7 64.9 65.1 107.0 811.8
37 2,000 202.8 180.5 85.8 22.2 7.1 11.2 16.0 14.1 15.5 51.5 40.5 124.9 772.1
38 2,001 197.2 68.5 127.9 23.3 13.5 0.0 12.6 23.1 4.0 48.5 81.1 69.3 669.0
39 2,002 114.3 158.8 107.2 87.4 15.0 7.2 32.2 12.9 27.4 102.4 108.0 106.0 878.8
40 2,003 168.2 239.4 198.1 43.3 6.3 2.0 0.4 14.4 16.6 33.0 112.0 145.2 978.9
41 2,004 123.8 134.1 67.7 47.9 17.9 4.0 14.5 12.4 43.9 46.9 75.8 171.7 760.6
42 2,005 77.7 97.9 102.0 28.7 2.7 0.0 11.2 2.1 8.7 31.5 66.9 141.4 570.8
43 2,006 183.0 116.9 103.9 77.0 0.0 6.0 0.0 1.1 3.1 51.4 87.8 120.0 750.2
44 2,007 75.6 109.1 151.6 19.1 7.0 0.0 7.8 1.5 1.0 47.4 84.1 174.5 678.7
45 2,008 107.5 127.3 109.4 18.7 8.4 6.3 0.0 11.1 5.5 48.2 93.3 152.1 687.8
46 2,009 84.2 185.3 86.5 47.1 13.0 0.0 4.6 5.8 13.5 36.5 119.4 95.0 690.9
47 2,010 175.5 91.2 114.6 35.5 7.1 0.0 1.8 4.6 9.9 43.0 58.7 113.0 654.9
48 2,011 145.9 185.7 152.4 52.2 5.4 1.5 6.8 18.6 33.4 37.2 90.4 102.8 832.3
49 2,012 135.9 240.1 209.5 69.2 0.0 5.6 7.0 1.4 16.8 57.0 36.8 233.6 1012.9
50 2,013 113.2 186.9 64.6 22.1 3.7 1.0 2.3 27.0 7.9 51.7 124.2 146.8 751.4
51 2,014 182.0 121.4 134.9 60.4 5.0 0.0 2.1 1.6 21.4 43.1 52.7 169.1 793.7
52 2,015 172.50 141.80 150.70 55.90 10.00 1.10 0.90 16.90 2.30 28.80 70.60 109.50 761.0
53 2,016 64.00 193.60 46.80 39.20 1.30 4.20 4.90 17.2 21.3 65.0 457.5
FUENTE: SENAMHI
105
Ilustración 30 Histograma Estación Curahuasi
250.0
200.0
150.0
100.0
50.0
0.0
1964
1965
1966
1967
1968
1969
1970
1971
1972
1973
1974
1975
1976
1977
1978
1979
1980
1981
1982
1983
1984
1985
1986
1987
1988
1989
1990
1991
1992
1993
1994
1995
1996
1997
1998
1999
2000
2001
2002
2003
2004
2005
2006
2007
2008
2009
2010
2011
2012
2013
2014
2015
2016
Fuente: PROPIA
La ilustración Nro. 32 muestra picos en los años 1967, 1969, 1995, 2003, 2012 y 2013, también muestra valores bajos en el periodo
1974-1975 y 1989-1992. Estos valores deben comprobarse mediante un análisis estadístico, si son fenómenos naturales que
efectivamente han ocurrido, o si son producto de errores sistemáticos.
106
3.6.5.2 Análisis estadístico de saltos y tendencias
Para realizar dicho análisis, dividimos los datos de precipitación, en dos series,
los cuales describimos a continuación:
Serie 1: La serie 1 está compuesta por todos los datos de precipitación mensual,
que van desde enero de 1964 hasta octubre de 1988. Esta serie cuenta con
281 datos.
Serie 2: La serie 2 está compuesta por todos los datos de precipitación mensual,
que van desde enero de 1990 hasta octubre de 2016. Esta serie cuenta con
319 datos.
107
Tabla 22 Precipitación mensual (mm) – Estación Curahuasi: Serie 1 y 2
N°
AÑO ENE FEB MAR ABR MAY JUN JUL AGO SET OCT NOV DIC TOTAL
REG
1 1,964 42.0 90.0 162.0 33.6 13.0 0.0 5.2 0.0 60.0 22.1 34.3 42.6 504.8
2 1,965 85.2 158.7 113.0 48.0 8.0 0.0 5.0 0.0 10.2 0.0 37.9 135.8 601.8
3 1,966 83.9 125.9 86.0 16.5 31.4 0.0 0.0 4.0 18.8 136.4 140.6 53.5 697.0
4 1,967 87.7 115.5 252.1 66.6 15.6 0.2 14.7 6.3 15.3 68.0 35.5 154.5 832.0
5 1,968 140.2 112.1 139.4 14.5 0.5 0.0 6.7 5.8 4.7 42.1 131.5 61.1 658.6
6 1,969 100.7 115.4 226.9 73.6 0.8 3.3 3.5 2.5 9.0 93.9 90.2 155.4 875.2
7 1,970 171.4 74.2 73.4 44.6 10.8 1.2 0.0 0.0 28.6 54.2 39.6 130.5 628.5
8 1,971 76.0 2.0 0.0 0.0 0.0 2.0 38.0 84.8 58.2 261.0
9 1,972 117.3 121.7 98.5 48.5 0.4 0.0 12.9 57.9 12.6 20.9 23.2 90.5 604.4
10 1,973 87.7 74.9 128.9 43.7 2.8 0.0 3.9 3.0 20.9 23.2 90.5 35.8 515.3
11 1,974 80.3 108.2 104.4 31.6 8.0 7.9 5.0 14.3 0.5 11.8 17.6 30.3 419.9
12 1,975 26.6 42.4 60.6 30.6 4.7 2.0 1.3 4.8 9.5 20.6 4.3 197.0 404.4
13 1,976 159.0 119.2 141.8 49.2 0.5 7.9 2.0 7.3 30.4 27.8 25.2 65.1 635.4
14 1,977 75.0 121.1 107.1 11.6 5.1 0.0 3.8 0.0 10.6 81.6 133.2 113.5 662.6
15 1,978 124.0 170.4 93.2 25.1 1.1 0.0 0.3 16.6 20.0 145.0 77.7 673.4
16 1,979 42.3 104.3 116.8 35.9 4.5 0.0 1.2 4.8 11.7 13.4 334.9
17 1,980 46.5 99.4 83.1 9.6 4.1 0.0 5.6 0.0 0.0 26.1 63.1 137.3 474.8
18 1,981 114.7 122.7 63.7 46.8 0.0 6.0 0.9 18.3 24.5 53.0 145.7 173.1 769.4
19 1,982 206.4 177.2 113.1 30.4 0.0 3.4 0.0 13.6 4.7 8.7 125.1 55.2 737.8
20 1,983 138.5 87.5 67.8 45.8 4.5 0.0 0.0 0.0 0.0 53.7 93.8 491.6
21 1,984 182.1 157.9 0.0 0.0 0.0 0.0 340.0
22 1,985 93.3 35.9 10.9 2.5 0.3 9.5 33.4 78.3 99.8 363.9
23 1,986 98.3 96.9 157.5 34.0 6.8 0.0 0.0 9.2 29.9 26.8 34.4 54.4 548.2
24 1,987 150.4 81.4 27.1 54.1 11.6 0.0 5.0 0.0 41.3 43.1 63.9 477.9
25 1,988 178.6 113.6 124.6 71.6 9.0 0.0 0.0 0.0 0.0 32.6 530.0
26 1,989
27 1,990 37.1 48.5 27.1 24.1 20.8 0.0 0.0 0.0 0.0 21.1 178.7
28 1,991 24.7 81.1 25.0 12.6 0.0 0.0 0.0 0.0 0.0 0.0 61.3 204.7
29 1,992 40.0 61.3 29.0 37.9 0.0 0.0 0.0 24.7 10.8 42.1 54.0 85.0 384.8
30 1,993 164.5 113.2 119.4 71.0 6.4 0.0 5.1 38.8 0.0 40.1 124.0 164.6 847.1
31 1,994 157.6 81.3 59.7 19.0 0.0 0.0 0.0 0.0 0.0 229.3 95.6 155.5 798.0
32 1,995 129.1 89.7 139.5 34.5 6.8 0.0 0.0 0.0 5.9 4.4 105.8 119.2 634.9
33 1,996 120.9 165.8 77.1 48.7 0.8 0.0 0.0 21.4 17.9 64.1 62.2 63.7 642.6
34 1,997 171.3 136.6 115.1 21.2 18.2 0.0 0.0 19.2 4.7 14.2 91.3 129.3 721.1
35 1,998 176.0 106.6 85.7 12.7 2.2 3.9 0.0 1.7 2.2 35.1 35.8 99.3 561.2
36 1,999 182.3 171.1 119.8 60.0 8.1 3.8 2.0 0.0 27.7 64.9 65.1 107.0 811.8
37 2,000 202.8 180.5 85.8 22.2 7.1 11.2 16.0 14.1 15.5 51.5 40.5 124.9 772.1
38 2,001 197.2 68.5 127.9 23.3 13.5 0.0 12.6 23.1 4.0 48.5 81.1 69.3 669.0
39 2,002 114.3 158.8 107.2 87.4 15.0 7.2 32.2 12.9 27.4 102.4 108.0 106.0 878.8
40 2,003 168.2 239.4 198.1 43.3 6.3 2.0 0.4 14.4 16.6 33.0 112.0 145.2 978.9
41 2,004 123.8 134.1 67.7 47.9 17.9 4.0 14.5 12.4 43.9 46.9 75.8 171.7 760.6
42 2,005 77.7 97.9 102.0 28.7 2.7 0.0 11.2 2.1 8.7 31.5 66.9 141.4 570.8
43 2,006 183.0 116.9 103.9 77.0 0.0 6.0 0.0 1.1 3.1 51.4 87.8 120.0 750.2
44 2,007 75.6 109.1 151.6 19.1 7.0 0.0 7.8 1.5 1.0 47.4 84.1 174.5 678.7
45 2,008 107.5 127.3 109.4 18.7 8.4 6.3 0.0 11.1 5.5 48.2 93.3 152.1 687.8
46 2,009 84.2 185.3 86.5 47.1 13.0 0.0 4.6 5.8 13.5 36.5 119.4 95.0 690.9
47 2,010 175.5 91.2 114.6 35.5 7.1 0.0 1.8 4.6 9.9 43.0 58.7 113.0 654.9
48 2,011 145.9 185.7 152.4 52.2 5.4 1.5 6.8 18.6 33.4 37.2 90.4 102.8 832.3
49 2,012 135.9 240.1 209.5 69.2 0.0 5.6 7.0 1.4 16.8 57.0 36.8 233.6 1012.9
50 2,013 113.2 186.9 64.6 22.1 3.7 1.0 2.3 27.0 7.9 51.7 124.2 146.8 751.4
51 2,014 182.0 121.4 134.9 60.4 5.0 0.0 2.1 1.6 21.4 43.1 52.7 169.1 793.7
52 2,015 172.50 141.80 150.70 55.90 10.00 1.10 0.90 16.90 2.30 28.80 70.60 109.50 761.0
53 2,016 64.00 193.60 46.80 39.20 1.30 4.20 4.90 17.2 21.3 65.0 457.5
Fuente: SENAMHI
108
Análisis de saltos
a) Consistencia en la media
Cálculo de la media y de la desviación estándar para las series 1 y 2
𝑛1 = 281
Xprom= x ̅= 49.97
S 1= S 1= 54.56
𝑛2 =319
(Xprom)= x ̅= 57.95
S 1= S 2= 60.48
Sp= 57.786
Sp= Sp=
109
Reemplazando los valores obtenemos el valor de tc:
tc = -1.688
Entonces:
V=281+319-2 V=598
α=0.050
Teniendo en cuenta estos datos procedemos a hallar el valor de t tabular en la
tabla Nro. 23 de t de Student.
110
Tabla 23 Tabla de distribución “t – Student”
111
Comparando tc con tt
I tc I < tt
1.6881 1.9626 Entonces S1(x) = S2(x) estadísticamente
2
1 = 2977.015
2
2 = 3558.150
2 2
= ; si 1 𝑥 2 (𝑥)
= = 1.1
; si 𝟐 𝟐
= 𝟐 𝟏( )
𝟐
Se verifica que
𝟐( )> 𝟐𝟏 ( ) por lo tanto
2 2
= ; si 1 𝑥 2 (𝑥)
= = 1.195
; si 𝟐 𝟐
= 𝟐 𝟏( )
112
Se calculó el valor de “Ft” en la tabla F-Fisher considerando que:
G.L.N= 𝑛1 − 1
2 2
sí 1 𝑥 2 (𝑥)
G.L.D = 𝑛2 − 1
𝑡 = 1.21
G.L.N= 𝑛2 − 1
𝟐 𝟐
𝟏( )
sí
𝟐
G.L.D = 𝑛1 −1
113
Tabla 24 Valores para la distribución F de Fischer (a = 0.05)
𝑡 = 1.2112
𝑡
1.195 < 1.2112 Entonces S1(x) = S2(x) estadísticamente
Por lo tanto, no es necesario corregir los datos
114
Análisis de tendencias
a) Tendencia en la Media
Forma general de la ecuación polinomial
𝑇𝑚 = 𝐴𝑚 + 𝐵𝑚 𝑡 + 𝐶𝑚 𝑡 2 + 𝐷𝑚 𝑡 3 + ⋯.
donde:
t: tiempo en años, tomando la variable independiente de la tendencia
t: 1, 2, 3, …n
Tm: Tendencia en la media, para este caso
Tm: Datos a usarse para el cálculo de parámetros
Am, Bm, Cm, Dm: coeficientes de los polinomios de regresión que deben
ser estimados con los datos
n = n 1 + n2
𝑇̅𝑚 = 54.215
𝑆 = 57.738
115
Tabla 25 Valores registrables entre la precipitación mensual (mm) Vs. el tiempo (1er bloque)
116
Tabla 26 Valores registrables entre la precipitación mensual (mm) Vs. el tiempo (2do bloque)
117
𝑡 ̅ = 300.5
𝑆𝑡 = 173.35
Se tiene que:
𝑡 ∗ 𝑇𝑚 =17096.59
R R
= 0.080
𝑡𝑐 = 1.9626135
118
A continuación, calculamos el valor de “t tabular” con el 95% de probabilidad el
cual se calcula mediante la tabla de Student con:
tt = 1.964472
Si 𝑡𝑐 𝑡𝑡 entonces R no es significativo
𝑡𝑐 𝑡𝑡 entonces R si es significativo
Se tiene que:
I tc I < tt
1.9626 1.9645
Entonces R no es significativo
por lo tanto, no es necesario corregir los datos
𝐵𝑚 = 0.0266457
𝐴𝑚 = 𝑇̅𝑚 − 𝑡̅ ∗ 𝐵𝑚 𝐴𝑚 = 46.208286
119
Ilustración 31 Histograma con línea de tendencia en la media – Estación Curahuasi
Fuente: PROPIA
120
b) Tendencia en la Desviación estándar
La tendencia en la desviación estándar es representada por la siguiente
ecuación:
121
Tabla 27 Desviaciones estándar de cada período
𝑆 = 56.751
𝑆𝑠 = 14.628
𝑡̅ = 26.5
122
𝑆𝑡 = 15.155
t Sp t*Sp t Sp t*Sp
1 46.001 46.0014 27 28.142 759.822
2 58.316 116.632 28 26.860 752.077
3 54.494 163.481 29 63.977 1855.325
4 74.969 299.876 30 78.742 2362.265
5 59.251 296.256 31 58.212 1804.562
6 72.451 434.705 32 51.191 1638.111
7 53.780 376.457 33 63.494 2095.302
8 35.772 286.179 34 57.209 1945.112
9 45.752 411.772 35 64.872 2270.532
10 42.829 428.288 36 69.182 2490.556
11 39.409 433.496 37 58.626 2169.166
12 54.671 656.058 38 51.054 1940.054
13 56.458 733.96 39 86.398 3369.530
14 54.560 763.837 40 53.962 2158.496
15 63.438 951.563 41 48.027 1969.125
16 43.133 690.13 42 61.627 2588.324
17 46.531 791.026 43 62.131 2671.615
18 60.378 1086.81 44 56.475 2484.891
19 74.678 1418.89 45 57.499 2587.455
20 48.234 964.683 46 57.213 2631.792
21 88.121 1850.53 47 64.426 3028.027
22 39.846 876.621 48 94.789 4549.890
23 48.483 1115.12 49 64.539 3162.387
24 44.666 1071.99 50 68.045 3402.243
25 59.167 1479.18 51 64.410 3284.930
26 17.411 452.692 52 57.146 2971.597
Fuente: PROPIA
𝑡. 𝑆𝑝 = 1520.374
123
R= R= 0.0743113
𝑡𝑐 = 1.8887423
t t= 2.009
Comparando los valores de t calculado y t tabular
I tc I < tt
1.8887 2.0090
Entonces R no es significativo
por lo que no es necesario corregir los datos
𝐵𝑠 = 0.0717287
𝐴𝑠 = 𝑆̅ − 𝑡 ̅ ∗ 𝐵𝑠 𝐴𝑠 = 54.850145
124
𝑇 = 𝐴𝑠 + 𝐵𝑠 ∗ 𝑡
𝑇 = 54.85 + 0.0717 ∗ 𝑡
125
Ilustración 32 Histograma con línea de tendencia en la Des. Estándar – Estación Curahuasi
Fuente: PROPIA
Del análisis estadístico realizado a la serie histórica de precipitación mensual de la estación Curahuasi, se concluye que, no
presenta saltos ni tendencias, siendo dicha información confiable.
126
3.6.6 COMPLETACIÓN DE DATOS DE PRECIPITACIÓN MENSUAL
CON EASY FIT
a) Distribución Normal
g) Distribución Gumbel
127
3.6.6.2 Completación del histórico de precipitaciones de la estación de
Curahuasi
128
Tabla 29 Datos de precipitación mensual de la estación Curahuasi
N°
AÑO ENE FEB MAR ABR MAY JUN JUL AGO SET OCT NOV DIC TOTAL
REG
1 1,964 42.0 90.0 162.0 33.6 13.0 0.0 5.2 0.0 60.0 22.1 34.3 42.6 504.8
2 1,965 85.2 158.7 113.0 48.0 8.0 0.0 5.0 0.0 10.2 0.0 37.9 135.8 601.8
3 1,966 83.9 125.9 86.0 16.5 31.4 0.0 0.0 4.0 18.8 136.4 140.6 53.5 697.0
4 1,967 87.7 115.5 252.1 66.6 15.6 0.2 14.7 6.3 15.3 68.0 35.5 154.5 832.0
5 1,968 140.2 112.1 139.4 14.5 0.5 0.0 6.7 5.8 4.7 42.1 131.5 61.1 658.6
6 1,969 100.7 115.4 226.9 73.6 0.8 3.3 3.5 2.5 9.0 93.9 90.2 155.4 875.2
7 1,970 171.4 74.2 73.4 44.6 10.8 1.2 0.0 0.0 28.6 54.2 39.6 130.5 628.5
8 1,971 131.8 77.8 32.1 76.0 2.0 0.0 0.0 0.0 2.0 38.0 84.8 58.2 502.6
9 1,972 117.3 121.7 98.5 48.5 0.4 0.0 12.9 57.9 12.6 20.9 23.2 90.5 604.4
10 1,973 87.7 74.9 128.9 43.7 2.8 0.0 3.9 3.0 20.9 23.2 90.5 35.8 515.3
11 1,974 80.3 108.2 104.4 31.6 8.0 7.9 5.0 14.3 0.5 11.8 17.6 30.3 419.9
12 1,975 26.6 42.4 60.6 30.6 4.7 2.0 1.3 4.8 9.5 20.6 4.3 197.0 404.4
13 1,976 159.0 119.2 141.8 49.2 0.5 7.9 2.0 7.3 30.4 27.8 25.2 65.1 635.4
14 1,977 75.0 121.1 107.1 11.6 5.1 0.0 3.8 0.0 10.6 81.6 133.2 113.5 662.6
15 1,978 124.0 170.4 93.2 25.1 1.1 2.8 0.0 0.3 16.6 20.0 145.0 77.7 676.2
16 1,979 42.3 104.3 116.8 35.9 4.5 0.0 1.2 4.8 11.7 13.4 111.9 109.2 556.0
17 1,980 46.5 99.4 83.1 9.6 4.1 0.0 5.6 0.0 0.0 26.1 63.1 137.3 474.8
18 1,981 114.7 122.7 63.7 46.8 0.0 6.0 0.9 18.3 24.5 53.0 145.7 173.1 769.4
19 1,982 206.4 177.2 113.1 30.4 0.0 3.4 0.0 13.6 4.7 8.7 125.1 55.2 737.8
20 1,983 138.5 87.5 67.8 45.8 4.5 1.1 0.0 0.0 0.0 0.0 53.7 93.8 492.7
21 1,984 182.1 157.9 147.3 74.3 0.0 0.0 0.0 0.0 12.5 72.4 120.6 90.2 857.2
22 1,985 129.1 98.1 145.5 93.3 35.9 10.9 2.5 0.3 9.5 33.4 78.3 99.8 736.5
23 1,986 98.3 96.9 157.5 34.0 6.8 0.0 0.0 9.2 29.9 26.8 34.4 54.4 548.2
24 1,987 150.4 81.4 27.1 54.1 11.6 0.0 5.0 0.0 5.0 41.3 43.1 63.9 482.9
25 1,988 178.6 113.6 124.6 71.6 9.0 0.0 0.0 0.0 0.0 32.6 82.7 97.5 710.2
26 1,989 86.8 138.0 75.2 39.0 6.3 3.0 0.3 15.1 30.0 86.8 89.3 124.0 693.7
27 1,990 37.1 48.5 27.1 24.1 20.8 0.0 0.0 0.0 0.0 21.1 50.1 77.3 306.1
28 1,991 24.7 81.1 25.0 12.6 0.0 0.0 0.0 0.0 0.0 0.0 61.3 41.8 246.5
29 1,992 40.0 61.3 29.0 37.9 0.0 0.0 0.0 24.7 10.8 42.1 54.0 85.0 384.8
30 1,993 164.5 113.2 119.4 71.0 6.4 0.0 5.1 38.8 0.0 40.1 124.0 164.6 847.1
31 1,994 157.6 81.3 59.7 19.0 0.0 0.0 0.0 0.0 0.0 229.3 95.6 155.5 798.0
32 1,995 129.1 89.7 139.5 34.5 6.8 0.0 0.0 0.0 5.9 4.4 105.8 119.2 634.9
33 1,996 120.9 165.8 77.1 48.7 0.8 0.0 0.0 21.4 17.9 64.1 62.2 63.7 642.6
34 1,997 171.3 136.6 115.1 21.2 18.2 0.0 0.0 19.2 4.7 14.2 91.3 129.3 721.1
35 1,998 176.0 106.6 85.7 12.7 2.2 3.9 0.0 1.7 2.2 35.1 35.8 99.3 561.2
36 1,999 182.3 171.1 119.8 60.0 8.1 3.8 2.0 0.0 27.7 64.9 65.1 107.0 811.8
37 2,000 202.8 180.5 85.8 22.2 7.1 11.2 16.0 14.1 15.5 51.5 40.5 124.9 772.1
38 2,001 197.2 68.5 127.9 23.3 13.5 0.0 12.6 23.1 4.0 48.5 81.1 69.3 669.0
39 2,002 114.3 158.8 107.2 87.4 15.0 7.2 32.2 12.9 27.4 102.4 108.0 106.0 878.8
40 2,003 168.2 239.4 198.1 43.3 6.3 2.0 0.4 14.4 16.6 33.0 112.0 145.2 978.9
41 2,004 123.8 134.1 67.7 47.9 17.9 4.0 14.5 12.4 43.9 46.9 75.8 171.7 760.6
42 2,005 77.7 97.9 102.0 28.7 2.7 0.0 11.2 2.1 8.7 31.5 66.9 141.4 570.8
43 2,006 183.0 116.9 103.9 77.0 0.0 6.0 0.0 1.1 3.1 51.4 87.8 120.0 750.2
44 2,007 75.6 109.1 151.6 19.1 7.0 0.0 7.8 1.5 1.0 47.4 84.1 174.5 678.7
45 2,008 107.5 127.3 109.4 18.7 8.4 6.3 0.0 11.1 5.5 48.2 93.3 152.1 687.8
46 2,009 84.2 185.3 86.5 47.1 13.0 0.0 4.6 5.8 13.5 36.5 119.4 95.0 690.9
47 2,010 175.5 91.2 114.6 35.5 7.1 0.0 1.8 4.6 9.9 43.0 58.7 113.0 654.9
48 2,011 145.9 185.7 152.4 52.2 5.4 1.5 6.8 18.6 33.4 37.2 90.4 102.8 832.3
49 2,012 135.9 240.1 209.5 69.2 0.0 5.6 7.0 1.4 16.8 57.0 36.8 233.6 1012.9
50 2,013 113.2 186.9 64.6 22.1 3.7 1.0 2.3 27.0 7.9 51.7 124.2 146.8 751.4
51 2,014 182.0 121.4 134.9 60.4 5.0 0.0 2.1 1.6 21.4 43.1 52.7 169.1 793.7
52 2,015 172.50 141.80 150.70 55.90 10.00 1.10 0.90 16.90 2.30 28.80 70.60 109.50 761.0
53 2,016 64.00 193.60 46.80 39.20 1.30 4.20 4.90 17.2 21.3 65.0 141.4 103.3 702.3
Fuente: SENAMHI
129
En el siguiente cuadro se Muestra los datos introducidos al software EasyFit 5.6
para su análisis.
Fuente: PROPIA-EASYFIT
130
Tabla 30 Resultados de la prueba de bondad de ajuste (enero)
131
Log Pearson 3
Normal
Log Normal
Fuente: PROPIA-EASYFIT
132
Tabla 31 Resultados de la prueba de bondad de ajuste (enero)
133
Log Pearson 3
Normal
Log Normal
Fuente: PROPIA-EASYFIT
Fuente: PROPIA-EASYFIT
134
Ilustración 37 Histograma de distribución log normal (enero)
Fuente: PROPIA-EASYFIT
Fuente: PROPIA-EASYFIT
135
Ilustración 39 Números aleatorios generados (enero)
Fuente: PROPIA-EASYFIT
136
Tabla 32 Precipitación mensual completada en mm (enero)
AÑO ENE
1,964 42.0
1,965 85.2
1,966 83.9
1,967 87.7
1,968 140.2
1,969 100.7
1,970 171.4
1,971 131.8
1,972 117.3
1,973 87.7
1,974 80.3
1,975 26.6
1,976 159.0
1,977 75.0
1,978 124.0
1,979 42.3
1,980 46.5
1,981 114.7
1,982 206.4
1,983 138.5
1,984 182.1
1,985 129.1
1,986 98.3
1,987 150.4
1,988 178.6
1,989 86.8
1,990 37.1
1,991 24.7
1,992 40.0
1,993 164.5
1,994 157.6
1,995 129.1
1,996 120.9
1,997 171.3
1,998 176.0
1,999 182.3
2,000 202.8
2,001 197.2
2,002 114.3
2,003 168.2
2,004 123.8
2,005 77.7
2,006 183.0
2,007 75.6
2,008 107.5
2,009 84.2
2,010 175.5
2,011 145.9
2,012 135.9
2,013 113.2
2,014 182.0
2,015 172.50
2,016 64.00
Fuente: PROPIA-EASYFIT
137
De la misma forma procedemos a completar los datos del mes de Febrero.
Fuente: PROPIA-EASYFIT
Fuente: PROPIA-EASYFIT
138
Ilustración 42 Histograma de distribución normal (febrero)
Fuente: PROPIA-EASYFIT
Fuente: PROPIA-EASYFIT
139
Ilustración 44 Números aleatorios generados (febrero)
Fuente: PROPIA-EASYFIT
140
Ilustración 45 Precipitación mensual completada en mm (febrero)
AÑO FEB
1,964 90.0
1,965 158.7
1,966 125.9
1,967 115.5
1,968 112.1
1,969 115.4
1,970 74.2
1,971 77.8
1,972 121.7
1,973 74.9
1,974 108.2
1,975 42.4
1,976 119.2
1,977 121.1
1,978 170.4
1,979 104.3
1,980 99.4
1,981 122.7
1,982 177.2
1,983 87.5
1,984 157.9
1,985 98.1
1,986 96.9
1,987 81.4
1,988 113.6
1,989 138.0
1,990 48.5
1,991 81.1
1,992 61.3
1,993 113.2
1,994 81.3
1,995 89.7
1,996 165.8
1,997 136.6
1,998 106.6
1,999 171.1
2,000 180.5
2,001 68.5
2,002 158.8
2,003 239.4
2,004 134.1
2,005 97.9
2,006 116.9
2,007 109.1
2,008 127.3
2,009 185.3
2,010 91.2
2,011 185.7
2,012 240.1
2,013 186.9
2,014 121.4
2,015 141.80
2,016 193.60
Fuente: PROPIA-EASYFIT
141
De la misma forma procedemos a completar los datos de precipitación mensual
para cada mes del año según la mejor distribución probabilística usando el
software EasyFit.
142
Tabla 33 Histórico de datos de precipitación mensual de la estación de Curahuasi completado
Fuente: PROPIA
143
3.6.7 GENERACIÓN DE LOS CAUDALES MENSUALES POR EL
MODELO LUTZ SCHOLZ
Fuente: PROPIA
144
3.6.7.2 Datos básicos de la cuenca de estudio
Este dato fue obtenido por medio del software ArcGIS y su valor es:
Hc = 3774 msnm.
145
Tabla 36 Cálculo de la radiación extraterrestre de la cuenca
Cuenca 13.60 1,030.80 971.40 865.19 728.99 604.98 543.18 566.38 665.38 794.79 919.20 1,007.20 1,039.41
TABLA N° : Numero de Horas de Sol Maxima Media Diaria Probable para diferentes Meses y Latitudes Sur (D)
TABLA N° : Radiación Extraterrestre Media Diaria (RMD) Expresada en Equivalente de Evapotranspiración (mm/día) para diferentes latitudes y meses
Fuente: PROPIA
Ra = 5272.01
146
Temperatura de la estación de referencia (Te)
Te = 16.20 °C.
He = 2763.0 msnm.
Según Plan Meris Inka, (1995) se tiene que la duración de la temporada seca
para esta zona comienza el 1 de Abril y termina el 31 de octubre, es decir
tiene una duración de:
TD = 214 días.
Pc = 663.70
Fuente: PROPIA
147
3.6.7.3 Cálculo de parámetros para la obtención de caudales mensuales
Donde:
148
EP = 0.0075 (0.075(5238.63)√(50)*(38.02)*(1 + 0.06( 3774.69/1000))
EP = 1029.66 mm/año
Fuente: PROPIA
Donde:
C = Coeficiente de escurrimiento
P = Precipitación total anual (mm/año)
EP = Evapotranspiración anual según Hargreaves (mm/año)
r = Coeficiente de correlación
C = 519
149
También se obtuvo el coeficiente de escurrimiento, calculándose previamente
el Déficit de escurrimiento.
Donde:
P = Precipitación total anual (mm/año)
EP = Evapotranspiración anual según Hargreaves (mm/año)
𝑫 = −𝟏𝟑𝟖𝟎 + 𝟎. 𝟖𝟕𝟐 ∗ (𝟔𝟔𝟑. 𝟕𝟎) + 𝟏. 𝟎𝟑𝟐 ∗ (𝟏𝟎𝟐𝟗. 𝟔𝟔)
D= 261.53
Donde:
P = Precipitación total anual (mm/año)
D = Déficit de escurrimiento
Aplicando la formula se tiene:
𝟔𝟔𝟑. 𝟕𝟎 − 𝟐𝟔𝟏. 𝟓𝟑
𝑪=
𝟔𝟔𝟑. 𝟕𝟎
C = 0.606
150
por Lutz Scholz, debido a que la cuenca de estudio, presenta una cobertura
vegetal representativa de la puna (vegetación poco desarrollada).
Donde:
a: Coeficiente de agotamiento
A: Área de la cuenca
R: Retención de la cuenca
𝒂 = 𝟎. 𝟎𝟎𝟐𝟔
A partir del valor del coeficiente de agotamiento “a” se determinó b0 mediante la
siguiente ecuación:
Que muestra que los caudales durante la estación seca, están relacionados al
mes anterior por la siguiente relación, cuando se cuenta el mes por 30 días.
𝑸𝒕
𝒃𝒐 =
𝑸𝒐
Reemplazando se tiene:
Dónde:
Qt = Descarga en el tiempo t
Qo = Descarga inicial
a = Coeficiente de agotamiento
151
t = tiempo
𝒃𝒐 = 𝒆−𝟎.𝟎𝟎𝟐𝟔(𝟑𝟎)
𝒃𝒐 = 𝟎. 𝟗𝟐𝟒𝟗𝟔
Fuente: PROPIA
Fuente: PROPIA
152
Para proseguir con el cálculo de los caudales por el modelo hidrológico Lutz
Scholz, asumimos el valor del coeficiente de agotamiento:
a = 0.000529
Según Scholz los suelos con baja pendiente acumulan mayor cantidad de
láminas de agua y las áreas que tienen pendiente mayor a 15% simplemento no
las considera, an el siguiente cuadro se muestra las láminas de agua que
acumula un suelo según la pendiente del suelo.
Lamina Area
Tipo Descripción Porcentaje (%) Total
(mm/año) (Km2)
Lagunas - Pantanos ---- 500 0.110 55
Pend. < 2% Pend. Cuenca 0 - 2 % 300 6.462 1938
Napa Freatica Pend. < 8% Pend. Cuenca 2 < 8 % 250 3.689 922
Pend. < 15% Pend. Cuenca 8 < 15 % 200 2.108 422
Nevados ---- 500 0 0
Σ= 3337.3
Fuente: PROPIA R= 112.15
153
Ilustración 47 Pendientes de la cuenca
Fuente: PROPIA
Las zonas de mayor retención de aguas son las lagunas, pantanos y los
nevados que tienen una retención de 500 mm/año.
154
Tabla 42 Áreas de los lagos y lagunas
Area de Lagos
Nombre
(km2)
Ccelhua ccocha 0.09148
Tinqui ccocha bajo 0.00434
Tinqui ccocha medio 0.00364
Tinqui ccocha alto 0.00663
Ronto ccocha 0.00351
Area total 0.10959
Fuente: PROPIA
Fuente: PROPIA
155
Se aprecia que los lagos en la cuenca son muy pequeños destacando el lago
Ccelhua ccocha que tiene la mayor extensión.
Lamina Area
Tipo Descripción Porcentaje (%) Total
(mm/año) (Km2)
Lagunas - Pantanos ---- 500 0.110 55
Pend. < 2% Pend. Cuenca 0 - 2 % 300 6.462 1938
Napa Freatica Pend. < 8% Pend. Cuenca 2 < 8 % 250 3.689 922
Pend. < 15% Pend. Cuenca 8 < 15 % 200 2.108 422
Nevados ---- 500 0 0
Σ= 3337.3
R= 112.15
Fuente: PROPIA
R = 112.15 mm/año.
Precipitación efectiva
156
Donde:
𝑪𝟏 + 𝑪𝟐 = 𝟏
Entre la curva I y II
C1 = -1.574 C2 = 2.574
157
C1 = -0.619 C2 = 1.619
C1 = 0.469 C2 = 0.531
Condición 2: C1 + C2 = 1
C1 y C2 fueron obtenidos a partir de las curvas III y IV (PEIII y PEIV), siendo por
tanto los valores para C1 y C2 los obtenidos anteriormente:
C1=0.469 C2=0.531
158
Tabla 44 Cálculo de las precipitaciones efectivas mensuales para el año 1985
PRECIPITACION MENSUAL
MES DIAS TOTAL P PE I PE I PE II PE III PE IV PE
Fuente: PROPIA
Fuente: PROPIA
Gasto de la retención
159
Para calcular el gasto de la retención, primero se determinó los coeficientes bi
para los meses correspondientes a la época seca (de abril a octubre) a partir de
la siguiente relación:
𝒃𝒊 = 𝒃𝒊𝟎
Entonces: 𝒃𝒊 = (𝒆−𝒂𝒕 )𝒊
𝒊
𝒃𝟎
𝑮𝒊 = ( 𝒊
)∗
∑𝒎
𝒊=𝟏 𝒃𝟎
Donde:
R= Retención de la cuenca.
MES i
ABR 1
MAY 2
JUN 3
JUL 4
AGO 5
SET 6
OCT 7
Fuente: PROPIA
160
Tabla 46 Cálculo del gasto mensual (1985)
GASTO (mm/mes)
MES i
Fuente: PROPIA
Abastecimiento de la retención
161
Hacemos una sumatoria de la precipitación mensual para cada año y
ponderamos la incidencia de cada mes, de esta manera obtenemos los
porcentajes para el coeficiente de abastecimiento para cada mes para 30 años
diferentes, los cuales mostramos el cuadro siguiente:
AÑO ENE FEB MAR ABR MAY JUN JUL AGO SET OCT NOV DIC TOTAL
Precipitación 1985 (mm) 129.1 98.1 145.5 93.3 35.9 10.9 2.5 0.3 9.5 33.4 78.3 99.8 736.5
ai (%) 0.175 0.133 0.197 0.127 0.049 0.015 0.003 0.000 0.013 0.045 0.106 0.136 1
Fuente: PROPIA
𝑨𝒊 = 𝒂𝒊 ∗
Donde:
162
ABASTECIMIENTO
Ai (mm/mes)
MES ai
Fuente: PROPIA
Una vez determinados los parámetros del balance hídrico: Precipitación efectiva
(PE), Gasto (G), y Abastecimiento (A) se determinó los caudales para cada mes
con la siguiente ecuación.
𝑪𝑴𝒊 = 𝑷𝑬𝒊 + 𝑮𝒊 − 𝑨𝒊
Dónde:
163
A continuación, se muestra el cálculo de los caudales mensuales del año 1985.
Tabla 50 Calculo de los caudales mensuales por el modelo hidrologico Lutz Scholz
CONTRIBUCION DE LA RETENCION
PRECIPITACION CAUDALES MENSUALES
GASTO ABASTECIMIENTO
MES TOTAL P PE III PE IV PE bi Gi (mm/mes) Ai (mm/mes)
mm/mes mm/mes mm/mes mm/mes ai mm/mes m3/s
JUL 2.5 0.6 9.468 5.312 0.938 16.001 0.003 0.381 20.933 0.233
AGO 0.3 0.1 10.130 5.429 0.923 15.746 0.000 0.046 21.129 0.235
SET 9.5 2.3 7.945 5.300 0.911 15.542 0.013 1.446 19.395 0.223
OCT 33.4 6.8 9.428 8.197 0.893 15.246 0.045 5.085 18.358 0.204
NOV 78.3 23.4 36.316 30.263 0.000 0.000 0.106 11.921 18.342 0.211
DIC 99.8 39.9 56.583 48.765 0.000 0.000 0.135 15.194 33.570 0.373
ENE 129.1 69.4 87.238 78.878 0.000 0.000 0.175 19.655 59.223 0.658
FEB 98.1 38.4 54.880 47.157 0.000 0.000 0.133 14.936 32.221 0.396
MAR 145.5 85.8 104.694 95.839 0.000 0.000 0.198 22.152 73.687 0.819
ABR 93.3 34.3 50.150 42.722 0.985 16.803 0.127 14.205 45.320 0.520
MAY 35.9 7.3 10.161 8.820 0.968 16.526 0.049 5.466 19.880 0.221
JUN 10.9 2.6 7.748 5.335 0.954 16.287 0.015 1.660 19.963 0.229
TOTAL 736.6 310.9 444.7 382 6.571 112.2 1 112.1 382 0.36
Fuente: PROPIA
164
En la siguiente figura se muestra el comportamiento del caudal mensual del rio
Asmayacu durante el año 1985, se puede apreciar que el comportamiento del
caudal es similar al de la precipitación, también se puede diferenciar que
presenta un caudal bajo durante los meses secos (Abril a Octubre) una crecida
considerable en los meses lluviosos con una caída de caudal en el mes de
febrero (como también pasa con la precipitación).
Fuente: PROPIA
Se obtuvieron los caudales para cada mes del año 1985 por el modelo
hidrológico Lutz Scholz, de igual manera se procede a calcular los caudales de
cada mes de los siguientes 30 años en la cuenca.
165
Tabla 51 Generación de caudales mensuales con el modelo Lutz Scholz desde 1985 hasta el 2016
AÑO ENE FEB MAR ABR MAY JUN JUL AGO SET OCT NOV DIC
1985 0.658 0.396 0.819 0.520 0.221 0.229 0.233 0.235 0.223 0.204 0.211 0.373
1986 0.329 0.351 0.876 0.213 0.237 0.266 0.254 0.222 0.197 0.188 0.021 0.060
1987 0.789 0.213 0.014 0.245 0.227 0.280 0.246 0.265 0.248 0.189 0.023 0.091
1988 1.069 0.491 0.554 0.322 0.209 0.218 0.208 0.205 0.209 0.193 0.194 0.297
1989 0.279 0.841 0.195 0.237 0.242 0.257 0.254 0.215 0.212 0.448 0.309 0.622
1990 0.041 0.065 0.036 0.227 0.214 0.187 0.178 0.175 0.178 0.203 0.066 0.222
1991 0.001 0.170 0.000 0.199 0.184 0.187 0.178 0.175 0.178 0.169 0.032 0.000
1992 0.025 0.089 0.020 0.217 0.184 0.187 0.178 0.196 0.192 0.197 0.055 0.236
1993 0.667 0.246 0.265 0.224 0.184 0.187 0.193 0.168 0.178 0.163 0.309 0.668
1994 0.765 0.139 0.046 0.203 0.184 0.187 0.178 0.175 0.178 1.619 0.219 0.745
1995 0.600 0.238 0.700 0.215 0.192 0.187 0.178 0.175 0.186 0.175 0.369 0.497
1996 0.663 1.240 0.195 0.263 0.185 0.187 0.178 0.208 0.208 0.289 0.115 0.117
1997 0.972 0.714 0.432 0.209 0.199 0.187 0.178 0.191 0.184 0.182 0.227 0.576
1998 1.051 0.417 0.190 0.206 0.187 0.192 0.178 0.177 0.182 0.189 0.021 0.304
1999 1.144 1.147 0.532 0.275 0.196 0.194 0.181 0.175 0.207 0.267 0.102 0.395
2000 1.272 1.173 0.184 0.211 0.194 0.200 0.193 0.189 0.194 0.213 0.029 0.530
2001 1.281 0.139 0.611 0.218 0.215 0.241 0.211 0.200 0.225 0.220 0.205 0.130
2002 0.535 1.090 0.464 0.491 0.222 0.235 0.219 0.215 0.217 0.587 0.488 0.452
2003 0.818 1.681 1.112 0.211 0.188 0.188 0.178 0.182 0.186 0.180 0.313 0.594
2004 0.620 0.800 0.149 0.260 0.219 0.245 0.217 0.217 0.235 0.232 0.204 1.093
2005 0.183 0.402 0.408 0.226 0.190 0.187 0.195 0.178 0.194 0.204 0.123 0.861
2006 1.167 0.571 0.376 0.356 0.184 0.197 0.178 0.177 0.184 0.228 0.241 0.551
2007 0.132 0.444 0.828 0.211 0.193 0.187 0.188 0.177 0.179 0.208 0.191 1.041
2008 0.409 0.694 0.429 0.215 0.195 0.197 0.178 0.189 0.187 0.216 0.282 0.868
2009 0.220 1.360 0.240 0.248 0.204 0.187 0.187 0.184 0.199 0.208 0.601 0.314
2010 1.102 0.290 0.500 0.226 0.194 0.187 0.182 0.183 0.194 0.210 0.076 0.481
2011 0.770 1.281 0.833 0.253 0.202 0.205 0.197 0.197 0.209 0.202 0.252 0.347
2012 0.454 1.629 1.169 0.244 0.184 0.190 0.180 0.176 0.182 0.193 0.006 1.407
2013 0.423 1.255 0.085 0.213 0.199 0.202 0.193 0.196 0.196 0.217 0.553 0.750
2014 1.076 0.548 0.626 0.261 0.191 0.188 0.181 0.178 0.197 0.201 0.048 0.953
2015 1.008 0.793 0.801 0.252 0.196 0.189 0.179 0.191 0.182 0.192 0.110 0.393
2016 0.098 1.418 0.049 0.233 0.185 0.195 0.185 0.197 0.203 0.270 0.816 0.381
PROMEDIO 0.644 0.698 0.429 0.253 0.200 0.205 0.195 0.193 0.198 0.271 0.213 0.511
Fuente: PROPIA
166
Con los datos de la tabla anterior se procede a graficar el cuadro que muestra las variaciones de los caudales mensuales generados por el modelo hidrológico Lutz Scholz, indicándose los años que ocurrieron fenómenos del Niño y de la Niña en los últimos 30 años.
Ilustración 51 Variación de los caudales mensuales según el modelo Lutz Scholz en los últimos 30 años
Fuente: PROPIA
En la ilustración 53 el color rojo representa los caudales generados con presencia del Fenómeno del Niño y el color azul representa a los caudales generados con la presencia del fenómeno del fenómeno de la Niña.
Fuente: PROPI
En la ilustración 54 se muestra el periodo donde se presentó un evento de sequía severa en los años 1990-1991-1992, presentándose en este evento de sequía una ocurrencia del Fenómeno del Niño, teniendo una incidencia del 8.33% en los últimos 30 años.
167
Ilustración 53 Incremento notable de caudales con fenómeno la Niña
Fuente: PROPIA
En la ilustración anterior mostramos la presencia del fenómeno de la Niña en los años (1999-2000) -(2011-2012) donde tuvo influencia este fenómeno, teniendo un porcentaje de incidencia del 33.33% en los últimos 30 años.
Fuente: PROPIA
En la ilustración 56 mostramos los máximos caudales pico sin la presencia del fenómeno de la Niña ocurridos en los años 1994 y 2003, siendo este último año el pico más alto de los caudales generados en los 30 años de estudio, demostrándose que estas alteraciones en
los caudales no necesariamente se presentan con la presencia de fenómenos naturales.
168
3.6.8 CAUDALES MAXIMOS
Coeficinte de escorrentia
Donde:
El valor de Ka es:
Donde:
kA : Factor reductor
169