IRD

(Instituto de investigaci6n para el desarrollo)

ANALISIS DE LOS COMPONENTES

NATURALES y ARTIFICIALES DEL VALLE
DEL ZONGO EN BOLIVIA

Vertiente Este dei Huayna Potosi vista dei valle deI rio Taypi KJ1UChu. Foto J.P. Chazarin

Por Igor REINHARDT

Memoria de posgrado de la Escuela Polytecnica Federal de

Lausanne (SUIZA)

Diciembre 1997
1
1 Agradecimientos

1 Quisiera agradecer a las personas siguientes por su contribuci6n en este trabajo:

Pierre Chevallier, mi dircctor de diploma, quien ha organizado mi curso de practica en el

1 ORSTOM y ha revisado (parcialmente) el francés deI manuscrito. Sin sus sugerencias y sus
consejos este informe no tendria su forma actual.

1 Victor Ticona de la COBEE, deI cual he obtenido la mayoria de los datos que son la base de
este trabajo, y quien ha dedicado mucho tiempo respondiendo mis preguntas. Sin su
contribuci6n este trabajo no habria podido realizarse.
1 Yves Arnaud, con quien he pasado horas en el analisis de los MNT. Su experiencia en la
materia y su ayuda generosa permitieron el esbozo de la modelizaci6n deI deshielo de los
1 glaciares.

Pablo Iturri de la COBEE, quien me asisti6 en la estimaci6n de los rendimientos de las plantas
1 hidroeléctricas.

Bernard Pouyaud, quien me recibi6 generosamente en La Paz y quien me dio buenos consejos
1 para este trabajo.

1 Bernard Francou, quien amablemente reley6 y corrigi6 una parte deI manuscrito.

A todas las personas que me ayudaron en los problemas cotidianos, en particular los

1 empleados deI ORSTOM y de la COBEE en Bolivia, quienes se mostraron muy amables y

hospitalarios.

1 Y final mente, a la Escuela Politécnica Feferal de Lausana, que financi6 este trabajo y mis
estudios con una beca.

1
1
1
1
1
1
1
1
1 Indice

1. Introducci6n
1 1.1. Objetivos deI estudio
1. 1.1. El marco deI proyecto
1.1.2. Objetivos generales y objetivos deI estudio
1 1.1.3. Metodologia deI estudio
1.2. Descripci6n de la zona de estudio

1 1.2.1. La cuenca de drenaje

1.2.2. Las plantas hidroeléctricas
2. Amilisis preliminar de las influencias deI sistema hidniulico artificial sobre los

1 caudales en Sainani
2.1. Los factores de influencia
2.1.1. El sistema hidniulico para la producci6n de electricidad

1 2.1.2. Las influencias posibles sobre los caudales diarios

2.2. La influencia de los cambios diarios deI stock de agua en Jas represas
2.2.1. La gesti6n de la producci6n planta hidroeléctrica -
1 2.2.2. Los datos disponibles
2.2.3. Estimaciones deI cambio de stock total diario
2.2.4. Estimaciones de los caudales turbinados en Sainani y de las
1 precipitaciones espaciales hasta Sainani
2.2.5. Presentaci6n de los resultados
2.3. Discusi6n y conclusiones respecto a la modelizaci6n hidrol6gica
1 2.3.1. Desalmacenamiento en estaci6n seca
2.3.2. Almacenamiento en estaci6n de I1uvias
3. Un método para reconstituir los caudales naturales provenientes deI deshielo de
1 los glaciares y deI escurrimiento de las vertientes
3.1. Principio deI método de reconstituci6n de los caudales naturales
3.1.1. La estimaci6n de Jos caudales na/urates sin aporte de los glaciares
1 3.1.2. La separaci6n de los periodos de modelizaci6n
3.1.3. El problema de la diferencia de horarios de los datos

1 3.1.4. El aporte hidrico deI valle Coscapa

3.2. Los métodos para la reconstituci6n de los caudales naturales
3.2.1. El método para la reconstituci6n de los caudales naturales en

1 estaci6n seca
3.2.2. El método para la reconstituci6n de los caudales naturaJes en
estaci6n de lIuvias

1 4. Mejoramiento y aplicaci6n deI método para la reconstituci6n de caudales

naturales en estaci6n seca
4.1. Reconstituci6n de caudales turbinados
1 4.1.1. Estimaci6n de los rendimientos de las plantas hidroeléctricas
4.1.2. Reconstituci6n de caudales turbinados
4.1.3. Discusi6n de los hidrogramas diarios medios de las plantas
1 hidroeléctricas
4.2. El papel de las pequefias cuencas de regulaci6n intermedia
4.2.1. Datos e informaci6n disponibles sobre las cuencas intermedias y sus
1 maneJos
4.2.2. Un método para la estimaci6n de la curva de almacenamiento dei

1 conjunto de las cuencas intermedias

 1
4.2.3. Simulaci6n de la curva de almacenamiento
1 4.2.4. La influencia deI cambio de stock en las cuencas intermedias sobre
los caudales naturales
4.3. Reconstituci6n de los caudales naturales en Sainani en estaci6n seca 1997
1 4.3.1. Reconstituci6n de los caudales naturales en estaci6n seca 1997
4.3.2. Discusi6n de los resultados
4.3.3. Estimaci6n cualitativa deI deshielo de los glaciares
1 5.1. Esbozo de una modelizaci6n deI agua de deshielo de todos los glaciares
5.1.1. Introducci6n
5.1.2. Las hip6tesis de un enfoque topografico
1 5.1.3. La importancia de la exposici6n
5.1.4. El concepto y la base fisica de un modelo a desarrollar

1 5.2. Comparaci6n de los caudales de deshielo y de las topografias de los

glaciares Charquini y Zongo
5.2.1. Los glaciares de referencia

1 5.3. Los datos obtenidos con ayuda deI :rvINT

5.3.1. Los MNT de los glaciares
5.3.2. Los parâmetros topognificos obtenidos con ayuda deI MNT para
1 cada malla
5.4. La estimaci6n de las diferencias de la radiaci6n solar directa durante el dia
y durante el ano en los glaciares Zongo y Charquini
1 5.4.1. Introducci6n
5.4.2. La posici6n deI sol en el transcurso deI ano y durante el dia
5.4.3. Calculos de los angulos de incidencias deI sol con las malIas de los
.1 glaciares
5.4.4. La insolaci6n media de las superficies totales de los glaciares
5.4.5. La insolaci6n media de las superficies bajo la isoterma cero
1 5.5. Simulaci6n de los caudales de los glaciares Charquini con dos modelos
simples y comparaci6n con la insolaci6n de los glaciares
5.5.1. Simulaci6n de los caudales de deshielo deI glaciar Charquini con dos
1 modelos simples
5.5.2. Interpretaci6n de los resultados de la modelizaci6n comparando las
diferencias de la insolaci6n
1 6. Conclusiones
7. Anexos

1 8. Referencias

1
1
1
1
1
1
1. Introducci6n
1 1.1. Objetivos dei estudio

1 1.1.1. El marco dei proyccto

Este estudio se inscribe en el marco deI programa « Nieve y glaciares Tropicales» (NGT), que
intenta evaluar el impacta de la variabilidad climatica sobre los sistemas hidrol6gicos en los
1 Andes tropicales. Estos sistemas hidrol6gicos son muy sensibles a la variabilidad climatica
(cf. FRANCOU, B. y RIBSTEIN, P. (1995)). Su funcionamiento esta vinculado al régimen de

1 los glaciares en las partes altas de los Andes. Si bien estos glaciares son pequefios en relaci6n
a las superficies de las cuencas, el papel que juegan es importante. El clima en estas regiones
tropicales se caracteriza pOl' una estaci6n de lluvias relativamente calida y una estaci6n seca

1 mas fria. En estaci6n de lluvias la contribuci6n de los glaciares al escurrimiento de los rios
puede sel' importante, pero es sobre todo en estaci6n seca 0 en las interestaciones que el
deshielo de los glaciares se vuelve el factor principal deI régimen hidrol6gico. Hasta ahora los

1 trabajos realizados en el marco del programa NGT han priorizado sobre todo el estudio dei
funcionamiento de los glaciares. Desde el inicio deI programa, en 1990, se lleg6 a comprender
bastante bien el funcionamiento de estos glaciares. Un modelo con base fisica para la
1 modelizaci6n de los caudales de deshielo se elabor6 y da resultados convincentes
(RIGAUDIERE, P. et al., 1995). Gracias a estos conocimientos deI funcionamiento de los
glaciares, ahora se puede extender la investigaci6n a los hidrosistemas directamente ligados a
1 sus regimenes. Frecuentemente, de estos hidrosistemas dependen las actividades del hombre,
las que a su vez influyen sobre eHos. Las consecuencias de la variabilidad climatica son de un
interés primordial para la continuaci6n de las actividades humanas.
1 Los glaciares estudiados se encuentran en la Cordillera Real de los Andes bolivianos. El
glacial' Zongo, cuyo funcionamiento ha sido mejor estudiado, cubre un costado deI macizo dei
1 Huayna Potosi que corona el valle de Zongo, y su escurrimiento alimenta el rio principal deI
valle, el rio Zongo. En los primeros cuarenta kil6metros de este rio estan instaladas varias
centrales plantas hidroeléctricas que explotan una caida de mas de 3000 m a corta distancia.
1 Es 10 que influy6 en la elecci6n de la cuenca de drenaje de Zongo para este cstudio; las
actividades deI hombre y los sistemas naturales se encontraban estrechamente vinculados.

1 1.1.2. Objetivos generales y objetivos dei estudio

A fin de poder estudiar los impactos de la variabilidad climâtica sobre el régi men hidrol6gico

1 deI valle de Zongo y las consecuencias sobre la explotaci6n deI recurso hidrico, hay que sel'
capaz de hacer predicciones sobre el comportamiento de los glaciares y sobre la producci6n de
las vertientes no glaciares bajo diferentes escenarios climaticos. Tanto el desarrollo de un

1 modelo hidrol6gico coma su aplicaci6n durante un estudio del funcionamiento deI conjunto
de los sistemas naturales y deI sistema hidraulico dedicado a la producci6n planta
hidroeléctrica, hacen necesaria una aclaraci6n sobre la vinculaci6n entre estos sistemas. Esta
1 ultima se manifiesta principalmente a través deI escurrimiento deI rio Zongo, donde
repercuten las influencias deI régimen de los glaciares, de la producci6n de las vertientes no
glaciares y deI sistema hidraulico dedicado a la producci6n de electricidad.
1 La meta principal de este trabajo es encontrar un método de deconvoluci6n, segun su origen,
de los escurrimientos deI Zongo a un paso de tiempo diario. Primero se intenta reconstituir los
1 caudales naturales deI rio Zongo, que no son inlluenciados pOl' los sistemas hidraulicos de

1
1
producci6n de electricidad y que no contienen los aportes deI deshielo de los glaciares. Estos
1 caudales naturales reconstituidos serviran luego para el perfeccionamiento de un modelo
hidrol6gico con el fin de modelizar la producci6n de las vertientes no glaciares.

1 La gesti6n sigue dos vias. En primer lugar, se trata de desarrollar un método de reconstituci6n
de los caudales naturales deI Zongo, incluyendo los aportes deI deshielo de los glaciares y, en
segundo lugar, mostrar c6mo se puede modelizar el deshielo deI conjunto de las superficies
1 glaciares en el valle a partir de los conocimientos adquiridos sobre su funcionamiento.

1 1.1.3. Metodologia deI estudio

El método general adoptado en este estudio es ir de 10 general a 10 particular. Es decir que se
analiza un componente deI sistema estudiando primero su influencia sobre el conjunto y,

1 segun el resultado obtenido, interesarse mas en el detalle. Este principio ha sido aplicado
sobre todo en el primer objetivo: la reconstituci6n de los caudales naturales incluyendo los
aportes de los glaciares. A partir de una estimaci6n de la importancia de la influencia deI

1 sistema hidraulico sobre los caudales naturales, se realiz6 un analisis preliminar de los
componentes dei sistema hidraulico cuya importancia parecia preponderante (Capitulo 2).
Interpretando el resuitado de este analisis, se estableci6 luego un método de reconstituci6n de
1 los caudales naturales deI conjunto de los glaciares y sus vertientes no glaciares (Capitulo 3).
Los conocimientos adquiridos durante el analisis preliminar y el desarrollo dei método
permitieron poner de relieve la influencia de otros componentes deI sistema hidraulico que
1 hubo que analizar antes de reconstituir los caudales naturales de un periodo de observaci6n en
1997 (Capitulo 4).

1 La estructura de este informe sigue la cronologia de nuestra investigaci6n, la cual fue

impuesta por la disponibilidad de datos. En cfecto, el sistema hidraulico ha sufrido cambios
importantes en los meses que precedieron a este trabajo. Como era esencial incluir estos
1 cambios en nuestro razonamiento, hubo que esperar el final de la estaci6n seca de 1997 para
disponer de suficientes datos para poder aplicar y evaluar nuestro método.

1 En cuanto al objetivo de modelizar el deshielo deI conjunto de los glaciares, se pudo

comenzar con el analisis en una escala mas detallada. La topografia de los glaciares es el

1 término que parece ser el factor principal a considerar en una modelizaci6n; fue analizado en
detaIIe (Capitulo 5).

1 1.2. Descripcion de la zona de estudio

1 1.2.1. La cuenca de drenaje

La cuenca de drenaje del rio Zongo esta esquematizada en la Figura 1. Su superficie es de
aproximadamente 400 km 2 a la salida de la tiltima planta hidroeléctrica en Cahua. A la salida
1 de la segunda planta hidroeléctrica subiendo el rio, en Chururaqui, la cuenca de drenaje
contiene todas las superficies glaciares cuya agua de deshielo alimenta los escurrimientos dei

1 rio Zongo. El esquema representa la cuenca de drenaje dei rio Zongo en Chururaqui. La planta
hidroeléctrica de Sainani, aguas arriba de Chururaqui, debia ser equipada con un limnigrafo a
finales deI ano 1997. La mayor parte deI agua que alimenta la planta hidroeléctrica de Sainani

1 es captada por una toma en el rio principal, justo aguas abajo de las plantas hidroeléctricas
Santa Rosa 1 y II. Suponiendo que no hay desbordamiento de la toma de Santa Rosa ni de sus
tomas en el canal de derivaci6n en la planta hidroeléctrica de Sainani, la cuenca de drenaje
1
1
que alimenta la planta hidroeJéctrica corresponde al de una salida al frente de la planta

1 hidroeléctrica Sainani, a excepci6n de una pequefia cuenca de aproximadamente 2 km cuyo

escurrimiento no pasa par la planta hidroeléctrica (Figura 1). Raz6n por la cual, en ausencia
de series de caudales medidos, se eligi6 la cuenca de drenaje de Sainani como unidad de base

1 de este estudio,

Anticipando un resultado de este trabajo, se estima el total de superficies glaciares en la

1 cuenca de drenaje de Zongo en 18.3 km 2 , 10 que hace un 5,8% de la superficie de la cuenca en
Chururaqui. Dejando de lado los aportes de la planta hidroeléctrica Santa Rosa II al
escurrimiento deI rio Zongo en Sainani que comprenden aguas de deshielo de los glaciares dei
1 valle de Coscapa, la cuenca de drenaje en Sainani comprende Il.5 km 2 de superficies
glaciares, es decir el 6.6.% de la superficie total de 172.6 km 2• El punta mas elevado de la
cuenca de drenaje cs la cima dei macizo dei Huayna Potosi, a 6050 m, mientras que la salida
1 de Sainami se cncuentra a una altura de 2210 m. Esto da globalmente pendientes
extremadamente inclinadas.

1 Las subcuencas prcsentadas en la Figura 1 corresponden al estado de la cuenca de drenaje no

influenciado por el sistema hidraulico dedicado a la producci6n de electricidad. Las tomas de

1 agua y los canales que alimentan las plantas hidroeléctricas alteran considerablemente esta
representaci6n, particularmente en estaci6n seca cuando los escurrimientos son minimos. En
efecto, una gran parte de los escurrimientos de las subcuencas es captada por el sistema

1 hidraulico y desviada a otras subcuencas, correspondiendo en general a subcuencas situadas

mas arriba de las subcuencas de origen. Ademas, una pequefia subcuenca de una superficie
total de 10 km 2 fue captada en el exterior de la cuenca de drenaje natural, en el valle de

1 Tiquimani. El sistema hidraulico de producci6n de electricidad se describe en detalle en el

segundo capitulo de este trabajo. El parrafo siguiente describe brevemente las plantas
hidroeléctricas.

1 1.2.2. Las plantas hidroeléctricas

En el valle dei Zongo, la Compafiia Boliviana de Energia Eléctrica (COBEE) tiene en
1 funcionamiento actualmente nueve plantas hidroeléctricas, con una capacidad total instalada
de 123.3 MW. La mas reciente de estas plantas hidroeléctricas es la de Tiquimani, que fue
puesta en servi cio en febrero de 1997. Toda la energia producida es utilizada para alimentar la
1 ciudad de La Paz. DUI'ante las horas de punta, a mediodia y por la tarde, la producci6n no es
suficiente para satisfacer la demanda de la ciudad, y la COBEE debe importar electricidad de
otra parte. Hasta 1998, la capacidad total va a ser aumentada a 42% de la capacidad actual por
1 las ampliaciones de las plantas hidroeléctricas existentes y por una nueva planta hidroeléctrica
(Huaji), llevandola a 173.6 MW (Cuadro 1). Este aumento se hace sobre todo para utilizar
mejor el potencial deI recurso hidrico en estaci6n de lluvias. Efectivamente, la capacidad
1 actual de las plantas hidroeléctricas esta adaptada al potencial dei recurso en estaci6n seca. En
estaci6n de lluvias, este potencial es mucha mas grande que la capacidad instalada de las

1 plantas hidroeléctricas que producen a la capacidad maxima y la mayoria de las tomas

desbordan. En estaci6n de lluvias, el aumento de la capacidad permitira no solamente
satisfacer la demanda de las horas punta, pero también exportaI' electricidad hacia otros

1 lugares de consumo fuera de la ciudad de La Paz.

1
1
-------------------
Cerro Ilarnpu Cerro Telata Saltuni Zongo Jisthàna
_ _ _ _ _---'J~isthana
5900
5600
5300
5000
4700

4400

4100

3800

3500

3200

2900

2600

2300

2000

1700 .-----
Légende: 0 chambre de charge - usine hydroéléctrique canal ~ surfaces glaciaires
L::.. prise d'eau = retenue, lac cours d'eau naturel C bassin extérieur!artiliciel
Figura 1. Esquema dei sys/ema hydraulico de produccion de elec/ricidad dei valle dei zongo. En el esquellla es/an indicadas las sub-cuencas na/urales entre cada toma
de agua sobre ellecho principal, con sus superficies individuales/acumuladas. Las lagunas/represas some/idas a operaciones hydraulicas es/an subrayados.
1
1 Planta hidroeléctrica: Capacidad instalada (MW)
hidroeléctr. Acumulada % total
Capacidad futura
hidroeléctr. Acumulada % total
Zongo 4.4. 4,4 3,6% 10.0 10.0 5,7%
1 Tiquimani
Botijlaca
9.6
3.5
14
17,5
7,8%
2,8%
9.6
6.6
19,6
26,2
5,5%
3,8%
Cuticucho 8.4 25,9 6,8% 20.0 46,2 Il,5%

1 Santa Rosa [
Santa Rosa II
3.2
6.8
29,1
35,9
2,6%
5,5%
6.2
6.8 -
52,4
59,2
3,6%
3,9%
Sainani 10.4 46,3 8,4% 10.4 69,9 6,0%

1 Chururaqui
Harca
Cahua
24.0
26.0
27.0
70,3
96,3
123,3
19,5%
21,1%
21,9%
24.0
26.0
27.0
93,6
119,6
146,6
13,7%
14,9%
15,5%

1 Huaji

Cuadro f
28.0 174,6 16,0%

Capacidades instaladas y capacidadesJuturas de las plantas hidroelécricas

1 en el valle de Zongo. En negrillas: capacidades aumentadas

Como las plantas hidroeléctricas estan organizadas en cascada, sus capacidades aumentan

1 al descender el valle. Hay dos plantas hidroeléctricas de pequefia capacidad aguas arriba
(Zongo y Botijlaca), seguidas de tres plantas hidroeléctricas de capacidad media
(Cuticucho, Santa Rosa 1 y II, Sainani) y de tres plantas hidroeléctricas de gran capacidad .
1 aguas abajo (Churaqui, Harca y Cahua). La nueva planta hidroeléctrica (Tiquimani), que es
alimentada por caudales captados deI exterior de la cuenca y que restituye sus
escurrimientos entre Zongo y Botijlaca, es de capacidad media. El aumento de la capacidad
1 se hace a través de la ampliaci6n de tres pequefias plantas hidroeléctricas y agregando una
nueva planta hidroeléctrica de gran capacidad mas abajo de Cahua en Huaji.

1 El aumento de la capacidad instalada no cambia intrinsecamente la influencia deI sistema

hidrauIico, dedicado a la producci6n de electricidad, sobre los caudales naturales. En
estaci6n seca, cuando esta influencia es la mas importante, como se mostrara de manera
1 cuantitativa en este trabajo, la capacidad adicional no es utilizada, y en estaci6n de lluvias,
las tomas continuaran rebalsando porque la capacidad instalada seguira siendo muy inferior

1 al potencial dei recurso disponible durante esta estaci6n. La raz6n por la cual no se utiliza
completamente este potencial es de orden econ6mico. Como los precios de la electricidad
son relativamente bajos en estaci6n de Iluvias, los costos marginales de la producci6n de

1 electricidad en determinada planta hidroeléctrica vienen a ser rapidamente superiores a las

ganancias marginales, 10 que hace no rentables los aumentos de su capacidad mas alla de
cierto limite. En consecuencia, los resuitados obtenidos en este estudio seran igualmente

1 validos en el futuro.

1
1
1
1
 1
1 [Link] preliminar de las influencias dei sistema hidraulico
artificial sobre los caudales en Sainani

1 2.1. Los factores de influencia

1 2.1.1. El sistema hidrâulico para la producci6n de electricidad

Las plantas hidroeléctricas en el valle de Zongo estan organizadas en cascada. La represa
Zongo se encuentra a la cabeza deI sistema. Es alimentada princîpalmente por el agua de

1 deshielo de los glaciares que cubren las cimas que coronan el valle, el Huayna Potosi y el
Charquini. El agua es captada en pequenos lagos al pie de los glaciares y llevada al lago
Zongo por cuatro canales. El canal Prado capta los cscurrimientos deI glaciar Charquini. El
1 canal Tubo Vertedero lleva las aguas de deshielo deI glacial' Zongo y los canales Alpaca y
Vicufia captan una parte de las aguas de deshielo des los otros glaciares deI Huayna Potosi.
La represa de Zongo alimenta directamente los conductos forzados de la planta
1 hidroeléctrica Zongo, la primera de la cascada. Luego el agua es captada en el rio principal
aguas abajo de cada planta hidroeléctrica y dirigida hacia la planta hidroeléctrica siguiente
pOl' medio de canales 0 de tuneles (cf. 1.2.1, Figura 1).
1 Las represas Hatilata, Taipicota (lagunas Tiquimani) y Guaraguarani constituyen una
segunda cabeza de sistema. Los lagunas Tiquimani se encuentran al exterior de la cuenca
1 de drenaje deI Zongo en una pequena cuenca deI valle Tiquimani (aproximadamente 6
km 2), que contiene un pequefio glacial' (Cerro Telata). El agua es llevada por el tunel
Taipicota y el canal Taipicota hasta el lago Guaraguarani, que se encuentra en el valle
1 Tiquimani al exterior de la cuenca natural de Zongo. Una toma de agua dellago redirige los
escurrimientos hacia el tunel Guaraguarani y el canal Tiquimani hasta la cuenca de carga de

1 la nueva planta hidroeléctrica Tiquimani. El pequeno canal Tiquimani capta dos afluentes
deI Zongo que drenan una superficie total de aproximadamente 2 km 2 Y que no fluyen en
estaci6n seca.

1 Las represas Livinosa y Sankayuni forman parte de dos sistcmas de toma que funcionan de
la misma manera que el de las represas Tiquimani/Guaraguarani, salvo que en estos casos
1 el agua proviene deI interior de la cuenca deI Zongo. Contrariamente al sistema deI lago
Zongo, las dos subcuencas drenadas (Livinosa 27 km2 y Sankayuni 75 km 2) son
esencialmente no glaciares, aunque contienen glaciares importantes, pero cuya influencia es
1 minima. El sistema de Livinosa comienza a alimentar la cascada de las plantas
hidroeléctricas en Botijlaca, el deI Sankayuni en Santa Rosa II. La planta hidroeléctrica
Santa Rosa II y la de Tiquimani no forman directamente parte de la cascada de plantas
1 hidroeléctricas, sino que son aflucntes de ésta.

2.1.2. Las influencias posibles sobre los caudales diarios -

1 El conjunto hidroeléctrico constituido por las represas, las tomas de agua, los canales y las
plantas hidroeléctricas influyen en el sistema natural de dos maneras. POl' una parte, el

1 escurrimiento es parcialmente encaminado hacîa aguas abajo a través deI sistema

hidraulico, 10 que influye en la forma de Jas crecidas, Por otra parte, cierta proporcion dei
escurrimiento es almacenado en las represas en estacion de Iluvias -periodo de abundancia-

1 .y desalmacenado durante la estaci6n seca -periodo de falta de agua-. En el parrafo

siguiente, la importancia cualitativa de estas dos influencias sera brevemente an~lizada.

1
1
Se puede admitir que la transferencia de escurrimientos a través deI sistema hidniulico no
1 tiene influencia sobre los caudales diarios en Sainani. Esto esta demostrado por la
comparacion de los tiempos de propagacion en el sistema hidraulico (canales / tuneles /
lechos naturales) y en el sistema natural (Cuadro 2). Las vclocidades de propagacion en las
1 transferencias naturales fueron medidas solo en ciertos tramos. Para establecer el cuadro, se
admite la hip6tesis simplicadora que éstas son constantes a 10 largo de los tramas referidos.
Consideramos por ejemplo la toma de agua Alpaca, aguas arriba de la cuenca de drenaje,
1 que constituye el camino mas largo que puede tomar el escurrimiento. Por via natural, el
agua llega 2 h 35 mas temprano a la represa Dique Intermedio que el agua que pasa por la
planta hidroeléctrica Zongo. En los casos de los otros canales, las diferencias de los
1 tiempos de propagaci6n son minimas, 20 minutos para la transferencia entre la toma
Livinosa y la planta hidroeléctrica Botijlaca, y 45 minutos para el encausamiento entre la
toma Coscapa y la planta hidroeléctrica Chururaqui.
1 N° De (toma): A (union) Distancias Velocidad Tiempo
lecho natural canales lecho natural lecho natural canales
1 1)
2)
Alpaca
Livinosa
Dique inlcrmctlio
Ilidroci. Botijlaca
3km
4 km
8.7 km
6.8 km
1.14 km/h
2 km/h
5 h 15
1 h 40
2 h 40
2h
3) Coscapa 1lidrocl. Chururaqui 9.5 km 15 km 3.5 km/h 1 h 55 2 h 40

1 Cuadro 2 Tiempo de propagacion entre algunas tomas de agua y la confluencia dei

rio correspondiente con el rio principal (rio Zongo). Estimacion de las velocidades de
1 transferencia natural, para 1): 3 h 30 para 4 km entre la planta hidroeléctrica Zongo y
Dique Intermedio; para 2): 1 h para 2 km entre la represa Livinosa y la toma Livinosa; 3):
2 h para 7 km entre la represa Sankayuni y la toma Coscapa
1
1 .., ,.,
Tl
r ••
450000 0
"'". "'
r
1 l
~
1

400000

1 350000
5

...::s 10 ...::s
1 0

"E
M
300000

250000 15
0

"EE

1 200000
Débit 20
150000

1 100000
'0 '0 \0 \0 '0 '0 \0 '0
~ \0 '0 '0 \0 '0 '0 '0 '0 '0
25

---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- ---- -- --g:- -- ------ ------
0\ 0- 0- 0- 0\ 0- 0- 0\ 0- 0- 0- 0- 0- 0- 0- 0- 0-
or) 'n V'l '0 '0 \0 r-- r-- r-- 00 00 00 0- 0; o 0 0

1 ~ ~ '<t
N
M M M
N
M M

date
r"l
N
N N N
N
- N
N

1 Figura 2 Precipitaciones espaciales hasta Sainani y caudales turbinados en la planta

hidroeléctrica Sainani
1
1
1
Para tener una idea de la importancia de estas diferencias hay que compararlas con la
1 duraci6n de las crecidas. Todavia no se dispone de mediciones directas dei caudal dei rio
Zongo en Sainani, pero es posible reconstruir los caudales turbinados en la planta
hidroeléctrica Sainani a través de la energia diaria producida (cf. 2.2.4). En estaci6n seca
1 casi todo el escurrimiento deI Zongo es turbinado. El caudal que no es captado en la planta
hidroeléctrica aguas arriba (Santa Rosa) es mu)' débil. Se puede observar en la Figura 2 que
la duraci6n de las crecidas en estaci6n seca de 1996 era siempre superior a 3-4 dias. La
1 diferencia de algunas horas de una parte dei escurrimiento por el sistema hidraulico
artificial no puede tener influencia significativa sobre las crecidas. En consecuencia, esta
influencia no sera estudiada de manera detallada en este trabajo.
1 A primera vista, la influencia dei almacenamiento de agua en las crecidas parece poco
importante; en efecto, la suma de los volumenes totales utiles de todas las represas es
1 pequena en relaci6n con las precipitaciones anuales. En relaci6n a la superficie de la cuenca
de drenaje en Sainani, solamente 44 mm de agua son almacenados y desalmacenados
durante un ano, 10 que s610 constituye un 4% de las precipitaciones anuales. Si se considera
1 que el desalmacenamiento de las represas tiene lugar en estaci6n seca, cuando casi no
llueve, la influencia a escala mensual 0 diaria puede ser mucho mas importante. Las

1 represas son vaciadas durante los meses mas secos, desde comienzos de mayo a fines de
agosto. Durante este periodo s610 llueve 191 mm (promedio ca1culado en 25 anos; ENDE
1992). Suponiendo que los caudales especificos no son muy diferentes de esta altura de

1 lluvia, la influencia dei desalmacenamiento seria superior al 23% del caudal en Sainani en
estaci6n seca. Es por ello que se tratara de cuantificar la influencia dei almacenamiento de
agua en las represas en relaci6n a los caudales diarios con un primer analisis de los datos.
1 Este amilisis servira primero para seleccionar las represas donde el cambio de stock influye
de manera significativa en los caudales medidos en Sainani. Esto conducira luego a la
elaboraci6n de un método para la reconstituci6n de los caudales naturales que provienen
1 des los escurrimientos de las vertientes y dei deshido de los glaciares.

1 2.2. La influencia de los cambios dia rios dei stock-de agua en las rcpresas

