Tesi Dal Molin Luca
Tesi Dal Molin Luca
Tesi Dal Molin Luca
INTRODUZIONE ............................................................................................................. 1
CAPITOLO 1 .................................................................................................................... 3
CAPITOLO 2 ................................................................................................................. 15
CAPITOLO 3 ................................................................................................................. 25
CAPITOLO 4 ................................................................................................................. 33
i
4.1 Sistemi di adduzione..............................................................................................................33
4.2 Turbine ..................................................................................................................................36
4.3 Funzione di trasferimento complessiva del sistema turbina e sistema di ammissione.........42
CAPITOLO 5 ................................................................................................................. 47
CAPITOLO 6 ................................................................................................................. 53
CAPITOLO 7 ................................................................................................................. 65
CAPITOLO 8 ................................................................................................................. 85
ii
CONCLUSIONI .............................................................................................................. 91
APPENDICE 1............................................................................................................... 93
BIBLIOGRAFIA ..........................................................................................................101
iii
Indice delle figure
Figura 1.1 Rete italiana a 380 kV - Fonte sito web Terna spa (2012). ............................................. 6
Figura 1.2 Rete italiana a 220 kV – Fonte sito web Terna spa (2012).............................................. 7
Figura 1.3 Diagramma di carico giornaliero dell'intera rete italiana – Fonte sito web Terna spa.... 8
Figura 1.4 Collegamento puramente induttivo. .............................................................................. 10
Figura 1.5 Caratteristica statica compensatore. .............................................................................. 11
Figura 2. 1 Schema a blocchi regolazione primaria ........................................................................ 19
Figura 2. 2 Schema regolazione primaria - Simulink ..................................................................... 21
Figura 2. 3 Grafico risultante simulazione - Simulink .................................................................... 21
Figura 2. 4 Schema a blocchi regolazione secondaria .................................................................... 22
Figura 2. 5 Schema a blocchi regolazione secondaria con Simulink .............................................. 23
Figura 2. 6 Grafico dell'andamento della frequenza con regolazione secondaria ........................... 23
Figura 3.1 Schema a blocchi princio della regolazione della velocità di un gruppo....................... 26
Figura 3.2 Schema a blocchi linearizzato della regolazione di velocità di un gruppo .................... 27
Figura 3. 3 Diagramma di Bode fdt 3.21 ........................................................................................ 30
Figura 3. 4 Caratteristica statica frequenza-potenza ....................................................................... 31
Figura 3. 5 Caratteristica statica frequenza-potenza in p.u. ............................................................ 31
Figura 4. 1 Sistema di adduzione di una centrale idroelettrica ....................................................... 33
Figura 4. 2 Condotta forzata di un impianto idroelettrico............................................................... 34
Figura 4. 3 Schema gruppo di generazioni di un impianto idroelettrico ......................................... 36
Figura 4. 4 Scelta della turbina in funzione di salto e portata disponibile. [14] ............................. 37
Figura 4. 5 Regolazione della portata per una turbina Pelton. Fonte www.energialternativa.org .. 37
Figura 4. 6 Caratteristica potenza-rendimento di una turbina Pelton. [14] ..................................... 38
Figura 4. 7 Girante di una turbina Pelton a cinque getti. Fonte www.atmosferis.com ................... 38
Figura 4. 8 Turbina Francis con sistema di regolazione del distributore. ....................................... 39
Fonte italian.hydropower-turbine.com ........................................................................................... 39
Figura 4. 9 Regoalazione inclinazione pale del distributore. Fonte [14] ........................................ 39
Figura 4. 10 Pala regolabile di un distributore di una turbina Francis. Fonte [14] ......................... 39
Figura 4. 11 Caratteristica generica di regolazione di una turbina Francis. Fonte sito www.cap2.it
.............................................................................................................................................................. 40
Figura 4. 12 Girante di una turbina Kaplan. Fonte sito www.tbhic.cn ........................................... 40
Figura 4. 13 Diagramma di Bode della funzione di trasferimento del sistema turbina e sistema di
ammissione........................................................................................................................................... 44
Figura 4. 14 Grafico caratteristiche potenza-apertura valvole per le diverse turbine ..................... 45
Figura 5. 1 Schema generale impianto termoelettrico .................................................................... 47
v
Figura 5. 2 Schema a blocchi linearizzato della caldaia ..................................................................48
Figura 5. 3 Caratteristica statica portata-apertura valvole per immissione ad arco totale ...............49
Figura 5. 4 Caratteristica statica portata-apertura valvole per immissione ad arco parziale ............49
Figura 5. 5 Schema a blocchi linearizzato della turbina e della caldaia di un impianto
termoelettrico ........................................................................................................................................52
Figura 6. 1 Schema a blocchi regolatore di velocità accelerometrico .............................................55
Figura 6. 2 Schema a blocchi regolatore di velocità con retrazione transitoria ...............................55
Figura 6. 3 Schema a blocchi ciclo di regolazione della frequenza .................................................57
Figura 6. 4 Diagramma di Bode (modulo) del ciclo di regolazione della frequenza .......................57
Figura 6. 5 Schema a blocchi linearizzato di un regolatore di velocità per gruppi termoelettrici ...59
Figura 6. 6 Schema di principio della regolazioni di un gruppo termoelettrico ..............................60
Figura 6. 7 Schema a blocchi semplificato per grandi variazioni del regolatore di velocità. della
turbina e del sistema di adduzione di un impianto idroelettrico ............................................................62
Figura 6. 8Schema a blocchi semplificato per grandi variazioni del regolatore di velocità. della
turbina e del sistema di adduzione di un gruppo termoelettrico ............................................................63
Figura 6. 9 Schema a blocchi semplificato per grandi variazioni per i gruppi con turbine a gas ....64
Figura 7. 1 Tipico regolatore di velocità per gruppi idroelettrici ....................................................66
Figura 7. 2 Caratterisitca statica rettilinea del gruppo di generazione .............................................68
Figura 7. 3 Pendolo di Watt. Fonte www.istitutomontani.it ............................................................69
Figura 7. 4 Regolatore di velocità oleodinamico della centrale idroelettrica del Furlo lungo il fiume
Candigliano. Fonte www.educambiente.it ............................................................................................70
Figura 7. 5 Regolatore di velocità oleodinamico. Fonte http://www.educambiente.it/ ...................71
Figura 7. 6 Regolatore di velocità per turbina a vapore. Fonte http://www.museibologna.it/ .........71
Figura 7. 7 Schema di un regolatore di velocità di un impianto idroelettrico..................................72
Figura 7. 8 Vista in sezione dell'impeller e dell'albero della turbina. Fonte [4] ..............................73
Figura 7. 9 Vista in sezione del servomotore di comando valvole. Fonte [4]..................................73
Figura 7. 10 Ruota fonica con sensori induttivi. Fonte [7] ..............................................................74
Figura 7. 11 Ruota fonica con sensori ottici. Fonte [7] ...................................................................75
Figura 7. 12 Principio di funzionamento di un LVDT. Fonte dito www.efunda.com .....................76
Figura 7. 13 LVDT installato sul distributore della turbina. Fonte [7] ............................................76
Figura 7. 14 Schema regolatore di velocità con attuatore ................................................................77
Figura 7. 15 Controllore a microprocessori Emerson Ovation per il controllo della regolazione di
velocità con moduli ingresso/uscita. Fonte [4]......................................................................................78
Figura 7. 16 Schema regolatore di velocità SIMO. Fonte [9] ..........................................................79
Figura 7. 17 Risposta in frequenza del sistema idraulico adottato per la simulazione. Fonte [9] ....80
Figura 7. 18 Risposta al gradino del sistema di regolazione SIMO.................................................81
Figura 7. 19 Risposta del regolatore alla variazione del carico elettrico .........................................81
Figura 7. 20 Schema dei regolatori pre e post revamping ...............................................................82
vi
Figura 7. 21 Risposta al gradino dei regolatori pre e post revamping ............................................ 82
Figura A1. 1 Libreria Simulink ...................................................................................................... 94
Figura A1. 2 Foglio Simulink ......................................................................................................... 95
Figura A1. 3 Finestra impostazione parametri blocco G1 .............................................................. 96
Figura A1. 4 Finestra impostazione parametri blocco disturbo deltaP ........................................... 96
Figura A1. 5 Finestra salvataggio risultati delle simulazione ......................................................... 97
Figura A1. 6 Finestra di configurazione della simulazione ............................................................ 97
Figura A1. 7 Barra Simulink .......................................................................................................... 98
Figura A1. 8 Grafico risultante della simulazione .......................................................................... 98
Figura A1. 9 Grafico risultante con comparazione di diverse simulazioni al variare di un parametro
.............................................................................................................................................................. 99
vii
Introduzione
Il lavoro di tesi svolto tratta i regolatori di velocità delle centrali termoelettriche e idroelettriche.
Questi dispositivi sono fondamentali per il corretto funzionamento di un qualsiasi sistema elettrico di
trasmissione. Il loro obiettivo è quello di mantenere la frequenza del sistema a un valore costante e il
più possibile prossimo al suo valore nominale.
In questo elaborato ho approfondito le parti che compongono l’anello di regolazione di velocità
comprendendo anche il sistema di adduzione (condotta forzata per impianti idroelettrici e caldaia per
impianti termoelettrici) e la relazione apertura valvole-potenza della turbina. Si sono quindi definiti
quelli che sono i requisiti richiesti al regolatore per avere un ciclo di regolazione stabile.
Successivamente si sono analizzati come erano costruiti nel passato i regolatori di velocità,
generalmente oleodinamici, descrivendo il loro funzionamento al variare della velocità di rotazione
del gruppo controllato. Successivamente vengono presentati alcuni degli interventi che vengono
eseguiti su impianti esistenti per modernizzare i regolatori di velocità e migliorarne le prestazioni.
Viene poi mostrata una tipologia di controllo innovativa basata sull’acquisizione di più segnali
dell’anello di regolazione. Questo comporta una complicazione nella progettazione delle funzioni di
trasferimento del regolatore di velocità, ma consente di ottenere una regolazione migliore rispetto ai
classici metodi di controllo con regolatori PI o PID.
Nell’ultimo capitolo sono stati riportati invece i requisiti che il gestore della rete impone ai
regolatori di velocità dei rispettivi gruppi di produzione.
1
In appendice è stato descritto l’utilizzo del programma Simulink che è stato usato per effettuare
alcune simulazioni presenti nell’elaborato.
2
Capitolo 1
3
Andiamo ora ad analizzare in dettaglio le varie parti del sistema.
Il compito di questa parte del sistema è la produzione di energia elettrica, ovvero, convertire una
fonte di energia primaria in energia elettrica. Questo processo avviene nelle centrali elettriche di
produzione. Prima che la fonte primaria si trasformi in energia elettrica, deve subire diverse
trasformazioni. Per esempio, in una centrale termoelettrica, il carbone viene usato per produrre il
calore necessario a riscaldare l'acqua della caldaia, che viene inviata in una turbina a vapore. Il
generatore elettrico è generalmente accoppiato direttamente con la turbina, la quale lo mantiene in
rotazione durante la produzione di energia elettrica. Le fonti maggiormente utilizzate per la
produzione di energia sono:
− carbone;
− nucleare;
− gas;
− petrolio.
Per quanto riguarda le fonti rinnovabili, abbiamo:
− eolico;
− solare;
− idroelettrico.
Tralasciando la parte termica e meccanica dell'impianto, specifichiamo che solitamente i
generatori producono energia elettrica in media tensione (dai 10 kV ai 20 kV). Questa viene poi
innalzata da un trasformatore fino a livelli di alta o altissima tensione (220 kV o 380 kV), per poter
essere così trasmessa attraverso le linee elettriche di trasmissione. Per quanto riguarda impianti di
generazione a energia rinnovabile, visto che le potenze potrebbero essere limitate, si hanno molti casi
in cui l'energia elettrica viene direttamente utilizzata oppure immessa in reti di media (20 kV) o bassa
tensione (400 V).
Questa parte del sistema non si limita alla sola produzione di energia, ma la sua importanza è
anche dovuta al fatto che è un punto nevralgico per la regolazione della frequenza e in parte anche
della tensione della rete. Infatti per avere una frequenza di rete costante, è necessario che in ogni
istante sia pareggiato il bilancio delle potenze attive prodotte e richieste (incluse le perdite di rete). In
poche parole, variando la potenza prodotta dalle centrali, si può seguire la variazione del carico e
quindi mantenere la frequenza costante. Regolando invece la tensione in uscita dai generatori, si
possono modificare i flussi di potenza, e quindi ottenere una più economica alimentazione del carico e
una tensione in arrivo adeguata (la regolazione della tensione avviene anche localmente in punti
strategici della rete).
4
1.1.3 Trasmissione[2]
Questa parte della rete si occupa della trasmissione di grandi quantità di energia elettrica a grandi
distanze (nell'ordine delle centinaia di chilometri). È costituito da tutte le linee di alta e altissima
tensione che attraversano il paese, e formano una rete magliata a cui vengono connesse le più grandi
centrali di produzione. I rami della rete possono essere soggetti a flussi bidirezionali di potenza,
permettendo una maggiore flessibilità del sistema. È sempre in questa rete che avvengono le
interconnessioni con i sistemi elettrici dei paesi confinanti come Austria, Svizzera e Francia, e anche
particolari carichi elettrici che richiedono potenze elevate (trazione ferroviaria ad alta velocità e
laboratori di ricerca).
Esiste poi il sistema di subtrasmissione, che opera a un livello di tensione di 132 kV, sempre
magliato, e assicura l'inter collegamento tra i siti di prelievo sulla rete di trasmissione e le cabine
primarie da cui originano le reti capillari a tensione inferiori. A questa rete sono anche connesse
grosse utenze industriali e altre centrali di produzione.
Le linee che operano alla tensione di 220 o 380 kV, sono le uniche in Italia a presentare i
conduttori a fascio (due o più conduttori in parallelo per fase), e sono sostenute da tralicci tronco-
piramidali. Sono presenti anche linee a doppia terna su un unico traliccio.
Nelle figure seguenti vengono rappresentate le linee elettriche di alta e altissima tensione presenti
sul territorio italiano.
5
Figura 1.1 Rete italiana a 380 kV - Fonte sito web Terna spa (2012).
6
Figura 1.2 Rete italiana a 220 kV – Fonte sito web Terna spa (2012).
7
1.1.4 Distribuzione[2]
Il sistema di distribuzione si divide su due livelli di tensione: media (20 kV) e bassa tensione (400
V). La rete è configurata in modo radiale o diramato, con possibilità di contro alimentazione dei vari
rami, ed è da questa che vengono alimentati la maggior parte delle utenze industriali e domestiche.
Negli ultimi anni, vista la grande diffusione della generatore distribuita, a questa rete sono
connessi anche gli impianti di autoproduzione residenziale e industriale, provocando diversi problemi
di proteggibilità e qualità della forma d'onda della tensione. La maggior parte delle linee sono in cavo,
con, soprattutto in media tensione, presenza di linee aeree.
Figura 1.3 Diagramma di carico giornaliero dell'intera rete italiana – Fonte sito web Terna spa.
Il diagramma di carico, o curva di carico, è dunque la curva della potenza attiva prelevata
dall'utenza in funzione del tempo.
8
In funzione del periodo di tempo preso in considerazione si possono avere diagrammi di carico:
− giornalieri;
− settimanali;
− annuali.
Sfruttando il fatto che questi diagrammi sono prevedibili e comunque molto simili tra loro
considerando intervalli di tempo uguali, è possibile prevedere in anticipo il fabbisogno di energia e la
richiesta massima di potenza. Questo viene fatto dal gestore della rete di trasmissione ogni giorno per
il giorno dopo, in modo da organizzare al meglio il sistema e assicurare la fornitura di energia e
potenza ai carichi. (Fellin & Benato, 2011)
La tensione della rete elettrica è uno di quei parametri che contribuiscono a definire la qualità
dell'energia fornita ai carichi. Infatti le varie utenze allacciate, funzionano correttamente e con le
migliori prestazioni se la tensione che le alimenta è quella nominale, cioè la tensione per cui sono
state progettate. Viceversa, se la tensione si discosta molto dal valore nominale, le prestazioni
degradano anche molto velocemente (si pensi per esempio alla coppia dei motori asincroni che varia
col quadrato della tensione, o la luminosità delle lampade, anch'essa proporzionale alla tensione). Le
variazioni del modulo della tensione sono conseguenza di disturbi che si vengono a verificare nella
rete, come variazione del carico (soprattutto potenza reattiva), e di azioni di controllo su sorgenti di
potenza reattiva (generatori sincroni, reattori, condensatori, variatori sotto carico).
La regolazione della tensione è l'insieme di tutte le azioni di controllo attuate in varie parti del
sistema con l'obiettivo di mantenere il livello di tensione a un livello idoneo per il funzionamento
corretto di impianti e dispositivi connessi.
9
Per capire come agire per regolare la tensione nella rete di trasmissione, si deve prima studiare il
regime delle tensioni e delle potenze in un collegamento puramente induttivo, che rappresenta una
approssimazione di un tratto di linea (si può ritenere una buona approssimazione in quanto la parte
resistiva della linea è trascurabile rispetto alla componente induttiva). Lo schema del tratto di linea
sarà dunque come segue:
$ % $
! "# % &'( ! "% (1.1)
Dove le grandezze sono espresse in per unità, e è l'angolo di sfasamento tra e . Separando
nella (1.1) la parte reale e immaginaria, si ottengono:
sin , (1.2)
cos
. (1.3)
Dalle relazioni ottenute si può osservare come la potenza attiva dipenda principalmente dallo
sfasamento delle due tensioni, mentre la potenza reattiva sia legata alla differenza dell'ampiezza delle
due. Infatti, considerando sufficientemente piccolo1:
≃ , (1.4)
≃ .2 (1.5)
Si può inoltre specificare che, nel collegamento preso in esame, la potenza reattiva positiva fluisce
dal nodo a tensione più alta al nodo a tensione più bassa.
Possiamo quindi concludere che, la caduta di tensione lungo una linea di trasmissione, è in gran
parte dovuta alla potenza reattiva circolante. Volendo quindi limitare la caduta di tensione nei diversi
nodi della rete, si deve evitare di trasmettere consistenti potenze reattive su grandi distanze. Ciò
significa che è necessario immettere in rete le potenze reattive in nodi non lontani da quelli in cui esse
vengono richieste, lasciando a gran parte della rete di trasmissione soprattutto il compito di
convogliare potenze prevalentemente attive: sarà perciò opportuno che i generatori delle centrali
1
cos ≃ 1 ; sin ≃
2
| |
10
lontane dai centri di consumo vengano scarsamente impegnati nella fornitura di potenza reattiva, che
dovrà essere invece prodotta possibilmente in vicinanza dei carichi (Paolucci, 1998).
Un primo provvedimento che va verso la riduzione della potenza reattiva circolante in rete è il
rifasamento dei carichi, "imposto" agli utenti attraverso la tariffazione del kvarh assorbito dalla rete. I
carichi rifasati presenteranno quindi un maggior fattore di potenza e quindi minor potenza reattiva
richiesta dalla rete. Questo unico provvedimento non è però sufficiente e spetta al gestore della rete
completare e regolare la compensazione delle potenze reattive, in modo da controllare il regime delle
tensioni nei nodi della rete. Infatti i carichi non vengono mai rifasati completamente, e anche singoli
elementi della rete assorbono potenza reattiva (per esempio trasformatori e linee a vuoto). Il gestore
deve quindi considerare questo aspetto e munirsi di apparecchiature che possono erogare o assorbire
potenza reattiva quando è necessario. I principali mezzi con cui è possibile immettere in rete potenze
reattive sono due:
− compensatori sincroni,
− condensatori,
− reattori.
I compensatori sincroni sono macchine sincrone che non forniscono potenza attiva alla rete e che
non sono accoppiate né a turbine, cioè motori primi, né cedono potenza meccanica ad utilizzatori. La
potenza attiva necessaria per far fronte alle perdite della parte elettromagnetica e nella parte
meccanica viene assorbita dalla rete. I compensatori sincroni possono essere considerati come dei
motori sincroni funzionanti, dato il piccolo valore delle perdite di potenza attiva, a vuoto. A regime
permanente, la corrente statorica è allora in quadratura con la tensione. A seconda dell'eccitazione
della macchina, questa può erogare in rete potenza reattiva positiva (funzionamento in
sovraeccitazione, comportamento da condensatore) o negativa (funzionamento in sottoeccitazione,
comportamento come reattore), come si può vedere nella figura 1.5.
La variazione di potenza reattiva viene ottenuta variando l'eccitazione della macchina mediante un
regolatore di tensione. In condizioni di regime, cioè v mantenuta costante dal regolatore, la potenza
reattiva q varia
linearmente con la
tensione di eccitazione vf
secondo una caratteristica
simile a quella di figura
1.5.
Se, ad esempio, la
tensione del nodo in cui è
collegato il compensatore,
o dei nodi vicini,
diminuisce per effetto di
Figura 1.5 Caratteristica statica compensatore. un aumento del carico, il
11
regolatore di tensione determina un aumento di vf e quindi di q. nelle nuove condizioni di regime la
tensione v ritorna al valore precedente con un'eccitazione più grande e quindi con una maggiore
erogazione di potenza reattiva. Il compensatore quindi svolge il compito di stabilizzazione della
tensione.
Queste speciali macchine sincorne sono molto efficaci per regolare la tensione, ma il loro costo
elevato ne sconsiglia l'uso nelle reti di distribuzione. Sono quindi presenti solo in particolari nodi della
rete di trasmissione AT.
Per le reti di distribuzione e subtrasmissione vengono usate batterie di condensatori statici in
derivazione. La potenza reattiva che queste batterie erogano in rete, è pari a:
)* 3, % -. (1.6)
con E tensione stellata della rete, ω pulsazione angolare, C capacità di una fase. La potenza
dipende quindi solo dalla tensione della rete, e varia con questa secondo una legge quadratica; può
quindi aumentare o diminuire, con effetto opposto a quello di un compensatore sincrono 3 . I
condensatori statici esercitano quindi un'azione destabilizzante sulla tensione. La capacità installata
viene suddivisa in diverse unità con possibilità di disinserzione mediante dispositivi di manovra
idonei.
I reattori, sempre inseriti in derivazione, assorbono una potenza reattiva:
3, %
)/ (1.7)
-0
Come i condensatori, vengono installati in diverse unità disinseribili, e sono utili in quelle situazioni
in cui la rete di trasmissione di grande estensione richiede a vuoto4, o con bassi carichi attivi, una
considerevole quantità di potenza reattiva dovuta alle capacità delle linee. In queste situazioni viene
definito un grado di compensazione trasversale della linea, che è il rapporto tra la suscettanza
capacitiva totale della linea e quella induttiva dei reattori. Anche questi reattori, come le capacità,
sono elementi passivi, ma hanno un'azione stabilizzante rispetto alla tensione della rete5.
Per quanto riguarda la rete di distribuzione, oltre al rifasamento, per mantenere la tensione a livelli
accettabili, si ricorre all'uso di commutatori sottocarico (nelle cabine primarie AT/MT) e di variatori a
vuoto (nelle cabine sencondarie MT/BT). Questi dispositivi consentono di variare il rapporto spire del
trasformatore e quindi modificare la tensione al secondario (questo richiede comunque un aumento o
una diminuzione della potenza reattiva transitante nella linea che alimenta il trasformatore). Un uso
non controllato di questi dispositivi, soprattutto variatori sottocarico (che intervengono
automaticamente), potrebbe portare instabilità nella rete, fino ad arrivare al collasso della tensione.
3
All'aumentare della tensione, assorbono una potenza proporzionale al quadrato di quest'ultima,
ovvero iniettano in rete una potenza superiore.
4
effetto Ferranti, specialmente in linee molto lunghe (l > 100 km).
5
All'aumentare della tensione, assorbono una potenza proporzionale al quadrato di quest'ultima,
ovvero iniettano in rete una potenza minore.
12
1.3.2 Cenni sulla regolazione della frequenza[3]
2
nota relazione:
.
60
1 (1.8)
6
Un aumento del carico può essere anche visto come la perdita di una centrale di produzione,
viceversa per la diminuzione del carico.
13
Capitolo 2
2.1 Premessa
Come già brevemente descritto nel capitolo precedente, la regolazione della frequenza è necessaria
nel momento in cui si viene ad avere uno squilibrio tra potenza attiva generata e potenza attiva
assorbita dai carichi. In seguito a questo fenomeno, si ha l’inizio di un transitorio in cui la frequenza
della rete varia in funzione dell’entità dello squilibrio, delle caratteristiche del parco generatori e dalle
caratteristiche dei carichi in quel momento connessi. In particolare possiamo dividere la regolazione
della frequenza in tre fasi:
− regolazione primaria,
− regolazione secondaria,
− regolazione terziaria.
Queste tre fasi si susseguono nel tempo nell’ordine scritto, e la loro durata totale può essere di
diversi minuti.
Le modalità e i tempi con cui i gruppi di produzione intervengono per la regolazione della
frequenza dipendono dalle caratteristiche delle unità di produzione, dal regolatore di velocità e dalle
disposizioni date dai gestori della rete di trasmissione (per l’Italia è Terna S.p.a.).
Nel seguito del capitolo prenderemo come disturbo che varia il regime della frequenza di rete un
aumento a gradino della potenza richiesta dei carichi. Questo disturbo verrà identificato con il nome
∆P.
15
2.2 Regolazione primaria[3]
La regolazione primaria ha inizio nel momento in cui in rete si verifica un disturbo che va a
rompere l’equilibrio tra potenza generata e assorbita, con conseguente variazione della frequenza. Con
la frequenza variano anche le potenze assorbite dai carichi, e grazie all’intervento del regolatore di
velocità, anche la potenza immessa in rete dalle unità di produzione. Si instaurerà quindi, alla fine
della regolazione primaria, una nuova condizione di regime in cui la potenza assorbita e generata si
equivalgono, ma la frequenza della rete sarà diversa da quella antecedente il disturbo.
Per determinare il modo in cui la frequenza varia in questo intervallo di tempo, è necessario
studiare il bilancio delle potenze attive, assumendo per ora la presenza di una sola unità di produzione
che alimenta un insieme di carichi, e per semplicità un rendimento unitario dell’insieme turbina-
alternatore (potenza immessa in turbina uguale a quella elettrica generata dall’alternatore). Con le
ipotesi assunte, prima che si verifichi il disturbo il bilancio delle potenze sarà:
56∗ = 58∗
Dove con 56∗ si intende la potenza attiva immessa in rete dal gruppo turbina alternatore in una data
condizione di regime e con 58∗ la potenza attiva assorbita dai carichi a regime.
Nel momento in cui si presenta in rete un disturbo ∆5 assunto come un aumento a gradino della
potenza assorbita dai carichi, la frequenza varia e quindi anche l’assorbimento dei carichi connessi e la
potenza immessa in rete dai generatori. Un altro fattore che riveste importanza nei transitori di
frequenza è l’energia cinetica delle macchine elettriche rotanti connesse (motori e generatori asincroni
e sincroni). A causa della variazione di frequenza, cambia anche l’energia cinetica di queste macchine
(poiché ad una variazione di frequenza corrisponde una variazione della loro velocità di rotazione).
L’energia cinetica delle macchine si oppone alla variazione della frequenza, poiché esse tendono a
mantenere la loro velocità. Di conseguenza si deve definire una grandezza che entra in gioco e che
viene chiamata potenza accelerante, definita come
=*6>
∆5 :;< =
;
(2.1)
si dovrà calcolare questo contributo dovuto da una parte dalle unità generatrici e dall’altra dai carichi
rotanti connessi e tenerne conto nel bilancio delle potenze.
A seguito del disturbo ∆5, tenendo conto dei cambiamenti delle potenze assorbite e generate da
parte dei carichi e delle unità di produzione e infine delle potenze acceleranti, il bilancio delle potenze
è il seguente:
56∗ + ∆56 :;< − ∆5 ? :;< = 58∗ + ∆58 :;< + ∆5 @ :;< + ∆5:;< (2.2)
con il seguente significato dei simboli:
∆56 variazione di potenza immessa dal gruppo di generazione a seguito della variazione di
frequenza,
∆5 ? variazione di potenza accelerante dovuta alle unità di produzione,
∆58 variazione di potenza assorbita dai carichi a seguito della variazione di frequenza,
∆5 @ variazione di potenza accelerante dovuta ai carichi connessi.
16
Il segno meno attribuito alla ∆5 ? è giustificato dal fatto che, per esempio, per un calo di frequenza
si ha una variazione di energia cinetica negativa, e quindi con il segno attribuito nel bilancio delle
potenze si somma a 56∗ + ∆56 (la potenza accelerante dà quindi un apporto effettivo di potenza a
quella immessa in rete).
Considerando che 56∗ = 58∗ , e unendo in un unico fattore la variazione di energia cinetica del sistema
(cioè accorpando ∆5 ? e ∆5 @ ), si ottiene
∆56 :;< − ∆5:;< = ∆58 :;< + ∆5 :;< (2.3)
con
=*6>
∆5 :;< = ∆5 ? :;< + ∆5 @ :;< =
;
(2.4)
Per arrivare a calcolare lo scarto Δf al termine della regolazione primaria, è necessario conoscere i
valori assunti da ciascun addendo dell’equazione.
Come già accennato in precedenza, l’assorbimento di potenza da parte dei carichi connessi subisce
una certa variazione in funzione dell’andamento della frequenza. Nell’equazione (2.4) questa
variazione è chiamata ∆5 @ ed è necessario quantificarla per tenerne conto durante i transitori di
frequenza. Considerando che normalmente le variazioni di frequenza sono molto limitate, possiamo
assumere che, la potenza assorbita dai carichi in funzione della frequenza sia:
1:;<
C
5@ :;< = 5@∗ A B
1∗
(2.5)
Vengono ovviamente trascurate nella formule la variazioni di tensione e altri fenomeni secondari.
L’esponente α è sempre superiore a uno e assume valori compresi tra 1,5 ÷ 2. Questo è dovuto al fatto
che la maggior parte dei carichi connessi è composto da motori asincroni, i quali, per le caratteristiche
dei carichi meccanici, devono fornire una coppia tanto maggiore all’aumentare della velocità (imposta
dalla frequenza di rete); per esempio, per pompe centrifughe, α vale 3 ÷ 4 e per ventilatori 2 ÷ 3. Se la
coppia richiesta non varia con la velocità di rotazione, si avrà α = 1. Altri carichi, per esempio di tipo
resistivo, presentano α = 0.
Possiamo ora calcolare la derivata rispetto al tempo della (2.5), che secondo le note regole di
derivazione, risulta
5@ 1 C&? 1
= 5@∗ D E ∗ F
; 1 1∗
(2.6)
Per rendere più semplice l’equazione e linearizzare il comportamento dei carichi nell’intorno di 1 ∗ ,
possiamo assumere che # ∗( = 1. Quindi per incrementi finiti della frequenza, si ha:
G
G
5@∗ D
∆5@ :;< = ∆1:;< = H@ ∆1:;<
1∗
(2.7)
17
Il parametro H@ ha le dimensioni di un’energia A = N ∗ O = P B. Ad esso viene dato il nome di
I
K
J M
L
La formula che esprime l’energia cinetica di un corpo in rotazione attorno ad un asse è la seguente:
1
=*6> = P -%
2
(2.8)
L’energia cinetica dipende quindi dal quadrato della velocità angolare -. La velocità di rotazione è
imposta dalla frequenza del sistema, quindi, l’energia cinetica in funzione della frequenza di rete si
può esprimere con la seguente equazione:
1:;<
%
=*6> = = ∗ A B
1∗
(2.9)
2 =∗
=*6> = ∆1
; 1∗ ;
(2.12)
Definendo HI =
% I∗
G∗
possiamo brevemente scrivere che
∆5 = HI ∆1
;
(2.14)
1 %
=W = PΩ
2 W
(2.17)
18
Dove con 5W si intende la potenza attiva del generatore, con 5Y quella del generico motore e con Q W e
Q Y i corrispondenti tempi di avviamento. Raccogliendo 5W nella (2.18), si può ottenere il tempo di
avviamento della rete Q , e quindi il termine 2 = ∗ risulta essere espresso dalla seguente equazione
1
2 =∗ 5W :Q Z 5Y Q Y < 5W Q
W
5W (2.19)
Y
Possiamo ora sostituire nella (2.3) le equazioni (2.7) e (2.14) ricavate ottenendo così l’equazione:
Per proseguire e ricavare la risposta del sistema è necessario passare al dominio della frequenza con
l’utilizzo della trasformata di Laplace. Quindi, la (2.20) diventa
∆56 ∆5 ∆[:H@ O HI < (2.21)
Quindi la variazione di frequenza ∆[ è data dalla (2.20)
∆[ :∆56 ∆5)
?
:\] ^_ \` <
(2.22)
2 =∗ 5W Q W Z 5Y Q Y (2.18)
Y
Per conoscere ∆[, manca soltanto il contributo ∆56 . Il suo valore nel tempo dipende dal regolatore di
velocità e dalle caratteristiche del gruppo di produzione. Come si capirà meglio nel capitolo 3, nella
regolazione primaria non si agisce sul variagiri del regolatore, ma si ha soltanto il contributo del
dispositivo accelerometrico che rileva la variazione di frequenza e va a modificare la posizione del
distributore della turbina, variando quindi la potenza immessa in rete del valore ∆56 . Con le ipotesi
fatte, la funzione di trasferimento che lega ∆56 a ∆[ è
1 O Q?
∆56 Ha : ∆[<
1 O Qa
(2.23)
19
Semplificando lo schema a blocchi e riducendolo ad un unico blocco, si ha
Per verificare la stabilità del sistema è necessario valutare le radici del’equazione caratteristica. Per
facilitare la comprensione, definiamo
Ha H@
-U c
Qa HI
(2.25)
1 H@ Qa HI Ha Q?
d
2 eQa HI :Ha H@ <
(2.26)
D?,% -U d f -U ed % 1 (2.29)
e presentano parte reale negativa (escludendo l’ipotesi di d 0). In base a queste considerazioni
possiamo dire che nel tempo la variazione di frequenza della rete potrà avere un andamento aperiodico
oppure oscillatorio smorzato (questo dipende dal valore che assume il coefficiente d).
In base allo schema a blocchi di figura, possiamo ricavare
1 -U% -U% Qa O
∆[ : ∆5<
Ha H@ O% 2d-U O -U%
(2.29)
20
Si può ricavare l’andamento della variazione della frequenza nel tempo antitrasformando
l’equazione (2.29), oppure simulare il sistema al calcolatore con programmi di calcolo appropriati, per
esempio utilizzando Simulink.
Nel nostro caso andremo a ricavare l’andamento di ∆1:;< con l’uso del programma di calcolo
appena citato.
21
2.3 Regolazione secondaria[3]
La regolazione secondaria della frequenza ha il compito finale di riportare la frequenza di rete al
valore nominale a seguito di un disturbo che ne ha causato la variazione. Abbiamo visto come a
seguito della regolazione primaria, viene ristabilito l’equilibrio tra potenza richiesta e potenza
immessa nel sistema, ma a una frequenza superiore o inferiore a quella prima del disturbo. Per
riportare la frequenza al valore nominale del sistema si dovrà quindi immettere o diminuire la potenza
generata dai vari gruppi collegati. Questa variazione di potenza sarà proporzionale allo scarto Δf che si
verifica in seguito alla regolazione primaria; dovrà inoltre intervenire solo a regolazione primaria
terminata.
Come abbiamo visto nel precedente capitolo, per variare la potenza immessa in rete dal generatore,
si deve agire sul variagiri del regolatore di velocità. Per variare la posizione del variagiri, nel passatto
si agiva manualmente, vista la presenza di personale in centrale. Ora l’azione è automatica; è
necessario quindi un segnale che indichi all’automazione di centrale di agire sul variagiri. Questo
segnale elaborato arriva alle centrali che effettuano la regolazione secondaria dal Centro Nazionale di
Controllo (CNC) di Roma. Il segnale è convenzionalmente compreso tra 0 e100 in relazione alla
potenza massima che ogni centrale è in grado di fornire.
In base alle considerazioni fatte, per ottenere il segnale di livello è sufficiente integrare nel tempo
lo scarto Δf moltiplicato per una costante (piccola in modo che il segnale di livello assuma un valore
consistente dopo il termine della regolazione primaria) e cambiato di segno (infatti a seguito di un
disturbo che provoca una diminuzione della frequenza in rete, quindi Δf negativo, il segnale di livello
deve essere positivo in modo da aumentare la potenza resa dai gruppi di generazione). Ne consegue
che lo schema a blocchi di fig. 2.1 si modifica e risulta:
22
Nello schema le funzioni dei blocchi G1 è GB dipendono dal regolatore di velocità e sono i due
membri che compongono la 7.10, G2 è la funzione di trasferimento presentata nel precedente
paragrafo.
Con Simulink possiamo simulare il conportamento del sistema a seguito di un disturbo.
Dai risultati della simulazione notiamo come nei primi venti secondi dopo il disturbo (che avviene
a cinque secondi dall’inizio) è la regolazione primaria ad intervenire per contenere la diminuzione
della frequenza. Una volta esaurita la regolazione primaria (a circa venti secondi), la regolazione
secondaria agisce sul distributore della turbina in modo da aumentare la potenza generata e quindi
ristabilire la frequenza nominale dopo un tempo che è di circa un centinaio di secondi.
Fin’ora abbiamo trattato il semplice caso di una sola unità di produzione. Nella realtà però i
generatori elettrici sono molteplici e tutti cooperano per il mantenimento della frequenza. Gli
andamenti delle curve della frequenza a seguito del disturbo non cambiano anche in caso di più
generatori.
Nel caso di più generatori avremo che il contributo totale ΔPi’ di fig 2.4 sarà dato dalla somma
della potenza fornita da ogni unità di produzione in regolazione primaria(l’entità di ogni contributo
23
dipende dalle caratteristiche dell’unità e del suo regolatore). Il segnale di livello sarà sempre unico,
ma sarà inviato a tutte le unità che effettuano la regolazione secondaria (solo alcuni tipi di centrali
svolgono questa funzione). Anche il contributo di regolazione secondaria sarà quindi dato dalla
somma di tutti i relativi contributi.
24
Capitolo 3
3.1 Premessa
La regolazione di velocità di una qualsiasi centrale di produzione è necessaria ai fini di mantenere
la frequenza di rete costante. Esistendo una stretta correlazione tra la velocità di rotazione e la
frequenza della tensione prodotta da un generatore asincrono
60 1
2 (3.1)
25
Velocità di rotazione
coppie polari
[g/min]
1 3000
2 1500
3 1000
4 750
5 600
Il problema della regolazione della velocità si può studiare direttamente controllando la frequenza
della rete che è la grandezza che alla fine si vuole mantenere costante.
Pe
Figura 3.1 Schema a blocchi princio della regolazione della velocità di un gruppo
Il blocco che riguarda l’inerzia del gruppo è descritto dall’equazione del bilancio meccanico, cioè
dall’equazione differenziale:
1>
1:;< i :;< j :;<
Q
(3.2)
Dove la frequenza è espressa in Hz, le potenza sono in p.u. della potenza apparente nominale
dell’alternatore, Ta è il tempo di avviamento in secondi del gruppo turbina alternatore riferito alla
26
velocità e alla potenza apparente nominale. Considerando che nella regolazione della frequenza si fa
riferimento alla potenza attiva nominale Pn, la definizione del tempo di avviamento assume la forma:
PΩ%>
Q ≜
5>
(3.3)
Il valore del tempo di avviamento per un gruppo turbina alternatore si aggira intorno alla decina di
secondi e l’inerzia M del gruppo è così definita:
Q 5>
k
Ω>
(3.4)
1>
Gr:s<
Δfrif Δεf Δθ ΔA ΔPm Δf
OQ 5>
Gv(s) Ga(s)
Gf:s<
Gc:s<
Nello schema i blocchi rappresentano in termini di trasformate di Laplace quelli dello schema di
fig. 3.1. Le grandezze Δf, ΔPm, ΔPc, rappresentano le trasformate di Laplace della variazione di
frequenza, della variazione di potenza meccanica o regolante e delle variazioni di carico, cioè potenza
attiva generata. Quest’ultima grandezza è quindi il disturbo o perturbazione che si manifesta in rete e a
seguito del quale occorre regolare la frequenza per riportarla al valore desiderato.
La funzione Δf si può quindi esprimere come segue:
1 1>
∆1 = ∆5i − ∆5j − b* :O<∆1
O Q 5>
(3.5)
La funzione di trasferimento b* :O< tiene conto della variazione dei carichi in funzione della
frequenza della rete e solitamente viene approssimata con un puro guadagno. Per simulare la presenza
di motori elettrici in rete si può anche aumentare corrispondentemente l’inerzia del gruppo.
Le restanti funzioni di trasferimento presenti nella figura 3.2 si esprimono come segue
∆PY
b :O< ≜
∆m (3.6)
27
∆m
bq :O< ≜
∆r
∆r
bY :O< ≜
∆gG
bG :O< ≜ b :O<bq :O<bY :O<
Ovviamente le funzioni di trasferimento b :O<, bq :O< dipendono dal tipo di impianto e dal punto
in cui si effettua la linearizzazione mentre b* :O< dipende dalla variazione dei carichi in funzione della
frequenza come già visto. Gli effetti di b* :O< vengono successivamente trascurati perché poco
influenti e per semplificare l’individuazione delle caratteristiche fondamentale richieste alla bG :O< e di
conseguenza alla bY :O<, cioè al regolatore di velocità.
Il nostro obiettivo è quello di ottenere a regime un errore di frequenza nullo (∆gG = 0) in ogni
possibile condizione di regime. Per ottenere un tale comportamento è necessario supporre che ci sia
un effetto integrale nella bG :O<, ovvero nel regolatore. Una bG :O< puramente integrale non sarebbe
accettabile in quanto, i poli dell’equazione caratteristica risulterebbero immaginari coniugati. Con tale
caratteristica il sistema risulterebbe a smorzamento nullo, con conseguente continua oscillazione del
valore di regime della frequenza in seguito a un disturbo.
Risulta essere adatta al nostro scopo una funzione bG :O< del tipo proporzionale-integrale, cioè:
1 + OQ%
bG :O< = HG
O
(3.7)
Dove:
1 HG Q%
y D=−
w 2 k
x HG HG Q% %
- = c − E F
(3.11)
w U
v k 2k
Ovvero
̅ = −ℒ-> ± -> e1 − ℒ %
s?,% (3.12)
Con
28
y H
w-> c G -U{
k
x 1 HG
(3.13)
w ℒ Q% c
v 2 k
I poli ottenuti con una funzione di trasferimento proporzionale integrale risultano quindi spostati verso
l’asse reale negativo, lungo un cerchio avente il centro nell’origine del piano e raggio pari a -> . Lo
smorzamento risulta positivo e quindi il ciclo di regolazione della frequenza risulta stabile
Rinunciando ad avere un errore nullo a regime, si può assumere come funzione di trasferimento
Gf(s) un’equazione del tipo:
1 + OQ%
bG :O< = HG Q?
1 + OQ?
(3.14)
Ponendo T1 abbastanza grande si riesce ad avere un errore Δεf abbastanza piccolo a regime, ovvero
Δf piccolo se si suppone nulla la variazione del riferimento di frequenza. Questo si vede in particolare
se poniamo T1 tendente all’infinito nella 3.14, la quale con questa ipotesi tende alla 3.7.
Ponendo nella 3.14 s=0, si può valutare il guadagno statico che risulta essere pari a
∆5Y ∆5Y
HG Q? = bG :0< = :0< = |− :0<}
∆gG ∆1 ∆G ~•€ •U
(3.15)
Risulta pari al rapporto in condizioni di regime, tra la generica variazione di potenza regolante
(ovvero di carico) e quella corrispondente di frequenza, cambiata di segno, grandezza che viene
chiamata energia regolante permanente del gruppo, indicata generalmente con Ep e misurata in
[MW/Hz].
Viene definito anche lo statismo permanente di regolazione, indicato con bp e calcolabile con la
seguente relazione:
−∆1/1> 5> 5>
‚ =| :0<} = =
∆5Y /5> ∆G
1> HG Q? 1> bG :0<
~•€ •U
(3.16)
−∆1/1> 5>
‚ =| :∞<} =
∆5Y /5> ∆G
1> HG Q%
~•€ •U
(3.19)
29
Con queste definizioni si può scrivere la 3.15 nelle seguenti forme alternative:
1 OQ% , ,„
bG :O< , ,„
1 OQ? 1 OQ?
(3.21)
1‰ 1‰
5> 1 1 OQ% 5> 1 ‚ ‚„
bG :O< † ∗ ‡ ˆ ∗ Š
(3.22)
Possiamo ora analizzare l’andamento qualitativo della 3.21 mediante la costruzione del diagramma
di Bode, considerando come detto che T1>T2. Il diagramma ricavato è il seguente:
Dal diagramma possiamo notare come al verificarsi del disturbo in linea (aumento o calo di potenza
richiesta), l’energia regolante è modesta e pari a quella transitoria; col svilupparsi del transitorio di
regolazione questa energia aumenta progressivamente fino a raggiungere il valore dell’energia
regolante permanente nelle nuove condizioni di regime permanente. Se avessimo mantenuto la
funzione di trasferimento 3.7 avremmo trovato come energia regolante permanente un valore infinito,
a cui corrisponde uno statismo nullo.
In conclusione possiamo ottenere il legame dinamico tra variazioni di potenza meccanica e variazioni
di frequenza, che sarà del tipo:
5Y 1 OQ%
,
1 1 OQ?
(3.23)
5Y ⁄5> 11 OQ%
1⁄1> ‚ 1 OQ?
(3.24)
30
Figura 3. 4 Caratteristica statica frequenza-potenza
L’angolo α è funzione dell’energia regolante permanente. Tanto più quest’ultima è grande (piccola)
tanto più l’angolo α tende a 90° (0°).
La caratteristica statica si può tracciare ovviamente in termini relativi, risulta sempre rettilinea come si
vede in figura:
31
La Gf(s) ottenuta nel corrente paragrafo ha carattere del tutto generale e prescinde dal tipo di impianto
di generazione (idroelettrico, nucleare, termico). Per progettare quindi il regolatore di velocità, ovvero
la sintesi di Gr(s), occorre conoscere Ga(s), cioè la funzione di trasferimento del sistema di adduzione
e della turbina, al fine di rispettare i requisiti prima visti della Gf(s).
32
Capitolo 4
Lo schema mostra la galleria in pressione che preleva l’acqua dal serbatoio o fiume e la indirizza
verso la condotta forzata che va in centrale alla turbina. Nel mezzo è presente il pozzo piezometrico
33
necessario per smorzare le onde di sovrappressione che si presentano durante transitori idraulici come
parzializzazione o chiusura del distributore della turbina che danno origine al colpo d’ariete.
Nel nostro caso non è necessario studiare i transitori che coinvolgono tutto il sistema, ma è
sufficiente prendere in considerazione la sola condotta forzata (gli effetti dell’impianto a monte della
condotta forzata sono modesti). Di conseguenza lo schema semplificato che analizzeremo sarà il
seguente:
Con le ipotesi fatte il fluido si può considerare perfetto e quindi per il nostro studio è sufficiente l’uso
della nota equazione di Bernoulli:
%
• ’“O; 2;!
_
• 2‘
(4.1)
Con g accelerazione di gravità, z altezza geometrica (energia di posizione), ps/γ altezza piezometrica
(energia di pressione) e v2/2g altezza cinetica (energia di velocità).
Supponiamo ora che l’impianto sia ad alta caduta e quindi possiamo trascurare l’energia di
velocità, piccola rispetto alle altre grandezze della formula. L’equazione 4.1 assume allora la forma
seguente se applicata tra la quota del pelo libero dell’acqua nel serbatoio e la sezione di imbocco della
condotta:
34
”S 0 sin Ž
_?
•
(4.2)
”≜
_%
•
−
La differenza
”S − ” = + 0 sin Ž
_? _%
•
(4.3)
Rappresenta una forza accelerante (Ht-H)γAc agente sulla colonna d’acqua, avente massa pari a
γLAc/g e velocità Q/Ac, contenuta nella condotta forzata
•0m* :) ⁄m* <
:”S − ”<•m* = ∗
‘ ;
(4.4)
La 4.4 diviene
”S − ”
)=
P*
(4.6)
L’energia H allo sbocco della condotta dipende quindi dalla portata Q. se trascuriamo le variazioni
di salto totale, dalla 4.4 si deduce:
” = −•– :O< ) (4.7)
Dove
0
•– :O< ≜ OP* = O
‘m*
(4.8)
Che in analogia ai fenomeni elettrici viene chiamata impedenza della condotta (H corrisponde alla
tensione, Q alla corrente).
Un’altra analogia tra i due fenomeni si ritrova nella propagazione delle onde di tensione nelle linee
di trasmissione con la propagazione delle onde elastiche in un fluido. Per ottenere relazioni analoghe è
necessario però tener conto dell’elasticità e degli attriti del sistema idraulico. Per tener conto di questi
fenomeni, si assume come impedenza d’onda della condotta la seguente definizione:
•– :O< ≜ •—* tanh r* (4.9)
dove le variabili sono così definite:
y ’ 2H)— ‘m*
w • = c1 +
—*
‘m* O0
x O0 2H)— ‘m*
(4.10)
w r* = c1 +
v ’ O0
dove K è il coefficiente delle perdite di attrito, proporzionali al quadrato della portata, Q0 è il valore
della portata allo sbocco della condotta e c è la velocità di propagazione delle perturbazioni elastiche
nella condotta (circa 1000 m/sec).
35
Trascurando i fenomeni elastici (c=∞) e gli attriti (K=0) ed elaborando le 4.10 si ritorna ad avere
come impedenza d’onda l’equazione 4.8.
4.2 Turbine[2]
Le turbine idrauliche sono il motore che converte l’energia idraulica in energia meccanica, che
sarà poi riconvertita in energia elettrica.
La potenza meccanica ottenibile in uscita dalla turbina dipende dunque dalla potenza idraulica in
ingresso. Tale potenza deve poter essere regolabile, e per questo è presente prima dell’ingresso
dell’acqua in turbina una valvola di ammissione, detta anche distributore (come mostrato nello
schema di figura 4.3). Agendo su questa valvola si va a modificare la sezione di passaggio dell’acqua
ottenendo come effetto la riduzione della portata d’acqua.
La relazione dinamica che ci interessa conoscere per i nostri studi, è quella che intercorre tra
l’apertura A della valvola e la potenza meccanica Pm sviluppata all’asse della turbina.
A seconda delle turbine presenti nel singolo impianto, la valvola di regolazione, o distributore,
varia a seconda della tipologia di turbina. Le turbine utilizzate maggiormente negli impianti idraulici
sono le Pelton, Francis e Kaplan. La scelta di una turbina rispetto ad un’altra è una conseguenza delle
caratteristiche dell’impianto di produzione. In linea di massima, per impianti che hanno un elevato
salto ma portata limitata, le Pelton sono le più indicate; per impianto con salto e portate medie, sono
da preferire le turbine Francis, mentre per salti modesti ma portate elevate sono più indicate le turbine
Kaplan. In figura 4.4 è presente un diagramma che mostra quanto appena enunciato.
36
Figura 4. 4 Scelta della turbina in funzione di salto e portata disponibile. [14]
A seconda della turbina, varia anche la valvola di ammissione, o distributore. Per le turbine Pelton,
la regolazione della portata viene ottenuta tramite la valvola a spina. Questa valvola è di forma
cilindrica e varia la sezione di efflusso con uno spostamento assiale in modo tale da non disturbare il
moto dell’acqua e minimizzare le perdite di carico.
Figura 4. 5 Regolazione della portata per una turbina Pelton. Fonte www.energialternativa.org
Le turbine Pelton possono avere diversi ugelli dai quali l’acqua esce e viene indirizzata sui
cucchiai. Solitamente non si trovano turbine con un numero maggiore a sei ugelli. Questo per far si di
ottenere una caratteristica di regolazione migliore (in termini di rendimento), in quanto per ottenere
una bassa potenza, è preferibile diminuire il numero di ugelli. Questo si vede bene nella figura
sottostante che mostra la caratteristica potenza all’asse – rendimento.
37
Figura 4. 6 Caratteristica potenza-rendimento di una turbina Pelton. [14]
Le turbine Francis sono di costruzione più complessa e gli organi di regolazione principale sono le
pale del distributore, e in alcune turbine anche le pale della girante sono ad inclinazione regolabile. È
possibile quindi avere due organi distinti su cui agire. In ingresso alla turbina Francis troviamo una
chiocciola che distribuisce la portata su tutta la circonferenza, di seguito troviamo il distributore a pale
mobili e quindi la girante (eventualmente anche quest’ultima a pale mobili); l’acqua poi abbandona la
turbina attraverso il diffusore. In figura si riesce a notare la chiocciola e colorato in arancione il
sistema di regolazione dell’inclinazione delle pale del distributore. Si vede come agendo sulla staffa
circola si varia contemporaneamente l’inclinazione di tutte le pale.
38
Figura 4. 8 Turbina Francis con sistema di regolazione del distributore.
Fonte italian.hydropower-turbine.com
Nelle due figure sottostanti sono mostrati una pala girevole e il movimento che si ottiene agendo
sul sistema di regolazione.
39
La caratteristica di regolazione di questa turbina in funzione della portata è mostrata nella figura
seguente:
Figura 4. 11 Caratteristica generica di regolazione di una turbina Francis. Fonte sito www.cap2.it
Si nota come la portata e la potenza varino quasi linearmente all’aumentare della sezione
disponibile per l’acqua, mentre il rendimento cala velocemente quando la sezione di efflusso scende al
di sotto del 50%.
Le turbine Kaplan sono macchine assiali a reazione con un limitato numero di pale (fino a sei)
disposte assialmente sul mozzo che interagiscono con un flusso il cui moto elicoidale, senza
componente radiale della velocità, è imposto dal distributore; sia il distributore che la girante possono
essere a pale fisse (turbine a elica) o orientabili in modo da avere un buon rendimento anche ai carichi
parziali. Sono caratterizzate da un elevato numero di giri caratteristico che supera di norma quello
delle Francis.
40
Tornando alla nostra analisi, per una turbina Pelton, la potenza entrante nella turbina è quella
dovuta all’energia cinetica del getto, in uscita dalla valvola, che investe le pale della girante della
turbina e vale
56 •)” (4.10)
Dove H e Q hanno le seguenti espressioni:
)%
”
2‘m%
(4.11)
) = me2‘” (4.12)
Essendo
5i∗ = •™ ∗ ) ∗ ” ∗ (4.16)
Come abbiamo visto il rendimento della turbina è funzione del grado di apertura A, sarà quindi:
™ m
= š›œ ∗
™∗ m
(4.17)
Dove
m∗ ™ ∗
š›œ ≜ E F
™∗ m
(4.18)
Che come si vede dipende dalla pendenza della funzione che lega il rendimento all’apertura della
valvola. Sostituendo nella 4.14 la 4.16 otteniamo
5i ∆) ∆” m
= ∗ + ∗ + š›œ ∗
5i∗ ) ” m
(4.19)
41
•) ∗ ∗ )∗
E F = J:2‘”< ‰% M = ∗
?
•m m
(4.21)
•) ∗ 1 ∗
1 )∗
E F = | m:2‘”<& ‰% 2‘} =
?
•” 2 2 ”∗
(4.22)
Abbiamo quindi ricavato con la 4.18 e la 4.22 il comportamento linearizzato della turbina nell’intorno
di un punto di funzionamento a regime. Andiamo ora ad analizzare la funzione di trasferimento del
complesso sistema di adduzione e turbina.
Questa equazione lega insieme le variazioni di energia allo sbocco della condotta con le variazioni
dell’apertura valvole. Sostituendo la 4.22 nella 4.18 ricaviamo
5i m 3 ”
= ž1 + š›œ Ÿ ∗ +
5i∗ m 2 ”∗
(4.26)
Utilizzando anche la 4.24 possiamo trovare la funzione di trasferimento cercata che, posto
š œ ≜ 1 + š›œ (4.27)
Risulta essere:
O 1 )∗
1−E − 1F • :O<
5i š œ 2 ”∗ – m
=š
5i∗ œ
1 )∗ m∗
1+ • :O<
(4.28)
2 ”∗ –
Quindi risulta essere
5i = b :O< m (4.29)
Con
)∗
” ∗ •– :O<
1−‚
b :O< = H
1) ∗
1+ • :O<
(4.30)
2 ”∗ –
42
5i∗ 5i∗
yH H =ž1 + š›œ Ÿ
m∗ œ
m∗
x 1 3
‚= E − 1F
(4.31)
v 2 Hœ
Analizzando la relazione ingresso uscita si nota che la risposta del sistema dipende dal punto di
funzionamento (soprattutto dall’andamento della caratteristica rendimento, apertura), e anche
dall’impedenza della condotta che può variare se si considerano o no i fenomeni elastici ed attriti.
Se trascuriamo l’elasticità della condotta e gli attriti, la Ga assume la forma seguente:
1 − ‚OQ–
b :O< = H
Q
1+O –
2
(4.32)
Anche detto tempo di avviamento della condotta. Si nota subito che questo tempo è proporzionale
alla portata secondo l’equazione:
)∗
Q– ≜ Q–>
”>
(4.34)
¢ ¢
e sfruttando la proprietà:
tanh # - ( = tan #- (
2 2
(4.38)
¢
la Ga(jw) diviene
1 − 2‚ £ tan #- (
b :O< = H 2
¢
1 + ‚£ tan #- (
2
(4.36)
Avendo quindi ora a disposizioni le funzioni di trasferimento possiamo ricavare il diagramma di Bode
e confrontare l’andamento nei due casi visti (uno senza attriti e fenomeni elastici, l’altro senza attriti
43
ma con fenomeni elastici). Dall’analisi fatta , risultano due andamento molto diversi da come si può
vedere in figura 4.12. I valori usati per ricavare le curve sono i seguenti:
b=1 μ=1
Tω=1 sec τ=1,57 sec
Ka=10 0,1<ω<9 rad/s
Figura 4. 13 Diagramma di Bode della funzione di trasferimento del sistema turbina e sistema di
ammissione
Si nota come i fenomeni elastici influiscano in maniera importante dopo una certa pulsazione, e come
la stessa fase cali maggiormente rispetto al primo caso semplificato. In particolare, la funzione di
trasferimento presenta delle frequenze di risonanza e anti risonanza e la fase cala velocemente verso -
∞ al crescere della frequenza. Questo costituisce un limite nella rapidità del ciclo di regolazione della
frequenza in quanto è pericoloso per la stabilità avvicinarsi con la banda di regolazione alla prima
frequenza di risonanza (nel caso in esempio uguale a 2 rad/s).
Possiamo comunque dire che per frequenze molto minori della prima frequenza di risonanza, le
due funzioni sono praticamente identiche. Ne consegue che per basse frequenze si può usare
l’approssimazione seguente per la Ga(jω):
1 ‚-Q–
b :O< H
Q
1 - –
2
(4.37)
44
1 − ‚? OQ–
b :O< H
1 + ‚% OQ–
(4.38)
avendo assunto come grandezze base la potenza nominale e l’apertura valvola nominale. Il guadagno
statico, cioè il valore di ga con s=0, è per definizione la pendenza della caratteristica statica nel punto
di funzionamento ed è espressa dalla:
:5i ⁄5> <
∗
‘ :0< ≜ t ¤
:m⁄m> <
(4.42)
Il valore del guadagno statico può essere ricavato dalle caratteristiche statiche di potenza-apertura
valvole per le diverse turbine. In figura 4.13 sono mostrati gli andamenti qualitativi di queste curve
per le diverse turbine.
0,8
Pn/Pm
0,6
Pelton
0,4 Francis
Kaplan
0,2
0
0 0,5 1
a(p.u.)
45
Tramite queste curve è possibile ricavare il valore di Ka. Per turbine Pelton, il suo valore va
crescendo e parte dal valore di 0,3 per arrivare fino a circa 2. Per turbine Francis invece il suo valore è
circa costante e pari a 1.
46
Capitolo 5
47
Come mostrato in figura 5.1, le fasi che compongono il ciclo sono le seguenti:
- riscaldamento dell’acqua a pressione costante sia nel preriscaldatore sia in caldaia;
- vaporizzazione a pressione e temperatura costante in caldaia;
- surriscaldamento del vapore, praticamente isobaro a pressione pS in caldaia;
- caduta di pressione, in parte sulla valvola di ammissione, in parte nello stadio di alta pressione
della turbina;
- risurriscaldamento del vapore, praticamente isobaro, a pressione pR in caldaia;
- espansione del vapore nello stadio di bassa pressione dalla turbina;
- trasformazione isoterma ed isobara (pressione p0) del vapore esausto in acqua nel
condensatore.
Senza entrare nei particolari, spiegati nel dettaglio in, lo schema a blocchi linearizzato della
caldaia di un impianto termoelettrico è il seguente:
48
È di nostro interessa conoscere la relazione che lega l’apertura della valvola A di ammissione del
vapore in turbina alla portata di vapore surriscaldato qs e alla sua pressione ps. Questa relazione è
particolarmente complessa perché deve tener conto del peso specifico del vapore e del salto di
pressione dovuto alle valvole. Possono essere di aiuto i grafici seguenti che mostrano le caratteristiche
statiche portata-apertura per i casi di ammissione ad arco parziale o totale (nell’ammissione ad arco
parziale le valvole non si aprono in modo sincronizzato, ma in modo sequenziale).
Vediamo dalle figure che la caratteristica è fortemente non lineare con variazioni di pendenza
molto elevate. Tali punti sono quelli in cui una valvola è completamente aperta e quella successiva
comincia ad aprirsi. Tale linearità è comunque compensata introducendo una non linearità nel sistema
di comando delle valvole tra l’apertura del servomotore e l’apertura delle valvole di ammissione.
49
Dalle caratteristiche è facile dedurre la legge di dipendenza di qs da ps e da a. Si ha.
_ 1: _ , < (5.1)
linearizzando si ricava:
_ =š _ +š (5.2)
con
y š ≜ E• _ F
∗
w • _
x • _ ∗
wš ≜ E F
(5.3)
v •
Per ricavare la potenza meccanica in uscita dalla turbina si deve effettuare il bilancio di energia per
tutti gli stadi della turbina, cioè l’insieme di quelli di alta e bassa pressione.
Tralasciando i calcoli intermedi, si ricava la relazione che esprime pm come:
Φℎ¥
= + = |ℎœ + }
i œ ¥ œ
1 + OQa
(5.4)
Seguendo il procedimento mostrato in Capitolo 9.4, si possono ricavare i contributi dello stadio AP e
BP, e si arriva ad esprimere la variazione di potenza meccanica con la seguente espressione:
Ž
= |D + }
i œ
1 + OQa
(5.7)
Solitamente α=0.3÷0.35, di conseguenza β=0.65÷0.7. È notare che la portata qAP del vapore negli
stadi di alta pressione è dinamicamente diversa dalla portata qS del vapore surriscaldato a causa
dell’accumulo di massa del vapore che si ha nella turbina stessa. Si tiene conto di ciò introducendo
50
una costante di tempo Tsc1 (del valore di 0.1÷0.3 secondi)di ritardo tra qs e qap’, ovvero, traducendo il
tutto in trasformate di Laplace:
1
œ
1 + OQ_*? _ (5.8)
La stessa cosa avviene negli stadi di media e bassa pressione della turbina. Si aggiunge quindi
un’ulteriore costante di tempo Tsc2 che varia tra 0.4÷0.6 secondi. in questo modo la 5.6 diviene:
Ž
= |D + }
i
:1 + OQa <:1 + OQ_*% < œ (5.9)
In conclusione si ha quindi:
=D œ œ
© Ž
=
¥
:1 + OQa <:1 + OQ_*% < œ
(5.10)
con
Q_*? + Q_*%
Qª« =
2
(5.13)
Si perviene ora facilmente allo schema a blocchi linearizzato della turbina e della caldaia di un
impianto termoelettrico.
51
Figura 5. 5 Schema a blocchi linearizzato della turbina e della caldaia di un impianto termoelettrico
i ‘ :O<
(5.14)
1 1 ODQa
‘ :O< š
:1 OQ_* < 1 OQa
(5.15)
%
in cui ka è il guadagno statico ed è pari alla pendenza della caratteristica in figura 5.4 ed è quindi
fortemente dipendente dal punto di funzionamento.
Se invece trascuriamo le variazioni Δpv della pressione del vapore in caldaia, ma ammettendo che
Δps possa variare (a causa di Δqs), si ricava dallo schema a blocchi:
_ š š š _
(5.16)
_ š{
(5.17)
dove
š
š{
1 š š
(5.18)
Con tali ipotesi, valgono ancora le 5.14 e 5.15, salvo sostituire ka con ka’. Tale guadagno, viste le
definizioni 5.3, è la derivata di qs rispetto a ps ad a costante, nel punto di funzionamento. Dalle
caratteristiche di figura 5.4, kp=1 p.u. / p.u. mentre k è al massimo 0.1÷0.2 p.u. / p.u. Si può ricavare
quindi che:
š{ •š
(5.19)
dove 0.85 ≤ γ ≤ 1. Risulta quindi che la funzione di trasferimento avrà un guadagno statico diminuito
al massimo del 15% rispetto al caso con ps costante.
52
Capitolo 6
Regolatori di velocità
m ∗
Hq ≜ E F
r
(6.1)
Hq
bq :O< ≜
1 + OQq
(6.2)
53
- la velocità del servoposizionatore e quindi quella del servomotore dθ/dt, devono essere
mantenute entro certi limiti in modo da non causare elevate sollecitazioni di pressione
nell’impianto idraulico
Supponiamo ora che sia in prima approssimazione
bq :O< ≜ Hq
(6.3)
e che la funzione di trasferimento della turbina e del sistema di adduzione sia anch’essa
assimilabile a un puro guadagno (vedi capitolo 4), cioè:
bœ :O< ≜ H
(6.4)
r> Hq H
‚{ ≜ ‚
5>
(6.7)
r> Hq H
‚„{ ≜ ‚„
5>
(6.8)
la 6.6 diviene:
una tale funzione di trasferimento si può ottenere con una delle due seguenti soluzioni:
- regolatore accelerometrico, mostrato in figura 6.1,
- regolatore a retroazione transitoria (o asservimento cedevole), mostrato in figura 6.2.
Si devono inoltre prendere in considerazione le nuove costanti di tempo dovute ai ritardi del
trasduttore di velocità Tb e la costante di tempo dell’accelerometro Tac.
54
Nel primo caso, figura 6.1, ricaviamo la funzione di trasferimento gr(s) pari a:
1 OQ * 1 1 1 OQ * 1
‘Y :O<
1 OQ¬ š O⁄šœ šª š 1 OQ¬ 1 O⁄šœ šª š
(6.11)
Considerando che Tb è solitamente di piccolo valore, concludiamo che la 6.11 è equivalente alla
6.10 con:
‚{ š
y
1
Q?
x šœ šª š
(6.12)
v Q% Q *
šœ šª :1 OQ- <
‘Y :O<
šœ šª š O1 šœ šª :š š S <Q- O % Q-
(6.13)
55
š š š ≫ 1⁄Q-
® œ ª S
šS ≫ š
(6.14)
si ricava
1 + OQ- 1
‘Y :O< =
š + O:š + š S <Q- 1 + OQ°
(6.15)
in cui T3 è dello stesso ordine di grandezza di Tb e può quindi essere trascurata. Si ottiene quindi
che gr(s) è della forma cercata, ponendo:
‚{ = š
y
š + šS
Q? = Q-
x š
(6.16)
v Q% = Q-
e
‚„{ = š + š S > š S
(6.17)
Dalla teoria dei controlli automatici sappiamo che il comportamento statico (a regime) di un ciclo
di regolazione è definito dalla precisione a regime; per avere una buona precisione si deve avere un
guadagno elevato. Nel nostro caso quindi questo sta ad indicare che lo statismo permanente dovrà
essere sufficientemente basso, cioè assumere un valore che sia:
56
Figura 6. 3 Schema a blocchi ciclo di regolazione della frequenza
11 OQ% 1
‚ 1 OQ? ‚„
(6.19)
si ricava per
1
ω±
Q%
(6.20)
1
|g:jω<|
© -‚„ Q
¸‘: -<
(6.21)
180° tan&? :-Q% <
57
Dal diagramma tracciato è evidente che per ottenere un buon margine di fase e un ciclo di
regolazione che non sia troppo lento, occorre scegliere ωt tra le pulsazioni che soddisfino la 6.20.
Una volta scelta ωt dalle 6.21 si ricava
1
b„
© -„ Q (6.22)
•„ = ; 2&? :-„ Q% <
Se per esempio assumiamo ωt=0.5 rad/sec e si desidera γt=60°, poiche Ta è dell’ordine di 8 sec si
ricavano rispettivamente bt=0.25 p.u./p.u. e T2=3sec. Se inoltre scegliamo bp=0.05 p.u./p.u. risulta
T1=(bt/bp)*T2 = 15sec.
Dobbiamo tenere conto prima di concludere, che il guadagno statico ka della funzione di
trasferimento della turbina e del sistema di adduzione dipende dal punto di funzionamento, come il
guadagno kv del sistema di comando valvole. Il guadagno 1/bp è proporzionale ad entrambi. Se
supponiamo che ka aumenta come nel caso delle turbine pelton (causa aumento del carico), il
diagramma di Bode trasla verso l’alto e la pulsazione di taglio aumenta avvicinandosi a quella di
risonanza della condotta. Di conseguenza il ciclo può diventare instabile anche se era stabile a bassi
carichi. Per ovviare a questo problema viene introdotta nei regolatori una non linearità di
compensazione della non linearità potenza-apertura. In alternativa si aumenta lo statismo permanente
bp’ di apertura all’aumentare del carico in modo tale che bp resti costante: ciò significa rendere il
regolatore adattativo alle varie condizioni di funzionamento dell’impianto.
1 1 + ODQa
‘ :O< = š
:1 + OQ_* < 1 + OQa
%
Se trascuriamo in questa relazione la costante di tempo Tsc capiamo subito che la desiderata forma
per la funzione gf(s) è già realizzata dalla ga(s) a men o di una costante moltiplicativa. Sarà quindi
richiesto al regolatore di velocità soltanto un’azione puramente proporzionale.
Supponiamo che valga, come ipotizzato per i gruppi idroelettrici, la 6.3, e indicando con 1/b il
guadagno incognito della gr(s), dalla
in valori relativi si ha
58
1 1 1
‚ ‚ šq š
(6.24)
T? Qa
®
T% D Qa
(6.25)
Ricordiamo inoltre che Tsc è dell’ordine dei 0,5 secondi, e quindi 1/Tsc è circa 2 rad/sec e
paragonabile alla prima frequenza di risonanza della condotta, e che tale costante di tempo è relativa
ad un polo doppio; si comprende facilmente come valgano anche in questo caso, per la scelta della
pulsazione di taglio ωt, le stesse considerazioni fatte nel paragrafo precedente.
La stessa cosa possiamo dirla per quanto riguarda il guadagno 1/b e la sua dipendenza da kv e ka,
che dipendono dal punto di funzionamento.
In base a quanto detto esiste quindi una forte analogia tra la funzione di trasferimento complessiva
di un gruppo termoelettrico e quella di un gruppo idroelettrico, anche in termini di costanti di tempo
T1 e T2 (risulta infatti TR pari a 10-15 sec, α pari a 0.3-0.35).
Si è prima detto che il regolatore deve avere un’azione puramente proporzionale con guadagno
definito dalla 6.24. Se si suppone che siano completamente compensate le non linearità che rendono kv
e ka variabili con il punto di funzionamento, si può ritenere che sia in ogni condizione di
funzionamento
šq š 1 . ´.⁄ . ´. (6.26)
e di conseguenza
1 1
‚ ‚
(6.27)
Nella realtà però occorre sempre un servomotore di regolazione per comandare l’apertura delle
valvole; per questo lo schema del regolatore di velocità è del tipo illustrato in figura 6.5 nella quale
1/Ts rappresenta il guadagno del servomotore dove Ts è dell’ordine di 0.02 sec.
59
Spesso esiste insieme al ciclo di regolazione della velocità, un ciclo per la regolazione della
potenza elettrica erogata, che a regime coincide con la potenza meccanica. La figura 6.6 rappresenta in
modo schematico questa situazione. Questo ciclo di regolazione è molto più lento di quello
riguardante la frequenza e i suoi obiettivi sono quelli di variare facilmente la potenza da produrre
tramite il regolatore di carico, introdurre dei limiti nel gradiente di presa di carico (per un buon
funzionamento dell’impianto), inviare alle regolazioni in caldaia i segnali adeguati per far fronte alla
richiesta di potenza. In tale ciclo è inserito il regolatore di potenza, di tipo proporzionale integrale per
garantire che, a regime, la potenza prodotta sia uguale a quella richiesta. La sua funzione di
trasferimento è del tipo:
1 OQ
‘Y :O<
OQ
(6.28)
Nella figura 6.6 sono raffigurati i due interruttori SW1 e SW2. Quando sono chiusi la regolazione
della frequenza avviene come nel modo indicato: tale modo di funzionamento è detto controllo
coordinato (della frequenza, della potenza e della caldaia). Agendo sugli interruttori citati si possono
ottenere diversi tipi di risposta a una perturbazione che avviene nel sistema. Se SW1 e SW2 vengono
aperti, si ottiene solamente una regolazione della frequenza e non della potenza, le regolazioni della
60
caldaia continuano a interagire con il ciclo di regolazione attraverso ps e qi. Tale modo di operare
viene chiamato caldaia segue nel senso che il ciclo di regolazione dominante è quello della frequenza
e la caldaia si deve adeguare a ciò che viene richiesto (le regolazioni di caldaia intervengono a seguito
della variazione di apertura delle valvole di ammissioni del vapore)
Un altro tipo di regolazione che si può ottenere è quello denominato turbina segue. In questo caso
la richiesta di potenza viene inviata alle regolazioni di caldaia (SW2 chiuso) che adegua la portata del
vapore al nuovo carico in modo da avere a regime l’apertura valvole uguale a quella precedente il
disturbo.
61
avviamento della condotta Tω dipendente dal punto di funzionamento iniziale, cioè dalle
condizioni di regime precedenti i transitori da esaminare. Si suppone cioè:
con turbine Pelton
1 ‚Q–>
o
∗
‘ : <
Q–> ∗
1
2
(6.29)
‘ : < 1 (6.31)
Nelle relazioni precedenti q* è la portata in p.u. nel punto di funzionamento, mentre b,b1 e b2
dipendono dal punto di funzionamento.
In base agli schemi precedenti e con le ultime ipotesi fatte possiamo costruire il diagramma a
blocchi della regolazioni nel caso di grandi variazioni, mostrato in figura 6.7:
Figura 6. 7 Schema a blocchi semplificato per grandi variazioni del regolatore di velocità. della turbina
e del sistema di adduzione di un impianto idroelettrico
I valori che assumono le diverse grandezze indicate sono solitamente simili ai seguenti:
- per regolatori accelerometrici
Per i gruppi termoelettrici vale il diagramma di figura 6.6 con tutte le considerazioni fatte nel
paragrafo precedente. Anche in questo caso si trascurano:
- la zona morta del tachimetro e i limiti di velocità del servomotore,
62
- la dinamica del comando valvole,
- le non linearità della caratteristica portata-apertura e del sistema del comando valvole, in
quanto entrambe conpensate.
Nel caso di funzionamento in caldaia segue e trascurando le variazioni della pressione del vapore
all’ammissione, ci si può ricondurre allo schema di figura in cui si è trascurata anche la dinamica,
relativamente veloce, del servomotore di regolazione, il che corrisponde ad assumere che la posizione
θ del servomotore coincida con la variabile di uscita dal blocco con guadagno 1/bp. Valori tipici dei
parametri dello schema a blocchi sono:
Si riporta inoltre lo schema semplificato relativo alla regolazione di velocità di gruppi con turbine
a gas in cui Tf e la costante di tempo del combustibile, Tm quella del motore, mentre α è il guadagno
transitorio del motore.
In questo caso i valori tipici delle grandezze sono i seguenti:
Figura 6. 8Schema a blocchi semplificato per grandi variazioni del regolatore di velocità. della turbina
e del sistema di adduzione di un gruppo termoelettrico
63
Figura 6. 9 Schema a blocchi semplificato per grandi variazioni per i gruppi con turbine a gas
64
Capitolo 7
7.1 Introduzione
Finora abbiamo visto soltanto le funzioni di trasferimento che i regolatori di velocità dovrebbero
avere per ottenere un ciclo di regolazione della frequenza che sia stabile e con errore a regime il più
possibile contenuto.
Nel seguente capitolo vedremo alcuni esempi di regolatori di velocità più comunemente utilizzati
negli impianti idroelettrici e termoelettrici
65
Figura 7. 1 Tipico regolatore di velocità per gruppi idroelettrici
Il regolatore di velocità mostrato ha come sensore di velocità un pendolo di Watt. Questo pendolo,
ruota alla stessa velocità dell’alternatore di cui deve regolare la velocità. Le due masse ruotano e
grazie alla forza centrifuga, al variare della velocità fanno variare la posizione del punto A. In
particolare, se la velocità aumenta, il punto A si muove verso l’alto, viceversa, se la velocità
diminuisce, il punto A si muove verso il basso. La molla limita e regola la variazione della posizione
del punto A in funzione della velocità. Il sistema che regola l’apertura delle valvole è invece costituito
da un insieme di leve e un sistema oleodinamico in cui l’olio è mantenuto in pressione da una pompa.
Vediamo ora il funzionamento qualitativo del sistema rappresentato.
Partendo da una certa condizione di regime, supponiamo che a causa di un disturbo, la velocità di
rotazione del gruppo di generazione aumenti. Supponiamo inoltre di mantenere per ora costante la
posizione del punto B. In seguito all’aumento della velocità, le due sfere del pendolo tendono ad
allontanarsi dal centro di rotazione e quindi a muovere il punto A nel senso indicato dalla freccia in
figura. Si muove di conseguenza tutta l’asta i cui estremi sono i punti A e B (ricordiamo che B è
mantenuto fisso). Con il suo movimento, l’asta provoca lo spostamento verso l’alto del punto C e di
conseguenza si aprono le luci 1 e 2 del cilindro. L’olio in pressione quindi fluisce dalla luce 1 nella
parte superiore del cilindro 3 e spinge l’olio presente nella parte inferiore a tornare nella vasca
attraverso la luce 2. Si ottiene quindi nel caso ipotizzato lo spostamento graduale verso il basso del
punto D che va di conseguenza a richiudere le luci 1 e 2 alla fine del transitorio (si hanno nel
66
transitorio delle oscillazioni attorno alla posizione di chiusura finale delle luci.). Abbiamo ottenuto
alla fine un movimento che va a chiudere parzialmente le valvole della turbina a seguito dell’aumento
della velocità (si riduce cioè la potenza prodotta dalla turbina).
Ipotizzando che il sistema sia lineare, possiamo calcolare il modello matematico attraverso le
trasformate di Laplace. Definiamo quindi le grandezze necessarie per questo studio: Δa, Δb, Δc, Δd
sono gli scostamenti a partire dalla posizione iniziale dei punti A, B, C, D; ΔA, ΔB, ΔC, ΔD sono le
rispettive trasformate. Indichiamo con ΔF e ΔP le trasformate degli incrementi di frequenza e potenza.
Considerando i versi delle frecce in figura, possiamo determinare le seguenti relazioni:
šG ∆1 m = šG ∆[ + O[š { (7.1)
’ = š- ∆ + š ∆1 − š¬ ∆‚ . = š- ∆À + š ∆[ + š¬ ∆Á (7.2)
∆.
= −š R ∆’ ; À = −š
O
(7.3)
Nella 7.3 è presente l’integrale in quanto il movimento del punto D è proporzionale alla portata
fluente attraverso le luci 1 e 2 che è proporzionale allo spostamento Δd.
Sostituendo nella 7.3 la 7.2 ed elaborando l’equazione otteniamo:
š¬
š #−∆[ + ∆Á(
š
À= O
š- +
(7.4)
š
ed essendo kd diverso da zero:
š
š- š¬
À= O E−∆[ + š ∆ÁF
1+
(7.5)
šš-
Supponendo proporzionale la relazione tra ΔD e ΔPi abbiamo di conseguenza:
56 = š6 ∆À (7.6)
Sostituendo si ha:
š6 š
š- š¬
56 = O E−∆[ + š ∆ÁF
1+
(7.7)
šš-
Introducendo ora le seguenti definizioni:
š6 š 1 š¬
Ha = Qa = H¥ =
š- šš- š
(7.8)
otteniamo:
Ha
56 = :−∆[ + HÁ ∆Á<
1 + O Qa
(7.9)
Questa è l’equazione che troviamo con la sola azione tachimetrica. Se come si fa comunemente
introduciamo anche un’azione accelerometrica del tipo ka’df/dt (migliora la risposta dinamica del
sistema), troviamo:
67
1 OQ? Ha
56 Ha : ∆[< :H ∆Á<
1 O Qa 1 O Qa Á
(7.10)
Analizziamo ora il comportamento del regolatore a regime a seguito di una variazione ΔF della
frequenza e con ΔB=0:
La 7.10 diventa:
1 O Q?
56 Ha : ∆[< b? : ∆[<
1 O Qa
(7.11)
Applicando il teorema del valore finale all’equazione 7.11 arriviamo ad avere per t tendente a
infinito:
lim ∆56 :;< Ha lim ∆1:;<
„→Å „→Å
(7.12)
Ciò sta ad indicare che per tutte le volte in cui Δf(t) assume un valore costante a regime, il valore
ΔPi assume il valore:
∆56 Ha ∆1 (7.13)
Se consideriamo unitario il rendimento del gruppo turbina alternatore allora la 7.13 rappresenta
anche la variazione di energia prodotta. Il fattore KR è chiamato energia regolante della macchina e il
suo valore dipende dai parametri del regolatore e dal parametro ki.
Considerando costanti questi parametri per ogni condizione di funzionamento della macchina,
possiamo disegnare la caratteristica statica rettilinea del gruppo di generazione:
La figura mostra come a regime permanente, ad una data potenza erogata corrisponde ad una
determinata frequenza.
Dalla caratteristica possiamo ricavare il valore dell’energia regolante della macchina tramite la
relazione:
68
5W 5W
Ha
1U 1? 1U 1?
1W
1W
(7.14)
Definendo nella formula il parametro σ=(f0-f1)/fN che viene chiamato grado di statismo del
regolatore, otteniamo:
5W
Ha
Æ1W
(7.15)
Per comprendere il significato di KR possiamo fare un piccolo esempio: solitamente lo statismo del
regolatore ha un valore che si aggira attorno al 5%; supponiamo che il gruppo di generazione abbia
una potenza di 100MW, troviamo quindi che KR vale:
100 k=
Ha 40
0.05 ∗ 50 ”•
(7.16)
Vediamo ora quale è l’effetto che ha il variagiri nella caratteristica statica del gruppo. Il variagiri è
quel dispositivo che consente la variazione della posizione del punto B di fig. 7.1 attraverso un’azione
manuale o con telecomando.
Ritornando alla relazione 7.10 e supponendo questa volta ΔF=0, troviamo:
Ha
56 H ∆Á
1 O Qa Á
(7.17)
Applicando anche in questo caso il teorema del valore finale possiamo calcolare la risposta a
regime a seguito di una variazione ΔB. In questo caso si ha, a regime:
56 Ha HÁ ∆‚ (7.18)
La variazione della posizione del punto B quindi comporta la variazione della potenza immessa in
turbina. Agendo sul variagiri si ottengono caratteristiche statiche che danno luogo alla stessa potenza
al variare della frequenza.
69
Anti trasformando la 7.17 possiamo ricavare qual è il suo andamento in funzione del tempo a
seguito di una sollecitazione a gradino. Il suo andamento sarà di tipo esponenziale con costante di
tempo TR. Sarà quindi:
.
&„‰S
56 Ha HÁ ∆‚ #1 ! Ç( (7.19)
Le caratteristiche del regolatore studiato sono in linea di massima valide anche per regolatori di altri
tipi di centrali (per esempio le termoelettriche) e anche per regolatori di più recente costruzione.
Figura 7. 4 Regolatore di velocità oleodinamico della centrale idroelettrica del Furlo lungo il fiume
Candigliano. Fonte www.educambiente.it
70
Figura 7. 5 Regolatore di velocità oleodinamico. Fonte http://www.educambiente.it/
71
7.2 Regolatore di velocità per impianto termoelettrico[4]
L’esempio che prendiamo in esame nel caso di impianto termoelettrico è presentato nel seguente
schema:
Il dispositivo che rileva la velocità di rotazione della turbina è l’impeller, o altrimenti detto anche
impulsore. Il suo funzionamento si basa sul fatto di convertire la rotazione dell’albero della turbina in
una pressione d’olio. questo dispositivo era costituito da un cilindro cavo coassiale all’albero della
turbina con fori di ingresso e di uscita calibrati, attraversato da olio in alta pressione proveniente dal
sistema alimentazione della turbina. Nel suo funzionamento e grazie alla particolare costruzione
consentiva di ottenere a turbina in rotazione una pressione dell’olio in uscita dall’impeller
praticamente proporzionale alla velocità di rotazione del gruppo di generazione.
72
Figura 7. 8 Vista in sezione dell'impeller e dell'albero della turbina. Fonte [4]
Descriviamo ora il funzionamento del regolatore di velocità: partendo da una situazione di regime
permanente (velocità costante), supponiamo che a seguito di un disturbo in rete, la velocità di
rotazione aumenti. A seguito dell’aumento di velocità, e in base al funzionamento dell’impeller, la
pressione in uscita da quest’ultimo aumenta. l’aumento di questa pressione provoca la compressione
del soffietto indicato con la lettera S in figura 7.4. Di conseguenza l’asta indicata con la lettera L si
alza e allontana la coppetta sottostante dalla sua sede. Si ha un aumento della sfuggita di olio dalla
camera e la conseguente diminuzione di pressione del collettore dell’olio di regolazione. Olio di
regolazione comanda direttamente il servomotore di comando delle valvole e una diminuzione della
sua pressione provoca chiusura delle valvole. Ne deriva quindi una diminuzione della potenza in
turbina in modo da contrastare l’aumento della velocità.
73
Nella figura 7.6 viene mostrato il servomotore di comando valvole che riceve l’olio modulato dal
regolatore di velocità. Le sue parti principali sono: un cassetto ricevitore di pressione (riceve l’olio
modulato dal regolatore), un cassetto distributore dotato di camicia mobile e una camma di
retroazione di posizione (sposta la camicia del cassetto distributore in funzione dello spostamento del
pistone motore).
Supponiamo di essere nel caso precedente. A seguito del disturbo si ha quindi la diminuzione
della pressione dell’olio modulato; questo provoca lo spostamento verso l’alto del pistoncino pilota.
Contemporaneamente al pistoncino pilota, viene spinto verso l’alto anche il cassetto distributore che
apre le luci di reflusso dell’olio di potenza consentendo al pistone motore di abbassarsi e chiudere
tramite il sistema di leve le valvole che regolano il flusso di vapore alla turbina. Si ha quindi la
diminuzione di potenza in turbina.
Come già detto, i regolatori di velocità sono composti da un sensore di velocità e da un sistema
che al variare della velocità va ad agire sul distributore della turbina.
Nei regolatori moderni, la rilevazione della velocità viene effettuata generalmente con la ruota
fonica. La ruota fonica è composta da una ruota dentata e da uno o più sensori induttivi o ottici.
74
I sensori induttivi sono costituiti da una bobina avvolta attorno a un magnete permanente. Grazie
alla conformazione della ruota, il sensore rileva la variazione di flusso in corrispondenza ad ogni
dente e invia un segnale alternato nel quale ogni picco corrisponde a un dente della ruota fonica. Si
può quindi elaborare il segnale del sensore in modo da conoscere la velocità angolare della ruota
fonica. Questo sistema di misura non presenta parti delicate in movimento ed è quindi molto più
affidabile e richiede minor manutenzione rispetto ad esempio al pendolo di Watt.
Il segnale dei sensori induttivi va quindi inviato a un sistema in grado di elaborarlo e ricavarne il
valore di velocità.
Esistono poi ruote foniche con sensori ottici (encoder); in questo caso, la ruota è costituita da un
laser e un sensore che lo rileva quando tra i due non è presente un ostacolo. Questi sensori vengono
posizionati in modo da rilevare la presenza dei denti della corona. Alla rotazione della ruota si avrà in
uscita del sensore una sequenza di impulsi la cui frequenza sarà proporzionale alla velocità di
rotazione.
75
seconda della sua posizione si hanno nei due secondari determinati valori della tensione di uscita in
quanto è diversa la forza elettromotrice indotta. Conoscendo la relazione tra tensione di uscita e
posizione del nucleo mobile, si può determinare la variazione di posizione dell’oggetto solidale a
quest’ultimo.
Il sistema oleodinamico viene modernizzato, oltre che con l’introduzione degli LVDT, anche con
l’utilizzo di attuatori lineari per la movimentazione dei cilindri comandati da un controllore a
microprocessori. Schematicamente il sistema oleodinamico si può così rappresentare:
76
Figura 7. 14 Schema regolatore di velocità con attuatore
Tutto il sistema come già detto viene controllato da un controllore a microprocessore. A questo
controllore giungono tutti i segnali necessari per la regolazione dell’impianto, a partire dalla velocità
di rotazione del gruppo di generazione. Nel controllore è possibile inserire tutti i modelli matematici
dei componenti dell’impianto che entrano in gioco nella regolazione della velocità come il modello
del distributore e della condotta forzata. In base ai segnali di ingresso e ai modelli matematici inseriti,
il controllore comanda in uscita gli attuatori e riceve i feedback di posizione dagli LVDT. Ovviamente
i parametri del regolatore così ottenuto devono soddisfare quelli richiesti dall’ente distributore, per
esempio lo statismo.
Con questi interventi di modernizzazione degli impianti si riesce a ottener un ciclo di regolazione
con un margine di fase maggiore rispetto al regolatore oleodinamico precedente.
È possibile inoltre monitorare da remoto in tempo reale le condizioni dell’impianto, sempre
tramite il controllore.
Il sistema deve presentare una grande affidabilità e continuità di esercizio, per questo si installano
due controllori di cui uno è il master e l’altro è lo slave. In questo modo è assicurata la continuità di
esercizio dell’impianto anche in caso di guasto a uno dei controllori. Stessa regola vale per i sensori di
misurazione della velocità: i sensori utilizzati sono due o più di due.
77
Figura 7. 15 Controllore a microprocessori Emerson Ovation per il controllo della regolazione di
velocità con moduli ingresso/uscita. Fonte [4]
78
Figura 7. 16 Schema regolatore di velocità SIMO. Fonte [9]
79
Figura 7. 17 Risposta in frequenza del sistema idraulico adottato per la simulazione. Fonte [9]
80
funzionare anche al variare, entro certi limiti, dei parametri principali dell’impianto. Il sistema così
ottenuto viene quindi confrontato con i classici sistemi di controllo che prevedono regolatori PI e PID.
La risposta dei tre diversi sistemi di controllo viene confrontata con un segnale di ingresso a
gradino. Per gli impianti di produzione è particolarmente importante che l’overshoot sia il più limitato
possibile (per evitare problemi alla turbina).
La risposta al gradino dei tre diversi sistemi presenta le seguenti caratteristiche:
Vediamo dalla figura che il sistema migliore risulta quello dimensionato con l’ottimizzazione H∞.
La risposta infatti presenta un minor tempo per arrivare a regime e un overshoot molto limitato
rispetto ai regolatori PI e PID.
Nella prossima figura possiamo vedere invece come risponde il sistema alla variazione del
carico elettrico:
81
Vediamo come nel caso studiato l’ottimizzazione H∞ sia decisamente la soluzione migliore per
l’impianto. La sua risposta non presenta oscillazioni e va a regime in un tempo minore.
L’impianto di produzione risulta quindi meno sollecitato e risponde prontamente alle variazioni
del carico.
Dalla simulazione si vede come le prestazioni siano migliorate dopo il revamping del regolatore. la
risposta al gradino del regolatore dopo il revamping presenta una sovra elongazione ridotta (minore
stress all’impianto di produzione) e una risposta più veloce alla variazione di frequenza.
82
La migliore risposta al gradino è una delle migliorie che si va ad aggiungere agli altri vantaggi
descritti in questo capitolo.
83
Capitolo 8
La regolazione primaria inizia negli istanti successivi al cambiamento della frequenza della rete.
Viene eseguita in maniera automatica dai regolatori di velocità dei singoli gruppi di produzione.
L’azione del regolatore di velocità, che rileva la variazione di frequenza, è quella di aumentare, a
fronte di una diminuzione della frequenza (supponendo una perdita di potenza di generazione), la
potenza erogata dal proprio motore primo.
85
Considerando la totalità dei generatori che effettuano la regolazione primaria, la potenza immessa
nel sistema viene progressivamente aumentata fino a compensare quella perduta. I gruppi citati
dovranno quindi avere una riserva di potenza da poter sfruttare; questa potenza viene definita riserva
di regolazione primaria.
L’azione autonoma dei regolatori cessa nel momento in cui l’equilibri della potenza in rete viene
ristabilito e quindi la frequenza assume un valore costante. Al termine della regolazione primaria, la
rete avrà una frequenza inferiore a quella di programma e la riserva complessiva di regolazione sarà
anch’essa inferiore a quella precedente il disturbo. Cambia in generale anche il valore della potenza
importata da altre reti, poiché anch’esse hanno contribuito a sopperire alla perdita di generazione
ipotizzata.
L’ente che definisce le unità di produzione che partecipano alla regolazione primaria è Terna SpA,
che gestisce la Rete di Trasmissione Nazionale. Terna prescrive che il servizio di regolazione primaria
è obbligatorio per tutte le unità di produzione con potenza efficiente (potenza massima che l’unità di
produzione può produrre con continuità, nel caso di impianti termoelettrici, o per un determinato
numero di ore per gli impianti idroelettrici) non inferiore a 10 MW, ad eccezione di impianti di fonti
rinnovabili non programmabili.
La regolazione del regolatore deve essere sempre automatica, e in particolare, per le unità
termoelettriche, la regolazione deve essere coordinata con il frequency-bias, mentre per gli impianti
idroelettrici non deve essere limitata né in gradiente né in banda, ma solo vincolata ai parametri
idraulici degli organi di regolazione e dalla disponibilità energetica. In ogni caso il gradiente di
regolazione di velocità deve essere superiore a quello della regolazione secondaria.
Terna prescrive anche quelle che devono essere le prestazioni funzionali minime dei regolatori di
velocità:
- precisione della misura di velocità migliore dello 0,02% in qualsiasi condizione di
funzionamento,
- insensibilità del regolatore di velocità, esclusa la parte di misura, non superiore a ±10
mHz,
- capacità di far funzionare l’unità di produzione in modo stabile a tempo indefinito, con
qualunque grado di statismo impostabile tra 2 e 8%, per qualunque punto di lavoro
corrispondente alle frequenze comprese fra 47,5 Hz e 51,5 Hz, e con qualunque carico
compreso fra il minimo tecnico e la potenza massima generabile dal gruppo
Sapendo che all’occorrenza le unità di produzione devono poter erogare una potenza
supplementare per riportare l’equilibrio in rete, queste non possono lavorare costantemente alla loro
potenza massima, ma devono avere un certo margine disponibile. Il margine che Terna chiede ai
distributori non deve essere inferiore a 1,5% della potenza efficiente dichiarata nel RUP (registro delle
unità di produzione), quando la potenza erogata è pari alla potenza massima erogabile o è pari alla
potenza minima erogabile. In base a questo, l’unità di produzione può essere esercita nel campo di
funzionamento ammissibile seguente:
- Pmin=PMT+1,5% Peff
86
- PMAX=Pmax erogabile-1,5% Peff.
con PMT potenza di minimo tecnico
Diversi requisiti hanno gli impianti situati in Sardegna e Sicilia: queste unità di produzione devono
avere a disposizione una riserva primaria non inferiore al 10% della propria potenza efficiente; da
questo valore derivano conseguentemente le potenze a cui questi impianti possono lavorare.
Nel caso di regolazione primaria della frequenza, si possono distinguere due casi:
- fornitura del servizio in condizioni normali di esercizio,
- fornitura del servizio in condizioni di emergenza.
In condizioni normali di esercizio, ogni unità di produzione che deve erogare una quota ∆Pe della
banda di riserva primaria disponibile tenendo conto della variazione di frequenza ∆f e del grado di
statismo permanente σp impostato nel regolatore secondo la relazione:
∆1 5jGG
∆5j = − 100
50 Æ
(7.1)
Di questa quota non è consentita nessuna limitazione in gradiente e ampiezza. Si devono in oltre
rispettare i seguenti requisiti:
1) entro 15 secondi dall’inizio della variazione di frequenza deve essere erogata almeno
metà della potenza espressa dalla 7.1.
2) entro 30 secondi dall’inizio della variazione di frequenza deve essere erogata tutta la
potenza richiesta dall 7.1.
Se il transitorio di frequenza si prolunga oltre i 30 secondi, le unità di produzione devono
proseguire con l’aumento o diminuzione di potenza erogata in funzione della variazione di frequenza.
Quando le unità di produzione hanno messo in gioco la totale ΔPe devono essere in gradi mantenere la
nuova potenza erogata per almeno 15 minuti consecutivi.
Per le condizioni di emergenza, la potenza da fornire è sempre quella espressa dalla 7.1. Le unità
di produzione devono fornire tutta la potenza disponibile fino al raggiungimento di uno dei limiti
tecnici (massimo e minimo), dichiarati nel RUP. La potenza deve essere fornita con il massimo
gradiente tollerabile dall’impianto, certificato con prove di tipo: non si deve arrivare a compromettere
il funzionamento continuativo dell’impianto sia per motivi elettrici, termici o meccanici. Tale
gradiente deve essere superiore al 3% al minuto della Peff. Anche in questo caso, la nuova potenza
erogata dall’unità di produzione deve essere mantenuta per almeno 15 minuti consecutivi.
A seconda del tipo di impianto, Terna prescrive che i regolatori di velocità siano tarati come
segue:
- impianti idroelettrici:
o grado di statismo 4%;
o banda morta intenzionale non superiore a ±10mHz
- impianti termoelettrici:
o grado di statismo 5%;
87
o banda morta intenzionale non superiore a ±10 mHz per le unità a vapore a ciclo
semplice, ±20 mHz per le unità turbogas o a ciclo combinato.
Per le unità di non recente costruzione sono ammesse delle deroghe: non si deve superare il valore
di ±30mHz. Per quanto riguarda le turbine a gas, che hanno una potenza prodotta che varia
sensibilmente con le condizioni ambientali, Terna impone che il profilo di potenza offerto sia
preventivamente calcolato in funzione della temperatura prevista, utilizzando i diagrammi di
prestazione fornite dal costruttore (ad ogni valore di temperatura corrisponde un certo valore massimo
di potenza erogabile). Questi impianti devono comunque essere in grado di fornire la banda di riserva
anche a temperature diverse da quelle previste.
88
8.1.3 Regolazione terziaria
Una volta terminata la fase di regolazione secondaria, supponendo a seguito di una diminuzione di
frequenza, sono diminuite le riserve di regolazione secondaria, che devono quindi essere ricostituite.
Nella regolazione terziaria il gestore della rete di trasmissione nazionale impartisce ai singoli
produttori delle disposizioni di esercizio, per esempio l’entrata in servizio di gruppi di generazione di
riserva o la variazione di potenza prodotta da quelli in servizio. Questa regolazione quindi non è
automatica ma viene eseguita manualmente su richiesta
89
Conclusioni
Nell’elaborato sono stati descritti come erano e come vengono ora realizzati i regolatori di velocità
per le centrali termoelettriche ed idroelettriche.
A conclusione di questo è interessante sottolineare quali siano i vantaggi di adeguare i regolatori di
velocità con l’inserimento di nuove tecnologie in diverse parti dell’anello di regolazione.
Con la sostituzione dei sensori di velocità da meccanici a induttivi è possibile avere più misure di
velocità da inviare al regolatore a microprocessore. Inoltre, la minor presenza di parti in movimento
consente di migliorare ulteriormente l’affidabilità di questa parte importante del regolatore.
L’innovazione più importante nei nuovi regolatori è sicuramente l’introduzione del regolatore a
microprocessore. Questo dispositivo, attraverso schede di acquisizione, può ricevere segnali di input
come la velocità di rotazione del gruppo o la posizione (e quindi l’apertura) della valvola di
immissione in turbina. Nel regolatore si possono implementare le funzioni di trasferimento
caratteristiche dell’impianto regolato in modo da ottenere come risultato la replica virtuale del sistema
da controllare. In base a questo è possibile raggiungere le migliori prestazioni di regolazione. Le
schede di output del regolatore vanno a comandare l’apertura o chiusura delle valvole in base
all’evoluzione del segnale di ingresso. I dispositivi che agiscono sulle valvole possono essere cilindri
idraulici o attuatori lineari controllati da circuiti idraulici o elettricamente. La loro posizione viene
monitorata da uno o più sensori di posizione lineare in modo da avere un segnale di feedback (non
svolge funzioni di regolazione).
91
Tutti questi provvedimenti portano in estrema sintesi i seguenti importanti vantaggi:
• minor manutenzione sull’impianto
• prestazioni (margine di fase e frequenza di taglio aumentati) e miglior sfruttamento
dell’impianto
• gestione remota per il controllo e il comando dell’impianto
• affidabilità del regolatore (ridondanza di sensori e controllori a microprocessore)
• adeguamento dei regolatori ai parametri stabiliti dal gestore della rete
Viene poi presentato un sistema di controllo in cui vengono considerati svariati segnali per il
controllo dell’impianto, come per esempio la posizione delle valvole di immissione, la potenza
erogata e l’angolo di carico. questi segnali si aggiungono a quello di velocità. Questa tipologia di
controllo viene chiamata SIMO (single input multi output) e per ottenere le massime prestazioni non è
sufficiente progettare il regolatore classico PI o PID. In questo caso la progettazione è molto più
complessa e richiede l’ottimizzazione H∞ che richiede l’uso di matrici e di una funzione peso (vedi
bibliografia). I risultati ottenuti con quest’ultima teoria di progettazione mostrano delle prestazioni
nettamente migliori sotto tutti gli aspetti rispetto alla teoria classica dei regolatori PI e PID.
92
Appendice 1
Introduzione a Simulink
A1.1 Matlab
Matlab, abbreviazione di Matrix Laboratory, è un linguaggio di programmazione per applicazioni
scientifiche e numeriche che offre un vasto set di funzioni predefinite, con la possibilità di scriverne di
nuove. L’elemento di base di Matlab è un array (vettore), questo consente un facile utilizzo ed
elaborazione delle matrici attraverso algoritmi più semplici di quelli che sarebbero necessari con l’uso
di un linguaggio diverso tipo Fortran o C. È usato in ambiente universitario come attrezzo didattico,
ma anche utile nell’ingegneria, nella matematica e nella scienza.
Nel programma sono presenti delle toolboxes che offrono delle soluzioni specifiche per risolvere
determinate categorie di problemi. Alcune di queste toolboxes sono:
• elaborazione dei segnali,
• sistemi di controllo;
• simulazione e molti altri.
Quello che viene utilizzato nelle nostre analisi è la simulazione, tramite l’uso di Simulink.
Simulink è un ambiente per la modellizzazione, l’analisi e la simulazione di sistemi matematici e
fisici. Offre un’interfaccia grafica per la costruzione di modelli tramite schemi a blocchi, presenti in
una libreria di blocchi predefiniti. La modellizzazione di un sistema risulta quindi essere veloce e
relativamente semplice e consente di visualizzare i risultati in tempo reale durante la simulazione.
93
Per avviare Simulink è sufficiente digitare ‘simulink’ nel command window di Matlab oppure fare
click sull’apposito pulsante in alto a sinistra.
A1.2 Simulink
Una volta aperto Simulink si presenta all’utente una finestra:
A sinistra sono presenti tutte le categorie che contengono blocchi specifici per determinate
applicazioni, a destra è possibile navigare nelle sottocategorie di ogni sezione. Da questa finestra è
possibile caricare oppure creare un nuovo file.
Cliccando sull’icona con il foglio bianco si inizia a operare su un nuovo file, dove sarà necessario
inserire vari blocchi e collegarli tra loro per costruire il nostro sistema da analizzare.
94
A1.3 Costruzione schema a blocchi
Per costruire il sistema sono stati inseriti i blocchi necessari presi dalle librerie presenti nel
programma. Si possono notare i blocchi con le funzioni di trasferimento, la sorgente di segnale a
gradino, i blocchi sommatori, e i blocchi che consentono la visualizzazione delle grandezze scelte nel
tempo. Per ogni blocco è possibile definire delle proprietà, per esempio i valori dei coefficienti delle
funzioni di trasferimento,
95
Figura A1. 3 Finestra impostazione parametri blocco G1
96
In funzione di analisi successive alle simulazioni, torna utile una proprietà dei blocchi che graficano la
soluzione del sistema: la possibilità di salvare su un vettore tutti i dati del grafico.
Una volta impostate le proprietà dei blocchi, è necessario impostare i parametri di simulazione:
97
A1.4 Simulazione
Una volta impostato i vari parametri di interesse, è possibile avviare la simulazione tramite il pulsante
play:
Se tutti i parametri sono stati impostati correttamente, il programma non presenta messaggi di errore
ed è possibile procedere alla visualizzazione dei risultati.
Se ci interessa visualizzare l’andamento di deltaF nel tempo sarà suffuciente fare doppio click
sull’icona
A seguito del doppio click, ci verrà mostrata una finestra con il grafico in funzione del tempo
98
Gli script non sono altro che m-file (file con estensione .m compatibili con Matlab) che
contengono all’interno diverse istruzioni che Matlab esegue nel momento in cui lo script viene
avviato. Tramite l’uso di questi file è possibile definire variabili e richiamare funzioni per effettuare le
operazioni necessarie. Tutte le variabili usate negli script vengono automaticamente salvate nel
workspace e quindi utilizzabili anche al di fuori dello script stesso.
In questo caso gli script possono essere usati per fare simulazioni ripetute al variare di una o più
variabili del sistema, per poi confrontare l’andamento delle grandezze d’uscita in unico grafico, come
nell’esempio seguente:
Figura A1. 9 Grafico risultante con comparazione di diverse simulazioni al variare di un parametro
Nella figura A1.9 è mostrato il grafico che confronta l’andamento di deltaF al variare dell’energia
regolante: è stato ottenuto tramite il seguente script:
for Kr=25000000:2000000:35000000
sim('untitled')
%delta(:,a)=deltaF(:,2);
figure(2)
hold on
99
plot(deltaF(:,1),deltaF(:,2))
title('Andamento deltaF in funzione di Kr (energia regolante)
[MW/Hz]')
xlabel('Tempo [s]')
ylabel('deltaF [Hz]')
end
100
Bibliografia
[1] Fellin, L., & Benato, R. (2011). Impianti Elettrici. Wolters Kluwer Italia.
[2] Marconato, R. (1985). Sistemi elettrici di potenza (Vol. 2). Milano: Cooperativa libraria
universitaria del politecnico.
[3] Paolucci, A. (1998). Lezioni di trasmissione dell'energia elettrica. Padova: CLEUP.
[4] Da Lio F. (2010). Analisi del revamping di un regolatore turbina di un gruppo di
produzione da 320 MW
[5] Terna Spa. Prescrizione tecnica - Partecipazione alla regolazione di frequenza e
frequenza-potenza
[6] Terna Spa. Specifiche tecniche per la verifica e valorizzazione del servizio di regolazione
primaria di frequenza
[7] Viani D. Revamping del sistema di controllo frequenza-potenza di una filiera di impianti
idroelettrici di taglia media
[8] Manunza A., Maccioni A. Analisi del comportamento dinamico delle centrali idroelettriche
[9] Eker I. Governors for hydro-turbine speed control in power generation: a SIMO robust
design approach, Energy conversion an magaement 45 (2004) 2207-2221
[10] Dulau M, Bica D., Simulation of speed steam turbine control, Procedia Technology 12
(2014) 716-722
[11] P. Nelson, S. Muthukumar, Frequency regulation by free governor mode of operation in
power station
[12] E. Malusardi, G. Tagliabue, La risposta delle centrali alla regolazione di frequenza:
impatti tecnici/economici e vincoli sulle prestazioni per i vari tipi di centrale di nuova
generazione e convenzionali
[13] ENTSOE, Documentation on controller tests in test grid configurations
[14] Paradiso B. Impianti idroelettrici e turbine idrauliche, Politecnico di Milano
[15] H. Kwakernaak Robust control and H∞ optimisation - Tutorial paper, Automatica, Vol. 9,
No. 2, pp. 255-273, 1993
.
101