2.3.1.

Oferta de agua

El objeto del presente Ítem es la estimación de los caudales naturales o escorrentía superficial
natural en los puntos de aprovechamiento involucrados en el esquema hidráulico del proyecto y
que corresponden básicamente al Proyecto “ELABORACIÓN DEL ESTUDIO DE PERFIL A
NIVEL DE FORMULACIÓN Y EVALUACIÓN DEL PROYECTO: "CREACIÓN DE FUENTE DE
AGUA POR EXTRACCIÓN Y RECARGA ARTIFICIAL INDUCIDA EN EL VALLE DEL RÍO
CHILLÓN DISTRITO DE CARABAYLLO - PROVINCIA DE LIMA - DEPARTAMENTO DE LIMA ”.
Como no se ha encontrado información hidrométrica alguna en la zona del proyecto y debido a
que el área de la cuenca de recolección es pequeña, se ha optado por utilizar un coeficiente
precipitación-escorrentía,
El agua, cuando consideramos como recurso, es necesario conocer las cantidades disponibles en
la cuenca de estudio. Las fuentes de agua disponibles se pueden dividir en superficiales: ríos,
lagos, embalses y subterráneas: acuíferos. Para este estudio hidrológico solo se tomó en cuenta
los cursos superficiales.
Para saber la cantidad de agua disponible en una cuenca se realiza el llamado balance hídrico,
que consiste en contabilizar todas las entrada y salidas de agua de la cuenca.
2.3.2 ANÁLISIS Y TRATAMIENTO DE LA ESCORRENTÍA SUPERFICIAL

Los recursos hídricos en los afluentes en el río Chillón, considerados para usos poblacional, y
agrícola, están constituidos por la escorrentía superficial que discurre por su cauce.
Se ha utilizado información de descargas mensuales y diarias en el río Chillón, para generar,
analizar y calibrar o ajustar caudales medios mensuales en el Chillón, en el punto de interés
(Regulación y Derivación o Sección de Control de Captación).
La metodología referida, puede aplicarse por la similitud de las características de las unidades
hidrográficas en lo relacionado a fisiografía, pisos altitudinales, geomorfología y cobertura
vegetal, aunque las áreas sean diferentes.

2.3.3. MODELO DETERMINISTICO – ESTOCÁSTICO LUTZ SCHOLZ

Actualmente sabemos que el agua es un recurso finito, pues la oferta es limitada y demanda va
en aumento en el tiempo; por esta razón para usar el agua racionalmente en las cuencas se debe
estimar las ofertas y demandas de aguas utilizando teorías adecuadas que sean las más
aproximadas a la realidad.
Bajo este principio, para estimar la oferta existen diferentes métodos: La hidrología estocástica
es uno de estos métodos, que se ha aplicado en diferentes partes del Perú y del mundo con
resultados satisfactorios. La modelación estocástica en la hidrología es una herramienta que
evalúa el comportamiento espacial y temporal del agua utilizando los registros de caudales
anuales; estacionales (mensuales) de varios años; por tanto, el análisis de las series de tiempo
hidrológico se realiza sobre base de los registros que ya han sucedido en el tiempo, y estas
sirven para proyectar el comportamiento en el futuro.
En la actualidad, la gestión del recurso hídrico en el Perú exige la evaluación de la OFERTA
HÍDRICA de la cuenca colectora y de la demanda multisectorial del agua. La evaluación del
recurso hídrico involucra el análisis e interpretación del comportamiento espacial y temporal de
los componentes del ciclo hidrológico, especialmente de las descargas de las aguas
superficiales.
Para pronosticar el comportamiento futuro de la oferta de agua es necesario conocer el
comportamiento histórico cuantitativo de este recurso, para lo cual se necesita de datos
históricos de aforo de varios años. En los diferentes puntos de la cuenca del río Chillón se
cuenta con datos de varios años, lo cual es una herramienta útil para formular modelos
estocásticos que permitirían pronosticar la OFERTA DE AGUA para varios años en cualquier
punto.
El primer aporte del presente estudio consiste en que es una herramienta para estimar la
OFERTA DE AGUA en puntos no aforados de la cuenca, mediante modelos estocásticos
regionales y estos caudales estimados servirían para realizar planificación hidrológica e
hidráulica.
El segundo aporte consiste en que se puede estimar las descargas anuales o mensuales para
períodos mayores a 25 años y con estos datos generados se puede determinar la oferta para
diferentes niveles de persistencia y para diferentes usos.