2.2.1. La gestion de la produccion planta hidrocléctrica

1 La gesti6n dei stock de agua en las represas depende de la poHtica economlca de la
empresa que administra la producci6n planta hidroeléctrica, la Compania Boliviana de

1 Energia Eléctrica (COBEE). Como la capacidad de producci6n en el valle no es suficiente

para satisfacer la demanda local de La Paz, la COBEE debe comprar electricidad a otros
operadores a un precio mas alto que el de su propia producci6n. En estaci6n de lluvias, hay

1 una abundancia de agua pero la capacidad de las plantas hidroeléctricas limita la

producci6n. Durante este periodo se escurre mucha agua sin ser utilizada para la
producci6n. De hecho, los aumentos de las plantas hidroeléctricas que se realizan

1 actualmente sirven para poder satisfacer la demanda de punta en estaci6n de lluvias (cf.
1.2.2). En estaci6n seca en cambio, hay una carencia de agua y la producci6n es limitada
por la disponibilidad de agua. El precio de la electricidad es mas alto en esta época. Como
1 se ha visto en el capitulo anterior, el volumen total util de las represas es relativamente
bajo, 10 que limita las posibilidades de regulaci6n de la producci6n de la planta
hidroeléctrica. El principio es desalmacenar las represas durante toda la estaci6n seca y
1 esperar el minimo de stock al comienzo de la estaci6n de lluvias. La producci6n en las
plantas hidroeléctricas principales dei sistema esta adaptada a la evoluci6n de la demanda

1
1
durante el dia. El escurrimiento aguas abajo de las represas entra inmediatamente en la
1 cascada de las plantas hidroeléctricas. La producci6n de las plantas hidroeléctricas aguas
abajo esta determinada, por una parte, por los caudales que lIegan de las plantas
hidroeléctricas principales deI sistema y, por otra parte, por los escurrimientos naturales de
1 las otras partes de la cuenca de drenaje que no son concentradas hacia las represas.

La gesti6n deI stock de agua en las represas es sobre todo importante en estaci6n seca. La
1 operaci6n de las valvulas es diferente en cada represa. No hay automatizaci6n; todas las
valvulas son reguladas manualmente. Las maniobras hidraulicas en la represa Zongo se
hacen durante todo el ano. Cada mes son previarnente establecidos niveles de alcance. El
1 objetivo es que la represa alcance su nivel maximo al final de la estaci6n de lIuvias, y que
sea vaciada hasta un minimo durante la estaci6n seca. Durante el periodo de vaciado, se

1 consideran las entradas medias (run-ins) de cada mes, es decir la acumulaci6n de energia
potencial en el sistema por medio de las precipitaciones y los aportes de los glaciares, de 10
cual sc sustrac una producci6n mensual que dcpcndc de la energia total almacenada en la

1 represa al final de la estaci6n de lIuvias y de los dias habiles deI mes en referencia. Esta
diferencia corresponde al descenso deI nivel de agua en términos de energia que hay que
alcanzar al final deI mes. Las vâlvulas de todas las otras represas son reguladas s610

1 durante la estaci6n seca en duraciones que varian segun el volumen tItil de diez dias (San
Pedro) a tres meses (Sankayuni). En estaci6n de lIuvias las valvulas de estas represas
permanecen cerradas (Cuadro 3).
1 2.2.2. Los datos disponibles
Durante el periodo de operaci6n de las represas hay una lectura diaria deI nivel de agua,
1 mientras que fuera de este periodo los niveles son leidos cada dos a siete dias, salvo
algunos casos en los que no se efectua, ninguna lectura (caso de la represa San Pedro). El
cuadro 3 muestra la frecuencia y las horas de lectura de los niveles de agua en las represas.
1 Nombre Volumen Relaci6n Operaci6n 96 Operaci6n/Jectura Lectura de nivel fuera
(10 3 m3) [Link]. dei al abierto cerrado periodo de operaci6n
1
1 1

Zongo 3256 todo el ano operaci6n continua 7h30 y 16h00

Mamankota 236 desconocida 26.7 2.9 op cont. 9h30 cada 2 dias, 9h30
Hatilata 2483 conocida 19.8 15.10 op cont. 9hOO cada 2 dias, 8hOO

1 Taipicota

Guaraguarani
3000
se~undo periodo de operacion
432
conocida
1996:
conocida
19.6
10.10
10.6
22.9
24.10
26.6
op cont. 9h30

op cont. 12h00
cada 2 dias, 8h30

cada 2 dias, Il hOO

1 A. Viscanchani desconocid

Livinosa
0
523
desconocida 5.6

dcsconocida 24.7
22.9

229
op cont. 5h45

8h30 20h30
cada 7 dias, ? h

cada 2 a 7 dias, ? h
1.7
1 Sankavuni 443 conocido 26.9 16h00 19h30 no lectura
San Pedro desconocid desconocida 18.6 29.6 13h30 17h30 no lectura
0

1 Cuadro 3 Volumenes de las represas mas importantes y caracteristicas de su

maniobra. Salvo indicacion contraria, la maniobra de la valvula y/a la lectura dei nivel se

1 hacen cada dia

Ademâs de las observaciones de los ni veles, se tienen mediciones diarias de las

1 precipitaciones leidas a las 7 de la mafiana y de la cantidad de energia producida en cada

planta hidroeléctrica anotada cada media hora.

1
1
2.2.3. Estimaciones dei cambio de stock total diario
1 Como los niveles de agua en las represas no son leidos cada dia y como las caracteristicas
dimensionales de las represas son generalmente desconocidas (Cuadro 3), la variaci6n
diaria de stock no es facil de estimar y es necesario hacer hip6tesis simplificadoras.
1 En las represas mas importantes, Zongo, Hatilalta, Taipicota, Guaraguarani y Sankayuni,
las relaciones entre el nivel de agua y el volumen son conocidas. En cambio, s610 se
1 conocen los volumenes totales utiles de las represas Livinosa y Mamankota, que son de una
importancia mediana (Cuadro 3). Para estimar las relaciones entre el nivel de agua y el
volumen de estas represas, se ha aproximado sus formas tridimensionales con casquetes
1 esféricos l . Conociendo el volumen total, se puede calcular el radio y la superficie de la
secci6n deI casquete. Esta superficie corresponde a la superficie de la represa cuando esta
Bena. Las superficies de las dos represas calculadas de esa manera son bastantes parecidas
1 a las superficies determinadas' por planimetro, 10 que muestra que la aproximaci6n es
aceptable (Anexo A). Respecto a las dos pequefias represas, Alto Viscanchani y San Pedro,
ni siquiera se conocen sus volumenes. Pero acercando sus formas a través deI método de
1 los casquetes esféricos, el calculo de los radios y de los volumenes puede hacerse a partir
de las superficies determinadas por planimetraje. De esta manera se obtienen los radios de
los casquetes, 10 que permite establecer las relaciones entre el nivel de agua y e! volumen.
1 Se pueden calcular los volumenes totales utiles de estas represas modelizadas, pero como
no se conocen los volumenes reales no se tiene base de validaci6n. No obstante, los
volumenes calculados (Alto Viscanchani: 19'500 m3, San Pedro: 9'000 1113) son probables.
1 Se ha ùtilizado esta aproximaci6n y las relaciones nivel / volumen conocidas para
transformar las lecturas de nivel de todos los lagos dei Cuadro 3 en volumenes

1 almacenados 2, para el periodo de enero 1996 a abril 1997.

Hay otros dos lagos y varias pequefias cuencas de regulaci6n cuyos niveles son observados
1 y los que no han sido tomados en cuenta en este estudio preliminar. Son de la dimension de
San Pedro, que fue elegido como represa Lipo. Las cuencas que se encuentran en el rio
principal sirven para la regulacion durante el dia y tienen cambios de nivel horaires. Como
1 las puntas de produccion planta hidroeléctrica son a mediodia y por la tarde, el agua es
soltada suficientemente temprano para que Begue antes de mediaoche a Sainani, hora en la
cuallos caudales son incluidos en el caudal diario (cf. 2.2.4).
1 Para obtener cronicas diarias completas de los volumenes almacenados se han estimado los
datos faltantes por interpolaciones lineales entre dos dias de los que se conocen los
1 volumenes3 . Se tuvo que interpolar en periodos comprendidos entre uno y siete dias (cf.
Cuadro 3). Cuando hal' pocas entradas en las represas, es decir en estacion seca, los datos
interpolados entre dos volumenes conocidos deben dar una aproximacion de los valores
1 reales. En efecto, en ausencia de lectura de nive!, la posicion de las valvulas no es

1 1 El casquete esférico es una de las partes de una esfera cortada en dos par una secci6n horizontal. Su
volumen y la superficie de la secci6n, que corresponde a la superficie de la represa, son funciones de la altura
de la parte 'cortada' y dei radio de la esfera 'entera'. La relaci6n nivel/volumen es simplemente V=f (radio,
nivel). El radio es un panimetro constante de la represa, que puede ser calculado a partir dei volumen total 0
1 de la superficie maxima, conociendo el nivel (=Ia altura) maxima (Anexo A).
2 A excepci6n dei caso de la aproximaci6n con el casquete esférico, las relaciones nivel/volumen no son
funciones analiticas. Por 10 tal1to hay que buscar en los cuadros el volumen correspondiente a cada nive!. A
1 fin de facilitar esta transforrnaci6n en cinco series de mas de 600 datos, se la ha automatizado con ayuda de
un programa en Visual Basic para Microsoft Excel, versi6n 7.0.
3 Estas interpolaciones también han sido realizadas con ayuda de un programa en Visual Basic.

1
1
modificada y las variaciones de nivel son bajas en relaci6n a los cambios de volumen. Sin
1 embargo, en estaci6n de lluvias, cuando las represas reciben aportes de agua importantes, la
interpolaci6n puede ser err6nea. Se tomani en cuenta [Link] durante la discusi6n de los
resultados (cf. 2.3).
1 El desfase de las horas de lectura es otro problema sobre los datos de nivel de agua en las
represas. Se puede verificar en el Cuadro 3 que las horas de lectura son distirbuidas en todo
1 el dia. Para obtener una primera aproximaci6n dei cambio de stock total se ha elegido la
primera lectura dei dia en los casos cn quc hay dos y sc han ignorado las diferencias.

1 2.2.4. Estimaciones de los caudales turbinados en Sainani y de las precipitaciones

espaciales hasta Sainani

1 Para conocer la influencia dei cambio dei stock de agua diario en las represas sobre los
caudales en Sainani, es necesario tener series diarias de estos caudales. Hasta el momento
no se tienen mediciones de caudal en Sainani. Se debe entonces reemplazar los caudales

1 reales por caudales turbinados, que pueden ser reconstituidos a partir de la energia
producida cada dia (cf. 1.2.1). No obstante, este enfoque tiene sus limites. Hay que conocer
los rendimientos de las maquinas para poder reconstruir los caudales turbinados 4 • El

1 rendimiento de la producci6n planta hidroeléctrica no ès lineal, depende de la energia

producida por unidad de tiempo. En la planta hidroeléctrica Sainani se midieron algunos
puntos de curva capacidad / rendimiento s. No obstante, por razones de simplicidad se ha

1 utilizado el factor de planta hidroeléctrica para el calculo, que se basa en el rendimiento

maximo, admitiendo un pequell0 error6 •

1 Las producciones diarias de electricidad de las plantas hidroeléctricas son calculadas de

medianoche a medianoche. La transformaci6n de estas producciones en caudales
turbinados da igualmente los caudales acumulados durantc el dia de medianoche a
1 medianoche. Son las cr6nicas de caudales que se utilizan para este estudio preliminar.
Como las producciones de las plantas hidroeléctricas son anotadas cada media hora, es
posible calcular los caudales turbinados diarios en otros intervalos. Esta posibilidad sera
1 utilizada en el método de reconstituci6n de los caudales naturales, como se vera en el
capitulo 3.1.3.

1 Los caudales turbinados son truncados en estaci6n de lluvias, cuando las turbinas alcanzan
su capacidad maxima (cf. 2.2.1). Este caudal maximo es de aproximadamente 2.8 m 3/s en

1 la planta hidroeléctrica Sainani, 10 que corresponde a un caudal especifico de 1,85 mm /

