MODELO COMBINADO ENTRE EL BALANCE HIDRICO Y UN

PROCESO MARKOVIANO DE PRIMER ORDEN

(ING. ABELARDO VELASQUEZ ALPACA) SETIEMBRE 2,010

 GENERACION DE CAUDALES EN LA SIERRA PERUANA
RESUMEN HANOCCA - TAYPITUNGA
La necesidad de contar con caudales medios mensuales en
pequeñas y medianas cuencas en la sierra peruana que nos
sirvan como recursos hídricos para elaborar proyectos de
irrigación es posible determinar utilizando el modelo
combinado entre el balance hídrico con un proceso
Markoviano de primer orden.
En un primer paso, se ha analizado la información
hidrometeorológica de 19 cuencas de Cusco y Cajamarca,
donde se obtuvieron los parámetros necesarios para el
análisis de los fenómenos de la escorrentía promedio.
En un segundo paso, se ha establecido un conjunto de
modelos estocásticos parciales de los parámetros para el
cálculo de caudales en estas cuencas que carecían de
información hidrométrica.
 GENERACION DE CAUDALES EN LA SIERRA PERUANA
HANOCCA - TAYPITUNGA
Empleando los datos meteorológicos regionalizados de la
cuenca y los modelos parciales, se pueden calcular los
caudales mensuales para el año promedio, todos ellos
generados mediante el balance hídrico.
Un tercer paso nos permite generar los caudales medios
mensuales para periodos extendidos mediante el cálculo
combinando entre la precipitación efectiva con las descargas
del mes anterior, mediante un proceso Markoviano I,
calibrando el modelo integral con los aforos realizados en la
estación de Pisac.
Posteriormente a los caudales generados utilizando el
modelo se les ha realizado un test estadístico de
comparación con los registros hidrométricos aforados, de
cuyos resultados se constato la existencia de una
correspondencia satisfactoria.
 GENERACION DE CAUDALES EN LA SIERRA PERUANA
HANOCCA
FINALIDAD DEL MODELO - TAYPITUNGA
E INFORMACION DISPONIBLE
1. La mayoría de
CARENCIA DE INFORMACION HIDROMETRICA.-
cuencas pequeñas y medianas en nuestra región no cuentan
con información hidrométrica, la necesidad de contar con
caudales medios mensuales que nos permitan contar como
recursos hídricos para elaborar proyectos de riego, es
posible determinar los caudales utilizando el modelo
combinado entre el balance hídrico con un proceso
Markoviano I.
2. GENERALIDADES DEL MODELO.- Para disponer de caudales en
cuencas, el modelo analiza la precipitación mensual
mediante el balance hídrico, considerando las características
físicas de la cuenca. En base a los datos disponibles se
calibraron varios modelos parciales utilizando variables
como: precipitación efectiva, déficit de escurrimiento,
retención de la cuenca y coeficiente de agotamiento.
 GENERACION DE CAUDALES EN LA SIERRA PERUANA
HANOCCA - TAYPITUNGA
RESTRICCIONES DEL MODELO
Se debería descuidar varios factores en el régimen de
escorrentía, por ejemplo, el proceso de la intercepción, la
variación anual del coeficiente de escurrimiento y del
coeficiente de agotamiento, la reacción específica de los
almacenes hídricos durante la estación lluviosa y la sequía
para llegar a un método simple y practico. Esta simplificación
implica algunas restricciones en el uso del modelo a
presentarse. La primera y la mas importante es su uso
únicamente dentro del rango de la calibración de los modelos
parciales. La segunda es que es permisible emplearlo
solamente para el cálculo de caudales mensuales.
DATOS DISPONIBLES
En la parte hidrológica se conto con los caudales en el
Río Vilcanota registrados entre 1965 y 1976 en la estación
de Pisac que controla una superficie de 7,145 Km2.
 GENERACION DE CAUDALES EN LA SIERRA PERUANA
