APROVECHAMIENTO DE LOS RECURSOS HIDRICOS

CHIMENEA DE EQUILIBRIO

Estabilidad en los sistemas de regulación de las turbinas

1 INTRODUCCION

El funcionamiento correcto del sistema de regulación de las turbinas requiere que estas puedan
seguir rápidamente las variaciones de potencia de la demanda. El propósito de lo que se expondrá
a continuación es de disponer de algunas herramientas que permitan tomar decisiones a la hora
de desarrollar el proyecto de las obras. Se analizará el efecto de los estados impermanentes en
los escurrimientos en tuberías provocados por modificaciones en la demanda de energía, es decir
durante aquellos períodos de carácter transitorio durante los cuales el sistema de generación pasa
de un estado de equilibrio a uno nuevo y diferente.

Todos los grupos de generación eléctrica se encuentran condicionados por la demanda. Un grupo
turbina-generador solo puede generar la potencia que el sistema de consumo le requiera en cada
instante. De existir una diferencia, es decir una potencia desequilibrada, esta acelerará o frenara
el grupo generador. Esto es consecuencia de que entre el centro de generación y la demanda no
existe ningún elemento que pueda actuar sobre las diferencias que se producen. Los cambios en
la velocidad de rotación de los grupos generadores así provocados son inadmisibles salvo dentro
de límites más o menos estrechos.

Por lo tanto, cualquier cambio de la potencia exigida o demandada por la red requerirá de una
pronta adecuación de la potencia producida por el grupo turbogenerador. La modificación de la
demanda puede ser pequeña o grande. La mayor variación imaginable se presentaría en el caso
de rotura de la línea de transmisión caso en el que la potencia demandada pasa a cero en forma
instantánea.

En términos generales cuando aparece un desequilibrio entre la potencia generada y la tomada

por la red la diferencia se manifiesta como una modificación en la velocidad de rotación del grupo
generador. Esta alteración se utiliza generalmente como elemento básico para controlar en cada
instante los requerimientos de potencia a generar.

Cuando la potencia generada es diferente a la absorbida por el sistema las fuerzas

desequilibradas aceleran o frenan el grupo generador y el sistema de regulación de la turbina se
encarga de reducir o aumentar el caudal de ingreso y así mantener la velocidad de giro
controlando la frecuencia de generación.

Los estados inicial y final de este proceso serían:

.........................................................Estado 1
N1 = η γ Q1 H1
N2 = η γ ( Q1 ± ∆Q ) H1 .........................................................Estado 2

Donde para ambos estados se mantiene el salto H1 ya que el embalse, en el pequeño lapso
dentro del cual se producen estas variaciones, se mantiene constante.

Planteado así el problema seria de solución relativamente sencilla sino fuera por que toda
modificación del caudal, y por consiguiente de la velocidad en la conducción, origina
paralelamente variaciones en la presión debidas al fenómeno del golpe de ariete provocado por
las aceleraciones del fluido en la transición de los estados 1 y 2.

Lo que sigue no se debe tomar como un planteo riguroso del problema pues tiene como propósito
principal tener una comprensión cualitativa del fenómeno que estamos tratando. Considerando los

CHIMENEA DE EQUILIBRIO Pág. 1

 APROVECHAMIENTO DE LOS RECURSOS HIDRICOS

estados 1 y 2 como los iniciales y finales, entre ambos se ubican los estados transitorios que
pueden expresarse como:

N 2* = η γ ( Q1 ± ∆Q ) ( H1 m ∆H ) ..................................Estado
transitorio

Correspondiendo valores negativos de DH para valores positivos de DQ y viceversa. Recordando

la expresión de Michaud para cierre brusco, es decir cuando: TC<2L/a :

L ∆V
∆H = 2
g ∆T
y poniendo en la expresión del estado transitorio la variación del caudal en función de la velocidad,
resulta:
2 L ∆V
N 2* = η γ ( Q1 ± Ω ∆V ) ( H1 m )
g ∆t
Resulta claro que se trata de lograr que la variación, (deseada), del caudal:
( Q1 ± Ω ∆V )
prime sobre la variación, (no deseada), del salto:
2 L ∆V
( H1 m )
g ∆t
Es necesario que se cumpla N 2* < N1 ante una disminución de la demanda o N 2* > N1 cuando
esta aumenta. De no ser así el regulador de la turbina seguirá cerrando los álabes en el primer
caso o abriendo en el segundo.

Resulta claro que frente a conducciones cortas, L pequeños, y valores elevados de Dt, la
variación del salto H1 no será suficiente como para alterar la operación del regulador. También se
puede advertir que frente a iguales valores de DV, Dt y L la situación es diferente según resulten
entre sí las magnitudes de Q1 y H1.

En general no es de esperar una

relación lineal entre las variaciones de la
Vmáx velocidad en la conducción y los
tiempos de cierre, pero de todas formas
V la magnitud de los valores Dt están
fuertemente ligados a los tiempos de
cierre TC .

El problema ahora se plantea en

términos cuantitativos. Si la distancia
desde la obra de toma hasta la central
t Tc es L, ¿qué valor debe tener el TC ?, y
complementariamente: ¿qué otros
elementos pueden restringir el valor de TC ?

Si el TC es elevado se trata de un regulador que actúa lentamente frente a un cambio de la

potencia. Esto implica que se mantendrán potencias desequilibradas durante períodos más o
menos importantes que darán como resultado la imposibilidad de mantener acotados los valores
de la velocidad angular de las masas rotantes y por ende de la frecuencia.

La influencia sobre la velocidad de rotación de estas potencias desequilibradas dependerá de la

inercia de las masas rotantes. Si estas son muy elevadas las variaciones de la velocidad angular

CHIMENEA DE EQUILIBRIO Pág. 2

 APROVECHAMIENTO DE LOS RECURSOS HIDRICOS

serán relativamente poco importantes. De aquí que la inercia de las masas rotantes resulte un
2 L ∆V
elemento equilibrante. Por el contrario las variaciones a que estará sujeto el salto, , frente
g ∆t
a una variación de la velocidad dan una idea de la magnitud del elemento desequilibrante del
sistema.