dia. Al final de la estaci6n de lluvias (1995/96 los caudales turbinados en la planta
hidroeléctrica Sainani comienzan a ser inferiores a este limite a partir dei 3 de mayo y

1 vuelven a pasarlo a partir dei 23 de octubre. Es posible que la capacidad maxima sea ya
alcanzada el 24 de septiembre, fecha de fin de la estaci6n seca, cuando los caudales
turbinados especificos comienzan a alcanzar valores alrededor de 1,7 mm: dia (Figura 3).

1
4 Conociendo este rendimiento, los caudales se calculan con la formula: Q (m 3/s) = E • 100/ (r • g • h • p),

1 donde E = energia producida por unidad de tiempo (kW), r = rendimienlos (%) h = altura de caida (m), p =
1000 kg!m 3 y g = 9.81 rn/s 2,
5 Como la produccion de electricidad nunca es constante, los aforos aguas abajo de las turbinas no son muy
exactos. Un método mas fiable utiliza mediciones de presion en el conducto asi como su diâmetro para
1 calcular un caudal instantaneo. Este método sera utilizado para calcular los caudales turbinados que se
utilizara para el perfeccionamiento dei modelo hidrologico.
l
6 El factor de planta hidroeléctrica es calculado con la formula: f(kWh/m ) = r· g • h/3600.

1
1
1 Los caudales turbinados no proporcionan ninguna informaci6n sobre los caudales totales en
estaci6n de lluvias. En estaci6n seca, en cambio, casi todo el caudal deI Zongo en Sainani

1 es turbinado. Se estima que mas deI 90% deI caudal deI rio Zongo es captado por la toma
aguas abajo de las plantas hidroeléctricas Santa kosa y que es llevado a la planta
hidroeléctrica Sainani. Hay dos pequeiias subcuencas aguas abajo de Santa Rosa donde el
agua no es turbinada en Sainani, que tienen una superficie de aproximadamente 2 km 2 (la
1 superficie total de la cuenca de drenaje hasta la salida de Sainani es de 173 km 2) (cf. 1.2.1,
Figura 1).

1 Por los tres canales, Coscapa y Tributario J y II, una parte deI agua de la subcuenca deI rio
Coscapa, un afluente deI Zongo, es captada y llevada a la planta hidroeléctrica Santa Rosa

1 II. Como este afluente se une al Zongo aguas abajo de Sainani, se trata de un aporte de agua
deI exterior de la cuenca de drenaje de Sainani. Los caudales que se vierten en las tomas de
Coscapa y de estos dos afluentes (Tributario 1 y II) son muy débiles en estaci6n seca. Hay
1 algunos arroyos que atraviesan los tres canales y donde un débil caudal es captado en
estaci6n de lluvias. En estaci6n seca no se pudo observar escurrimiento en estos arroyos, ni
durante un evento lluvioso bastante fuerte. Se volvera a hablar sobre la importancia de la
1 toma de agua en el valle Coscapa en el transcurso deI capitulo 3.1.4.

A fin de tener una base de comparaci6n para el almacenamiento en las represas en estaci6n
1 de lluvias, también se ha calculado la cronica de las precipitaciones espaciales en la cuenca
de Sainani sirviendose de las series de las seis estaciones pluviométricas entre la
Plataforma Zongo y la planta hidroeléctrica Sainani y utilizando el método dt: isohietas
1 (obtenidas de ENDE 1992). Para atenuar la variabilidad de las lluvias diarias y para
comparar con los caudales, se ha aproximado las precipitaciones por medias m6viles en 4

1 dias (Figura 3).

2.2.5. Presentacion de los resultados

1 Los caudales turbinados en Sainani y las precipitaciones diarias de 4 dias son mostrados en
la Figura 3 por todo el periodo de observaci6n de enero 1996 a abril 1997. Se nota
claramente el comienzo y el final de la estaci6n seca 1996, tanto por la disminuci6n de las

1 precipitaciones como de los caudales. El desalmacenamiento deI sistema, calculado por la

suma de los cambios de stock de cada represa,· es comparado a los caudales turbinados
durante la estaci6n seca en la Figura 4. En la Figura 5, la influencia dei desalmacenamiento
1 sobre los caudales esta expresada por el porcentaje dei caudat que proviene de las represas.
Las influencias de cada represa se encuentran en el Anexo. En estaci6n de lluvias una parte
dei escurrimiento es utilizada para llenar las represas. Como no se tienen mediciones deI
1 caudal en estaci6n de lluvias (cf. 2.2.4), el cambio de stock total sera comparado con las
precipitaciones medias (Figuras 6, 7 y 8).

1
1
1
1
1
1 10
9
,. '1· ~
::1
o
~
1 ...
8
7
,
+10 :~
E
....,
~
(;j
::1
E
1 6
0 !
'.:::" ::1
E o
E 5 15 ';;;
u
li (;j
:ë
~
4 , .~

1 "0

3
2
saison sèche 96
~
20 ~
'"c
,2
~
25 '5.
1 ! • 30
'ü
...
·U
0.

° '-0
0-
.....
'-0
~
'-0
~
'-0
~
'-0
~
'-0
~
'-0
~
'-0
~
'-0
~
'-0
~
'-0
~
'-0
~
'-0
~
l"-
~
l"-
~
l"-
~
l"-
~

1 ..... .....
r'"l
~
-.....
r'"l

M
~
0
M
.....
If')

0-
N
.....
'-0
00
roi
r-
.....
00
roi
.....
00
r-
N
~
'-0
roi
0
::::
~
roi
.....
M
N
roi
::::
r'"l
roi
.....
roi
roi
.....
N

N
" -
.....
M
r'"l
roi
~
roi
roi
~- ".

Figura 3 Medias movi/es de las precipitaciones diar-ias en la cuenca de Sainani y

1 caudales turbinados en Sainani

1 1,0

0,9 1
--Déstockage
1,8

1,6
- - Débit turbiné à Saimilli

1 ,...., 0,7
0,8

u
CJl
1,2

1 ~

~
~
..loi:
u
B
'"
0,6
0,5 1,0

... 0,4 0,8

1
i
::1
0
'.:::" 0,3
E 0,6
E
0,2

1 0,1

0,0
0,4

0,2

1 -0,1
3/5/96 18/5/% 2/6/96 17/6/96 2/7/96 17/7/96 1/8/% 1618/96 31/8/96 15/9/96 30/9/96 15/10/96
. f· _. ,. 0,0

1 Figura 4 Desalmacenamiento de las represas dU/'ante la estacion seca 1996 y caudal

turbinado en la planta hidroeléctrica Sainani

1
1
1
1
1
1 o • A A, ._''''. ."' ........... , .- ...... " . .. ',-
........... " '" "',

','
. ,
."
'
~

.
:",
..," , . .
;
,
,
'
70"10
60"/0
2 ,.
: 50"10
1
'

-
,g'"
c
4
-.
, ' "
.'
"

:!
40"10

!:! 30"10
ï5.

1 'u
.~
E:
...
:;l
6
20"10
10"10
o ................................................................................ ~, .. 0%
i
1 E
Il

10
....... Ilrécipitations moycnnes sur 4 jours
-10%
-20%
- - - Déstockagc Cil % du débit

1 12
1

3/5/96 1815/96 2/6/96 17/6/96 2/7/96 17/7/96 1/8/96 1618/96 31/8/96 15/9/96 30/9/96 15/10/96
-30%
-40% •

1 Figura 5 Variacion de almacenamienlo de agua en las represas en porcenlaje dei

caudal turbinado en Sainani y medias movi/es de las precipitaciones de 4 dias durante la
estacion seca 1996
1
40
1
1200 Rctcnucs plt"ines .....

35
1000

1 Ê
g
800
, - - Précipitations cumu lécs
'"co
1 .~
ï5.
'u
600 .
--Stockage cumulé

.~ 400

1 200

1 0
o.D o.D o.D o.D o.D r-- r-- r-- r-- r-- r--
~ ~ ~ ~ ~ 0'0
-.. 0'0
-.. 0'0
-.. ~ ~ ~

1
0 N N N N M
::::: ::::: -.. ::::: -..
r--
-..
N
-..
o.D
-.. !:!
00
-..
M
M
""
N
00
N
00 M
N
N N N

Figura 6 Precipitaciones medias en la cuenca de Sainani y almacenamiento de agua

1 de depositos, valores acumulados

1
1
1
1
 1,3
o' ""
.,. '"" ..
,
, , ,"
:..
.......
.. ... 1
.
l''" . . . .',
1,2
l' ......
......
" ,

:
.
1 ::0."
. " .. :
"
" ,
,
.'
"
,',
.. '" l.l
5 ..... . " " 1,0
",' .
.'
~.
,
',0'
:. ... "

. 0,9
. '
. ,,,"
,",

.
t' ....
,
,
...., ','
.. 0,8
~IO
o
"

"
,
"
,
'" ., 0,7 1
"
.~ ',f'

0,6
ï5.
'ü 0,5
, .~ 15
.=::
... 0,4
::l
o 0,3
'fE 20 0.2
0,1
0,0
25 ..... ,' Précipitations moyennes sur 4 -0,1
jours
- - - Stockage -0,2
-0,3
30 1

24/10/968/11/9623/11/968/12/9623/12/96 7/1/97 22/1/97 6/2/97 21/2/97 8/3/97 23/3/97 7/4/97 22/4/97

Figura 7 lv/edias movi/es de las precipitaciones de 4 dias y almacenamiento de agua

en las represas

30% 1

III 25% 1
c:
0
.~ 20%
ï5.
'ü /5%
...0..
'Q)

III
Q) 1()o1o
"0
-;g,
0 5%
1
l,
0% ..... ..... ..... ..... ..... ..... ..... .....
--
\0 \0 \0 \0 \0
~ ~ ~ ~ ~ ~ ~ ~ -~
~
~
t:'
~
t:'
~
:t
'"
~
0 N N
::::: ~ .....
::::: ..... N N
:;:r
N
Do --
<'l
N
Do <'"l
N
N
\0
N
00 <'l
N N

dUle

Figura 8 Almacenamiento de agua en los depositos en porcentaje des las medias

movi/es de las precipitaciones de 4 dias

1
1
Discusi6n y conclusiones respecto a la modelizaci6n hidrol6gica
1 2.3.1. Dcsalmacenamicnto en estaci6n scca
La influencia deI desalmacenamiento de las represas en estaci6n seca es muy importante,
1 como se puede observar en la Figura 5. La parte deI caudal turbinado en Sainani, que
proviene deI desalmacenamiento, comienza a alcanzar valores superiores a 30% a
comienzos de junio. Luego, la influencia aumenta hasta un maximo de 70% a mediados de
1 agosto y disminuye rapidamente después. Pero sigue siendo importante hasta fines de
octubre. El maximo es alcanzado después dei periodo sin lluvia mas largo dei verano. Se
observa que después de los fuertes eventos de precipitaciones la influencia dei
1 desalmacenamiento disminuye (Figura 5). Es el caso después dei evento de mediados de
julio y sobre todo después de los fuertes eventos de fin de agosto y de fin de septiembre. La
Figura 4 muestra que esta disminuci6n de la intluencia durante las crecidas no se debe
1 solamente al aumento dei caudal. El desalmacenamiento total disminuye durante las
crecidas asociadas a estas precipitaciones. Estas observaciones van de acuerdo con la
estrategia de la gesti6n de las represas (cf. 2.2.1).
1 La disminuci6n dei desalmacenamiento durante las crecidas podria hacer pensaI' que la

1 toma en consideraci6n dei desalmacenamiento durante una modelizaci6n hidrol6gica no es

tan importante. De hecho, durante los fuertes eventos de fin de agosto y de fin de
septiembre la influencia es inferior a 10%. La conclusi6n es seductora pero falsa. Aunque

1 es posible que la influencia dei desalmacenamiento tenga poca importancia durante las
puntas de las crecidas, no 10 es durante las otras fases. El procedimiento de
desalmacenamiento paralelamente a la subida de la crecida la hace menos aguda y aumenta
1 el volumen pOl' comparaci6n con una crecida no influenciada. Ademas, casi no hay
disminuci6n dei desalmacenamiento durante las débiles crecidas, que se espera poder
modelizar. El desalmacenamiento de las represas tiene pOl' 10 tanto una influencia sobre los
1 caudales en Sainani que no es insignificante para la modelizaci6n hidrol6gica.

Los aportes hidricos de las represas a los caudales de Sainani son variables. Expresado en
1 porcentaje dei caudal turbinado en Sainani, el desalmacenamiento deI lago Zongo varia
entre 15 y 25% Y el de los lagos Tiquimani entre 20 y 25%. Los lagos Tiquimani son
vaciados uno después dei otro, de manera que contribuyen a los caudales durante toda la
1 estaci6n seca, como en el caso dei lago Zongo. Los aportes de los otros lagos son menos
importantes. La influencia de la represa Mamankota, que alimenta al lago Zongo pOl' el

1 canal Prado, es normalmente inferior a 5%, pero aIcanza un mâximo de 10% durante
algunos dias en agosto. Los lagos Livinosa y Sankayuni contribuyen con aproximadamente
un 5% a los caudales, con maximos de 22% y 12%. La contribuci6n de San Pedro no pasa

1 de 2% (ver Anexo). Esto confirma que se pueden ignoraI' los cambios de stock de las
pequefias represas como San Pedro, muchas de las cuales no son consideradas en este
estudio preliminar (cf. 2.2.3).

1 Los aportes maximos dei lago Sankayuni tienen lugar en octubre, cuando los caudales
turbinados en Sainani estan cerca dei limite de capacidad (Figura 3), 10 que es muy extrafio
1 pues segun las informaciones de la COBEE las valvulas de esta represa ya no eran
maniobradas en octubre 1996 (Cuadro 3). Como muehos valores de nivel dei lago
Sankayuni son interpolados fuera dei periodo de operaci6n, basta que algunas lecturas sean
1 falsas para llegar a estos desalmacenamientos diarios extremos. Dado que las valvulas ya
no estân en este periodo, es mas probable que estos valores sean err6neos. Si es el caso, los

1
1
grandes aportes totales de desalmacenamiento en octubre (Figuras 4 y 5) serian también
1 sobrestimados. En cuanlo al almacenamiento duranle· los dias anteriores, se trata
probablemente de un artefaclo causado por los mismos errores de lecturas.

1 La contribuci6n de la represa Alto Viscanchani es muy irregular. Alcanza hasta un 15% de

los caudales turbinados en Sainani. Es sorprendente para una represa relativamente
pequcfia. Se puede observar, en la figura en Anexo, que aportes positivos y negativos se
1 siguen muchas veces durante la estaci6n seca, 10 que no es el casa de las otras represas. De
hecho, en 1996 la represa Alto Viscanchani era alimentada por las aguas de los lagos

1 Tiquimani y Guaraguarani; dej6 de serlo desde la puesta en servicio de la planta

hidroeléctrica Tiquimani en febrero de 1997 (cf. 2.1.1). Hasta ese momento, servia de
cuenca de regulaci6n diaria, coma las pequefias reprcsas dei rio principal. Tai regulaci6n es

1 necesaria porque las valvulas de los lagos Tiquimani y Guaraguarani s610 son maniobradas
cada dos dias (cf. 2.2.2)7. Probablemente es por eso que los aportes relativamente
importantes de este pequefio lago sean artificiales a causa de las diferencias de horas de

1 lectura (cf. 2.2.3). Sea 10 que sea, en 1997 las valvulas dellago Alto Viscanchani ya no son
tocadas, y coma la modelizaci6n hidrol6gica sera realizada con una cr6nica de la estaci6n
seca 1997, no nos preocuparemos.
1 Como conclusi6n de 10 expuesto, se puede decir que son las represas Zongo, Hatilata,
Taipicota, Guaraguarani, Livinosa, Mamankota y eventualmente Sankayuni cuyo
1 desalmacenamiento en estaci6n seca debe ser tomado en cuenta en una modelizaci6n
hidraulica.

1 2.3.2. Almacenamiento en cstacion de lluvias

Las entradas de agua en las represas provienen de las precipitaciones y dei deshielo de los
glaciares. No se conoce la importancia respectiva de estos dos origenes difcrentes. En la
1 Figura 7 se observa que la variaci6n dei almaccnamiento diario sigue a la de las
precipitaciones diarias (promediadas en 4 dias). El deshielo de los glaciares disminuye

1 durante los eventos de precipitaciones a causa de la baja de energia consecutiva a la

intercepci6n de la radiaci6n por las nubes, la baja de las temperaturas y el aumento dei
albedo de los glaciares (cf. RIGAUDIERE, P. et al.). En -consecuencia, las variaciones

1 paralelas dei almacenamiento y de las precipitaciones muestran que una gran parte dei agua
almacenada durante las precipitaciones proviene de las precipitaciones. No obstante, se
observa que en noviembre de 1996 el almacenamiento es mayor en comparaci6n a las
1 precipitaciones que en pleno invierno (Figura 7). En la figura de los valores acumulados
(Figura 6) este hecho es atm mas marcado. La curva de almacenamiento acumulado cambi6
bruscamente de inclinaci6n, mientras que la inclinaci6n de las precipitaciones acumuladas
1 no cambia mucho. Se constata que este evento se produjo entre el 8 y el 25 de noviembre
de 1996. Durante este periodo, los aportes dei deshielo de los glaciares en el llenado de las
represas deben haber sido importantes. Justamente, en el mes de noviembre el deshielo de
1 los glaciares es generalmente el mas elevado. Asimismo, entre el 8 y el 25 de noviembre
los valores dei almacenamiento diario en porcentaje de las precipitaciones diarias son los
mas altos (Figura 8). No obstante, se debe ser prudente si se quiere atribuir los picos de la
1 Figura 8 a eventos de fuerte deshielo de los glaciares. Un pico asi puede ser también el
efecto de un escurrimiento de base que continua llenando las represas cuando las

1 7 Hay en dia, la cuenca de carga de la hidroeléctrica Tiquimani tiene esta funci6n de regulaci6n. Por esta
raz6n su valumen es excepcionalmente grande (30'000 ml).
1
1
precipitaciém ya cesaron s, 0 de una reparticion desequilibrada de las precipitaciones en las
1 partes altas de la cuenca de drenaje, que son drenadas en las represas.

Hacia fines de febrero de 1997 las represas Zongo y Taipicota alcanzaron sus niveles
1 mâximos. La represa Hatilata se Heno a fines de marzo. En marzo, hay un corto periodo de
desalmacenamiento que es seguido de dos picos de almacenamiento, el primero de los
cuales (23 de marzo) es precedido por precipitaciones menos fuertes. En abri!, hay cinco
1 eventos muy pequefios de almacenamiento; el primero (4 abril) y el ultimo (26 abril) siguen
a dias sin precipitaciones (Figura 7), de donde resultan los dos grandes picos, en abril, de
desalmacenamiento en porcentaje de las precipitaciones (Figura 8). El mes de abril es el
1 segundo periodo de deshielo elevado en el ano hidrologico, y el almacenamiento durante
estas dos fechas se debe principalmente a los aportes de agua de los glaciares.

1 _ Précipitations spatiales sur le bassin versant jusqu'à )'éxutoire de Sainani

1 20%
--%des précipitations sur Ic bassin vcrsantjusqu'à l'éxutoirc dc Tiquimani
0
18% "fl~
1 16%
14%
5

10

1 12%
10010
15

20
8%
1 6% i
25
4%

1 2%
0010
moyenne
-- l- 1·
30

35
-.0 t- t- t- t- t- t- t-
a.. a.. a.. a..

1 a..
-. -. a..
-. a..
-. -. a..
-. -.
~ N ~
...., v v
-.
...., -.
-.0
-.
on 0 N -.
-.0
N 00 M t::: N t::: M N N N N
1 N N N date
1. . __ ~ __ .... ... . . . .... _

1 Figura 9 Parte de las precipitaciones medias en la cuenca de Sainani caida en la

cuenca de Tiquimani

1
1
1
1
1
8 Es justamente para disminuir este efecto que se ha promediado las precipitaciones diarias de 4 dfas.

1
1
En cuanto a la modelizaci6n hidrol6gica, es importante saber si hay que tomar en cuenta el
1 almacenamiento de agua en las represas durante la estaci6n de lluvias. Mas precisamente,
si hay que tomar en cuenta el almacenamiento deI escurrimiento de las precipitaciones y de
la fusi6n de nieve. El volumen total de agua almacenada durante la estaci6n de lluvias 1996
1 s610 alcanza un 3.9% de las precipitaciones acumuladas dm"ante este periodo (Figura 6). En
caso de que la parte deI escurrimiento de las precipitaciones (escurrimiento directo y

1 deshielo de nieve) almacenada en las represas sea también minima, se la puede dejar de
lado durante una modelizaci6n hidrol6gica, como se vera mas tarde (cf. 3.2.1). Este podria
ser el" caso si los picos en la Figura 8 se debieran al aporte de agua de deshie10 de los

1 glaciares mas que a precipitaciones particularmente abundantes en las partes altas de la

cuenca de drenaje que son drenadas en las represas.

1 La Figura 9 muestra la parte de las precipitaciones medias diarias en la cuenca de Sainani

que cay6 sobre la cuenca de drenaje Tiquimani. La cuenca de drenaje de Tiquimani
contiene todas las [Link] que pueden contribuir al llenado de las represas mas

1 importantes. Las precipitaciones mas fuertes alcanzan las zonas altas y aportan
generalmente un mayor volumen de agua a las partes altas de la cuenca. El valor medio de
la parte de las precipitaciones sobre las partes altas seria de 8.8% en estaci6n de lluvias
1 1996, el valor maximo, 12.7% Y el valor minimo 0%. Esto significa que la parte deI
escurrimiento proveniente de las precipitaciones que ha sido almacenada durante este
evento maximo es al maximo 44% mas grande que la parte almacenada en promedio
1 durante toda la estaci6n9 . Es decir que la parte de los escurrimientos acumulados de las
precipitaciones que fue almacenada en las represas durante la estaci6n de lluvias 1996 es
deI 4%, la parte almacenada de los escurrimientos diarios no pudo pasar de 5.8%.
1 No obstante, a fin de saber qué parte deI escurrimiento de las precipitaciones fue
almacenada en las represas en estaci6n de lluvias 1996, habra que tomar en cuenta las
1 pérdidas de agua (evapotranspiraci6n y sublimaci6n). Para el desarrollo del método de
reconstituci6n de los caudales en estaci6n de lluvias (cf. 3.2:-1), se parte dc la hip6tesis que
esta parte deI escurrimiento es insignificante, sin demostrarJa lo .
1
1
1
1
1
1
1 944% = (12.7%-8.8%) / 8.8%.
10 Se podrâ veriticar esta hip6tcsis mas facilmentc cuando se tengan mediciones dei caudal en estaci6n de
lIuvias en Sainani.
1
1
3. Un método para reconstituir los caudales naturales
1 provenientcs dei deshielo de los glaciares y de los escurrimientos
de las vertientes
1 3.1. Principio dei método de reconstitucion de los caudales naturales

1 3.3.1. La estimacion de las caudales Ilaturales sin aporte de los glaciarcs

Un modelo hidrol6gico perfeccionado utiliza en entrada las precipitaciones y proporciona
en salida los caudales en el cierre de la cuenca de drenaje. Los caudales no deben ser
1 influenciados por otros datos diarios no tomados en cuenta en el modelo. Es por esto que se
debe hacer abstracci6n, en los caudales medidos en Sainani,-de la influencia dei cambio de
stock de agua en las represas y de aportes de agua de deshielo de los glaciares.
1 En una primera etapa, se intenta utilizar los datos disponibles sobre la variaci6n de los
niveles de agua en las represas y en los caudales turbinados en las plantas hidroe1éctricas
1 para estimar los caudales na/urales. Los caudales provienen ûnicamente de las
precipitaciones, de la fusi6n de nieve y dei deshielo de los glaciares al interior de la cuenca

1 de drenaje. En ese sentido, el desalmacenamiento de las represas en estaci6n seca significa

un aporte exterior de agua que hay que sustraer de los caudales totales medidos. En
cambio, su llenado en estaci6n de lluvias es una sustracci6n de agua deI sistema que hay

1 que sumar a los caudales medidos.

Càrlte:;daicsdi ~io Qulie:;daicstutinrl:B

1 ..._di~Zaw
--~--

. mSinri
1

I--=~:~=
Rniénàtua/munn
1
enfiniénœSl tqxgafia Plrnwurie1d •
1 • ds:l1rnnnuicrto
œ1a;n:peiE

1 [Link]:; ~ei,
~ T •
Càrlt œ1a; \t11Î61es
àtuayeqrsicién ~é4Dtfs œ1a;lqmxE
1 œlamu.1a; yœ la; gme;

1 [Link] ~
l\1J..U~~
ll'1A&1U._~ ....
1

m\a;JlItfsm~ 1
1

1 œlaaJmlœ~e

Figura la Principio de la modelizaci6n hidrol6gica dei valle dei Zongo

1 La modelizaci6n precipitaci6n!caudales no toma en cuenta el deshielo de los glaciares, sin

1 embargo constituye un aporte anual. La sustracci6n de agua de precipitaciones en las zonas

de acumulaci6n de los glaciares tampoco es tomada en cuenta. Las superficies glaciares

1
1 deben ser tratatas aparte y se considera el agua de deshielo proveniente de esas superficies

1 como un aporte exterior de agua. El principio de la modelizaci6n hidrol6gica deI valle deI
Zongo se presenta en la Figura 10.

1
3.1.2. La sepraci6n de los periodos de modelizaci6n

1 En los capitulos anteriores, se vio que la influencia deI sistema artificial, es decir deI
cambio de stock en las represas, no es la misma en estaci6n de lluvias y en estaci6n seca.
Por 10 tanto hay que tomarla en cuenta de una manera diferente durante las dos estaciones

1 (cf. 2.3.1 y 2.3.2). De hecho, hay que distinguir cuatro periodos diferentes durante el afio.
Primero se tiene la estaci6n de lluvias y la estaci6n seca, que contiene el periodo de
operaci6n de retenidas. Enseguida, hay que distinguir los dos periodos de transici6n entre
1 las estaciones. El mas importante cs el fin de la estaci6n de lluvias.

En 1997, este periodo cubria el mes de abril (Figura 3, 2.2.5), cuando las precipitaciones ya
1 habian disminuido, pero los caudales continuaban siendo importantes, es decir superiores a
la capacidad maxima de la planta hidroeléctrica de Sainani(ct: 2.2.4), bajo la influencia deI
deshielo elevado de los glaciares. Durante este periodo, solamente la represa Zongo es
1 controlada y los dias de almacenamiento y de desahnacenamiento se siguen (Figura 7). No
obstante, el desalmaccnamiento durante este periodo de transici6n no parece ser
importante. Al final de la estaci6n de Buvias en 1996 y 1997 constituia aproximadamente el
1 2% de los caudales turbinados en Sainani, con maximas de 8.7% en 1996 y de 8.0% en
1997 (resultados que no se muestran en las figuras). Si se considera que los caudales

1 turbinados durante este periodo de transici6n estân al limite de capacidad de la planta

hidroeléctrica de Sainani, y por consecuencia son inferiores a los caudales totales, la
influencia real del desalmacenamiento es atm nuis débil. Entonces se puede suponer que es

1 insignificante.

Durante el periodo de transici6n entre la estaci6n de l1uvias y la estaci6n seca en 1996 se

1 observan algunos dias de almacenamiento alto alrededor del 30 de septiembre (Figuras 4 y

5, 2.2.1) después de las fuertes precipitaciones. El almacenameinto alcanza el 30% deI
caudal turbinado en Sainani. Es evidente que hubo almacenamiento durante esos dias, pero
1 el porcentaje parece exagerado, tanto asi que durante este periodo el deshielo de glaciares
es minimo. Se trata muy probablemente de un arteficio causado por errores de lectura (cf.
2.3.1). Sin embargo, coma no se toma en cuenta la parte deI almacenamiento que proviene
1 deI escurrimiento de las precipitaciones, y coma el deshielo de los glaciares no es
importante durante este periode de transici6n, se puede utilizar el método para la
reconstituci6n de los caudales de estaci6n seca (cf. 3.2.1).
1 Basta entonces con diferenciar dosperiodos de modelizaci6n. El primera es de la estaci6n
seca. Comienza cuando el nivel de la represa Zongo empieza a descender de manera
1 continua, después deI periodo de transici6n, y se termina al final de la estaci6n seca,
cuando las valvulas son cerradas para Benar las represas- y cuando los valores deI

1 almacenamiento se vuelven definitivamente positivos. En 1996, esto sucedi6 el 24 de

octubre (Cuadro 2). Durante el resto dei ano, el desalmacenamiento puede no tomarse en
cuenta y es el método que se utiliza para la reconstituci6n de los caudales naturales en

1 estaci6n de lluvias (cf. 3.2.1).

1
1
3.1.3. El problema de la diferencia horaria de los datos
1 Para facilitar el calculo deI cambio de stock de agua total en las represas, no se ha tomado
en cuenta las diferencias de las horas de lectura de los niveles de las represas y de los
caudales turbinados en Sainani (cf. 2.2.3 y Cuadro 3, 2.2.2). A fin de obtener una primera
1 estimaci6n de la influencia deI cambio de stock en los caudales, esta simplificaci6n era
aceptable. Sin embargo, en una modelizaci6n hidrol6gica precipitaciones / caudales durante
la estaci6n seca, el defase de las horas de lectura ya no es insignificante, sabiendo que las
1 lecturas son repartidas en casi todo un dia (Cuadro 3, 2.2.2).

Lugar Tiempo de Hora de lectura Lectura 6ptima en Diferencia

1 transferencia a la
planta hidroeléctrica
(nive 1, caudal) relaci6n al lago
Zongo,
r

Sainani 16.00 h

1 Zongo
Mamankota
7 h 35 min.
7 h 55 min.
16h00
09h30
16h00
14h40
-
- 5h 10 min.
LivÎnosa 5 h 10 min. 19h30 18h25 + 1 h 5 min.

1 Sankayuni
Hatilala
Taipicota
3hlOmÎn.
3 h 50 min.
4h
20h30
08h00-09h00
08h30-09h30
20h25
19h45
19h55
+ 5 min.
-9 h 45 (-lOh 45)
-lOh 25 (-llh 25)

1 Guaraguaran i
Hidroel. Tiquimani
Hidroel. Santa Rosa
3 h 10 min.
2 h 50 min.
10min.
11h00-12h00
21h00
23h30
20h25
20h45
23h25
-811 25 (-911 25)
+ 15 min.
- 5 min.
Hidroel. Sainani - 23h30 23h35 + 5 min.
1 Cuadro 4 Tiempo de transferencia entre las represas (las plantas hidroeléctricas) y la
planta hidroeléctrica Sainani, horas de lectura, horas Ideales de lecfura en relaelon a la
1 lectura dei nivel de la represa Zongo

1 Con ayuda de los tiempos de transferencia entre las represas y la planta hidroeléctrica
Sainani se puede calcular las horas de lectura ideales en relacién a una hora de lectura dada
en Sainani (0 en otra parte), Para minimizar las diferencias entre estas horas ideales y las

1 horas de lectura reales, se elegira la segunda lectura dei dia en el caso en que haya dos pOl'
dia. Los resultados se presentan en el Cuadro 4. Se ha calculado la diferencia en relaci6n a
la segunda lectura de nivel de la represa Zongo a las 16:00 h. Se ve que en los casos de las
1 represas Livinosa y Sankayuni las diferencias entre la hora ideal y real son minimas. No es
coincidencia. Las lecturas se hacen durante la apertura y cierre de las valvulas de estas
represas. Como se intenta satisfacer las horas de demanda de punta de la tarde con su
1 operaci6n, el interés es que el agua soltada llegue durante estas horas a las plantas
hidroeléctricas. Para el lago Zongo, la mayoria dei agua es también soltada antes de las
17:00 h, de manera que llegue a las horas punta a las plantas hidroeléctricas.
1 Lamentablemente, no es el caso de las represas Hatilata, Taipicota (lagos Tiquimani) y
Guaraguarani. Como estas represas son localizadas a medio dia de caminata, sus valvulas
s610 son reguladas cada dos dias (Cuadro 3, 2.2.2). Alli es donde se encuentran las
1 diferencias entre las horas de lectura ideaJ y real mas importantes; estas alcanzan la mitad
de un dia (Cuadro 4). La influencia sobre los caudales de estas represas es demasiado

1 importante para no dar importancia a estas diferencias.

La soluci6n para la modelizaci6n en estaci6n seca consiste en reemplazar el cambio de

1 stock diario de estas tres represas pOl' los caudales diarios turbinados en la planta
hidroeléctrica Tiquimani. Ademas deI desalmacenamiento de las represas Hatilata,
Taipicota y Guaraguarani, estos caudales contienen una parte de los escurrimientos de las
1
1 vertientes que atraviesan los canales de derivaci6n. Mientras estos aportes

1 complementarios de agua provengan de las vertientes deI valle Tiquimani, es decir deI
exterior de la cuenca de drenaje deI Zongo, la aproximaci6n elegida sera mejor que si s610
se considerara el cambio de stock, ya que se considera la totalidad de los aportes de agua

1 dei exterior dei sistema (cf. 3.1.1). En cambio, al exterior de la cuenca de drenaje deI
Zongo, los escurrimientos que lIegan a los canales l'nlsean la estimaci6n de los aportes
externos. Felizmente, al interior de la cuenca de drenaje el agua es principalmente

1 transferida dentro de tuneles. Cerca de la cuenca de regulacl6n de la planta hidroeléctrica

Tiquimani, hay dos tomas de agua en afluentes deI Zongo que estan activas en estaci6n de
lluvias. Las subcuencas captadas tienen una superficie de aproximadamente 3 km 2, que es

1 poco en relaci6n a la superficie total de la cuenca de drenaje de 237 km 2• Las pérdidas en

los canales y tuneles de toma de la planta hidroeléctrica Tiquimani son minimas, es decir
que casi todo el aporte externo de agua es turbinado.
1 Otra diferencia horaria importante de lectura de nivel se refiere al lago Mamankota. Esta
represa de dimensi6n mediana tiene una intluencia sobre los caudales en Sainani que esta
1 normalmente por debajo de 5%, pero que puede a1canzar el 10% durante algunos dias deI
ano (cf. 2.3.1 y Anexo). Como las valvulas dellago Mamankota s610 son maniobradas cada
dos dias, el caudal que sale no cambia bruscamente. Dada la pequena influencia de la
1 represa, es admisible una interpolaci6n lineal deI volumen almacenado entre dos lecturas
de niveles para obtener el volumen almacenado a la hora ideal.

1 3.1.4. El aporte hidrico deI valle Coscapa

Una parte deI escurrimiento dei valle Coscapa es captada por tomas de agua en el rio

1 Coscapa y en dos de sus afluentes. El agua es llevada a la planta hidroeléctrica Santa Rosa
II, mas arriba de Sainani (cf. 2.2.4). Como el rio Coscapa se une al rio Zongo mas abajo de
Sainani, este apOlte hidrico proviene deI exterior de la cuenca de drena:Îe en Sainani. Si

1 todo caudal es captado en las tres tomas, los aportes de agua de deshielo de todos los
gla~iares dei valle dei Zongo pasan a Sainani. Es casi el caso en estaci6n seca, cuando los
caudales que rebalsan de las tomas son muy débiles y probablemente insignificantes en

1 comparaci6n con los caudales captados. En estaci6n de Iluvias, la situaci6n es inversa, el

desagüe excede con mucho 10 que es captado. La cuesti6n que se plantea es saber si se
quiere inc1uir estos aportes de agua externos en la modelizaci6n en estaci6n seca. Si se los
1 incluye, hay que tomar en cuenta el desalmacenamiento de la represa Sankayuni que se
encuentra aguas arriba de las tomas. Esto tendria la ventaja de efectuar la modelizaci6n
hidrol6gica tomando en cuenta aportes de todas las superficies glaciares. Pero el problema
1 de esta inclusi6n es que los limites de cuenca de drenaje son vagos durante las pérdidas en
los canales, 0 los desagües de las tomas de agua, no son insignificantes. Si no se los
inc1uye, hay que sustraer de los caudales en Sainani los caudales que son transferidos al
1 Zongo par la planta hidroeléctrica Santa Rosa II. Como toda el agua que es turbinada en
Santa Rosa II proviene deI exterior de la cuenca de drenaje, esto tiene la ventaja de una
delimitaci6n exacta de la cuenca de drenaje.
1
1
1
1
1
Los métodos para la reconstitucion de los caudales naturales
1 3.2.1. El método para la reconstitucion de los caudales naturales en estacion seca
En estacion seca hay que tener en cuenta el desalmacenamiento de las represas. Los
1 cambios de nivel de las represas Zongo, Mamankota, Livinosa y eventualmente Sankayuni
son directamente utilizados en el método, mientras que los desalmacenamientos de los
lagos Tiquimani y Guaraguarani son asimilados a los caudales turbinados en la planta
1 hidroeléctrica Tiquimani. Si no se incIuye la modelizacion hidrologica las zonas deI valle
de Coscapa que son drenadas por los canales de toma de la planta hidroeléctrica Santa Rosa
II, hay que tomar en cuenta los aportes externos de esta planta hidroeléetrica en los
1 caudales dei Zongo.

1 El esquema de los flujos de agua (Figura Il) muestra como se pueden estimar los caudales
en la salida de Sainani que provienen unicamente dei escurrimiento de las precipitaciones y
dei deshielo de la nieve (cf. 3.1.1). Estos caudales corresponden a los flujos de agua (2) en
el esquema ll . Si se parte dei caudal turbinado en la planta hidroeléctrica Sainani (8), que
1 corresponde al caudal total dei Zongo en estacion seca, se ve que esta compuesto de otros
flujos de agua, dos de los cuales provienen dei exterior de la cuenca de drenaje. Son los

1 aportes de agua de las cuencas externas Tiquimani y Coscapa. El caudal turbinado en Santa
Rosa II (6) corresponde exactamente al aporte externo, mientras que el caudal turbinado en
Tiquimani (7b) contiene una pequeiia parte de agua proveniente dei exterior de la cuenca
1 de drenaje (3a) (cf. 3.1.3 y 2.1.1). Si se sustrae (6) y (7b) de (8), se obtiene una buena
estimacion dei caudal proveniente dei interior dei sistema (la suma de los flujos 3b, 4a, 4b
5 Y7a). El otro componente importante dei caudal turbinado en Sainani (8) es el agua dei
1 desalmacenamiento de las represas al interior de la cuenca de drenaje (5). Este flujo no es
conocido, ya que los cambios dei stock de agua incIuyen también las entradas en las
represas, que comprenden los aportes de los glaciares (4c), de las precipitaciones y de la
1 fusion de nieve (3e). El cambio dei stock de agua es la diferencia entre el flujo (5) y los dos
flujos (3c) y (4c). En estacion seca, estos ultimos son débiles en relacion al flujo (5). Si se
sustrae este cambio de stock de los caudales que provienen dei interior dei sistema, se
1 obtiene el escurrimiento natural compuesto dei agua de deshielo de los glaciares (4a, b, c)
y de los escurrimientos de las precipitaciones y de la fusion de nieve (3b,c). La suma de

1 estos flujos corresponde al escurrimiento natural de la cuenca de drenaje (suma de los

flujos 1 y 2) disminuido por los aportes interiores en los caudales turbinados en la planta
hidroeléctrica Tiquimani (3a) y aumentado por las pérdidas de los canales de derivacion de

1 esta planta hidroeléctrica al interior de la cuenca de drenaje (7a). Como estos dos ultimos
flujos son muy débiles en relacion al total, se obtiene una buena estimacion dei
escurrimiento natural. Basta con sustraer los aportes de los glaciares (1) para obtener el

1 flujo buscado, el escurrimiento de las preeipitaciones y de la fusion de nieve de la cuenca

de drenaje (2). Resumiendo :

1 (2) = (3a,bc,) - (3b,c) + (7a) = [(3b) + (4a,b) + (7a) + (5)] - [(5) - (3c) - (4c)] - (4a,b,c)

donde: [(4a,b) + (3b) + (7a) + (5)] = [(8) - (6) - Ob)] = Diferencias de los caudales
1 turbinados

1 Il Se parte de la hip6tesis, para la modelizaci6n hidrol6gica, que no hay pérdidas de agua durante la
transferencia en el rio y en los canales. En la realidad, el ancho de las fechas de los f1ujos disminuye hacia
abajo a causa de estas pérdidas y la estimaci6n de los f1ujos desconocidos es mucho mas complicada.

1
1
1 [(5) - (3e) - (4e)] = Cambio de stock de las represas
(4a,b,e) = (1) = Aguas de deshielo de los glaciares

1 El esquema detallado deI método para la reconstitucién de los caudales naturales en

estacién seca se presenta en la Figura 12.

1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
Agua originaria dei deshielo Agua originaria dei

1 de glaeiares, las precipitaciones

y la nieve
esnlelo de los glaciales
...........................! ~)

........................... .: .

1
Escurrimiento al exterior de la cuenca Escurrimiento al interior de la cuenca

1 Sub-euenca dei Sub-cuenca dc Sub-cuenca de Partcs de la cuenca Partcs de la euenca que

1
valle Coscapa
(Iago Sankayuni),
al exlerior de la
cuenca Sainani
lagos Tiqllimani.
al cXlcrior dei
valle dei
Zongo
al exterior dcl
vallée du Zongo
.
lago Guaraguarar i. que no son drenadas
en las represas
:
. . . !.... ..
3~L ..1.... • 3~) .. . :-
:
no son drenadas en las
represas Zongo, Maman-
kola, Livinosa (Sankayuni)

!
!...
:::
:
1 13 c)
.;. •.•
1 ••.. ....... ~ ••.
····················r··· : :
1 !! , 4a)
.
1 ,:.!~b~:::" J,o.
Prcsas Coscapa, Lago Guaraguara i
i •
1 Tributario 1 ct Il
Agua originari'
l;l

Agua originaria ! 1
Mamankota, Livi-
nosa (Sankayuni)
dei déstockage dei déstockage
1
1 ~
1

Il
Agua originaria
dei déstockage

1
1
Canales Cos-
capa, Tribu-
tario r et Il
Canal TiqUimanii

G",rngo'rn~ \.V
1

j 1
Rcbalsc de 10

1::~~~1 ••
::::::

1 1---
1 c~" Gor'~:o~:=~·:_1 i, f. i : <::~,v:~~:'
T;qo;m'~~ c,'!", d, :do"~.o
Rebalse
1 R;o
I(ruera
!ii.! 1 1 Arno""",." ';0 Zoo,o, do 1" pl.o'" 1

:?v ! 1
~
system '"

1 PlantaIlSanta
Rosa 1
'" " ' .
"'''"."''".." ,, " ,....... """."""""""""" "..""""..,, ,, l:llll
:,,:.~, "..•

c:: ~~6) p'~1?~:~ L..\y5.~~~(~~~T~~;~-l'-·

1
1 1 Rio Zongo, agllas arriba planta Sainani

u
1

1
v
1 Rio Zongo,
aguas abajo de
1

L----.-
..- -..-..-..- -..-..- -
X
Rio Zongo, aguas abajo de planta Sainani
-)-------::....--=---....:.....-....:.....-)- - - - - - - - - - - - - '
1

Chururaqui = Flujo desconocido C :::>-x = f1ujo conocido

1
Figura Il Esquema de los flujos de agua en el valle dei Zongo en estacion seca. Las
1 explicaciones se encuentran en el texto.

1
-------------------
Figura 12 Esquema dei principio de reconstifucion de los caudales a partir de los caudales totales en estacion seca