HANOCCA - TAYPITUNGA
En la misma cuenca del Río Vilcanota se ha registrado
caudales de los Río Pitumarca, Pomacanchi y Tungasuca
durante un solo año. Fuera de la cuenca existe un registro de
los caudales del Río Paucartambo entre 1965 y 1972.
Además, existen aforos aislados en la cuenca ejecutados a
partir del mes de agosto de 1978 y algunos aforos
esporádicos ejecutados anteriormente.
En la parte meteorológica se conto con los registros
de la precipitación total mensual, temperatura media y
humedad relativa para cinco estaciones ubicadas en la
cuenca alta del Vilcanota aguas arriba de Pisac y otras 15
estaciones más situadas alrededor de la cuenca misma. Se
ha generado los datos principales como valores mensuales
para un periodo uniforme de 20 años desde 1957 hasta
1976.
 GENERACION DE CAUDALES EN LA SIERRA PERUANA
HANOCCA - TAYPITUNGA
SISTESIS DEL MODELO
La generación de caudales mensuales mediante el
modelo hidrológico para el año promedio se lleva a cabo
combinando los factores de la ecuación fundamental del
balance hídrico:
CMi = Pi – Di + Gi - Ai
Donde:
CMi = Caudal mensual (mm/mes)
Pi = Precipitación total mensual sobre la cuenca (mm/mes)
Di = Déficit de escurrimiento (mm/mes)
Gi = Gasto de la retención en la cuenca (mm/mes)
Ai = Abastecimiento de la retención (mm/mes)
Posterior a la explicación se presentaran notas
prácticas para la preparación de los datos básicos y la
aplicación del conjunto de modelos parciales.
 GENERACION DE CAUDALES EN LA SIERRA PERUANA
HANOCCA
PRECIPITACION SOBRE LA CUENCA- TAYPITUNGA
Se determino la precipitación total mensual sobre la
cuenca utilizando el método de Thiessen por polígonos de
influencia de las estaciones meteorológicas dentro y
alrededor de la cuenca, teniendo en cuenta la gradiente de la
precipitación calculada a partir de datos regionalizados.
COEFICIENTE DE ESCURRIMIENTO
Para la estimación del coeficiente del escurrimiento C,
en la cuenca del río Vilcanota se utiliza las ecuaciones
siguientes que fueron desarrolladas sobre la base de
observaciones en la región de Cusco y Huancavelica.
Ĉ = 3.16 E12 * P. -.571 * EP -3.686 ; r = 0.96
= -1380 + 0.872 * P + 1.032 * EP ; r = 0.96
Donde:
MAPA TOPOGRAFICO (BASE) DE LA CUENCA DEL RÍO MALA

M A P A D E U B IC A C IO N UBICACIÓN DE SECCIONES TRANSVERSALES

DE LA CUENCA DEL RÍO MALA
N
C
CO
L OM
O R BIA
U AD W E
EC

S
B

MALA La g. Ch um picoch a

C
BR A SIL La g. Tot oral

A l

a
Ma
O

ío
B
C

R
E
A
N
O

Santiago de Anchucaya
P
A

%
C
IF
IC

%
San Juan de T antar anche
O

San Pedro
BO LI VIA %
Huar ochiri % % San Lorenzo de Quinti
% A
Sangallaya

CHIL E

Mar iatana
%

00
30
% San Joaquin
Cochas % %
Huañec
%
Quinches
0 La g. Huascaco ch a
200

Quinocay Huampar a%
%

% Ayaviri
4 00 0

0
100
Calango Mal a
O %
C MI NI STER IO DE AGR ICU LTU RA
E INS TIT UTO NA CION AL DE R ECU RS OS N ATU RA LES
A
N D IR ECC IÓN GENER AL D E A GUA S Y SUE LOS
O ADMINISTRACIÓN TÉCNICA DEL DISTRITO DE RIEGO MALA-OMAS-CAÑETE
P PROYEC TO :
LE Y E NDA "EVALUAC IÓN DE LOS R ECU RSOS HÍDR IC OS EN LA CU ENC A D EL R ÍO MALA
A

ESTU DIO H ID ROLÓGICO E INV ENTARIO D E FU ENTES D E AGUA SUPER FIC IALES"
C
IF

Curvas de niv el_500 PLANO :

% Santa Cr uz de Flores "MAPA T OPO G RÁF ICO
IC

Rio princ ipal DE LA CUENC A DEL RÍO MALA"

O

% San Antonio
%
Nevados FE CHA : ESCALA : 1:6 00 000 PLANO N°
6000 0 6000 Meter s
NOVIEM BR E DEL 20 0 6