Se tienen así los dos elementos básicos, equilibrante la inercia de las masas rotantes,
desequilibrante la sobrepresión generada en cambios de velocidad en la conducción, que
interactuan y permiten establecer si un sistema de regulación que opera sobre los grupos de
generación alimentados por una conducción resulta estable o inestable. Debemos aclarar que
interpretamos como estable aquel regulador que logra, luego de transitar un corto período dentro
de un régimen impermanente, converger en la nueva potencia demandada N2.

Para cuantificar estos parámetros definiremos los siguientes elementos:

Tiempo característico de las masas rotantes, TM : Las masas rotantes del sistema de
generación se integran con el rotor de la turbina, el eje que vincula el rotor con el generador y este
último. De la inercia total el generador provee entre el 85 y el 90% de total. Se define como tiempo
característico de las masas rotantes al tiempo que requiere una potencia desequilibrada P, igual a
la potencia instalada, en acelerar el sistema rotante desde una velocidad 0 hasta la nominal n. El
cálculo de TM es relativamente complejo. Una expresión empírica aplicable para niveles de
anteproyecto es la siguiente:

TM = 0.242 P 1 / 4 n 1 / 8

Donde P es la potencia expresada en kW y n la velocidad de rotación en rpm, quedando el valor

de TM en segundos. Resulta claro advertir que la inercia de los grupos es directamente
proporcional al valor de TM.

Tiempo característico de la tubería, TW : Es el equivalente al TM pero aplicado a la masa de

agua dentro de la tubería. Se define como el tiempo que necesita una carga H en acelerar a la
masa de agua contenida dentro de la tubería de longitud L desde una velocidad V=0 hasta V=V :

LV
TW =
gH

Novelli establece, a partir de experiencias propias, que para lograr una buena regulación es
necesario que:

∆H
< 0.40
H

Donde DH es la sobrepresión por golpe de ariete para un cierre total del flujo en un tiempo igual a
TC, que es el tiempo que demora el regulador en operar un cierre completo. Aplicando Michaud se
tendrá:
2 LV ∆H 2 LV
∆H = → = < 0.40 ........................................................1
g TC H g TC H

Considerando que los valores de TC resultan en la mayoría de los casos similares a los
correspondientes TM se puede escribir:

CHIMENEA DE EQUILIBRIO Pág. 3

 APROVECHAMIENTO DE LOS RECURSOS HIDRICOS

∆H 2 LV LV 1 TW
= < 0.40 → < 0.20 → < 0.20 ................2
H g TM H g H TM TM
Cuando se trata de grupos generadores medianos o chicos, cuya potencia no supere los 10.000
kW, es posible, aun en los casos en que la relación TW/TM resulte mayor que 0.20 lograr que el
sistema sea estable mediante el aumento de las masas rotantes agregando un volante de inercia.
Esto es válido dentro del rango dado por:

TW
0.20 < < 0.40 ..............................................................................................3
TM

Para valores superiores a 0.40 se considera imposible la regulación.

2 FORMAS DE ASEGURAR LA ESTABILIDAD

Volviendo a la expresión que establecía la condición general de estabilidad (2) se llega a:

TW LV
< 0.20 → TW < 0.20 TM → < 0.20 TM ..................................4
TM gH

Cuando el proyecto incluya grupos de generación de potencias superiores a los 10.000 kW no

resulta económicamente factible el agregado de volantes de inercia por lo tanto solo queda actuar
sobre el primero de los términos es decir tratando de reducir TW.

De las variables que definen TW tanto g como H no se pueden modificar, quedando solo analizar si
es posible actuar sobre L ó V.

2.1 Modificación de la velocidad

Para alterar, disminuir en este caso, el valor de la V se puede recurrir a dos procedimientos: uno,
el más directo se orienta a un aumento del diámetro, en tanto que el restante es haciendo uso de
los llamados descargadores sincrónicos.

2.1.1 Modificación del diámetro

A partir de la ecuación (4) se puede establecer:

LV L Q L 4Q
< 0.20 TM → = < 0.20 TM y finalmente:
gH g H Ω g H π D2

4 LQ 20 L Q
D> es decir: D> ...............................5
0.20 TM g H π π TM g H

Este sería el diámetro que provoca una disminución de la velocidad que hace estables a los
sistemas de regulación de las turbinas. Dado que el diámetro, sobre todo en conducciones
relativamente largas, tiene un peso importante desde el punto de vista económico, las
posibilidades de modificarlo son generalmente escasas. En todo caso esta alternativa deberá ser
comparada con otra que, manteniendo el diámetro económico, agregue algún dispositivo
destinado a solucionar el problema de la estabilidad.

2.1.2 Descargadores sincrónicos

CHIMENEA DE EQUILIBRIO Pág. 4

 APROVECHAMIENTO DE LOS RECURSOS HIDRICOS

Los descargadores sincrónicos permiten una lenta reducción de las velocidad dentro del conducto.
En realidad su función es solo la de prevenir las grandes sobrepresiones producto de una salida
simultánea y no programada de todos
los grupos de generación haciendo que
TURBINA
el TC resulte muy lento y por lo tanto el
valor de DH pueda mantenerse dentro QC QT
de los rangos dispuestos por el
proyectista.
DESCARGADOR
SINCRONICO
La operación de los descargadores
QD
sincrónicos se realiza cumpliéndose:

QC = QT + QD = cte

Como los descargadores solo pueden dejar o no escapar un cierto caudal, los caudales
correspondientes a la turbina tienen como opción la de disminuir su caudal. Esto por lo tanto se
vincula solamente a una caída en la demanda. Si la demanda aumentara, por lo tanto Q , para
que el descargador sincrónico colaborara en la estabilidad de la regulación debería estar T en
condiciones de dejar entrar un caudal Q cosa imposible. Por lo tanto los descargadores
sincrónicos tienen como función básica la deD controlar las sobrepresiones ante cierres parciales o
totales de las turbinas permitiendo que estas si reduzcan rápidamente su caudal de alimentación
impidiendo que se produzcan embalamientos peligrosos de los equipos generadores y evitar
sobredimensionar la conducción para que resista las sobrepresiones a que estaría sometida de no
disponerse del descargador. Las máquinas así equipadas no pueden seguir las variaciones de la
demanda y deben generar a potencia constante.
Los descargadores sincrónicos presentan dos formas características que están relacionadas con
el tipo de turbina a las que se hayan asociados.