Caudales diarios Caudales diarios Nivel dei agua Niveldel agua Niveldelagua Energia producida Energia producida
dei aftluente dei dei canal Prado, en el lago Ma- en el lago en el lago Zongo diariamente en diariamente en
glacial Zongo affluente dei glacial man kota, medido Livinosa , medido medido cada planta Sant Rosa Il, planta Tiquimani,
(exposiciôn sud-est Charquini (exposiciér Cada dos cada dia Dia a las 16:00 Medida a las 23:30 h Medida a las 21 :00 h
nor-oeste) y dei dias a las 9.301 20:30 h 1
lago Mamankota
Fotos satelitales [" ..1 1. 1 ,
De las partes glaciales Volumenes y Estimacion de la relacion nivel/volumen ~ 1
superficies max - simplificando la forma de los (agos coma ~
Rendimiento en funcion de la energia
De la cuenca vertiente
mas de los lago L . .
i~ ~.~.~~.~~ ~~:.~.~!.~~ ..1 L -_ _ ----:,p,o_r_u_n_id_a_d_d_e_ti_e_m-:-p_o_'_a-rlt_u_ra_d_e_la_ca_i_d

Malla dei MNT rC-;;b~de st~k--l

de superficies gla- r----- . . . ---... Idellago Mamankota 1
cia1es. parâmetros 1caudales dei afluente 1 I(cada dos dias) __ -' ,~
topogrâficos !2!:.!.~alci~ Ch~r~u~LJ ~::::.:I....:.............. rAj,;;;ï;sd~gu;ck1 Aport~&-;g~;P~lrAI

(inclinacion. exposici' ) 'U!.1.~:D??~~.~.~~.~..~.i.~:.~! ..! lias tomas en el 1 111: changement de 1

r~ï~d~ï~d~"d~~hi~ïo'd~I""l ~ca~bi;-d;~Ck----~ !~':i~x1:r~~~c~pa

!! ~::~i~a~~.I~:i~~~~~a, i
l conjunto de los glaciales i 1 diario (Mamankota) __ ..! pistema) 1 n Guaraguarani, 1
~ que se basa sobre dif- ~ ------r-- ---- .... ---- 1';' entradas naturales l
,· · · · . ·}· · ·.r· · · ,
i férencias topogrâficas ~ --~-- ...
C~~t~s~~,~~~;n~o~:n;~~v~o~~~~~I:
_ __ ...
J
1a los canaleS
:(s~::.~1exte~~-
rp~~sd;-agu;de deili~l~de gi~G:i~-ï
~':!:..--------...,------7-..
..j..(-)
1
(-) (-) (-)
..

r-------------------------------------------·-----------------------
1 Caudales de la cuenca vertiente hasta Sainani (sin los aportes de dos pequerias sub-cuencas cuyo escurrimiento no es
1 turbinado en Sainani), proveniente de las precipitaciones sobre las partes no glaciales y dei derretimiento de la nieve
--~-------------------------------------r---------- ---------------.

~~{~F~)i~:r:sfi;OS, ~-{.~·~i;j;.:~~~~:~!~~~.~i.l.~f..~.~;;;t~~~·~.:··.·.:::.·'.:: > c 2 e l o hydrologico :>< 1precipitaciones diarias

L...- ....JI = Série de datos primârios L- I= Dtros datos
r-----
I J= Série de datos tratados L.·.·.·.·.·.·.·.·.·.·.·.·.·.~·.·.·.·.·.·.·.·.]= Modelos y algoritmos para el tratamiento de los dé
1
3.2.2. El método para la reconstitucion de los caudales naturales en estacion de Iluvias
1 En estaci6n de Iluvias, la situaci6n es completamente diferepte a la estacion seca. Los
caudales dei Zongo y de sus afluentes son mucho mas importantes y la influencia deI sistema
artificial es mas débil (cf. 2.3.2). Las cuencas externas Tiquimani y Guaraguarani no aportan
1 agua al sistema sino a partir deI momento en que sus represas estan llenas. En cuanto a las
represas al interior de la cuenca de drenaje, no es su desalmacenamiento sino su llenado que
hay que tomar en cuenta. Una parte del escurrimiento producido en la cuenca de drenaje no
1 llega al cierre porque esta almacenada en las represas. Se plantea la hipotesis que la parte dei
agua que es utilizada para llenar las represas, y que proviene de las precipitaciones 0 de la
fusi6n de nieve, es insignificante en relacion a los escurrimientos totales que provienen de
1 estos dos origenes (cf. 2.3.2). En estacion de lluvias, no es posible elegir si se incluye 0 no, en
la modelizacion hidrologica, las partes dei valle Coscapa que son captadas pOl' los canales de
derivacion de la planta hidroeléctrica Santa Rosa II (cf. 3.1.4). El desagüe de las tomas es
1 dcmansiado importante.

La Figura 13 de los flujos de agua muestra como se pueden estimaI' los caudales al cierre de
1 Sainani que provienen unicamente de las precipitaciones y de la fusion de nieve (cf. 3.1.1).
Estos caudales corresponden a los flujos de agua (2) en el esquema. En estaci6n de lIuvias, se

1 debe disponer de mediciones de caudal total cn Sainani, pues los caudales turbinados no
tienen relacion con los caudales reales (cf. 2.2.4). Este caudal total corresponde al f1ujo (8) en
el esquema. Si se resta a éste los aportes externos dei valle de Coscapa, conocidos pOl' los

1 caudales turbinados en la planta hidroeléctrica Santa Rosa II (7), se obtiene el escurrimiento

proveniente dei interior de la cuenca de drenaje, con un pequefio aporte extemo de agua dei
sistema Tiquimani/Guaraguarani (6). Este aporte no se produce sino cuando los lagos
1 Tiquimani y Guaraguarani se desbordan, al final de la estaci6n de lIuvias. Se parte de la
hip6tesis que este flujo (6) es pequefio en comparaci6n al flujo total de agua considerando la
pequefia superficie dei sistema externo en relacion a la superficie total de la cuenca de drenaje.
1 Mientras que las represas Mamankota y Livinosa se llenan, no hay rebalse, es decir que los
escurrimientos de las subcuencas de estas represas (4a, 3a) son enteramente almaccnados y no
contribuyen a los caudales en Sainani. No es el caso de la represa Zongo, cuyas valvulas son
1 maniobradas durante todo el afio. El agua soltada corresponde al caudal turbinado en la planta
hidroeléctrica Zongo (5). Si se sustrae este flujo (5) dei escurrimiento proveniente dei interior
de la cuenca de drenaje ((8)-(7», se obtiene el escurrimiento natural de la cuenca de drenaje
1 menos la parte que es almacenada en las represas ((3a,b) y (4a,b».

1 En la hipotesis que los f1ujos de agua en las represas que provienen de las precipitaciones y de
la fusion de nieve (3a,b) son insignificantes (cf. arriba y 0), solo hay que tomar en cuenta la
parte dei agua de almacenamiento que proviene dei dcshielo de los glaciares (4a,b). Dicho de

1 otra manera, basta con sustraer de los escurrimientos naturales ((8)-(7)-(5» el escurrimiento
que proviene de los glaciares y que no esta almacenado en las represas (4c,d) para obtener el
escurrimiento des las precipitaciones y de la fusion de nieve. A diferencia de las otras

1 represas, el lago Zongo es alimentado pOl' canales. Como estos canales y las tomas de agua
que los alimentan pueden desaguar, no toda el agua de las subcuencas drenadas llega al lago.
A partir de las superficies glaciares, se puede estimaI' el total deI escurrimiento de agua de
1 deshielo en estas subcuencas «4b)+(4c», pero no la parte de 10 que llega allago Zongo. En la
hipotesis que el desagüe de agua de deshielo (4c) es insignificante, el escurrimiento de las
precipitaciones y de la fusion de nieve (2) es obtenido pOl' la diferencia entre los
1 escurrimientos naturales ((8)-(7)-(5» y el escurrimiento de las superficies glaciares que no es

1
1
1 transferido hacia las represas (4d)12. Los aportes de agua de deshielo de los glaciares deben ser
calculados en funci6n dei desbordamiento de las represas Livinosa y Mamankota. En el
momento que una de las dos represas, Livinosa 0 Mamankota, desagüe, es necesario aumentar

1 las superficies glaciares calculando su aporte a los caudales naturales (4d) por las superficies
glaciares de subcuenca drenada por esta represa (una parte dei flujo (4a». No obstante, las
superficies glaciares drenadas pOl' las represas Mamankota y Livinosa son bastante pequefias

1 en relaci6n a las superficies glaciares totales de la cuenca de drenaje para que los aportes
adicionales de agua de deshielo sean insigllificalltcs. Resumiendo:

1 1) Sin desbordamiento de los lagos Livinosa 0 Mamankota:

(2) = (3a,b,c) = (8) - (7) - (6) - (5) - (4d) - (4c) ~ (8) - (7) - (5) - (4d)
1 2) Con desbordamiento:

1 (2) = (3a,b.c.) = (8) - (7) - (6) - (5) - (4d) - (4c) (4a) ~ (8) - (7) - (5) - (4d) - (4a)

donde: (8) = Caudales totales medidos en Sainani

1 (7)
(5)
= Caudales turbinados en la planta hidroeléctrica Santa Rosa II
= Caudales turbinados en la planta hidroeléctrica Zongo
(4d) = Aguas de deshielo de los glaciares que no son drenadas por las
1 (4a)
represas
= Aguas de deshielo de los glaciares que son drenadas por las represas

1 Mamankota y/o Livinosa

1
1
1
1
1
1
1
12 Si esta hip6tesis es falsa, se podrian utilizar los caudales de los canales Prado, Tubo Vertcdero y Alpaca para

1 caJcular el aporte hidrico al lago ((4b+(3b». Agregando este flujo a los escurrimientos naturales ((8)-(7)-(5» y
sustrayendo el aporte de las superficies glaciares que no son drenadas por las represas Mamankota y Livinosa
((4d)+(4b», se obtendria una estimaci6n dei escurrillliento de las precipitaciones que toma en consideraci6n el

1 desbordamiento de los canales de derivaci6n allago Zongo.

1
1 Dren~ie de las
Agua originaria dc
La nievc

1 .....
...... .
....................................... ~
:
, 2)
.

1 scurrimiento al exterior de la cuenca

Sub-cuenca de
Sub-cuenca dei Sub-cuenca de Partes dc la cuenca verticn e

1
lagos Tiquimani. que no son drcnadas cn la
valle Coscapa lago .Guaraguarani.
al exterior dei represa Zongo
(/ago Sankayuni). al exterior dei
valle dei 1
al exterior de la vallée du Zongo

••
Zongo 1
cuenca Sainani 1 1
1 4a)[Link]
1
1
1 '
.' ..··
i 3')
Ô..·
,...
1 l
1·'
1"

4b)".. J.....
"31»
;..
.

rC)..:.J: •
3c) 1
1 4d)
[Link]j:::>

1 ~ t "\f"
l':stockage en las
~ 1 ••
••
[Link].
1
1
1
1 .o,,(j,;oo",
represas Maman-
1 •
•• 1
1
1 •• 1
1 \
Canales Cos-
Rébalse cuando la
1
1

••
1
1
1
capa, Tribu-
Represa esta lIena •
1 tario 1 et Il
1
•• 1
1
...........,6)
..........
1
1
1
•• 1

1 1
1 • 1

~
1
~I

1 Afluentes dcl rio Zongo, Canalcs deaducci6n de las plantas

............................................ •_ 1 "

1 l-----T~·~~~"i~y~~~'l----~--j
Planta Santa
Rosa Il

1 Rio Zongo. aguas arriba planta Sainani

1
1
1 1RIO Zongo, aguas a6ajo de Chururaqui Il

:::·
t··
· ·::•• 1) .. Fluio dCSCOllllCido
1 •• •• ~
RIO Zongo. aguas a6aJo planta Samani

C --....1)
~ - 1"1'
' uuo conoci'd0

1 Figura 13 Esquema de los flujos de agua en el valle dei Zongo en estacion de lluvias. Las
explicaciones se encuentran en el texto

1
1
4. Mejoramiento y aplicaci6n deI método para la reconstituci6n de
1 caudales naturales en estaci6n seca

1 4.1. Reconstituci6n de caudales turbinados

1 4.1.1. Estimaci6n de los rendimientos de las plantas hidroelé.ctricas

Para el estudio preliminar de la inf1uencia deI cambio de stock sobre los caudales medidos en
Sainani (cf. i), nos hemos basado en rendimientos aproximados de las plantas hidroeléctricas,

1 es decir en los factores de planta hidroeléctrica utilizados pOl' la COBEE (cf. 2.2.4). Dichos
factores se basan en rendimientos de turbinas que se conocen desde la instalaci6n de las
plantas hidroeléctricas. Ahora bien, dichas turbinas han sido reemplazadas pOl' modelos

1 nuevos, que tienen mejores rendimientos. Para el estudio preliminar, los factores de plantas
hidroeléctricas utilizados eran accptables, pero para calculaI' los caudales naturales a aplicar
con los métodos de la reconstituci6n de los caudales en estaci6n seca y en estaci6n de lluvias
1 (cf. 3.2.1, 3.2.2), se necesitan valores mas precisos. Como los rendimientos, y pOl'
consecuencia los factores de plantas hidroeléctricas, varian con la potencia instantanea (en
relaci6n con la potencia maxima 0 capacidad), conviene hacer variaI' estos factores de planta
1 hidroeléctrica. Se estima que son de 25, 50, 75 Y 100% de la pOlencia instantanea en relaci6n a
su capacidad. Para obtener un punto cada 5% de la curva, se interpolan linealmente.

1 Los factores de planta hidroeléctrica fueron estimados diferenciando los rendimientos totales
de las plantas hidroeléctricas. Se han considerado las pérdidas en el conducto forzado y en el
inyector de la turbina (rendimiento deI conducto), las pérdidas en la meda hidniulica
1 (rendimiento de la turbina) y las pérdidas deI generador de elcctricidad (rendimiento deI
generador) (cf. Cuadro 2). Generalmente, los rendimientos de los conductos varian mas de una

1 planta hidroeléctrica a otra. Las pérdidas en el conducto dependen de diferentes factores,

siendo los principales el largo dei conducto y su mgosidad. Se dispone de valores calculados
de los rendimientos deI conducto forzado de la planta hidroeléctrica Santa Rosa Il. El método

1 de calculo se basa en mediciones de presi6n en el inyector y potencias instantaneas (cf. Anexo

C). Este permite obtener separadamente las pérdidas (los rendimientos) en el conducto, en el
inyector y en el conjunto turbina/generador. Respecto a las otras plantas hidroeléctricas, no se
1 dispone de mediciones necesarias para la aplicaci6n deI método. La caida deI conducto
forzado de la planta hidroeléctrica Sainani es la mas corta (291 m), 10 que lleva a rendimientos
pr6ximos a 1. Se han elegido valores experimentales de la literatura para los rendimientos deI
1 conducto de esta planta hidroeléctrica. Los rendimientos de los otros conductos fueron
calculados pOl' interpolaci6n lineal. Las pérdidas en el inyector, que se incluyen en el
rendimiento deI conducto, pueden sel' estimados para la planta hidroeléctrica Santa Rosa II
1 segun el mismo calculo. Se hace la hip6tesis que no varian mucho de una planta hidroeléctrica
a la otra.

1 El rendimiento de la turbina depende principalmente de la edad de la meda hidraulica. Como

las medas de todas las plantas hidroeléctricas datan de los ultimos afios, se hace la hip6tesis

1 que los rendimientos de las turbinas corresponden a los rendimientos indicados pOl' los
fabricantes, sin consideraI' las otras pérdidas en las turbinas, que son pequefias en relaci6n a
las de la meda. Estas son las informaciones que tenemos para las plantas hidroeléctricas

1 Tiquimani y Sainani. Se ha adoptado para la planta hidroeléctrica Zongo, mas vieja que la
planta hidroeléctrica de Tiquimani que es nueva, el mismo rendimiento de las medas
hidraulicas que en Sainani. Para Santa Rosa II, no es necesario estimaI' los rendimientos, ya
1
1
que el calculo mcncionado antes da dircctamente los rendimientos dei conjunto
1 turbina/generador.

Los rendimientos dc los gcneradorcs varian pnicticamente de una planta hidroeléctrica a la

1 olra y no disminuycn con cl licmpo. POl' csa raZ<ln nos basul110s cn los valorcs indicados por
los l~lbrical1lcs, suponicndo lJUC los valorcs son los mismos cn las planlas hidroeléctricas de
las cuales no se tienc informaci6n.
1
Tiquimani 25 50 75 100 Origen:
1 Turbina
Générador
0,893
0,932
0,902
0,961
0,906
0,969
0,899
0,972
datos dei fabricante
datos de Santa Rosa II

1 Conducto
Total
0,921
0,767
0,875
0,759
0,841
0,738
0,762
0,666
interpol. lin. caidas Sain./SR II
caida bruta de agua (m) : 487
F. P. (kWh/m3) 1,017 1,008 0,980 0,883 capacidad max. (MW) : 9.4

1 Santa Rosa II
Turbina
25
0,730
50
0,832
75
0,846
100
0,843 calculo (Turbina y Gén./Gén.)
Générador 0,932 0,961 0,969 0,972 datos deI fabricante
1 Conducto
Turbina y Gén.
0,921
0,680
0,864
0,800
0,818
0,820
0,731
0,820
calculo °1)
calculo °1)
Total 0,626 0,691 0,671 0,599 caida bruta de agua(m) : 830
1 F. P. (kWh/m3)
Sainani
1,416
25
1,563
50
1,518
75
1,356
100
capacidad max. (MW) : 6.8

1 Turbina
Générador
0,884
0,925
0,895
0,956
0,901
0,967
0,897
0,9715
0,779
datos dei fabricante
datos dei fabricante
expérience °2)
Conducto 0,921 0,882 0,854

1 Total
F. P. (kWh/m3)
0,753
0,597
0,755
0,598
0,744
0,590
0,679
0,538
caida bruta de agua(m) : 291
capacidad max. (MW) : 10.4
Zon2;o 25 50 75 100
1 Turbina
Générador
0,884
0,932
0,895
0,961
0,901
0,969
0,897
0,972
como datos de Sainani
coma datos de S.R. II
Conducto 0,921 0,878 0,846 0,769 Interpol. lin. chutes Sain./SR
1 Total 0,759 0,756 0,739 0,671
II
caida bruta de agua(m) : 377
F. P. (kWh/m3) 0,779 0,776 0,759 0,689 capacidad max.(MW):4.4/l 0.5
1 Cuadro 5 Rendimientos y lactores de planta hidroeléctrica en lunci6n dei porcentaje de

1 la potencia instantanea en relaci6n con la potencia nominal (capacidad). 0') Calculos dei
rendimiento a partir de series de mediciones de la presi6n en la salida dei conduclo forzado
y la palencia inslanlanea correspondiente (paquete de la COBEE). °2) Valores aproximados

1 dados par el responsable de la producci6n de electricidad.

El resumen de calculos de los rendimientos y de los factores de planta hidroeléctrica se

1 presenta en el Cuadro 5. Se ve que es principal mente la altura de la calda que determina la
variabilidad de los factores de planta hidroeléctrica. Los rendimientos totales son
practicamente los mismos para todas las plantas hidroeléctricas, a excepci6n de Santa Rosa II,
1 cuyo rendimiento total es debilitado por la aHura excepcional de la calda. Como los
rendimientos varian poco entre elIos, la correcci6n repercute en los caudales de manera

1
1
1 sensiblemente proporcional; es decir que todos los caudales aumentan 0 disminuyen deI
mismo factor. Para la reconstituci6n de los caudales naturales se trabaja con diferencias de

1 caudales turbinados y de stocks de agua en las represas (3.2.1 Y 1.1.1). POl' ello es que una
correcci6n de los rendimientos, y pOl' consecuencia, de los caudales turbinados, puede influir
mucho en el hidrograma de los caudales reconstituidos. Incluso una correcci6n proporcionai

1 de los caudales turbinados, no conduce a un aumento 0 a una disminuci6n, de un factor

constante, de los caudales reconstituidos. Para poder reconstituir correctamente los caudales
naturales de la cuenca, es absolutamente necesario que la estimaci6n de los rendimientos de

1 las plantas hidroeléctricas sea correcta. Habria que mejorar estas estimaciones con datos
complementarios como aforos 0 la aplicaci6n deI método de estimaci6n de la COBEE a otras
plantas hidroeléctricas.
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
4.1.2. Reconstitucion de caudales turbinados

1 Los caudales turbinados de las plantas hidroeléctricas Zongo, Tiquimani, Santa Rosa II y
Sainani, que interviencn cn la detcrminaci6n de los caudales naturalcs, fucron reconstituidos a
partir de potencias instantâncas, regisLradas a inLcrvaJos de un hora 0 de una media hora

1 (durante las horas punta). Se han determinado los factores de .plantas hidroeléctricas y los
caudales turbinados en cada intervalo de medici6n. Los caudales diarios son obtenidos por
suma de esos caudales entre una hora elegida y la misma hora al dia siguiente. Como se ha
1 visto en el capitulo anterior, es necesario desplazar los caudales turbinados por el tiempo de
propagaci6n entre las plantas hidroeléctricas para poder reconstituir los caudales naturales (cf.
3.1.3). En la figura 14 los caudales diarios de las cuatro plantas hidroeléctricas se presentan
1 por el periodo 2 de abril al 30 de septiembre 1997.

Se ve claramente que hasta fines de abril los caudales en Sainani son constantes, es decir al
1 limite de la capacidad maxima de la planta hidroeléctrica. Este limite ya no es superado sino
hasta fines de septiembre. Este es el periodo durante el cual los caudales turbinados
corresponden a los caudales totales en Sainani y con el cual se va a trabajar mas tarde (cf.
1 2.2.4).

1 50‫סס‬oo , Déborderœlll Tout le débit naturel passe dans l'usine Sainani

.... ! ~
Débits turbinés 1997

45‫סס‬oo --Zongo

1
r·'·:,'········ ", \

1·' " --Tiquimani :'. ,

400000;',' i ~'." ~ •
i
- - Santa Rosa "
350000 ...' ",
:~:"

1
"

" /.
.... , .. Sainuni
\ '
300000 . , .'
g .'
~:
~ 25‫סס‬oo
e
....... ' ; ••• . : ..... " ....... ',:':.",': l' • .\. :.', ,

.
._..';.. . i
1 20‫סס‬oo

150000 i
'0' \: ••••• '".: ." •••••: ... '. ''.: ' ..... \:'
.,,,:
.....: ...:'.
",:
1

100000 r - - - - \ - - - .... .,,''''.....\

1 50000 :
',,-/:;-.lJ.J.~-~~~~~~~~nr~~
o'
1 Figura 14
2-avr 17-avr 2-mai 17-mai I·juin

Caudales turbinados reconstituidos para las plantas hidroeléctricas Zongo,

16-juin I~iuil 16-juil 31-juil 15-août 30·août 14-sept 29-sept

1 Tiquimani, Santa Rosa II y Sainani, que intervienen en la determinacion de los caudales

naturales
A partir de los caudales horarios se pueden calcular los caudales medios para cada intervalo

1 de tiempo, es decir para los caudales medios diarios a partir de horas fijas. Los hidrogramas
diarios medios de las plantas hidroeléctricas obtenidos asi son presentados y discutidos en el
parrafo siguiente.
1 4.1.3. Discusion de los hidrogramas diarios medios de las plantas hidroeléctricas
Se puede observar en la Figura 15 que los caudales turbinados de las plantas hidroeléctricas
1 Zongo y Tiquimani varian sensiblemente durante el dia. Estas plantas hidroeléctricas se
encuentran al frente deI sistema, aguas abajo de las represas importantes (lago Zongo y lagos
Tiquimani) que drenan de las subcuencas relativamente pequefias. Por 10 tanto no son plantas
1 hidroeléctricas de corriente de agua y sus producciones plantas hidroeléctricas en estaci6n
seca estan completamente adaptadas a la demanda durante el dia (cf. 0). En consecuenca, los

1
1
caudales turbinados de estas dos plantas hidroeléctricas son los mas elevados durante las horas
1 punta de la demanua, a mediodia y por la tarde. Ese no es el casa de Santa Rosa II. Aunque
aguas abajo deI sistema de represas donde la mas importante es el lago Sankayuni, la
posibilidad de regular la producci6n es limitada. Ya que la subcuenca deI rio Coscapa,
1 particularmente drenada por la planta hidroeléctrica, es mucho mas grande ante los volumenes
de regulaci6n que en el casa de las plantas hidroeléctricas Zongo y Tiquimani. Por 10 tanto, el
hidrograma es menos variable durante el dia. Se observa, a pesar de todo, la influencia de la
1 apertura de las vâlvulas deI lago Sankayuni entre las 16h00 y 19h30, 10 que aumenta la
producci6n y los caudales turbinados de la planta hidroeléctrica durante las horas de punta.
Los hidrogramas medios de las otras tres plantas hidroeléctricas se encuentran finalmente en
1 los caudales turbinados de la planta hidroeléctrica Sainani. Por esta raz6n el hidrograma es
también variable coma en las planta hidroeléctricas Zongo y Tiquimani, con la diferencia que
los caudales turbinados no se anulan jamas, a causa de los escurrimientos que no dependen de
1 maniobras de las represas.

1 La organizaci6n deI sistema hidraulico explica claramente los hidrogramas diarios medios
encontrados en las cuatro plantas hidroeléctricas. No obstante, se observa una anomalia que, a
primera vista, parece dificil explicar. i,Por qué las puntas deI hidrograma de la planta

1 hidroeléctrica Sainani no son desfasadas por el tiempo de propagaci6n desde Zongo y

Tiquimani? S610 hay una explicaci6n posible. Una parte deI agua turbinada por las plantas
hidroeléctricas Zongo y Tiquimani es retenida en cuencas de regulaci6n intermedia y soltada
1 al dia siguiente, a fin de aumentar la producci6n de Sainani durante las horas punta. Ahora
bien, no se ha considerado la variaci6n de stock de estas cuencas intermedias durante el
calculo deI cambio de stock total de las represas (cf. 2.2.3), estimando sus volumenes utiles
1 muy pequefios para tener una influencia significativa. El resultado del estudio preliminar
parecia confirmar esta hip6tesis, porque la influencia sobre los caudales turbinados en Sainani
de la represa mas pequefia estudiada, San Pedro, deI tamafio de las cuencas intermedias, era
1 insignificante (cf. 2.3.1). Pero el agua deI lago San Pedro es soltada de una manera regular.
Durante un periodo de déficit, se abren las valvulas durante ciertas horas deI dia (cf. 2.2.2), sin
soltar toda el agua de una sola vez.
1 Los stocks de agua de las cuencas intermedias, en cambio, deben variar mucho mas, porque se
intenta a través de su gesti6n adaptar mejor la producci6n de las plantas hidroeléctricas aguas
abajo a la demanda de electricidad durante el dia, es decir desplazar por algunas horas los
1 caudales que llegan de las plantas hidroeléctricas aguas arriba. i,Esto hace necesario tomar en
cuenta los cambios del stock de las represas intermedias en el método para la reconstituci6n

1 de los caudales naturales en estaci6n seca (cf. 3.2.1)? Mas adelante se buscara una respuesta a
esta pregunta.
8000 !

1 <E
6000
..... 4000 !
Usine Zongo .- ........

E 1

1 2000
o
!
1

1
---i---:---";";'"';";""..;;.".;..;'"..;.,".;..;""~.
;...1

2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
heure de lajoumée
1
1
1
1 10000 t - .... - .......
Usine Tiquimani

1 8000·.
~ 6000
E
~
... .
'. ''',_
. .•..•. .

1
4000 :

2 ~ : _ -_ _- - - - -

2 3 4 5 6 7 8
'~:_:r._::~-'.', '~'.'.'.'.'.'.'.'.'.'" .-~ .. ~:>;
9 10 Il 12 13 14 15 16 17 18
f'

. ." '.: ~~~: '.

19 20
..

21 22 23
..

24
heure de lajournée

1 5000
4()()()
Usine Santa Rosa"
---- .
....................... . .
o:E 3 000 ".- ............... _

1 E 2000
1000 ï· .... ·· .... ·····....·......·........ ···.. ··· ....... -'" 1
i
~ l"~ ,~ ,~ ;';3~-;~
1 ~----_.
0: ,

__
3

_-
4 ,

. ----------._--
6 7 ,

'-~
9

.--..__
JO Il 12 13 14
heure de la journée
_.... .
" 16 1

_--_.- ._-_ .. _ - - - - -
,', 1

i
20000

1 Us ine Sainan i
15000 r············· ··························· ._ .. '

':: t....... :=....... .~ .

1

1 ~
1

o,
1 1 2 3 4 5 6 7 8 9 10 Il 12 13
heure de la journée
14 15 16 17 18 19 20 21 22 23 24

Figura 15 Hidrogramas diarios medios de las plantas hidroeléctricas (+/- diferencia-

1 tipo), determinados a partir de los caudales horarios entre el 2 de maya y el 30 de septiembre
de 1997

1 4.2. El papel de las pequeiias cuencas intermedias de regulacion

1 4.2.1. Datos e inform~lcion disponibles sobre las cuencas interrnedias y su gestion

Hay tres cuencas intermedias, dos grandes y una mas pequefia. Las dos grandes se encuentran
justo aguas abajo de la planta hidroeléctrica Tiquimani (Dique Botijlaca) y en frente de la

1 toma deI canal Cuticucho, mas abajo de la planta hidroeléctrica Botijlaca (Dique Cuticucho).
La tercera (Dique intermedio) esta situada a tres cuartos de distancia entre las plantas
hidroeléctricas Zongo y Tiquimani. Para la localizaci6n de estas cuencas se puede consultar el
1 esquema de la cuenca de drenaje (cf. 1.2.2. Figura 1). Los niveles de las cuencas intermedias
Botijlaca y Cuticucho son registrados por la COBEE a intervalos de media hora a una hora, asi
como para las potencias instantaneas (cf. 4.1.2). Se conocen los volumenes totales de estas
1 dos cuencas, pero las relaciones altura de agua/volumen son desconocidas. Sobre la tercera
cuenca intermedia (Dique Intermedio) no se conoce cl volumen y los niveles no son
apuntados. A causa de la ausencia de datos, se ha preferido estimar la variaci6n dei stock de
1 agua en las cuencas en el periodo de la estaci6n seca de 1997, con ayuda de los hidrogramas
horarios de las cuatro plantas hidroeléctricas Zongo, Tiquimani, Santa Rosa II y Sainani,
cuyos datos de base (potencias instantaneas) son disponibles. La variaci6n de stock y su
1 influencia sobre los caudales naturales en Sainani es enseguida ca1culada e interpretada.

Los volumenes de las represas Botijlaca y Cuticucho son respectivamente de 16 880 m3 y de 9

1 390 m 3• Los volumenes de la represa Intermedia y de la represa de carga de la planta

1
1
hidroeléctrica Sainani deben ser mucho nuis bajos. Se dispone de informaci6n cualitativa
1 sobre la gesti6n de las represas intermedias, brindada por el responsable de la gesti6n de
producci6n: se intenta maximizar la producci6n durante las horas punta; es decir, soltar el
maximo de agua de las represas de manera que este llegue a las horas punta a las plantas
1 hidroeléctricas con la corriente de agua de gran capacidad, siendo Sainani la primera aguas
abajo. Esta tipo de maniobra perturbe la producci6n de electricidad en las plantas
hidroeléctricas mas pequefias como Botijlaca, situada entre las represas de regulacion de
1 Botijlaca y Cuticucho, y Cuticucho, aguas abajo de la represa Cuticucho. Asi, cuando la
represa Botijlaca esta llena, el agua no solamente es redirigida por el canal Botijlaca, que
alimenta la planta hidroeléctrica, sino que también es principalmente desaguada en el rio
1 Zongo. El rio es utilizado coma volumen de almacenamiento, ya que los caudales toman mas
tiempo para llegar aguas abajo de la planta hidroeléctrica Botijlaca cuando se propaguan en el

1 lecho natura!. El mismo principio se aplica para el casa de la represa de regulacion y de la

planta hidroeléctrica Cuticucho. Otro principio de la gesti6n de las represas intermedias
consiste en tratar de alcanzar los niveles de llenados maximos hacia las seis de la mafiana. De

1 esta manera, estan listas para ser vaciadas a fin de aumentar la producci6n eléctrica durante las
horas de punta. En estacion seca, siempre segun las informaciones dei responsable de la
gestion de producclon, casi todos los caudales provenientes de las plantas hidroeléctricas

1 Zongo y Tiquimani pueden ser desplazados hacia las horas punta de la producci6n de las
grandes planta hidroeléctricas aguas abajo, utilizando los vohimenes de stock disponible. Es
decir que no hay caudal turbinado en Zongo y Tiquimani que no pase directamente a Sainani,
1 sin ser desplazados en el tiempo por el almacenamiento en las represas intermedias.

4.2.2. Un método para la estimacion de la curva de almacenamiento dei conjunto de las

1 cuencas intermedias

Se pueden hacer cinco hip6tesis simplificadoras para reconstituir los vohimenes almacenados
1 en las represas intermedias:

1) Los caudales en el lecho natural dei rio Zongo se propagan aproximadamente como una
1 onda cinemâtica.
2) La velocidad de propagaci6n de crecida es mas 0 menos constante en todo el trama entre la

1 planta hidroeléctrica Zongo y la planta hidroeléctrica Sainani.

3) En estaci6n seca, todos los caudales que provienen de las plantas hidroeléctricas Zongo y
Tiquimani son desplazados en el tiempo por el almacenamiento en las represas.

1 4) Los caudales almacenados en las represas son mas bien propagados por ellecho natural que
por los sistemas artificiales canal Botijlaca / planta hidroeléctrica y canal Cuticucho /
planta hidroeléctrica Cuticucho (cf. 4.2.1).

1 5) Los caudales que provienen dei desalmacenamiento de la represa Livinosa estan delante de
los caudales turbinados en las plantas hidroeléctricas Tiquimani y Zongo.

1 La primera hipotesis se justifica porque no se necesita de una gran precision para el analisis de
la evoluci6n dei stock total estimado. La segunda se basa en el hecho que los declives no
varian mucha en cl canal considerado. La tercera y la cuarta hip6tesis se basan en la
1 experiencia dei responsable de la gestion de produccion (cf. 4.2.1). La quinta hipotesis es
necesaria porque los caudales que provienen dei desalmacenamiento de la represa Livinosa

1 son almacenados en la represa Cuticucho. En el periodo considerado ha sido verificada porque

las maniobras de esta represa apenas comenzaron a fines de septiembre de 1997.

1
1
Las hip6tesis 2) y 3) permiten asimilar las tres represas de regulaci6n y la represa de carga de
1 la planta hidroeléctrica Sainani a una sola represa de retenci6n. El volumen lltil de esta represa
ficticia contiene los volumenes de las tres represas de regulaci6n y de la represa de carga de la
planta hidroeléctrica Sainani, asi como la capacidad de stock de la secci6n deI rio Zongo entre
1 la primera cuenca (Dique Intermedio) y la represa de carga de la planta hidroeléctrica Sainani
(cf. 0).

1 Bajo las hip6tesis 1) Y 2) la posici6n de esta represa ficticia, en la secci6n deI rio entre la
primera y la ultima represa, no tiene influencia sobre la variaci6n deI stock de agua total
durante el dia (curva de almacenamiento). Como, bajo las dos hip6tesis, una cantidad de agua
1 se propaga independientemente de su volumen 13 y con una velocidad constante en toda la
secci6n hacia aguas abajo, entonces cualquiera sea la posici6n de la represa, los desfases en el
tiempo de los hidrogramas de los caudales de entrada y de salida son los mismos y la curva de
1 variaci6n de stock no cambia de forma. Lo que cambia con la localizaci6n de la represa
ficticia es la posici6n absoluta en el eje de los tiempos de esta curva. Es decir que si la represa
esta aguas arriba, se alcanza un punto singular de la curva (por ejemplo el stock maximo) mas
1 que si se encontrara aguas abajo. El desfase en la escala de tiempo corresponde al tiempo de
propagaci6n deI agua en el rio entre dos posiciones de la represa ficticia. A continuaci6n (yen
los graficos) vamos a suponer que la represa ficticia se encuentra en el lugar de la cuenca de
1 carga de la planta hidroeléctrica Sainani.

1 A fin de poder estimar la variaci6n deI stock de agua durante el dia, hay que conocer los
hidrogramas instantaneos de entrada y de salida de la represa ficticia. El hidrograma de
entrada es facil de determinar. Basta con retardar los hidrogramas instantaneos de la plantas

1 hidroeléctricas Zongo y Tiquimani, de la duraci6n de transferencia entres estas plantas

hidroeléctricas y la planta hidroeléctrica Sainani (Figura 16). Las dos plantas hidroeléctricas
producen unicamente durante las horas punta y los caudales totales constituyen el volumen
1 que hay que desplazar hacia las horas punta de las grandes plantas hidroeléctricas aguas abajo,
a través de la gesti6n de las cuencas intermedias (cf. 4.1.3 y 4.2.1). Las entradas de caudales
naturales en la represa son tomadas en cuenta en la estimaci6n de los caudales de salidas, ya
1 que la retenci6n de un caudal natural instantaneo disminuye mucho el hidrograma de salida,
como se vera a continuaci6n.