%
Lagunas
Capital de Dis trit o
FUENTE :
Cartas N acionales 1:100,000 IGN - Año 2002
DAT UM : WGS 84 / Proyección UTM - Zona 18.
08
 GENERACION DE CAUDALES EN LA SIERRA PERUANA
HANOCCA - TAYPITUNGA
PRECIPITACION EFECTIVA
Suponiendo que los caudales promedios observados
pertenezcan a un estado de equilibrio entre gasto y el
abastecimiento de la retención de la cuenca estudiada, la
precipitación efectiva se calculó para que el coeficiente de
escurrimiento promedio este de manera que la relación entre
precipitación efectiva y total sea igual al coeficiente de
escurrimiento. Para este cálculo se adopto el método del
United States Bureauof Reclamation (USBR) para la
determinación de la porción de lluvias que es aprovechable
para cultivos.
PEi = C * P
Donde:
PEi = precipitación efectiva mensual ( mm/mes )
P = precipitación total anual ( mm/años)
C = coeficiente de escurrimiento ( =1 )
 GENERACION DE CAUDALES EN LA SIERRA PERUANA
HANOCCA - TAYPITUNGA
El elemento constitutivo del método del USBR para el cálculo
de la precipitación efectiva para cultivos es el principio que
cuando aumenta la precipitación total mensual se toma un
porcentaje disminuyendo del incremento de la lluvia con
aumento de la precipitación efectiva para los cultivos se
mantenga constante.
Para la Hidrología se toma como precipitación efectiva esta
parte de la precipitación total mensual que sale como el
deficit según el método original del USBR. Para facilitar el
cálculo de la precipitación efectiva se ha determinado un
polinomio de quinto grado:
PE = a0 + a1 * P + a2 * P2 + a3 * P3 + a4 * P4 + a5 * P5
Donde:
PE = Precipitación efectiva (mm/mes)
P = Precipitación total mensual (mm/mes)
a1 = Coeficientes del Polinomio
 GENERACION DE CAUDALES EN LA SIERRA PERUANA
HANOCCA - TAYPITUNGA
En el cuadro siguiente se presentan tres juegos de los
coeficientes a1 que permiten alcanzar cada valor C
predeterminado entre 0.15 y 0.45 por interpolación.
COEFICIENTES PARA EL CALCULO SEGÚN
CONSTANTES CURVA I CURVA II CURVA III
a0 (-0.018) (-0.021) ( -0.028 )
a1 -0.0185 +0.1358 +0.2756
a2 +0.001105 -0.002296 -0.004103
a3 -1204 E – 8 + 4349 E - 8 + 5534 – 8
a4 +144 E – 9 - 89.0 E - 9 + 124 E – 9
a5 - 285 E - 12 - 879 E - 13 - 142 E - 11

RETENCION DE LA CUENCA
DIMENSION TOTAL
Se mide la extensión de los almacenes naturales –
acuíferos potenciales, lagunas, pantanos y nevados en el
mapa 1:100,000 y se multiplica su área por la lámina de
agua respectiva, según el cuadro siguiente
 GENERACION DE CAUDALES EN LA SIERRA PERUANA
HANOCCA - TAYPITUNGA
LAMINA ACUMULADA (mm / año)
TIPO
Pendiente de la cuenca

Napa 2% 8% 15%
Freática
300 250 200
Lagunas y
500
Pantanos
Nevados 500

La suma de los volúmenes particulares se transforma

a una lámina de agua sobre toda la cuenca. El resultado
indica la retención total de la cuenca para el año promedio
en (mm/año) que por definición es igual al gasto total G
durante la estación seca y además igual al abastecimiento A
durante la estación lluviosa.
El análisis ejecutado para 19 cuencas de la Sierra muestra
que la retención varía entre 43 y 188 (mm/año), de los
resultados se indican en el cuadro siguiente:
 GENERACION DE CAUDALES EN LA SIERRA PERUANA