El objetivo del presente estudio es determinar modelos estocásticos que permitan generar
descargas medias anuales o descargas medias mensuales en sitios no aforados en la cuenca del
río Chillón y estos servirán para una planificación hidrológica
Este modelo hidrológico es combinado porque cuenta con una estructura determinística para el
cálculo de los caudales medios mensuales para el año promedio (Balance Hídrico – Modelo
determinístico) y una estructura estocástica para la generación de series extendidas de caudal
(Proceso Markoviano – Modelo Estocástico). Fue desarrollado por experto en hidrología, Lutz
Sholtz para cuencas de la sierra peruana, entre los años 1979 – 1980, en el marco de
Cooperación Técnica de la República de Alemania a través del Plan Meris II.
Determinado la situación de ausencia de registros de caudal en la Sierra Peruana, el modelo se
desarrolló tomando en consideración parámetros físicos y meteorológicos de las cuencas, que
puedan ser obtenidos a través de mediciones cartográficas y de campo. Los parámetros más
importantes del modelo son los coeficientes para la determinación de la precipitación efectiva,
déficit de escurrimiento, retención y agotamiento de las unidades hidrográficas. Los
procedimientos que se han seguido en la implementación del modelo son:
 Computo de los parámetros necesarios para la descripción de los fenómenos de escorrentía
promedio.
 Establecimiento de un conjunto de modelos parciales de los parámetros para el cómputo de
los caudales en unidades hidrográficas sin información hidrométrica. Con base en lo
anterior se realiza el cálculo de los caudales necesarios.
 Calibración del modelo y generación de caudales extendidos por un proceso Markoviano
combinado de precipitación efectiva del mes con el caudal del mes anterior.
Este modelo fue implementado con fines de pronosticar caudales a escala mensual,
teniendo una utilización inicial en estudios de proyectos de riegos y posteriormente
extendiéndose el uso del mismo a estudios hidrológicos con prácticamente cualquier
finalidad (abastecimiento del agua, hidroelectricidad, presas, otros). Los resultados de la
aplicación del modelo a las cuencas hidrográficas de la Sierra Peruana, han producido
una correspondencia satisfactoria respecto a los valores medidos en los puntos de
interés. Este método se aplicó para microcuenca del río Iruro.

a. Ecuación del balance hídrico

La ecuación fundamental que describe el balance hídrico mensual en mm/mes es la siguiente:

Qmi = Pi – Di + Ci – Ai

Donde:

Qi : Caudal mensual en mm/mes.

Pi : Precipitación total mensual sobre la cuenca en mm/mes

Di : Déficit de escurrimiento en mm/mes

Ci : Gasto de la retención en la cuenca (Aporte)

Ai : Abastecimiento de la retención (Retención).

Asumiendo que para períodos largos (en este caso 1 año) el Gasto y Abastecimiento de la
retención tienen el mismo valor, es decir G = A y que para el año promedio, una parte de la
precipitación retorna a la atmósfera por evaporación; luego reemplazando (P-D) por (C*P), y
tomando en cuenta la transformación de unidades (mm/mes a m3/s) la ecuación se convierte en:

Q = c'*C * P * AR

Que es la expresión básica del método racional, donde:

Q = Caudal (m3/s)

c1 = Coeficiente de conversión del tiempo (mes/s)

C = Coeficiente de escurrimiento

P = Precipitación total mensual (mm/mes)

AR = Área de la cuenca (m2).

b) Coeficiente de escurrimiento

Se ha considerado el uso de la fórmula propuesta por L. Turc:

P−D
C=
P
Donde:

C = Coeficiente de escurrimiento (mm/año)

P = Precipitación Total anual (mm/año)

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

Para la determinación de D se utiliza la expresión:

P
D=

( )
1/ 2
P2
0 . 9+ 2
L

Donde:

L = Coeficiente de Temperatura

T = Temperatura media anual (°C)

L = 300 + 25(T) + 0.05 (T)3

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

C=3.16 E 12( P−0.571 )( EP− 3.686 );r=0.96

D=−1380+0.872(P )+1.032( EP);r=0.96
Donde;

C = Coeficiente de escurrimiento

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

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

EP= Evapotranspiración anual según Hargreaves (mm/año)

r = Coeficiente de correlación

La evapotranspiración potencial, se ha determinado por MÉTODO DE THORNTHWAITE