1 El hidrograma de salida puede estimarse utilizando los hidrogramas instantâneos de las cuatro
plantas hidroeléctricas Zongo, Tiquimani, Santa Rosa II y Sainani. El hidrograma instantâneo
de los caudales turbinados en Sainani se compone deI hidrograma instantaneo de salida de la
1 represa ficticia, deI hidrograma instantaneo de Santa Rosa Il que constituye un aporte lateral
(cf. 3.1.4), y de los caudales naturales de la cuenca de drenaje hasta la planta hidroeléctrica
Sainani, sin los aportes de las subcuencas captadas por las plantas hidroeléctricas Tiquimani y
1 Zongo. En efecto,bajo la hip6tesis 3), la totalidad de los caudales de estas cuencas (los
caudales turbinados) es retenida en las represas intermedias y por 10 tanto contenidas en el

1 hidrograma de salida. El hidrograma de salida de las represas corresponde final mente a la

diferencia entre el hidrograma de los caudales turbinados en Sainani, sin los aportes de Santa
Rosa 11 14 , y los caudales naturales hasta Sainani, sin los aportes de las cuencas captadas por
1 las plantas hidroeléctricas Tiquimani y Zongo (Figura 16).

1 13
14
Esta condici6n es verificada por la ausencia de término de difusion para la onda cinematica.
No es necesario un desplazamiento hacia atras de los caudales de Santa Rosa II porque el tiempo de
propagacion entre las dos hidroeléctricas solo es de diez minutos.

1
1
1
1 :Débits instantanés
i: Debits de ('usine Santa Rosil'I!, moycn- dc Santa Rosa Il
~ Ë ~ne mobile sur:24 h~ de.çalé~ 5 nin,
,~
0lJ ,

1 ~ ...
~' ...
o
5~ E
6'.~ E
E"O .~.-
oU
Débits de l'liSme--Tlquunam moyenne§
: mobile sw'·[Link];'; ée ée 2h 50 min. .~
(/)
1':t'C
;ë '~§ e ~ ~
1 l-

A 'éli'iiiilie l,uSine Zongo, moyerme 1

,5 ~. ~. ~ oliile sûÎ:Î24 heures, décalée 7h 30 min. ~

>. ...

1 Débits instantanés
usine Zongo
Débits instantanés
usine Tiquimani
j
~ ~
0lJ
c~
.~~
<l>
··············································i-(13
'"
·U
~
~ '; §
1
Débits naturels du bassin versant

1
"0 '" jusqu'à Sainani, sans les apports ;;;
1< 'fenêtre' de 24 heures
:Ê ~ des sous-bassins captés par les .;
oE usines Tiquimani et Zongo ~
o

1 :< 'fcnêtrc' uc 24 hcurcs >

neeonslitution de l'hydrogramme d'entrée Reconslilulion de l'hydrogramme de sortie
1

1
" J

Figura 16 Reconstitucion de los hidrogramas de entradas y de salida de la represa jicticia

(,C6mo se pueden estimar los caudales naturales hasta Sainani, sin los aportes de las
1 subcuencas captadas por las plantas hidroeléctricas Zongo y Tiquimani? Se parte de la
suposici6n que el hidrograma diario medio de estos caudales naturales es constante. Es decir
que los caudales instantaneos varian en funci6n dei tiempo, pero no dependen de la hora dei
1 dia. Esta claro que el deshielo de los glaciares no es constante durante el dia, 10 que deberia
repercutir en la forma dei hidrograma diario medio. No obstante, si se mira el hidrograma

1 diario medio de la planta hidroeléctrica Santa Rosa II (Figura 15), se constata que los caudales
instantaneos medios son mas 0 menos constantes durante el dia, con una pequefia ondulacion
hacia la tarde, debido al desalmacenamiento deI lago Sankayuni (cf. 4.1.3). Ahora bien, la

1 planta hidroeléctrica Santa Rosa II drena la parte superior deI valle Coscapa, que contiene
superficies glaciares importantes (cf. 1.2.2, Figura 1). Se puede objetar que es justamente la
gesti6n de la represa Sankayuni que amortigua el hidrograma diarios medio de Santa Rosa II.

1 Pero incluso si fuera el caso, nos damos cuenta que el pequefio volumen de agua contenido en
la ondulacion dei desalmacenamiento (Figura 15), repartido en el dia, no puede cambiar
mucho la forma general dei hidrograma.
1 El hidrograma instantaneo de los caudales naturales de la cuenca de drenaje hasta Sainani, sin
aportes de subcuencas captadas pOl' las plantas hidroeléctricas Zongo y Tiquimani, puede ser
1 estimado sustrayendo, de la media m6vil de 24 horas de los caudales turbinados en Sainani,
las medias m6viles de los caudales turbinados en las plantas hidroeléctricas Santa Rosa II,
Zongo y Tiquimani, desplazados por su tiempo de propagaci6n hasta Sainani. Evidentemente,
1 los caudales naturales estimados asi no se adaptan a la realidad cuando hay cambios rapidos
de caudales naturales instantaneos reales. En estaci6n seca, en ausencia de fuertes

1 precipitaciones, eso no sucede. La diferencia entre los caudales naturales instantaneos

estimados asi y los caudales instantaneos turbinados en Sainani, de los que sustrayeron los
caudales instantaneos de Santa Rosa II, corresponde al hidrograma buscado de los caudales de

1 salida de la represa ficticia (Figura 16). En la Figura 17, se presentan los hidrogramas diarios
medios de los caudales naturales y de los caudales de Sainani, excepto los de Santa Rosa Il,

1
1
calculados en un periodo de observaciém de la estaci6n seca 1997, deI 1 de maya al 27 de

1 septiembre.

1 18000 i
1

--Débits turbinés instantanés à l'usine Sainani, moins ies

1 16000 ~

14000 .
-
débits mstantanés de Santa Rosa JI

- Moyennes mobiles sur 24 heures des débits turbinés à

l'usme Saman l, sans les moyennes mobiles des déblls
de Santa Rosa Il. T1quimalll et Zongo

1
"
12000 ' ," "
"

oE 10000 •', ••• , ••• , '" . . . _.. '" ...... • .:Hydrograrrme

M .....

1 E
8000 ~

6000
1

!i -
- ' " .. .. • • • .. ....
_. . . . . . . . . . . . . . . . . . . . . . . . .

_=-=::-::_=--__0::""::_ -
.

-::.=-- - - - -
. . .
~.
, .
_. . . . . .

.......
;' moyenne de sortie
o. . .•. :.: ••••• _••.•.•••.••.•••••••

• .a.:

1 4000 ~.

2000 :
~ :: ::::::::~:::::::; .. - - _ .

1 o i
234 6 7 8 9 10 Il 12 13 14 15 16 17 18 19 20 21 22 23 24
Heure de la joumée

1 Figura 17 Hidrograma diario medio de los caudales instantcineos turbinados en la planta

hidroeléctrica Sainani sin los caudales instantaneos de la planta hidroeléctrica Santa Rosa Il.
0- ,

1 Yel hidrograma medio de los caudales naturales de la cuenca hasta Sainani, sin los aportes
de las subcuencas captados pOl' las plantas hidroeléctricas Zongo y Tiquimani. Media de los
valores horario dei periodo de observaci6n, dei JO de maya al 27 de septiembre de 1997. La
1 diferencia de estos dos hidrogramas es el hidrograma diario medio de salida de la represa
ficticia. La escala de la hora dei dia corresponde a la posici6n de la represa ficticia en el
lugar de la cuenca de carga de la planta hidroeléctrica Sainani
1
1 En los caudales de salida calculados de esta manera se incluyen las entradas naturales de la
represa ficticia. La retenci6n de un caudal natural instantaneo en la represa ficticia, a una hora
dada deI dia, disminuye el caudal turbinado instantaneo en Sainani un tiempo después. Se
1 puede suponer que los caudales naturales en Sainani, que son de medias de 24 horas, casi no
disminuyen a causa de la retenci6n. En consecuencia, el caudal instantaneo retenido es
incluido en la diferencia entre el caudal natural instantâneo y el caudal turbinado instantaneo,
1 es decir, en el caudal de salida instantaneo calculado para esta h~ra. Por eHo es que se puede
admitir que las salidas instantaneas toman valores negativos.

1 4.2.3. Simulacion de la curva de almacenamiento

Conociendo los hidrogramas instantâneos de los caudales de entrada y de salida de la represa

1 ficticia, la variaci6n del stock de agua en esta represa puede fâcilmente ser reconstituida. Entre
dos pasos de tiempo, el cambio en la represa es la diferencia entre los caudales instantâneos de
entrada y de salida. El volumen de almacenamiento titil que se necesita para no tener
1 desbordamiento (hip6tesis 3), es calculado pOl' la diferencia entre el maximo y el minimo deI
stock alcanzado durante el periodo de observaci6n. Si se considera todo el periodo de

1
1
1 observacion, deI 1ero de maya al 30 de septiembre de 1997, el volumen util de la represa debe
ser de aproximadamente 60 000 m3. Parece un volumen un poco exagerado. Este volumen de

1 almacenamiento solo es necesario durante dos dias, el 2 y el 16 de septiembre de 1997. En

esos dias los caudales turbinados en la planta hidroeléctrica Zongo eran excepcionalmente
altos a causa de pruebas de funcionamiento de la nueva planta hidroeléctrica Zongo II (cf.

1 1.2.2). Estos caudales adicionales no pudieron ser retenidos en las cuencas de drenaje
interrnedias. Si no se consideran estos dos dias en la modelizaci6n de la curva de
almacenamiento, el volumen util necesario de la represa ficticia es de 42 500 m 3 . Si se

1 consideran los volumenes utiles conocidos de las dos cuencas Botijlaca y Cuticucho (26 000
m 3), el conjunto de los volumenes utiles de la cuenca Intermedia y de la cuenca de carga de
Sainani tendria un volumen de 17 500 m 3 . Este volumen adicional parece razonable y

1 confirma la informacion deI responsable de la gestion, que todos los caudales turbinados en
las plantas hidroeléctricas Zongo y Botijlaca son almacenados y desplazados en el tiempo
hacia horas punta de las grandes plantas hidroeléctricas en estacion seca (cf. 4.2.1). La
1 hipotesis 3) deI método es justificada. La Figura 18 muestra la cvolucion diaria media de la
curva de almacenamiento asi como los hidrogramas diarios medios de las entradas y de las
salidas.
1 Se observa en la Figura 18 que el stock de agua alcanza su maximo en promedio hacia las
6hOO, en el casa que la represa ficticia se encuentre frente a la represa de carga de la planta
1 hidroeléctrica Sainani (Figura 18). Si se supone que la reprega esta situada entre la represas
mas grandes, Botijlaca y Cuticucho, se debe retardar la curva de almacenamiento coma
maximo 2h30, que corresponde al tiempo de propagacion de la cuenca Botijlaca en la planta
1 hidroeléctrica Sainani. El maximo deI stock es alcanzado entre 4h 15 Y 6hOO. Esto coincide
con el principio de gestion de llenar todas las cuencas intermedias antes de las 6hOO (cf.

1 4.2.1).

1
1
1
1
1
1
1
1
1
1 12000 .
,,
.
.- .... - .
,
.
....... Enlrees auxrelel\ues inlerm:dialres

--SortieS des retenues II\lennediuircs

IO()()() , ' .. - - Stockage dans les retenues

1 8000 !' •• ' •

intennedlaires ••••• ~

1 ~E 6000·

_ • • tII
.,
:
,

1 4000 .

., 10000

. ..
,
2000 i ,

1 o. 1 j

0 0 0 0
~ ~ ~ ~ Cl 8 ~
0 0 0
86 8~ 8~ 8~ 8~ 8~ 0 Cl .,.
Cl Cl
or, 0 8 0
0
0
0
Cl
0
Cl
0
Cl ...,
- N M ~
- -- -r-.i '"
6 \ci r:..: 00 b; 0 N

1 1le ure <le la jOlllllèe

N N N N

_ J
Figura 18 Entradas y saUdas medias en las cuencas intermedias durante el dia y
1 variacion media dei stock de agua en las cuencas (+/- diferencia fipo). La escala de la hora
dei dia corrresponde a una posicion de la represaficticia en le lugar de la cuenca de carga
de la planta hidroeléctrica Sainani.
1 Hay que hacer una observaci6n importante respecta a la validez deI amilisis de la curva de
almacenamiento, que se hani en el parrafo siguiente. La puesLa en servicio de la nueva planta
1 hidroeléctrica no cambiara la gesti6n deI lago Zongo ni de las cuencas intermedias en estaci6n
seca. Ellago Zongo sera vaciado de una manera lineal como en el pasado (cf. 2.2.1), ya que la
capacidad adicional de la nueva planta hidroeléctrica no es utilizada sino para aprovechar los
1 escurrimientos abtindantes en estaci6n de lluvias (1.2.2). La respuesta a la pregunta de si la
gesti6n de las cuencas intermedias debe ser tomada en cuenta en la reconstituci6n de los

1 caudales naturales de la cuenca de drenaje, que sera dada a continuaci6n, es igualmente valida
para las futuras estaciones secas.

1 4.2.4. La influencia de la variaci6n de stock en las cuencas intermedias sobre los

caudales naturales
La curva de almacenamiento durante el periodo de observaci6n en estaci6n seca 1997 (2 de
1 maya al 27 de septiembre) permite calcular la variaci6n de stock de un dia al siguiente, para
una hora de observaci6n dada. La Figura 20 da un ejemplo de calcula entre el 27 de agosto a
las 23h30 y el 28 de agosto a la misma hora. Para cada hora deI dia se puede calcular el
1 promedio en toda la estaci6n seca y la variaci6n maxima alcanzada durante el periodo de
observaci6n. Lo que nos interesa es la inf1uencia de estas variaciones de stock sobre los
caudales naturales diarios en Sainani que seran reconstituidos en el pr6ximo capitulo. Para
1 estimar esta inf1uencia, nos servimos de los caudales naturales reconstituidos en un parrafo
ahterior (cf. 4.2.2), 10 que sobrestima un poco la influencia. Las variaciones de stock de las
cuencas intermedias para cada hora deI dia, promediadas por el periodo de observaci6n y
1 expresadas en porcentaje de los caudales naturales en Sainani, estan representadas en la Figura
20.

1
1
1
1 14000 . .
- - Entrees à la retenue
SortieS de la relenu
--Stock d'cau
Changement de stock: 10'585 m3
En % du débit naturels Journalier
(Sainan 1sons SR Il. Zongo, Tlq.): 17%
· 40000

12000 . · 35000

1 :?
....
10000 .
: 30000

.§. · 25000 ~

1 .!l
~
II
8000'

6000·
,
l
: 20000 ~
:>
'"

-""
g
.~ 15000 cil

1 ë
lJJ 4000

2000
10000

............

1 o 1 1 1 li
888 8
1 1 1

~
f ! ~ Il j

8N ë8 $ 888 888 8
1 t ! ! • • 1 :..

~
-T ~

~
• t

8N
; , .. , -,
000
000
• l'----...-......,...J·
88 8
o
- - - -
N ~ ~ 00 o
N N
N.::t~OO~~:!~ 00
-
0
N N
ë ~ N 00 ë
":j

1 Figura 19
Heure de la journée

Calculos de la variaci6n de slock de agua en las cuencas inlermedias enlre el

27 y el 28 de agoslo de 1997 a las 23h30. Esla variaci6n corresponde al 17% dei caudal
1 natural en Sainani, sin los apories de las subcuencas de las plantas hidroeléctricas Zongo y
Tiquimani, escurridos entre esas horas.

1 Hay que recordar que a las 23h30 se deben determinar los caudales turbinados en la planta
hidroeléctrica Sainani, para reconstituir los caudales naturales seglin el concepto desarrollado

1 en el capitulo anterior (cf. 3.2.1, Figura 12). Pero también hay que observar que la escala de
tiempo de la curva de stock es relativa y variable seglin el lugar de la represa ficticia en el
Zongo (cf. 4.2.2). Si se supone primero que la represa ficticia se encuentra en lugar de la
1 cuenca de carga de Sainani, no es necesario desplazar la escala de tiempo. La hora de
comparaci6n de los vollimenes almacenados es 23h30, pues los caudales deI
desalmacenamiento no fueron transferidos para alcanzar Sainani. Si se desplaza la represa
1 ficticia hacia aguas arriba, la escala de tiempo debe ser avanzada por una duraci6n igual al
tiempo de propagaci6n entre la nueva posici6n de la represa y Sainani, pues los caudales
deben ser soltados antes desde la represa (cf. 4.2.2). Pero como la represa se encuentra ahora a
1 cierta distancia mas arriba de Sainani, el cambio de stock debe ser determinado a una hora de
referencia anterior, porque los caudales deI desalmacenamiento toman cierto tiempo para
llegar a Sainani. Evidentemente, estas dos duraciones son las mismas y las dos acciones se
1 anulan. La posici6n de la represa ficticia en el rio Zongo no tiene pues influencia en la
determinaci6n de la variaci6n de stock diario para una hora de reœrencia de medici6n dada de
caudales en Sainani.
1 übservemos la Figura 20, en la cual la represa ficticia esta situada en Sainani, las variaciones

1 de stock en porcentaje de los caudales naturales de Sainani a las 23h30. La intluencia media
sobre estos caudales durante todo el periodo de observaci6n es de 2.7%, la influencia maxima
es de 17%. S610 hay 4 dias (3% de los dias) en los cuales la influencia deI desalmacenamiento
1 / almacenamiento sobre los caudales naturales en Sainani supera el 10%.

1
1
1
hecho podria hacer necesaria la determinaci6n de la vanaClOn de stock de las represas

1 intermedias durante eventos de crecida. Esto puede realizarse ya sea por medio de la
utilizacî6n de medidas de nivel disponibles para las dos represas, Cuticucho y Botijlaca, con
una estimaci6n de las curvas nivel/volumen de las represas, es decîr por la utilizaci6n deI

1 método desarrollado. Podria preverse también una combinaci6n de los métodos.

A continuaci6n se reconstituiran los caudales naturales en Sainani, con ayuda deI principio
1 expuesto anteriormente (cf. 3.2.1, Figura 12). Como todavia no se tiene una modelizacî6n de
los caudales de deshielo deI conjunto de los glaciares de la cuenca de drenaje hasta Sainani,
aun no es posible un analisis fino de los caudales reconstituidos a partir de las precipitaciones.
1 Por e110 es que en esta etapa de los trabajos, no se considerara la influencia de la gesti6n de las
represas intermedias sobre los caudales reconstituidos.

1 4.3. Reconstituci6n de los caudales naturalcs en Sainani en estacion seca

1997
1
4.3.1. Reconstitucion de los caudales naturales en estacion seca 1997
1 Recordemos que, para aplicar un modelo hidrol6gico a los caudales de Sainani, hay que
reconstituir los caudales naturales no influenciados por las operaciones hidraulicas
(almacenamiento y desalmacenamiento de las represas) y sin los aportes externos a la cuenca
1 y de deshielo de los glaciares (cf. 3.1.1). Para este fin se han desarrollado métodos para la
reconstituci6n de los caudales en estaci6n seca y en estaci6n de lluvias (cf. 3.2). Como no se
dispone por el momento de datos hidrométricos, debcmos conformarnos con una
1 reconstituci6n de los caudales en estaci6n seca (cf. 3.2.1) comparando los caudales totales en
Sainani con los caudales turbinados en la planta hidroeléctrica Sainani (cf. 2.2.4). En un
primer tiempo, une se limita a tratar de reconstituir los caudales naturales de la cuenca de
1 drenaje no influido por las operaciones hidraulicas y sin los aportes dei exterior de la cuenca
de drenaje, pero con los aportes deI deshielo de los glaciares.

1 El principio de la reconstituci6n de estos caudales naturales en estaci6n seca fue desarrollado

en un capitula anterior (cf. Figura 12,3.2.1). Recordemos que para reconstituir los caudales

1 naturales en Sainani hay que sustraer de los caudales diarios de la planta hidroeléctrica
Sainani los de las plantas hidroeléctricas Tiquimani y Santa Rosa II, determinados con
desplazamientos en el tiempo que corresponden al tiempo de transito de los caudales entre

1 estas plantas hidroeléctricas y Sainani. A los caudales obtenidos hay que agregar las
variaciones de stock (el desalmacenamiento siendo negativo) de las represas Zongo,
Mamankota, Livinosa y Sankayuni, cuyos volumenes son calculados a la hora de lecturas de
1 nive!. Se eligi6 la hora de la determinaci6n de los caudales en Sainani de manera que el
tiempo de transito de los caudales dei desalmacenamiento de estas represas hasta Sainani sean
tomadas en cuenta 10 mejor posible (cf. 3.1.3, Cuadro 4). En el caso de la represa Mamankota
1 la hora de lectura de nivel es la mafiana (cf. 2.2.2 y Cuadro 3), mientras que se debe conocer la
variaci6n de stock diario a las 14h40 de un dia al dia siguiente. Por ello es que se hace una
estimaci6n dei volumen almacenado a las 14h40 de un dia dado por interpolaci6n lineal de los
1 stocks de agua a las horas de lectura de niveles. Las variaciones de stock diarios de las
represas fueron estimadas de la misma manera que en el estudio preliminar (cf. 2.2.3). Es

1 decir, que se han calculado los volùmenes instantâneos de agua almacenados en las represas
para cada dia a las horas de lectura de nivel con ayuda de las curvas niveles/volumenes. Paras

1
1
las represas mas importantes estas curvas fueron establccidas por la COGEE basandosc en
1 mediones topogrMicas. Para olras represas, se tuvieron que estimar las curvas (cf. 2.2.3 y
Anexo A).

1 En la Figura 21, se presentan los resultados de la reconstituciûn de los caudales naturales en

Sainani, con el aporte deI agua de deshielo de los glaciares, por el periodo de observaci6n dei
1cm de mayo al 30 de septiembre 1997.
1
1 450000 o Dést6ekage des relenues cl apports des
0

bassins c~érieurs 2

1 400000

350000
[] Débits nalurels. y inclus l'apport des glaciers
4

....
1 ~
g
300000

250000
fenêlre Il ju illel - 28 août
6 ~
0
r')
'tE
1 E
200000
8

150000 10

1 100000 IOCY.> ''J i. r.

12
50000

1 o
2-mai 17-mai I-ju in 16-juin I-ju il 16-juil 3 J -ju il 15-aoû t 30-aoû t 14-sept
14
29-sept

1 Figura 21 Precipitaciones medias sobre la cuenca de Sainani, caudales naturales

reconstituidos en Sainani y variaciones de stock diario de las represas dl/l'ante el perfodo de

1 observacion en estacion seca 1997 Los caudales naturales lueron calcu/adus segun el
método de reconstitucion de los caudales naturales en estacion seca (cf 3.2.1). Para las
precipitaciones medias sc ha utilizado el método de isohietas.

1 4.3.2. Discusi6n de los resultados

La suma de los caudales naturales en Sainani y de los caudales dei desalmacenamiento de las

1 represas, con los aportes externos (caudales de Santa Rosa II y una parte de los caudales de
Tiquimani), da el hidrograma de los caudales turbinados en la planta hidroeléctrica Sainani.
Se observa en la Figura 21 que los caudales turbinados y los caudales naturales diarios
1 evolucionan de manera diferente. Dicho de otra manera, la toma en consideraci6n deI
desalmacenamiento de las represas y de los aportes externos de la cuenca cambian el aspecto
dei hidrograma. Observando solamente una parte dei hidrograma (Figura 22) se ve que este
1 cambio se debe sobre todo a la toma en consideraci6n deI desalmacenamiento de las represas.
La evoluc6n de los caudales de la planta hidroeléctrica Santa Rosa II, en cambio, es casi
paralela a la de los caudales turbinados de Sainani, que es un poco menos el caso de los
1 caudales de la planta hidraeléctrica Tiquimani. Esta se explica facilmente. Los caudales
turbinados de Santa Rosa II reflejan principal mente la producci6n de escurrimiento de la parte

1 superior dei valle Coscapa, cuyo comportamicnto debe ser similar al de la cuenca de drenaje
hasta Sainani, de donde resultan las formas paralelas de estos dos hidrogramas. La influencia
de las maniobras deI Ingo Sankayuni sobre el hidrograma de Santa Rosa II es minima, a causa

1
1
1 dei poco volumcn de esta represa en comparaClOn a la superficie drenada (cf. 4.1.3). En
cambio, la cuenca de drenaje externa drenada par los canales de derivaci6n de la planta
hidroeléctrica Tiquimani contiencn los tres lagos Tiquimani, cuya manipulaci6n en estaci6n

1 seca tiene una inf1uencia mas importante sobre los caudales (cf. 2.3.1 y 1.2.1, Figura 1). Coma
la variaci6n de stock de las represas disminuye mas bien cuando las producciones de las
vertientes aumenta, los caudales de Tiquimani evolucionan de manera menas paralela a los

1 caudales de Sainani que los de Santa Rosa II. Par ello es que finalmente los hidrogramas dei
desalmacenameinto de las represas no se parecen al hidrograma de Sainani, influenciando asi
no soIamente los volumenes de las crecidas, sino también sus formas.

1
1 1 •
450000 0
I- 1--· -. ..'-'·-1"1-;,

400000 .. . . ,

.... ... ,
,
1 35000()

300000
. ,
"
"
.' ,
, '
: ..
,
2
' '
... "
. ,

1 0
....
~

'::'
......
250000

E 200000
, ,
Saison sèche 1997
3
....
~
0

'FE
1 150000 .• 4

100000 .
1 50000
5

1 0
1J-juil 15-juil 19-juil 23-juil 27-juil
- 6
31-juil 4-août S-aoûl 12-aoûl 16-août 2û-août 24-août 28-août

i [Link];:;::::J DéSlockage lac Zongo _ Déslockage lac Mamankota

1 1
[Link]J Déstockagc lac Livinos<l
[Link]J Débits Santa R. Il
[Link]J Débits lïquirnani
_ Débits naturels, y inclus l'apport des glaciers
, _ Précipitations ....... Ablation des glaciers (qualitatif)

1 Figura 22 Sumas de los hidrogramas de los caudales naturales, de los caudales

turbinados en Santa Rosa 11 y en Tiquimani y de los caudales dei desalmacenamiento de las

1 represas Zongo, Livinosa, Sakayuni y Mamankota por el periodo comprendido entre el Il de

Julio y el 28 de agosto 1997. Esta suma corresponde al hidrograma de los caudales
turbinados en Sainani. lndicacion de la tendencia de deshielo de los glaciares en punteado
1 (cf 2.2.4).

1 Se observa que, durante los largos periodos sin preclpltaciones, los caudales naturales
reconstituidos varian mucha (Figuras 21 y 22). Es un poco extrafio, pues en ausencia de
precipitaci6n se espera una forma deI hidrograma mas fina. Debe existir otra influencia
1 diferente a las precipitacioncs que hace variar los caudales. Esta influencia viene
evidentemente de los glaciares. POl' el momento, no se tiene ninguna posibilidad de verificar si

1 la reconstituci6n de los caudales naturaJes corresponde a la realidad.

La pregunta es de saber si los caudales que provicnen deI deshielo de los glaciares pueden
explicar la variaci6n dei hidrograma en ausencia de eventos de precipitaci6n. Se observa en

1
1
30% . .Sainani 23:30 h: 16

1 fi
~
25% .
--Moycnne dc l'Inllucncc
- •• - Moycnnc + eeart-typc
....... MaXImum
Moycnnc' 2,7%
Moycnnc +
eeart-ty pc: 5.5%
Maximum: 17%
14 ~
....
E 12 .~
1 --Valcurs au dcssus de 10"10 Valcurs au dcss us
dc 10%: 4
.:. ':
"
......',
, . \0 ~
"li::>
l!!
!l
~

1
8 ]

.... l'!
.. "
" . 6 ,,;
~
g
:s
1 4

2
J:)

z
g
~

1 0"/.
0
0
0
0
0
- :5 O
N ..,
0
' 0
0
v
0
.,.,
0
0
0
-<:i
0
0
,:...: ~
00
:5
a;
0
0
.::;
$ 0
0
N
$
~
0
0
~
0
0
V"l
0
0
-<:i
:5
!::
0
0
00
0
0
a;
:5 ~ ~ ..,
.::;
N
N
~
N N N
0 0

Hcure dc la jouméc

1 Figura 20 Influencia de la variacion de stock de agua en las cuencas intermedias sobre

... -- -

los caudales naturales en Sainani, expresada como la variacion dei stock entre una hora
1 dada y la misma hora dei dia siguiente (hora de observacion sobre el eje x), en proporcion
con los caudales naturales escurridos entre las horas de observacion (caudales naturales de
Sainani, sin aporte de subcuencas de las plantas hidroeléctricas Tiquimani y Zongo, cf
1 Figura 17), En el grafico se muestran las medias de todo el periodo de observacion (2 de
maya al 27 de septiembre), las diferencias tipo, el maxima y el numero de difèrencias
superiores a 10%. La escala de la hora dei dia corre~ponde a una posici6n de la represa
1 jicticia en la ubicacion de la represa de carga de la planta hidroeléctrica Sainani.

1 Para dar un ejemplo deI calculo de la variaci6n de stock a una hora dei dia dada (Figura 20),
se eligiô la parte de la curva de almacenamiento que muestra la variaciôn de stock maxima a
las 23h30. Se ve en la figura que la variaciôn elevada de stock e..ntre el 27 y el 28 de agosto a

1 .23h30 se debe a una diferencia marcada de los picos de caudales de salidas que sirven para
satisfacer las demandas de las horas punta. Si se compara con la evoluciôn media de la curva
de almacenamiento durante el dia (Figura 18), se constata que mas bien el pico deI 27

1 disminuyô en relaciôn al deI 28. Parece que después de algunos dias de escasas
precipitaciones entre el 24 y el 27 de agosto, la crecida causada pOl' este evento fue utilizada
para lIenar las represas intermedias durante una hora deI dia en el cual se encuentran
1 normalmente a nivel minimo. Recordemos que, segun el modo de calculo adoptado, los
caudales de salida incluyen las entradas en las represas intermedias de los caudales naturales
(cf. 4.2.2). Al dia siguiente, el 28 ùe agosto, las aguas almacenadas fueron utilizadas para
1 aumentar la producciôn de electricidad durante las horas punta. Nos encontrabamos ya en
decrecida, pues el stock alcanza su bajo nivel habituaI. El mismo efecto deberia observarse al
comienzo de la crecida, cuando el almacenamiento de los caudales naturales comienza. En
1 cambio, durante la crecida, no hay que esperar fuertes variaciones de stock, pues las represas
deben estaI' lIenadas la mayor parte dei tiempo.

1 Como conclusiôn de este capitulo, se puedc decir que la influencia de la gestiôn de las
represas intermedias sobre los caudales en Sainani es generalmente minima. Pero puede sin

1 embargo alcanzar valores importantes durante los eventos de crecida, sobre todo al comienzo
yal fin, cuando los caudales naturales excedentes son utilizados para lIenar las represas, Este