HANOCCA - TAYPITUNGA
DATOS BASICOS PROVENIENTES DE LA SIERRA PERUANA
Lat. S. AR Q P C Ai T EP R Qo Qt t a
CUENCA
gº m´ Km2 mm/año mm/año (1) m.s.n.m. ºC mm/año m3/s mm/s mm/s días (1)
Vilcanota 14º 10´ 7,145 294 923 0.32 4,290 6.3 1134 80 92.4 25.1 150 0.0084
Paucartambo 13º 30´ 2,468 240 665 0.36 4,110 7.2 1184 57 20.5 8.7 120 0.007
Pitumarca 14º 00´ 694 403 919 0.44 4,650 4.4 1077 96 6.4 3.5 96 0.0063
Pomacanchi 14º 05´ 259 256 882 0.29 3,970 8 1218 90 2.4 0.55 150 0.0096
Sta. Elena 13º 02´ 595 370 916 0.4 4,370 5.8 1134 71 4.4 1.3 120 0.0102
Huapa 13º 00´ 965 311 902 0.34 4,300 6.2 1149 56 6.3 2.1 120 0.0092
Orcococha 13º 00´ 100 389 932 0.42 4,600 4.6 1087 188 1.4 1.1 210 0.0011
Quilión 12º 23´ 1,325 269 690 +) 0.39 4,300 6.2 1149 43 7.1 4.5 90 0.0051
Pachacayo 11º 49´ 720 412 1,168 0.35 4,525 5 1108 99 7.2 3.3 120 0.0065
Pinascocha 11º 51´ 195 250 899 0.28 4,470 5.3 1120 66 2.9 0.4 150 0.0132
Cochas 11º 45´ 490 398 1,162 0.34 4,600 4.6 1093 88 4.3 1.7 190 0.0103
Huari 11º 39´ 430 527 1,249 0.42 4,650 4.4 1086 129 6.6 1.6 120 0.0118
Casaracra 11º 07´ 370 170 610 0.28 4,255 6.4 1163 48 2.6 0.8 150 0.0079
Namora B 7º 16´ 461 368 968n 0.38 3,290 11.6 1364 82 5.2 1.3 90 0.0154
Matar 7º 14´ 140 407 963 0.42 3,260 11.7 1366 88 1.6 0.5 90 0.0129
Michiquillay 7º 02´ 6.4 731 1,191 0.61 3,825 8.7 1263 780 0.114 0.046 90 0.0101
Namora G 7º 12´ 249 335 970 0.35 3,369 11.2 1350 66 2.7 0.5 90 0.0187
Cajamarca 7º 15´ 789 249 797 0.31 3,140 12.4 1390 50 5.6 0.7 120 0.0173
Crisnejas 7º 18´ 4,150 323 883 0.37 3,300 11.5 1360 57 41.1 6.6 120 0.0152

La información muestra además los datos básicos de las

cuencas de los ríos Vilcanota, Paucartambo, Pitumarca, y
Pomacanchi.
 GENERACION DE CAUDALES EN LA SIERRA PERUANA
HANOCCA - TAYPITUNGA
GASTO DE LA RETENCIÓN
Para la estimación del coeficiente de agotamiento se utiliza
resultados de aforos ejecutados durante la estación seca y/o
el juego de ecuaciones empíricas siguientes:
Agotamiento muy rápido por temperatura elevada mayor de 10° C y
retención reducida (50 mm/año) hasta retención mediana (80 mm/año) ;
a = -0.00252 * Ln AR + 0.034
Agotamiento rápido por la retención entre 50 y 80 mm/año y vegetación
poco desarrollada (puna) ;
a = -0.00252 * Ln AR + 0.030
Agotamiento mediano por retención mediana (alrededor de 80 mm/año) y
vegetación mezclada (pastos, bosques y terrenos cultivados);
a = -0.00252 * Ln AR + 0.026
Agotamiento reducido por alta retención (arriba de 100 mm/año) y
vegetación mezclada.
a = -0.00252 * Ln AR + 0.023
 GENERACION DE CAUDALES EN LA SIERRA PERUANA
HANOCCA - TAYPITUNGA
En el caso que interese en particular calcular
descargas mínimas seguras durante la estación seca, se
puede elegir un valor “a” un poco más alto. La ecuación ( 11 )
indica que los caudales durante la estación seca están
relacionados a los del mes anterior por la siguiente relación
cuando se cuenta el mes por 30 días.
bo = e –a x 30
Donde:
bo = relación entre la descarga del mes actual y del mes anterior
a = Coeficiente de agotamiento
Durante la estación seca de m meses de duración, el caudal
disminuye en la relación :
Qj / Qo = b0j
Donde: Qj = Descarga del mes j
Qo = Descarga inicial
Boj = Relación entre las descargas Qj y Qo
 GENERACION DE CAUDALES EN LA SIERRA PERUANA