T
ETP=16∗(10∗ )a
I
ETP=f∗ETP '
Evapotranspiración potencial mensual (mm), por mes de 30
ETP´ : días y 12 horas de duración

ETP : Evapotranspiración potencial mensual corregida (mm)

T : Temperatura media mensual (ºC)

I : Índice térmico anual

i : Índice térmico mensual

Exponente que varía con el índice anual de calor de la

a : localidad

f : Factor de corrección

X|

ETP=0 . 003∗( RS )2. 5 +0 . 16∗(T )0 .88

 ETP : Evapotranspiración potencial

RS : Radiación solar extraterrestre, expresada en equivalentes

de evaporación

T : Temperatura media mensual

c) Precipitación efectiva (PE)

Para el cálculo de la Precipitación Efectiva, se supone que los caudales promedios observados
en la cuenca pertenecen a un estado de equilibrio entre gasto y abastecimiento de la retención.
La precipitación efectiva se calculó para el coeficiente de escurrimiento promedio, de tal forma
que la relación entre precipitación efectiva y precipitación total resulta igual al coeficiente de
escorrentía.

A fin de facilitar el cálculo de la precipitación efectiva se ha determinado el polinomio de

quinto grado:

PE=a0 + a1 P+a 2 P 2 +a 3 P3 + a4 P 4 + a5 P5

Donde:

PE = Precipitación efectiva (mm/mes)

P = Precipitación total mensual (mm/mes)

ai = Coeficiente del polinomio

Se muestra los valores límite de la precipitación efectiva y el cuadro 8 muestra los tres juegos de
coeficientes para, ai que permiten alcanzar por interpolación.

Límite superior para la Precipitación Efectiva.

Curva I: PE = P- 120.6 para P > 177.8 mm/mes

Curva II: PE = P-86.4 para P> 152.4 mm/mes

Curva III: PE = P -59.7 para P> 127.0 mm /mes

Cuadro 8: Coeficientes para el Cálculo de la Precipitación Efectiva

Valor del Coeficiente

Descripción
Curva I Curva II Curva III

ao -0.01800 -0.02130 -0.0280

 a1 -0.01850 0.13580 0.2756

a2 0.001105 -0.002296 -0.004103

a3 -1.20E-05 4.35E-05 5.53E-05

a4 1.44E-07 -8.90E-08 1.24E-07

a5 -2.85E-10 -8.79E-11 -1.42E-09

De esta forma es posible llegar a la relación entre la precipitación efectiva y precipitación total:

12 PE
Q
C= =∑ i
P i=1 P

Donde:

C = Coeficiente de escurrimiento

Q = Caudal anual

P = Precipitación Total anual

12 PEi
∑ P = Suma de la precipitación efectiva mensual
i=1

d.Retención de la cuenca

Bajo la suposición de que para un año promedio exista un equilibrio entre el gasto y el
abastecimiento de la reserva de la cuenca y además que el caudal total sea igual a la
precipitación efectiva anual, la contribución de la reserva hídrica al caudal se puede calcular
según las fórmulas:

Ri=CM i −Pi
CM i=PEi +G i − Ai

Donde:

CMi = Caudal mensual (mm/mes)

PEi = Precipitación Efectiva Mensual (mm/mes)

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

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

A¡ = Abastecimiento de la retención (mm/mes)

 R¡ = Gi para valores mayores que cero (mm/mes)

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

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

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

e. Relación entre descargas y retención

Durante la estación seca, el gasto de la retención alimenta los ríos, constituyendo el caudal o
descarga básica. La reserva o retención de la cuenca se agota al final de la estación seca; durante
esta estación la descarga se puede calcular en base a la ecuación:

Qt =Q0 e−a t

Donde:

Qt : Descarga en el tiempo t.

Q0 : Descarga inicial.

a : Coeficiente de agotamiento.

t : Tiempo.

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

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

Coeficiente de Agotamiento

Mediante la fórmula anterior se puede calcular el coeficiente de agotamiento "a", en base a

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

Con fines prácticos se puede despreciar la variación del coeficiente "a" durante la estación seca
empleando un valor promedio.

El coeficiente de agotamiento de la cuenca tiene una dependencia logarítmica del área de la

cuenca.

a=f ( Ln A )
El análisis de las observaciones disponibles muestra, además, cierta influencia del clima, de la
geología y de la cobertura vegetal. Se ha desarrollado una ecuación empírica para la sierra
peruana:

a=3.1249 x10 67 A−0 .144 EP−19 . 336 T −3 . 369 R−1. 429 ; r=0.86
Donde:

a : Coeficiente de agotamiento por día.

A : Área de la cuenca en Km2.

EP : Evapotranspiración potencial anual en mm/año.

T : duración de la temporada seca en días.

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

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

en el río durante la estación seca; sin embargo cuando no sea posible ello, se puede recurrir a las
ecuaciones desarrolladas para la determinación del coeficiente "a" para cuatro clases de cuen -
cas:

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

a=−0.00252 Ln A+0.034
Cuencas con agotamiento rápido, cuya retención varía entre 50 y 80 mm/año y vegetación poco
desarrollada (puna).

a=−0.00252 Ln A+0.030
Cuencas con agotamiento mediano, cuya retención es alrededor de 80 mm/año y vegetación
mezclada (pastos, bosques y terrenos cultivados).

a=−0.00252 Ln A+0.026
Cuencas con agotamiento reducido por alta retención (>100mm/año) y vegetación mezclada.

a=−0.00252 Ln A+0.023
Donde:

a = coeficiente dé agotamiento por día

AR = área de la cuenca (km2)

EP = evapotranspiración potencia! anual (mm/año)

T = duración de la temporada seca (días)

R = retención total de la cuenca (mm/año)

Almacenamiento Hídrico

Tres tipos de almacenes hídricos naturales que inciden en la retención de la cuenca son
considerados:

• Acuíferos

• Lagunas y pantanos

• Nevados

La determinación de la lámina "L" que almacena cada tipo de estos almacenes está dado por:

• Acuíferos

Siendo:

L A=−750 I +3 . 15 (mm /año )

LA = lámina específica de acuíferos

I = pendiente de desagüe: / <= 15 %

• Lagunas y Pantanos

LL=500 (mm/año).

Siendo:

LL = Lámina específica de lagunas y pantanos

• Nevados

LN = 500 (mm/año)

Siendo:

LN = lámina específica de nevados

Las respectivas extensiones o áreas son determinadas de los mapas o aerofotografías. Los
almacenamientos de corto plazo no son considerados para este caso, estando los mismos
incluidos en las ecuaciones de la precipitación efectiva.

Abastecimiento de la Retención

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

región climática. En la región del Cuzco el abastecimiento comienza en el mes de noviembre
con 5%, alcanzando hasta enero el valor del 80 % del volumen final. Las precipitaciones altas
del mes de febrero completan el 20% restante, y las precipitaciones efectivas del mes de marzo
escurren directamente sin contribuir a la retención. Los coeficientes mensuales expresados en
porcentaje del almacenamiento total anual se muestran en el cuadro siguiente.

Cuadro 9: Almacenamiento hídrico durante la época de lluvias - (valores –a i %)

Mes Cuzco Huancavelica Junín Cajamarca

1 5 3 4

Ene 40.0 30.0 30.0 20.0

Feb 20.0 20.0 30.0 25.0

Mar - 5.0 5.0 35.0

Abr

May

Jun

Jul

Ago

Sep

Oct - 10.0 10.0 25.0

Nov 5.0 - - - 5.0

Dic 35.0 35.0 25.0 -

Total 100.0 100.0 100.0 100.0

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

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

Ai =ai ( )
R
100

Siendo:
Ai= abastecimiento mensual déficit de la precipitación efectiva (mm/mes)

ai= coeficiente de abastecimiento (%)

R = retención de la cuenca (mm/año)

DETERMINACIÓN DEL CAUDAL MENSUAL PARA EL AÑO PROMEDIO