1
1
1 las Figuras 21 y 22 que los caudales turbinados en Sainani varian tanlbién en ausencia de
precipitaci6n. Si la reconstituci6n de los caudales naturales es buena, el hidrograma de estos

1 caudales debe adaptarse mejor a la evoluci6n dei deshielo de los glaciares que el hidrograma
de los caudales turbinados. En el parrafo siguiente, se va a tratar de verificar si esto es verdad
con ayuda de una estimaci6n cualitativa de la evoluci6n dei deshielo de los glaciares.
1 4.3.3. Estimacion cualitativa dei deshielo de los glaciares
La variaci6n diaria deI stock de agua en el lago Zongo corresponde a la suma de los caudales
1 de entrada y de salida de la represa. Como en estaci6n seca no hay desbordamiento de la
represa, los caudales turbinados en la planta hidroeléctrica Zongo constituyen las unicas
salidas exceptuando las pérdidas por evaporaci6n. Como las variaciones diarias de stock
1 pueden ser estimadas a partir de las lecturas diarias de nivel, së puede calcular el total de los
caudales de entrada, si no se consideran las pérdidas por evaporaci6n. Estas entradas estan
constituidas por los aportes dei canal Tubo que capta los aguas dei emisario dei glaciar Zongo,
1 de los canales Mamankota / Prado que drenan la vertiente norot:ste deI Charquini y de los
canales Alpaca / Vicuiia que son alimentados por la vertiente y dei Huayna Potosi (cf. 1.2,
Figura 1 5.2.1, Figura 25). La mayoria de las superficies drenadas por estos canales es glaciar.
1 Por ello es que en ausencia de precipitaci6n, los caudales de entrada en la represa provienen
principalmente deI deshielo de los glaciares. La Figura 23 muestra estos caudales de entrada

1 durante el periodo de observaci6n en estaci6n seca 1997 con las variaciones de stock de la
represa Zongo (desalmacenamiento negativo) y los caudales turbinados en la planta
hidroeléctrica Zongo.

1 Se nota que los caudales de entrada toman a veces valores negativos. Esto significaria que
durante estos dias hay otras salidas de la represa que los caudales turbinados en el planta
1 hidroeléctrica Zongo. Como no hay desbordamiento, sc podria atribuir estas pérdidas a la
evaporaci6n. Sin embargo, hay otra explicaci6n mas plausible que los valores negativos. Se
observa que estos valores coinciden con dias durante los cuales las variaciones de stock de la
1 represa y los caudales diarios turbinados en Sainani son altos (Figura 23). Esto se explica por
el hecho que, cuando hay pocos caudales que provienen dei deshielo de los glaciares, se debe
soltar mas agua de la represa Zongo para mantener la producci6n de e1ectricidad, sobre todo
1 en ausencia de precipitaci6n. Pero si los caudales de salida y la variaci6n de sotck son
elevados y las entradas disminuidas, un minimo error relativo en el calculo de la variaci6n de
stock puede conducir a un valor negativo en entrada. Un error como este se puede producir
1 por ejemplo cuando la medici6n de nivel se realiza mas temprano 0 mas tarde que la hora
prevista.

1
1
1
1
1
1
.. ,.
1 60000
40000 .
l'
0

11
1 0 Il' "'
1· 011
1 r Il
. 0

1 20000 . 4

o: ...
::l

1 8
o
'ê'E

1 -40000 '

·60000
00

" ...
10

12
00

1
~o

-80000 1 " 14
2-mai J7-mai I-juin 16-juin I-juil l6-juil 31-juil 15-août 30-août 14-sept 29-sept

1 - - Entrées à la retenue Zongo

--Chungem.:nt de stock du lac Zongo
... _... Débits turbinés à J'usine Zongo (sorties)
_Précipitations

1 Figura 23 Caudales de en/rada, desalmacenamien/o de la represa Zongo (signo posilivo),

caudales turbinados en la planla hidroeléc/rica Zongo y precipilaciones. Valores en mm / dia.

1 Las entradas en la represa Zongo pueden ser un excelente indicador de la producci6n de agua
de deshielo de todos los glaciares en la cuenca de drenaje, pero no contienen informaci6n
sobre los valores absolutos deI deshielo. Esto debido a que las superficies glaciares, cuya agua
1 de deshielo alimenta la represa Zongo, no contienen sino una pequena parte deI total de las
superficies glaciares contenidas en la cuenca de drenaje de Sainani (4.3.1, Figura 1). En
cambio, se puede esperar que la evoluci6n de la producci6n de agua de deshielo de esta parte
1 de las superficies glaciares siga mas 0 menos la deI conjunto de los glaciares, porque los
glaciares que alimentan los canales de derivaci6n de la represa cubren una arnplia gama de
exposiciones y de alturas. Para utilizar como un indicador las entradas en la represa Zongo, se
1 puede atribuir un signo positivo a un aumento de los caudales de entrada de un dia al otro y de
un signo negativo a una disminuci6n. Si se hace la misma cosa con los caudales en Sainani, se
pueden comparar cualitativamente las evoluciones deI deshielo (las entradas en la represa
1 Zongo) y de los caudales en Sainani. El resultado de la comparaci6n entre los caudales
naturales en Sainani y el indicador, y entre los caudales turbinados y el indicador, para los

1 periodos en estaci6n seca 1997 sin precipitaci6n, se presenta en la Figura 24.

Se puede observar que, en ausencia de precipitaci6n, los caudales naturales reconstituidos en

1 Sainani (cf. 4.3.1) siguen mejor el indicador deI deshielo total de los glaciares que los
caudales turbinados. Lo que uno observa visualmente es confirmado por las cifras. Los dias de
coincidencia, es decir los dias en los cuales las variaciones del caudal diario y deI indicador
1 deI deshielo estan en el mismo sentido (positivo 0 negativo), se elevan al 78% de los dias
considerados para los caudales naturales reconsituidos y solamente al 56% para los caudales
t\Jrbinados. De 42 dias en los que la evoluci6n de los caudales turbinados no coincide con la
1 deI deshielo, 27 dias tienen una coincidencia con los caudales reconstituidos (64%). En
cambio, de 53 dias durante los cuales los caudales turbinados coinciden, solamente 6 dias no
muestran coincidencia con los caudales naturales (Il %). 78% de coincidencia para la
1 evoluci6n de los caudales naturales reconstituidos es una proporci6n significativa. En fecto,
hay que tener en cuenta el hecho que no se puede esperar una coincidencia perfecta de las

1
1
vanaClOnes durante todo el periodo de observaciôn: es totalmente posible que el
1 comportamiento dei deshielo de todos los glaciares difiera dei de las superficies glaciares que
proporcionan las entradas a la represa Zongo.

1 Périodo 26 mayo à 13 junio 1997 (un dia por pla a)

Caudales de - - + - + + + - + - + - + - - - + - -
1 entrada
Caudales - - + - + - + - + - + - + - + - + - +
naturales
1 Caudales
turbinados
hi:, - - - + - - + + - - + + - + :', + '
+

Périodo 14 'unio à 2 'ulio 1997

1 Caudales de
entrada
+ - + - + - - - + + - - - + - + - - +

Caudales - - + - - - - - + - - - - + - + - - +
1 naturales
Caudales - + .+ - + - - - - + + - - - - + + ~:
+
turbinados
1 Périodo
Caudales de
3 Julio à 21 'ulio 1997
- + - + - + + + - + - - + - + - + - +
entrada
1 Caudales
, naturales
+ - - + - + - + - - - + + - - - + - +

Caudales + + - - + - - + + - - + - - +< -
,~),
+
1
"

turbinados ,
'

Périodo 22 julio à 3 agosto 1997 15 agosto al 20 agosto

1997
1 Caudales de
entrada
- - + - + - + + - - + + - + - + - + +

Caudales - - - - + - - + - + - + - + - + - + +
1 naturales
Caudales - - + + - - + + - + - - - + - - + + ,
-
turbinados

1 Périodo
Caudales de
21 agosto à 8 sept 1997
+ - + + - - + + - + - + + - + + + - +
entrada

1 Caudales
naturales
+ + - - -
+ + + - + - +
-
- - + - + - +

Caudales - + - - + - + - - - - + - + - - + - .+

1 turbinados
Figura 24 Aumento y disminucion de un dia al otro de los caudales turbinados en
Sainani, de los caudales naturales reconstituidos y de los caudales de entrada a la represa
1 Zongo. Si el caudal diario ha aumentado de una dia al otro, el signo es positivo, si ha
disminuido, negativo. En gris: coincidencias entre el caudal y el deshielo. Para la
comparacion se han elegido los dos largos periodos sin precipitacion de la estacion seca
1 1997.

1 Los resuitados deI anâlisis anterior son alentadores. Confirman que nos encontramos en una
buena via para el método de reconstituciôn de los caudales naturales en Sainani. A fin de
poder reconstituir los caudales provenenientes dei escurrimiento de las vertientes, que podrân

1
1
1 ser utilizados en una modelizacion hidrologica, se debell conocer los aportes de agua de
deshielo de los glaciares. Una modelizacion de estas aguas de deshielo constituye la proxima
etapa.
1 Se vera en el proximo capitulo que esta tarea no es simple. Se presenta un esbozo de la
modelizacion dei deshielo total de los glaciares, el cual se basa en las diferencias topograficas
1 entre los glaciares. Para este estudio de la inf1uencia de la topografia sobre el deshielo de los
glaciares se ha elegido el glaciar Zongo y el glaciar Charquini, para los cuales se tienen

1 mediciones de los caudales de deshielo.

1
1
1
1
1
1
1
1
1
1
1
1
1
1
,-
1 .

1
5. Esbozo de una modelizacion dei agua de deshielo de todos los
1 glaciares

1 5.1. La modelizaci6n dei desbielo de los glaciares

5.1.1. Introduccion
1 Recordemos que el objetivo deI programa al cual este trabajo da su contribuci6n es el estudio
deI funcionamiento hidrol6gico de la cuenca de drenaje dei Zongo y el perfeccionamiento de
un modelo hidrol6gico. Este modelo debe servir para hacer predicciones sobre los cambios deI
1 régimen hidrol6gico en el valle, bajo diferentes escenarios climaticos. Evidentemente, el
cambio climâtico no influye solamente sobre el régimen de las vertientes, sino también sobre
el régimen de los glaciares. El estudio deI funcionamiento de los glaciares en el valle deI
1 Zongo comenz6 en 1991 con la instalaci6n de instrumentos de medici6n sobre el glacial'
Zongo, que cubre una ladera deI Huayna PotOS!, la montafia mas alta de la cuenca de drenaje

1 (Figura l, 1.2.1). Las mediciones de diferentes parametros climatol6gicos en el paso de

tiempo horario permitieron establecer y perfeccionar modelos deterministas deI
funcionamiento deI glacial', algunos de los cuales ya marchan bastante bien (cf.

1 RIGAUDIERE, P. el al., 1995). Se deberia disponcr dentro de poco tiempo de un modelo que
permita hacer predicciones deI funcionamiento dcl glacial' Zongo bajo diferentes escenarios
climaticos. No obstante, a fin de poder reconstituir los caudales naturales que provienen deI
1 escurrimiento de las vertientcs, la contribuci6n deI agua de deshielo de todos los glaciares de
la cuenca de drenaja debe sel' conocida (cf. Figura 10, 3.1. 1).

1 Suponiendo que el comportamiento deI glacial' Zongo podra sel' modelizado bajo diferentes
escenarios climaticos, el problema consiste en extender este comportamiento a otros glaciares
de la cuenca de drenaje. En este capitulo, se presenta un enfoque posible para llegar a este fin.
1 5.1.2. Las hipotesis de un enfoque topografico
Los caudales reales de un glacial' dependen de la interacci6n de la superficie deI glaciar con la
1 atm6sfera y de su funcionamiento fisico interno. Un modela determinista intenta describir los
procesos fisicos reales de la manera mas exacta posible. Se utiliza un numero limitado de

1 datos de entrada para describir las condiciones climatol6gicas en la superficie dei glacial' y de
parâmetros fisicos y conceptuales para simulaI' los procesos fisicos. En el glacial' Zongo,
diversos factores climatol6gicos, como la temperatura, la humedad relativa, la veIocidad dei

1 viento, etc., son medidos a diferentes distancias de la superficie deI glacial' (cf. BERTON P. el
al.). El modelo determinista que se ha desarrollado para el glacial' Zongo utiliza estas
mediciones climatol6gicas, asi como mediciones de temperatura en la capa superior de hielo.

1 Como datos de entrada (cf. RIGAUDIERE, P. et al., 1995). A fin de poder simulaI' la
producci6n de agua de deshielo deI conjunto de los glaciares a partir de los caudales medidos
o modelizados deI glacial' Zongo, es necesario conocer la repartici6n espacial de estos datos de
1 entrada deI conjunto de los glaciares y los parametros fisicos y conceptuales de su
funcionamiento interno.

1 La ausencia de mediciones climatoligicas en la mayoria de las superficies glaciares hace

necesaria la estimaci6n de las condiciones en la interfase de estas superficies y de la
atm6sfera. El unico enfoque posible consiste en explicar la variabilidad de estas condiciones
1 pOl' la topografia. Esto quiere decir que se debe hacer la hip6tesis que las condiciones
climatol6gicas en la superficie de los glaciares s610 dependen de la topografia, siendo la

1
1
1 situaci6n climatica idéntica en todos los glaciares. En el modelo de Riguaudière et al., el
estado de la capa superior de hielo también influye en el deshielo. El estado de esta capa se

1 caracteriza, por una parte, por su pertil de temperatura deI cual depende su recalentamiento 0
su enfriamiento, y, por otra parte, por su edad que influye en la absorci6n de la radiaci6n solar
(albedo). En ausencia de mediciones sobre los otros glaciares, se debe igualmente hacer la

1 hip6tesis que este estado s610 depende de la situaci6n topografica de la unidad de superficie y
deI tiempo, que tiene el mismo efecto sobre el conjunto de los glaciares.

1 La extensi6n de los parametros deI funcionamiento interno deI glaciar Zongo al conjunto de
los otros glaciarcs cs todavia mas dificil, inclusa imposible. Esto sc debe a que, en el modelo
de Rigaudière, los caudales producidos en la superficie dei glaciar son transferidos en el
1 cierre, a través de un modelo conceptual, a dep6sitos, cuyos parametros deben ser regulados
con los caudales no disponibles para los otros glaciares. En el modelo de Riguadière et al., la
sectorizaci6n deI glaciar y la regulaci6n de los parâmetros de los dep6sitos de cada sector
1 conducen a tomar en consideraci6n el glaciar en forma tridimensional, es decir la posici6n
relativa deI sector en el glaciar (arriba, al medio, abajo, etc.) y la profundidad deI hielo
(mientras mas espeso es el hielo debajo de un sector, mayor es el dep6sito). Para un modelo
1 de escurrimiento deI conjunto de glaciares, la toma en consideraci6n tridimensional de los
glaciares probablemente no es necesaria. En un primer enfoque, se puede admitir la hip6tesis
que el transito de las aguas de deshielo producidas en las superficies deI conjunto de los
1 glaciares es en promedio similar al transito observado en el glaciar Zongo. En consecuencia,
se puede dejar de lado la posici6n relativa deI lugar de producci6n de agua de deshielo al

1 interior deI glaciar.

Las hip6tesis para un modelo preliminar son:

1 1) La situaci6n climatol6gica, principalmente la cobe11ura nublosa deI cielo y la temperatura a

gran escala, es idéntica en todos los glaciares de la cuenca.

1 2) La topografia es el unico factor que determina el microclima en los glaciares en una

situaci6n climatol6gica dada.
3) La topografia es el unico factor que determina el estado de la capa superior de hielo (nieve)
1 deI glaciar y el efecto temporal sobre el estado deI hielo es idéntico en todos glaciares.
4) El tiempo de transita de las aguas de deshielo desde el lugar de producci6n hasta la salida
es, en promedio, el mismo en todos los glaciares y es el que fue observado en el casa deI
1 glaciar de Zongo.

Resumiendo estas cuatro hip6tesis en una sola:

1 • La· relaci6n entre los escurrimiento deI glaciar Zongo (medidos 0 simulados) y el
escurrimiento deI conjunto de los glaciares es unicamente determinada por las condiciones
1 climaticas de superficie de los glaciares que no dependen sino de la topografia local de las
mallas de superficie.

1
1
1
1 5.1.3. La importancia de la exposicion
Se observa que, en el hemisferio sur, los deshielos de los glaciares expuestos al norte se
1 encuentran generalmcnte a una altura mas devada que los frentes de los glaciares expuestos al
sur. Sus dimensiones son generalmcnte mas pequenas y retroceden mas rapido (cf.
FRANCOU, B. et RlBSTEIN, P., 1995). Es la consecuencia de la ablaci6n elevada de las
1 superficies glaciares expuestas al norte en relaci6n a las expuestas al sur, debido a una
intensidad de la radiaci6n solar globalmente mas elevada. La posici6n deI sol a cierta altura
varia en el transcurso deI dia y en el transcurso deI ano.
1 En la latitud de 16° l' sur donde estan situadas las superficies glaciares de la cuenca de drenaje
deI Zongo, el sol se encuentra casi todo el ano al norte. Por consiguiente, las laderas expuestas
1 al norte reciben en promedio mas radiaci6n solaro Pero s610 la insolaci6n media influye en el
deshielo de los glaciares. La distribuci6n de esta energia, segun las horas deI dia, también
debe ser importante. Asi, es probable que la insolaci6n de la tarde produzca mas agua de
1 deshielo que el de la manana, la humedad aumenta durante el dia a causa deI ascenso de masas
de aire humedo que vienen deI valle bajo: esta tiene como consecuencia una disminuci6n de la

1 sublimaci6n deI hielo en provecho deI deshielo.

La observaci6n de que los glaciares de cierta exposici6n retroceden mas rapido hace pensar

1 que la exposici6n debe ser une de los factores mas importantes de la topografia de los
glaciares que influye en su deshielo. En un modelo basado en las hip6tesis hechas en el
parrafo anterior, las influencias de las caracteristicas topograficas en el deshielo pueden ser
1 estudiadas en cada malla independientemente de las otras. El deshielo total en este tipo de
modelo corresponde al integral deI deshielo de una malla sobre el conjunto de mallas que
contribuyen agua de deshielo. Es decir, si se llega a describir la influencia de la topografia
1 sobre el deshielo de una malla, se puede integrar la relaci6n encontrada en una superficie y en
una duraci6n determinadas a fin de obtener un factor topografico para el total deI agua de
deshielo de un glaciar. Esto se hara a continuaci6n.
1 5.1.4. El concepto y la base fisica de un modelo a desarrollar
En el modelo de Rigaudière et al. (cf. RIGAUDIERE, P. et al., 1995) la producci6n de las
1 aguas de deshielo en la superficie deI glaciar Zongo se base en el balance de energia. Los
términos de este balance son presentados en la Ecuaci6n 1). Dependen deI microclima en el

1 glaciar que es la resultante de la situaci6n climatol6gica general y de la topografia, asi coma

deI estado de la capa superior de hielo (0 de nieve). Con las hip6tesis hechas en el pârrafo
anterior, basta conocer la influencia de la topografia de los glaciar en estos términos para

1 modelizar los escurrimientos deI conjunto de los glaciares a paI1ir de los caudales medidos 0
modelizados deI glaciar Zongo.

1 La topografia de los glaciares posee dos aspectos. Estos son, por una parte, los parametros
topograficos de cada malla descriptiva deI glaciar, es decir la altitud, la exposici6n y la
inclinaci6n, y, por otra parte, la topografia a gran escala. Esta ultima concierne sobre todo la
1 circulaci6n atmosférica y por consecuencia el término de intercambio de calor sensitivo, pero
también la radiaci6n solar directa tomando en cuenta las sombras proyectadas sobre partes
mas 0 menos grandes deI glaciar. En un primer enfoque se puede hacer la hip6tesis que los
1 movimientos de aire son en promedio idénticos en todas Jas mallas de los glaciares y en
consecuencia la evoluci6n deI término deI calor sensitivo es la misma en todos los glaciares.
A excepci6n deI efecto de las sombras sobre la radiaci611 solar directa, los otros términos dei
1
1
1 balance no son influidos por la topografia a gran escala y se puede lIevar el estudio a la

1 influencia de los panimetros topogrâficos de cada malla.

Para nuestro futuro modelo, se puede separar la suma de los términos de derecha deI balance

1 (ECUACION 1), que corresponde a la energia disponible para el deshielo y para la

evaporaci6n, en dos grupos.

1
q Caudal especifico de deshielo (m J/s/m 2 )

1 HI
P
Calor latente de fusion (2.5 * 10 6 J/kg)
densidad dei agua (kg/m 3 )
H. Calor latente de sublimacion de la nieve/hielo (3.3 *1 0 5 J/kg)

1 e
a
[Link]'
Masa de agua evaporada por unidad de tiempo y de superficie (kg/s/m 2)
Albédo, sin unidad
Radiacion solar directa (W1m 2)
Rdif Radiacion solar difusa (W1m 2 )
1 Ralon
Rsul
Radiacion atmosferica
Radiacion deI suelo (gran longilud de onda) (W1m 2 )
Qscns Flujo de calor sensible (W1m 2)

1 Qprof
Qprec
Q'ech
Calor aportado por las precipitaciones (W1m ) 2
2
Flujo de calor creada por las capas illferiores (W/m )

Calientamiento/Enfriamiento dei manto de nieve (W/m 2)

1 Ecuacion 1 Balance de energia en la superficie dei glaciar (segun LIBOUTRY, 1964)

El primer grupo comprende los términos de energia que son influidos par la situaci6n

1 topografica de la malIa de calculo. Contiene la expresi6n de la radiaci6n solar (l-a)* (RJir +

RJif), de la radiaci6n atmosférica (Ra1m ) Y el término de los aportes de las capas inferiores
(Qpror). El segundo grupo contiene todos los términos cuya evoluci6n en el tiempo es idéntica

1 en todas las mallas. Los términos pueden ser influidos por las condiciones climaticas
generales 0 por el estado general de la capa superficial de hielo, pero depcnden poco 0 nada de
la altura 0 de la exposici6n de las mallas. Este es mas 0 menos el caso de los términos de la
1 radiaci6n deI suelo (Rsuelo), deI calor sensible (QSCIlS), deI aporte de calor de las precipitaciones
(Qprec) y deI recalentamiento/enfriamiento de la capa superficial de hielo (nicve) (Qrcc) (cf.
RIGAUDIERE, P. et al., 1995). Como en la aplicaci6n deI modelo de Rigaudière et al. Se
1 observan variaciones de R a1m Y de Qprof relativamente pequefias en comparaci6n a las de la
energia total, se las puede asignar al segundo grupo. De esta manera, se puede describir la
energia disponible para el deshielo y la evaporacion para una malla dada de la siguiente
1 manera:

1 Etat = q* Hj*p + Hv*e = LE topo + LE prop (2)

1 y
LE topo = (l-a)*Rdir(exp) + (l-a)*RJif
LEprop = Qscns + Qpréc - Qrcch - Rso1 + R atm + Qprof
1
1
1
1 d'onde:

1 a Albedo, que depende de la altura (ait), dei angulo incidente dei sol
(exp), de la edad dei mante nevoso (t) y de la cobertura nublosa dei
cielo
1 Rdir(exp) la radiaci6n solar directa, que depende dei ângulo incidente dei sol

1 Como se trata de estimar finalmente la cantidad de agua de deshielo producida en la superficie

de los glaciares, s610 hay interés en la parte de la energia total disponible que no es utilizada
para la sublimaci6n. Arriba de la isoterma cero grade de los glaciares, no hay deshielo: toda la

1 energia disponible es utilizada para la sublimaci6n. En cambio, bajo la isoterma cero grado, la
mayor parte de esta energia es utilizada para la fusi6n de nieve. Ya que si el deshielo es
posible (temperatura superior a cero grado), es utilizada diez veces menos energia para el

1 deshielo que para la sublimaci6n de una misma cantidad de hielo/nieve. Basta pues con tomar
en cuenta las mallas abajo de esta isoterma cero grado, que se desplaza en el transcurso deI
dia. Abajo de esta isoterma, la cantidad de hielo sublimado depende dei gradiente de la
1 humedad relativa, que es mas elevado y presenta una mayor variabilidad en estaci6n seca que
estaci6n de lIuvias (cf. RIGAUDIERE, P. el al., 1995). Un gradiente elevado favorece el
alejamiento de masas de aire producidas por la sublimaci6n, 10 que aumenta el fIujo deI calor
1 latente (la sublimaci6n). En estaci6n seca, las masas de aire humedo, que en la tarde suben de
la Amazonia hacia la alta cuenca de drenaje, disminuyen los gradientes de humedad y en
consecuencia la sublimaci6n. Se puede esperar pues que en la tarde es utilizada una mayor
1 parte de la energia total disponible para el deshiclo que en la mafiana. Aunque el fIujo de calor
latente se atenua un poco con la altura, se puede suponer que estas variaciones son pequeiias.
Esto permite aceptar la hip6tesis que la parte de energia utilizada para el deshielo depende
1 unicamente de la hora deI dia.

1 F(h) = (q*Hv*p)/(q*Hlp+ Hv*e) (3)

El caudal instantaneo constituido por el agua de deshielo de una malIa puede escribirse a

1 partir de las ecuaciones 2 y 3 :

q*Hr*P= F(h)* (l-a)*Rdir(exp) + F(h)*(1-a)*Rdir+ F(h)*LEprop (4)

1
1 Planteando para el albedo un valor correspondiente a una altura media, se integra (4) en la
superficie de ablaci6n variable :
1
1 Q, . Htp = (1- am(h»· F(h)[.f Rdi,(exp)· dS + Rd!! . Sr (h)] + F(h)·
.~ r
~E prop ·Sf(h) (5)

1
donde:
1 Qs = caudales de deshielo instantaneo deI glaciar entera

1
1
1 St{h)
am(h)
Superficie bajo la isoterma cera grado en una hora dada del dia
= Albedo que corresponde a la altura deI centra de gravedad de la
1 superficie bajo la isoterma cera grado, en una hora dada deI dia y un
éingulo de incidencia nulo

1 Veamos el primer término de esta ecuacÎôn. Tanto la parte de la radiaciôn solar reflejada (el
albedo) coma la radiaciôn solar directa dependen deI ângulo de incidencia deI sol (ex.), medido
en relaciôn a la normal de la malIa (cf. Figura 30). Mientras que la radiaciôn disminuye en

1 proporciôn al coseno de este ângulo, la absorciôn de esta radiaciôn aumenta, es decir que hay
menos energia reflejada por la superficie (1-a aumenta). Si se pone la = cos (ex.) y se toma en
consideraciôn la disminuciôn deI albedo se puede escribir para la radiaciôn solar directa
1 absorbida:

(1- a", (h)) . Rdir(exp) = (1- a",(h))· R{~~;'X f / (6)

1 <.lan<.le:

1 R;::X = Radiaciôn solar directa para un ângulo de incidencia nulo,

constante en el tiempo (aproximadamente 25% de la radiaciôn
1 x =
extraterrestre, cf. MUSY, A. y SOUTTER, M., 1991)
Factor de potencia, entre 0 yi, constante

1 Para obtener el caudal diario, se debe integrar la ecuaciôn (5) en el dia. Teniendo en cuenta la
variaciôn horaria de la energia disponible, se puede integrar en dos medio dias, la mafiana
hasta mediodia y la tarde hasta la noche. Esto dada:
1
1 Q' HJP = LF;[O- am,;}' R~~:x JfaxdS + (1- am,;)' RdirSr(i)] + l F;(L EpropSJ(i))
1 Sfli) 1
(7)

1 donde:
Q = caudal diario
= indice deun medio dia: 1= mafiana a mediodia, 2 = mediodia a la
1 am,i
tarde
= Albedo que corresponde al angulo de incidencia nulo y a la altura deI
centro de gravedad de la superficie media bajo la isoterma cero
1 Sr(i)
durante un medio dia.
= Superficie media bajo la isoterma cero grado durante un medio dia.

1
1
1
1
1 Partiendo de la ecuaci6n (2) se han reemplazado los términos que dependen de la exposici6n y
de la altura dei pixel, la radiaci6n solar directa (Rdir) Y el albedo (a), por funciones dei coseno
1 dei angulo de incidencia dei sol (la) y valores medios (am,.) 0 maximos (R ) de estos
términos. Estos ultimos ya no dependen de la topografia ùe los glaciares y su evoluci6n en el
tiempo es idéntica en el conjunto de los glaciares, como los términos Eprop . Esto permite
1 reemplazar estas dos energias por un factor que exprese la parte de la energia de la radiaci6n
directa de la energia total:

1
1
Si se reemplaza esta ecuaci6n en la ecuaci6n (7), y suponiendo que la influencia dei angulo de
1 inciùencia sobre el albedo es insignificante y que el albedo medio es el mismo por la manana
y por la tarde, se puede asimilar la radiaci6n difusa (Rdir) a los términos que no dependen de la
topografia y escribir:
1 m j)·E tot [f .J
Q'= (l-aH '.. '~Fj Ri' ladS+(l-Rj),Sf(i) (9)
1 y
f p 1 Sf{i)

1 fladS ::.: I(Sr,i) (10)

SI (1)

de 10 cual:
1 I(Sr,i) ::= Sumas de cosenos dei angulo de incidencia dei sol sobre las mallas
de la superficie bajo la isoterma cero. Es una medici6n de la
1 radiaci6n solar directa recibida por el glaciar en funci6n de su
topografia

1 La ecuaci6n de los caudales (9) derivada dei balance de energia Cl) es un modelo simple que
toma en consideraci6n la influencia de la topografia sobre los caudales de los glaciares. Esta
muestra que un enfoque 'topografico' no puede funcionar con factores de exposici6n. En
1 efecto, I(St;i) es el unico término influido por la topografia dei glaciar. Es una medici6n de la
radiaci6n recibida por la superficie bajo la isoterma cero grado durante un paso de tiempo de
la integraci6n en el dia (por ejemplo durante un medio dia) y en funci6n de la topografia de las
1 mallas de superficie a las cuales unD integra. Se ve en la ecuaci6n (8) que la influencia de este
término sobre el caudal no es constante, pero que depende de la energia total recibida por el

1 glaciar y de la parte de la radiaci6n solar directa maxima (R ) de este total (Ri).

El término I(Sf,i) hace intervenir la topografia de malias de superficies de los glaciares en la

1 producci6n de agua de deshielo, Corresponde a la integral en una superficie variable en

funci6n .dei tiempo dei coseno dei angulo de incidencia dei sol. Mide la radiaci6n directa
recibida por la superficie variable bajo la isoterma cero durante un paso de tiempo variable. Se
1 ha elegido el dia entera para este paso de tiempo. Pero con el método que se ha desarrollado
para su calculo, presentado a continuaci6n, es fâcilmente posible calcular la integral para un
paso de tiempo cualquiera. Hay que observar que cada paso de tiempo en un dia hace
1 necesario la fijaci6n de un parametro que determine la parte de la energia disponible que es

1
1
1 utilizada para el deshielo (Fi). La variante mas simple deI modelo desarrollado en este parrafo
divide el dia en dos pasos de tiempo, como ya se ha mencionado, para tener cuenta de la

1 conversion de energia disponible de agua de deshielo que es mas elevada por la tarde que por
la manana.
5.2. Comparaci6n de caudales de deshielo y de topografias de los glaciares
1 Charquini y Zongo

5.2.1. Los glaciares de referencia

1 Por el momento no se dispone de medicion de caudales sino para los glaciares Zongo y
Charquini. Un limnigrafo acaba de ser instalado al pie deI glaciar de la vertiente norte deI
Huayna Potosi, pero atm no se dispone de mediciones explotables. La situacion de ambos
1 glaciares, Zongo y Charquini, asi como las instalaciones limnimétricas y meteorol6gicas se
ubican en la Figura 25. Los caudales deI emisario deI glaciar Zongo son medidos justo al pie

1 deI glaciar y la parte de la cuenca de drenaje no glaciar solo comprcnde las morrenas,
correspondiendo al 30% de la superficie total captada (0.9 km 2 de 3.05 km 2). En cambio, la
parte no glaciar de la cuenca de drenaje deI Charquini es mucho mas importante, abarca casi el

1 70% de la superficie total captada. Las dos tomas de agua instaladas en los emisarios deI
Charquini que alimentan el canal Prado, en el cual son medidos los caudales, no captan la
totalidad de los escurrimientos glaciares, en particular en estaci6n de lluvias. Ademas, los

1 caudales medidos en el canal Prado incluyen los caudales provenientes deI

desalmacenamiento deI lago Mamankota. Esto limita las posibilidades de comparacion de los
caudales deI canal Prado con los deI emisario deI Zongo en periodos sin lluvia y sin
1 operaciones hidraulicas sobre el lago Mamankota.

1
1
1
1
1
1
1
1
1
1
1
1
. ",,'

1 ...................................:; :/ ;" /~. / /

1 (
.
.0
-o ,/)
]l5200