HANOCCA - TAYPITUNGA
La suma de los factores mensuales b j durante la o
estación seca de m meses :
bm = Σ
i=1
boi
Corresponde al gasto total durante el periodo de
meses secos o el agotamiento total de la reserva hídrica de
la cuenca. En consecuencia, la contribución mensual de la
retención durante la estación seca se puede indicar por la
relación siguiente:
Ci / R = boj / Σi=1 boj
Donde:

bo j = relación entre el caudal del mes actual y anterior

bo i = relación entre el caudal del mes i y el caudal inicial
Gi = gasto mensual de la retención (mm / mes)
R = retención de la cuenca (mm / mes)
 GENERACION DE CAUDALES EN LA SIERRA PERUANA
HANOCCA
ABASTECIMIENTO DE - TAYPITUNGA
LA RETENCION
La lámina de agua Ai que entra en la reserva de la
cuenca se muestra en forma de un déficit mensual de la
precipitación efectiva mensual PEi. Se calcula mediante la
ecuación :
Ai = ai * R / 100
Donde:
Ai = Abastecimiento mensual déficit de la PE (mm/mes)
ai = coeficiente de abastecimiento (%)
R = retención de la cuenca (mm/año)

CAUDAL MENSUAL PROMEDIO

La lámina de agua que corresponde al caudal
mensual para el año promedio se calcula según la ecuación
básica del balance hídrico a partir de la información básicas
de la cuenca:
 GENERACION DE CAUDALES EN LA SIERRA PERUANA
CMi = PEi + Gi - HANOCCA
Ai - TAYPITUNGA
CMi = caudal del mes i (mm/mes)
PEi = precipitación efectiva del mes i (mm/mes)
Gi = gasto de la retención en el mes i (mm/mes)
Ai = abastecimiento en el mes i (mm/mes)
GENERACION DE CAUDALES MENSUALES PARA PERIODOS EXTENDIDOS
METODO
El modelo hidrológico permite determinar los caudales
medios mensuales se requiere otros parámetros estadísticos,
sobre todo la desviación tipo que necesita para el cálculo de
caudales y consiste en una combinación de un proceso
Markoviano de primer orden, con una variable de impulso,
que para este caso es la precipitación efectiva cuyas
ecuaciones son las siguientes.
Qt = f ( Qt -1), y Qt = g ( PEt)
 GENERACION DE CAUDALES EN LA SIERRA PERUANA
HANOCCA - TAYPITUNGA
Para aumentar el rango de los valores generados y
obtener una aproximación óptima 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 combinando los componentes anteriormente
citados, se escribe:
Qt = B1 + B2 * Q t-1 + B3 * PEt + z * S* (1 – r 2 ) ½
Donde:
Qt = caudal del mes t
Q t-1 = caudal del mes anterior
PEt = precipitación efectiva del mes t
B1 = factor constante = caudal básico
S = Desviación de los residuos
Z = Números aleatorios (0,1)
r, r2 = Coeficiente de correlación y coeficiente de determinación
 GENERACION DE CAUDALES EN LA SIERRA PERUANA
HANOCCA - TAYPITUNGA
Los parámetros B1 , B2, B3, r y S se calculan sobre la base
de los resultados del modelo para el año promedio, mediante
un regresión triple con Qt como valor independiente y Q t -1 y
PEt , como valores independientes. En lo que se refiere a los
valores aleatorios, se puede emplear los valores publicados
en tablas o mediante la computadora.
Existen varias posibilidades para elegir el valor inicial
del proceso de generación, entre ellos se tiene:
- Empezar el cálculo en el mes para el cual se dispone de un aforo.
- Empezar con el caudal promedio de cualquier mes, y/o
- Empezar con el caudal cero, calcular un año y tomar el último valor
como Qo sin considerar estos valores en el cálculo de los parámetros
estadísticos del periodo generado.

La calidad de la coincidencia de los caudales generados con los

observados, se lleva a cabo mediante comparación de los promedios y
desviaciones mediante las pruebas del test de Student y de Fisher.
 GENERACION DE CAUDALES EN LA SIERRA PERUANA