Diseño Hidrológico-Hidraúlico - Pucp
Diseño Hidrológico-Hidraúlico - Pucp
Diseño Hidrológico-Hidraúlico - Pucp
En el segundo capítulo se procede a exponer las bases teóricas sobre las que se constituye la
tesis, pasando por recapitular algunos conceptos básicos en ingeniería hidrológica, explicar la
estructura y los parámetros descriptores del modelo mensual de transformación precipitación
– escurrimiento de Témez, así como presentar los principales tópicos para el diseño hidráulico
y balance hídrico en pequeñas presas de materiales sueltos. Prosiguiendo, se explican los
factores que producen los fenómenos de erosión y depositación de sedimentos en embalses,
así como las metodologías para su cálculo (USLE). Finalmente, se presentan fundamentos del
geoprocesamiento en sistemas de información geográfica (GIS) para fines de ingeniería
hidrológica, y modelamiento numérico de filtraciones y redes de flujo en presas.
El tercer capítulo desarrolla los procedimientos necesarios en las distintas etapas para calcular
y dimensionar la presa del caso de estudio. Se comienza con una breve descripción general del
proyecto, para luego pasar a la caracterización de la cuenca Champaccocha mediante la
definición de las distintas variables en el entorno, sea geomorfometría, geología, edafología,
etc., y el análisis geofísico y estratigráfico.
5 BIBLIOGRAFÍA ................................................................................................ 97
ANEXOS
I. Apéndices auxiliares
II. Gráficos totales
III. Mapas nacionales
IV. Memorias de cálculo
Índice de Figuras
Índice de Cuadros
Índice de Tablas
Índice de Imágenes
Índice de Ecuaciones
En años recientes, diversas instituciones del estado han realizado publicaciones de mapas
oficiales de gran valor ingenieril, lo cual sugiere un avance importante al momento de
caracterizar cuencas, sin embargo, estos adelantos no están esquematizados y plasmados en
una guía de actualización de procedimientos y fuentes de información.
Por otro lado, en el ámbito actual de la ingeniería hidráulica nacional no se cuenta con un nivel
de consenso suficientemente alto sobre un modelo de simulación hidrológica continua, para
que este pueda ser dispuesto como la directriz oficial que determine el proceder en una
infinidad de proyectos de disponibilidad hídrica. En este contexto, no es posible concentrar
esfuerzos para el calibrado de los parámetros descriptores en las distintas zonas del país, que
conlleve a la creación de mapas regionales que permitan estandarizar los resultados, tales como
los que ya se han generado en otros países iberoamericanos como Uruguay o España.
Cuadro 1.1-1 Inventario de presas en el Perú por finalidad o uso (ANA, 2015)
-1-
La principal finalidad de uso de las presas en el Perú (cuadro 1.1-1) es con fines de riego (442
de 743 presas), por lo que, en la actualidad, instituciones del estado como el Proyecto
Subsectorial de Irrigaciones (PSI) o el MINAGRI, gestionan múltiples proyectos sociales de
ayuda para la tecnificación del riego de parcelas agrarias. Los cientos de proyectos de pequeñas
presas que se requieren y se desean poner en marcha para apoyar al sector agro, se encuentran
en su gran mayoría en el ámbito altoandino del país. Este escenario supone modelar la
capacidad hídrica del embalse en cabeceras de cuenca, donde no existe aforo de caudales.
Por otro lado, desde que ha cobrado relevancia el problema de sedimentación en embalses, se
han desarrollado múltiples métodos para la remoción de sedimentos (como el flushing), sin
embargo, las alternativas convencionales no alcanzan altos grados de eficiencia. En el Perú, el
Instituto de Hidráulica, Hidrología e Ingeniería Sanitaria (IHHS) de la Universidad de Piura,
ha instalado un sistema de hidrosucción en el embalse Cirato de la central hidroeléctrica de
Carhuaquero, el cual se encuentra en la actualidad en funcionamiento.
No obstante, si bien este tipo de técnicas avanzadas1 da buenos resultados, solo se justifica su
uso en grandes proyectos, por lo que su aplicación queda fuera del alcance de las pequeñas
presas, ya que el nivel de impacto que representan no es comparable al de una obra principal.
Los principales lineamientos sobre los que se desarrolla el documento son los siguientes:
1
Para más detalle en técnicas de remoción de sedimentos, se puede consultar Ref. [59].
-2-
el diseño de pequeñas presas, donde por la magnitud y menor complejidad de la estructura,
es aceptable y necesaria una estandarización del procedimiento de ingeniería a utilizar.
Finalmente, se exponen las bondades tecnológicas de varios de los softwares actuales más
potentes con aplicaciones para la ingeniería hidrológica e hidrogeológica, demostrando su
mayor exactitud y precisión con los métodos computacionales modernos.
Alcances de la investigación
Al contrario de un proyecto concebido desde sus inicios, esta tesis retoma y usa la información
de varios estudios previos a nivel de perfil realizados en años anteriores. Por lo tanto, ya que
el lugar del emplazamiento de la presa es fijo, por la información de ingeniería de detalle
disponible, ya no se ejecuta la tradicional evaluación previa de ubicación óptima de la presa.
Luego, se evaluará si es posible satisfacer la demanda agrícola proyectada para esta cerrada.
Otros procesos intermedios como el análisis de filtraciones y el cálculo del volumen muerto
de sedimentos también se desarrollan de forma conveniente.
En específico, se busca que la tesis llegue a ser considerada un precedente importante para el
establecimiento de una norma oficial para la caracterización de cuencas y el cálculo de la
disponibilidad hídrica de pequeñas presas en cuencas no aforadas en Perú. Se espera que ante
la futura publicación del referido manual en Perú, más países opten por tomar el mismo
camino, estandarizando y enriqueciendo el modelamiento hídrico, tan necesario para poder
producir a gran escala proyectos hidráulicos importantes a nivel social en el continente.
-3-
2 MARCO TEÓRICO
2.1 TÓPICOS DE INGENIERÍA DE RECURSOS HÍDRICOS
La ciencia hidrológica ha sido ampliamente estudiada y documentada. Este acápite solo se
limita a especificar ciertos aspectos puntuales sobre evapotranspiración y clases texturales de
suelos. Para temas más elementales se recomienda revisar ref. [18] y [19].
Evapotranspiración de referencia
El concepto fue desarrollado haciendo uso de la definición del cultivo de referencia como un
cultivo hipotético con una altura de 0,12 m, con una resistencia superficial de 70 s.m -1 y un
albedo de 0,23. Este cultivo ideal representa a la evapotranspiración de una superficie extensa
de pasto verde de altura uniforme, creciendo activamente y adecuadamente regado.
A medida que el cultivo se desarrolla, tanto su altura como el área foliar varían
progresivamente, por lo que existen diferencias en la ET y, por lo tanto, en el valor de Kc. Este
período de crecimiento puede ser dividido en cuatro etapas: inicial (Kc ini), mediados de
temporada (Kc med), final de temporada (Kc fin), y las de desarrollo intermedio (fig 2.1-1).
-4-
Debido a que las superficies de las cuencas altoandinas están fuertemente dominadas por
pastizales andinos y lagunas, usualmente basta solo con conocer estos coeficientes de cultivo,
para lo cual se hace uso de las tablas referenciales proporcionadas por la FAO (cuadro 2.1-1).
Cuadro 2.1-1 Coeficiente de cultivo Kc para pastos y cuerpos de agua superficial (FAO, 2006)
Por otro lado, existe un coeficiente final de ajuste (Kc aj), dado que bajo condiciones de campo
se puede disminuir el crecimiento y densidad de las plantas, reduciendo la tasa de ET. Esto
sucede principalmente debido al estrés hídrico, y a limitaciones ambientales y de prácticas de
manejo (fig. 2.1-2). Sin embargo, dado que solo con un control sistematizado de agricultura
es posible obtener ese tipo de información tan específica, se opta por descartar el uso de Kc aj,
y emplear de forma inicial la ET potencial que desarrolla el cultivo (ETc), pudiendo este valor
preliminar ser distinto de la ET real, la cual se determinará por balance hídrico en el cap. 2.2.1.
Método de Hargreaves
Si bien la FAO proporciona fórmulas simplificadas para el cálculo de las variables faltantes,
puede ser conveniente utilizar criterios basados solo en la temperatura del aire, no obstante,
estos métodos al ser empíricos requieren cuidadosas calibraciones locales para alcanzar
resultados satisfactorios. Una excepción posible es el método de Hargreaves, que ha producido
valores razonables de ETo con cierta validez global (FAO, 2006).
-5-
donde:
ETo = Evapotranspiración de referencia diaria, mm/día
t med, máx y min = Temperatura media, máxima y mínima diaria, oC
R0 = Radiación solar extraterrestre, convertida en mm/día
A diferencia del sistema SUCS, este sistema de clasificación establece grupos de suelos con
comportamiento análogo frente al agua y las plantas (Porta et al, 2003).
La textura USDA expresa al suelo con las proporciones de las partículas minerales inferiores
a 2 mm. Las combinaciones posibles en los porcentajes de arcilla, limo y arena se agrupan en
unas pocas clases de tamaño de partículas o clases texturales2, y su determinación se realiza
mediante un diagrama conocido como el triángulo textural de suelos (fig. 2.1-3).
Figura 2.1-3 Relación de los constituyentes de suelo por tamaño, definiendo las clases
texturales y subclases de arena (USDA, 1967, adaptado por FAO, 2015)
2
El calificativo “Franco” (Loam) se puede encontrar también bajo la denominación “Marga” o “Migajón”
-6-
2.2 MODELO MENSUAL DE TRANSFORMACIÓN
PRECIPITACIÓN - ESCURRIMIENTO
El sistema hidrológico
Desde la unidad de estudio que se define en la cuenca hidrológica, se precisan los componentes
y sus interacciones, modelando la realidad hacia un esquema teórico que pueda ser usado como
base para realizar procesamientos posteriores que estimen las cantidades de agua necesarias
para proyectar las grandes obras de la ingeniería hidráulica.
En base a este modelo conceptual del comportamiento del sistema hidrológico, es posible crear
un modelo de simulación continua que, mediante algoritmos especialmente diseñados, busca
representar y cuantificar los procesos físicos que ocurren en la cuenca hidrológica.
Debido a los largos períodos de evaluación de datos (que usualmente son años, sino décadas
de procesamiento de datos), y a diferencia de los modelos de eventos extremos, estos toman
en consideración los procesos lentos de transferencia hídrica de la cuenca, tales como el flujo
en acuíferos o la evapotranspiración, que pueden marcar diferencias substanciales mes a mes
con sus flujos másicos totales acumulados.
-7-
En definitiva, los modelos continuos aplican simulaciones de balance hídrico por largos
períodos de tiempo, con el fin de determinar la denominada oferta hídrica de la cuenca,
expresada en el volumen de agua procedente de los sistemas de drenaje superficial.
Por otro lado, si bien los efectos de trasvase de flujos subterráneos pueden ser determinantes
en un estudio de balance hídrico de cuencas pequeñas, se asume que los gradientes de carga
hidráulica en la superficie freática siguen la tendencia topográfica, ya que, como se ve luego
en el cap. 3.2, no hay evidencia de factores que puedan desencadenar diferencias importantes
en los niveles piezométricos esperados, tales como fallas geológicas, influencias zonales de
pozos de bombeo (fig. 2.2-2), entre otros. Por lo tanto, para el modelo conceptual se considera
que la cuenca hidrogeológica es equivalente a la cuenca hidrográfica superficial.
Se muestra dicho efecto comparando los registros de la humedad del suelo de dos simulaciones
con diferentes condiciones iniciales. Los dos histogramas mostrados (fig. 2.2-3) convergen
gradualmente y se vuelven idénticos después de algunos meses. Por lo tanto, se recomienda
excluir los resultados del primer año de simulación, dado que los contenidos de humedad (y
por tanto, la escorrentía) se ven significativamente influenciados por la condición inicial.
-8-
El modelo realiza una valoración global del área de estudio, ya que no considera la distribución
espacial de las variables y parámetros que intervienen en los cálculos, sustituyéndolos por
valores medios, por lo que se limita su aplicación a cuencas pequeñas en las que existe cierta
homogeneidad climática, edafológica y geológica (Paredes et Al, 2014).
Todo el proceso está gobernado por los principios de continuidad y de balance de masas, y
regulado por leyes específicas de reparto y transferencia entre los distintos términos del
balance (fig. 2.2-4). El modelo considera dos capas de suelo: una superior no saturada, y una
inferior saturada (acuífero), generadoras de los flujos superficiales rápido y lento,
respectivamente, las cuales se consideran como embalses (fig. 2.2-5).
Figura 2.2-4 Esquema de entradas y salidas en Figura 2.2-5 Sistema de flujos y almacenamiento
el Modelo de Témez (DINAGUA, 2011) del Modelo de Témez (Paredes et al, 2014)
Se considera que una fracción del agua que precipita sobre el terreno (P) es almacenada en la
zona superior del suelo (H), y que el resto, el excedente (T), se distribuye entre el aporte al
escurrimiento de origen superficial (Asup) y la infiltración hacia el acuífero (I) o recarga.
siendo:
-9-
y los parámetros:
De este modo, cuando la cantidad de lluvia es inferior a P0 toda el agua de lluvia se almacena
como humedad del suelo (siendo susceptible de convertirse en evapotranspiración) y el
excedente total es nulo. A partir del valor de P0, la precipitación incrementa simultáneamente
el excedente y la componente de humedad del suelo. La curva que expresa los excedentes no
es lineal, parte de 0 para P=P0 y tiende a P-𝛿 para lluvias altas (fig. 2.2-6).
Figura 2.2-6 Ley de excedentes Figura 2.2-7 Ecuación que relaciona la infiltración
(Paredes et al, 2014) con el excedente existente (Paredes et al, 2014)
Lo cual representa el cierre del balance entre el agua que queda en el suelo y la ETP. Si el
resultado es negativo la humedad del suelo al final del periodo es nula, y la evapotranspiración
potencial no se habrá desarrollado en su totalidad.
- 10 -
Siempre que exista suficiente agua en el suelo, la evapotranspiración se desarrolla hasta
alcanzar la ETP. En caso de no haberla, la humedad en el suelo al final de mes será nula.
Asimismo, la infiltración hacia el acuífero adopta una ley empírica (fig. 2.2-7), que es función
del excedente (T) y de la infiltración máxima Imáx -mm/mes- (parám. 3), que es la máxima
cantidad que puede infiltrarse en el terreno.
Para simular el comportamiento del acuífero se plantea un modelo unicelular basado en una
función exponencial que representa la curva de agotamiento de un manantial.
- 11 -
Este modelo subterráneo presenta una ventaja de simplicidad de cálculo, pero también supone
una simplificación excesiva del sistema, por lo que queda reducida su utilidad a pequeños
acuíferos. Por ejemplo, no es válido para modelar acuíferos cársticos, donde se ha identificado
más de una rama de descarga en el agotamiento del acuífero (Estrela & Sahuquillo, 1985,
citados por Paredes et al, 2014).
De esta forma, el aporte subterráneo Asub es: Finalmente, la escorrentía total AT se calcula:
Cuadro 2.2-1 Valores calibrados de los parámetros del modelo de Témez (DINAGUA, 2011)
Donde CAD es una constante propuesta para relacionar proporcionalmente Hmáx y el Agua
Disponible Total (ADT) en los suelos del Uruguay.
Cuadro 2.2-2 Rangos de los parámetros del modelo de Témez (Paredes et al, 2014)
- 12 -
Métodos y rangos específicos de estimación
Existen varios criterios para determinar este parámetro. A nivel de Latinoamérica, dos
metodologías requieren especial mención: el estudio de “Agua disponible en suelos del
Uruguay” del DINAGUA y la “Metodología de Procedimientos metodológicos de la
zonificación agroecológica y socioeconómica” del ZONISIG.
Debido a esta situación, Hmáx se estima mediante el criterio propuesto por el GIRH de la
UPV3, que realiza una correlación con el uso del suelo y la cobertura superficial (cuadro
2.2-3). Asimismo, de tener referencia de niveles freáticos superficiales, estos se deben
considerar aparte, ya que el parámetro es distinto para cuerpos de agua.
3
También hubiese sido posible hallar Hmáx mediante su correlación con texturas asociadas a clases
edafológicas (USDA-MAG, 1985, citado por Barberà, 2007)
- 13 -
Nuevamente, debido a limitaciones en la información de campo, se estima Imáx con el
criterio del GIRH de la UPV (cuadro 2.2-4), el cual propone valores promedio de
infiltración máxima para un estrato o litología específica.
El rango típico se encuentra entre 0.10 – 0.40, aunque su valor más frecuente es 0.30
(Témez, 1977, citado por Molina & Espinoza, 2005).
El ámbito de aplicación del proyecto está limitado solo a pequeñas presas de materiales sueltos,
es decir, aquellas que no puedan encuadrarse en la definición de “Gran Presa”, la cual hace
referencia a toda presa con altura mayor a 15 m., o bien la que teniendo entre 10 - 15 m. de
altura cumpla con alguna de las siguientes condiciones (ICOLD, citado por DINAGUA, 2011):
- 14 -
Componentes del volumen embalsado
Figura 2.3-1 Componentes del volumen de almacenamiento de un reservorio (Loucks & Van Beek, 2017)
Capacidad útil o activa: utilizada para la regulación del flujo aguas abajo y para el
suministro de agua, desarrollo recreativo o producción hidroeléctrica.
Volumen muerto: reservado principalmente para la recolección de sedimentos.
Adicionalmente, incluye un Volumen inactivo, dado que el punto de toma debe estar
suficientemente por debajo del almacenamiento útil para vaciar totalmente el embalse.
Almacenamiento para control de inundaciones: reservada para la laminación de
caudales durante eventos de avenidas extremas. Junto al volumen comprendido en el
borde libre, constituyen la Capacidad de contingencia, la cual controla adicionalmente
otros fenómenos como el oleaje producido por vientos dominantes, etc.
Tipología de presa
Figura 2.3-2 Ejemplo de presa de escollera con núcleo RMA (SPANCOLD, 1994)
- 15 -
Núcleo impermeable: con el objetivo de mantener la estanqueidad del represamiento,
esta capa actúa como una barrera hidráulica contra las filtraciones excesivas. Se conforma
de materiales con muy baja conductividad hidráulica, tales como las arcillas.
Espaldones: son las zonas del cuerpo de presa que soportan al núcleo. Aportan la
estabilidad y la resistencia al conjunto, por lo que definen en gran medida los taludes de
la presa. Para alcanzar la estabilidad óptima, es necesario evitar la formación de presiones
intersticiales excesivas, ya sea por niveles freáticos altos en el talud aguas abajo, como
por efectos de desembalse rápido en el talud aguas arriba, por lo que se usan materiales
granulares suficientemente drenantes, como gravas arenosas. Asimismo, para lograr
taludes más inclinados, se puede mejorar la estabilidad mecánica con el uso de geomallas.
Transiciones o filtros: este material intermedio tiene como fin asegurar la separación
entre capas diferenciadas, y, por lo tanto, reducir los problemas de tubificación interna
por el arrastre de finos. Para material de filtro se usan tradicionalmente arenas, sin
embargo, en tiempos recientes se usan geotextiles.
Drenes: con el propósito de captar y redirigir flujos internos hacia zonas adecuadas para
su control y desfogue, se diseñan drenes horizontales e inclinados, si el material de cuerpo
de presa no es lo suficientemente autodrenante. Se usan materiales granulares con buena
conductividad hidráulica como gravas bien graduadas. También es posible usar geonets.
Coraza protectora: es necesario recubrir los taludes con un material especialmente
grueso que resista la erosión producida por el oleaje, la lluvia, etc. Usualmente, se utiliza
escollera por encima de geomembrana, la cual otorga impermeabilidad extra a la presa.
Figura 2.3-3 Presa zonificada con zanja de impermeabilización total en núcleo (Novak et Al., 2001)
Una presa de materiales sueltos se puede cimentar sobre roca o suelo, sin embargo, existen
estratos especialmente desfavorables que se deben evitar, por lo que como mínimo se debe
retirar la primera capa superficial de suelo orgánico -top soil-, y las capas debajo de la zona de
actuación de la presa con capacidad portante baja y de conductividades hidráulicas altas.
- 16 -
En el escenario especial en que el cimiento sea permeable y poco profundo4 (como en el caso
de estudio), el medio más adecuado para controlar las filtraciones y asegurar que no se
produzca una falla por sifonamiento a través de los cimientos, o subpresiones en el pie del
paramento de aguas abajo, es a través de una zanja de impermeabilización que llegue hasta el
lecho rocoso u otro estrato impermeable (fig. 2.3-3). Este elemento actúa como barrera
hidráulica en la cimentación (debe ubicarse ligeramente aguas arriba del eje de la presa),
además aumenta la estabilidad del núcleo impermeable, dado que su dentellón va unido a este.
Pre-dimensionamiento
Ancho de corona
Esta medida está en función de la altura de la presa y los análisis de estabilidad, aunque se
recomienda un ancho mínimo de 3 m. (fig. 2.3-4). Existen normativas de varios países que
proporcionan fórmulas para su determinación, entre las más reconocidas se encuentra la
expresión recomendada por el Bureau of Reclamation:
La inclinación de los taludes está en función de la estabilidad de la presa, para lo cual se deben
ejecutar diversos análisis que definen los valores adecuados, sin embargo, de forma preliminar,
es posible usar valores referenciales. El Bureau of Reclamation, con base en la posibilidad de
desembalse rápido, recomienda taludes y materiales para la construcción de los terraplenes
(cuadro 2.3-1), así como la proporción geométrica del núcleo proyectado en función del grado
de permeabilidad del cimiento (fig. 2.3-4). Del mismo modo, la base mayor del núcleo
impermeable se fija de una longitud mínima igual a la altura de la presa, actuando además
como umbral superior de la zanja de impermeabilización. Por último, también se recomienda
un espesor y granulometría específica para la escollera protectora (cuadro 2.3-2).
4
En pequeñas presas este término está asociado a profundidades menores a la altura de la presa.
- 17 -
Figura 2.3-4 Proporción de núcleos impermeables en presas de tierra (Bureau of Reclamation, 2007)
Cuadro 2.3-1 Taludes y materiales recomendados para presas de tierra (Bureau of Reclamation, 2007)
Cuadro 2.3-2 Espesor y granulometría recomendados para la escollera (Bureau of Reclamation, 2007)
- 18 -
Resguardo y borde libre
Se entiende por resguardo a la distancia que existe entre el nivel superior del dique y el máximo
del embalse. El "borde libre normal" se define como la diferencia de altura entre la coronación
de la presa y el nivel normal de embalse, fijado por las necesidades del proyecto. El "borde
libre mínimo" se define como la diferencia de altura entre la coronación y el máximo nivel de
embalse que puede producirse durante la avenida de proyecto, asumiendo que el aliviadero y
los desagües funcionan conforme a lo previsto (Bureau of Reclamation, 2007).
Por lo tanto, según el criterio del Bureau of Reclamation (citado por DINAGUA, 2011):
a) Es borde libre normal para considerar el b) Es borde libre mínimo para considerar el
efecto de la ola máxima que se desarrolla efecto de la ola máxima que se desarrolla
con el embalse en la cota de inicio de en el embalse cuando se produce la lámina
vertido, por lo que: máxima de vertido (E), por lo que:
Resguardo = BL Normal Resguardo = E + BL Mínimo
Figura 2.3-5 Determinación del Fetch máximo Cuadro 2.3-3 Borde libre según Fetch (Bureau of
en espejo de agua irregular (elab. propia) Reclamation, 1987, citado por DINAGUA, 2011)
Como se menciona en el capítulo inicial (cap. 1.2), por estar fuera del alcance de la tesis no se
realiza laminación en embalses, por lo que se toma directamente el valor del BL normal.
- 19 -
2.3.2 BALANCE HÍDRICO EN RESERVORIOS
Como se sabe, la altura de una presa pequeña tiene un límite, 15 metros. Esta restricción
condiciona de manera significativa el rango de elevaciones que puede tomar la estructura, por
lo que análisis de pre-dimensionado de embalses que usan la capacidad útil que satisface el
100% de la demanda como función objetivo (p. ej. método de los picos secuenciales), no son
óptimos. En su lugar, se usa un modelo de simulación para calcular las demandas posibles a
ser atendidas, ya que permite acotar mejor los valores de las alturas evaluadas.
Figura 2.3-6 Elementos del balance hídrico en el reservorio (modificado de DINAGUA, 2011)
Aportaciones: Pérdidas:
Escorrentía de la cuenca (Vol esc) Vol. entregado por riego agrícola (Vr)
Precipitación directa en el embalse (P) Evaporación en el embalse (Ev)
Vertido por encima del aliviadero (Vv)
Salidas hídricas varias (Vs):
filtraciones, caudal ecológico, etc.
Estas ecuaciones representativas del embalse se ajustan mejor a una escala logarítmica
(DINAGUA, 2011), por lo que con este fin se calcula una cota especial (H*), utilizando las
áreas medidas (Ai) y las cotas (Hi) en el rango de niveles entre la obra de toma y el vertedero:
- 20 -
Con las alturas virtuales calculadas por las diferencias Hi - H* y sus áreas correspondientes,
se ajusta la función Aemb, calculando los coeficientes α y β por mínimos cuadrados de las
funciones logarítmicas. Luego, la función Vemb resulta de integrar Aemb.
Finalmente, las funciones de la geometría del embalse presentan las siguientes distribuciones:
Tal como se vio previamente, para lagunas someras (h < 2 m.) el coeficiente de cultivo de 1.05
es constante (cuadro 2.1-1), sin embargo, en cuerpos de agua superficial profundos (h > 5 m.)
existe una distorsión en la distribución estacional de la ETP producida, debido al efecto de
retardo térmico originado por la absorción y retención temporal de energía de radiación dentro
del embalse. Sin embargo, este ciclo de amortiguación y retraso solo se completa en grandes
embalses, dado que operan en buena parte de su tiempo en niveles por encima de los 5 m.
En el caso de embalses pequeños, las alternancias entre niveles menores y mayores a este
umbral son más frecuentes, por lo que el ciclo se rompe y se libera la radiación retenida a
penas el embalse se vuelve somero. Debido a esta particularidad y para simplificar los cálculos,
se toma el valor medio entre los coeficientes correspondientes a verano e invierno, por lo tanto,
se fija para profundidades del embalse mayores a 5 m. un Kc = 0.95. Del mismo modo, para el
rango intermedio (2<h<5 m.), se precisa el valor promedio entre los otros dos grupos (Kc = 1).
Finalmente, y aunque los valores son muy cercanos y podría tomarse uno solo sin mayores
implicancias, se programa el cálculo para tomar en consideración las diversas situaciones.
ℎ < 2 𝑚 → 1.05
𝐾𝑐 = [ 2 < ℎ < 5 𝑚 → 1
ℎ ≥ 5 𝑚 → 0.95
Una vez hallados todos los valores de los parámetros antes descritos, y definida el área
cultivada, con la cual queda determinada la demanda agrícola total, ya es posible simular la
operación del embalse para determinar los caudales ofertados por la represa.
- 21 -
A priori, no es posible determinar el volumen de aporte de la precipitación y la evaporación
dentro del embalse, ya que no se dispone de un espejo de agua representativo del mes, por lo
que, en primera instancia, se calcula un volumen auxiliar despreciando su diferencia (Evi - Pi
= 0), y suponiendo además que no existe vertido (Vvi = 0) y se satisface la demanda (Vri =
Vdi). Esto se realiza con el propósito de aproximar el espejo de agua representativo, para luego
hallar el volumen bruto en el embalse ahora sí con todas las variables involucradas y con los
reajustes respectivos, tal como se sigue a continuación:
No obstante, este volumen bruto Bi (en el esquema anterior aparece de forma general como
Vembi) debe compararse con la capacidad del embalse, y en caso sea mayor que el volumen
máximo, o menor que el mínimo, entonces se hacen las correcciones correspondientes. Si es
el caso, también se determina el volumen vertido, según las siguientes relaciones lógicas:
Finalmente, según se haya procedido en los apartados anteriores, se determina el volumen real
regado (respecto a la demanda inicial) en el mes analizado, según el siguiente algoritmo:
- 22 -
Este índice varía en base a la cota de vertido (Hv) y el área regada (Ar), por lo que para una
mejor elección de los alcances finales del proyecto es conveniente realizar un análisis de
sensibilidad alterando porcentualmente dichas variables.
Entre los procesos externos que han moldeado la geomorfología de la Tierra, el concepto más
relacionado con la erosión es el fenómeno de intemperización. Es necesario puntualizar este
tema para tener un entendimiento más amplio de la dinámica erosiva.
Figura 2.4-1 Formación del suelo a través de la intemperización (según Cooper, Soil Basics,
U. Minnesota, http://www.soils.umn.edu/academics/classes/soil2125/index.html)
La mayoría de la superficie seca del planeta está compuesta por materiales no consolidados
provenientes de rocas. Aquellas partículas sólidas alcanzaron su estado actual a través de un
proceso de interacción química, física y/o biológica denominado intemperización, que consiste
en la alteración progresiva de la estructura litológica original, dando como resultado un
material particulado (Ritter et al., 2002, citado por García-Chevesich, 2015).
El fenómeno erosivo
- 23 -
Figura 2.4-2 Proceso de erosión (Suárez, 2001)
Erosión hídrica
Erosión y sedimentación son dos conceptos que se encuentran ligados cuando la causante de
la pérdida de suelo es el agua. Existe una relación clara entre la precipitación media anual y
las tasas de sedimentación producidas (fig. 2.4-3). La dinámica erosiva inicia en los desiertos
(lluvias escasas), en donde mientras más intensa se torna la precipitación, mayor es también
el nivel de erosión. Esta relación se mantiene lineal hasta superar el umbral en el que comienza
a formarse vegetación, la cual amortigua la caída de la lluvia, con lo cual a partir de este punto
de inflexión la erosión unitaria se ve reducida con el aumento de la precipitación, ya que
produce a su vez más vegetación.
Por otro lado, típicamente, las tasas más altas de erosión hídrica se concentran en un reducido
número de eventos, representados por las grandes tormentas (alto período de retorno). No
obstante, en el largo plazo, la mayor proporción de la erosión total sucede debido a los eventos
de menor magnitud (bajo período de retorno), los que en cantidad superan el suelo perdido en
los eventos más catastróficos (Morgan, 2005).
Erosión eólica
Evidentemente, el factor más importante es la velocidad del viento, la cual se puede asumir
como cero a una altura 𝑧0 determinada en función de la rugosidad superficial producida por
obstáculos, y que a partir de allí aumenta con la altura en forma exponencial.
Si bien existen métodos relativamente confiables para medir la erosión eólica, las técnicas de
medición están mucho menos desarrolladas que las utilizadas en erosión hídrica (Morgan,
2005). No obstante, a fin de tener una idea del orden de magnitud de la velocidad de los
vientos, se puede consultar la escala de Beaufort, la cual consta de 13 rangos (cuadro 2.4-1).
- 24 -
Con el fin de evaluar el potencial eólico en el país con fines energéticos, el MINEM publicó
en 2016 un atlas eólico que muestra las velocidades del viento en el territorio nacional,
medidas a distintas alturas, desde los 10 m hasta los 150 m. (fig. 2.4-4). Para poder medir la
influencia local de la fuerza eólica, basta con clasificar dentro del correspondiente rango el
valor de la máxima velocidad media del viento esperada.
Este apartado no requiere mayor detalle, debido a que, como se demuestra más adelante, no se
desarrollan vientos importantes en la zona de estudio.
Modelación de la erosión
Existen diversas metodologías desarrolladas para distintos procesos erosivos. Entre las más
conocidas y mundialmente utilizadas se encuentran los siguientes:
Universal Soil Loss Equation (USLE) Soil Loss Estimation for Southern Africa
Revised Universal Soil Loss Equation (SLEMSA)
(RUSLE) Morgan, Morgan and Finney Model (MMF)
Water Erosion Prediction Project (WEPP) European Erosion Model (EUROSEM)
Si bien el modelo WEPP es el más reciente e implica una nueva generación de tecnología de
predicción de la erosión hídrica, superando limitaciones de variabilidad espacial y temporal
en topografía, propiedades y uso del suelo, etc.; requiere información mucho más precisa y,
sobre todo, en cantidades considerablemente mayores.
Debido a lo anterior, y dado el alcance trazado, se elige realizar la predicción de la erosión con
la Ecuación Universal de Pérdida de Suelo (USLE), desarrollada por el USDA.
- 25 -
2.4.2 MODELO USLE
La Ecuación Universal de Pérdida de Suelo (USLE)
La pérdida de suelo promedio anual por área unitaria está descrita y cuantificada por la
Universal Soil Loss Equation – USLE (Wischmeier y Smith, 1978 citado por Morgan, 2005):
𝐸 = 𝑅 x 𝐾 x 𝐿𝑆 x 𝐶 x 𝑃
donde:
E = Pérdida de suelo promedio anual (t ha-1 año-1)
R = Erosividad de la lluvia y escorrentía (MJ mm ha-1 h-1)
K = Erodabilidad del suelo (t ha h ha-1 MJ-1 mm-1)
LS = Factor de longitud y pendiente de la cuenca (adimensional)
C = Factor del manejo de la cobertura vegetal (adimensional)
P = Factor de prácticas de conservación de suelos (adimensional)
Una de las fórmulas de validez global con mayor aceptación es la propuesta por Morgan:
𝑅 = (𝐸𝐶 x 𝑃𝑎 x 𝐼30)/100
Ecuación 2.4-1 Fórmula de Morgan (Morgan, 2005)
donde:
R = Erosividad anual de la lluvia y escorrentía (MJ mm ha-1 h-1)
EC = Energía cinética de la lluvia (MJ ha-1 mm-1)
Pa = Precipitación anual promedio (mm)
𝐼30 = Intensidad de lluvia en 30 minutos, para un TR=2 años (mm h-1)
Si bien es razonable usar un tiempo de retorno (TR) que caracterice a la media poblacional
(50% de recurrencia), en la actualidad existe gran controversia en torno al TR asociado a 𝐼30.
Sin embargo, aun así, se adopta el período de retorno de 2 años, debido a su representatividad
estadística de la lluvia en el tiempo (USDA, 1997, citado por Mendoza, 2014).
Para poder obtener 𝐼30, se hace uso de las curvas IDF IILA–SENAMHI–UNI, una familia de
curvas de Intensidad–Duración–Frecuencia desarrolladas para las distintas regiones del Perú,
- 26 -
por el convenio de cooperación conformado por el Instituto Ítalo Latinoamericano -IILA-, el
Servicio Nacional de Meteorología e Hidrología -SENAMHI-, y la Universidad Nacional de
Ingeniería -UNI-5. En el cap. 3.6.1 se comenta la idoneidad de esta fórmula en la actualidad.
En específico, la materia orgánica puede mejorar casi todas las propiedades del suelo, pues
su presencia aumenta la aireación e infiltración, así como la cohesividad interparticular.
5
Para un mayor alcance sobre esta metodología revisar Ref. [20].
- 27 -
Asimismo, la capacidad de infiltración está en directa relación con la generación de
escurrimiento superficial, responsable de la erosión laminar y otros procesos de erosión
hídrica. Por último, es interesante indicar que las partículas medianas son las que más fácil
se erosionan, pues tanto las arcillas con su alta cohesión, como las arenas con su alta fricción,
resisten mejor, siendo los limos la clase textural más vulnerable (García-Chevesich, 2015).
No obstante, el nivel de procesamiento es alto, por lo que, para el caso de pequeñas cuencas,
se toma el criterio simplificado del Servicio de Investigación Agrícola (ARS) del USDA, el
cual relaciona el contenido de materia orgánica con características texturales (Stewart et al,
1975, citado por León, 2001), indicando medias estimadas de escalas ampliadas de valores
específicos del suelo, especificadas para la determinación de K (cuadro 2.4-2).
Como regla general, mientras más densa y homogénea sea la cubierta vegetal, mayor es su
efectividad en la disminución de la erosión laminar (García-Chevesich, 2015).
- 28 -
Figura 2.4-6 Ciclos estacionales de lluvia,
cobertura vegetal y erosión en un clima semi-
húmedo (Kirkby, 1980, citado por Morgan, 2005)
Por estas razones, el factor C (Cover-Management Factor) es uno de los agentes más
influyentes en la determinación de la erosión, tal como queda evidenciado en la fig. 2.4-6.
En el cuadro 2.4-3 se definen los valores del factor C para los tipos más comunes de cultivos
y coberturas (Wischmeier y Smith, 1978, citado por Portuguez, 2015).
La longitud de la pendiente (L) se define como la distancia horizontal desde el origen del
flujo superficial hasta el punto donde:
- 29 -
En la fig. 2.4-7 se ilustran situaciones hipotéticas para las cuales existen variaciones en las
tasas de erosión y depositación, a lo largo de cinco tipos de configuraciones de pendiente.
Figura 2.4-7 Variaciones en las tasas de erosión (E) Cuadro 2.4-4 Valores del Factor LS para
y depositación (D), según cambios topográficos diferentes pendientes (López Cadenas de Llano
(Terrence et al., 2002 citado por BREA, 2010) y Blanco, 1968, citados por Portuguez, 2015)
Sin embargo, los nomogramas y las ecuaciones fundamentales para determinar el factor LS
(Wischmeier & Smith 1978, citado por USDA, 1997) son complicados de implementar en
un entorno de geoprocesamiento GIS, debido principalmente al trabajo en formato ráster (ver
cap. 2.5). En la actualidad, el procedimiento se ha complejizado con el desarrollo de
algoritmos de cálculo, por lo que para simplificar el cálculo se utiliza el criterio sencillo del
Departamento de Hidráulica e Hidrología de la E.T.S.I. de Montes de la UPM (cuadro 2.4-
4), que permite asignar directamente el valor del factor LS en función de la pendiente de una
ladera (López Cadenas de Llano y Blanco, 1968, citados por Portuguez, 2015).
Este parámetro recoge el efecto que tienen las diferentes prácticas de conservación de suelos
en las pérdidas por erosión (León, 2001). Se pueden mencionar la siembra en curvas de nivel,
y el cultivo en fajas o terrazas (fig. 2.4-8). Estas prácticas afectan a la erosión modificando
el patrón de flujo, el grado o la dirección del escurrimiento superficial, y reduciendo la
cantidad y tasa de escorrentía (Renard y Foster, 1983 citados por USDA, 1997).
- 30 -
Finalmente, mediante análisis geoespacial, se procede a la superposición de los mapas
elaborados con los factores R, K, LS y C, y así finalmente se calcula la pérdida total de suelos
en la cuenca de estudio. En el cuadro 2.4-5 se muestra una escala de valores en t/ha/año para
poder desarrollar una idea general del orden de magnitud del monto total de erosión.
El coeficiente de entrega de sedimentos (CES) surge para ajustar la erosión total estimada y
compensarla por la deposición intermedia de sedimentos, para así, finalmente, estimar el
aporte de sedimentos en un punto específico de la cuenca.
𝐴𝑆 = 𝑇𝐸 ∗ 𝐶𝐸𝑆
donde:
AS = Aporte neto de Sedimentos de una cuenca [ton/ha/año]
TE = Erosión potencial bruta de la cuenca fuente, estimada por USLE [ton/ha/año]
CES = Coeficiente de entrega de sedimento [%]
El valor del CES depende de la distancia del área de análisis al cuerpo principal de agua.
Mientras más alejadas se encuentren las zonas de producción de sedimentos de los cursos de
agua, menor será el CES; asimismo, su valor es disminuido en grandes áreas de drenaje, suelos
de textura gruesa y topografía suave (Hairston, 2001, citado por Mendoza, 2014).
- 31 -
Su predicción se realiza mediante ecuaciones desarrolladas para relacionar las características
de la cuenca con la estimación de la erosión total. Entre los criterios simplificados con mayor
aceptación, la ASCE propone el uso del ábaco de la fig. 2.4-9 para evaluar el CES.
1
𝜂𝐴 = 1 −
𝑉
1 + 𝐾 ∗ (2.1 ∗ 103 ) ∗ ( )
𝐴
Ecuación 2.4-3 Fórmula de Brown (Brown, 1943, citado por Bodoque et Al., 2001)
siendo:
ηA = Eficiencia de atrape ; A = Área de la cuenca [km2]
V = Capacidad volumétrica total del embalse [MM3]
K = Coeficiente del efecto de las reglas de operación del embalse. Puede tomar valores
entre 0.046 y 1, dependiendo de la frecuencia de vaciado y el tipo de sedimento.
Para diseño se suele usar un valor característico de 0.10.
Por lo tanto, si un embalse retiene el agua por varios meses tendrá más sedimentación que otro
que en pocos días descarga, asimismo, aquel que sea pequeño en comparación con la magnitud
de la escorrentía anual dejará pasar gran parte del sedimento, mientras que uno con capacidad
relativamente significativa tendrá una alta eficiencia de atrape (Fatorelli, 2011).
- 32 -
El criterio de Miller
𝑡
𝛿𝑡 = 𝛿𝑜 + 0.4343 ∗ 𝐾 ∗ [ ∗ ln 𝑡 − 1]
𝑡−1
Ecuación 2.4-4 Fórmula de Miller (Miller, 1953, citado por Sánchez, 2014)
siendo:
δt = Peso específico de los sedimentos después de t años de compactación [kg/m3]
δo = Peso específico inicial, es decir, al final del primer año (cuadro 2.4-6) [kg/m3]
K = Factor de compactación que depende de las reglas de operación del embalse, y
del tamaño y textura del sedimento (cuadro 2.4-6)
Cabe realizar dos acotaciones: (1) el tipo de material está definido por el sistema USDA, y (2)
cuando se presentan distintos materiales, δo y K se obtienen mediante promedios ponderados.
Tal como fue comentado en el cap. 1.1.2, no es factible ejecutar procedimientos de remoción
de sedimentos de alta eficiencia en pequeñas presas. En este contexto, es fundamental
determinar la ubicación del sedimento dentro del embalse, ya que permite definir las zonas del
vaso donde se tendrá una disminución de capacidad. Esta parte del embalse que durante la vida
- 33 -
útil se colmata de sedimentos se conoce como volumen muerto, y no se tiene en cuenta en la
capacidad útil usada en la simulación de la operación del embalse.
La forma de distribución del material sólido en un embalse depende principalmente del tipo
de sedimento, las reglas de operación y las características del vaso y las avenidas. Según la
clasificación de Morris y Fan (fig. 2.4-10), se presentan cuatro esquemas típicos de depósitos
de sedimentos. En específico, los depósitos en cuña son característicos de embalses pequeños,
y están formados por materiales gruesos constituidos a pie de presa, que se van haciendo más
finos hacia aguas arriba (Morris y Fan, 1997, citados por Sánchez, 2014).
Figura 2.5-1 Concepto de capa de información Figura 2.5-2 Extensión de la fórmula de la USLE para la
geográfica dentro de un GIS (Olaya, 2014) combinación de conjunto de capas GIS (Olaya, 2014)
Un GIS representa la información mediante dos formatos, el vectorial, con formas definidas
como líneas, puntos o polígonos; y el ráster, con una grilla estructurada de celdas con valores
de la variable medida (fig. 2.5-3). Mediante la combinación de diversas capas (fig. 2.5-1) se
pueden obtener nuevas capas unificadas con nuevos parámetros – resultado, tal como en la
adecuación GIS del método USLE (fig. 2.5-2), modelo visto en el capítulo anterior.
Figura 2.5-3 Representación vectorial -línea roja- Figura 2.5-4 Procesos de acumulación según valores
vs ráster -cuadros negros- (Olaya, 2014) de curvaturas horizontal y vertical (Olaya, 2014)
- 34 -
Conceptos de geomorfometría y análisis del terreno
- 35 -
Figura 2.5-7 El área acumulada de una cuenca Figura 2.5-8 División de una cuenca en
con el Modelo D8 (Olaya, 2014) múltiples subcuencas (Olaya, 2014)
Lo anterior permite establecer relaciones entre las celdas y proceder así a un estudio regional.
A esta escala la medida principal es el área aportante, la cual se emplea como parámetro
auxiliar para la extracción de las redes de drenaje de la cuenca (fig. 2.5-7).
Esta extensión contiene un conjunto de herramientas creado para desarrollar de forma rápida
en ArcGIS los geoprocesamientos necesarios para llevar a cabo las simulaciones hidrológicas.
Usando el MDE como insumo, el geoprocesamiento del terreno recopila una serie de
subrutinas que son ejecutadas en secuencia y que describen colectivamente el patrón de
drenaje, además de conducir a la demarcación de subcuencas y cauces.
Demarcación de subcuencas
−𝑪𝒂𝒕𝒄𝒉𝒎𝒆𝒏𝒕 𝑮𝒓𝒊𝒅 𝑫𝒆𝒍𝒊𝒏𝒆𝒂𝒕𝒊𝒐𝒏 −
𝑫𝒆𝒍𝒊𝒎𝒊𝒕𝒂𝒄𝒊ó𝒏 𝒅𝒆 𝒔𝒖𝒃𝒄𝒖𝒆𝒏𝒄𝒂𝒔 Vectorización de polígonos de cuencas
=
(𝑷𝒓𝒐𝒄𝒆𝒔𝒂𝒅𝒐 𝒈𝒆𝒏𝒆𝒓𝒂𝒍 𝒅𝒆𝒍 𝑴𝑫𝑬) −𝑪𝒂𝒕𝒄𝒉𝒎𝒆𝒏𝒕 𝑷𝒐𝒍𝒚𝒈𝒐𝒏 𝑷𝒓𝒐𝒄𝒆𝒔𝒔𝒊𝒏𝒈 −
Agregación de subcuencas adyacentes
[ −𝑨𝒅𝒋𝒐𝒊𝒏𝒕 𝑪𝒂𝒕𝒄𝒉𝒎𝒆𝒏𝒕 𝑷𝒓𝒐𝒄𝒆𝒔𝒔𝒊𝒏𝒈 −
- 36 -
Finalmente, se procede al recorte de todas las capas con la máscara del perímetro de la cuenca
obtenida en el último paso. En caso la cuenca sea muy grande también es posible sectorizar en
varias subcuencas. Para el caso de pequeñas cuencas, basta con los pasos mencionados para
desarrollar el caso de estudio.
Todo este proceso es esquematizado por el constructor de modelos de ArcGIS, Model Builder,
el cual compila esta serie de rutinas en un solo algoritmo de cálculo, permitiendo la
automatización del análisis morfométrico (Collado & Navarro, 2013).
Figura 2.5-9 Red de flujo por debajo de una presa de concreto (GEO-SLOPE International Ltd., 2012)
SEEP/W es un modelo numérico que puede simular el proceso físico real del agua que fluye a
través de un medio poroso, siguiendo la formulación de la ley de Darcy, por lo que es
ampliamente utilizado en ingeniería para cuantificar las filtraciones (fig. 2.5-9) y describir su
comportamiento (redes de flujo) a través de distintos tipos de presas y sus cimientos.
Tipos de materiales
Son cuatro los materiales usados por SEEP/W en base a su comportamiento hidrogeológico:
a) Modelo completamente saturado: solo se aplica en el caso que se tenga la total certeza de
que el estrato asignado se encontrará por completo sumergido, por debajo de la superficie
freática. Básicamente, solo es necesario el dato de la conductividad hidráulica saturada.
b) Modelo saturado/no saturado: cuando existe duda respecto al régimen hidrogeológico
del estrato, o si se sabe con anticipación que se encontrará parcialmente sumergido,
entonces es necesario utilizar este modelo, en donde se debe definir una función de
conductividad hidráulica que varíe con la humedad y el grado de succión del estrato.
- 37 -
c) Modelo de interface: son elementos lineales usados para representar objetos de espesor
despreciable. En el caso de redes de flujo, se usa para modelar geosintéticos.
d) Ninguno: se usa para definir zonas en donde se elige no simular el flujo, debido a que el
estrato no contribuye a la disipación de la carga hidráulica, y, por lo tanto, su definición
no aporta a la solución, sino más bien ralentiza el procesador.
Las dos condiciones de contorno de interés son las de carga hidráulica total (altura de posición
más altura de presión), y la de carga de presión cero, en donde internamente SEEP/W asigna
una carga total correspondiente a la altura de posición del elemento. A modo ilustrativo, en la
fig. 2.5-10 se muestra la conceptualización típica que se sigue para la modelización numérica
de las filtraciones en una presa zonificada. Como se observa, en aras de definir condiciones de
frontera que propicien un modelo ligero, pero eficaz, se aprovecha la gran diferencia de las
conductividades hidráulicas entre el núcleo y su interacción con el material granular permeable
y los cimientos de roca madre, para poder idealizar a los últimos como estratos con carga de
presión cero (totalmente permeable) y flujo cero (totalmente impermeable), respectivamente.
Por otro lado, la opción de cuadrángulos & triángulos como elementos base para la
construcción del mallado es la más usual y la más optimizada. El tamaño de malla puede ser
distinto en cada estrato, y es posible definir discretizaciones focalizadas.
- 38 -
3 CASO DE ESTUDIO: PROYECTO PRESA
CHAMPACCOCHA
3.1 DESCRIPCIÓN GENERAL DEL PROYECTO
El proyecto fue desarrollado como una versión optimizada de diversos estudios previos que
fueron destinados a darle viabilidad al embalsamiento con represa en la cerrada de la laguna
Champaccocha (cap. 1.2). Por lo tanto, se procede al análisis en el emplazamiento ya elegido,
donde se cuenta con información de ensayos de campo y topografía de detalle.
Ubicación y acceso
Generalidades
- 39 -
Figura 3.1-3 Proyecto en GIS del sistema de abastecimiento de riego (elab. propia)
- 40 -
En décadas pasadas, la población lugareña construyó una pequeña presa “artesanal” de tierra
(imagen 3.1-2) para represar el volumen de la laguna, por lo que durante muchos años las
inmediaciones aguas arriba de la presa existente han recibido los aportes de sedimentos
transportados por los eventos de precipitación que han transcurrido en todo este tiempo. Esta
situación generó un estrato aluvial levemente consolidado que a simple inspección visual se
concluye que muy probablemente no constituye una capa adecuada para la cimentación de la
presa proyectada, y posiblemente tendrá que ser removido de la zona de actuación.
Figura 3.2-1 Algoritmo unificado de Arc Hydro Tools creado en Model Builder (elab. propia)
Posteriormente, con el software ArcGIS se realiza el geoprocesamiento del total del MDE
creando en Model Builder un algoritmo unificado de cálculo para análisis hidrológicos con
Arc Hydro Tools (fig. 3.2-1).
- 41 -
Para concluir, se ejecuta por separado el “Marcado interactivo del punto de salida”, y así se
delimita la cuenca Champaccocha y se obtiene la red de drenaje, además de los mapas
fisiográficos intermedios (fig. 3.2-2), mencionados en el cap. 2.5.1.
Figura 3.2-2 Mapas generados desde Arc Hydro Tools superpuestos con la cuenca Champaccocha
a) MDE original, b) MDE corregido sin depresiones, c) Direcciones del flujo, d) Acumulación del
flujo, e) Demarcación de subcuencas, y f) Agregación de subcuencas (elab. propia)
Los archivos se exportan a QGIS para su representación (fig. 3.2-3), en donde además se
genera un modelo pseudo 3D de la cuenca ya demarcada, superpuesto con una máscara de
Google Satellite para una mejor visualización del objeto de estudio (fig. 3.2-4).
Figura 3.2-3 MDE de la microcuenca superpuesto con Figura 3.2-4 Visualización pseudo 3D de la
curvas de nivel y red de drenaje (elab. propia) microcuenca en QGIS (elab. propia)
- 42 -
INFORMACIÓN
MORFOMETRÍA UND VALOR
REQUERIDA
Área Km2 2.69
Perímetro Km 8.64
CUENCA Kc (Coef. Gravelius) --- 1.49
Pend. Cuenca % 18.9%
Altitud media Cuenca msnm 4530
Long. Cauce Principal km 2.41
Cota Mayor msnm 4600
CAUCES Cota Menor msnm 4413
Pend. Cauce Principal % 7.76%
Densidad de drenaje km/km2 1.71
Curva hipsométrica
Las curvas hipsométricas representan la relación entre la altitud y el área acumulada a cada
elevación. Realizado el cálculo tabulado (ver anexo IV), se crea la curva hipsométrica de la
cuenca (gráfico 3.2-1/ fig. 3.2-5), y desde esta se calcula que la altura media es 4532 msnm.
- 43 -
3.2.2 BIO-CLIMATOLOGÍA, HIDROGRAFÍA Y GEOLOGÍA
Dado que la zona de estudio se encuentra en cabecera de cuenca, una caracterización precisa
es inmensamente más importante y esencial, ya que no existe información para calibrar.
Bio-Climatología
Según el Mapa Ecológico Del Perú 6, que establece la distribución geográfica de Zonas de
Vida basadas en el sistema de clasificación del Dr. Leslie R. Holdridge, el área de estudio
pertenece a la categoría de “Páramo muy húmedo Subalpino Subtropical”.
En esta categoría la temperatura media anual varía entre 6.9 a 4.6 Co, la precipitación total
anual promedio, entre 1,088.5 y 513.4 mm., y se ubica en el tipo de provincia de humedad
denominado PerHúmedo. La configuración topográfica está definida por áreas con laderas de
moderado a fuerte pendiente hasta presentar afloramientos rocosos. Por otro lado, los suelos
son relativamente profundos, generalmente con influencia volcánica (Andosols), y el escenario
vegetal está constituido por una abundante mezcla de gramíneas y otros pastos naturales.
Hidrografía
El otrora INRENA, siguiendo la metodología Pfafstetter para la delimitación, codificación y
clasificación de cuencas hidrográficas, elaboró el mapa actual de cuencas hidrográficas del
Perú (INRENA, 2005, citado por ANA, 2012). Según esta clasificación, el área de estudio es
parte de la red hidrográfica de la quebrada Chapccara, en la cuenca norte del río Pampas, y
según la demarcación administrativa realizada por la Autoridad Nacional del Agua (ANA),
pertenece a la Autoridad Local del Agua (ALA) de Bajo Apurímac – Pampas (fig. 3.2-6).
6
Ver Ref. [31]
- 44 -
Geología y cobertura superficial
Los resultados de la caracterización de los estratos se contrastan con fotografías tomadas en
campo en la zona de la laguna, y con imágenes satelitales en planta del área de estudio,
descargadas en alta resolución mediante el aplicativo SAS-Planet (ver anexo I). Estos dos
elementos auxiliares se utilizan para comprobar y evaluar los resultados obtenidos.
Siguiendo con el procedimiento, desde los Web Map Service (WPS) del INGEMMET,
MINAGRI y del MINAM se obtuvo acceso a diversas capas correspondientes a geología, uso
del suelo, fisiografía, entre otros (fig. 3.2-7), las cuales fueron revisadas en línea desde
ArcGIS, con el fin de seguir caracterizando de forma más específica la zona de actuación.
Finalmente, se descartaron mapas con información redundante y/o de poca utilidad, y solo se
emplearon cinco (5) mapas nacionales y un (1) mapa departamental7:
Figura 3.2-7 Mapas nacionales. 1. Cobertura vegetal (Fuente: MINAM) 2. Clase de suelo FAO
(Fuente: ONERN) 3. Ecología – zonas de vida (Fuente: ONERN) 4. Rugosidad superficial (Fuente:
MINEM) 5. Geología nacional (Fuente: INGEMMET) 6. Geología de Ayacucho (Fuente: MINEM)
El Mapa Ecológico Del Perú se usa para poder caracterizar la cuenca como la Zona de Vida
“Páramo muy húmedo Subalpino Subtropical”, que se vio en el apartado de bioclimatología.
Geología
Dado el mayor nivel de discretización del mapa (1:50,000), la cuenca presenta diversas capas
geológicas en su superficie. Un nivel de extensión amplio (fig. 3.2-8) es ideal para una mejor
comprensión de la historia geológica de la zona, que ultima en un mejor entendimiento de su
estratigrafía, sin embargo, para poder caracterizar la geología superficial de la cuenca, es
necesario discretizar el espacio de trabajo para un análisis más focalizado (fig. 3.2-9).
7
Todos los mapas se pueden revisar con mayor detalle en el anexo III.
- 45 -
Figura 3.2-8 Distribución referencial de la
geología regional en los alrededores de la Figura 3.2-9 Distribución geológica superficial
cuenca de estudio (elab. propia) en la cuenca Champaccocha (elab. propia)
En la leyenda de la carta geológica nacional 28n-I, los símbolos son clasificados respecto a su
edad y formación geológica, y presentan breves descripciones. La identificación de los estratos
se hizo sobre la imagen de la leyenda del archivo original escaneado del MINEM (por esto la
escala multicromática no es la misma que en el mapa acotado en GIS), y se presenta en el
cuadro 3.2-1 (Np-po/s está con la nomenclatura Np-ms/tbl).
Para la descripción de los estratos, dado que la memoria descriptiva de Paras no contiene un
inventario completo de todos sus componentes, fue necesario apoyarse en los sectores
adyacentes Huachocolpa y Arequipa para completar la caracterización.
A continuación se describen las capas superficiales asociadas, por orden de edad geológica:
La formación Labra es una de las cinco formaciones geológicas del Grupo Yura, tiene un
espesor medio de aprox. 800 m. y es una unidad de amplitud regional que se extiende desde
el sur de Arequipa hasta el suroeste de Ayacucho y sur de Huancavelica, encontrándose en
los bordes de la cuenca estratos volcánicos. Presenta una secuencia litológica bastante
variable, predominando las areniscas cuarcíticas y las cuarcitas (Vargas, 1970).
En general, se puede decir que a groso modo la Formación Labra es una arenisca
cuarzosa, una roca sedimentaria de tipo detrítico casi impermeable.
El volcanismo explosivo del Grupo Huachocolpa se manifiesta en varias etapas, entre ellas
las Formaciones Caudalosa y Portugueza. Se define como Fm. Caudalosa a un conjunto de
rocas volcánicas, esencialmente lavas andesíticas y flujos de brechas (Morche et al., 1996).
Mediante el diagrama TAS -Total Alkali Silica- (que permite definir una roca volcánica); y
el diagrama K20 vs. Si02 (para clasificación de Andesitas), se determina que el 40% de las
muestras analizadas son Andesitas, y el 60% restante, Andesitas Basálticas.
- 46 -
Cuadro 3.2-1 Descripción y clasificación de capas geológicas en la carta 28n-I (INGEMMET, 1994)
Se describe como Fm. Portugueza a las rocas volcánicas cuya litología está constituida por
brechas andesíticas y dacíticas en las unidades inferiores; y domos riodacíticos en las
unidades superiores (Morche et al., 1996). Su espesor oscila entre 20 - 50 m.
Nuevamente, mediante diagrama TAS y diagrama K20 vs. Si02, se determina que el 67% de
las muestras analizadas son Dacitas, mientras que las Andesitas y las Riolitas se dividen en
proporciones iguales el 33% restante. Sin embargo, la descripción litológica en la leyenda
detalla tobas riolíticas de feldespatos, aunque también contiene biotitas, minerales propios
de las andesitas, y, por lo tanto, de las dacitas, ya que estas últimas son rocas de formación
intermedia entre las dos anteriores (Bureau of Reclamation, 2007).
- 47 -
Por otro lado, dado que esta zona presenta vegetación escasa, es muy notorio el cambio que
sufre el entorno hacia un tono gris oscuro cuando se acerca a los depósitos piroclásticos de
la Fm. Portugueza, tal como se ve en las fotos de campo (imagen 3.1-1) y la imagen satelital
(imagen 3.2-1). Esto induce a considerar rocas de composición céntrica, ya que la Riolita es
la roca ígnea más clara (cuadro 3.2-2). En consecuencia, se opta por la opción intermedia
típica del miembro superior, es decir, se decide que la Formación Portugueza es Riodacita.
Son depósitos cuaternarios constituidos por los materiales arrastrados por la escorrentía y
depositados en su trayecto formando lechos, terrazas y llanuras de capas delgadas. Los
materiales aluviales están constituidos por cantos, arena y cascajo (Palacios, 1994).
Este mapa contiene información que define el tipo de suelo según el sistema de clasificación
de suelos de la FAO8. Al identificar el área de estudio, el 100% de la cuenca cae dentro de la
denominación “Leptosol dístrico - Andosol úmbrico - Afloramiento lítico”:
Leptosols: Incluye los litosoles (se ve en bioclimatología). Son suelos jóvenes muy delgados
sobre roca continua, y que son extremadamente ricos en fragmentos gruesos (pedregosos,
menos de 20% de finos). Los Leptosols conforman el Grupo de Suelos de Referencia (GSR)
más extenso sobre la tierra y son particularmente comunes en regiones montañosas.
8
Ver Ref. [34]
- 48 -
Andosols: Son suelos porosos de origen volcánico de color oscuro (se corrobora lo expuesto
en bioclimatología). Usualmente, se desarrollan en eyecciones volcánicas ricas en vidrio,
pero también se pueden desarrollar en otros materiales ricos en silicato bajo meteorización
ácida en climas húmedos y perhúmedos (FAO IUSS Working Group WRB., 2015).
Este escenario edafológico se correlaciona muy bien con el apartado geológico: estratos de
origen volcánico (Fm. Portugueza y Caudalosa), y el afloramiento rocoso (Fm. Labra) con su
correspondiente leptosol superficial. Finalmente, aunque no aparece en la edafología regional,
los depósitos aluviales pertenecen a los Fluvisols 9.
Otro de los productos que se liberaron con el Atlas eólico del Perú 10 (MINEM, 2016) es el
mapa de rugosidades, que mide la longitud de rugosidad superficial Zo. Usando la consulta
espacial en GIS, se determina que la cuenca tiene un Zo promedio de 0.11 m.
Cobertura vegetal
9
Se infiere a partir de Ref. [34], por caer en la acepción de sedimentos aluviales y lacustres.
10
Revisar Erosión eólica del cap. 2.4.
- 49 -
En esta gran unidad de cobertura vegetal se registran tres subclases: a. Césped, hierbas de hasta
de 15 cm de alto; b. Pajonal, hierbas de hasta 80 cm; y c. Tolar: arbustos de hasta 1,20 m.
(Flores et al., 2005, citado por MINAM, 2015). No obstante, por el apartado anterior, se sabe
que la cuenca tiene una cobertura de 0.11 m de altura promedio, por lo que es evidente que la
subclase específica es césped. Además, en la imagen de referencia de este ejemplar que
aparece en el documento (imagen 3.2-2) se puede verificar su enorme parecido a los pastos
andinos fotografiados en la zona de estudio.
El denominado “césped de puna”, con alturas de hasta 15 cm, está dominado por gramíneas y
graminoides (coincide con lo visto en bioclimatología).
Figura 3.2-10 Perfil típico geo-sísmico de la zona del vaso del embalse (Romero, 2016)
Se realizaron 20 líneas sísmicas de refracción con el fin de conocer de manera específica los
estratos en el sector del vaso, donde se proyecta la cimentación de la presa Champaccocha.
Para cada una de las líneas sísmicas se estableció un perfil estratigráfico con valores de
velocidad de propagación de ondas longitudinales, basándose en las dromocrónicas calculadas
de los tiempos de llegada de las ondas P.
11
Para más información sobre métodos sísmicos se recomienda revisar Ref. [14] y [24]
- 50 -
En la fig. 3.2-10 se generó un modelo geo-sísmico de tres estratos, el cual se correlaciona con
los estratos tipificados en el apartado de geología, y se subdivide en la siguiente secuencia:
Las velocidades de onda longitudinales (Vp) fluctúan entre rangos específicos dependiendo
del tipo de estrato que atraviesan, por lo que es necesario comparar estos rangos con escalas
de velocidades tipificadas (cuadro 3.2-3), y así poder deducir el tipo de estrato analizado.
Cuadro 3.2-3 Valores promedios de Vp, según la norma ASTM-D5777 (ASTM, 2011)
El primer y último estrato caen sin problema dentro de los rangos esperados de “Suelo
intemperizado” y “Arenisca”, que, como es evidente, corresponden a las capas de Depósitos
Aluviales Recientes y Arenisca cuarzosa, determinadas previamente. Sin embargo, el segundo
estrato, que previamente fue establecido como Andesita, está muy por debajo de los umbrales
inferiores de los rangos de rocas, aunque también se encuentra por encima de los valores
correspondientes a suelos granulares. Esta situación evidencia que la roca no es sana.
Continuando con el análisis, cabe señalar que por inferencia geológica, las líneas de refracción
sísmica cortan ambos miembros de la Fm. Caudalosa, produciendo resultados igualmente
dispersos en el rango antes indicado, sin embargo, debido a su explícita descripción en la
leyenda de la carta geológica, el grado de alteración del miembro inferior se asume mayor.
Por lo tanto, teniendo en cuenta los resultados geofísicos, ahora se considera el Miembro
Medio (Nm-ca_m) como Andesita fracturada y alterada, mientras que el Miembro Inferior
(Nm-ca_i) recibe la denominación de Andesita muy fracturada y alterada.
- 51 -
Análisis estratigráfico
En los perfiles estratigráficos de las calicatas elegidas (CC-2, CC-3, CC-5 y CC-12) se aprecia
que la caracterización geotécnica (SUCS y otros ensayos) se realiza a partir del segundo
estrato, dado que la capa orgánica superficial no es relevante para la ingeniería de
cimentaciones, dado su pésimo comportamiento como suelo de fundación.
Como se explica en el cap. 2.4, mientras mayor sea el grado de desarrollo del perfil, los estratos
comienzan a presentar características cada vez más diferenciadas. En el caso de las calicatas
analizadas, se observa que los perfiles no se encuentran totalmente desarrollados, pero sí
presentan un nivel de intemperización intermedio, lo que implica que si bien las propiedades
de las dos primeras capas superficiales pueden ser similares, la composición del estrato
inorgánico no es totalmente representativa de la capa superior. En este contexto, no es posible
usar la granulometría de las calicatas para representar los estratos superficiales con fines de
cálculo de la erosión (cap. 2.4) o de determinación de parámetros del modelo de Témez (cap.
2.2). Por esta razón, se optó por utilizar los métodos alternativos ya explicados en cada acápite.
12
Ver Ref. [47]
- 52 -
Figura 3.2-12 Ficha descriptiva de perfil estratigráfico típico (Municipalidad de Vilcanchos, 2011)
Sin embargo, en base a las clases de suelo descritas en el acápite edafológico, se obtienen
valores típicos de las proporciones de arena, limo y arcilla, así como el porcentaje de materia
orgánica (%MO) del estrato orgánico -top soil-, usando la tabla anexa del mapa mundial de
suelos (FAO, 2007), creado en base a la clasificación edafológica FAO del 2006, por lo que
algunos calificativos de suelos varían con respecto a la versión actual (2015). En esta versión,
los Leptosols están representados por los Litosols, y el estrato aluvial se asimila mejor al
subgrupo de los Fluvisols dístricos. En el caso del suelo correspondiente a la Fm. Caudalosa,
el subgrupo más cercano al horizonte úmbrico que aparece en esta versión, es el de los
Andosols móllicos. Por otro lado, dada su composición de roca inferior, la mayor probabilidad
está en que la capa superior de la Fm. Portugueza pertenezca al subgrupo de Andosols vítricos
(su bajo %MO explica la escasa cobertura vegetal que se observa).
Sin embargo, siendo estos datos generalizaciones de propiedades de suelo a nivel mundial,
deben validarse de alguna forma con información de campo, por lo que se contrasta la clase
textural determinada con la deducida desde las calicatas, dado que se asume que las dos capas
superficiales son por lo menos similares, esperando que el estrato vegetal presente una
composición ligeramente más fina que el inorgánico.
Por consiguiente, se realiza la comparación de las clases texturales en cada capa principal
(tabla 3.2-3), excepto en la Fm. Portugueza donde no hay calicata. Finalmente, se observa que
en todos los casos la “textura – generalización FAO” presenta una composición similar, pero
marcadamente más fina que su par, tal cual había sido previsto, dado que esta describe la capa
más superficial (la “textura – calicatas” es del estrato inorgánico).
- 53 -
Textura USDA
# Calicatas Generalización FAO
1 Franco limosa Franco
3 Franco Franco arcillosa
4 Franco arenosa Franco arcillosa
5 Areno franca Franco arcillo - arenosa
Tabla 3.2-3 Clases texturales por zona principal desde las calicatas y tabla FAO (elab. propia)
Finalmente, habiendo validado los datos extraídos desde la tabla FAO, se presentan las
características finales por capa principal (tabla 3.2-4).
PARÁMETROS COORDENADAS
ESTACIÓN PRECIPIT. TEMPERAT. TEMPERAT. TEMPERAT. ELEV.
DEPARTAM. PROV. DISTRITO LONG. LAT.
MEDIA MÁX MIN MED (msnm)
Paras 1965 - 1973
NA NA NA AYACUCHO CANGALLO PARAS 74° 38'1 13° 33'1 3,330
(156203) 1975 - 1981
Huancapi 1965 - 1982 1994 - 1999 VICTOR
1994 - 2014 1994 - 2014 AYACUCHO HUANCAPI 74° 04'14 13° 45'1 3,120
(000665) 1994 - 2014 2002 - 2014 FAJARDO
1965 - 1974 1970 - 1976
La Quinua 1965 - 1979 1966 - 1979
1977 - 1979 1998 - 2005 AYACUCHO HUAMANGA QUINUA 74° 08'7 13° 02'2 3,316
(000663) 1988 - 2014 1988 - 2014
1988 - 2014 2011 - 2014
Los Libertadores
1983 - 1991 1983 - 1991 1983 - 1991 1983 - 1991 HUANCAVELICA HUAYTARA PILPICHACA 74° 55'1 13° 24'1 4,200
(000706)
1967 - 1982
Túnel Cero
1984 - 1990 NA NA NA HUANCAVELICA HUAYTARA PILPICHACA 75° 5'10 13° 15'34 4,512
(000647)
1993 - 2015
- 54 -
Figura 3.3-1 Ubicación geográfica de las 5 estaciones y la cuenca Champaccocha (elab. propia)
Para un mejor análisis de la climatología regional se hace uso de los mapas nacionales 13 de
variación espacial de Precipitación Promedio Anual -Pp-, y de Temperatura Promedio
Máxima -Tmáx- y Mínima -Tmin- (SENAMHI, 2010). Los valores absolutos de estos mapas
(fig. 3.3-2) no son muy precisos, sin embargo, marcan tendencias de rangos regionales.
Figura 3.3-2 Mapas de variación espacial de Pp, Tmáx y Tmin en el período 1965 - 2006 (SENAMHI, 2010)
Es complejo el análisis para la elección de las estaciones a utilizar para representar el clima de
la zona de estudio. Usualmente, los patrones de temperatura se mantienen en grandes
extensiones, por lo que se usan las tres estaciones disponibles para extrapolar la información.
Al contrario, los patrones de precipitación pueden variar considerablemente entre lugares
relativamente poco distantes, por lo que solo se usan las estaciones que coinciden con la cuenca
de estudio en los mismos rangos de Pp, Tmáx y Tmin (fig. 3.3-3).
13
Mapas a escala ampliada en anexo III
- 55 -
Naturalmente, este análisis a priori debe ser verificado luego con coeficientes de correlación
representativos (mayores a 0.80) en gráficas de la variable analizada vs altitud.
Debido a este análisis, aunque la estación Túnel Cero tiene la elevación más parecida, además
del período de datos más extenso, esta no puede representar directamente la climatología del
proyecto, ya que si bien cae dentro del mismo rango de T min y Pp, no coincide en Tmáx. Por lo
tanto, se escogen las estaciones Los Libertadores y Paras para caracterizar la precipitación de
la zona de estudio, las cuales sí convergen dentro de los tres rangos visualizados.
- 56 -
3.3.2 ANÁLISIS EXPLORATORIO DE DATOS (AED)
El AED consiste en el análisis de los datos empleando una variedad de técnicas para lograr un
mayor entendimiento sobre ellos y las relaciones existentes entre las variables analizadas
(CAZALAC, 2015). Para el procesamiento de los datos en esta sección se usan principalmente
los softwares Hydrognomon y CHAC.
Evaluación preliminar
Un ciclo completo de fluctuaciones hidrológicas dura un período de tiempo que se denomina
año hidrológico, y su aplicación es fundamental para lograr la mejor independencia estadística
entre los valores anuales, por lo que se usa como estándar al momento de procesar las series
de tiempo. En Perú, comienza el 1 de Septiembre y termina el 31 de Agosto del año siguiente.
Figura 3.3-4 Cronograma inicial de Pp de las cinco estaciones analizadas (elab. propia)
Por un lado, es útil analizar los histogramas superpuestos de todas las estaciones (fig. 3.3-5)
para obtener una idea general de la distribución mensual de todas las variables involucradas.
Asimismo, la dispersión de los datos se estudia mediante su representación en gráficos de caja
-Box Plots- (realizados en el software estadístico SPSS), en los cuales se visualiza de forma
preliminar algunos outliers (datos fuera de rango) que pueden significar errores de medición
en la muestra (fig. 3.3-6).
- 57 -
Como se observa en ambas técnicas gráficas14, los valores de temperatura se mantienen con
una variabilidad mínima, al contrario de la precipitación, que presenta diferencias importantes.
En consecuencia, en adelante solo se comprueba la validez de los datos de precipitación.
14
Para revisar los registros completos de cronogramas, histogramas y box plots, ver anexo II.
- 58 -
Las siguientes etapas se realizan con Hydrognomon, el cual ofrece para la presentación de
series de tiempo un muy útil conmutador de vista lista/calendario. Introducidos los datos en la
vista calendario (fig. 3.3-7) son ordenados por año gregoriano o hidrólogico, y además se
muestran directamente los descriptores estadísticos (media, la desviación estándar, etc.).
Figura 3.3-7 Vista de calendario por año hidrológico con descriptores estadísticos (elab. propia)
Análisis de consistencia
Por lo tanto, en base a estas valoraciones, se elige a La Quinua como estación patrón.
- 59 -
ESTACIÓN CONDICIONES EN CAMPO VALORACIÓN
Probabilidad
Estación ubicada en mitad de valle en zona
Huancapi aceptable de
interior alejada de poblados importantes.
buen control
Mayor
Estacion a media ladera ubicada en la
La Quinua Probabilidad de
periferia de la ciudad de Huamanga.
buen control
Control
Estación ubicada en zona alta con
Túnel Cero especialmente
condiciones agresivas (heladas).
difícil
a) Left line slope l1=1.078, right line slope l2=.723, reduction factor m=l1/l2=1.492
1965/09
26,000
24,000
22,000
20,000
18,000
1995/09
16,000
14,000
12,000
10,000
8,000
6,000
4,000
2013/09
2,000
2,000 4,000 6,000 8,000 10,000 12,000 14,000 16,000 18,000 20,000 22,000 24,000 26,000
2,000 4,000 6,000 8,000 10,000 12,000 14,000 16,000 18,000 20,000 22,000 24,000 26,000
Sin embargo, para que la muestra sea homogénea, la CDM debe presentar una sola pendiente,
por lo que se debe elegir una de las dos partes para que quede intacta y sirva como referencia
para corregir el período contiguo. Debido a que los períodos más antiguos tienen mayor
probabilidad de estar influenciados por malas condiciones de medición, ante una
- 60 -
diferenciación evidente de períodos, y ante la falta de información histórica de las condiciones
de mantenimiento y control de la estación, solo queda asumir para restaurar la homogeneidad
que la información más reciente es la más confiable.
En el ajuste gráfico realizado se aprecia la corrección de pendiente con un factor igual a 1.492
(1.078/0.723, por las pendientes de las rectas ajustadas del período más reciente y más antiguo,
respectivamente). Finalmente, Hydrognomon automáticamente extiende la corrección para los
registros a escala mensual y presenta la CDM homogeneizada (gráfico 3.3-1b).
La CDM de La Quinua vs Túnel Cero presenta similares condiciones y se puede revisar, junto
a otra información relevante de este apartado, en el anexo II.
Asimismo, el módulo también detecta automáticamente outliers, pero en ninguno de los casos
se mostraron notificaciones al respecto. Por esta razón, y ya que se ha comprobado la
consistencia regional de las series, ya no es necesario ejecutar una prueba de datos extremos15.
CORMUL crea una ecuación de regresión para cada par de estaciones correlacionadas con la
estación a completar. Para elegir la ecuación que otorga el dato con mejor performance, se
genera una matriz de priorización, la cual es función de los coeficientes de correlación múltiple
y del número de datos comunes entre las estaciones. La influencia de este último se ajusta e
incrementa de forma proporcional con el exponente de priorización.
Por consiguiente, se elige el par de estaciones con el coeficiente de priorización más alto, y
este va decreciendo de forma gradual hasta llegar al par de estaciones que contengan el mes
que se quiere completar. Este descenso es controlable a través del umbral de priorización, el
cual establece un límite inferior, lo cual permite controlar la acumulación de errores, ya que
15
El análisis de datos atípicos mediante el módulo Range Check, se puede revisar en el anexo I
- 61 -
un mayor umbral de coeficiente de priorización implica una mejor calidad de datos predichos,
sin embargo, en contraparte, restringe la cantidad de datos completados (Álvarez, 2011).
Valores usuales en la bibliografía fijan al primero en 0.10, y admiten un umbral hasta de 0.70
(Álvarez, 2011), por lo que se usan estos valores para ejecutar la simulación.
Figura 3.3-8 Series de tiempo de Pp luego del completado por curva de doble masa (elab. propia)
Por consiguiente, la regresión se ejecuta (fig. 3.3-9) bajo el criterio de estacionalidad (la
muestra se subdivide en 12 series, por cada mes), y se rellenan los pocos meses faltantes de la
estación Huancapi. Luego, en CHAC se termina la completación total de las series de tiempo.
Por otro lado, en el caso de las temperaturas, no es posible en este momento aplicar el
algoritmo CORMUL, dado que en ningún período se cuenta con dos estaciones de referencia.
- 62 -
Dada la naturaleza de esta variable, es aceptable rellenar series de datos mediante regresión
lineal simple, sin embargo, existen casos en que se detectan lapsos de tiempo sin datos en todas
las estaciones, especialmente en la variable de Tmin (fig. 3.3-11).
Debido a la situación anteriormente comentada, solo queda recurrir a una técnica de rellenado
que no requiera estaciones de referencia, lo cual restringe bastante su alcance.
Por lo tanto, se ejecuta la predicción con las opciones anteriores, pero sin restringir el uso de
valores negativos dado que es temperatura, e incluyendo además un término aleatorio como
recurso para contrarrestar la implícita subestimación de la desviación estándar (fig. 3.3-10).
280
260
240
220
200
180
160
mm
140
120
100
80
60
40
20
0
dic 1965
dic 1967
dic 1969
dic 1971
dic 1973
dic 1975
dic 1977
dic 1979
dic 1981
dic 1983
dic 1985
dic 1987
dic 1989
dic 1991
dic 1993
dic 1995
dic 1997
dic 1999
dic 2001
dic 2003
dic 2005
dic 2007
dic 2009
dic 2011
dic 2013
9.5
9
8.5
8
7.5
7
6.5
6
5.5
°C
5
4.5
4
3.5
3
2.5
2
1.5
dic 1982
dic 1983
dic 1984
dic 1985
dic 1986
dic 1987
dic 1988
dic 1989
dic 1990
dic 1991
dic 1992
dic 1993
dic 1994
dic 1995
dic 1996
dic 1997
dic 1998
dic 1999
dic 2000
dic 2001
dic 2002
dic 2003
dic 2004
dic 2005
dic 2006
dic 2007
dic 2008
dic 2009
dic 2010
dic 2011
dic 2012
dic 2013
Gráfico 3.3-2 Histogramas completados de la estación Túnel Cero de a) Pp y b) Tmin (elab. propia)
Como se observa, los datos completados siguen una tendencia lógica respecto a la serie base,
tanto en los histogramas de precipitación como de temperatura (gráfico 3.3-2).
- 63 -
Por último, para una mejor representatividad de la muestra, debido a la mayor confiabilidad
de las series de tiempo más recientes, solo se toman los últimos 30 años, más un año previo a
esta serie (el primer año tiene una utilidad especial en el cap. 3.4.3), por lo que, finalmente,
los 31 años transcurrirán desde el año hidrológico 1982/83 hasta el 2013/14.
Con todos los datos rellenados, Hydrognomon calcula directamente el valor mensual promedio
histórico de cada estación, y luego realiza una regresión lineal para estimar el cambio con la
altitud de la variable analizada. Como resultado se presenta un gráfico de precipitación vs
altitud en donde se muestran la altitud media de la cuenca, la altitud ponderada de las
estaciones, y la recta de regresión, junto a sus estadísticos descriptores, con los cuales se
calcula el factor de corrección por altitud. Finalmente, se calcula la lluvia areal, creando la
nueva serie ponderada de la cuenca de estudio.
Tal como explicó previamente, para el caso de la precipitación se usan solo las estaciones
Paras y Los Libertadores, mientras que para las temperaturas se usan todas las estaciones
disponibles, es decir, Los Libertadores, La Quinua y Huancapi.
- 64 -
En el caso de las precipitaciones (fig. 3.3-12), los resultados son influenciados casi de forma
absoluta por la estación Paras, dada su mayor cercanía a la cuenca de estudio (cálculo de pesos
por distancia inversa), determinando un factor de corrección por altitud de 1.2. Finalmente, se
obtiene la nueva serie de precipitación mensual de la cuenca Champaccocha (gráfico 3.3-3a).
Precipitación - Champaccocha
350
300
250
200
mm
150
100
50
0
dic 1982
dic 1983
dic 1984
dic 1985
dic 1986
dic 1987
dic 1988
dic 1989
dic 1990
dic 1991
dic 1992
dic 1993
dic 1994
dic 1995
dic 1996
dic 1997
dic 1998
dic 1999
dic 2000
dic 2001
dic 2002
dic 2003
dic 2004
dic 2005
dic 2006
dic 2007
dic 2008
dic 2009
dic 2010
dic 2011
dic 2012
dic 2013
Temperatura media - Champaccocha
6.4
6.2
6
5.8
5.6
5.4
5.2
°C
5
4.8
4.6
4.4
4.2
4
3.8
dic 1982
dic 1983
dic 1984
dic 1985
dic 1986
dic 1987
dic 1988
dic 1989
dic 1990
dic 1991
dic 1992
dic 1993
dic 1994
dic 1995
dic 1996
dic 1997
dic 1998
dic 1999
dic 2000
dic 2001
dic 2002
dic 2003
dic 2004
dic 2005
dic 2006
dic 2007
dic 2008
dic 2009
dic 2010
dic 2011
dic 2012
dic 2013
Del mismo modo, se procesan las series de temperaturas (ver anexo II), hallando coeficientes
de correlación altos (|𝑟| > 0.93) en todos los escenarios. Lógicamente, al contrario de la
precipitación, los factores de corrección son menores que la unidad dado que la relación de la
altitud con la temperatura es inversamente proporcional (coeficientes r negativos).
Una vez creadas todas las series temporales de la cuenca Champaccocha, se determina que la
precipitación media anual es 1033.01 mm y que la temperatura media anual es 5.15 Co, valores
que encajan perfectamente en los rangos de Pmed (1,088.5 - 513.4 mm) y Tmed (6. 9 - 4.6 Co) de
la zona de vida “Páramo muy húmedo Subalpino Subtropical” (estimado en cap. 3.2).
- 65 -
Asimismo, es conveniente contrastar los valores promedio de las series obtenidas con estudios
zonales que los validen. En el caso de las precipitaciones, dado que presentan variaciones
espaciales importantes, no es adecuada una comparación con valores generalizados regionales,
sin embargo, en el caso de las temperaturas, sí puede ser un indicador útil.
Como se aprecia, las diferencias son despreciables, comprobando la validez del cálculo.
Tmáx (Co) ENE. FEB MAR. ABR. MAY. JUN. JUL. AGO. SET. OCT. NOV. DIC. ANUAL
ANA 11.30 11.30 12.30 12.80 13.00 13.00 12.30 12.50 12.30 12.30 12.80 12.80 12.39
Calculado 14.17 13.68 13.72 13.93 14.20 13.80 13.89 14.40 14.92 15.57 15.87 15.09 14.44
%Var. (Ko) -1.0% -0.8% -0.5% -0.4% -0.4% -0.3% -0.5% -0.7% -0.9% -1.1% -1.1% -0.8% -0.7%
Tmin (Co) ENE. FEB MAR. ABR. MAY. JUN. JUL. AGO. SET. OCT. NOV. DIC. ANUAL
ANA -0.30 0.00 -0.30 -1.00 -2.80 -4.80 -5.00 -4.30 -2.50 -1.80 -1.30 -0.80 -2.07
Calculado -3.51 -3.44 -3.42 -4.08 -5.89 -7.08 -7.57 -6.80 -5.33 -4.48 -4.28 -3.89 -4.98
%Var. (Ko) 1.2% 1.2% 1.1% 1.1% 1.1% 0.8% 0.9% 0.9% 1.0% 1.0% 1.1% 1.1% 1.1%
Tabla 3.3-2 Comparación de resultados de temperatura con valores regionales (elab. propia)
Tomando como base las descripciones precisadas por la FAO16, se infiere que el cultivo
definido previamente como césped de puna se adapta mejor a la definición de Heno -Rye
Grass-, y las lagunas al tipo Agua libre somera en climas sub-húmedos o trópicos, detallados
en el cuadro 2.1-1, desde el cual son tomados los valores de Kc. Ya que este tipo de pasto
presenta valores muy cercanos durante su ciclo de vida, se usa el valor promedio de Kc = 1.
Este coeficiente resulta consistente con la realidad, ya que sus características y altura (0.11 m)
son prácticamente las mismas que las del cultivo de referencia.
16
Ver Ref. [25]
- 66 -
En el caso de los cuerpos de agua, aparte de la laguna Champaccocha, existe una adicional en
el perímetro de la cuenca (fig. 3.4-1). Estas lagunas solo liberan agua mediante evaporación,
y dado que son poco profundas su Kc es 1.05. Sin embargo, no representan ni el 1% de la
superficie total, por lo que dadas estas condiciones se fija un único Kc = 1 para toda la cuenca.
Evapotranspiration (Hargreaves)
125
120
115
110
105
mm
100
95
90
85
80
75
dic 1982
dic 1983
dic 1984
dic 1985
dic 1986
dic 1987
dic 1988
dic 1989
dic 1990
dic 1991
dic 1992
dic 1993
dic 1994
dic 1995
dic 1996
dic 1997
dic 1998
dic 1999
dic 2000
dic 2001
dic 2002
dic 2003
dic 2004
dic 2005
dic 2006
dic 2007
dic 2008
dic 2009
dic 2010
dic 2011
dic 2012
dic 2013
Ya que esta sección se desarrolla con la fórmula de Hargreaves (cap. 2.1), se utilizan las series
temporales de Tmed, Tmáx y Tmin, calculadas en el capítulo anterior. Esta operación se ejecuta
mediante el módulo de Evapotranspiración de Hydrognomon (fig. 3.4-2), en el cual, una vez
incorporadas las series mencionadas, se introduce la latitud del centroide de la cuenca de
estudio y se indica que para los cálculos de determinación de la radiación solar extraterrestre
el día representativo se fije a mitad del mes.
- 67 -
Finalmente, se obtiene la serie temporal de ETo (gráfico 3.4-1), que es igual a la ETP (ya que
Kc es 1), desde la cual se calcula un valor promedio anual de 1207.29 mm.
Tal como se realizó con las temperaturas, se comparan los resultados con los valores de la
cuenca del Alto Pampas (ANA, 2010). Si bien la evapotranspiración potencial es una variable
más compleja y no es tan constante a escala regional como la temperatura, la variación a escala
anual entre ambos registros no debería ser tan significativa. Como se aprecia en la tabla 3.4-1,
la variación anual calculada respecto del promedio regional es aceptable (menos del 10%), por
lo que se confirma la validez de los cálculos.
ETP (mm) ENE. FEB MAR. ABR. MAY. JUN. JUL. AGO. SET. OCT. NOV. DIC. ANUAL
ANA 99.10 95.10 92.90 89.10 82.70 79.70 85.20 85.20 91.00 105.50 105.90 106.40 1117.8
Calculado 113.93 100.56 104.57 93.16 88.20 78.95 84.27 96.21 103.98 117.74 119.18 118.66 1219.4
%Var. (Ko) -15.0% -5.7% -12.6% -4.6% -6.7% 0.9% 1.1% -12.9% -14.3% -11.6% -12.5% -11.5% -9.1%
Tabla 3.4-1 Comparación de resultados de ETP con valores regionales (elab. propia)
Cuerpos de agua superficial: las lagunas son medidas en base a la imagen satelital,
y caen dentro de capas aluviales y volcánicas (fig. 3.4-1).
Espacios con poca vegetación: debido a lo observado en las fotografías de campo
e imágenes satelitales, es evidente que la zona de la Fm. Portugueza presenta
cobertura vegetal escasa debido a su rango superior de elevación, dejando notar la
coloración oscura de su capa formadora.
Praderas y pastizales naturales: como se vio previamente, la vegetación se definió
como césped de puna, por lo que el resto de la cuenca se asigna a este grupo.
Por consiguiente, con esta clasificación, y haciendo uso del cuadro 2.2-3, se obtiene un
valor promedio de Hmáx = 147 mm (tabla 3.4-2), el cual se encuentra dentro del rango
típico de 50 – 250 mm (cuadro 2.2-2), por lo que se valida el resultado obtenido.
- 68 -
b) Capacidad máxima de infiltración, Imáx (mm/mes)
Por último, haciendo uso del cuadro 2.2-4, se definen los valores Imáx (tabla 3.4-3).
El resultado cae dentro del rango de 100 – 400 mm/mes, especificado por Témez, por
lo que se acepta el valor de Imáx = 335 mm/mes.
Se calcula con la fórmula de Lutz Scholz para cuencas de puna (ec. 2.2-1):
Según lo indicado en el cap. 2.2.2 por Témez, para propósitos de diseño, en la práctica
se considera que C = 0.30.
17
Se puede revisar el original en https://mvotma.gub.uy/component/k2/item/10003003-manual-de-
diseno-y-construccion-de-pequenas-presas
- 69 -
dependiente de las condiciones iniciales (ver cap. 2.2), se le excluye de la serie histórica del
escurrimiento (tabla 3.4-4). Finalmente, con los resultados restantes, se determina que el
volumen de escurrimiento anual promedio de la cuenca es 0.85 Hm3 (Qprom = 27.5 l/s).
Gráfico 3.4-2 Comparación entre precipitación y escurrimiento a nivel mensual (elab. propia)
- 70 -
Gráfico 3.4-3 Proporción del escurrimiento por épocas del año (elab. propia)
La determinación del volumen óptimo del embalse es un proceso iterativo (cap. 2.3.2), por lo
que inicialmente se comienza el análisis para la altura máxima de una pequeña presa: 15 m.
desde su cimentación. Posteriormente, se analiza la sensibilidad del modelo con la altura.
Luego, se procede al modelado BIM -Building Information Model- de la represa con Civil 3D.
Antes de iniciar la conformación de la presa, se deben retirar las capas que no cumplen con las
condiciones de ingeniería propicias para servir como estratos de cimentación de la presa.
- 71 -
De acuerdo a los perfiles correlacionados (tabla 3.2-4), y tomando en cuenta los lineamientos
del cap. 2.3.1, en el primer estrato se encuentran depósitos aluviales recientes, los cuales al
estar semi-consolidados, presentan malas propiedades geomecánicas e hidrogeológicas,
además de albergar en su superficie capas de suelo orgánico, por lo cual, dada su relativa poca
profundidad, esta capa será removida en su totalidad de la zona de actuación de la presa.
Por otro lado, el siguiente estrato de roca fracturada se considera adecuado con fines de
capacidad portante, sin embargo, aun permite filtraciones moderadas hacia el exterior de la
presa, por lo que se instalará una zanja de impermeabilización total que penetrará en la capa
hasta encontrar la arenisca y asegurar la estanqueidad del embalse.
Con motivo de demarcar el área de actuación, se realiza el trazado (fig. 3.5-1) del eje de la
futura represa (L = 258 m.), en donde mediante una breve inspección visual se puede encontrar
intuitivamente la ubicación de la presa artesanal actual (PK 180 - 200), mencionada en el cap.
3.1. En este lugar en especial se ha presentado el mayor espesor de la capa aluvial respecto de
las demás regiones del alineamiento, debido a que es la zona de fondo de valle por donde se
concentra el escurrimiento, transportando y depositando (por el efecto de la presa artesanal)
mayores cantidades de sedimentos, razón por la cual se realiza una distinción entre los sectores
existentes, y se les separa entre “zona de fondo” al emplazamiento de la presa actual, y “zona
de ladera” a todas las demás regiones. En base a los perfiles geofísicos, la distribución de las
profundidades de los estratos según el tipo de zona se resume en la tabla 3.5-1.
Por otro lado, como es evidente, en la zona de fondo se encuentra el mayor desnivel
topográfico entre ambos lados del eje, además de la menor cota a la que podría proyectarse la
nueva presa. Precisamente, la altura de una represa se considera desde la cota más baja de la
cimentación, por lo que este sitio en específico es crítico para el diseño de ingeniería.
Mediante el uso del software Autocad Civil 3D e importando como data reference la superficie
topográfica usada en el cap. 3.2, se realiza un corte transversal al eje de la actual presa artesanal
(fig. 3.5-2) para tener referencia de los niveles que condicionarán el pre-dimensionamiento.
Como se observa, la cota más baja en la zona de actuación es poco mayor a 4410 msnm, por
lo que considerando además los 2 m. de profundidad que se remueven del estrato aluvial por
los motivos antes expuestos, entonces la cota de cimentación de la presa se define en 4408
- 72 -
msnm. Asimismo, asumida la máxima altura que puede tener una pequeña presa (15 m.), la
cota de coronamiento queda definida en 4423 msnm.
Por otro lado, realizando el proceso iterativo mediante el cual se aproximan los derrames de
la presa para definir la sección de fondo, se estima una base de 30 m. a cada lado, con un talud
de corte hacia la superficie de 3:1.
Las secciones típicas, así como los perfiles longitudinales de los cortes de la excavación previa
y el alineamiento longitudinal de la futura represa, superpuestos con el terreno natural actual,
se observan en la fig. 3.5-3 y fig. 3.5-4, respectivamente. Una vez creado el TIN de la represa,
se calcula un volumen de obra construida de aprox. 30,000 m3. Finalmente, la topografía con
la presa en la zona de influencia del embalse se muestra en vista 3D (fig. 3.5-5).
- 73 -
Figura 3.5-4 Perfiles longitudinales de excavación y corona de represa en Civil 3D (elab. propia)
Figura 3.5-5 Visualización 3D de terreno con la presa en su zona de influencia (elab. propia)
Para determinar las dimensiones del embalse, se limita la topografía del vaso hasta la cota de
coronamiento (4423 msnm), por el lado donde se produce el almacenamiento. Luego, se
realiza una clasificación de elevaciones cada 1 m (fig. 3.5-6).
Figura 3.5-6 Clasificación por elevaciones del vaso de la represa (elab. propia)
- 74 -
A partir de lo anterior, se obtiene una tabla con las áreas parciales entre curvas de nivel (ver
anexo IV), con la cual obtener las curvas cota-área-volumen, y así finalmente definir la curva
de capacidad del embalse en Millones de Métros Cúbicos -MM3- (gráfico 3.5-1).
Sin embargo, la curva hace referencia al volumen total del embalse, por lo que para obtener el
almacenamiento útil es necesario descontar el volumen muerto y el volumen de contingencias.
El primero se desarrolla en el siguiente capítulo, mientras que el segundo se determina
mediante el criterio del Bureau of Reclamation descrito en el cap. 2.3.1.
Así entonces, dado que el espejo de agua en la cota de coronamiento no pasa de los 400 m., se
determina que para un Fetch máximo menor de 1.6 km., según el cuadro 2.3-3 corresponde un
resguardo de 1.20 m, lo que posiciona la cota de inicio de vertido en los 4421.8 msnm,
formando un espejo de agua máximo en el embalse de 6.64 Ha. Esto último define hasta estas
instancias la capacidad máxima del embalse en 0.293 MM3 (gráfico 3.5-2).
Gráfico 3.5-2 Determinación del volumen del embalse por debajo del borde libre (elab. propia)
- 75 -
3.6 EVALUACIÓN DE SEDIMENTOS EN EL EMBALSE
Erosión eólica
Previamente, se presentó el mapa eólico del Perú para velocidades medias anuales del viento
a distintas alturas (fig. 2.4-4), sin embargo, los valores de fuerza eólica más desfavorables se
presentan en el mes de Julio. Por lo tanto, mediante el uso del aplicativo web-GIS del MINEM
(fig. 3.6-1), se determina que la máxima velocidad media mensual erosiva del viento, oscila
entre los 4.25 m/s (h=10 m) y los 7.13 m/s (h=150 m). Según la escala de Beufort (cuadro 2.4-
1), dichas velocidades califican como viento leve y moderado, lo que expresa que no hay
potencial erosivo significativo, por lo tanto, no es necesario el cálculo de erosión eólica.
Figura 3.6-1 Vista del aplicativo que muestra el aumento de la velocidad media del
viento con la altura (http://mapas.minem.gob.pe/map_eolico/, MINEM, 2016)
La fórmula de Morgan está compuesta por la Pmedia anual, la I30 y la energía cinética (EC). La
primera se obtuvo en el cap. 3.3.4 (Pmedia = 1033.01 mm). Por otro lado, para determinar la
intensidad de lluvia en 30 minutos, se hizo uso de las curvas IDF IILA-SENAMHI-UNI (fig.
3.6-2), con lo cual se determina que la I30 para un TR = 2 años de la zona de estudio es 18.1
mm/h. Los cálculos realizados para obtener la curva IDF se presentan en el anexo I.
Tal como se explica con mayor detalle en el cap. 3.9, el SENAMHI concluye que para la sierra
central (zona del proyecto) no se han presentado variaciones significativas en los últimos 40
años, ni se proyectan grandes cambios con miras hacia el año 2030. Por lo tanto, al menos para
el caso del área de estudio, el método IILA-SENAMHI-UNI está totalmente vigente en la
- 76 -
actualidad y para las proyecciones futuras, en donde incluso se puede tomar como un criterio
conservador, dado que sobreestimaría las intensidades máximas.
Figura 3.6-2 Curvas IDF Champaccocha por método IILA-SENAMHI-UNI (elab. propia)
El parámetro restante, la energía cinética, se calcula con la fórmula de Van Dijk (ec. 2.4-2):
Factor K
Haciendo uso de la tabla 3.2-4, y mediante el criterio del ARS (cuadro 2.4-2), se determina la
erodabilidad K en base a la clase textural USDA y el porcentaje de materia orgánica, del estrato
superficial de cada zona principal (tabla 3.6-1). Cabe resaltar que se consideran aparte las
superficies de agua (donde K=0, ya que por obvias razones no hay erosión sobre agua).
Por otro lado, aunque de forma exacta el espejo de agua de la laguna Champaccocha debe ser
reemplazado por el del embalse (que es variable), esta modificación no tiene mayor relevancia,
por lo que se decide mantener la definición inicial de los cuerpos de agua superficial.
Tabla 3.6-1 Resumen del cálculo del factor K de cada estrato superficial (elab. propia)
- 77 -
Factor C
Este valor se estima mediante el cuadro 2.4-3, y se distribuye de forma muy similar al
parámetro Hmáx, del apartado de simulación del modelo de Témez. Se define a la extensión
sobre la capa geológica No. 2 (Fm. Portugueza) como territorio con cubierta escasa, y al resto
del área de suelo como pradera herbácea (nuevamente, para cuerpos de agua superficial C=0).
Cultivo y Práctica C
Cubierta escasa (60%) 0.10
Pradera herbácea en buenas condiciones 0.01
Embalse y lagunas 0
Factor LS
Figura 3.6-3 Mapa de pendientes por Figura 3.6-4 Mapas parciales USLE a) Factor R,
rangos en la cuenca (elab. propia) b) Factor K, c) Factor C, y d) Factor LS (elab. propia)
Una vez definidos todos los parámetros, se crean los mapas de los cuatro factores (no es
necesario un mapa del factor P, ya que se había definido previamente como 1, dado que es una
cuenca natural), los cuales se muestran en la fig. 3.6-4.
- 78 -
De este modo, ya que todos los valores son menores a 10 t ha-1 año-1, según la escala de
clasificación FAO-PNUMA-UNESCO (cuadro 2.4-5), la erosión producida en la cuenca es de
nula a ligera. Asimismo, se obtiene el valor total de la pérdida de suelo acumulando los
productos parciales de las pérdidas de suelo, multiplicadas por sus áreas en hectáreas.
Finalmente, la erosión bruta de la cuenca se calcula en 66.96 t/año.
Tal como se definió en el cap. 3.5.2, el volumen restante del embalse sin el borde libre queda
en 0.293 MM3. Por lo tanto, la proporción neta de los sedimentos que retiene el embalse,
representada en la eficiencia de atrape, es calculada mediante la fórmula de Brown (ec. 2.4-3):
1
𝜂𝐴 = 1 − = 0.958
0.293
1 + 0.1 ∗ (2.1 ∗ 103 ) ∗ ( )
2.69
Se considera la vida útil de la presa en 50 años (cap. 2.4.3), y una operación del embalse con
descensos de normales a moderados. Luego, se caracteriza el sedimento arrastrado, para lo
cual se ponderan los porcentajes de arena, limo y arcilla en los cinco estratos superficiales de
- 79 -
la cuenca, obtenidos desde la tabla 3.2-4. Posteriormente, con estos porcentajes promedio se
ponderan los parámetros de la fórmula de Miller, y, haciendo uso del cuadro 2.4-6, se obtiene
que δo = 1241 y K = 41 (tabla 3.6-3).
50 𝑘𝑔 𝑡
𝛿𝑡 = 1241 + 0.4343 ∗ 41 ∗ [ ∗ ln 50 − 1] = 1294 = 1.294
50−1 𝑚3 𝑚3
𝑊 34.6 ∗ 50
𝑉= = = 1337 𝑚3
𝛿 1.294
Gráfico 3.6-1 Determinación de cota teórica de volumen muerto por colmatación (elab. propia)
La gráfica marca una cota situada aprox. en 4412.4 msnm, sin embargo, como se indicó en el
cap. 2.3.1, el volumen muerto también incluye un pequeño volumen inactivo necesario para
que el punto de toma tenga una carga hidráulica mínima que garantice la correcta operación
del embalse. Esta altura es mutuamente dependiente del máximo caudal descargado, y la
sección del conducto de la toma. Ya que el diseño de estructuras de toma no es objeto de
- 80 -
estudio, se estima apropiado establecer el nivel superior del volumen muerto total en la cota
4415 msnm., y el inicio del punto de toma en la cota 4413 msnm.
Gráfico 3.6-2 Curva de capacidad efectiva o útil del embalse (elab. propia)
Como es evidente, debido al efecto del embalse se consideran las capas geológicas totalmente
sumergidas, actuando como acuíferos, mientras que la presa funcionará parcialmente saturada.
Se define el modelo para la máxima capacidad del embalse, según los lineamientos definidos
en el cuadro 2.3-1 (materiales constituyentes) y la fig. 2.5-10 (modelamiento numérico).
Núcleo impermeable: se usa como material constituyente Arcilla de baja plasticidad (CL).
Espaldones: en el caso del cuerpo de la presa, se elige Arena Limosa (SM).
Coraza protectora: dado que el enrocado no aporta a la disipación de la carga hidráulica,
y debido a su reducido espesor en comparación al cuerpo de presa (cuadro 2.3-2), se decide
con propósitos de reducir la complejidad numérica no incluir al enrocado en el modelo. Sin
embargo, se usarán elementos de interface para simular la geomembrana del talud aguas
arriba, dada su relevancia como barrera hidráulica.
18
Ver anexo II para revisar las curvas efectivas cota-área-volumen
- 81 -
Drenes: se considera que el material de cuerpo de presa no es lo suficientemente
autodrenante, por lo que se definen geonets -chimenea y horizontal-, los cuales serán
modelados como sumideros, con una condición de borde especial de presión cero. Debido
al uso de drenes, el espaldón del talud aguas abajo ya no se modela (no contribuye a la
disipación de energía), por lo que se le asigna un material tipo nulo (ver fig. 2.5-10).
Filtros: ya que se usan geotextiles para controlar la erosión interna, y estos no aportan a la
disipación de la energía, tampoco se modelarán.
Teniendo caracterizadas las capas modeladas (estratos geológicos y materiales que componen
la presa), con el uso del ábaco de la fig. 3.7-1 se procede a estimar a groso modo la
conductividad hidráulica (tabla 3.7-1). En el caso del elemento de interface, se considera de
forma conservadora para una geomembrana un K = 10-11 m/s, aunque este material puede
llegar a tomar valores tan bajos como 10-14 m/s (Koerner, 2005).
-5
Espaldones Cuerpo 10 25%
Dépositos Aluviales -4 35%
10
Capas
- 82 -
No se considera plano preferencial de flujo (no hay restricciones por estratos impermeables
inclinados), ni anisotropía (sí existe, pero si se especifica un menor valor para K y, el caudal
sería menor, por lo que para modelar de forma conservadora se establece que Kx = Ky).
Importadas las dimensiones se proceden a definir los materiales en la presa (fig. 3.7-2). Las
capas geológicas adoptan modelos totalmente saturados, por lo que solo son necesarias la
porosidad y la conductividad hidráulica saturada de cada estrato, sin embargo, en contraste,
- 83 -
los materiales de la presa se definen con el modelo saturado/no saturado, por lo que se deben
definir funciones de contenido volumétrico de agua (W/C) vs succión -presión de poros
negativa-, mediante la función aproximada por tipo de material que ofrece SEEP/W, para
luego, en base a la forma de estas, y con el método de Van Genuchten, obtener las funciones
de conductividad hidráulica vs succión (fig. 3.7-3).
Siguiendo los lineamientos de la fig. 2.5-10, dado que la coronación de la presa se define en
coordenadas relativas en la cota 38 m. del modelo, se asigna una condición de contorno de
carga hidráulica total constante del NAMO de la presa de 36.8 m (1.20 m de borde libre),
aplicada al fondo del vaso y al talud aguas arriba de la presa (fig. 3.7-4). Por otro lado, se
asigna una condición de sumidero (P=0) al talud aguas abajo del núcleo, y a la base del
espaldón aguas abajo, simulando el arreglo de geonets, con el cual se asegura un drenaje total
del agua captada (se supone un diseño previo óptimo de las geonets por capacidad drenante).
La misma condición de presión cero recibe la superficie aguas abajo de la presa.
Figura 3.7-5 Distribución de la discretización por tamaño de malla en el modelo (elab. propia)
- 84 -
Resultados del análisis de filtraciones
Figura 3.7-6 Resultados de SEEP/W: líneas equipotenciales, superficie freática y caudal circundante (elab. propia)
Figura 3.7-7 Simulación en SEEP/W de los patrones de flujo dentro de la presa (elab. propia)
Por otro lado, mediante un acercamiento de pantalla (fig. 3.7-7) se observan los patrones que
siguen las filtraciones dentro de la presa, en donde se advierte que el flujo entra desde la parte
inferior del talud aguas arriba, lo cual demuestra el buen funcionamiento de la geomembrana
como cubierta impermeabilizante, obligando al flujo a bordearla para poder seguir su camino
hacia el núcleo impermeable, el cual disipa casi de forma íntegra la carga hidráulica, para
finalmente drenar el caudal hacia las geonets (línea roja). Además, a comparación de las demás
zonas de la presa, los vectores de velocidad son casi imperceptibles en el estrato de roca sana,
lo cual pone en evidencia la gran barrera hidráulica que supone esta capa geológica.
Finalmente, SEEP/W calcula un caudal circundante de 1.4 x 10-8 m3/s por metro lineal. Dado
que la presa tiene 258 m. de longitud (fig. 3.5-1), entonces el caudal de filtración total es 3.6
x 10-6 m3/s. Esto supone un volumen anual perdido por filtraciones de aprox. 113 m3, lo que
significa el 0.01% del escurrimiento anual de la cuenca, valor que es extremadamente bajo.
Esto corrobora la alta efectividad del sistema de impermeabilización empleado.
- 85 -
3.8 SIMULACIÓN DE OPERACIÓN DEL RESERVORIO
3.8.1 DETERMINACIÓN DEL CAUDAL ECOLÓGICO
Un caudal ecológico es aquel que se diseña para mantener un nivel de funcionalidad biológica
suficiente que asegure la continuidad de los procesos y las comunidades naturales existentes
tanto en la ribera como en el medio acuático (Almagro et al., 2003, citado por Alcázar, 2007).
Cuadro 3.8-1 Extracto de la metodología de para determinar caudales ecológicos (ANA, 2016)
Figura 3.8-1 Representación de la distribución de Weibull para la serie histórica total (elab. propia)
- 86 -
Debido a la pequeña magnitud de los caudales de escurrimiento generados, se trabaja en l/s.
Exportando estas cantidades a Hydrognomon, se usa su módulo Pythia para aplicar el análisis
estadístico. De forma automática la serie mensual de escurrimientos es procesada y se pueden
elegir tanto las funciones de distribuciones teóricas (Gumbel, Log-Pearson tipo III, etc.), como
los puntos generados de las distribuciones empíricas (Blom, Gringorten, etc.), con los que son
comparados. Naturalmente, se elige la distribución de Weibull en los dos grupos y se observa
gráficamente el muy buen ajuste que presenta (fig. 3.8-1).
Las gráficas previas corresponden al uso de toda la data existente histórica, sin embargo, para
cumplir el criterio de la ANA se ajusta una curva de distribución de Weibull para cada mes
del año. Para verificar la idoneidad del ajuste se ejecuta la prueba de Kolmogorov - Smirnov
a los 3 niveles de significancia estándar (1%, 5% y 10%), en cada uno de los doce grupos. En
todos los casos la hipótesis nula es aceptada, incluso en algunos meses con niveles de
significancia de ajuste alcanzados de casi 100% (tabla 3.8-1).
Una vez confirmada la validez del ajuste, se procede a obtener el valor correspondiente al 95%
de excedencia de cada mes con la opción de pronósticos -Forecasts- (tabla 3.8-2).
Tabla 3.8-1 Aplicación de prueba de Kolmogorov – Smirnov Tabla 3.8-2 Predicción del caudal al
en el mes de Febrero (elab. propia) 95% de excedencia (elab. propia)
- 87 -
Gráfico 3.8-1 Comparación de caudales ecológicos vs escorrentía (elab. propia)
En general, la demanda es una variable que se expresa con un ciclo anual que se repite durante
todo el período simulado. Para uso agrícola, este volumen depende del comportamiento del
tipo de cultivo durante el año, el área a regar y la eficiencia del sistema de riego. Dado que
este apartado corresponde a la ingeniería agronómica, su desarrollo está fuera del alcance de
esta tesis. Debido a esto, los coeficientes usados de la demanda agrícola, en mm/ha, se obtienen
desde un documento adjunto del antiguo estudio del proyecto19. Este apartado considera
19
Ver Ref. [47]
- 88 -
alternancia entre diversos cultivos (papa, choclo, alverja, etc.), y calcula sus propios ETo, Kc
y Precipitación efectiva (Pe) por métodos distintos a los aplicados en los capítulos anteriores.
Como resultado, se obtiene la demanda estacional de riego agrícola por hectárea cultivada
(gráfico 3.8-2). Para cuantificar el total de demanda hídrica agrícola, se consideran de forma
inicial 100 hectáreas cultivadas en el proyecto. Posteriormente, se analiza su sensibilidad
respecto a su variación porcentual.
Figura 3.8-2 Ajuste de funciones probabilísticas de la geometría del embalse (elab. propia)
Las relaciones de las alturas en la presa con las áreas y capacidades del embalse deben ser
representadas por funciones continuas. Aplicando las operaciones explicadas en el cap. 2.3.2,
se obtiene una función de área con un ajuste -R2- del 99% respecto de los pares discretos (fig.
3.8-2), por lo que se acepta la función, y se procede a generar las cuatro ecuaciones
representativas de la geometría del embalse, las cuales se muestran en la tabla 3.8-3.
- 89 -
Simulación del balance hídrico
Una vez obtenidas todas las componentes descritas en el cap. 2.3.1, y asumiendo un volumen
medio inicial del embalse, se ejecutan los Algoritmos de operación del embalse, descritos en
el cap. 2.3.2, descartando nuevamente el primer año de simulación (tabla 3.8-4).
Tabla 3.8-4 Extracto de planilla de la simulación del balance hídrico en el reservorio (elab. propia)
Con los resultados del balance hídrico, se compara el volumen total de riego agrícola
demandado contra el suministrado (cap. 2.3.2), calculando un Índice de Satisfacción de la
Demanda -ISD- de 83.18%. En el anexo IV se muestra el registro temporal de la variación
de los niveles y volúmenes del embalse, así como la comparación simbólica en proporción de
los déficits y superávits mediante una representación en barras. Con estos datos se elabora el
gráfico 3.8-3, donde se observa de manera esquemática los resultados de la simulación.
- 90 -
3.8.3 ANÁLISIS DE SATISFACCIÓN DE LA DEMANDA
Para una mejor evaluación del proyecto con fines de toma de decisiones para la ingeniería
final, se realiza un análisis de sensibilidad de los parámetros relevantes, variando la cota de
vertido (Hv) en la altura del borde libre (que es hasta donde llegan las curvas área-cota-
volumen procesadas del embalse), es decir, en ±1.20 m., y el área regada (Ar), en un valor
porcentual arbitrariamente adoptado de ±20%. Para este fin, mediante una macro en Excel que
automatiza los cálculos, se genera finalmente el gráfico 3.8-4.
Para la elección de la cota definitiva de inicio de vertido (que condiciona directamente la altura
de la presa), se debe definir el margen de confiabilidad que se permite el proyecto según su
nivel de importancia. En el Perú, a diferencia de otros países de la región, no existe una norma
que determine el grado mínimo de ISD para abastecimiento hídrico (p. ej. en España, para usos
agrícolas, se definen garantías mínimas de 85-90%), sin embargo, en la revista oficial de la
ANA, se recomienda asegurar una confiabilidad del volumen atendido mayor al 90%
(González-Otoya, 2016).
- 91 -
Sin embargo, se recuerda que el área total inicialmente proyectada para ser abastecida con
propósitos de riego tecnificado era de 534 Ha (ver cap. 3.1). En este escenario se obtiene un
ISD sin proyecto de 20% y con proyecto de 21.4%, valores muy lejanos de porcentajes de
eficiencia óptimos, lo cual fuerza a reducir de forma considerable la magnitud de la demanda
hídrica agrícola y así llegar a un ISD aceptable.
Como se aprecia, la influencia del caudal ecológico provoca variaciones mínimas, del orden
de 5-10%, lo que evidencia que casi no tiene dependencia del tipo de escenario analizado, por
otro lado, se observa que, mientras el escenario final presenta aumentos significativos del
orden de 30% en condiciones con proyecto, el escenario original prácticamente no se ve
influenciado por la infraestructura.
Debido a este minucioso análisis, se decide mantener las dimensiones calculadas para la altura
máxima de una pequeña presa, por lo que para 15 m. de altura, las dimensiones y niveles
finales de la represa Champaccocha son las dispuestas en la fig. 3.8-3.
- 92 -
3.9 VULNERABILIDAD ANTE EL CAMBIO CLIMÁTICO
Aunque no es usual que se le incluya en un diseño habitual de disponibilidad hídrica, a fin de
obtener una visión global de la vigencia en el tiempo de los cálculos realizados, se estima
conveniente introducir un análisis general de cambio climático en el área de estudio.
En el año 2010, como parte de los documentos preparados para el Proyecto Segunda
Comunicación Nacional del Perú a la Convención Marco de las Naciones Unidas sobre
Cambio Climático (CMNUCC), el SENAMHI publica un informe sobre la variabilidad
meteorológica en el territorio nacional en los últimos 40 años, y realiza proyecciones de
precipitaciones y temperaturas con miras hacia décadas futuras.
Con base en el análisis de las precipitaciones ocurridas en el período 1965 – 2006, se concluye
que la región de la sierra central es la que a través de los años se convierte en la menos propicia
para la ocurrencia de inundaciones por la progresiva disminución de sus lluvias extremas, así
como de la intensidad de las mismas (SENAMHI, 2010). En el caso de las temperaturas, estas
se han mantenido relativamente constantes en la sierra central.
Las lluvias extremas estimadas mostrarían un probable decrecimiento en los próximos 30 años
en gran parte del territorio. Por otro lado, los períodos secos (CDD) están incrementándose en
mayor grado que los períodos húmedos (CWD) a nivel nacional. La tendencia del CDD
muestra incrementos en la sierra central occidental entre 2 a 3 días/año.
Figura 3.9-1 Mapas de proyecciones al 2030 de las variaciones de Pp, Tmáx y Tmin (SENAMHI, 2010)
El SENAMHI ha creado mapas nacionales de las variaciones porcentuales de las tres variables
meteorológicas principales para el escenario de emisión A2 (altas emisiones) con proyecciones
al año 2030 (fig. 3.9-1), determinando para la zona de estudio incrementos de Tmáx y Tmin de
- 93 -
+0.6% y 0.7%, respectivamente, y una disminución del 9% para Pp, valores que comparados
con los estimados en otras zonas del país (con alteraciones en algunos casos de hasta 30%) no
son significativos en términos de cambio climático.
Los valores de precipitación media anual (1033 mm) y temperatura media anual (5.15 Co),
se encuentran dentro de los rangos característicos de Pmed (1,088.5 – 513.4 mm) y Tmed (6.9
– 4.6 Co) de la zona de vida “Páramo muy húmedo Subalpino Subtropical”. Asimismo, los
rangos de temperatura calculados y los valores regionales de la ANA son muy cercanos.
Estas dos situaciones validan la metodología utilizada en el procesamiento hidrológico.
Los gráficos de variables meteorológicas vs altitud (fig. 3.3-12) corroboran las buenas
elecciones de estaciones escogidas para caracterizar la climatología de la cuenca de estudio,
tomando como base los mapas nacionales de rangos climáticos de SENAMHI. En el caso
de las precipitaciones, evidentemente siendo solo dos estaciones el coeficiente de
correlación es 1, sin embargo, si se regresiona con todas las estaciones el r (coeficiente de
Pearson) está alrededor de 0.5 (valor inaceptable), demostrando la buena selección actual
de estaciones. Asimismo, en el caso de las temperaturas, se presentan coeficientes de
correlación muy altos en todos los escenarios (mayores a 0.93), lo cual demuestra la
correcta decisión de utilizar todas las estaciones disponibles.
Los resultados observados en la gráfica de volúmenes escurridos en las épocas húmeda y
seca del año (gráfico 3.4-3), determinan que se concentra casi el 80% del escurrimiento
total en solo cuatro meses, lo que manifiesta la gran variabilidad de los caudales generados
en la cuenca Champaccocha, y evidencia su gran necesidad de regulación para uniformizar
la oferta hídrica mensual.
El volumen teórico ocupado por los sedimentos retenidos en el embalse representa poco
más del 0.37% de su capacidad útil total. Si bien existe bastante variabilidad en los métodos
para cálculo de la erosión, por lo cual esta cantidad podría ser hasta el doble, aun así seguiría
siendo una tasa relativamente baja. Esta ínfima cantidad se explica principalmente por su
cubierta vegetal, la cual propicia un adecuado control de la erosión.
Las filtraciones que traspasan la presa por año significan el 0.01% del volumen de
escurrimiento anual de la cuenca, valor casi nulo que demuestra la alta efectividad del
sistema empleado mediante zanja de impermeabilización total en los cimientos.
La simulación de la operación de presa con la demanda hídrica original (534 Ha.) implica
un aumento tan solo de 1.5% en el ISD con el proyecto ejecutado, respecto a las condiciones
naturales actuales, lo que evidencia que la presa prácticamente no genera influencia alguna
- 94 -
ante una solicitación de demanda de riego agrícola tan alta. Este escenario sería
estadísticamente imposible de satisfacer con las condiciones actuales proyectadas, por lo
que una garantía hídrica proporcionada para 80 Ha. se considera adecuada.
Esto último refleja la gran restricción hídrica de proyectar la presa en una cuenca tan
pequeña (menos de 3 Km2), y de la forma tan poco favorable del vaso desarrollado por el
represamiento (muy empinado, se gana muy poco volumen por metro de curva de nivel).
Por restricciones en la información recibida, ya se tenía fija la cerrada de la presa, sin
embargo, en las inmediaciones del proyecto se pueden identificar mejores emplazamientos
potenciales para superar la oferta hídrica alcanzada con la presa Champaccocha. Por otro
lado, otra opción es recrecer la presa Champaccocha a alturas mucho mayores (por lo
menos el doble), aunque el proyecto ya pasaría al dominio de las grandes presas, y su
viabilidad se complicaría dada la naturaleza social de este tipo de iniciativas agrícolas.
Se estima que las variaciones de temperatura y de precipitación no son significativas en
comparación con los cambios proyectados en otras zonas del país, siendo Ayacucho uno
de los departamentos más estables y que menos se verán influidos por los efectos del
cambio climático. Por lo tanto, según las proyecciones descritas, se concluye que el
presente estudio tiene validez técnica hasta el año 2030.
Comentarios finales
Sobre las fuentes de información WEB-GIS
Se evidencia la alta capacidad de geoprocesamiento del terreno digital que producen las
herramientas hidrológicas de software GIS, ya sea por ArcHydroTools en ArcGIS, o
GRASS en QGIS. La determinación de direcciones de flujo, áreas acumuladas y otros
parámetros morfológicos para cuencas y cauces, son de gran ayuda para la simplificación
del trabajo y posterior interpretación del ingeniero.
SEEP/W de Geostudio es un software muy potente para el cálculo y post-visualización de
filtraciones a través de presas. El nivel de personalización que presenta permite construir
- 95 -
modelos complejos con diversos materiales y geometrías, simulando desde estratos rocosos
y suelos, hasta materiales recientes como geosintéticos.
El algoritmo CORMUL del programa CHAC ha demostrado una potencia enorme de
predicción de datos faltantes, el cual mediante sus matrices de priorización y parámetros
de ajuste completa de manera muy razonable los períodos sin información.
El gran kit de herramientas que compila Hydrognomon implica un gran alivio en la carga
de trabajo para análisis hidrológico de series de tiempo para las distintas variables
meteorológicas, lo cual reduce en gran medida el tiempo en esta parte de la ingeniería.
Si bien en la actualidad los modelos distribuidos, que calculan de forma mallada los
resultados (como el modelo Tetis de la UPV), son parte de la siguiente generación en
modelamiento hidrológico continuo, ha quedado evidenciada la gran versatilidad del
modelo de Témez en su aplicación para pequeñas cuencas, el cual calcula con relativa
simpleza y facilidad la escorrentía superficial mensual, lo que lo hace muy útil en el medio
aplicativo de cabeceras de cuenca altoandinas.
Se considera lo realizado como un primer paso para la validación del modelo de Témez.
Lo desarrollado debe ser complementado, como cualquier modelo de balance hídrico, con
calibración con aforos de caudales. Esto puede ser realizado mediante el módulo
especializado de CHAC para calibración del coeficiente de cultivo (Kc) y los parámetros
del modelo Témez, que incluso permite modelar al acuífero con dos ramas de descarga.
Sin duda, justamente los valores de referencia con los cuales determinar los parámetros del
modelo (Hmáx, C, Imax y α) son escasos, circunstancia que resta precisión al modelo, sin
embargo, se espera que de adoptarse el modelo de Témez en el territorio nacional como
método oficial para el cálculo de la disponibilidad hídrica en pequeñas presas, se
concentren esfuerzos de calibrado para poder generar una gran base de datos en el país con
la cual potenciar la aplicación y alcance de este modelo, llegando así a resultados muy
satisfactorios, tal y como ya sucede en países vecinos como Uruguay o España.
Asimismo, concluir que para el diseño de pequeñas presas es más que suficiente este
sencillo modelo, y si se decide adoptarlo como directriz, la publicación de mapas de sus
parámetros descriptores agilizaría los cálculos de disponibilidad hídrica en gran medida.
Finalmente, resaltar que si bien es posible que con posteriores estudios técnico – económicos
se determine que, según las eficiencias de abastecimiento mostradas, la represa no es rentable,
esta conclusión es solo un resultado específico de este proyecto en particular, siendo la
principal intención y legado de esta tesis la exposición y sistematización del uso en conjunto
de los modelos y herramientas informáticas empleadas para el cálculo de la disponibilidad
- 96 -
hídrica en pequeñas presas, por lo cual son perfectamente aplicables en cientos de casos de la
mismo índole, donde muchos de ellos probablemente resultarán con evaluaciones positivas
que le den la viabilidad deseada.
Por lo tanto, se espera que esta tesis sea un documento de consulta útil para futuros expedientes
técnicos de disponibilidad hídrica de pequeñas presas en zonas altoandinas, con lo cual puedan
llevarse a cabo correctos estudios de ingeniería de forma más práctica y productiva, que luego
se traduzcan en versiones corregidas y mejoradas que potencien los cientos de proyectos
hidráulicos que necesita nuestro país.
5 BIBLIOGRAFÍA
1. Aguirre, M. (1992). Análisis y aplicación de modelos matemáticos para la generación
de caudales mensuales en cuencas de la región. Cusco: UNSAAC.
2. Alcázar, J. (2007). El Método del Caudal Básico para la determinación de Caudales
de Mantenimiento. Aplicación a la Cuenca del Ebro. Lleida: Universidad de Lleida.
3. Álvarez, J. (2011). Estimación de la distribución espacial de la precipitación en zonas
montañosas mediante métodos geoestadísticos. Universidad Politécnica de Madrid.
4. Álvarez, J. (2013). Completado de datos meteorológicos: Modelo CORMUL. Madrid:
CEDEX .
5. Autoridad Nacional del Agua (ANA). (2010). Evaluación de recursos hídricos
superficiales en la cuenca del Río Pampas. Lima: MINAGRI.
6. Autoridad Nacional del Agua (ANA). (2012). Delimitación y codificación de unidades
hidrográficas del Perú, aplicando la metodología Pfafstetter. Memoria Descriptiva.
Lima: MINAGRI.
7. Autoridad Nacional del Agua (ANA). (2015). Inventario de presas en el Perú -
Primera parte. Lima: MINAGRI.
8. Autoridad Nacional del Agua (ANA). (2016). Metodología para determinar caudales
ecológicos. Resolución Jefatural No. 98-2016-ANA (págs. 3-4). Lima: Diario Oficial
El Peruano.
9. Barberà, C. (2007). Plan director de abastecimiento y saneamiento de las cuencas
internas del departamento de la Libertad - El Salvador. Barcelona: Universitat
Politècnica de Catalunya.
10. Barreda, A., & Llasat, M.-d. C. (2004). Evolución regional de la precipitación en
España en los últimos 100 años. Madrid: CEDEX.
11. Bodoque, J., & et al. (2001). Evaluación de la degradación específica en la cuenca
vertiente al embalse de Puente Alta (Segovia), mediante métodos de estimación
directos e indirectos. C&G, 3-4(15), 21-36.
12. Bureau of Reclamation. (2007). Diseño de pequeñas presas 3ra ed. Español. Madrid:
Bellisco Ediciones.
- 97 -
13. Campos Cedeno, A. (2016). Hidráulica e Hidrología para Ingeniería. Moscú:
Universidad de la Amistad de los Pueblos de Rusia.
14. Cantos Figuerola, J. (1978). Tratado de Geofísica Aplicada 2da ed. Madrid: Litoprint.
15. CAZALAC. (2015). Guía metodológica para la elaboración del mapa de zonas
áridas, semiáridas y subhúmedas secas de américa latina y el caribe. Chile:
UNESCO.
16. CEDEX. (2010). Técnicas Geofísicas. Madrid.
17. CEDEX. (2013). Cálculo Hidrometeorológico de Aportaciones y Crecidas: Manual
CHAC. Madrid.
18. Chereque, W. (1989). Hidrología para estudiantes de ingeniería civil. Lima: PUCP.
19. Chow, V. T. (1994). Hidrología Aplicada. Bogotá: McGraw-Hill.
20. Convenio de cooperación técnica IILA & SENAMHI & UNI. (1982-1983). Estudio
de la Hidrología del Perú. Lima: Publicación del Instituto Italo-Latino Americano -
SENAMHI.
21. DINAGUA. (2011). Manual de diseño y construcción de pequeñas presas 2da ed.
Volumen l: Diseño hidrológico - hidráulico. Montevideo: MVOTMA.
22. Drake, K., & Hogan, M. (2013). Watershed management guidebook. USA: Integrated
Environmental Restoration Services, Inc.
23. ESRI. (2011). Arc Hydro Geoprocessing Tools Version 2.0 - Tutorial. New York.
24. Estrada, L. (2008). Prospección sísmica. Tucumán: Universidad Nacional de
Tucumán.
25. FAO Riego y Drenaje. (2006). Evapotranspiración del cultivo: Guías para la
determinación de los requerimientos de agua de los cultivos. Roma: FAO.
26. Fattorelli, S., & Fernández, P. (2011). Diseño Hidrológico 2da ed. Zaragoza: Water
Assessment & Advisory Global Network.
27. GEO-SLOPE International Ltd. (2012). Seepage Modeling with SEEP/W: An
Engineering Methodology. Obtenido de https://www.geo-slope.com/products/seep-w
2012
28. González-Otoya, V. (2016). Los balances hídricos y los indicadores de atención de las
demandas. Agua y más - Revista de la Autoridad Nacional del Agua. Edición No. 05,
42-47.
29. Gracia, J. (1997). Manual de Ingeniería de Ríos. México: UNAM-CONAGUA.
30. IBM Corporation. (2012). Manual del usuario del sistema básico de IBM SPSS
Statistics 21. USA.
31. INRENA. (1995). Mapa Ecológico del Perú – guía explicativa 2da ed. Lima:
MINAGRI.
32. ITIA. (2010). Documentación teórica para el software Hydrognomon versión 4.
Grecia: Universidad Técnica Nacional de Atenas.
- 98 -
33. IUSS Working Group WRB. (2007). Base referencial mundial del recurso suelo.
Primera Actualización 2007. Informes sobre recursos mundiales de suelos No. 103.
Roma: FAO.
34. IUSS Working Group WRB. (2015). Base referencial mundial del recurso suelo 2014.
Actualización 2015. Informes sobre recursos mundiales de suelos No. 106. Roma:
FAO.
35. Koerner, R. M. (2005). Designing with geosynthetics - 5th ed. New Jersey, USA:
Pearson Prentice Hall.
36. Kozani, S., & Markonis, I. (2009). Hydrognomon versión 4. Manual de usuario. ITIA.
Grecia: Universidad Técnica Nacional de Atenas.
37. Land & Water Division CBL. (2007). Digital soil map of the world v.3.6. -
Generalized soil unit information on soil characteristics. Roma: FAO.
38. López, G. (2011). Introducción a la modelización de redes de flujo mediante el uso
del programa SEEP/W 2007. España: Universidad de Córdoba.
39. Martínez, R., & Otros. (2011). SEDAH: Servidor de datos para el estudio de la
alteración hidrológica. Revista Ingeniería Civil (164), pág. 78-90.
40. McEnroe, B. (2010). Guidelines for Continuous Simulation of Streamflow in Johnson
County, Kansas, with HEC-HMS. USA: Universidad de Kansas.
41. Mendoza, F. (2014). Modelación del transporte de sedimentos en cauces de alta
pendiente, aplicación al río Las Cruces. Toluca: Universidad Autónoma del Estado
de México.
42. MINAM. (2015). Mapa Nacional de Cobertura Vegetal – Memoria descriptiva. Lima.
43. Molfino, J., & Califra, A. (2001). Agua disponible de las tierras del Uruguay.
Montevideo: DINAGUA.
44. Molina, J., & Espinoza, D. (2005). Balance Hídrico de la Cuenca Alta del rio
Pilcomayo. La Paz - Bolivia: Universidad Mayor de San Andrés.
45. Morche, W., La Torre, O., De la Cruz, N., & Cerrón, F. (1996). Geología del
cuadrángulo de Huachocolpa. Boletín, Serie A: Carta Geológica Nacional, n° 63.
Lima: INGEMMET.
46. Morgan, R. P. (2005). Soil Erosion and Conservation, 3rd edition. Oxford: Blackwell
Publishing.
47. Municipalidad distrital de Vilcanchos. (2011). Perfil Técnico: “Instalación del
servicio de agua del sistema de riego Vilcanchos, Espite, San Jacinto y San Miguel”.
Ayacucho.
48. Novak, P., Moffat, A., & Nalluri, C. (2001). Estructuras hidráulicas 2da ed.
University of Newcastle, U.K.: Chapman & Hall.
49. Ocampo, L. (2012). Análisis de Vulnerabilidad de la cuenca del río Chinchiná.
Manizales: Universidad Nacional de Colombia.
50. Olaya, V. (2004). Hidrología Computacional y Modelos Digitales del Terreno.
Madrid.
- 99 -
51. Olaya, V. (2014). Sistemas de Información Geográfica. Madrid.
52. Organización Meteorológica Mundial (OMM). (2011). Guía de prácticas
climatológicas (Vols. OMM-Nº 100). Suiza.
53. Palacios, O. (1994). Geología del cuadrángulo de Santiago de Chocorvos y Paras.
Boletín, Serie A: Carta Geológica Nacional, n° 49. Lima: INGEMMET.
54. Paredes, J., Solera, A., Andreu, J., & Lerma, N. (2014). Herramienta EvalHid para la
evaluación de recursos hídricos: Manual Técnico v1.1. Valencia: Universidad
Politécnica de Valencia.
55. PHI-LAC. (2010). Procesos de erosión - sedimentación en cauces y cuencas Vol. l.
Montevideo: UNESCO.
56. Porta, J., López-Acevedo, M., & Roquero, C. (2003). Edafología para la agricultura
y el medio ambiente 3ra ed. Madrid, España: Ediciones MundiPrensa.
57. Portuguez, D. (2015). Estimación de la pérdida de suelos por erosión hídrica en la
cuenca del Río Siguas utilizando Geoinformática. Universidad Nacional Agraria La
Molina (UNALM). Revista Anales Científicos, 76 (2), 324-329.
58. Romero, W. (2016). Estudio geofísico de refracción sísmica de la presa
champaccocha. Ayacucho: SM&W Investigaciones Geológico Geotécnico SRL.
59. Sánchez, E. (2014). Análisis de métodos para la limpieza hidráulica de sedimentos
desde embalses. Quito: Universidad Central del Ecuador.
60. Sánchez, F. J. (2013). Hidrología e Hidrogeología. Universidad de Salamanca.
Obtenido de http://hidrologia.usal.es
61. SENAMHI. (2010). Escenarios Climáticos en el Perú para el año 2030. Resumen
técnico. Segunda Comunicación Nacional de Cambio Climático. Lima: MINAM.
62. SPANCOLD. (1994). Selección de presas españolas, 1973-1993. Madrid: Centro de
Publicaciones, Ministerio de Obras Públicas, Transportes y Medio Ambiente.
63. Suárez, J. (2001). Control de erosión en zonas tropicales. Bucaramanga, Colombia:
Universidad Industrial de Santander.
64. USDA. (1997). Predicting soil erosion by water: A guide to conservation planning
with the Revised Universal Soil Loss Equation (RUSLE). Tucson: Agriculture
Handbook No. 703.
65. Van Dijk, A., Bruijnzeel, L., & Rosewell, C. (2002). Rainfall intensity–kinetic energy
relationships: a critical literature appraisal. Journal of Hydrology, vol. 261, 1-23.
66. Vargas, L. (1970). Geología del cuadrángulo de Arequipa. Boletín, Serie A: Carta
Geológica Nacional, n° 24. Lima: INGEMMET.
67. Vortex FdC SL & Barlovento Renovables Latinoamérica SAC. (2016). Atlas eólico
del Perú. Estudio para aprovechamiento eólico con fines energéticos y actualización
del atlas eólico del Perú. Lima: MINEM.
68. ZONISIG. (2001). Procedimientos metodológicos de la zonificación agroecológica y
socioeconómica. La Paz: MDP.
- 100 -
ANEXOS
I. Apéndices auxiliares
II. Gráficos
III. Mapas nacionales
IV. Memorias de cálculo
ANEXO I:
APÉNDICES AUXILIARES
Este método alternativo para realizar la curva I-D-F ha sido definido por el estudio “Hidrología
del Perú”, mediante el convenio la cooperación técnica I.I.L.A – SENAMHI – UNI en 1983.
Parámetro de intensidad
El proyecto se encuentra en subzona 12310 (fig. A1-1), con Altitud media (Y) = 4532 msnm.
Champaccocha
Fig. A1-1 Ubicación del proyecto en plano n.2-c del estudio IILA-SENAMHI-UNI (elab. propia)
Parámetro de tiempo
Parámetro de duración
Parámetro de frecuencia
Luego se procede conformar la curva IDF con los parámetros hallados y se reemplazan las
variables para el caso de interés.
Cuadro. A1-1 Intensidades para distintas duraciones y tiempos de retorno (elab. propia)
Valores de los parámetros a y n que junto con K, definen las curvas de probabilidad
pluviométrica en cada punto de las subzonas.
N° TOTAL
SUB VALOR DE
ESTACIÓN DE VALOR DE a
ZONA n
ESTACIONES
1231 321-385 2 0.357 32.2
1233 384-787-805 3 0.405 a = 37.85 – 0.0083 Y
12313 244-193 2 0.432
1235 850-903 2 0.353 9.2
1236 840-913-918-958 4 0.380 11.0
654-674-679-709-713
1238 714-732-745-752 9 0.232 14.0
Y= Altitud -msnm-, Dc= Distancia a la cordillera -Km-, y Dm= Distancia al mar -Km-
TABLA A1-2
Subdivisión del territorio en zonas y subzonas pluviométricas y valores de los parámetros K’g
y ɛg que definen la distribución de probabilidades de hg en cada punto.
El USWRC (1982) define outlier como “valor extremo o dato puntual que se separa de la
tendencia general de la muestra”. Antes de realizar cualquier tratamiento estadístico a la
información, se procede a realizar un análisis de datos dudosos, para determinar aquellos datos
de la información que se alejan significativamente de la tendencia de la información restante.
La inclusión de estos datos sin un análisis previo puede conllevar a una curva de frecuencias
distorsionada con relación a la que la muestra podría indicar.
Este análisis de extremos se realiza mediante la aplicación de umbrales que marcan un límite
máximo desde el cual valores superiores se consideran outliers. Estos umbrales se determinan
considerando la muestra con distribución normal y definiendo un nivel de confiabilidad.
Fig. A2-1 Vista de lista donde se aprecia la eliminación de la gran mayoría de etiquetas
de valores sospechosos cuando se trabaja por etapa lluviosa del año (elab. propia)
Se considera para la evaluación de la muestra marcar como extremos los valores que superen
al denominado valor máximo probable (µ + 2σ), correspondiente aproximadamente al 95% de
confiabilidad. Mediante la opción Range Check se realiza el procedimiento, por el cual al lado
derecho de las celdas que superan el umbral aparece la etiqueta “sospechoso” -suspect-.
Primero se realiza la corrida con todos los meses del año, con la cual se aprecia que los pocos
valores marcados como outliers se concentran entre los meses de Enero a Marzo (tan solo un
par muestra desfase temporal, apareciendo en Diciembre o Abril), es decir, la época húmeda,
lo cual puede indicar que dado que las lluvias de esta parte del año hidrológico concentran la
mayor proporción del total anual, se pudo otorgar la condición de extremo de forma errónea.
Para corroborar esto, se realiza un proceso de agregación de los meses Enero, Febrero y Marzo,
y sobre este valor se vuelve a correr el Range Check. En las estaciones Huancapi, Paras y Los
Libertadores, que presentan menor variabilidad entre sus valores en época lluviosa, ya no
aparece ninguna etiqueta de valor sospechoso, sin embargo, en las estaciones La Quinua y
Túnel Cero, aún se detectan un par de puntos remanentes como valores extremos (fig. A2-1).
Sin embargo, como se aprecia en el gráfico superpuesto de todos los histogramas (fig. 3.3-5a),
todas las estaciones -con datos- presentan picos en los años en que el test ha detectado outliers,
por lo que aunque estadísticamente estos pocos valores sean considerados fuera rango, es
notorio que no se debe a malas mediciones u otros problemas externos, sino a condiciones
climáticas medianamente atípicas que se desarrollaron en dichos años aislados, por lo que se
decide mantener los valores para continuar con su procesamiento en el AED.
Uso de data satelital
El estudio se apoya en la descarga o consulta online de data satelital, la cual es publicada por
múltiples organismos y está a libre disposición para ser usada con fines científicos.
Para poder visualizar y/o descarga información de mapas nacionales de las distintas variables
analizadas, primero se debe obtener la dirección url del servidor WMS -Web Map Service-
que se va a consultar. Esto se puede hacer de dos formas.
Fig. A3-1 Conexión a ArcGIS Online para visualizar mapa geológico (elab. propia)
Para poder descargar imágenes georreferenciadas de altas resolución, se hace uso del SAS-
Planet (fig. A3-4), el cual es un aplicativo ruso ampliamente difundido en la actualidad, que
cuenta con conexiones a decenas de servidores públicos, tales como Google Satellite, Bing
Roads, Yahoo Maps, entre otros. Para poder descargar la imagen que se usó en la tesis, se
delimita el área de influencia de la cuenca Champaccocha, y una vez en las opciones se clickea
el servidor Google Satellite, y se designa el formato vectorial ecw para el archivo de salida.
Antes de terminar se elige el archivo georreferenciador con extensión .map en la ventana final,
y se inicia con la descarga.
2,000 4,000 6,000 8,000 10,000 12,000 14,000 16,000 18,000 20,000 22,000 24,000 26,000
2,000 4,000 6,000 8,000 10,000 12,000 14,000 16,000 18,000 20,000 22,000
20,000
18,000
16,000
14,000
12,000 2001/09
10,000
8,000
6,000
4,000
2013/09
2,000
2,000 4,000 6,000 8,000 10,000 12,000 14,000 16,000 18,000 20,000 22,000 24,000
La Quinua vs Huancapi HMG
Homogeneous line, slope l=1.078 Determination factor=.999
1965/09
30,000
28,000
26,000
24,000
22,000
20,000
18,000
16,000
14,000
12,000
10,000
8,000
6,000
4,000
2013/09
2,000
2,000 4,000 6,000 8,000 10,000 12,000 14,000 16,000 18,000 20,000 22,000 24,000 26,000
26,000
1967/09
24,000
22,000
20,000
18,000
16,000
14,000
12,000
10,000
8,000
6,000
4,000
2013/09
2,000
2,000 4,000 6,000 8,000 10,000 12,000 14,000 16,000 18,000 20,000 22,000
2,000 4,000 6,000 8,000 10,000 12,000 14,000 16,000 18,000 20,000 22,000 24,000
CONFIGURACIONES Y SERIES DE
TIEMPO FINALES POR
INTEGRACIÓN ESPACIAL
Cuadro de integración espacial de la Temperatura Máxima
-8.5
-8
-7.5
-7
-6.5
-6
-5.5
-5
-4.5
-4
-3.5
-3
-2.5
-2
12
12.5
13
13.5
14
14.5
15
15.5
16
16.5
dic 1982 dic 1982 17
dic 1983 dic 1983
dic 1984 dic 1984
dic 1985 dic 1985
dic 1986 dic 1986
dic 1987 dic 1987
dic 1988 dic 1988
dic 1989 dic 1989
dic 1990 dic 1990
dic 1991 dic 1991
dic 1992 dic 1992
dic 1993 dic 1993
dic 1994 dic 1994
dic 1995 dic 1995
dic 1996 dic 1996
dic 1997 dic 1997
dic 1998 dic 1998
dic 1999 dic 1999
dic 2000 dic 2000
dic 2001 dic 2001
dic 2002 dic 2002
dic 2003 dic 2003
dic 2004 dic 2004
dic 2005 dic 2005
dic 2010
Serie temporal final de Temperatura Máxima
dic 2010
dic 2011 dic 2011
dic 2012 dic 2012
dic 2013 dic 2013
CURVAS EFECTIVAS
COTA-ÁREA-VOLUMEN
CURVA DE PERSISTENCIA
Curvas Cota – Área – Volumen de la
capacidad activa del embalse Champaccocha
Curva Hipsométrica
Curva de Capacidad Bruta del Embalse
Parametrización de la geometría del vaso
Modelo de transformación de Témez
Simulación de operación del reservorio
Variación de niveles en el embalse –
Déficits y Superávits
Curva Hipsométrica
Cota (msnm) Area Parcial (m2) Área Total (m2) Vol (m3) Vol acum (m3)
4411 0 0 0.00 0.00
4412 42.25 42.25 21.13 21.13
4413 482.69 524.94 283.60 304.72
4414 3436.57 3961.51 2243.23 2547.95
4415 17770.44 21731.95 12846.73 15394.68
4416 4753.01 26484.96 24108.46 39503.13
4417 5809.1 32294.06 29389.51 68892.64
4418 4835.1 37129.16 34711.61 103604.25
4419 4861.17 41990.33 39559.75 143164.00
4420 6534.03 48524.36 45257.35 188421.34
4421 7172.29 55696.65 52110.51 240531.85
4422 13422.42 69119.07 62407.86 302939.71
4423 14645.61 83764.68 76441.88 379381.58
Área de cuenca Champaccocha: Ac 269 Ha
PARÁMETROS FIJOS PARA EL MODELO DE TEMEZ
(Hm3) (mm) (mm) (mm) (Hm3)/Há (Hm3) (Hm3) (Hm3) (Ha) (Hm3) (Ha) (Hm3) (Hm3) (Hm3)
0.147 5.22
1982 9 0.03 142.16 106.62 25.65 0.0004 0.026 0.000 0.155 5.39 5.31 0.156 0.156 5.42 0.000 0.026 0.035
1982 10 0.05 111.48 107.03 0.00 0.0000 0.000 0.000 0.203 6.37 5.90 0.202 0.202 6.36 0.000 0.000 0.046
1982 11 0.05 127.34 102.67 0.00 0.0002 0.000 0.000 0.255 7.36 6.86 0.256 0.256 7.36 0.000 0.000 0.054
1982 12 0.03 4.91 117.18 35.14 -0.0011 0.035 0.004 0.251 7.28 7.32 0.242 0.242 7.12 0.000 0.035 0.026
1983 1 0.07 173.83 113.93 25.50 0.0006 0.025 0.020 0.266 7.55 7.33 0.269 0.269 7.59 0.000 0.025 0.072
1983 2 0.09 150.25 99.82 55.78 0.0005 0.056 0.052 0.250 7.26 7.42 0.251 0.251 7.28 0.000 0.056 0.091
1983 3 0.08 123.91 102.57 78.51 0.0002 0.079 0.053 0.199 6.31 6.79 0.199 0.199 6.29 0.000 0.079 0.080
1983 4 0.06 106.87 88.78 99.06 0.0002 0.099 0.029 0.134 4.91 5.60 0.133 0.133 4.90 0.000 0.099 0.063
1983 5 0.04 24.47 85.08 98.36 -0.0006 0.098 0.013 0.058 2.91 3.91 0.055 0.055 2.82 0.000 0.098 0.033
1983 6 0.02 0.68 79.79 84.60 -0.0008 0.085 0.006 -0.020 0.00 1.41 -0.021 0.001 0.16 0.000 0.063 0.015
1983 7 0.01 0.16 90.99 59.49 -0.0009 0.059 0.003 -0.055 0.00 0.08 -0.055 0.001 0.16 0.000 0.004 0.007
1983 8 0.00 19.1 102.69 19.59 -0.0008 0.020 0.001 -0.017 0.00 0.08 -0.017 0.001 0.16 0.000 0.002 0.003
1983 9 0.01 96.35 109.53 25.65 -0.0001 0.026 0.000 -0.014 0.00 0.08 -0.014 0.001 0.16 0.000 0.011 0.011
1983 10 0.01 28.62 124.98 0.00 -0.0010 0.000 0.000 0.011 1.05 0.61 0.011 0.011 1.01 0.000 0.000 0.010
1983 11 0.00 22.43 126.50 0.00 -0.0010 0.000 0.000 0.015 1.25 1.13 0.014 0.014 1.19 0.000 0.000 0.004
1983 12 0.10 218.17 117.46 35.14 0.0010 0.035 0.004 0.072 3.33 2.26 0.073 0.073 3.38 0.000 0.035 0.099
1984 1 0.08 110.57 106.05 25.50 0.0000 0.025 0.020 0.112 4.40 3.89 0.111 0.111 4.38 0.000 0.025 0.083
1984 2 0.35 355.7 90.83 55.78 0.0026 0.056 0.052 0.354 9.02 6.70 0.363 0.294 8.03 0.069 0.056 0.360
1984 3 0.42 299.39 94.42 78.51 0.0020 0.079 0.053 0.581 12.28 10.16 0.586 0.294 8.03 0.292 0.079 0.424
1984 4 0.22 108.54 84.91 99.06 0.0002 0.099 0.029 0.388 9.56 8.79 0.383 0.294 8.03 0.089 0.099 0.218
1984 5 0.11 7.06 81.65 98.36 -0.0007 0.098 0.013 0.296 8.07 8.05 0.287 0.287 7.91 0.000 0.098 0.105
1984 6 0.05 11.54 71.74 84.60 -0.0006 0.085 0.006 0.246 7.19 7.55 0.240 0.240 7.08 0.000 0.085 0.044
1984 7 0.02 9.3 80.16 59.49 -0.0007 0.059 0.003 0.200 6.32 6.70 0.195 0.195 6.22 0.000 0.059 0.017
1984 8 0.01 21.92 89.65 19.59 -0.0007 0.020 0.001 0.184 5.99 6.10 0.179 0.179 5.90 0.000 0.020 0.005
1984 9 0.00 1.1 99.08 25.65 -0.0010 0.026 0.000 0.158 5.45 5.68 0.152 0.152 5.33 0.000 0.026 -0.001
1984 10 0.00 36.71 106.47 0.00 -0.0007 0.000 0.000 0.154 5.36 5.34 0.150 0.150 5.28 0.000 0.000 -0.002
1984 11 0.06 174.61 107.16 0.00 0.0007 0.000 0.000 0.206 6.44 5.86 0.209 0.209 6.50 0.000 0.000 0.059
1984 12 0.14 207.48 107.61 35.14 0.0010 0.035 0.004 0.305 8.22 7.36 0.309 0.294 8.03 0.015 0.035 0.139
1985 1 0.08 58.73 106.04 25.50 -0.0005 0.025 0.020 0.331 8.65 8.34 0.325 0.294 8.03 0.031 0.025 0.076
1985 2 0.06 125.42 90.09 55.78 0.0004 0.056 0.052 0.248 7.22 7.63 0.249 0.249 7.24 0.000 0.056 0.063
1985 3 0.06 120.52 97.12 78.51 0.0002 0.079 0.053 0.178 5.87 6.56 0.178 0.178 5.87 0.000 0.079 0.061
1985 4 0.07 130.6 85.52 99.06 0.0005 0.099 0.029 0.116 4.50 5.19 0.117 0.117 4.52 0.000 0.099 0.067
1985 5 0.04 41.93 80.51 98.36 -0.0004 0.098 0.013 0.047 2.57 3.55 0.046 0.046 2.51 0.000 0.098 0.040
1985 6 0.02 6.44 73.81 84.60 -0.0007 0.085 0.006 -0.026 0.00 1.25 -0.027 0.001 0.16 0.000 0.057 0.017
1985 7 0.01 1.87 85.72 59.49 -0.0008 0.059 0.003 -0.053 0.00 0.08 -0.053 0.001 0.16 0.000 0.005 0.008
1985 8 0.00 26.66 99.92 19.59 -0.0007 0.020 0.001 -0.017 0.00 0.08 -0.017 0.001 0.16 0.000 0.002 0.003
1985 9 0.01 88.71 105.96 25.65 -0.0002 0.026 0.000 -0.016 0.00 0.08 -0.016 0.001 0.16 0.000 0.009 0.009
1985 10 0.01 3.66 123.28 0.00 -0.0012 0.000 0.000 0.009 0.91 0.54 0.008 0.008 0.87 0.000 0.000 0.008
1985 11 0.00 27.37 120.38 0.00 -0.0009 0.000 0.000 0.012 1.07 0.97 0.011 0.011 1.02 0.000 0.000 0.003
1985 12 0.04 152.41 120.30 35.14 0.0003 0.035 0.004 0.012 1.09 1.05 0.012 0.012 1.10 0.000 0.035 0.041
1986 1 0.22 289.3 106.67 25.50 0.0018 0.025 0.020 0.188 6.09 3.59 0.192 0.192 6.16 0.000 0.025 0.225
1986 2 0.39 312.74 89.74 55.78 0.0022 0.056 0.052 0.472 10.79 8.47 0.478 0.294 8.03 0.184 0.056 0.394
1986 3 0.26 168.14 93.99 78.51 0.0007 0.079 0.053 0.425 10.11 9.07 0.423 0.294 8.03 0.129 0.079 0.261
1986 4 0.14 60.56 86.07 99.06 -0.0003 0.099 0.029 0.310 8.31 8.17 0.304 0.294 8.03 0.010 0.099 0.138
1986 5 0.07 53.28 80.16 98.36 -0.0003 0.098 0.013 0.253 7.32 7.67 0.249 0.249 7.24 0.000 0.098 0.067
1986 6 0.03 3.93 76.22 84.60 -0.0007 0.085 0.006 0.191 6.14 6.69 0.185 0.185 6.03 0.000 0.085 0.027
1986 7 0.01 4.8 77.98 59.49 -0.0007 0.059 0.003 0.138 5.00 5.52 0.133 0.133 4.90 0.000 0.059 0.010
1986 8 0.01 12.99 87.55 19.59 -0.0007 0.020 0.001 0.119 4.57 4.74 0.115 0.115 4.48 0.000 0.020 0.003
1986 9 0.01 89.8 95.80 25.65 -0.0001 0.026 0.000 0.100 4.10 4.29 0.100 0.100 4.09 0.000 0.026 0.010
1986 10 0.01 49.01 111.81 0.00 -0.0006 0.000 0.000 0.109 4.33 4.21 0.106 0.106 4.26 0.000 0.000 0.007
1986 11 0.00 4.5 109.72 0.00 -0.0011 0.000 0.000 0.110 4.35 4.31 0.106 0.106 4.24 0.000 0.000 0.000
1986 12 0.00 44.85 110.95 35.14 -0.0007 0.035 0.004 0.068 3.22 3.73 0.065 0.065 3.15 0.000 0.035 -0.001
1987 1 0.11 230.99 111.43 25.50 0.0012 0.025 0.020 0.132 4.88 4.01 0.135 0.135 4.95 0.000 0.025 0.115
1987 2 0.23 249.32 92.93 55.78 0.0016 0.056 0.052 0.258 7.41 6.18 0.262 0.262 7.48 0.000 0.056 0.235
1987 3 0.12 56.84 97.55 78.51 -0.0004 0.079 0.053 0.253 7.32 7.40 0.247 0.247 7.20 0.000 0.079 0.116
1987 4 0.06 76.13 86.53 99.06 -0.0001 0.099 0.029 0.181 5.94 6.57 0.179 0.179 5.89 0.000 0.099 0.060
1987 5 0.03 13.91 84.51 98.36 -0.0007 0.098 0.013 0.098 4.05 4.97 0.094 0.094 3.95 0.000 0.098 0.027
1987 6 0.01 8.64 76.78 84.60 -0.0007 0.085 0.006 0.017 1.37 2.66 0.015 0.015 1.27 0.000 0.085 0.012
1987 7 0.01 15.53 84.24 59.49 -0.0007 0.059 0.003 -0.041 0.00 0.64 -0.041 0.001 0.16 0.000 0.018 0.005
1987 8 0.00 14.68 101.92 19.59 -0.0009 0.020 0.001 -0.018 0.00 0.08 -0.018 0.001 0.16 0.000 0.001 0.003
1987 9 0.00 0.33 112.00 25.65 -0.0011 0.026 0.000 -0.024 0.00 0.08 -0.024 0.001 0.16 0.000 0.001 0.001
1987 10 0.00 63.02 124.90 0.00 -0.0006 0.000 0.000 0.002 0.39 0.28 0.002 0.002 0.37 0.000 0.000 0.002
1987 11 0.00 65.27 121.53 0.00 -0.0006 0.000 0.000 0.005 0.65 0.51 0.005 0.005 0.62 0.000 0.000 0.003
1987 12 0.00 2.73 129.30 35.14 -0.0013 0.035 0.004 -0.032 0.00 0.31 -0.032 0.001 0.16 0.000 0.002 0.002
1988 1 0.07 193.21 118.95 25.50 0.0007 0.025 0.020 0.029 1.87 1.02 0.029 0.029 1.89 0.000 0.025 0.074
1988 2 0.26 294.13 104.59 55.78 0.0019 0.056 0.052 0.183 5.99 3.94 0.187 0.187 6.06 0.000 0.056 0.266
1988 3 0.14 88.52 99.93 78.51 -0.0001 0.079 0.053 0.197 6.26 6.16 0.193 0.193 6.18 0.000 0.079 0.138
1988 4 0.07 61.75 85.38 99.06 -0.0002 0.099 0.029 0.139 5.04 5.61 0.136 0.136 4.97 0.000 0.099 0.072
1988 5 0.03 34.88 81.57 98.36 -0.0005 0.098 0.013 0.059 2.96 3.96 0.057 0.057 2.88 0.000 0.098 0.032
1988 6 0.02 9.65 78.58 84.60 -0.0007 0.085 0.006 -0.018 0.00 1.44 -0.019 0.001 0.16 0.000 0.065 0.014
BALANCE HÍDRICO EN EL EMBALSE ENTRE LA OFERTA Y DEMANDA DE AGUA
Presa Champaccocha
Ht Hv Ar VolHt VolHv ISD Índice de Satisfacción de la Demanda 84.07 %
4413.0 4421.8 100.0 0.00 0.29
Vemb Aemba Aporte
año mes i Vesci Prec. Evap. Dem. Pi-Evi Vdi Vsi Ai Bi Vembi Aembi Vvi Vri
auxi ux total
(Hm3) (mm) (mm) (mm) (Hm3)/Há (Hm3) (Hm3) (Hm3) (Ha) (Hm3) (Ha) (Hm3) (Hm3) (Hm3)
1988 7 0.01 4.04 88.66 59.49 -0.0008 0.059 0.003 -0.055 0.00 0.08 -0.055 0.001 0.16 0.000 0.004 0.007
1988 8 0.00 9.2 104.74 19.59 -0.0010 0.020 0.001 -0.017 0.00 0.08 -0.017 0.001 0.16 0.000 0.002 0.003
1988 9 0.00 58.62 112.16 25.65 -0.0005 0.026 0.000 -0.023 0.00 0.08 -0.023 0.001 0.16 0.000 0.002 0.002
1988 10 0.00 37.04 128.01 0.00 -0.0009 0.000 0.000 0.002 0.35 0.26 0.002 0.002 0.32 0.000 0.000 0.001
1988 11 0.00 48.79 130.62 0.00 -0.0008 0.000 0.000 0.002 0.37 0.34 0.002 0.002 0.33 0.000 0.000 0.000
1988 12 0.02 120.61 123.62 35.14 0.0000 0.035 0.004 -0.017 0.00 0.17 -0.017 0.001 0.16 0.000 0.017 0.020
1989 1 0.07 172.02 112.37 25.50 0.0006 0.025 0.020 0.029 1.90 1.03 0.029 0.029 1.91 0.000 0.025 0.074
1989 2 0.08 134.33 96.49 55.78 0.0004 0.056 0.052 0.006 0.70 1.31 0.006 0.006 0.71 0.000 0.056 0.085
1989 3 0.19 235.4 104.32 78.51 0.0013 0.079 0.053 0.061 3.00 1.85 0.062 0.062 3.03 0.000 0.079 0.188
1989 4 0.13 112.5 88.13 99.06 0.0002 0.099 0.029 0.067 3.18 3.11 0.066 0.066 3.16 0.000 0.099 0.133
1989 5 0.07 18.81 85.01 98.36 -0.0007 0.098 0.013 0.026 1.79 2.47 0.024 0.024 1.69 0.000 0.098 0.070
1989 6 0.03 3.52 75.05 84.60 -0.0007 0.085 0.006 -0.035 0.00 0.84 -0.035 0.001 0.16 0.000 0.049 0.031
1989 7 0.01 0.41 84.04 59.49 -0.0008 0.059 0.003 -0.048 0.00 0.08 -0.048 0.001 0.16 0.000 0.011 0.014
1989 8 0.01 27.22 94.96 19.59 -0.0007 0.020 0.001 -0.014 0.00 0.08 -0.014 0.001 0.16 0.000 0.005 0.006
1989 9 0.00 40.88 110.08 25.65 -0.0007 0.026 0.000 -0.023 0.00 0.08 -0.023 0.001 0.16 0.000 0.002 0.003
1989 10 0.00 30.16 117.60 0.00 -0.0009 0.000 0.000 0.002 0.31 0.24 0.001 0.001 0.28 0.000 0.000 0.001
1989 11 0.00 37.38 128.63 0.00 -0.0009 0.000 0.000 0.002 0.30 0.29 0.001 0.001 0.27 0.000 0.000 0.000
1989 12 0.00 6.4 134.38 35.14 -0.0013 0.035 0.004 -0.038 0.00 0.13 -0.038 0.001 0.16 0.000 0.000 0.000
1990 1 0.07 186.35 118.38 25.50 0.0007 0.025 0.020 0.021 1.57 0.87 0.022 0.022 1.59 0.000 0.025 0.066
1990 2 0.07 117.25 103.35 55.78 0.0001 0.056 0.052 -0.015 0.00 0.79 -0.015 0.001 0.16 0.000 0.040 0.071
1990 3 0.06 113.78 111.36 78.51 0.0000 0.079 0.053 -0.070 0.00 0.08 -0.070 0.001 0.16 0.000 0.008 0.061
1990 4 0.04 78 99.30 99.06 -0.0002 0.099 0.029 -0.087 0.00 0.08 -0.087 0.001 0.16 0.000 0.012 0.041
1990 5 0.02 34.51 89.87 98.36 -0.0006 0.098 0.013 -0.090 0.00 0.08 -0.090 0.001 0.16 0.000 0.008 0.021
1990 6 0.01 26.34 75.93 84.60 -0.0005 0.085 0.006 -0.081 0.00 0.08 -0.081 0.001 0.16 0.000 0.003 0.009
1990 7 0.00 2.69 85.54 59.49 -0.0008 0.059 0.003 -0.057 0.00 0.08 -0.058 0.001 0.16 0.000 0.001 0.004
1990 8 0.00 27.59 97.14 19.59 -0.0007 0.020 0.001 -0.018 0.00 0.08 -0.018 0.001 0.16 0.000 0.001 0.002
1990 9 0.01 77.04 106.30 25.65 -0.0003 0.026 0.000 -0.020 0.00 0.08 -0.020 0.001 0.16 0.000 0.005 0.005
1990 10 0.01 91.24 117.95 0.00 -0.0003 0.000 0.000 0.013 1.17 0.67 0.013 0.013 1.16 0.000 0.000 0.013
1990 11 0.08 190.37 109.64 0.00 0.0008 0.000 0.000 0.093 3.93 2.54 0.095 0.095 3.96 0.000 0.000 0.082
1990 12 0.09 130.46 107.86 35.14 0.0002 0.035 0.004 0.141 5.08 4.52 0.141 0.141 5.07 0.000 0.035 0.085
1991 1 0.13 197.89 109.78 25.50 0.0009 0.025 0.020 0.227 6.83 5.95 0.229 0.229 6.87 0.000 0.025 0.133
1991 2 0.21 228.4 92.57 55.78 0.0014 0.056 0.052 0.327 8.59 7.73 0.332 0.294 8.03 0.038 0.056 0.211
1991 3 0.14 110.62 95.23 78.51 0.0002 0.079 0.053 0.302 8.17 8.10 0.299 0.294 8.03 0.005 0.079 0.137
1991 4 0.07 20.43 84.84 99.06 -0.0006 0.099 0.029 0.240 7.08 7.55 0.233 0.233 6.95 0.000 0.099 0.068
1991 5 0.03 42.06 81.88 98.36 -0.0004 0.098 0.013 0.154 5.37 6.16 0.151 0.151 5.30 0.000 0.098 0.029
1991 6 0.01 7.14 76.61 84.60 -0.0007 0.085 0.006 0.075 3.42 4.36 0.071 0.071 3.32 0.000 0.085 0.011
1991 7 0.01 0 86.39 59.49 -0.0009 0.059 0.003 0.016 1.29 2.31 0.014 0.014 1.18 0.000 0.059 0.004
1991 8 0.00 0 97.07 19.59 -0.0010 0.020 0.001 -0.004 0.00 0.59 -0.005 0.001 0.16 0.000 0.014 0.002
1991 9 0.00 65.84 108.95 25.65 -0.0004 0.026 0.000 -0.022 0.00 0.08 -0.022 0.001 0.16 0.000 0.003 0.003
1991 10 0.00 57.85 121.09 0.00 -0.0006 0.000 0.000 0.004 0.53 0.35 0.004 0.004 0.51 0.000 0.000 0.003
1991 11 0.00 9.35 124.06 0.00 -0.0011 0.000 0.000 0.005 0.64 0.58 0.004 0.004 0.59 0.000 0.000 0.001
1991 12 0.00 3.37 131.84 35.14 -0.0013 0.035 0.004 -0.034 0.00 0.29 -0.034 0.001 0.16 0.000 0.000 0.000
1992 1 0.03 131.06 122.21 25.50 0.0001 0.025 0.020 -0.019 0.00 0.08 -0.019 0.001 0.16 0.000 0.006 0.026
1992 2 0.17 255.41 108.22 55.78 0.0015 0.056 0.052 0.059 2.95 1.56 0.060 0.060 2.99 0.000 0.056 0.168
1992 3 0.12 117.1 105.37 78.51 0.0001 0.079 0.053 0.050 2.64 2.82 0.049 0.049 2.61 0.000 0.079 0.120
1992 4 0.08 82.65 91.10 99.06 -0.0001 0.099 0.029 -0.004 0.00 1.31 -0.004 0.001 0.16 0.000 0.094 0.075
1992 5 0.04 2 96.11 98.36 -0.0009 0.098 0.013 -0.072 0.00 0.08 -0.072 0.001 0.16 0.000 0.025 0.039
1992 6 0.02 6.36 82.62 84.60 -0.0008 0.085 0.006 -0.073 0.00 0.08 -0.073 0.001 0.16 0.000 0.011 0.017
1992 7 0.01 3.34 88.13 59.49 -0.0008 0.059 0.003 -0.054 0.00 0.08 -0.054 0.001 0.16 0.000 0.005 0.007
1992 8 0.00 17.79 97.58 19.59 -0.0008 0.020 0.001 -0.017 0.00 0.08 -0.017 0.001 0.16 0.000 0.002 0.003
1992 9 0.00 1.69 109.54 25.65 -0.0011 0.026 0.000 -0.024 0.00 0.08 -0.024 0.001 0.16 0.000 0.001 0.001
1992 10 0.00 70.07 120.70 0.00 -0.0005 0.000 0.000 0.004 0.52 0.34 0.003 0.003 0.50 0.000 0.000 0.003
1992 11 0.01 90.04 126.82 0.00 -0.0004 0.000 0.000 0.014 1.18 0.84 0.013 0.013 1.16 0.000 0.000 0.010
1992 12 0.02 91.83 123.65 35.14 -0.0003 0.035 0.004 -0.009 0.00 0.58 -0.010 0.001 0.16 0.000 0.025 0.017
1993 1 0.31 365.45 122.33 25.50 0.0024 0.025 0.020 0.265 7.52 3.84 0.269 0.269 7.61 0.000 0.025 0.314
1993 2 0.30 244.58 95.35 55.78 0.0015 0.056 0.052 0.462 10.65 9.13 0.465 0.294 8.03 0.172 0.056 0.304
1993 3 0.36 261.77 99.04 78.51 0.0016 0.079 0.053 0.523 11.50 9.77 0.525 0.294 8.03 0.232 0.079 0.364
1993 4 0.30 196.84 86.11 99.06 0.0011 0.099 0.029 0.469 10.75 9.39 0.469 0.294 8.03 0.175 0.099 0.303
1993 5 0.15 1.04 83.46 98.36 -0.0008 0.098 0.013 0.334 8.70 8.37 0.322 0.294 8.03 0.029 0.098 0.140
1993 6 0.07 1.84 74.86 84.60 -0.0007 0.085 0.006 0.270 7.62 7.82 0.262 0.262 7.48 0.000 0.085 0.059
1993 7 0.03 6.48 79.37 59.49 -0.0007 0.059 0.003 0.229 6.88 7.18 0.223 0.223 6.77 0.000 0.059 0.023
1993 8 0.01 15.16 85.96 19.59 -0.0007 0.020 0.001 0.215 6.62 6.69 0.210 0.210 6.52 0.000 0.020 0.008
1993 9 0.02 108.31 95.29 25.65 0.0001 0.026 0.000 0.205 6.42 6.47 0.206 0.206 6.43 0.000 0.026 0.021
1993 10 0.05 137.35 109.59 0.00 0.0003 0.000 0.000 0.252 7.30 6.86 0.253 0.253 7.31 0.000 0.000 0.047
1993 11 0.04 96.99 107.72 0.00 -0.0001 0.000 0.000 0.296 8.06 7.69 0.294 0.294 8.03 0.000 0.000 0.041
1993 12 0.17 258.24 108.59 35.14 0.0015 0.035 0.004 0.421 10.05 9.04 0.429 0.294 8.03 0.135 0.035 0.175
1994 1 0.17 181.65 108.65 25.50 0.0007 0.025 0.020 0.418 10.01 9.02 0.419 0.294 8.03 0.125 0.025 0.171
1994 2 0.32 291.15 96.87 55.78 0.0019 0.056 0.052 0.506 11.28 9.66 0.514 0.294 8.03 0.220 0.056 0.328
1994 3 0.20 134.18 101.55 78.51 0.0003 0.079 0.053 0.362 9.14 8.59 0.358 0.294 8.03 0.064 0.079 0.196
1994 4 0.12 85.84 86.31 99.06 0.0000 0.099 0.029 0.287 7.91 7.97 0.283 0.283 7.85 0.000 0.099 0.118
1994 5 0.06 23.44 84.40 98.36 -0.0006 0.098 0.013 0.233 6.95 7.40 0.227 0.227 6.84 0.000 0.098 0.055
1994 6 0.03 1.01 76.13 84.60 -0.0008 0.085 0.006 0.164 5.57 6.20 0.158 0.158 5.46 0.000 0.085 0.022
BALANCE HÍDRICO EN EL EMBALSE ENTRE LA OFERTA Y DEMANDA DE AGUA
Presa Champaccocha
Ht Hv Ar VolHt VolHv ISD Índice de Satisfacción de la Demanda 84.07 %
4413.0 4421.8 100.0 0.00 0.29
Vemb Aemba Aporte
año mes i Vesci Prec. Evap. Dem. Pi-Evi Vdi Vsi Ai Bi Vembi Aembi Vvi Vri
auxi ux total
(Hm3) (mm) (mm) (mm) (Hm3)/Há (Hm3) (Hm3) (Hm3) (Ha) (Hm3) (Ha) (Hm3) (Hm3) (Hm3)
1994 7 0.01 15.85 82.21 59.49 -0.0007 0.059 0.003 0.108 4.30 4.88 0.105 0.105 4.22 0.000 0.059 0.008
1994 8 0.01 15.59 93.58 19.59 -0.0008 0.020 0.001 0.089 3.81 4.01 0.086 0.086 3.73 0.000 0.020 0.002
1994 9 0.03 128.99 99.10 25.65 0.0003 0.026 0.000 0.088 3.78 3.75 0.089 0.089 3.80 0.000 0.026 0.028
1994 10 0.02 43.67 116.39 0.00 -0.0007 0.000 0.000 0.112 4.41 4.10 0.109 0.109 4.33 0.000 0.000 0.021
1994 11 0.05 147.63 114.41 0.00 0.0003 0.000 0.000 0.154 5.37 4.85 0.155 0.155 5.39 0.000 0.000 0.046
1994 12 0.04 73.51 117.83 35.14 -0.0004 0.035 0.004 0.153 5.35 5.37 0.150 0.150 5.28 0.000 0.035 0.034
1995 1 0.08 181.15 113.79 25.50 0.0007 0.025 0.020 0.183 5.98 5.63 0.185 0.185 6.03 0.000 0.025 0.081
1995 2 0.18 240.59 99.89 55.78 0.0014 0.056 0.052 0.261 7.46 6.74 0.266 0.266 7.54 0.000 0.056 0.188
1995 3 0.30 272.97 101.80 78.51 0.0017 0.079 0.053 0.434 10.24 8.89 0.439 0.294 8.03 0.145 0.079 0.306
1995 4 0.16 87.98 95.69 99.06 -0.0001 0.099 0.029 0.328 8.59 8.31 0.322 0.294 8.03 0.028 0.099 0.156
1995 5 0.08 14.07 86.62 98.36 -0.0007 0.098 0.013 0.263 7.50 7.76 0.255 0.255 7.35 0.000 0.098 0.073
1995 6 0.04 2.69 76.44 84.60 -0.0007 0.085 0.006 0.200 6.32 6.84 0.194 0.194 6.20 0.000 0.085 0.029
1995 7 0.02 0.01 83.09 59.49 -0.0008 0.059 0.003 0.148 5.23 5.72 0.143 0.143 5.12 0.000 0.059 0.010
1995 8 0.01 1.32 99.45 19.59 -0.0010 0.020 0.001 0.129 4.80 4.96 0.124 0.124 4.68 0.000 0.020 0.002
1995 9 0.00 14.25 99.47 25.65 -0.0009 0.026 0.000 0.101 4.12 4.40 0.097 0.097 4.02 0.000 0.026 -0.001
1995 10 0.00 48.92 117.74 0.00 -0.0007 0.000 0.000 0.098 4.06 4.04 0.096 0.096 3.98 0.000 0.000 -0.001
1995 11 0.01 92.75 111.67 0.00 -0.0002 0.000 0.000 0.104 4.21 4.10 0.103 0.103 4.18 0.000 0.000 0.008
1995 12 0.01 76.01 116.94 35.14 -0.0004 0.035 0.004 0.077 3.47 3.83 0.075 0.075 3.42 0.000 0.035 0.011
1996 1 0.17 274.81 116.40 25.50 0.0016 0.025 0.020 0.200 6.31 4.87 0.205 0.205 6.41 0.000 0.025 0.175
1996 2 0.34 305.37 102.32 55.78 0.0020 0.056 0.052 0.432 10.22 8.31 0.439 0.294 8.03 0.145 0.056 0.342
1996 3 0.21 147.95 101.66 78.51 0.0005 0.079 0.053 0.377 9.38 8.71 0.374 0.294 8.03 0.080 0.079 0.212
1996 4 0.12 62.77 87.57 99.06 -0.0002 0.099 0.029 0.286 7.89 7.96 0.280 0.280 7.80 0.000 0.099 0.115
1996 5 0.06 8.19 84.86 98.36 -0.0008 0.098 0.013 0.225 6.80 7.30 0.218 0.218 6.67 0.000 0.098 0.049
1996 6 0.02 0 76.96 84.60 -0.0008 0.085 0.006 0.152 5.33 6.00 0.147 0.147 5.22 0.000 0.085 0.020
1996 7 0.01 0.22 81.78 59.49 -0.0008 0.059 0.003 0.096 4.00 4.61 0.092 0.092 3.89 0.000 0.059 0.007
1996 8 0.00 12.28 91.86 19.59 -0.0008 0.020 0.001 0.076 3.46 3.68 0.073 0.073 3.37 0.000 0.020 0.002
1996 9 0.00 31.99 106.05 25.65 -0.0007 0.026 0.000 0.049 2.64 3.01 0.047 0.047 2.57 0.000 0.026 0.000
1996 10 0.00 26.8 122.75 0.00 -0.0010 0.000 0.000 0.048 2.59 2.58 0.046 0.046 2.51 0.000 0.000 -0.002
1996 11 0.01 94.57 122.73 0.00 -0.0003 0.000 0.000 0.054 2.80 2.66 0.054 0.054 2.78 0.000 0.000 0.008
1996 12 0.03 128.42 120.59 35.14 0.0001 0.035 0.004 0.047 2.56 2.67 0.047 0.047 2.55 0.000 0.035 0.033
1997 1 0.14 232.39 115.30 25.50 0.0012 0.025 0.020 0.139 5.03 3.79 0.141 0.141 5.08 0.000 0.025 0.140
1997 2 0.13 152.53 97.91 55.78 0.0005 0.056 0.052 0.165 5.59 5.34 0.165 0.165 5.60 0.000 0.056 0.132
1997 3 0.16 189.95 105.81 78.51 0.0008 0.079 0.053 0.194 6.19 5.90 0.195 0.195 6.22 0.000 0.079 0.162
1997 4 0.09 11.41 95.78 99.06 -0.0008 0.099 0.029 0.158 5.45 5.84 0.151 0.151 5.30 0.000 0.099 0.084
1997 5 0.04 26.15 85.93 98.36 -0.0006 0.098 0.013 0.079 3.55 4.43 0.076 0.076 3.46 0.000 0.098 0.037
1997 6 0.02 2.14 83.36 84.60 -0.0008 0.085 0.006 0.003 0.47 1.97 0.001 0.001 0.29 0.000 0.085 0.016
1997 7 0.01 7.08 90.58 59.49 -0.0008 0.059 0.003 -0.053 0.00 0.15 -0.053 0.001 0.16 0.000 0.006 0.008
1997 8 0.00 43.84 96.84 19.59 -0.0005 0.020 0.001 -0.017 0.00 0.08 -0.017 0.001 0.16 0.000 0.002 0.003
1997 9 0.00 61 110.57 25.65 -0.0005 0.026 0.000 -0.023 0.00 0.08 -0.023 0.001 0.16 0.000 0.002 0.003
1997 10 0.01 85.85 129.10 0.00 -0.0004 0.000 0.000 0.008 0.88 0.52 0.008 0.008 0.86 0.000 0.000 0.008
1997 11 0.01 87.14 123.55 0.00 -0.0004 0.000 0.000 0.022 1.58 1.22 0.021 0.021 1.56 0.000 0.000 0.013
1997 12 0.03 126.3 121.27 35.14 0.0001 0.035 0.004 0.014 1.22 1.39 0.014 0.014 1.22 0.000 0.035 0.032
1998 1 0.10 197.03 112.21 25.50 0.0008 0.025 0.020 0.070 3.28 2.25 0.071 0.071 3.31 0.000 0.025 0.102
1998 2 0.15 195.29 96.44 55.78 0.0010 0.056 0.052 0.116 4.51 3.91 0.118 0.118 4.55 0.000 0.056 0.155
1998 3 0.13 126.67 99.12 78.51 0.0003 0.079 0.053 0.111 4.38 4.46 0.110 0.110 4.36 0.000 0.079 0.124
1998 4 0.07 42.09 89.21 99.06 -0.0005 0.099 0.029 0.052 2.74 3.55 0.050 0.050 2.65 0.000 0.099 0.068
1998 5 0.03 1.64 92.16 98.36 -0.0009 0.098 0.013 -0.031 0.00 1.33 -0.032 0.001 0.16 0.000 0.066 0.030
1998 6 0.01 9.2 81.96 84.60 -0.0007 0.085 0.006 -0.076 0.00 0.08 -0.076 0.001 0.16 0.000 0.008 0.014
1998 7 0.01 3.03 89.44 59.49 -0.0009 0.059 0.003 -0.056 0.00 0.08 -0.056 0.001 0.16 0.000 0.003 0.006
1998 8 0.00 0.18 102.61 19.59 -0.0010 0.020 0.001 -0.018 0.00 0.08 -0.018 0.001 0.16 0.000 0.001 0.003
1998 9 0.00 11.58 115.47 25.65 -0.0010 0.026 0.000 -0.024 0.00 0.08 -0.024 0.001 0.16 0.000 0.001 0.001
1998 10 0.00 22.37 125.78 0.00 -0.0010 0.000 0.000 0.001 0.22 0.19 0.001 0.001 0.19 0.000 0.000 0.000
1998 11 0.00 40.59 129.85 0.00 -0.0009 0.000 0.000 0.001 0.17 0.18 0.000 0.001 0.16 0.000 0.000 0.000
1998 12 0.00 6.28 123.45 35.14 -0.0012 0.035 0.004 -0.039 0.00 0.08 -0.039 0.001 0.16 0.000 0.000 0.000
1999 1 0.01 84.78 119.71 25.50 -0.0003 0.025 0.020 -0.039 0.00 0.08 -0.039 0.001 0.16 0.000 0.000 0.006
1999 2 0.29 345.83 100.51 55.78 0.0025 0.056 0.052 0.178 5.88 3.02 0.182 0.182 5.96 0.000 0.056 0.290
1999 3 0.26 214.03 98.41 78.51 0.0012 0.079 0.053 0.309 8.29 7.13 0.311 0.294 8.03 0.017 0.079 0.260
1999 4 0.20 140.98 86.82 99.06 0.0005 0.099 0.029 0.361 9.13 8.58 0.359 0.294 8.03 0.065 0.099 0.194
1999 5 0.11 54.7 82.09 98.36 -0.0003 0.098 0.013 0.292 8.01 8.02 0.287 0.287 7.91 0.000 0.098 0.105
1999 6 0.05 7.29 77.06 84.60 -0.0007 0.085 0.006 0.248 7.23 7.57 0.242 0.242 7.11 0.000 0.085 0.045
1999 7 0.02 0 79.70 59.49 -0.0008 0.059 0.003 0.202 6.36 6.73 0.196 0.196 6.24 0.000 0.059 0.017
1999 8 0.01 14.25 94.17 19.59 -0.0008 0.020 0.001 0.185 6.02 6.13 0.180 0.180 5.92 0.000 0.020 0.005
1999 9 0.00 22.46 96.80 25.65 -0.0007 0.026 0.000 0.159 5.47 5.69 0.154 0.154 5.38 0.000 0.026 0.000
1999 10 0.01 99.71 106.66 0.00 -0.0001 0.000 0.000 0.167 5.65 5.51 0.167 0.167 5.64 0.000 0.000 0.012
1999 11 0.02 79.05 115.55 0.00 -0.0004 0.000 0.000 0.183 5.97 5.80 0.180 0.180 5.92 0.000 0.000 0.014
1999 12 0.15 258.66 110.89 35.14 0.0015 0.035 0.004 0.292 8.00 6.96 0.299 0.294 8.03 0.005 0.035 0.158
2000 1 0.27 271.3 105.32 25.50 0.0017 0.025 0.020 0.521 11.49 9.76 0.528 0.294 8.03 0.234 0.025 0.279
2000 2 0.55 387.04 97.50 55.78 0.0029 0.056 0.052 0.737 14.25 11.14 0.746 0.294 8.03 0.452 0.056 0.561
2000 3 0.44 265.32 97.51 78.51 0.0017 0.079 0.053 0.604 12.59 10.31 0.604 0.294 8.03 0.311 0.079 0.443
2000 4 0.20 27.37 86.34 99.06 -0.0006 0.099 0.029 0.364 9.18 8.60 0.353 0.294 8.03 0.059 0.099 0.187
2000 5 0.09 24.05 82.55 98.36 -0.0006 0.098 0.013 0.270 7.62 7.82 0.263 0.263 7.49 0.000 0.098 0.081
2000 6 0.04 1.12 74.61 84.60 -0.0007 0.085 0.006 0.211 6.53 7.01 0.205 0.205 6.41 0.000 0.085 0.032
2000 7 0.02 2.74 78.54 59.49 -0.0008 0.059 0.003 0.159 5.48 5.95 0.155 0.155 5.38 0.000 0.059 0.012
BALANCE HÍDRICO EN EL EMBALSE ENTRE LA OFERTA Y DEMANDA DE AGUA
Presa Champaccocha
Ht Hv Ar VolHt VolHv ISD Índice de Satisfacción de la Demanda 84.07 %
4413.0 4421.8 100.0 0.00 0.29
Vemb Aemba Aporte
año mes i Vesci Prec. Evap. Dem. Pi-Evi Vdi Vsi Ai Bi Vembi Aembi Vvi Vri
auxi ux total
(Hm3) (mm) (mm) (mm) (Hm3)/Há (Hm3) (Hm3) (Hm3) (Ha) (Hm3) (Ha) (Hm3) (Hm3) (Hm3)
2000 8 0.01 6.71 88.16 19.59 -0.0008 0.020 0.001 0.141 5.08 5.23 0.137 0.137 4.98 0.000 0.020 0.003
2000 9 0.00 55.62 101.60 25.65 -0.0005 0.026 0.000 0.115 4.47 4.73 0.113 0.113 4.41 0.000 0.026 0.002
2000 10 0.04 152.13 107.18 0.00 0.0004 0.000 0.000 0.153 5.35 4.88 0.155 0.155 5.39 0.000 0.000 0.042
2000 11 0.03 57.35 120.37 0.00 -0.0006 0.000 0.000 0.189 6.09 5.74 0.184 0.184 6.01 0.000 0.000 0.030
2000 12 0.04 129.87 113.57 35.14 0.0002 0.035 0.004 0.185 6.02 6.01 0.185 0.185 6.02 0.000 0.035 0.040
2001 1 0.11 204.87 107.64 25.50 0.0010 0.025 0.020 0.252 7.30 6.66 0.256 0.256 7.37 0.000 0.025 0.116
2001 2 0.09 117.13 98.33 55.78 0.0002 0.056 0.052 0.242 7.11 7.24 0.240 0.240 7.09 0.000 0.056 0.093
2001 3 0.11 174.63 101.61 78.51 0.0007 0.079 0.053 0.223 6.77 6.93 0.226 0.226 6.81 0.000 0.079 0.117
2001 4 0.07 6.98 90.17 99.06 -0.0008 0.099 0.029 0.167 5.64 6.22 0.160 0.160 5.49 0.000 0.099 0.063
2001 5 0.03 10.35 83.21 98.36 -0.0007 0.098 0.013 0.079 3.53 4.51 0.075 0.075 3.42 0.000 0.098 0.027
2001 6 0.01 5.15 76.53 84.60 -0.0007 0.085 0.006 -0.002 0.00 1.71 -0.003 0.001 0.16 0.000 0.081 0.012
2001 7 0.01 7.5 85.85 59.49 -0.0008 0.059 0.003 -0.056 0.00 0.08 -0.056 0.001 0.16 0.000 0.003 0.006
2001 8 0.00 23.19 97.62 19.59 -0.0007 0.020 0.001 -0.018 0.00 0.08 -0.018 0.001 0.16 0.000 0.001 0.003
2001 9 0.00 72.91 107.04 25.65 -0.0003 0.026 0.000 -0.021 0.00 0.08 -0.021 0.001 0.16 0.000 0.004 0.004
2001 10 0.01 85.2 122.62 0.00 -0.0004 0.000 0.000 0.011 1.01 0.59 0.010 0.010 0.99 0.000 0.000 0.010
2001 11 0.01 76.45 124.65 0.00 -0.0005 0.000 0.000 0.022 1.59 1.29 0.021 0.021 1.56 0.000 0.000 0.011
2001 12 0.03 123.95 122.34 35.14 0.0000 0.035 0.004 0.011 1.01 1.29 0.011 0.011 1.01 0.000 0.035 0.029
2002 1 0.02 58.99 125.53 25.50 -0.0007 0.025 0.020 -0.011 0.00 0.50 -0.011 0.001 0.16 0.000 0.014 0.023
2002 2 0.19 282.77 104.95 55.78 0.0018 0.056 0.052 0.085 3.69 1.93 0.087 0.087 3.75 0.000 0.056 0.194
2002 3 0.21 203.69 99.34 78.51 0.0010 0.079 0.053 0.167 5.65 4.70 0.169 0.169 5.68 0.000 0.079 0.214
2002 4 0.14 98.06 89.87 99.06 0.0001 0.099 0.029 0.178 5.88 5.78 0.176 0.176 5.82 0.000 0.099 0.135
2002 5 0.07 53.21 83.54 98.36 -0.0003 0.098 0.013 0.138 5.02 5.42 0.135 0.135 4.94 0.000 0.098 0.071
2002 6 0.03 6.48 74.33 84.60 -0.0007 0.085 0.006 0.079 3.53 4.24 0.075 0.075 3.44 0.000 0.085 0.031
2002 7 0.01 24.85 82.33 59.49 -0.0006 0.059 0.003 0.028 1.86 2.65 0.027 0.027 1.79 0.000 0.059 0.013
2002 8 0.01 5.9 93.10 19.59 -0.0009 0.020 0.001 0.012 1.12 1.46 0.011 0.011 1.04 0.000 0.020 0.005
2002 9 0.01 68.39 105.17 25.65 -0.0004 0.026 0.000 -0.009 0.00 0.52 -0.010 0.001 0.16 0.000 0.016 0.005
2002 10 0.00 30.47 122.14 0.00 -0.0009 0.000 0.000 0.004 0.57 0.37 0.004 0.004 0.54 0.000 0.000 0.003
2002 11 0.01 74.85 121.65 0.00 -0.0005 0.000 0.000 0.009 0.90 0.72 0.008 0.008 0.87 0.000 0.000 0.005
2002 12 0.01 83.46 128.43 35.14 -0.0004 0.035 0.004 -0.021 0.00 0.44 -0.021 0.001 0.16 0.000 0.013 0.010
2003 1 0.04 137.4 124.37 25.50 0.0001 0.025 0.020 -0.008 0.00 0.08 -0.008 0.001 0.16 0.000 0.017 0.037
2003 2 0.18 261.5 107.43 55.78 0.0015 0.056 0.052 0.073 3.37 1.76 0.074 0.074 3.41 0.000 0.056 0.182
2003 3 0.23 222.35 106.47 78.51 0.0012 0.079 0.053 0.172 5.74 4.58 0.173 0.173 5.77 0.000 0.079 0.231
2003 4 0.13 75.01 92.33 99.06 -0.0002 0.099 0.029 0.176 5.83 5.80 0.172 0.172 5.75 0.000 0.099 0.127
2003 5 0.06 1.06 84.08 98.36 -0.0008 0.098 0.013 0.125 4.70 5.23 0.119 0.119 4.57 0.000 0.098 0.058
2003 6 0.03 0 76.93 84.60 -0.0008 0.085 0.006 0.057 2.87 3.72 0.053 0.053 2.77 0.000 0.085 0.025
2003 7 0.01 0 85.69 59.49 -0.0009 0.059 0.003 0.004 0.52 1.64 0.002 0.002 0.37 0.000 0.059 0.011
2003 8 0.01 0.18 99.09 19.59 -0.0010 0.020 0.001 -0.013 0.00 0.19 -0.013 0.001 0.16 0.000 0.006 0.005
2003 9 0.00 1.15 107.89 25.65 -0.0011 0.026 0.000 -0.023 0.00 0.08 -0.023 0.001 0.16 0.000 0.002 0.002
2003 10 0.00 49.76 129.32 0.00 -0.0008 0.000 0.000 0.002 0.31 0.24 0.001 0.001 0.29 0.000 0.000 0.001
2003 11 0.00 55.38 130.18 0.00 -0.0007 0.000 0.000 0.002 0.37 0.33 0.002 0.002 0.34 0.000 0.000 0.001
2003 12 0.06 185.51 123.57 35.14 0.0006 0.035 0.004 0.027 1.83 1.09 0.028 0.028 1.85 0.000 0.035 0.065
2004 1 0.08 141.37 116.83 25.50 0.0002 0.025 0.020 0.065 3.14 2.49 0.065 0.065 3.14 0.000 0.025 0.083
2004 2 0.26 302.2 102.31 55.78 0.0020 0.056 0.052 0.219 6.69 4.91 0.224 0.224 6.78 0.000 0.056 0.267
2004 3 0.32 259.47 101.03 78.51 0.0016 0.079 0.053 0.408 9.86 8.32 0.412 0.294 8.03 0.118 0.079 0.319
2004 4 0.16 62.21 91.06 99.06 -0.0003 0.099 0.029 0.325 8.55 8.29 0.318 0.294 8.03 0.024 0.099 0.152
2004 5 0.07 10.77 85.57 98.36 -0.0007 0.098 0.013 0.257 7.39 7.71 0.249 0.249 7.25 0.000 0.098 0.067
2004 6 0.03 5.2 75.16 84.60 -0.0007 0.085 0.006 0.192 6.15 6.70 0.186 0.186 6.04 0.000 0.085 0.027
2004 7 0.01 3.55 78.43 59.49 -0.0007 0.059 0.003 0.138 5.02 5.53 0.134 0.134 4.92 0.000 0.059 0.010
2004 8 0.01 19.34 85.49 19.59 -0.0007 0.020 0.001 0.119 4.58 4.75 0.116 0.116 4.50 0.000 0.020 0.003
2004 9 0.04 144.28 93.60 25.65 0.0005 0.026 0.000 0.129 4.80 4.65 0.131 0.131 4.84 0.000 0.026 0.040
2004 10 0.03 59.57 108.23 0.00 -0.0005 0.000 0.000 0.163 5.57 5.20 0.160 0.160 5.50 0.000 0.000 0.030
2004 11 0.02 93.55 114.85 0.00 -0.0002 0.000 0.000 0.184 5.99 5.75 0.182 0.182 5.96 0.000 0.000 0.022
2004 12 0.12 225.93 105.73 35.14 0.0012 0.035 0.004 0.264 7.51 6.73 0.269 0.269 7.60 0.000 0.035 0.126
2005 1 0.08 68.69 111.25 25.50 -0.0004 0.025 0.020 0.302 8.17 7.88 0.296 0.294 8.03 0.002 0.025 0.073
2005 2 0.06 124.32 94.24 55.78 0.0003 0.056 0.052 0.246 7.19 7.61 0.247 0.247 7.20 0.000 0.056 0.061
2005 3 0.16 240.38 99.47 78.51 0.0014 0.079 0.053 0.278 7.76 7.48 0.284 0.284 7.87 0.000 0.079 0.169
2005 4 0.09 11.47 89.29 99.06 -0.0008 0.099 0.029 0.247 7.21 7.54 0.239 0.239 7.06 0.000 0.099 0.083
2005 5 0.04 2.77 87.17 98.36 -0.0008 0.098 0.013 0.167 5.65 6.35 0.161 0.161 5.51 0.000 0.098 0.034
2005 6 0.02 1.36 77.35 84.60 -0.0008 0.085 0.006 0.088 3.78 4.65 0.084 0.084 3.68 0.000 0.085 0.014
2005 7 0.01 7.06 81.32 59.49 -0.0007 0.059 0.003 0.030 1.93 2.80 0.028 0.028 1.84 0.000 0.059 0.006
2005 8 0.00 12.8 101.65 19.59 -0.0009 0.020 0.001 0.010 0.99 1.41 0.009 0.009 0.91 0.000 0.020 0.002
2005 9 0.00 3.54 109.85 25.65 -0.0011 0.026 0.000 -0.015 0.00 0.46 -0.016 0.001 0.16 0.000 0.009 0.001
2005 10 0.00 1.45 123.68 0.00 -0.0012 0.000 0.000 0.001 0.24 0.20 0.001 0.001 0.21 0.000 0.000 0.000
2005 11 0.00 18.41 131.93 0.00 -0.0011 0.000 0.000 0.001 0.19 0.20 0.001 0.001 0.16 0.000 0.000 0.000
2005 12 0.02 111.65 123.80 35.14 -0.0001 0.035 0.004 -0.023 0.00 0.08 -0.023 0.001 0.16 0.000 0.012 0.016
2006 1 0.02 62.61 120.28 25.50 -0.0006 0.025 0.020 -0.028 0.00 0.08 -0.028 0.001 0.16 0.000 0.000 0.016
2006 2 0.07 180.75 104.22 55.78 0.0008 0.056 0.052 -0.036 0.00 0.08 -0.035 0.001 0.16 0.000 0.020 0.072
2006 3 0.10 152.26 108.21 78.51 0.0004 0.079 0.053 -0.031 0.00 0.08 -0.031 0.001 0.16 0.000 0.047 0.100
2006 4 0.06 62.08 97.31 99.06 -0.0004 0.099 0.029 -0.063 0.00 0.08 -0.063 0.001 0.16 0.000 0.036 0.065
2006 5 0.03 14.92 94.56 98.36 -0.0008 0.098 0.013 -0.080 0.00 0.08 -0.081 0.001 0.16 0.000 0.017 0.031
2006 6 0.01 10.5 84.35 84.60 -0.0007 0.085 0.006 -0.076 0.00 0.08 -0.076 0.001 0.16 0.000 0.008 0.013
2006 7 0.01 7.29 90.46 59.49 -0.0008 0.059 0.003 -0.056 0.00 0.08 -0.056 0.001 0.16 0.000 0.003 0.006
BALANCE HÍDRICO EN EL EMBALSE ENTRE LA OFERTA Y DEMANDA DE AGUA
Presa Champaccocha
Ht Hv Ar VolHt VolHv ISD Índice de Satisfacción de la Demanda 84.07 %
4413.0 4421.8 100.0 0.00 0.29
Vemb Aemba Aporte
año mes i Vesci Prec. Evap. Dem. Pi-Evi Vdi Vsi Ai Bi Vembi Aembi Vvi Vri
auxi ux total
(Hm3) (mm) (mm) (mm) (Hm3)/Há (Hm3) (Hm3) (Hm3) (Ha) (Hm3) (Ha) (Hm3) (Hm3) (Hm3)
2006 8 0.00 10.52 99.52 19.59 -0.0009 0.020 0.001 -0.018 0.00 0.08 -0.018 0.001 0.16 0.000 0.001 0.003
2006 9 0.01 87.69 112.55 25.65 -0.0002 0.026 0.000 -0.017 0.00 0.08 -0.017 0.001 0.16 0.000 0.008 0.008
2006 10 0.01 37.12 126.37 0.00 -0.0009 0.000 0.000 0.008 0.86 0.51 0.008 0.008 0.83 0.000 0.000 0.007
2006 11 0.02 106.22 120.05 0.00 -0.0001 0.000 0.000 0.024 1.70 1.27 0.024 0.024 1.69 0.000 0.000 0.017
2006 12 0.05 146.49 118.48 35.14 0.0003 0.035 0.004 0.034 2.09 1.89 0.034 0.034 2.10 0.000 0.035 0.049
2007 1 0.18 261.62 114.59 25.50 0.0015 0.025 0.020 0.171 5.74 3.92 0.175 0.175 5.80 0.000 0.025 0.185
2007 2 0.32 290.39 99.52 55.78 0.0019 0.056 0.052 0.386 9.52 7.66 0.392 0.294 8.03 0.098 0.056 0.325
2007 3 0.35 257.04 98.72 78.51 0.0016 0.079 0.053 0.512 11.35 9.69 0.514 0.294 8.03 0.220 0.079 0.352
2007 4 0.22 128.42 86.76 99.06 0.0004 0.099 0.029 0.385 9.51 8.77 0.382 0.294 8.03 0.088 0.099 0.216
2007 5 0.11 8.89 84.56 98.36 -0.0008 0.098 0.013 0.297 8.08 8.06 0.287 0.287 7.92 0.000 0.098 0.105
2007 6 0.05 2.55 75.61 84.60 -0.0007 0.085 0.006 0.247 7.21 7.57 0.240 0.240 7.08 0.000 0.085 0.043
2007 7 0.02 0.28 80.70 59.49 -0.0008 0.059 0.003 0.200 6.32 6.70 0.194 0.194 6.20 0.000 0.059 0.016
2007 8 0.01 6.12 94.87 19.59 -0.0009 0.020 0.001 0.183 5.98 6.09 0.178 0.178 5.87 0.000 0.020 0.004
2007 9 0.00 19.21 94.26 25.65 -0.0008 0.026 0.000 0.156 5.41 5.64 0.152 0.152 5.32 0.000 0.026 0.000
2007 10 0.00 50.67 114.25 0.00 -0.0006 0.000 0.000 0.154 5.36 5.34 0.150 0.150 5.28 0.000 0.000 -0.001
2007 11 0.00 76.03 117.48 0.00 -0.0004 0.000 0.000 0.155 5.38 5.33 0.152 0.152 5.33 0.000 0.000 0.002
2007 12 0.06 175.64 112.39 35.14 0.0006 0.035 0.004 0.173 5.77 5.55 0.175 0.175 5.82 0.000 0.035 0.062
2008 1 0.19 254.85 105.88 25.50 0.0015 0.025 0.020 0.319 8.46 7.14 0.325 0.294 8.03 0.031 0.025 0.195
2008 2 0.14 126.1 98.10 55.78 0.0003 0.056 0.052 0.323 8.51 8.27 0.321 0.294 8.03 0.027 0.056 0.135
2008 3 0.09 94.02 101.60 78.51 -0.0001 0.079 0.053 0.251 7.28 7.65 0.248 0.248 7.22 0.000 0.079 0.086
2008 4 0.05 6.04 92.95 99.06 -0.0009 0.099 0.029 0.166 5.63 6.43 0.160 0.160 5.49 0.000 0.099 0.040
2008 5 0.02 0.21 84.83 98.36 -0.0008 0.098 0.013 0.068 3.23 4.36 0.064 0.064 3.11 0.000 0.098 0.016
2008 6 0.01 0.37 79.33 84.60 -0.0008 0.085 0.006 -0.017 0.00 1.56 -0.018 0.001 0.16 0.000 0.066 0.008
2008 7 0.00 0.15 90.10 59.49 -0.0009 0.059 0.003 -0.058 0.00 0.08 -0.058 0.001 0.16 0.000 0.001 0.004
2008 8 0.00 13.78 106.88 19.59 -0.0009 0.020 0.001 -0.018 0.00 0.08 -0.019 0.001 0.16 0.000 0.001 0.002
2008 9 0.00 23.35 112.96 25.65 -0.0009 0.026 0.000 -0.024 0.00 0.08 -0.025 0.001 0.16 0.000 0.001 0.001
2008 10 0.01 95.18 123.45 0.00 -0.0003 0.000 0.000 0.010 0.97 0.57 0.010 0.010 0.96 0.000 0.000 0.009
2008 11 0.01 80.36 128.35 0.00 -0.0005 0.000 0.000 0.024 1.67 1.31 0.023 0.023 1.64 0.000 0.000 0.013
2008 12 0.10 211.45 119.20 35.14 0.0009 0.035 0.004 0.081 3.58 2.61 0.082 0.082 3.63 0.000 0.035 0.099
2009 1 0.08 106.53 105.24 25.50 0.0000 0.025 0.020 0.118 4.54 4.08 0.116 0.116 4.51 0.000 0.025 0.080
2009 2 0.22 280.91 93.06 55.78 0.0019 0.056 0.052 0.232 6.94 5.72 0.238 0.238 7.05 0.000 0.056 0.230
2009 3 0.12 80.05 100.21 78.51 -0.0002 0.079 0.053 0.230 6.90 6.97 0.226 0.226 6.82 0.000 0.079 0.120
2009 4 0.06 65.27 85.72 99.06 -0.0002 0.099 0.029 0.162 5.53 6.17 0.159 0.159 5.47 0.000 0.099 0.061
2009 5 0.03 36.93 83.43 98.36 -0.0005 0.098 0.013 0.077 3.49 4.48 0.075 0.075 3.42 0.000 0.098 0.027
2009 6 0.01 9.18 81.14 84.60 -0.0007 0.085 0.006 -0.002 0.00 1.71 -0.004 0.001 0.16 0.000 0.080 0.012
2009 7 0.01 10.72 88.68 59.49 -0.0008 0.059 0.003 -0.056 0.00 0.08 -0.056 0.001 0.16 0.000 0.003 0.006
2009 8 0.00 14.22 104.91 19.59 -0.0009 0.020 0.001 -0.018 0.00 0.08 -0.018 0.001 0.16 0.000 0.001 0.002
2009 9 0.00 53.61 113.12 25.65 -0.0006 0.026 0.000 -0.024 0.00 0.08 -0.024 0.001 0.16 0.000 0.001 0.001
2009 10 0.00 50.06 127.17 0.00 -0.0008 0.000 0.000 0.001 0.29 0.23 0.001 0.001 0.27 0.000 0.000 0.001
2009 11 0.07 189.09 122.28 0.00 0.0007 0.000 0.000 0.069 3.25 1.76 0.070 0.070 3.27 0.000 0.000 0.069
2009 12 0.09 144.32 118.54 35.14 0.0003 0.035 0.004 0.117 4.53 3.90 0.117 0.117 4.52 0.000 0.035 0.087
2010 1 0.11 175.18 106.51 25.50 0.0007 0.025 0.020 0.184 5.99 5.26 0.185 0.185 6.02 0.000 0.025 0.113
2010 2 0.18 223.67 96.93 55.78 0.0013 0.056 0.052 0.262 7.47 6.75 0.266 0.266 7.54 0.000 0.056 0.188
2010 3 0.17 164.11 100.98 78.51 0.0006 0.079 0.053 0.303 8.19 7.86 0.303 0.294 8.03 0.009 0.079 0.170
2010 4 0.11 81.25 90.01 99.06 -0.0001 0.099 0.029 0.272 7.65 7.84 0.268 0.268 7.58 0.000 0.099 0.102
2010 5 0.05 13.13 84.86 98.36 -0.0007 0.098 0.013 0.209 6.50 7.04 0.203 0.203 6.37 0.000 0.098 0.047
2010 6 0.02 0 75.64 84.60 -0.0008 0.085 0.006 0.136 4.96 5.67 0.131 0.131 4.85 0.000 0.085 0.019
2010 7 0.01 0 83.43 59.49 -0.0008 0.059 0.003 0.079 3.54 4.20 0.076 0.076 3.44 0.000 0.059 0.007
2010 8 0.00 0.1 100.94 19.59 -0.0010 0.020 0.001 0.059 2.96 3.20 0.056 0.056 2.85 0.000 0.020 0.001
2010 9 0.00 19.16 105.86 25.65 -0.0009 0.026 0.000 0.032 2.02 2.44 0.030 0.030 1.94 0.000 0.026 0.000
2010 10 0.00 3.1 119.95 0.00 -0.0012 0.000 0.000 0.031 1.97 1.95 0.028 0.028 1.87 0.000 0.000 -0.001
2010 11 0.00 40.21 124.59 0.00 -0.0008 0.000 0.000 0.029 1.87 1.87 0.027 0.027 1.81 0.000 0.000 -0.001
2010 12 0.06 185.4 115.54 35.14 0.0007 0.035 0.004 0.052 2.73 2.27 0.053 0.053 2.77 0.000 0.035 0.066
2011 1 0.18 245.16 113.07 25.50 0.0013 0.025 0.020 0.192 6.15 4.46 0.195 0.195 6.21 0.000 0.025 0.187
2011 2 0.42 346.51 93.95 55.78 0.0025 0.056 0.052 0.510 11.33 8.77 0.519 0.294 8.03 0.225 0.056 0.432
2011 3 0.35 231.81 97.17 78.51 0.0013 0.079 0.053 0.511 11.35 9.69 0.512 0.294 8.03 0.218 0.079 0.350
2011 4 0.26 161.2 86.60 99.06 0.0007 0.099 0.029 0.425 10.11 9.07 0.423 0.294 8.03 0.129 0.099 0.258
2011 5 0.13 2.7 82.69 98.36 -0.0008 0.098 0.013 0.317 8.42 8.22 0.306 0.294 8.03 0.012 0.098 0.124
2011 6 0.06 5.27 75.25 84.60 -0.0007 0.085 0.006 0.262 7.48 7.76 0.255 0.255 7.36 0.000 0.085 0.052
2011 7 0.03 0 79.49 59.49 -0.0008 0.059 0.003 0.219 6.69 7.02 0.213 0.213 6.57 0.000 0.059 0.020
2011 8 0.01 10.03 92.32 19.59 -0.0008 0.020 0.001 0.203 6.38 6.48 0.198 0.198 6.27 0.000 0.020 0.006
2011 9 0.00 46.78 97.62 25.65 -0.0005 0.026 0.000 0.177 5.85 6.06 0.174 0.174 5.79 0.000 0.026 0.002
2011 10 0.00 62.32 110.99 0.00 -0.0005 0.000 0.000 0.177 5.86 5.82 0.174 0.174 5.80 0.000 0.000 0.001
2011 11 0.03 129.99 113.55 0.00 0.0002 0.000 0.000 0.201 6.33 6.07 0.201 0.201 6.34 0.000 0.000 0.027
2011 12 0.07 166.47 108.54 35.14 0.0006 0.035 0.004 0.235 6.98 6.66 0.237 0.237 7.01 0.000 0.035 0.075
2012 1 0.12 190.78 110.41 25.50 0.0008 0.025 0.020 0.313 8.36 7.69 0.316 0.294 8.03 0.022 0.025 0.125
2012 2 0.30 300.33 99.71 55.78 0.0020 0.056 0.052 0.481 10.92 9.47 0.489 0.294 8.03 0.195 0.056 0.304
2012 3 0.23 174.87 99.42 78.51 0.0008 0.079 0.053 0.387 9.54 8.78 0.386 0.294 8.03 0.093 0.079 0.225
2012 4 0.15 100.85 90.00 99.06 0.0001 0.099 0.029 0.311 8.32 8.17 0.307 0.294 8.03 0.013 0.099 0.142
2012 5 0.08 1.64 83.89 98.36 -0.0008 0.098 0.013 0.258 7.40 7.72 0.249 0.249 7.25 0.000 0.098 0.067
2012 6 0.03 7.67 73.63 84.60 -0.0007 0.085 0.006 0.192 6.16 6.70 0.187 0.187 6.05 0.000 0.085 0.028
2012 7 0.01 3.26 80.32 59.49 -0.0008 0.059 0.003 0.139 5.04 5.55 0.135 0.135 4.94 0.000 0.059 0.010
2012 8 0.01 10.89 93.56 19.59 -0.0008 0.020 0.001 0.120 4.60 4.77 0.116 0.116 4.50 0.000 0.020 0.002
2012 9 0.00 54.95 99.42 25.65 -0.0004 0.026 0.000 0.094 3.94 4.22 0.092 0.092 3.89 0.000 0.026 0.001
2012 10 0.00 54.44 115.50 0.00 -0.0006 0.000 0.000 0.094 3.94 3.91 0.091 0.091 3.88 0.000 0.000 0.000
2012 11 0.00 28.07 113.62 0.00 -0.0009 0.000 0.000 0.092 3.90 3.89 0.089 0.089 3.81 0.000 0.000 -0.002
2012 12 0.08 199.34 105.95 35.14 0.0009 0.035 0.004 0.128 4.79 4.30 0.131 0.131 4.85 0.000 0.035 0.081
2013 1 0.06 82.25 110.42 25.50 -0.0003 0.025 0.020 0.148 5.24 5.05 0.146 0.146 5.19 0.000 0.025 0.060
2013 2 0.10 183.52 91.30 55.78 0.0009 0.056 0.052 0.137 4.99 5.09 0.140 0.140 5.05 0.000 0.056 0.102
2013 3 0.06 30.87 101.35 78.51 -0.0007 0.079 0.053 0.071 3.30 4.18 0.067 0.067 3.19 0.000 0.079 0.059
2013 4 0.03 7.75 99.42 99.06 -0.0009 0.099 0.029 -0.034 0.00 1.59 -0.036 0.001 0.16 0.000 0.063 0.026
2013 5 0.01 25.53 91.47 98.36 -0.0007 0.098 0.013 -0.099 0.00 0.08 -0.099 0.001 0.16 0.000 0.000 0.012
BALANCE HÍDRICO EN EL EMBALSE ENTRE LA OFERTA Y DEMANDA DE AGUA
Presa Champaccocha
Ht Hv Ar VolHt VolHv ISD Índice de Satisfacción de la Demanda 84.07 %
4413.0 4421.8 100.0 0.00 0.29
Vemb Aemba Aporte
año mes i Vesci Prec. Evap. Dem. Pi-Evi Vdi Vsi Ai Bi Vembi Aembi Vvi Vri
auxi ux total
(Hm3) (mm) (mm) (mm) (Hm3)/Há (Hm3) (Hm3) (Hm3) (Ha) (Hm3) (Ha) (Hm3) (Hm3) (Hm3)
2013 6 0.01 8.37 80.57 84.60 -0.0007 0.085 0.006 -0.085 0.00 0.08 -0.085 0.001 0.16 0.000 0.000 0.005
2013 7 0.00 12.41 85.00 59.49 -0.0007 0.059 0.003 -0.059 0.00 0.08 -0.059 0.001 0.16 0.000 0.000 0.002
2013 8 0.00 27.4 101.23 19.59 -0.0007 0.020 0.001 -0.019 0.00 0.08 -0.019 0.001 0.16 0.000 0.000 0.001
2013 9 0.00 62.05 110.53 25.65 -0.0005 0.026 0.000 -0.023 0.00 0.08 -0.023 0.001 0.16 0.000 0.002 0.002
2013 10 0.00 46.42 121.16 0.00 -0.0007 0.000 0.000 0.002 0.36 0.26 0.002 0.002 0.34 0.000 0.000 0.001
2013 11 0.01 88.78 127.93 0.00 -0.0004 0.000 0.000 0.009 0.93 0.64 0.009 0.009 0.91 0.000 0.000 0.007
2013 12 0.03 125.07 121.31 35.14 0.0000 0.035 0.004 -0.001 0.00 0.46 -0.001 0.001 0.16 0.000 0.034 0.030
2014 1 0.10 194.42 116.88 25.50 0.0008 0.025 0.020 0.053 2.75 1.46 0.053 0.053 2.77 0.000 0.025 0.098
2014 2 0.16 203.45 97.68 55.78 0.0011 0.056 0.052 0.105 4.24 3.51 0.107 0.107 4.28 0.000 0.056 0.162
2014 3 0.32 285.77 98.08 78.51 0.0019 0.079 0.053 0.291 7.99 6.13 0.295 0.294 8.03 0.002 0.079 0.320
2014 4 0.16 67.09 88.98 99.06 -0.0002 0.099 0.029 0.322 8.50 8.26 0.315 0.294 8.03 0.021 0.099 0.150
2014 5 0.07 24.17 83.16 98.36 -0.0006 0.098 0.013 0.257 7.38 7.71 0.250 0.250 7.26 0.000 0.098 0.068
2014 6 0.03 5.08 75.77 84.60 -0.0007 0.085 0.006 0.192 6.17 6.72 0.187 0.187 6.06 0.000 0.085 0.027
2014 7 0.01 8.07 79.19 59.49 -0.0007 0.059 0.003 0.139 5.04 5.55 0.135 0.135 4.94 0.000 0.059 0.010
2014 8 0.01 20.57 90.98 19.59 -0.0007 0.020 0.001 0.120 4.60 4.77 0.117 0.117 4.52 0.000 0.020 0.003
Variación de Niveles en el Embalse
Volúmenes de déficit, entregados, y vertidos.
Nivel de Nivel de Volumen Vol entreg. para
año mes i Déficit Vri‐Vdi Volumen vertido Vvi
Embalse Emb. (%) disp. (%) riego Vri
1982 9 4419.5 |||||| ||||| ||
1982 10 4420.4 ||||||| ||||||
1982 11 4421.3 ||||||||| ||||||||
1982 12 4421.1 |||||||| |||||||| |||
1983 1 4421.5 ||||||||| ||||||||| ||
1983 2 4421.2 ||||||||| |||||||| |||||
1983 3 4420.3 ||||||| |||||| |||||||
1983 4 4419.0 ||||| |||| |||||||||
1983 5 4416.9 || | |||||||||
1983 6 4415.0 || |||||
1983 7 4415.0 |||||
1983 8 4415.0 |
1983 9 4415.0 | |
1983 10 4415.4
1983 11 4415.5
1983 12 4417.5 ||| || |||
1984 1 4418.5 ||||| ||| ||
1984 2 4421.8 |||||||||| |||||||||| ||||| |||||||
1984 3 4421.8 |||||||||| |||||||||| ||||||| |||||||||||||||||||||||||||||
1984 4 4421.8 |||||||||| |||||||||| ||||||||| |||||||||
1984 5 4421.7 ||||||||| ||||||||| |||||||||
1984 6 4421.0 |||||||| |||||||| ||||||||
1984 7 4420.2 ||||||| |||||| |||||
1984 8 4419.9 ||||||| |||||| |
1984 9 4419.4 |||||| ||||| ||
1984 10 4419.3 |||||| |||||
1984 11 4420.5 |||||||| |||||||
1984 12 4421.8 |||||||||| |||||||||| ||| |
1985 1 4421.8 |||||||||| |||||||||| || |||
1985 2 4421.1 |||||||| |||||||| |||||
1985 3 4419.8 ||||||| ||||| |||||||
1985 4 4418.4 ||||| ||| |||||||||
1985 5 4416.3 | | |||||||||
1985 6 4415.0 |||| ||||
1985 7 4415.0 |||||
1985 8 4415.0 |
1985 9 4415.0 |
1985 10 4415.3
1985 11 4415.3
1985 12 4415.4 |||
1986 1 4420.2 ||||||| |||||| ||
1986 2 4421.8 |||||||||| |||||||||| ||||| |||||||||||||||||||
1986 3 4421.8 |||||||||| |||||||||| ||||||| ||||||||||||
1986 4 4421.8 |||||||||| |||||||||| ||||||||| |
1986 5 4421.1 |||||||| |||||||| |||||||||
1986 6 4420.0 ||||||| |||||| ||||||||
1986 7 4418.9 ||||| |||| |||||
1986 8 4418.4 ||||| ||| |
1986 9 4418.0 |||| ||| ||
1986 10 4418.2 |||| |||
1986 11 4418.2 |||| |||
1986 12 4417.1 ||| | |||
1987 1 4418.9 ||||| |||| ||
1987 2 4421.3 ||||||||| |||||||| |||||
1987 3 4421.1 |||||||| |||||||| |||||||
1987 4 4419.9 ||||||| ||||| |||||||||
1987 5 4417.9 |||| || |||||||||
1987 6 4415.2 ||||||||
1987 7 4415.0 |||||
1987 8 4415.0 |
1987 9 4415.0 ||
1987 10 4415.0
1987 11 4415.1
1987 12 4415.0 |||
1988 1 4416.2 | | ||
1988 2 4420.2 ||||||| |||||| |||||
1988 3 4420.3 ||||||| |||||| |||||||
1988 4 4419.2 |||||| |||| |||||||||
1988 5 4417.2 ||| || |||||||||
1988 6 4415.0 | ||||||
‐ 1 ‐
Nivel de Nivel de Volumen Vol entreg. para
año mes i Déficit Vri‐Vdi Volumen vertido Vvi
Embalse Emb. (%) disp. (%) riego Vri
1988 7 4415.0 |||||
1988 8 4415.0 |
1988 9 4415.0 ||
1988 10 4415.0
1988 11 4415.0
1988 12 4415.0 | |
1989 1 4416.2 | | ||
1989 2 4415.3 |||||
1989 3 4417.3 ||| || |||||||
1989 4 4417.4 ||| || |||||||||
1989 5 4416.0 | |||||||||
1989 6 4415.0 ||| ||||
1989 7 4415.0 |||| |
1989 8 4415.0 |
1989 9 4415.0 ||
1989 10 4415.0
1989 11 4415.0
1989 12 4415.0 |||
1990 1 4415.9 | ||
1990 2 4415.0 | ||||
1990 3 4415.0 |||||||
1990 4 4415.0 |||||||| |
1990 5 4415.0 |||||||||
1990 6 4415.0 ||||||||
1990 7 4415.0 |||||
1990 8 4415.0 |
1990 9 4415.0 ||
1990 10 4415.5
1990 11 4418.2 |||| |||
1990 12 4419.3 |||||| |||| |||
1991 1 4420.9 |||||||| |||||||| ||
1991 2 4421.8 |||||||||| |||||||||| ||||| ||||
1991 3 4421.8 |||||||||| |||||||||| |||||||
1991 4 4420.8 |||||||| ||||||| |||||||||
1991 5 4419.3 |||||| |||| |||||||||
1991 6 4417.2 ||| || ||||||||
1991 7 4415.1 |||||
1991 8 4415.0 |
1991 9 4415.0 ||
1991 10 4415.1
1991 11 4415.0
1991 12 4415.0 |||
1992 1 4415.0 |
1992 2 4417.3 ||| || |||||
1992 3 4416.9 || | |||||||
1992 4 4415.0 |||||||||
1992 5 4415.0 ||||||| ||
1992 6 4415.0 ||||||| |
1992 7 4415.0 |||||
1992 8 4415.0 |
1992 9 4415.0 ||
1992 10 4415.1
1992 11 4415.5
1992 12 4415.0 | ||
1993 1 4421.6 ||||||||| ||||||||| ||
1993 2 4421.8 |||||||||| |||||||||| ||||| ||||||||||||||||||
1993 3 4421.8 |||||||||| |||||||||| ||||||| |||||||||||||||||||||||
1993 4 4421.8 |||||||||| |||||||||| ||||||||| |||||||||||||||||
1993 5 4421.8 |||||||||| |||||||||| ||||||||| |||
1993 6 4421.3 ||||||||| |||||||| ||||||||
1993 7 4420.7 |||||||| ||||||| |||||
1993 8 4420.5 |||||||| ||||||| |
1993 9 4420.4 ||||||| |||||| ||
1993 10 4421.2 ||||||||| ||||||||
1993 11 4421.8 |||||||||| ||||||||||
1993 12 4421.8 |||||||||| |||||||||| ||| |||||||||||||
1994 1 4421.8 |||||||||| |||||||||| || ||||||||||||
1994 2 4421.8 |||||||||| |||||||||| ||||| |||||||||||||||||||||
1994 3 4421.8 |||||||||| |||||||||| ||||||| ||||||
1994 4 4421.6 ||||||||| ||||||||| |||||||||
1994 5 4420.7 |||||||| ||||||| |||||||||
1994 6 4419.4 |||||| ||||| ||||||||
1994 7 4418.2 |||| ||| |||||
‐ 2 ‐
Nivel de Nivel de Volumen Vol entreg. para
año mes i Déficit Vri‐Vdi Volumen vertido Vvi
Embalse Emb. (%) disp. (%) riego Vri
1994 8 4417.7 ||| || |
1994 9 4417.8 |||| || ||
1994 10 4418.3 |||| |||
1994 11 4419.4 |||||| |||||
1994 12 4419.3 |||||| |||| |||
1995 1 4420.0 ||||||| |||||| ||
1995 2 4421.4 ||||||||| ||||||||| |||||
1995 3 4421.8 |||||||||| |||||||||| ||||||| ||||||||||||||
1995 4 4421.8 |||||||||| |||||||||| ||||||||| ||
1995 5 4421.2 ||||||||| |||||||| |||||||||
1995 6 4420.2 ||||||| |||||| ||||||||
1995 7 4419.1 |||||| |||| |||||
1995 8 4418.7 ||||| |||| |
1995 9 4418.0 |||| ||| ||
1995 10 4418.0 |||| |||
1995 11 4418.2 |||| |||
1995 12 4417.4 ||| || |||
1996 1 4420.4 ||||||| |||||| ||
1996 2 4421.8 |||||||||| |||||||||| ||||| ||||||||||||||
1996 3 4421.8 |||||||||| |||||||||| ||||||| ||||||||
1996 4 4421.6 ||||||||| ||||||||| |||||||||
1996 5 4420.6 |||||||| ||||||| |||||||||
1996 6 4419.2 |||||| |||| ||||||||
1996 7 4417.9 |||| || |||||
1996 8 4417.3 ||| || |
1996 9 4416.5 || | ||
1996 10 4416.4 || |
1996 11 4416.7 || |
1996 12 4416.4 || | |||
1997 1 4419.0 ||||| |||| ||
1997 2 4419.6 |||||| ||||| |||||
1997 3 4420.1 ||||||| |||||| |||||||
1997 4 4419.3 |||||| |||| |||||||||
1997 5 4417.4 ||| || |||||||||
1997 6 4415.0 | |||||||
1997 7 4415.0 |||||
1997 8 4415.0 |
1997 9 4415.0 ||
1997 10 4415.3
1997 11 4415.8 |
1997 12 4415.5 |||
1998 1 4417.5 ||| || ||
1998 2 4418.7 ||||| |||| |||||
1998 3 4418.5 ||||| ||| |||||||
1998 4 4416.8 || | |||||||||
1998 5 4415.0 ||| ||||||
1998 6 4415.0 |||||||
1998 7 4415.0 |||||
1998 8 4415.0 |
1998 9 4415.0 ||
1998 10 4415.0
1998 11 4415.0
1998 12 4415.0 |||
1999 1 4415.0 ||
1999 2 4420.1 ||||||| |||||| |||||
1999 3 4421.8 |||||||||| |||||||||| ||||||| |||
1999 4 4421.8 |||||||||| |||||||||| ||||||||| ||||||
1999 5 4421.7 ||||||||| ||||||||| |||||||||
1999 6 4421.0 |||||||| |||||||| ||||||||
1999 7 4420.2 ||||||| |||||| |||||
1999 8 4419.9 ||||||| |||||| |
1999 9 4419.4 |||||| ||||| ||
1999 10 4419.6 |||||| |||||
1999 11 4419.9 ||||||| ||||||
1999 12 4421.8 |||||||||| |||||||||| |||
2000 1 4421.8 |||||||||| |||||||||| || |||||||||||||||||||||||
2000 2 4421.8 |||||||||| |||||||||| ||||| |||||||||||||||||||||||||||||||||||||||||||||
2000 3 4421.8 |||||||||| |||||||||| ||||||| |||||||||||||||||||||||||||||||
2000 4 4421.8 |||||||||| |||||||||| ||||||||| ||||||
2000 5 4421.3 ||||||||| |||||||| |||||||||
2000 6 4420.3 ||||||| |||||| ||||||||
2000 7 4419.4 |||||| ||||| |||||
2000 8 4419.0 ||||| |||| |
‐ 3 ‐
Nivel de Nivel de Volumen Vol entreg. para
año mes i Déficit Vri‐Vdi Volumen vertido Vvi
Embalse Emb. (%) disp. (%) riego Vri
2000 9 4418.4 ||||| ||| ||
2000 10 4419.4 |||||| |||||
2000 11 4420.0 ||||||| ||||||
2000 12 4420.0 ||||||| |||||| |||
2001 1 4421.2 ||||||||| |||||||| ||
2001 2 4421.0 |||||||| |||||||| |||||
2001 3 4420.7 |||||||| ||||||| |||||||
2001 4 4419.5 |||||| ||||| |||||||||
2001 5 4417.4 ||| || |||||||||
2001 6 4415.0 | |||||||
2001 7 4415.0 |||||
2001 8 4415.0 |
2001 9 4415.0 ||
2001 10 4415.4
2001 11 4415.8 |
2001 12 4415.4 |||
2002 1 4415.0 | |
2002 2 4418.0 |||| ||| |||||
2002 3 4419.9 ||||||| |||||| |||||||
2002 4 4420.0 ||||||| |||||| |||||||||
2002 5 4419.2 |||||| |||| |||||||||
2002 6 4417.7 ||| || ||||||||
2002 7 4416.1 | |||||
2002 8 4415.5 |
2002 9 4415.0 | |
2002 10 4415.1
2002 11 4415.3
2002 12 4415.0 || |
2003 1 4415.0 |
2003 2 4417.7 ||| || |||||
2003 3 4420.0 ||||||| |||||| |||||||
2003 4 4419.9 ||||||| |||||| |||||||||
2003 5 4418.8 ||||| |||| |||||||||
2003 6 4417.0 ||| | ||||||||
2003 7 4415.1 |||||
2003 8 4415.0 |
2003 9 4415.0 ||
2003 10 4415.0
2003 11 4415.0
2003 12 4416.1 | |||
2004 1 4417.3 ||| || ||
2004 2 4420.9 |||||||| ||||||| |||||
2004 3 4421.8 |||||||||| |||||||||| ||||||| ||||||||||||
2004 4 4421.8 |||||||||| |||||||||| ||||||||| ||
2004 5 4421.1 |||||||| |||||||| |||||||||
2004 6 4420.0 ||||||| |||||| ||||||||
2004 7 4418.9 ||||| |||| |||||
2004 8 4418.5 ||||| ||| |
2004 9 4418.8 ||||| |||| ||
2004 10 4419.5 |||||| |||||
2004 11 4419.9 ||||||| ||||||
2004 12 4421.4 ||||||||| ||||||||| |||
2005 1 4421.8 |||||||||| |||||||||| ||
2005 2 4421.0 |||||||| |||||||| |||||
2005 3 4421.6 ||||||||| ||||||||| |||||||
2005 4 4420.9 |||||||| |||||||| |||||||||
2005 5 4419.5 |||||| ||||| |||||||||
2005 6 4417.6 ||| || ||||||||
2005 7 4415.7 | |||||
2005 8 4415.0 |
2005 9 4415.0 ||
2005 10 4415.0
2005 11 4415.0
2005 12 4415.0 || |
2006 1 4415.0 ||
2006 2 4415.0 ||| ||
2006 3 4415.0 ||| ||||
2006 4 4415.0 |||||| |||
2006 5 4415.0 |||||||| |
2006 6 4415.0 |||||||
2006 7 4415.0 |||||
2006 8 4415.0 |
2006 9 4415.0 |
‐ 4 ‐
Nivel de Nivel de Volumen Vol entreg. para
año mes i Déficit Vri‐Vdi Volumen vertido Vvi
Embalse Emb. (%) disp. (%) riego Vri
2006 10 4415.3
2006 11 4415.9 |
2006 12 4416.3 | | |||
2007 1 4419.9 ||||||| |||||| ||
2007 2 4421.8 |||||||||| |||||||||| ||||| ||||||||||
2007 3 4421.8 |||||||||| |||||||||| ||||||| |||||||||||||||||||||
2007 4 4421.8 |||||||||| |||||||||| ||||||||| ||||||||
2007 5 4421.7 ||||||||| ||||||||| |||||||||
2007 6 4421.0 |||||||| |||||||| ||||||||
2007 7 4420.2 ||||||| |||||| |||||
2007 8 4419.9 ||||||| |||||| |
2007 9 4419.3 |||||| ||||| ||
2007 10 4419.3 |||||| |||||
2007 11 4419.4 |||||| |||||
2007 12 4419.8 ||||||| ||||| |||
2008 1 4421.8 |||||||||| |||||||||| || |||
2008 2 4421.8 |||||||||| |||||||||| ||||| ||
2008 3 4421.1 |||||||| |||||||| |||||||
2008 4 4419.4 |||||| ||||| |||||||||
2008 5 4417.0 || | |||||||||
2008 6 4415.0 ||| |||||
2008 7 4415.0 |||||
2008 8 4415.0 |
2008 9 4415.0 ||
2008 10 4415.4
2008 11 4415.9 |
2008 12 4417.8 |||| || |||
2009 1 4418.7 ||||| |||| ||
2009 2 4421.1 |||||||| |||||||| |||||
2009 3 4420.9 |||||||| ||||||| |||||||
2009 4 4419.6 |||||| ||||| |||||||||
2009 5 4417.6 ||| || |||||||||
2009 6 4415.0 |||||||
2009 7 4415.0 |||||
2009 8 4415.0 |
2009 9 4415.0 ||
2009 10 4415.0
2009 11 4417.5 ||| ||
2009 12 4418.7 ||||| |||| |||
2010 1 4420.2 ||||||| |||||| ||
2010 2 4421.5 ||||||||| ||||||||| |||||
2010 3 4421.8 |||||||||| |||||||||| ||||||| ||
2010 4 4421.4 ||||||||| ||||||||| |||||||||
2010 5 4420.3 ||||||| |||||| |||||||||
2010 6 4418.8 ||||| |||| ||||||||
2010 7 4417.4 ||| || |||||
2010 8 4416.8 || | |
2010 9 4415.8 | ||
2010 10 4415.7 |
2010 11 4415.6
2010 12 4416.6 || | |||
2011 1 4420.1 ||||||| |||||| ||
2011 2 4421.8 |||||||||| |||||||||| ||||| ||||||||||||||||||||||
2011 3 4421.8 |||||||||| |||||||||| ||||||| |||||||||||||||||||||
2011 4 4421.8 |||||||||| |||||||||| ||||||||| ||||||||||||
2011 5 4421.8 |||||||||| |||||||||| ||||||||| |
2011 6 4421.2 ||||||||| |||||||| ||||||||
2011 7 4420.5 |||||||| ||||||| |||||
2011 8 4420.2 ||||||| |||||| |
2011 9 4419.8 ||||||| ||||| ||
2011 10 4419.8 ||||||| |||||
2011 11 4420.3 ||||||| ||||||
2011 12 4420.9 |||||||| |||||||| |||
2012 1 4421.8 |||||||||| |||||||||| || ||
2012 2 4421.8 |||||||||| |||||||||| ||||| |||||||||||||||||||
2012 3 4421.8 |||||||||| |||||||||| ||||||| |||||||||
2012 4 4421.8 |||||||||| |||||||||| ||||||||| |
2012 5 4421.1 |||||||| |||||||| |||||||||
2012 6 4420.0 ||||||| |||||| ||||||||
2012 7 4418.9 ||||| |||| |||||
2012 8 4418.5 ||||| ||| |
‐ 5 ‐
Nivel de Nivel de Volumen Vol entreg. para
año mes i Déficit Vri‐Vdi Volumen vertido Vvi
Embalse Emb. (%) disp. (%) riego Vri
2012 9 4417.9 |||| || ||
2012 10 4417.9 |||| ||
2012 11 4417.8 |||| ||
2012 12 4418.8 ||||| |||| |||
2013 1 4419.2 |||||| |||| ||
2013 2 4419.0 ||||| |||| |||||
2013 3 4417.1 ||| | |||||||
2013 4 4415.0 |||| |||||
2013 5 4415.0 |||||||||
2013 6 4415.0 ||||||||
2013 7 4415.0 |||||
2013 8 4415.0 |
2013 9 4415.0 ||
2013 10 4415.0
2013 11 4415.3
2013 12 4415.0 |||
2014 1 4417.0 || | ||
2014 2 4418.5 ||||| ||| |||||
2014 3 4421.8 |||||||||| |||||||||| ||||||| |
2014 4 4421.8 |||||||||| |||||||||| ||||||||| ||
2014 5 4421.1 ||||||||| |||||||| |||||||||
2014 6 4420.0 ||||||| |||||| ||||||||
2014 7 4418.9 ||||| |||| |||||
2014 8 4418.5 ||||| ||| |
‐ 6 ‐