.........,...- ~.:::..'" 'f';

1
Or)

1
1 Glacier Zongo
'Station limnigraphique

1 115200 Station météorologique (altitude)

~.•••.........
.,, Limites des bassins versants des lilllnigraphes
1 Limite du bassin versant du Zongo

1 Glacier, surfaces au dcsous d'une altitude

Cours d'eau naturel ~ Canal, prise d'eau

Huayna Potosi
6050 m.s.m.

1 Figura 25
Direction principale de l'exposition

Plan de siluacion de los glaciares Zongo y Charquini con las inslalaciones

limnimélricas y melcorol6gicas
1 Traduccion de la leyenda :

1 Estaci6n limnigrâfica
Estaci6n meteorol6gica (altura)
Limites de las cuencas de drenaje de los limnigrafas
1 Limite de la cuenca de drenaje deI Zongo
Glaciar, superficies par debajo de una altura
Rio natural Canal, toma de agua
1 Direccion principal de la exposici6n

1
1
1
1
1 5.3. Los datos obtenidos con ayuda dei MNT

1 5.3.1. Los MNT de los glaciarcs

Para obtener los datos topogrâficos de los glaciares, se sobrepuso una imagen satelital SPOT
de la estaci6n seca 1989 con el modelo numérico de terreno (MNT) de la cuenca de drenaje,
1 obtenido con ayuda de mapas geogrâficos militares de escala 1:50'000. Como los frentes de
los glaciares tropicales retroceden râpidamentc, se ha tratado de utilizar la imagen mas
reciente posible. La sobreposici6n fue realizada con cl programa de tratamiento de imâgenes
1 ERDAS 'Image'. En la imagen 'SPOT' las superficies glaciares son visibles, con pixeles de
un color azul claro u obscuro. Por superposici6n se pudo delimitar los contornos de los
glaciares. Se obtiene de esta manera los M}\fT deI glaciar Zongo y deI Charquini, asi como, a
1 manera de comparaci6n, deI conjunto de los otros glaciares. A modo de ejemplo, el MNT
obtenido para el glaciar Zongo se presenta en la Figura 26.

1 Un examen minucioso de los mapas geogràficos militares, gracias a los cuales se obtuvo el
MNT de la cuenca, muestra que las curvas de nivel de estos mapas probablemente fueron
trazadas por interpolaci6n entre puùtos de altura conocida. Por 10 tanto, el MNT de los
1 glaciares no es muy preciso y la topografia real es mucho mas accidentada.

1
1 6100
-5900

1 5700
5500 g
Cl)
5300
1
"Q

5100 'E-.;

1
1 axe csl-oucsl
200m

allC nord-sud

1
Figura 26 Modelo numérico de terreno (MNT) dei glaciar Zongo. La resoluci6n
1 elemental es de 100 m.

1
1
1
1
1 5.3.2. Los parametros topograficos obtenidos con ayuda dei MNT para cada malla
La superficie de una malla dei MNT proyectada sobre el piano horizontal es de 0.01 km 2

1 (lOO*IOOm). La suma de todas las mallas dei glaciar Zongo corresponde a una superficie de
2.15 km 2, 10 que coincide con la superficie total determinada por B. FRANCOU (FRANCOU,
B., et RIBSTEIN, P., 1995). La superficie dei glaciar Charquini es de 0.46 km 2 Y la dei
conjunto de los glaciares de 18.3 km 2. Sin considerar el valle de Coscapa, la superficie total de
1 los glaciares es de 11.5 km 2• Hay que observar que el 37% de las superficies glaciares de la
cuenca dei Zongo pertenecen al valle de Coscapa, cuyos escurrimientos no pasan por la planta

1 hidroeléctrica de Sainani. A partir de alturas de cada malla de un glaciar, se pueden estab1ecer

las curvas hipsométricas de los glaciares. Las curvas hipsométricas en los glaciares Charquini
y Zongo yen el conjunto de los glaciares (sin Coscapa) se presentan en la figura 27. La forma

1 de la curva dei conjunto de los glaciares es la mas tipica: una parte empinada justo al frente
dei glaciar, luego una pendiente menos importante y de nuevo un aumento hasta los valores
maximos observados en las zonas de las cumbres. El glaciar Zongo tiene sensiblemente el
1 mismo perfil, mientras que el glaciar Charquini tiene un perfil particular. La inflexion en la
curva dei glaciar Charquini, entre 4900 y 5100 m de altitud se debe a algunas mallas cubiertas
de hieIo aisladas y no contiguas al glaciar principal.
1
1 6000 1

i
--Glacier Charquini

1 5800
- - Glacier Zongo
. ". "... Tous les glaciers sans Cos capa 1
1

1<Ii 5600
..
1 É
--- 5400
.g
deux parties non adjacentes
du glacier Charquini
oE
:( r--- ~~" .. ~ -~

1 5200 : ~_~~~:::::::::::::::....__

5000 . ".
1 4800 l'

1 100% 90% 80% 70% 60% 50% 40% 30% 20% 10% 0%

Figura 27 Curvas hipsométricas de los glaciares Charquini, Zongo y dei conjunto de los
1 glaciares de la cuenca de drenaje (sin los glaciares dei valle Coscapa)

1
1
1
1
1
1 6000 1

1
... t : : : : : ~ : : : :

...
~ : ~
: ..: .. .:.: ,,".,_

...........
...... . - Glacier Charquini
5800 ' - - Glacier Zongo

1 ""
E 5600
....... :-
....... Tous les glaciers sans Coscapa ;

<Ii ..... : : : : : : ; ~ .
. ". ~. -. ". -..........
1 ~
E
'0'
·Ë
54001 1

'
". :::: : ..
......... .- .... -
-:: ..
< 5200
1 .

....... • 1
. ...
.:
5000 ........
1 4800 1 _. - --~~-
!

1 Figura 28
nord nord-est est sud-est sud

Exposic.:i6n media de las mal/as de lOOm*lOOm para c/ases de 25 m de

sud-ouest ouest nord-ouest nord

diferencia de a/tura
1 El programa de tratamiento de imagenes utilizado permite calcular las exposiciones y las
inclinaciones maximas de las mallas. La exposici6n de las mallas corresponde a la direcci6n
1 de la inclinaci6n maxima. Las exposiciones medias de los glaciares se presentan en la Figura
28. Se observa que el glaciar longe pasa de una exposici6n este en su parte inferior, hacia una

1 exposici6n sureste en su parte superior. Un cambio comparable de la exposici6n media de la

parte inferior hacia la parte superior deI glaciar puede observarse en el Charquini. Pero la
exposici6n deI Charquini es exactamente opuesta a la deI longe (cf. igualmente Figura 25).

1 La curva de exposici6n media deI conjunto de glaciares se encuentra entre las deI Charquini y
deI longo. Oscila irregularmente entre el sur-sureste y suroeste. Esta tendencia deI conjunto
de glaciares hacia una expasici6n al sur se debe a que los glaciares expuestos al norte son
1 generalmente reducidos a causa de la radiaci6n solar directa, que es en promedio mas fuerte
en las pendientes expuestas al norte (cf. 5.1.3).

1 5.4. La estimaci6n de las diferencias de la radiaci6n solar directa en el

transcurso dei dia y durante el afio en los glaciares Zongo y Charquini
1 5.4.1. Introducci6n

1 En el pIano de la situaci6n de los glaciares Charquini y longe (Figura 25) 0 en la Figura 28 de

las exposiciones medias, se observa que las partes inferiores de estos dos glaciares estân
orientadas de oeste a suroeste (longo) y de sureste a noreste (Charquini); es decir mas 0

1 menos opuestas en una lîna de este a oeste. En consecuencia, la diferencia entre las
insolaciones de estas exposiciones opuestas no debe variar mucha durante eI aiio. Los dos
glaciares difieren no obstante por las horas deI dia donde se encuentran directamente
1 expuestos al soL Las partes bajas deI Zongo reciben el sol por la mafiana, las deI Charquini
por la tarde. Si se consideran los glaciares en su conjunto, la variaci6nde la intensidad solar
difiere en los dos glaciares. Las partes medias y altas deI glaciar longe estan expuestas mas
1 bien hacia eI sur, 10 que hace que reciban una insolaci6n maxima en verano, cuando eI sol se
encuentra al sur de la latitud de los glaciares. A mediodia, es el casa entre el 26 de enero y el

1
1
1 24 de abriI; en las primeras horas de la manana y al final dei dia, entre el 12 de diciembre y el
8 de junio. Las partes medias y altas dei glaciar estan expuestas al sol durante la mayor parte

1 dei dia en estaci6n de Iluvias y particularmente en estaci6n seca, entre principios de enero y
mediados de mayo,

1 Parece que que la posici6n dei sol durante el ano es de considerarse si se quiere estimar las
diferencias de radiaci6n solar directa sobre los glaciares Charquini y Zongo. Por esta raz6n a
continuaci6n se toma en consideraci6n no solamente la evoluci6n de la insolaci6n durante el

1 dia sino también dei ano.

5.4.2. La posici6n dei sol en el transcurso dei aiio y durante el dia

1 La posici6n dei sol en relaci6n con un punto terrrestre es determinada por dos angulos. El
cenit es el angulo incidente sobre la horizontal dei punto, es decir, el angulo entre la !inea
solltierra y la vertical en el punto considerado. El acimut corrresponde a la declinaci6n dei sol
1 respecto al norte (Figura 30), Las f6rmulas de calculo de estos dos angulos en funci6n de las
coordenadas, deI dia y de la fecha, fueron establecidas por Robinson (cf. ROBINSON, N.,
1996). La posici6n dei sol durante el dia y el ano se presenta en la Figura 29.
1 La Figura 29 presenta las trayectorias dei sol durante el dia en los dias mas largos y los mas
cortos de la estaci6n seca, Se ve que al comienzo de la estaci6n el sol sube mas alla de la
1 altura de los glaciares, mientras que hacia el final de la estaci6n, su trayectoria pasa mucho
mas al norte.

1 n

1 ';j est .... ,

' ........' " ..
ë Azimut. 1 mai
'"
'
.
1 :;]
"0
'5
E
,
. \ ......
" ... ' "
........ ....
,
-5
E sud ., ,
. -,
,,,<;:> r.
.- 1/2n '2
~
..,
1 <l.l
"0
c:
.2
Ü
Zénith, 20 sept. - ,
.. ' bO
c
-<X:
e "
ëS
1
ouest .. ,.

.- Zénith, 1 mai

1 nord
o 1 2 3 4 5 6 7 8 9
-"-1

10 Il
-, ••
',-

12 13 14 15 16 17 18 19 20 21 22 23 24
Heure de la journée

1 Figura 29 Posicion dei sol duranle el dia mas largo (l de mayo) y mas corlo (20 sept.) de
-

la eSlacion seca. En gris oscuro: horas de noche dei 1 de mayo,' en gris clar%scuro: horas
1 de noche dei 20 de sepl.

1
1
1
5.4.3. Calculos de los angulos de incidencias dei sol con las mallas de los glaciares
1 La intensidad de iluminaci6n de una malia depende deI ângulo incidente deI sol sobre esta
malla. El poder de la insolaci6n es maximo cuando el angulo incidente, medido en relaci6n a
la nonnal de la superficie, es 0°. Es nula en un ângulo incidente de 90°. Mas alla de 90°, la
1 malia se encuentra en la sombra. La malia puede encontrarse también en la sombra a causa de
la elevaci6n de las malias vecinas. Para saber si una malia se encuentra a la sombra de las
mallas vecinas, hay que conocer los ângulos 'topograficos' entre la malla considerada y las
1 mallas que se encuentran en la alineaci6n deI acimut deI sol. Si el complemento a 90° dei cenit
deI sol es mas grande que el angulo 'topografico' maximo, la malla se encuentra al sol, en el
caso contrario, esta a la sombra. El esquema de la Figura 30 presenta una explicaci6n grafica.
1 La media de los cosenos de los ângulos de incidencia (Ia) en todas las mallas de superficie de
un glaciar es proporcional a la radiaci6n solar directa (medio) sobre esta superficie (cf. 2.2.4).

1 A fin de poder calcular estas medias deI coseno deI angulo de incidencia durante el dia y el
aiio, hay que calcular las caracteristicas de todas las mallas en cada paso de tiempo. Para
simplificar y permitir un calculo rapido, se ha elegido otro acercamiento. En lugar de calcular
1 la media deI coseno de los glaciares a cada paso de tiempo, se calculan las medias de los tres
conjuntos glaciares (Charquini, Zongo y el conjunto de glaciares) de una gama limitada de

1 valores de los angulos (cenit y acimut). Para una hora y fecha dadas, se puede facilmente
detenninar primero los angulos de la posici6n deI sol (con ayuda deI cuadro de los cenit y
acimut en funci6n de la hora y de la fecha, cf. 5.4.2), y luego leer los cosenos en el cuadro

1 establecido l5 • Esto puede parecer complicado pero para detenninar las medias de los cosenos
deI angulo de incidencia en relaci6n a la superficie de los glaciares abajo de cierta altura, este
procedimiento es el unico posible, como se vera en el pr6ximo parrafo. El cuadro fue
1 establecido para 8 angulos de acimut y 16 angulos de cenit.

1
1
1
1
1
1
1
1 IS Como el nûmero de valores de acimut y de cenit en el cuadro es necesariamente limitado, interpolaciones

Iineales son efectuadas entre los valores dei cuadro.

1
1
1 Normale de la maille 'Normale de la terre
Angle maximal
ensoleilé pour la

1 maille considéré

.......
.....

1 ......
........
....

1
1
1
1 ..';<-.-+-------'- , .
..: / ;
l.·······

1 Légende:

1 ! Maille du glacier

Maille du glacier à l'ombre

1 Maille hors du glacier / Limite du glacier

i ··..··! Maille ensolciléc pour laquelle l'angle entre

1 , ,
:.. .: le soleil et la normale est determiné
Maille pour laquelle l'angle zénit maximal
SUD I:.. ..J pour être ensoleilé est determiné

1 Figura 30 Ccilculo de la inlensidad de insolacion de las mallas de los glaciares para una
posicion dada (acimul/cenit) del sol. Delerminacion si una malla se encuenlra en la sombra y
1 ccilculo del cingulo incidenle del sol para una malla (cingulo enlre la normal de la malla y la
/inea sol/lierra)

1 5.4.4. La insolacion media sobre las superficies totales de los glaciares

El resultado dei calculo de las medias de los cosenos dei angulo de incidencia sobre las

1 superficies enteras de los glaciares se presenta en las Figuras 31 y 32. Se puede observaI', pOl'
una parte, la evoluci6n de las medias diarias (medias de los cosenos sobre las mallas de
superficie y las horas dei dia) y , pOl' otra parte, las variaciones durante el dia en los dias de

1 insolaci6n maxima (1 mayo) y minima (10 septiembre). La cone! usi6n de la comparaci6n

cualitativa hecha en el parrafo 0 es confirmada pOl' los datos. Se observa que la insolaci6n
media dei glaciar Zongo esta mas pr6xima que la de Charquini al comienzo de la estaci6n
1 seca, cuando la parte de la cumbre dei glaciar Zongo, expuesta hacia el sur, es insolada,
porque el sol sube mas alla de la latitud de los glaciares durante el dia (cf. Figura 25,5.2.1).
La insolaci6n dei glacial' Zongo disminuye enseguida porque la parte de la cumbre es cada vez
1 menos ilwninada cuando el sol se desplaza hacia latitudes norte durante la estaci6n seca. Es
1
1 interesante observar que el desarrollo deI conjunto de los glaciares (medias de sol en las
mallas de todos los glaciares deI valle) sigue mas 0 menos al deI glaciar Charquini, con un

1 desplazamiento hacia valores un poco bajos. Las mallas de superficie deI conjunto de
glaciares estân en promedio expuestas hacia el sureste, mientras que la exposici6n deI glaciar
Charquini esta al suroeste (cf. Figura 28, 5.3.2). Estas dos exposiciones son simétricas de une

1 y otro lado de una !inea sur-norte, en la cual la inclinaci6n deI sol se desplaza durante el afio.
Esto explica la evoluci6n paralela de la insolaci6n de ambos glaciares durante el ano.

1 La variaci6n de la insolaci6n durante el dia ref1eja la trayectoria deI sol durante el dia de este a
oeste. Como el glaciar Zongo esta expuesto mas bien hacia el este y el glaciar Charquini hacia
el oeste, el maxime de insolaci6n deI glaciar Zongo se desplaza hacia la mafiana y el deI
1 Charquini hacia la tarde. Se ha observado anteriormente que la insolaci6n de la tarde tiene
probablemente un efecto mas importante sobre los escurrimientos que la insolaci6n de la
mafiana, ya que una minima energia es utilizada para la sublimaci6n (0). En un enfoque
1 modelizador, las series de caudales de ambos glaciares son ideales para fijar los factores que
describen la parte de la energia disponible utilizada para la fusi6n deI hielo en el transcurso
deI dia.
1 La variaci6n de la insolaci6n durante el dia sobre el conjunto de glaciares se asemeja mas bien
a la deI glaciar Zongo que a la deI glaciar Charquini, porque el componente este de la
1 exposici6n tiene la misma importancia que en el caso de Zongo (cf. Figura 28, 5.3.2). Las
diferencias de sol durante el dia entre los glaciares Zongo y el conjunto de glaciares se

1 observan sobre todo durante las horas de la tarde y hacia el final de la estaci6n seca (IOde
septiembre), cuando el conjunto de glaciares es mas insolado que el glaciar Zongo. Al
comienzo de la estaci6n, en cambio, la insolaci6n durante el dia es casi idéntica en las dos

1 superficies glaciares. Esto se explica en parte por las sombras proyectadas. Se puede observar
en la Figura 32 que las curvas deI glaciar Zongo estan un poco aplanadas por la tarde, sobre
todo al comienzo de la estaci6n seca (1 mayo). Como las partes de la cumbre deI glaciar
1 Zongo estân expuestas al sureste, entran en la sombra de la cresta sur deI Huayna PotOS! (cf.
Figura 25, 5.2.1). Al comienzo de la estaci6n seca, la curva deI glaciar Zongo esta desplazada
hacia valores mas elevados, a causa deI aumento general de la insolaci6n sobre este glaciar
1 durante el verano. En consecuencia, la insolaci6n matinal deI glaciar Zongo es mas elevada
que en el conjunto de glaciares, mientras que durante la tarde, los valores de las dos
superficies glaciares se vuelven mas pr6ximos a causa deI efecto de las sombras presentes. En
1 invierno, hacia el final de la estaci6n seca, no hay efecto de sombras proyectadas en la parte de
la cumbre deI glaciar Zongo, pues esta zona ya se encuentra en la sombra por la manana. En
cambio, la morrena norte deI glaciar Zongo podria, por la tarde, disminuir la insolaci6n de las
1 partes bajas de este glaciar (cf. Figura 25,5.2.1). El efecto es aumentado por el hecho que s610
estas partes bajas reciben la mayoria de la insolaci6n en invierno. Esta podria ser la
explicaci6n de los valores cada vez mas elevados de la insolaci6n deI conjunto de los glaciares
1 el IOde septiembre en comparaci6n con el glaciar Zongo durante el ano.

1
1
1
1
1
0,65
Comparaison de l'évolution
journahère pour le 1 mai et le
0.60
10 septembre dans la ligure 31
0.55

1 '"
.......
" . . -,
~ _ _~_----/. --: 0,50

'"
- ... -. ; 0,45

1
........ - .... - ••• - •• - ..... - .. - - • • °1

- - Glacier Charqulni ; 0.40

,1
~ -+-__- ,, 0.35
--Glacier longo

1 ,
l,
Tous les glaciers (sans les glaciers du 1
sous-bassin Coscapa) j
0,30

0.25

1 ,. ,. 1
';;; ';;; ';;;
1
';;; c
1
c c
1
c c
1

':; ':; ':; ':; ':; ." ."0 ."0 ."0 .:;
1 0.20

"'. ',.~ ',.... ....',. .;.... .. 0. 0. 0. 0. 0.

~ E
,.:.
E
.:. ~ If ~
':; :; :; ':;
',. ',."" ',. ',.
.... ',.
:;
',. ',. 0
1 1
0
~ il: 1;l u
or ~
a- V>
N .... 'Cl
-
00

""
,.., 0 'Cl
""
0

- .....- ,..,1
N
do.
"" " 0 .;, N
N
00
N

1 Figura 31 Medias sobre el conjunto de mal/as de los glaciares deI coseno dei émgulo de
incidencia deI sol (IoJ. Evolucion de los valores medios diarios

1 ..!!
1.0

-8 0,9

1
u

~ ~ 0.8
g ~ ,
..!! 3 0,7 ;
1 ~ ~ .................

1 , Ji
~

'"
~ ~ ~
u-
~ ,~
~ 0,6 :

05:
, f
--Charquini, 1 mai
....

'.
~

~u
u
~c 0.4 t '
- •• Charquini, 10 septembre

1
"
DO fi ' --longo. 1. mai
c ~0.3t,
"
:g." E
~ 02 j., .. ., ..... lango. 10 septembre
~ . Tous Ics glaciers, 1 mai "

1
.5 "

1 8 Tous les glaciers. 10 septembre' '.

: '_0 ... "1

8 9 10 Il 12 13 14 15 16 17

1
1 Heure de la journée
1 _ __

Figura 32 Medias sobre conjunto de mal/as de los glaciares deI coseno dei émgulo de
incidencia deI sol (IoJ. Evolucion durante los dias de sol maximo.{1 mayo) y minimo (lO sept.)
1 de la estacion seca

Resumiendo y coma conclusion de este parrafo, se puede decir que la evoluci6n de la

1 insolaci6n en el conjunto de los glaciares durante el afio no es concomitante a la de Zongo. En
consecuencia, la estacionalidad es un factor de influencia que deberia ser tomado en cuenta en

1 el desarrollo de un modelo para estimar los caudales hidricos de deshielo dei conjunto de
glaciares a partir de los deI glaciar Zongo. La variabilidad de la insolacion durante el dia, en
cambio, es mas comparable sobre todo al comienzo de la estacion seca. Pero la insolacion deI

1 conjunto de glaciares durante la tarde es con todo elevada comparada con la deI Zongo. Por
esto, los factores de conversion de energias en caudal son probablemente importantes en la
modelizaci6n, sobre todo al final de la estaci6n seca.
1
1
1 5.45. La insolaci6n media de las superficies bajo la isoterma cero

1 En un futuro modelo de deshielo deI conjunto de glaciares, podria ser interesante separar
superficies de glaciares en dos zonas abajo 0 arriba de la isoterma cero grado, 0 incluso tomar
en cuenta diferentes clases de altura coma en el modelo de RIGAUDIERE et al. (1995) (cf.

1 2.2.4). Para ello, se debe poder calcular la media de insolaci6n de las mallas, es decir la media
de los cosenos deI angulo de incidencia, para una superficie de glaciarres hasta una curva de
nivel dada. En este caso, la insolaci6n media no depende solamente deI cenit ni deI acimut deI
1 sol sino también de. la altura. El procedimiento es idéntico al que se presenta en el pârrafo
5.4.3. En este casa es absolutamente indispensable calcular las medias de insolaci6n para
gamas de valores de cenit, de acimut y, ademâs, de altura hasta el limite en el cual se desea
1 tratar las superficies glaciares. En efecto, un calculo triple en todas las mallas (cenit, acimut,
altura) de la cuenca tomaria mucha tiempo.

1 Se eligen los mismos valores de cenit y acimut que en el calculo de insolaci6n en las
superficies enteras de los glaciares (cf. 5.4.4). La gama de curvas de nivel esta comprendida
entre 4800 y 6050 m, con pasos de 25m, 10 que da 50 valores de altura. Se calcul6 la
1 insolaci6n media de las mallas bajo estas alturas, haciendo variar los 8 angulos de acimut y los
16 angulos de cenit para los glaciares Zongo y Charquini. Con ayuda de los cuadros
resultantes, de dimensi6n 50 x 8 x 16, las insolaciones medias de ambos glaciares sobre las
1 superficies, bajo una isoterma a una hora y fecha dadas, puede ser facilmente calculadas.

1 Para la comparaci6n de glaciares se han elegido tres periodos de observaci6n durante los
cuales coesxisten datos de caudales dei glaciar Zongo y deI canal Prado y que no contienen
dias de fuertes precipitaciones ni dias de operaci6n de valvulas dei Jago Mamankota (cf.

1 5.2.1). El primer periodo de observaci6n es del 1 de maya al 25 de julio de 1995; el segundo

dei 16 de agosto al 3 de septiembre dei mismo ano y el tercero dei IOde junio al 19 de julio de
1996 (Anexo C). Todos los periodos elegidos se encuentran en estaci6n seca porque en
1 estaci6n de lluvias los caudales deI canal Prado son muy influenciados por las precipitaciones.
Para la estaci6n seca 1997 no disponemos todavia de valores tratados de caudales en el
momento de la realizaci6n dei estudio.
1 La estimaci6n de la altura de la isoterma cero grade cada media hora durante los periodos de
observaciones fue realizada con ayuda de temperaturas medidas cada media hora en la
1 estaci6n Campbell 21 X en el glaciar Zongo a 5200 m de altura (cf. Figura 25, 5.4.4). A fin de
atenuar los errores de mediciones y la oscilaci6n de la altura de la isoterma cero grado, se ha
preferido establecer la curva de medias de temperaturas en 30 minutos, centradas y reducidas
1 para cada dia y sobre el periodo de las dos estaciones secas 1995 y 1996 (Anexo D). Esta
curva de la evoluci6n media de la temperatura permite recalcular las temperaturas en 30
minutos, a partir de medias diarias y de diferencias tipo de valors en 30 minutos durante el
1 dia. La altura de la isoterma cero grade fue calculada suponiendo un gradiente constante de -
0.671/100 m.

1 Con ayuda dei cuadro en tres dimensiones, descrito arriba, las insolaciones medias (por
unidad de superficie) de superficies bajo la isoterma cero grade estimada (medias sobre las
1 mallas bajo la isoterma cero grade de los cosenos de los ângulos de incidencia dei sol) fueron
reconstituidas al paso de tiempo de una media hora para los glaciares Charquini y Zongo.
Luego se calcularon las medias diarias que se presentan en la Figura 33.
1
1
1
La Figura 33 corresponde de alguna manera a la de las evoluciones medias de sol de
1 superficies enteras de los glaciares (cf. 5.4.4, Figura 31). Asimismo, se observa una
disminuci6n de las medias diarias dei Zongo en relaci6n al Charquini durante la estaci6n seca.
Esta disminuci6n es, no obstante, menos pronunciada que en cl caso de medias sobre
1 superficies enteras. Esto tiene su explicaci6n en cl hecho que cl peso de la parte superior dei
glaciar Zongo, que provoca un efecto de estacionalidad sobre la diferencia de insolaci6n de los
glaciares Charquini y Zongo, es reducido, ya que se encuentra- con menos frecuencia bajo la
1 isotenna cero grado que las partes inferiores dei glaciar. Los valores nulos de la Figura 33 son,
por otra parte, el resultado de dias muy frios, durante los cuales la isoterma cero grado se
encuentra bajo glaciares a 10 largo deI dia.
1
0,90 ; 1 mai - 25 juillet 1995 16.8.-3.10.95 10 ju in - 19 ju illet 1996 i 0,30

1 ~ ~
~ ~
::>
~

c
0,80 .
~
~.~!

,
1

0,25
-0 ~ 0,70 .

1 u
g
u

~5 Ë
e•
·U
] ~ 0,60
,

"
.
,
0,20
~
-0 U
~ ~ 0.50 , -""

1
, ~

là=:: ...... ::
, "
",
,~1 0,15 :~
=;; -8 ..

::1"
1 t7'
"'0
En
En
::1
0,40
'.
~.'
.
, .:
" ,''.
",
. l;l
..c
U
.~ ~
'

..
1 !
0.30 "
0,10
0-0 \
u ::>

1
ri 8
-0
'" 0,20
~~::> ....... Zongo
0,05

1
·U
oS .;
..
~
0,10
--Charquini
0,00 0,00
l-mai-95 16-l1lai·95 31·mai-95 15-juin-95 3O-jum-95 IS-juil-95 Zo-aoùt-95 lo-juin-96 25-juin-96 lo-juil-96

1 Figura 33 Insolaci6n de las superficies bajo la isoterma cero grado (integrales en mal/as
. bajo la isoterma de los cosenos de los angulos de incidencia dei sol) en los glaciares Zongo y
1 Charquini. Valores integrados en el dia (sin integraci6n media hora)

En la Figura 34 se presenta una comparaci6n entre la insolaci6n en toda la superficie de los

1 glaciares y en las superficies bajo la isoterma cero grado. Como la Figura 33, ésta muestra la
integral de insolaci6n en la superficie (la suma de cosenos deI angulo de incidencia deI sol en
cada pixel). En la figura se muestra igualemnte la evoluci6n de la superficie bajo la isoterma
1 cero grado. Como el valor mâximo de insolaci6n de una malla es 1 (cos (0°», el valor mâximo
de la integral en las superficies bajo la isoterma cero corresponde a esta superficie (en km 2), 10
mismo para el valor mâximo de la integral en toda la superficie dei glaciar. Evidentemente,
1 estos valores maximos jamâs son alcanzados, ya que esto significaria que todas las malIas de
la superficie considerada estarian en cada instante expuestas al sol.

1 Es interesante comparar ahora la insolaci6n con los caudales de deshielo de los glaciares
Zongo y Charquini. Este el objetivo deI pr6ximo capitulo.

1
1
1
1
g
c:=::::::J Inlégrak au dessus de 1"solhemlC ;a;ro
....16.8 -110,95 Glacier Zongo

1 ...
'0
ô
.S
'c
1.8
1,6 .
,
1.4 .
....
Intégrale uu dessous de l'isolhennc zào
....... Surlàee au dessous de l'isolhennc ;a;ro degré. Zongo

1 mai - 25 Juillel 1995

10 JUill' 19 juillet 1996

'"... L---;-....:....:~ .' .'..

1 '&>
<:
os
...'"
'0
1.2

1,0 ,.
0,8 ..
i":
1
,

'"<:
: ;)

1
';;; 0,6 .
o
u
'" 0,4 -:
"
'0

~
0.2 -
0,0

1 o
t/)
l-mai-95 l6-mai-95 31-lTlJi-95 15-juin-95 3Û'juin-95 15-juil-95 2Û'aoûl-95 lÛ'juin-96 25-juin-96 lÛ'juit-96

0,5 1 Glacier Charqulnl

1
'.
"
"

0,4
"
:'
"
. ..
'.,
"

,'.

- .,
"
.
','

.',:;; .... . . , .'

,
"
. .'
1
"
:.
' ,", "
. " , "
. : .. : •... :: "
"
0.3 " " ; ", ""

0,2

1 0,1 '

1 0,0 '
1-lmi-95 16-rroi-95 31-rroi-95 15·juin-95 3Û'juin-95 15-juil-95 20-aoùt-95 IO-ju in-96 25-ju in-96 lÛ'juil·96

-- '

1 Figura 34 Sol de las superficies bajo la isoterma cero grado (integrales en las mal/as
bajo la isoterma de los cosenos de los angulos de incidencia dei sol) en los glaciares Zongo y
Charquini. Valores entegrados en el dia (sin integraci6n media hora). La unidad es el km 2.

1 En punteado: Superficies bajo la isoterma cero dei glacial' Zongo

1 5.5. Simulaci6n de los caudales deI glaciar Charquini con dos modelos
simples y comparaci6n con la insolaci6n de los glaciares

1 5.5.1. Simulacion de los caudales de deshielo dei glaciar Charquini con dos modelos
simples

1 El modelo mas simple para simular los caudales de una superficie glaciar, a partir de caudales
conocidos de otra superficie, consiste en admitir la hip6tesis que los caudales espécificos
llevados a las supeficies totales de los glaciares son iguales. Este modelo lineal simple se
1 escribe: -

1 donde:
Qch Caudal diario dei glaciar Charquini (simulado)

1 Qz
Qbase
Caudal diario dei glaciar Zongo
Caudal adicional debido a las diferencias entre los glaciares
SIOI
Superficie total dei glaciar Charquini
1 d,

Superficie total dei glaciar Zongo

1
1
Este modela se basa en las hip6tesis hechas en el capitulo 0 postulando que la producci6n de
1 agua de deshielo en la superficie de las malIas de los glaciares influye en la diferencia de los
caudales. Pero este modelo, extremadamente simple. supone ademas que la producci6n de
agua de deshielo es en promedio igual en toda la superficie de un glaciar y esta ligada de una
1 manera lineal a las producciones de los otros glaciares. A fin oe tomar en cuenta diferencias
de producci6n entre los glaciares, se introduce un caudal de base constante. Un modelo un
poco mas sofisticado adapta la uItima hip6tesis deI modela lineal simple planteando que las
1 malIas de la superficie de ablaci6n son las que producen el caudal. Este modelo se basa en la
idea que arriba de la isoterma cero grado, no hay deshielo y toda la energia es utilizada para la
sublimaci6n. Este segundo modelo puede escribirse:
1 s':w
Qch = SiS" Q
ch
Z + Qbase (11 )
1 donde:
Z

1 st'" = Superficie bajo la isoterma cero grado deI glaciar (i) (media deI dia
de las superficies bajo la isoterma cada media hora)