2.1.2.1 Turbinas Pelton

En términos generales las turbinas Pelton
actúan bajo saltos importantes. El rango se
sitúa entre los 100 m y los 1 500 m. Cuando
los saltos son de varios centenares de metros
de longitud las tuberías de conducción resultan
notoriamente largas. Esta situación que resulta
bastante común en este tipo de turbinas ha
hecho necesario encontrar una solución
económica a las sobrepresiones por golpe de
ariete. De no existir los descargadores
sincrónicos las sobrepresiones demandarían el
uso de tuberías de espesores importantes y
por lo tanto elevados costos de la misma. INYECTOR

El esquema adjunto muestra la posición del

descargador sincrónico, deflector en el caso
de las Pelton, en relación al inyector y las DEFLECTOR
palas de la turbina.

Cuando es necesario reducir el caudal que

INYECTOR incide sobre las palas, el deflector, girando
QC QT
sobre su eje, desvía parcial o totalmente el
chorro. Esta acción del deflector puede
realizarse en forma tan rápida como se desee ya
DEFLECTOR
QD

CHIMENEA DE EQUILIBRIO Pág. 5

 APROVECHAMIENTO DE LOS RECURSOS HIDRICOS

que no provoca ninguna modificación del escurrimiento dentro de la conducción.

Una vez que el deflector encuentra su ubicación final, para evitar una pérdida continua de agua, la
aguja del inyector comienza a cerrar, en un proceso que, este si, se realiza muy lentamente, hasta
reemplazar al deflector.
Esta facilidad para resolver el problema de controlar las sobrepresiones elevadas con una rápida
adecuación de los caudales esta ligada al hecho particular de que en las Pelton existe un recorrido
del circuito hidráulico, entre el inyector y la pala, en que el escurrimiento se produce fuera de los
contornos cerrados facilitando así la incorporación de este elemento mecánico.

2.1.2.2 Turbinas Francis. Turbinas de reacción.

Por sus características este tipo de turbinas resulta más compleja para intercalar dentro de sus
circuito hidráulico un sistema sincrónico de descarga. Dado que se trata de un circuito cerrado
resulta imposible aplicar elementos tan sencillos como los deflectores de las turbinas de acción.

Un descargador sincrónico para turbinas de reacciones el que se muestra esquematizado en la

figura. El gráfico complementario representa el giro de los álabes del distribuidor en una operación
de cierre.

Durante la operación de un cierre parcial es necesario que, como se muestra en el esquema

general adjunto que el punto A se desplace a A’. Ante este desplazamiento giran los álabes
provocando una reducción del caudal
de QC a QT.

Simultáneamente el descargador se
abre dejando escapar un caudal QD =
QC - QT

Esta operación permite una

reducción del caudal que alimenta la
turbina sin que se haya verificado
ninguna alteración del caudal que se
desplaza por tubería.

Para impedir que el caudal QD se

pierda en forma permanente, el eje
que comanda la válvula del
descargador lleva intercalada un
dispositivo de acción retardada,
(dash pot), constituido por un pistón y
un cilindro con las cámaras A y B del
mismo llenas de aceite. Cuando el
regulador fuerza un rápido desplazamiento del eje que comanda el
descargador, el pistón y el cilindro se mueven solidariamente, aumentando
la presión en B y disminuyéndola en A.

Ambas cámaras están unidas por secciones de pasaje muy pequeñas.

A Esta vinculación entre ambas, ante la diferencia de presión permite un flujo
lento del aceite de la cámara B hacia la A permitiendo un lento cierre del
descargador. Esta modificación, por su lentitud, es tolerada por el sistema
B de regulación que adecua convenientemente las palas del distribuidor.
Esta ultima fase del proceso hace variar el caudal QC a QT y el QD a 0.
Estos dispositivos ocupan dentro de la central importantes espacios y
tienen un peso económico relativamente importante sobre todo en el caso

CHIMENEA DE EQUILIBRIO Pág. 6

 APROVECHAMIENTO DE LOS RECURSOS HIDRICOS

de las turbinas Francis. Como todo sistema mecánico agrega al sistema un factor adicional de
riesgo.

Resulta evidente que los descargadores sincrónicos solo presentan utilidad cuando se presenta
una reducción de la potencia demandada o en otras palabras cuando el proceso de regulación
exige un cierre total o parcial del caudal. Cuando la potencia debe aumentar y por lo tanto el
caudal, los descargadores no son de utilidad alguna. Su uso se restringe a las centrales de altos
saltos y largas tuberías que equipadas con estos dispositivos permiten ahorros importantes en los
costos de las tuberías por reducción de los espesores que se requerirían para soportar las
sobrepresiones por golpe de ariete. Estos estados se presentan ante la necesidad de cerrar
totalmente las turbinas generalmente en situaciones de emergencia.

Si se está ante un aprovechamiento integrado por una conducción de longitud importante y que a
su vez no tiene como función seguir las variaciones de la demanda de potencia, el uso de
descargadores sincrónicos puede resultar una solución a tener en cuenta.

Cuando no es posible satisfacer los requerimientos de la expresión 2 :

TW LV
< 0.20 → < 0.20 TM
TM gH

se debe actuar sobre la longitud L de tal manera que cumpla:

LV gH
< 0.20 TM → L < 0.20 TM
gH V

Resulta evidente que en un aprovechamiento no es posible en general modificar la longitud de la

conducción sin alterar profundamente los parámetros que determinaron su factibilidad económica.
En casi todos los casos la longitud L esta fuertemente ligada al salto H.

Veamos una obra que permite lograr un acortamiento de L en esta última ecuación sin por ello
modificar el salto.

3 CHIMENEA DE EQUILIBRIO

La chimenea de equilibrio permite una modificación total o parcial de los caudales que ingresan a
las turbinas con una mínima alteración del que escurre por la tubería de aducción. Esto lo logra sin
dispositivo mecánico alguno.

En el esquema adjunto se Presa