a. Caudales medios mensuales para el año promedio río Iruro
Previamente se ha procesado y analizado la información hidrométrica de la cuenca hidrográfica
vecina del río Chillón, con el cálculo de la proporción de lluvia que genera escorrentía, los
coeficientes para el computo de la precipitación efectiva, los coeficientes anuales y mensuales
para el cálculo de la contribución de la retención, la evapotranspiración potencial mensual y
determinación de caudales medios mensuales por el Método Lutz – Sholtz
Igualmente, la determinación del caudal mensual para el año promedio en la cuenca del río
Chillón en los puntos de interés, está también basado en la ecuación fundamental que describe
el balance hídrico mensual hídrico a partir de los componentes descritos anteriormente:
CMi = PEi + Gi - Ai
Donde:
CMi = Caudal del mes í (mm/mes)
PEi = Precipitación efectiva del mes i (mm/mes)
Gi = Gasto de la retención del mes i (mm/mes)
Ai = Abastecimiento del mes i (mm/mes)
a. Ingreso de Datos
Los datos de ingreso son la Precipitación Áreal (para el presente estudio obtenido por la
metodología de regionalización), los porcentajes mensuales de almacenamiento del acuífero
(retención), los aportes del acuífero y los caudales de comparación.
Los caudales aforados de la estación Puente Magdalena, las retenciones y aportes serán
obtenidos por calibración.
b. Calibración del modelo
Consiste en el ajuste de los parámetros del modelo durante el período fijado, por comparación
entre valores simulados y los caudales reales medidos en estaciones de aforo en el punto de
interés. Es la comprobación de la bondad del modelo.
Considerando la estación de aforo más cercana (Estación Puente Magdalena), se ha visto por
conveniente realizar la calibración en este punto.
Con estos valores se elige el coeficiente C y la retención de las unidades hidrográficas en
estudio.
c. Simulación
Una vez realizada la calibración, podemos realizar la simulación de caudales en el punto de
interés, utilizando los valores de los parámetros ajustados previamente.
La simulación, para nuestro caso fue en el punto de almacenamiento de las aguas del río Chillón
con área de cuenca total de 1999.9 km2 y tiene una altitud media de 2682 m.s.n.m.
Como resultado de la Simulación, se obtienen los caudales medios mensuales, para el período
1987 – 2022,.

Generación de Caudales Mensuales para Períodos Extendidos

A fin de generar una serie sintética de caudales para períodos extendidos, se ha implementado
un modelo estocástico que consiste en una combinación de un proceso Markoviano de primer
orden, según la ecuación siguiente, con una variable de impulso, que en este caso es la
precipitación efectiva en la ecuación siguiente:
Qti = f (Qt-1)
Q = g (PEt)
Con la finalidad de aumentar el rango de Valores generados y obtener una óptima aproximación
a la realidad, se utiliza además una variable aleatoria.
Z = z(S) (1 – r2)
La ecuación integral para la generación de caudales mensuales es:
Qt =B1 + B2 (Qt-1) + B3 (PEt) + z(Z) (1 – r2)0.5
Donde:
Qt = Caudal del mes t
Qt-1 = Caudal del mes anterior
PEt = Precipitación efectiva del mes
B1 = Factor constante o caudal básico.
Se calcula los parámetros B1, B2, B3, r y S sobre la base de los resultados del modelo para el
año promedio por un cómputo de regresión lineal con Q t como valor dependiente y Q t-1 y PEt
como valores independientes. Para el cálculo se recomienda el uso de software comercial (hojas
electrónicas) o de uso específico (programas elaborados tales como el SIH).
El proceso de generación requiere de un valor inicial, el cual puede ser obtenido en una de las
siguientes formas:
 Empezar el cómputo en el mes para el cual se dispone de un aforo.
 Tomar como valor inicial el caudal promedio de cualquier mes.
 Empezar con un caudal cero, calcular un año y tomar el último valor como valor Q 0 sin
considerar estos valores en el cálculo de los parámetros estadísticos del período generado.
Con los datos generados podemos obtener el caudal de diseño al 75% de los ríos en estudio los
caudales generados se presentan en el Anexo.

d. Restricciones del Modelo

El modelo presenta ciertas restricciones de uso o aplicación tales como:
 El uso de los modelos parciales, únicamente dentro del rango de calibración establecido.
 Su uso es únicamente para el cálculo de caudales mensuales promedio.
 Los registros generados en el período de secas presentan una mayor confiabilidad que los
valores generados para la época lluviosa.
 La Aplicación del modelo se restringe a las cuencas en las que se ha calibrado sus
parámetros (Sierra Peruana: Apurímac, Ayacucho, Cusco, Huancavelica, Junín, Cajamarca).
Pero puede ser adaptado por similitud de cuencas y periodos de precipitación
Es importante tener en cuenta las mencionadas restricciones a fin de garantizar una buena
performance del modelo.
Resultados
Cuenca del Río Chillón
Como se mencionó en apartados anteriores en la cuenca Chillón se proyecta la defensa rivereña
con enrocados y pantallas para recargas inducidas, respecto a su extracción se realizara a partir
de pozos ubicados a través de un estudio hidrogeológico y geofísico que nos servirá para
determinar la ubicación de los pozos propuestos .
Se ha analizado, evaluado el modelo que ha permitido la determinación de caudales medios para
la cuenca en estudio. Obteniendo resultados presentados a continuación.