1 Para adaptar mejor los términos deI caudal de base de estos dos modelos, se utiliza la funci6n
criterio de Nash. El resultado de las simulaciones con los dos modelos se presenta en la Figura

1 35. Se realiza una aplicaci6n de cada modelo sobre el conjunto de periodos de observaci6n en
estaci6n seca 1995 y, separadamente, sobre el periodo de observaci6n en estaci6n seca 1996,
que es mas corto. Evidentemente, estos modelos simples no estan adaptados a la descripci6n

1 de procesos complejos en su origen de las diferencias entre los caudales de los glaciares, pero
estas simulaciones permiten la interpretaci6n de diferencias topograficas entre ambos
glaciares, Charquini y Zongo, detalladas en el capitulo anterior.

1
1
100

90

80
1

r
1
1

1
Imai-25juillet 1995

il mai 1995-3scpt 1995:

1modèle 1: Qbase = 25 Vs;fonction de critère: 0.47
16août-3sept. 1995

- 10 juin· 19 ju iIIet 1996

imodèle 1:
iQbase = 18 Vs, f.c. =0,86
-~

imodèle 2: Qbase = 19 Vs;fonction de critère: 0,32 imodèle 2'

1 ~
70

60
1
1
iQbase = 14 Vs, f.c. =0,84

1
.l!l 50
15
~
40

1 30

20
--débits observès
10 --débits simulés avec les surfaces tolales (modèle 1)

1 o'
l-mai-95
....... débits simulés avec les surfaces au dessous de l'isothennc zéro (modèle 2)
1

l6-mai-95
1

2-ju in-95
1

l7·ju in-95
l

2-ju il-95
,

17-ju il-95
•

22-août-95
1

12-ju in-96
i·

28-ju În-96
--- i

14-ju il-96
j

1 Figura 35 Simulaci6n de los caudales dei glaciar Zongo con modelos que se basan en los
caudales especificos relacionados con las superficies totales de los glaciares (modelo 1) y en
.. _--~-_.- -_ .. --~-----

1 relaci6n a las superficies bajo la isoterma cero grado (modelo 2). Apuntalamiento dei
parametro dei caudal de base con el criterio de Nash.

1
1 5.5.2. Interpretaci6n de los resultados de la modelizaci6n comparando las diferentes
insolaciones
1 Se observa en las simulaciones que los dos modelos utilizados reconstituyen bastante bien los
caudales deI glaciar Charquini en estaci6n seca 1996, mientras que en estaci6n seca 1995, en
un periodo mas largo, la adaptaci6n de los caudales simulados es menos buena (cf. Figura 35).
1 Evidentemente, se debe poner en relaci6n con la diferencia de la duraci6n de los periodos de
observaci6n. Observemos en la Figura 35 el resultado de la simulaci6n en los dos primeros

1 periodos. Los caudales al comienzo de la estaci6n seca frecuentemente son sobrestimados en

reIaci6n a los caudales reales, mientras que hacia el final de la estaci6n se encuentran
generalmente subestimados. Esto parece estar vinculado a las diferentes evoluciones de

1 insolaci6n durante la estaci6n en los glaciares Charquini y Zongo.

Si se retoman la Figura 33 y la 34, se observa una disminuci6n mas fuerte de la insolaci6n en

1 el glaciar Zongo que en el glaciar Charquini durante el avance deI invierno (cf. 5.4.4 y 0).
Estas evoluciones diferentes conducen muy probablemente a una disminuci6n mas importante
de los caudales deI Zongo en relaci6n a los caudales dei Charquini durante la estaci6n seca.
1 Ahora bien, en el casa deI modelo de superficies totales, los caudales simulados dei Charquini
estân vinculados, por un factor constante en el tiempo, a los caudales dei glaciar Zongo. POl'
consecuencia, cuando los caudales dei glaciar Zongo disminuyen mas fuertemente que los
1 caudales reales deI glaciar Charquini durante la estaci6n seca, los caudales simulados, que al
principio son sobrestimados, poco a poco son subestimados. Puede que la subestimaci6n de
los caudales deI Charquini entre el 16 de agosto y eI 3 de septiembre se debe a la incapacidad
1 deI modelo para simular una no linealidad de los caudales que no depende de la
estacionalidad, ya que los caudales durante este periodo de observaci6n son excepcionalmente
elevados. No obstante, la disminuci6n de los caudales estimados en relaci6n a los caudales
1 reales se observa también si s610 se mira el primer periodo de observaci6n en 1995 y durante
el periodo de observaci6n en estaci6n seca 1996. Durante estos dos periodos de observaci6n,

1 las subestimaciones no estân vinculadas a caudales excepcionalmente elevados. Ademas, las

subestimaciones de los caudales reales por los dos modelos corresponden a dias durante los
cuales la insolaci6n media de las superficies, bajo la isoterma cero grado dei glaciar Zongo,

1 era fuertemente inferior a la deI glaciar Charquini.

En el casa dei modelo de las superficies bajo la isoterma cero grado, los caudales simulados
1 no dependen de un factor constante en el tiempo de los caudales deI glaciar Zongo. Se podria
incluso pensar que este factor, el cociente de las superficies dc los glaciares bajo la isoterma
cero grado, introduce cierta estacionalidad en los caudales simulados, ya que las temperaturas
1 disminuyen al avanzar el invierno y, con ellas, las superficies..bajo el isoterma cero grado.
Recordando las curvas hipsométricas de los glaciares (cf. 5.3.2. Figura 27), se observa que la
dei Charquini es mucha mas plana que la dei Zongo. En consecuencia, un descenso de
1 temperatura conduce a una disminuci6n mas grande de las superficies bajo la isoterma cero
grado dei glaciar Charquini que deI glaciar Zongo. Como estas superficies son generalmente
mas pequeiias en el casa dei Charquini que en el casa dei Zongo, su disminuci6n relativa es
1 mas grande que la dei Zongo. y por 10 tanto los cocientes deberian tender mas bien a
aumentar. Esto introduciria una estacionalidad en el otro sentido. No obstante, se observa un

1 cociente de las superficies bajo la isoterma cero grado que no parece tener tendencia
estacional (Figura 36).

1
1
1
1 mai - 25 juillet 1995 16.8. - 3.9.95 10 juin - 19 juillet 1996 1
0,50·

1 .~
o
50
0,45
i
T
!.....-- ---
1

0,401
.~
1 -8
e
'1:! 0
0,35,
!

... ~ 0,30
~
1 05
Q.5
~~
,
0
~ 025
50
-;; ~ 020 ;
l'
1

"0 u ' ,

1 ~~
"0
~
0,15

0,10

1
<Il
~
tS 0,05
~
CI)

0,00 1 i ·1

1 Figura 36
l·mai-95 l6-mai-95 2-juÎn-95 17-juin-95 2-juil-95 17-juil-95 22-août-95 12-juin·% 28-juin-%

Cociente de las superficies de ablacion bajo la isoterma cero grado entre el

14-juil·96

glaciar Charquini y el glaciar Zongo

1 Un modelo a desarrollar podria basarse en la ecuaClOn (9) (cf. 5.1.4). El factor de la

1 multiplicaci6n corresponderia al cociente de esta ecuaci6n para el glaciar Charquini y para el

glaciar Zongo. El calculo de este cociente hace necesario conocer la parte de la energia total
que es originada por la radiaci6n solar directa en cada paso de tiempo de la integraciém (Ri).

1 Se podria escribir el modelo parapasos de tiempo de esta manera:

Jl dS + (1- Ri)' Sr har

(i»
1 2
Fi' (Ri
sChar
a '

Q Char = " f , iJ
LJ Zongo'
Q
Zongo
(12)
(i»
1 i=l Fi ·(Ri la ·dS+(I-Ri)·Sf
sZongo
f,i

1 de 10 cual :

1 Qchar
Qzongo
=
=
Caudal simulado deI glaciar Charquini
Caudal medido 0 simulado con un modelo determinista dei glaciar
Zongo

1 =
=
Indicador de la mafiana (= 1) Yde la tarde (=2)
Factor de la conversi6n de la energia disponible de agua de deshielo,
independiente deI glaciar
1 = Parte de la energia total recibida originada por la radiaci6n solar
directa, independiente deI glaciar
Superficie media bajo el isoterma cero grado deI glaciar Charquini
1 =

durante un medio dia

= Superficie media bajo el isoterma cero grado deI glaciar Zongo

1 durante un medio dia

1
1
1 Se ve que, contrariamente a la ecuacion (9), el término de la energia total desapareci6 en la

1 ecuacion de este modelo. Como la energia total disponible por unidad de superficie no
depende de la topogratla de la malla (se trata de la energia total hipotética disponible para una
malla siempre expuesta al sol, ElOI = E prop + R ), los términos de superficies de glaciares bajo

1 la isoterma cero son los que toman en consideracion la diferencia de esta energia total entre
los glaciares. Estas superficies, asi como las integrales de los cosenos de los angulos de
incidencia deI sol, pueden ser detemlinados con los métodos desarrollados en este capitulo a

1 un paso de tiempo de la integraci6n en el transcurso de un dia cualquiera.

Para verificar ,el funcionamiento de esta proposicion para el modelo habrâ que encontrar un

1 método para determinar la parte de la radiacion solar directa de la energia total en cada paso
de tiempo de la integracion en cl dia (Ri).

1
1
1
1
1
1
1
1
1
1
1
1
1
1 6. Conclusiones
1 El objetivo principal de este trabajo era encontrar un método para reconstituir los caudales
naturales deI rio Zongo, que no sean influenciados por el sistema hidniulico dedicado a la
producci6n de electricidad y que no contengan los aportes dei deshielo de los glaciares. El

1 paso para alcanzar esta meta sigui6 dos vias. La primera era encontrar un método de
reconstituci6n de los caudales naturales deI conjunto de escurrimientos de las vertientes y deI
deshielo de los glaciares. La segunda consistia en desarrollar el esbozo de una modelizaci6n

1 para el deshielo de los glaciares.

Se ha estudiado el funcionamiento deI sistema hidraulico y realizado un analisis cualitativo y

1 cuantitativo de los componentes que influyen en los caudales observados a la salida de la
cuenca de drenaje. Los resultados de este analisis mostraron que hay que tomar en cuenta
principalmente los cambios de stock de agua en las represas y los aportes externos de la
1 cuenca, para la reconstituci6n de los caudales naturales; la influencia de llenado de las
represas en estaci6n de lluvias es probablemente insignificante. Cuantificando la influencia
dei cambio de stock en diferentes represas, se pudieron distinguir represas cuyos volumenes
1 utiles son suficientemente grandes coma para ser tomados en consideraci6n en la
reconstituci6n de caudales. A partir de este analisis preliminar se desarro1l6 un método para la
reconstituci6n de los caudales naturales al cierre de Sainani en estaci6n seca. Se desarro1l6
1 también un método para la reconstituci6n en estaci6n de lluvias. Estos métodos utilizan las
mediciones diarias de nivel en las represas y caudales diarios, en ciertas plantas
hidroeIéctricas, para estimar los caudales que provienen dei desalmacenamiento de represas y
1 de aportes externos.

1 Una cstimaci6n de los rendimientos de las plantas hidroeléctricas Zongo, Tiquimnai, Santa
Rosa II y Sainani, permiti6 reconstituir sus caudales turbinados al paso de tiempo horario
durante la estaci6n seca 1997, Yestudiar la influencia de la gestion de las pequeiias cuencas de

1 regulaci6h en el rio principal que no habian sido tomadas en cuenta en el estudio preliminar.
Para analizar esta influencia, hubo que desarrollar un método de estimacion de la curva de
llenado deI conjunto de estas cuencas intermedias. El método desarrollado utiliza caudales de
1 las tres plantas hidroeIéctricas, Zongo, Tiquimani y Santa Rosa II, asi coma los caudales
totales en Sainani para reconstituir la curva de lIenado. No disponiendo de los caudales totales
en Sainani, se tuvo que estimarlos a partir de los caudales turbinados en la planta
1 hidroeléctrica.

Los caudales naturales provenientes deI deshielo de los glaciares y deI escurrimiento de las
1 vertientes fueron reconstituidos en 5 meses, de maya a septiembre 1997, con el método
desarrollado para la estaci6n seca. Finalmente, se utiliz6 una comparaci6n entre la evoluci6n
dei hidrograma diario reconstituido y la evoluci6n estimada deI deshielo de una pequeiia parte
1 de las superficies glaciares de la cuenca de drenaje para discutir la calidad deI hidrograma de
los caudales naturales obtenidos.

1 La calidad de los caudales naturales reconstituidos es buena. El método desarrollado para

calcular los caudales naturales en estaci6n seca toma en consideraci6n influencias principales

1 deI sistema hidraulico de la producci6n de electricidad sobre los caudales medidos en Sainani.
Podria ser mejorado tomando en cuenta la variaci6n deI stock en las cuencas intermedias. Esto
podria hacerse utilizando el método desarrollado para estimar la curva de llenado deI conjunto
1 de estas cuencas, 0 estimando los cambios de stock en las dos cuencas intermedias mas
importantes, con ayuda de mediciones de nivel. La variacion deI stock de estas cuencas es

1
1 importante sobre· todo al comienzo y al final de las crecidas. Pero s610 durante la puesta a
punta de un modelo hidrol6gico se sabra si hay que mejorar la calidad de los caudales
1 reconstituidos por medio de un refinamiento deI método de la reconstituci6n de los caudales
naturales.

1 En ausencia de datos de los caudales totales en Sainani, no era posible reconstituir los
caudales naturales en estaci6n de lluvias utilizando el método adaptado a esta estaci6n. En

1 consecuencia, la calidad deI método no pudo ser probada. Este trabajo se hara cuando se
disponga de series de datos de los caudales totales en estaci6n de l1uvias, probablemente en
1998.

1 La modelizaci6n deI deshielo deI conjunto de superficies glaciares en el valle es primordial.

Este modelo permitira reconstituir los caudales naturales en Sainani sin aportes de los
1 glaciares y también perfeccionar un modelo hidrol6gico. El esbozo de este modelo fue
desarrollado durante este trabajo. Se basa en una estimaci6n deI deshielo deI conjunto de
superficies glaciares a partir de los caudales deI emisario deI glaciar Zongo, tomando en
1 cuenta las diferencias topograficas entre los glaciares. Una consideraci6n te6rica de la
influencia de la topografia sobre la producci6n de agua de deshielo en la superficie de los
glaciares, ha permitido conc1uir que esta influencia no es lineal y que un enfoque con factores
1 correctivos probablemente no funciona. Haciendo algunas hip6tesis, para simplificar, sobre la
dependencia de los términos deI balance de energia frente a la topografia, se puede reducir
esta influencia a un factor que mide la energia de la radiaci6n solar directa recibida sobre la
1 superficie bajo la isoterma cero grado durante el dia. Esta medida es la doble integral sobre la
superficie bajo la isoterma cero grado, cuya posici6n varia durante cl dia, y en horas del dia,

1 del coseno deI angulo de incidencia deI sol.

Con ayuda de la superposici6n de una imagen satelital y de un modelo numérico de terreno

1 MNT) de la cuenca de drenaje, se obtuvo los MNT de los glaciares Zongo y Charquini, y deI
conjunto de superficies glaciares. A partir de estos MNT, y de la determinaci6n de la posici6n
deI sol durante el dia y el afio, se pudo calcular tanto la energia de la radiaci6n solar directa

1 . recibida en la totalidad de la superficie de los glaciares coma la recibida en la superficie bajo

la isoterma cero grado al paso de tiempo de una media hora. Se calcularon las integrales deI
dia de estas energia para los glaciares Charquini y Zongo. Estas integrales muestran una fuerte
1 estacionalidad, que es diferente en los dos glaciares. De tres periodos de observaci6n en
estaci6n seca 1995 y 1996, se modelizaron finalmente los caudales diarios deI glaciar
Charquini a partir de caudales deI glaciar Zongo, con ayuda de dos modelos simples que no
1 toman en cuenta la estacionalidad.

La comparaci6n entre la diferencia de los caudales simulados y observados y las integrales de

1 la energia solar directa recibida, conduce a la conc1usion de que la estacionalidad es un factor
que no puede despreciarse en la simulaci6n dei deshielo dei conjunto de superficies glaciares a
partir de los caudales deI glaciar Zongo. Se pudo mostrar que las diferencias entre las
1 evoluciones de la energia solar directa recibida por los glaciares comparados, son las que
determinan esta estacionalidad. Pero consideraciones teoricas mostraron también que la

1 influencia de la diferencia de estas energias recibidas por los glaciares, agregada a la

diferencia de sus caudales, depende de la relacion entre la energia solar directa mâxima (en un
ângulo de incidencia nulo) y la energia total recibida par los glaciares. Segtm las hipotesis
1 admitidas para el esbozo deI modelo a desarrollar, estas enfoques son independientes de la
topografia y, por consecuenciam iguales por unidad de superficie de los glaciares.

1
1
1 Parece final mente posible desarrollar un modelo que se basa en la energia de la radiaci6n solar

1 directa recibida por la superficie bajo la isoterma cero grado deJ gJaciar de referencia y de la
superficie glaciar deI cual se hace una simulaci6n deI deshielo. Pero es nccesario incluir en
este modelo la evoluci6n de la relaci6n entre la energia soJar directa maxima (en un ângulo de

1 incidencia nuJo) y la energia total. La evoluci6n de esta relaci6n puede ser probablemente
calculada a partir de mediciones climaticas en el glaciar Zongo.

1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1 7. Anexos

1 Anexo A Estimaci6n de las relaciones entre el nivel de agua y el volumen para las
represas Mamankota, Livinosa, Alto Viscanchani y San Pedro, con la hip6tesis que sus
formas son asimiladas a casquetes esféricos
1 Valeurs mésurées Valeurs calculées
Volumen Superficie Nivel Radius 1*) Superficie Volumen
1 maximo
(V max )
max.
(planimétria)
maximal
(h max ) (r)
calculada calculado
(Smax) (Vmax)
(Smax)

1 Hatilata
1000 m
2483,33
m..!
j
m m m..! j
1000 m

Taipicota 3000,00

1 Guaraguaran i
Alto Viscanchani ?
432,73
40000 8,83 725 40000 19462
Zongo 3256,59
1 Marnankota
Livinosa
236.00
523,00
67010
70352
7,19
7,72
1456
2796
65593
135430
236,00
523,00
Sankayuni 433,37
1 San Pedro ? 15000 2,7 886 15000 8869

([Link] + hmax /3)

1 1*) r =
7r' hmax
2 ou

Relaci6n nivel/volumen utilizada:

1
1
1
1 represa al nivel
maximal
1
1 1
1
1
1
1
1 Anexo B Desalmacenamjenlo de los lagos de regulacion en porcenlaje dei caudal
turbinado, en estadon seca 1996

1
25,0%

1 20,0010
15,00/0
.
1- -
L
. ·············_··1
i--Zongo
Mamankota
~ _ _~ .._ .

1
~

:ô 10,00/0
'd.I
"Cl
=s 5,00/0
"Cl
tt. 0,0%
1 -5,0%
-10,00/0

1 -15,00/0
'-0
~
'-0
~
1
'-0
~
1
'-0
~
1
'-0
1
'-0
~
'-0 '-0 '-0 '-0 '-0
i
'-0
1

.,... .,... ~ ~ ~ ~ ~
,.,...
~
0\
~
......
M
......
00
'-0
......
N
~
l"
l"
......
N
l"
......
l"
-
00
...... 00
......
'-0
-
00
...... .,... ~
0
0

1
M M

date

1 40,0010 . . -. . ·-1
- - - HatiJata

1 35,0%
30,0010
25,0010
- - - Taipicota
: Gua~~u.~~ani 1
1

] 20,0%
1 :g
.g
15,00/0
10,0%
tt. 5,00/0
·+------4b......J.-."':'-'\'·~II~'.~~....-:-:~~Io.A.N/
1 0,0010
-5,0010
-10,0010
"
"l
..........

,
1 -15,0010
'-0
~
~
M
'-0

~
N
'-0
~
~
l"
'-0
~
l"
......
N
\0
~
~
l"
\0
~
00
:::
\0
~
--
00
\0
'-0
~
00
......
\0

~
.,... 0
\0

~
\0

,.,...
~
0
M M

1 date

1
1
1
1
1
1
1 25,0%
--Alto Viscanchani!
20,0% ,
1 15,0%
--Livinosa

-
:0 10,0010

1
'Q)

":::l 5,00/0
"'$. 0,00/0

1 -5,0%
-10,00/0
-15,0010 1 1 1 1 1 t"

1 \0

----
0-
V)

.....
\0

----
0-
V)

00
\0

--
0-
::e
N
\0

--......
0-
::e
\0

----
0-
......
N
\0

--
0-
......
~
\0

----
0-
00
\0

----
0-
00
\0
\0

----
0-
00

.....
\0

----
0-
0-
V)
\0

--
0-
0-
25
.....
\0

--
0-
0

--
V)

1 datc

1 10,00/0 -Sankayuni
--San Pedro
5,0% i

1 :0
''1)
0,0%
"
1
:::l

"'$. -5,0010

-10,00/0

1 -15,00/0 t ""-+-"-"---i-
\0 \0 \0 \0 \0 \0 \0 \0 \0

1 ----
0-
V)

..... ----
0-
......
N
----
0-
......
......
----
0-
00 ----
0-
00
\0
----
0-
00 ----
0-
0-
V)
0-
0;
--
0
0-
25
~
..... ..... V)

date

1
1
1
1
1
1
1
1
1 Anexo C Delerm;nac;on de los per;odos de observacion de los caudales dei glaciar
Zongo y dei canal Prado

1 1· --~-l

J' ,.1
Il 0

~.,.
Ti 1

...... i
1 ..
,',',
'
,

. ,
'-. :: .
:;11: \._ ..
, :1,/ :;",
~ ~.,'
2

4
0-

-2·
i1
1
1

........ f
( 16août à
\. Û
0

1 .... 3 sepl 95 6 -4.

~~
~~

'. 1·
'"
2- 200
~
~I
,, -6. ,g E
....
:0
Irrui à 25 juillet 95
8 !9 ...
1i.. .2

1
Cl 150 1 '- f!
1 10 .~ '8-
Q" E
100 1 ~
12

1 50
1

14

1 o i
l-rrui-95 l6-rrui-95 JI-rrui-95 15-ju in-95 JO-juin-95 15-juil-95 2D-aoûl-95 4-sept-95 2D-sept-95
16

- - Débits de l'émissaire du glacier Zongo ---DébilS du canal Prado

1 _ Précipitation à la Platafonm ....... Tel11'ératurc moyenne à 5200 m (glacier Zongo)

1 J50

1
-i 0

-~
!
1
(0

JOO
." "
1 1

2 0-

1 250 '..
"
". ",
,
:'
,

;
:
: .. 4 tO~
Û
.'
'. E C

1
"
: ~5
~ 200
:
" '"
C 0>.
,~ .
\' 6 :j. E
.&>
0 150
"
.. .
"
'-
os. ~
'ü ~
• i
1
, '
" : .~
10juin à 19juillel 1996 "
8 .%oE
100 ~

1 10

- -1 " 1__ [ 12

1 lo-juin-96 2O-juin-96

- - Débils de l'émissaire du glacier Zongo

JO-juin-96 lo-juil-96 2D-juil-96 JO-juil-96

---Débits du canal Prado

9-aoùl-96 19-aoùt-96 29-aoûl-96

_ Précipilalion à la Platafoml3

1
....... Tel11'éralure lTlJyenne à 5200 m(glacier Zongo)

1
1
1
1
1 AnexoD Evolucion media durante el dia de las temperaturas centradas reducidas

1 ---
u
2,0
.... ,

B
0
'-'
1,5
Station Campbell 21 X .: Séries température des

1 .:;
..
, "1::1
'v
sur le gbicier Zongo,
à 5200rn
:
.. "
saison sèche 1995 et 1996

..
v
'4)

ë
1,0
. .. ,
.
1 v
u
li>
~
0,5
.
-...
..
::s
1;;
0,0
..
".
1 'v
0-
E
B
li>
v
-0,5
1 0
. '.
-. ..
. '. ...... ,
"1::1

1 V
C
c -1,0
v
>.
0
.... " . - ....
.... - . ..
: ...... - .. oo

"

1 -1,5
0:00 2:00 4:00 6:00 8:00 10:00 12:00
Heure de la journée
14:00 16:00 18:00 20:00 22:00 0:00

1
1
1
1
1
1
1
1
1
1
1
1
1
8. Référencias
1 ENDE, (1992). Hidrologia deI valle deI Zongo, Empresa National de Electricidad SA

1 EROUTA., (1996). Rapport de stage, version provisoire

ALEGRE W T, (1996). Influencia de los glaciares en el comportamiento hidrolologico de

1 cuencas de alta montafia, estudio de casos en Perù y Bolivia, Thèse d'ingénieur à
l'université nationale d'Ancash, faculté d'ingénierie civile, Pérou

1 BOUVIER, C. ET HINGRAY, B. (1995). Une méthode hydro-géographique pour la

caractérisation du risque pluvial en milieu urbain Dans: Hydrologie urbaine et
modélisation spatialisée. Mémoire de DEA de B. Hingray, Université de Montpellier 2
1 RIGAUDJERE, P., RIBSTEIN, p., FRA NCOU, B., POUYAUD, B. ETSARAVIA, R. (1995). Un modèle
hydrologique pour le glacier Zongo, Informe No. 44, ORSTOM-Bo/ivie
1 RIGAUDJERE, P., RIBSTEIN, P., FRANCOU, B., POUYAUD, B. ETSARAVIA, R. (1995). Mesures
météorologiques, hydrologiques et glaciologiques sur le glacier Zongo, année
1 hydrologiques 1993-94, Informe No. 43, ORSTOM-Bo/ivie

WAGNON, P., RIBSTEIN, P., FRA NCOU, B., POUYAUD, B. ET VALDIVIESO, H (1995). Mesures
1 météorologiques, hydrologiques et glaciologiques sur le glacier Zongo, année
hydrologiques 1994-95, Informe No. 49, ORSTOM-Bo/ivie

1 BERTON, P., WAGNON, P., FRA NCOU, B., POUYAUD. B. ET VALDJVIESO, H (1997). Mesures
météorologiques, hydrologiques et glaciologiques sur le glacier Zongo, année
hydrologiques 1995-96, Informe No. 55, ORSTOM-Bo/ivie
1 FRANCOU, B. ET RIBSTEIN, P. (1995). Glaciers et évolution climatiques dans les Andes
boliviennes, glacier de Zongo et glacier de Chacaltaya, Cordillère Royale, 16° S, Bull. Insl.
1 fr. études andines 24(1), 23-36

LIBOUTRY, L. (1964). Traité de glaciologie,p.364, Masson et Cie, Paris 1964

1 ROBINSON, N (1966). Solar Radiation, Elsevier Publishing Company, New York 1966

1
1
1
1
1