muestra un aprovecha-
Chimenea
miento dotado de una chi-
de Equilibrio
menea de equilibrio.

La chimenea divide la
con-ducción en dos L1
partes:
Central
L2
• Una aguas arriba del L
emplazamiento de la
chimenea que queda a
cubierto de las sobre-
presiones por golpe de
ariete.

CHIMENEA DE EQUILIBRIO Pág. 7

 APROVECHAMIENTO DE LOS RECURSOS HIDRICOS

• Una segunda, aguas abajo de la misma, donde las sobrepresiones por golpe de ariete de
desarrollan libremente.
En el esquema anterior la longitud total es L en tanto que L1 es la correspondiente al sector aguas
arriba de la chimenea y L2 al de aguas abajo.

El principio de funcionamiento de la chimenea es la de disponer de

un receptáculo que permita el ingreso de los caudales que
exceden la demanda de la Central. Se debe cumplir:

QC = QT + QD ..................................................6

Cuando el estado del escurrimiento es permanente QD = 0 , siendo

por lo tanto iguales los caudales de las tuberías L1 y L2. Al
modificarse la potencia demandada el regulador actúa sobre el
QD caudal modificando el valor de QT. Para que se modifique QC
deben cambiar las fuerzas que actúan sobre la masa de fluido que
QC se encuentra dentro de la tubería L1, y para que esto ocurra debe
modificarse la presión bajo la chimenea. Esto solo puede ocurrir si
cambia la altura de agua dentro de la misma, resultando que QT se
modifica, QC permanece inicialmente igual y por lo tanto QD
resulta distinto de cero debiendo cumplir 6.
QT
La variación de presión que se origina en la Central se propaga
hacia aguas arriba y al llegar a la sección de la chimenea deja de
progresar ya que el fluido tiene la posibilidad de ingresar o salir de ella evitando de allí en más
modificaciones bruscas del caudal.

Entre la Central y la chimenea de equilibrio el tramo L2 estará sujeto al fenómeno de golpe de

ariete y por lo tanto será esa la longitud a considerar en la condición de estabilidad buscada
pudiendo ésta expresarse:
gH
L2 < 0.20 TM ...................................................7
V

Finalmente cuando la condición de estabilidad no se cumple la solución es colocar una chimenea

de equilibrio a una distancia menor que la establecida por la expresión 7.

3.1 Dimensiones de la Chimenea de Equilibrio

Corresponde de aquí en más definir las dimensiones que debe tener la chimenea para que cumpla
correctamente sus funciones. Para disponer de alguna herramienta de cálculo que permita
correlacionar las
dimensiones con el
comportamiento de
la chimenea se
debe analizar el
comportamiento
hidráulico del
sistema a partir del
momento en que
se produce una
variación de la
potencia dando
como resultado

CHIMENEA DE EQUILIBRIO Pág. 8

 APROVECHAMIENTO DE LOS RECURSOS HIDRICOS

una modificación de los caudales.

Se analizará el comportamiento de la masa de agua dentro de la conducción comprendida entre la

embocadura y la base de la chimenea y la ubicada en el interior de la chimenea. La magnitud Y
mide la distancia entre el nivel estático de la superficie del embalse y el nivel libre de la chimenea.
Los valores Y se consideran positivos cuando la superficie libre dentro de la chimenea se
posiciona sobre el nivel del embalse.

En el tramo comprendido entre la embocadura y la base de la chimenea se pueden determinar las

condiciones generales del escurrimiento a partir de la segunda ecuación de Newton.:

dV
∑F = m a = ρ LA .................................................................................8
dt
Siendo la ΣF la resultante de todas las fuerzas externas que actúan sobre la masa líquida m. Se
partirá considerando solo la masa contenida dentro de la conducción. Esta simplificación, ignorar
la masa de agua dentro de la chimenea, en la generalidad de los casos es válida ya que dicha
masa es pequeña en relación a la de la conducción y las aceleraciones dentro de ella son también
de un orden menor. Como un paso posterior se considerará la influencia de la masa de agua
dentro de la chimenea.

Un detalle de estas fuerzas comprende:

• La fuerza generada por la presión sobre la cara en correspondencia con el plano de la

embocadura. La fuerza que actúa será: p1.A1, presión por área del conducto. Aplicando
Bernoulli en la embocadura se tendrá:
V2 V2
p1γ + + J e = H1 donde: Je = Ke son las pérdidas en la
2g 2g
embocadura expresadas en función de la velocidad.

Tanto las fuerzas como las velocidades se consideran positivas cuando dichos vectores se
dirigen hacia aguas abajo, es decir desde el embalse hacia la chimenea. De lo anterior
resulta que la fuerza actuante en la sección de entrada será:

⎡ V2⎤
F1 = γ A ⎢ H1 − (1 − K e ) ⎥ ................................................................................9
⎣ 2g ⎦

• Corresponde analizar la fuerza originada por la presión en la

sección ubicada en la base de la chimenea. Igual que en la
sección anterior, aplicando Bernoulli entre la superficie de la
chimenea y el centro del conducto en la base de la misma resulta:

p2 V2 V2
H +Y = + − Ks S
γ 2g 2g

Finalmente la fuerza en 2, base de la chimenea, actuando hacia

aguas arriba sobre la masa líquida de la tubería queda expresada
como:

⎡ V2 V2 ⎤
F2 = − γ A ⎢(H + Y ) − + KS S ⎥ .....................................................................10
⎢⎣ 2g 2 g ⎥⎦

CHIMENEA DE EQUILIBRIO Pág. 9

 APROVECHAMIENTO DE LOS RECURSOS HIDRICOS

• La tercera fuerza a considerar es la originada por la fricción de las paredes del conducto
sobre el fluido. Corresponde recordar la vinculación entre las fuerzas de fricción y la
pérdida de energía. Analizaremos un tramo de tubería en régimen permanente:

Aplicando Bernoulli entre las caras

extremas de un elemento de fluido de
longitud L se tiene:

V12 V2
L senα + R + H1 + = R + H 2 + 2 + J 1− 2
2g 2g

Tratándose de un escurrimiento en
régimen permanente se tiene que:

V1 = V2
Simplificando y reordenando la ecuación de Bernoulli anterior resulta:

J 1−2 = L senα + R + H1 − H 2

Siendo el escurrimiento permanente, la segunda ley de Newton planteada en (8)se reduce

a ΣF = 0. Siendo las fuerzas externas que actúan sobre el elemento considerado:

- F1 y F2 fuerzas originadas en la presión y por lo tanto iguales a (H1 y (H2.

- La fuerza de fricción sobre la superficie longitudinal originada por el esfuerzo

cortante ζ, que, siendo el perímetro de la sección Π resulta ζ Π L.

- La componente del peso en el sentido del escurrimiento: p sen ∀.

Por lo tanto:

ΣF = F1 − F2 + p senα − ζℵL = ( γH1 − γH 2 ) A + γALsenα − ζℵL = 0

Pudiéndose expresar las fuerzas de corte en función de las pérdidas como:

ζℵL = γA( H1 − H 2 + Lsenα ) = γAJ

Resultando finalmente:

⎡ V2 ⎤
F3 = −γAJ = −γA⎢K R ⎥ ........................................................................................11
⎣ 2g ⎦

• La última fuerza a ser considerada es la componente del peso de la masa de agua

contenida dentro de la conducción en el sentido del escurrimiento que como se acaba de
ver resulta:

F4 = P senα = γAL senα = γAJ ..............................................................................12

Teniendo en cuenta los valores definidos para las fuerzas externas actuantes resulta para
la ΣF:

CHIMENEA DE EQUILIBRIO Pág. 10

 APROVECHAMIENTO DE LOS RECURSOS HIDRICOS

⎡ V2⎤ ⎡ V2 V2⎤ ⎡ V2⎤

ΣF = γ A ⎢ H 1 − (1 − K e ) ⎥ − ⎢(H + Y ) − + K S S ⎥ − ⎢K R ⎥ + H2
⎣ 2 g ⎦ ⎢⎣ 2g 2 g ⎥⎦ ⎣ 2g ⎦

Que agrupada convenientemente se puede expresar como:

⎡ V2 V2⎤
ΣF = γ A ⎢ ( H 1 + H 2 − H − Y ) − (1 + K e − 1 + K R ) − KS S ⎥
⎣⎢ 2g 2 g ⎦⎥

Siendo: H = H1+H2 resulta: H1+H2-H-Y = -Y

y llamando a Ke+KR = ΣK que representan todas las pérdidas entre la embocadura y la

chimenea se tiene:

⎡ V2 V2⎤
ΣF = −γ A ⎢Y + ΣK + KS S ⎥ ..........................................................................................13
⎢⎣ 2g 2 g ⎥⎦

Reemplazando este valor en la (8) resulta:

dV ⎡ V2 V2⎤
ΣF = ρ LA = −γ A ⎢Y + ΣK + KS S ⎥
dt ⎢⎣ 2g 2 g ⎦⎥

Para finalmente obtener:

dV g⎡ V2 VS2 ⎤
= − ⎢Y + ΣK + KS ⎥ ............................................................................14
dt L ⎢⎣ 2g 2 g ⎥⎦

Expresión que establece las condiciones dinámicas que gobiernan las oscilaciones de masa
entre la conducción y la chimenea de equilibrio. Para lo que antecede se ha supuesto que la
conducción mantiene la sección en toda la longitud del tramo y además no se ha tomado en
cuenta la masa en el interior de la chimenea.

Las variables de la ecuación (14) son dos: V y Y. Para encontrar una

solución a dicha ecuación será necesario una segunda condición.
Esta surge de plantear las condiciones de continuidad en la sección
de la conducción en coincidencia con la chimenea como se muestra
en el gráfico adjunto. Se deberá cumplir:

Q = QCH + QT , y siendo A, ACH y AT las secciones en la

conducción, la chimenea y la conducción a la Central se tendrá como
expresión de la condición de continuidad:

V × A = VCH × ACH + VT × AT ...............................................15

La velocidad dentro de la chimenea se puede expresar como:

dY
VCH = que combinándola con la (15) resulta la:
dt

CHIMENEA DE EQUILIBRIO Pág. 11

 APROVECHAMIENTO DE LOS RECURSOS HIDRICOS

dY ACH A
V = + VT T ..............................................................................................16
dt A A

La resolución de las ecuaciones (14) y (16), que establecen las condiciones dinámicas y de
continuidad, permite analizar las relaciones entre V e Y frente a variaciones del QT y los
parámetros físicos A y ACH. En todos los casos se supone que el caudal que demandan las
turbinas es una función conocida, por lo tanto VT = f(t) resulta la condición que termina de definir el
problema.

3.2 Incorporación de la masa de agua de la chimenea en las ecuaciones

Como paso previo a la resolución de las ecuaciones planteadas consideraremos la influencia de la

masa de agua contenida entro de la chimenea. Para ello se deben plantear las condiciones
dinámicas vinculadas a la masa contenida dentro de la
misma.

dVCH
ΣF = ma = ρ ( H' +Y ) ACH .................................17
dt

Tratándose la chimenea de equilibrio de un elemento

relativamente corto y de diámetro importante se pueden
despreciar las fuerzas de fricción generadas en las paredes
de la misma. Por ello la ΕF se reduce a las presiones en la
base de la chimenea, (sección bb), y el peso de la masa
líquida contenida en su interior.

• Fuerza vertical en la sección bb. La fuerza de

presión en la base resulta de considerar Bernoulli
entre la superficie libre y la sección bb donde la presión pb resulta:

2
pb VCH
+ = H' +Y
γ 2g

⎛ V2 ⎞
de donde la pb resulta: p b = γ ⎜⎜ H' +Y − CH ⎟ y la Fbb será:
⎟
⎝ 2g ⎠

⎛ V2 ⎞
Fbb = p b ACH = γ ⎜⎜ H' +Y − CH ⎟ ACH
⎟
⎝ 2g ⎠

• Peso del líquido. El peso del líquido contenido dentro de la chimenea resulta:

P = −γACH ( H' +Y )

2
VCH
Resultando el valor de ΕF = Fbb – P = − γACH ......................................................................18
2g

Si se suma a las fuerzas de la ecuación (13) las correspondientes al líquido en la chimenea dada
por la (18) se tiene:

CHIMENEA DE EQUILIBRIO Pág. 12

 APROVECHAMIENTO DE LOS RECURSOS HIDRICOS

⎡ V2 V2⎤ V2
ΣFTotales = −γ A ⎢Y + ΣK + K S S ⎥ − γACH CH siendo AS y VS el área y la velocidad en la
⎢⎣ 2g 2 g ⎥⎦ 2g
sección de entrada a la chimenea resulta por lo tanto VCHACH=VSAS

y la ecuación anterior resulta finalmente:

⎡ V 2 ⎛⎜ AS2 ⎞ VS2 ⎤
ΣFTotales = −γ A ⎢Y + ΣK + ⎜KS + ⎟ ⎥ ..................................................................19
⎢⎣ 2g ⎝ ACH A ⎟⎠ 2 g ⎥⎦

Sumando ahora las ecuaciones (8) y (17) se tiene:

dV dVCH
Σma = ρ L A + ρ ( H' +Y ) ACH , ecuación que puede expresarse como:
dt dt

⎡ dV A dV ⎤
Σma = ρ A⎢ L + ( H' +Y ) S S
⎥ ...................................................................................20
⎣ dt A dt ⎦

igualando las ecuaciones (19) y (20) se llega a:

⎡ V 2 ⎛⎜ AS2 ⎞ VS2 ⎤ ⎡ dV A dVS ⎤

− γ A ⎢Y + ΣK + ⎜KS + ⎟ ⎥ + ρ A⎢ L + ( H' +Y ) S
⎣⎢ 2g ⎝ ⎟
ACH A ⎠ 2 g ⎦⎥ ⎣ dt A dt ⎥⎦

para llegar finalmente a:

dV ⎡ ( H' +Y ) AS ⎤ dVS g⎡ V 2 ⎛⎜ AS2 ⎞ VS2 ⎤

⎟
+⎢ = − ⎢Y + Σ K + K + ⎥ .......................................21
A ⎥⎦ dt 2 g ⎜⎝ ACH A ⎟⎠ 2 g ⎦⎥
S
dt ⎣ L L ⎢⎣

La ecuación (21) resume las condiciones dinámicas correspondientes al sistema completo tubería-
chimenea en tanto que la (14) lo hace solamente sobre la tubería. Estas condiciones, combinadas
con la (16), permiten hallar la solución del sistema.

A continuación se verá una resolución para las ecuaciones diferenciales que gobiernan las
oscilaciones en una chimenea de equilibrio.

3.3 Solución de las ecuaciones diferenciales

Una vez establecida la necesidad de intercalar una chimenea de equilibrio en la conducción

corresponde analizar cuales son los principales elementos que intervienen en el diseño estructural
y funcional de la misma para seleccionar adecuadamente las hipótesis de cálculo.

a- Dimensiones estructurales. Los requerimientos del dimensionado estructural hacen

necesario conocer los niveles extremos dentro de los cuales oscilará la superficie dentro de
la chimenea para toda alternativa operacional.

Tanto la masa de la chimenea como las pérdidas de energía son elementos que tienden a
reducir la amplitud de las oscilaciones. En un cálculo conservador resulta aceptable no
tenerlas en cuenta. En este caso las hipótesis que se adoptan son: no considerar la
influencia de las masas fuerzas en la chimenea ni las pérdidas de energía.

b- Analizada ahora desde un punto de vista funcional la chimenea, además de evitar la

propagación hacia aguas arriba de su posición de las sobrepresiones por golpe de ariete

CHIMENEA DE EQUILIBRIO Pág. 13

 APROVECHAMIENTO DE LOS RECURSOS HIDRICOS

debe producir un paulatino amortiguamiento de las oscilaciones para evitar el

mantenimiento prolongado de los estados impermanentes que se generan. Para ello es
fundamental tener en cuenta la influencia de las pérdidas ya que, como se verá, son las
únicas que posibilitan la disminución de las mismas. En estos casos las hipótesis de
cálculo pueden prescindir de la masa de la chimenea pero no de las pérdidas.

3.3.1 Solución de las ecuaciones para el caso de un fluido ideal.

Para encarar una solución con solo los requerimientos que se plantearon en el punto a es
necesario resolver las ecuaciones (14) y (16) despreciando el término correspondiente a las
pérdidas. Considerar un escurrimiento sin pérdidas es tratar con un fluido ideal, es decir un fluido
no viscoso e incompresible.

Bajo estas consideraciones la (14), ecuación que determina las condiciones dinámicas del
escurrimiento, se reduce a:

dV g
= − Y ........................................................................................................................22
dt L

Que juntamente con las condiciones de continuidad contenidas en la (16):

dY ACH A
V = + VT T
dt A A

permiten encontrar la solución buscada. Para ello es necesario determinar cual es la condición de
borde adicional que permita conocer el escurrimiento hacia las turbinas, VT = f(t), que origina el
desequilibrio de las fuerzas. Para ello basta considerar que el caudal varía de un valor QT1 a otro
QT2 en forma instantánea.

Admitiendo que las secciones A y AT son iguales la ecuación (16) se puede escribir:

dY ACH dY A
V = + VT 2 ∴ = VCH = ( V − VT 2 ) cuya derivada respecto del tiempo es:
dt A dt ACH

d 2 Y dV A dV d 2 Y ACH
= y por lo tanto: = que combinándola con la (22) resulta:
dt 2 dt ACH dt dt 2 A

d 2Y g A
+ Y = 0 ...............................................................................................................23
dt 2 L ACH

Ecuación de segundo grado en Y, incompleta de la forma:

d 2Y g A
+ mY = 0 donde m = y cuya solución es:
dt 2 L ACH

Y = a sen m t + b cos m t ..................................................................................................24

El valor de las constantes a y b se calcula a partir de las condiciones para el instante inicial t = 0.
Para dicho instante el valor de Y, nivel dentro de la chimenea medido desde el nivel del embalse,
está dado por la energía cinética en la conducción:

CHIMENEA DE EQUILIBRIO Pág. 14

 APROVECHAMIENTO DE LOS RECURSOS HIDRICOS

VT21
Y =− ................................................................................................................Condición 1
2g

Recordando que la velocidad en la chimenea aumenta cuando disminuye la velocidad hacia las
turbinas y que el caudal en la chimenea será el necesario para mantener inalterado en ese
instante el caudal de la conducción, se tiene:

dY A
= −( VT 2 − VT 1 ) .........................................................................................Condición 2
dt ACH

VT21
Para t = 0 la Condición 1 permite hallar el valor: b = −
2g

Aplicando para el mismo tiempo t = 0 la Condición 2 e la derivada de la (24) se obtiene:

A g A L A
− ( VT 2 − VT 1 ) =a m =a resultando: a = ( VT 2 − VT 1 )
ACH L ACH g ACH

Resultando finalmente:

L A g A V2 g A
Y = ( VT 1 − VT 2 ) sen t − T 1 cos t
g ACH L ACH 2g L ACH

Expresión que permite conocer las variaciones de Y en función del tiempo, de la velocidad final
hacia las turbinas, VT2, la longitud de la tubería y de las áreas de la conducción y de la chimenea.

Ocurre normalmente que la necesidad de implantar una chimenea está asociada a longitudes
importantes de la tubería y por lo tanto para evitar que las pérdidas resulten importantes las
velocidades son, en general, bajas. Esta situación hace que el coeficiente b, de energía cinética,
pueda, bajo estas condiciones, no ser tenido en cuenta.

Si se trata de analizar el comportamiento de la chimenea para un cierre total se tendrá:

VT1 = V y VT2 = 0 quedando como expresión de Y para este caso:

L A g A
Y =V sen t
g ACH L ACH

Donde las variaciones de Y resultan las una sinusoide con un valor máximo de:

L A
Y MAX = V resultando como período T de la onda:
g ACH

g A L ACH
T = 2π ∴ T = 2π
L ACH g A

3.3.2 Análisis de cierres parciales y condiciones de inestabilidad

En el punto 3.3, b-, se planteó la condición funcional de la chimenea que establecía que las
oscilaciones que en ella se desarrollaban debían amortiguarse lo más rápidamente posible. Ante

CHIMENEA DE EQUILIBRIO Pág. 15

 APROVECHAMIENTO DE LOS RECURSOS HIDRICOS

una variación de la demanda de potencia actuará el regulador de la turbina tratando de adecuar el

caudal. En un primer momento, y hasta tanto las presiones en la base de la chimenea se
modifiquen lo suficiente no se producirán modificaciones en la velocidad de la conducción y por lo
tanto se cumplirá: QCH = )QT, es decir que el caudal hacia la chimenea o desde ella será igual a la
modificación del caudal de las turbinas provocado por el regulador.

Cualquier variación en la demanda producirá una adecuación del caudal que alimenta las turbinas
hecho que generará una modificación del nivel en la chimenea. Si pensamos en alteraciones
pequeñas de la potencia tendremos dos casos extremos:

• Chimenea de gran diámetro. En este caso las pequeñas variaciones de caudal no

generarán variaciones apreciables del nivel en la chimenea y las fuerzas que actúan sobre
el fluido ubicado dentro de la conducción serán también de poca entidad y suavemente se
llegará a una nueva posición de equilibrio estable.

• Chimenea de pequeño diámetro. Aquí se presenta el caso de que pequeñas alteraciones

del caudal de generación se transforman en importantes modificaciones del nivel dentro de
la chimenea y por lo tanto se podrá producir un desequilibrio de la fuerzas que actúan
sobre el fluido de tal magnitud que este se traduzca en modificaciones del caudal de la
conducción mayor aún que el requerido por las turbinas. Esta situación impide que el
sistema converja hacia una nueva posición de equilibrio resultando por lo tanto inestable.

Una adecuada modelación matemática del estado impermanente que se desarrolla a partir de un
pequeño cambio en el requerimiento de potencia permite determinar fehacientemente si las
oscilaciones resultan convergentes o no, además de permitir apreciar claramente la duración de
este proceso transitorio. Con una herramienta de cálculo como la indicada es posible buscar la
dimensión más conveniente para la chimenea.

Sin necesidad de entrar en detalles de cálculo como los planteados siempre será necesario
conocer con que diámetro de la chimenea de equilibrio es posible asegurar el amortiguamiento de
las oscilaciones.

A continuación desarrollaremos el planteo realizado por el Dr. Dieter Thoma que permite definir el
límite mínimo del área de la chimenea para asegurar un sistema convergente de oscilaciones.

Para ello se parte de las ecuaciones (14) y (16) que definen las condiciones dinámicas y de
continuidad:

dV g⎡ V2 V2⎤
= − ⎢Y + ΣK + K S S ⎥ ....................ecuación dinámica y
dt L ⎢⎣ 2g 2 g ⎥⎦

dY ACH A
V = + VT T ....................................ecuación de continuidad
dt A A

A estas ecuaciones se les sumará el hecho de que el fenómeno a ser analizado genera valores de
Y muy pequeños frente a los valores del salto H, de tal manera que en todo lo que sigue podrán
despreciarse términos tales como:
2 2
⎛Y ⎞ Y dY ⎛ dY ⎞
⎜ ⎟ ; ; ⎜ ⎟ ; etc.
⎝H ⎠ H dt ⎝ dt ⎠

Al pasar la demanda de un valor P0 a P entre los cuales hay poca diferencia, apenas de 1 a 3%,
se deberá llegar a una posición estable tal que: P = ηgQ F H , pero inicialmente el regulador
modifica el caudal QT, que alimenta las turbinas, para adecuarlo a la demanda provocando una

CHIMENEA DE EQUILIBRIO Pág. 16

 APROVECHAMIENTO DE LOS RECURSOS HIDRICOS

alteración Y en el nivel de la chimenea y por lo tanto: P = ηgQT ( H + Y ) . Igualando ambas

expresiones de la potencia se tiene:

Q F H = QT ( H + Y ) y por lo tanto V F AT H = VT AT ( H + Y ) ∴ V F H = VT ( H + Y ) Siendo en

todos los casos P, QF y VF constantes en tanto que Y, QT y VT variables. Esto permite expresar:

H 1 1 1
VT = V F = VF expresión donde el factor es del tipo cuya solución es:
H +Y Y Y 1+α
1+ 1+
H H

1
= 1 − α + α 2 − α 3 + α 4 − K para valores de α < 1 . En este caso tratándose de valores de
1+α
Y
α = << 1 resulta que los términos α 2 , α 3 , α 4 , etc. pueden ser despreciados y por lo tanto
H
resulta:

⎛ Y⎞
VT = V F ⎜ 1 − ⎟ ......................................................................................................................25
⎝ H⎠

Si este valor se aplica en la ecuación de continuidad (16) se tiene:

dY ACH ⎛ Y⎞A
V = + VF ⎜1 − ⎟ T ...............................................................................................26
dt A ⎝ H⎠ A

Derivando esta ecuación respecto del tiempo e igualándola con la ecuación dinámica (14) se
dispondrá de una única ecuación que represente las oscilaciones de pequeña amplitud.

dV d 2 Y ACH dY V F AT
= − .................................................................................................27
dt dt 2 A dt H A

Poniendo ahora los valores de las velocidades V 2 y V S2 en función de Y en la (14) se tendrá la

ecuación general del movimiento como una ecuación diferencial en Y.

2 2 2 2 2
2 ⎛ dY ⎞ ⎛ ACH ⎞ ⎛A ⎞ ⎛Y ⎞ ⎛ A ⎞ Y AT dY ACH ⎛ Y⎞A
V =⎜ ⎟ ⎜ ⎟ + V F2 ⎜ T ⎟ + V F2 ⎜ ⎟ ⎜ T ⎟ − 2V F +2 VF ⎜1 − ⎟ T
⎝ dt ⎠ ⎝ A ⎠ ⎝ A ⎠ ⎝H ⎠ ⎝ A ⎠ H A dt A ⎝ H⎠ A

2 2
⎛ dY ⎞ ⎛ Y ⎞ dY Y
despreciando los términos donde figuren: ⎜ ⎟ ,⎜ ⎟ y se tiene:
⎝ dt ⎠ ⎝ H ⎠ dt H

2
2 ⎛
AT ⎞ Y AT dY A A
V = V F2 ⎜ ⎟ − 2V F +2 V F CH T ..................................................................28
⎝ A ⎠ H A dt A A

y como valor de V S2 resultará:

2
dY ⎛ dY ⎞ ⎛⎜ ACH ⎞⎟
VS AS = VCH ACH = ACH ∴ VS2 = ⎜ ⎟ ≅0
dt ⎝ dt ⎠ ⎜⎝ AS ⎟⎠

Aplicados a la (14) esta queda:

CHIMENEA DE EQUILIBRIO Pág. 17

 APROVECHAMIENTO DE LOS RECURSOS HIDRICOS

g⎡ ΣK ⎛⎜ 2 ⎛ AT ⎞ A A ⎞⎤
2
dV Y AT dY
= − ⎢Y + VF ⎜ ⎟ − 2VF +2 VF CH T ⎟ ⎥ .......................................29
dt L⎢ 2 g ⎜⎝ ⎝ A ⎠ H A dt A A ⎟⎠ ⎥
⎣ ⎦

Igualando las ecuaciones (27) y(29) y presentándolas en forma explícita se tiene:

1 AT ⎤ ⎡ A ⎤
2
d 2 Y dY ⎡ V F AT ΣK AT ⎤ ⎡g ΣK 1 2 ⎛ AT ⎞ A
2
− ⎢ − VF ⎥ +Y ⎢ − g VF ⎥+⎢ + ΣKV F ⎜ ⎟ ⎥=0
dt dt ⎣ H ACH L A ⎦ ⎣L L H ACH ⎦ ⎢⎣ ACH 2 L ⎝ A ⎠ ACH ⎥⎦

Expresión de la forma:

d 2Y dY
2
+β + γ Y +θ = 0
dt dt

Ecuación diferencial de segundo grado en Y, completa, cuya solución resulta una expresión
sinusoidal cuyas oscilaciones son de amplitud creciente si el valor de β resulta negativo y
decreciente si este es positivo. El valor de β es:

⎛V A ΣK AT ⎞ ⎛ ΣK 1 1 ⎞
β = −⎜⎜ F T − VF ⎟⎟ = V F AT ⎜⎜ − ⎟⎟
⎝ H ACH L A ⎠ ⎝ LA H ACH ⎠

Por lo tanto dado que es necesario que las oscilaciones se amortigüen deberá cumplirse que

β>0 es decir:

⎛ ΣK 1 ⎞
β = V F AT ⎜⎜ − ⎟⎟ > 0 desigualdad que puede expresarse:
⎝ LA ACH H ⎠

ΣK 1 LA
> y por lo tanto: ACH H > resultando como condición límite para el área de la
LA ACH H ΣK

chimenea, ACH , que:

LA
ACH > ........................Condición de Thoma
ΣK H

Si el ACH resulta menor que el establecido como valor límite por la condición de Thoma una
pequeña perturbación de Y, es decir un cambio mínimo en el valor de la potencia demandada,
generará una serie de oscilaciones de amplitud
creciente como las que se muestran en el gráfico
adjunto. En estas condiciones el regulador de la
turbina no podrá controlar la situación y él, o los
grupos deberán salir de servicio.

En la realidad la amplificación de las oscilaciones

evolucionan como se muestran en el gráfico hasta
Y
que se llega a un valor de Y tal que deje de ser
H
lo suficientemente pequeño como para que la
validez de las simplificaciones que dieron lugar a
esta solución tengan validez. Cuando se llega a

CHIMENEA DE EQUILIBRIO Pág. 18

 APROVECHAMIENTO DE LOS RECURSOS HIDRICOS

estas circunstancias las oscilaciones dejan de crecer y se mantienen dentro de un entorno más o
menos constante.

Dadas las simplificaciones que se han introducido para poder llegar a esta solución resulta
recomendable que, para chimeneas cilíndricas simples, se adopten áreas mínimas que superen
en no menos de un 50% el valor dado por Thoma.

CHIMENEA DE EQUILIBRIO Pág. 19