Programmazione idraulico di breve termine –Obiettivo della programmazione di breve termine...

Programmazione idraulico Programmazione idraulico di breve terminedi breve termine

– Obiettivo della programmazione di breve termine dell’idraulico (programmazione giornaliera) è quella di utilizzare al meglio la risorsa idraulica disponibile minimizzando la produzione termoelettrica.

– I dati del problema, relativi a ciascun impianto idroelettrico, riguardano:

– la conoscenza dei diagrammi di carico del giorno successivo (deterministici)

– la conoscenza degli apporti idrici del giorno successivo (deterministici)

– la conoscenza (dalla programmazione di medio termine) del livello iniziale del bacino (ri,0) e dei volumi di acqua da utilizzare nel giorno successivo (V i)

Programmazione idraulicoProgrammazione idraulico

– La principale ipotesi relative allo sviluppo del problema sono:

– il rendimento idraulico degli impianti viene assunto costante; in effetti esso dipende dalla potenza erogata dalla singola turbina nonché dalla portata complessiva di ciascuna condotta che alimenta generalmente più turbine. Tuttavia è possibile verificare che sopra il minimo tecnico e fino alla potenza massima della centrale il rendimento complessivo si mantiene abbastanza costante;

– Il carico che deve essere coperto è quello totale idro-termoelettrico diminuito della quota parte relativa ad impianti non modulati (acqua fluente, impianti termoelettrici di base, impianti contrattualizzati, ecc.); cioè si considera solo la quota parte di carico che può essere coperta con centrali di modulazione e di regolazione.


– Occorre innanzitutto osservare che:

l

k

l

k

n

jjkk

n

jjkk Pfqe

1 1 1,

1, )(

– dove:ek = costo specifico del combustibile del gruppo k-esimoqk,j= portata di combustibile al gruppo k-esimo nel j-esimo intervallofk = funzione di costo del gruppo k-esimo (termoelettrico)Pk,j = potenza erogata dal gruppo k-esimo nel j-esimo intervallo

– D’altra parte poiché le funzioni di costo dipendono dal quadrato della potenza si può scrivere:

n

j

l

kjk

l

k

l

k

n

jjkk

n

jjkk PPfqe

1

2

1,

1 1 1,

1, )(


– Di conseguenza trovare il minimo della funzione di costo

l

k

n

jjkk qe

1 1,

– equivale a trovare il minimo della seguente funzione

n

j

l

kjkP

1

2

1,

Programmazione idraulicoProgrammazione idraulico– Se si indica con

sj = il carico idro-termoelettrico da coprire nell’intervallo j-esimoxi,j = la potenza erogata dalla centrale i-esima nel periodo j-esimo

– Si ottiene

m

i

l

kjkjij Pxs

1 1,, da cui

m

ijij

l

kjk xsP

1,

1,

E ancora

2

1 1,

1

2

1,

n

j

m

ijij

n

j

l

kjk xsP

Quindi il minimo della funzione costo di produzione coincide con quello della funzione

min2

1 1,

n

j

m

ijij xs

Programmazione idraulicoProgrammazione idraulico– I vincoli a cui deve sottostare la funzione obiettivo prima definita

sono i seguenti:

i

j

h i

hihiii

ijii

i

n

j i

ji

Rc

xarr

xxx

Vc

x

1

,,0,

,

1

,

– doveci = coefficiente energetico della centrale i-esimaai,j = apporti al bacino della centrale i-esima nel periodo j-esimori = limite di svuotamento del bacino i-esimoRi = limite di sfioro del bacino i-esimo

Ottimazione economicaOttimazione economica– Consideriamo la seguente funzione di costo da minimizzare sotto

opportuni vincoli:

t

ng

kk

ng

kkk

CP

PC

min)(1

– doveCk = costo orario di produzione della centrale k-esimaPk = potenza oraria media erogata della centrale k-esimaCt = carico da coprire con la produzione termicang = numero dei generatori “in giri”

Ottimazione economicaOttimazione economica– Applicando il metodo dei moltiplicatori di Lagrange () si ha:

)()(1

ng

kkt

ng

kkk PCPC

– Il cui minimo, che coincide con quello della funzione di costo prima definita, può essere calcolato annullando le derivate parziali prime rispetto alle variabili Pk:

..........

0

.........

01

1

1

j

j

j P

C

P

P

C

P

Ottimazione economicaOttimazione economica

– In definitiva la condizione di minimo della funzione costo si ottiene quando tutti i generatori funzionano in corrispondenza dello stesso valore di , che rappresenta il costo incrementale (detto anche costo marginale).

– Il criterio prima esposto va sotto il nome di criterio degli uguali costi incrementali (uguali costi marginali).

– L’unica condizione da imporre perché il sistema sia convergente è che la funzione di costo sia monotona non decrescente.


Criterio degli uguali costi incrementali

Ottimazione economicaOttimazione economica– Se introduciamo un vincolo che tiene conto della potenza

massima e del minimo tecnico di ciascun gruppo il problema del dispacciamento a minimo costo assume la seguente forma:

kkk

t

ng

kk

ng

kkk

PPP

CP

PC

min)(1

– e la funzione lagrangiana diventa (condizioni di Kuhn-Tucker)

ng

k

kkk

ng

kkkk

ng

kkt

ng

kkk PPPPPCPC

1111

)()(


– Dato un insieme di generatori è possibile definire una curva monotona non decrescente che definisce un legame univoco (Pk; ); tale curva identifica il costo incrementale del sistema al variare della potenza complessivamente prodotta.


Costo incrementale di sistema

Dispacciamento generatori Dispacciamento generatori termicitermici– La funzione di costo di produzione che deve essere minimizzata

in funzione delle potenze erogate dai diversi generatori ed i relativi vincoli è la seguente:

IuxI

uuu

xxx

uxG

uxF

uxE

),(

0),(

0),(

min),(

– dove F e G rappresentano le equazioni di load-flow attivo e reattivo della rete.

Dispacciamento generatori Dispacciamento generatori termicitermici– I costi che si considerano sono tutti quelli variabili legati al

funzionamento di una centrale, e cioè sostanzialmente i costi di combustibile (gli altri costi possono essere considerati fissi e quindi indipendenti dalla potenza erogata dalla centrale).

– Il problema così come è stato posto in maniera generale è estremamente complesso; infatti si deve tenere conto di

– costi di funzionamento a regime

– costi di avviamento dei gruppi

– dispacciamento della potenza attiva e di quella reattiva

– Di conseguenza per trovare una soluzione è necessario effettuare alcune semplificazioni.

Dispacciamento generatori Dispacciamento generatori termicitermici– Poiché la soluzione del problema è funzione di un elevato numero

di variabili e ciascuna di esse ha un peso paragonabile alle altre, la funzione obiettivo presenta un minimo molto piatto; in tali condizioni si può dimostrare che la soluzione trovata separando il problema in più sottoproblemi e determinando il minimo di ciascuno di questi praticamente coincide con la soluzione globale del problema.

– Di conseguenza per trovare la soluzione del problema generale verranno determinate le soluzioni di alcuni sottoproblemi nei quali è scindibile il problema originario:

– determinazione dello U.C.;

– separazione del dispacciamento attivo e reattivo.

Unit CommitmentUnit Commitment– La prima semplificazione che viene fatta riguarda la

determinazione dei gruppi termici che devono rimanere in giri durante il periodo considerato; tale operazione prende il nome di Unit Commitment.

– Nello U.C. il periodo temporale che tipicamente viene preso in considerazione è la settimana.

– L’obiettivo che ci si pone è quello di determinare, per ciascun intervallo (generalmente l’ora) del periodo considerato, quali sono i gruppi che devono rimanere in giri garantendo il minimo costo complessivo di esercizio.

– La funzione obiettivo è costituita dalla somma dei costi di funzionamento a regime e di quelli di avviamento di tutti i gruppi del parco.

Unit Commitment:Unit Commitment:programmazione delle fermateprogrammazione delle fermate

– Questa tecnica può essere utilizzata quando il parco macchine sia ridotto (es.: i gruppi di una centrale); si stabilisce in anticipo quale gruppo si vuole fermare.

– Per la soluzione del problema occorre conoscere i costi di funzionamento e di avviamento di ciascuno dei gruppi oltre che, ovviamente, il diagramma di carico da coprire.

)1()( , eCPPC kakkk

– Se la condizione è verificata viene arrestata la macchina k-esima.

Unit Commitment:Unit Commitment:programmazione completaprogrammazione completa

– In un parco composto da m macchine il problema di minimo vincolato posto è il seguente:

k

n

hhk

hkhkhkhkhk

h

l

khk

l

k

n

h

l

kkkahkhk

PPP

LP

CPC

1,

,,,,,

1,

1 1 1,,,

1

min)(

Separazione Separazione attivo/reattivoattivo/reattivo

– Dalle equazioni di load-flow è nota la seguente relazione (con ovvia notazione di simboli):

VJJ

JJ

Q

P

VQQ

VPP

,,

,,

– Per le reti di trasmissione con buona approssimazione si può ipotizzare che:

VJ

J

Q

P

VQ

P

,

,

0

0

Separazione Separazione attivo/reattivoattivo/reattivo– Quindi le variazioni di transito della potenza attiva sono in prima

approssimazione indipendenti dalla tensione e quelle di potenza reattiva dagli angoli; in altre parole la potenza attiva e quella reattiva sono debolmente legate.

– Di conseguenza il problema può essere risolto in maniera iterativa determinando prima i transiti di potenza attiva a tensioni fissate e successivamente quelli di potenza reattiva a fissati; in questo modo il problema generale viene declassato alla risoluzione di due sottoproblemi più semplici e più agevoli da risolvere:

– Dispacciamento della potenza attiva

– Dispacciamento della potenza reattiva.

Considerazioni sulle Considerazioni sulle equazioni di l.f.equazioni di l.f.

– L’espressione in forma polare delle equazioni di load-flow consente di esplicitarle applicando il teorema del Dini.

– Tale teorema afferma che data una funzione implicita di variabili di stato e di controllo, sotto opportune ipotesi riguardanti la continuità della funzione e la sua derivabilità rispetto alle suddette variabili, è possibile definire due intorni, uno nello spazio delle variabili di stato ed uno in quello delle variabili di controllo, tali per cui se le variabili appartengono ai rispettivi intorni è possibile esplicitare le variabili di stato in funzione di quelle di controllo.

Considerazioni sulle Considerazioni sulle equazioni di l.f.equazioni di l.f.– L’espressione in forma polare delle equazioni di load-flow è la

seguente:

0)(

0)cos(

1,,

1,,

n

kkikikikii

n

kkikikikii

senYVVQ

YVVP

– Ossia, in forma implicita (indicando con x le variabili di stato e con u quelle di controllo):

0),(

0),(

uxG

uxF F = equazioni di l.f. della potenza attiva

G = equazioni di l.f. della potenza reattiva

Considerazioni sulle Considerazioni sulle equazioni di l.f.equazioni di l.f.– Dalla applicazione del teorema del Dini, limitandoci alle

equazioni della potenza attiva, possiamo scrivere:

uJx F

– Il cui Jacobiano JF vale:

u

F

x

FJF

1

Considerazioni sulle Considerazioni sulle equazioni di l.f.equazioni di l.f.– Verifica: poniamo

nn

nn

PPu

Px

,,...,

,,...,

11

11

– L’insieme delle F(x,u)=0 risulta

0),...,(

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

0),...,(

0),...,(

11

1122

1111

nnn

n

n

fP

fP

fP


1..

0......

0......

0..

11

1

1

1

1

n

nn

n

ff

ff

x

F

1

01

11

JfJ

x

Fn

Ix

F

x

F

1


00

0I

u

F

0

0

00

0

1

0

1

1

1

11

JfJ

I

JfJ

u

F

x

FJ

n

nF

– In base alla definizione:

Considerazioni sulle Considerazioni sulle equazioni di l.f.equazioni di l.f.– Particolare interesse riveste l’ultima riga che scritta in forma

esplicita diventa:

11,22,11, ....... nnnnnn PAPAPAP

– Dove le costanti A rappresentano gli elementi del vettore:

1

Jfn


– Tali costanti costituiscono però anche le e cioè le variazioni della potenza richiesta nel nodo di saldo in funzione delle iniezioni negli altri nodi della rete.

– Da un punto di vista fisico i coefficienti A rappresentano quindi i contributi alle perdite di rete per effetto delle iniezioni di potenza negli altri nodi.

i

n

P

P

Dispacciamento potenza Dispacciamento potenza attivaattiva– Il problema posto è quello di definire il dispacciamento della

potenza attiva dei gruppi termoelettrici che sono “in giri” (a valle dello U.C.) avendo definito come funzione obiettivo la minimizzazione dei costi variabili conoscendo il carico da coprire con il termico e tenendo conto delle curve di capability dei generatori.

– Il problema posto può essere risolto con diverse approssimazioni:

– I approssimazione: non si tiene conto delle perdite di rete (modello sbarra)

– II approssimazione: si tiene conto delle perdite di rete adottando differenti modelli

Dispacciamento potenza Dispacciamento potenza attivaattiva– I approssimazione: modello sbarra

G1 GngGjG2 Ct

– Viene applicata la tecnica degli uguali costi incrementali già dimostrata in precedenza

Dispacciamento potenza Dispacciamento potenza attivaattiva– II approssimazione: modello sbarra ad uguali costi incrementali

corretti

G1 GngGjG2 Ct pr

– Anche in questo caso viene applicata la tecnica degli uguali costi incrementali introducendo però un ulteriore condizione.

Dispacciamento potenza Dispacciamento potenza attivaattiva– Il problema di minimo costo è così definito

– dove p rappresenta le perdite nella rete di trasmissione; dalla applicazione delle condizioni di Kuhn-Tucker si ottiene

kkk

t

ng

kk

ng

kkk

PPP

pCP

PC

min)(1

ng

k

kkk

ng

kkkk

ng

kkt

ng

kkk PPPPPpCPC

1111

)()(

Dispacciamento potenza Dispacciamento potenza attivaattiva– Annullando le derivate prime rispetto alle potenze iniettate si ha

– In assenza dei vincoli la soluzione porta alla seguente espressione

0)1(

kkkk

k

P

p

P

C

)1(k

k

k

Pp

PC

– Dove la correzione tiene conto della “distanza” elettrica delle iniezioni.

Dispacciamento potenza Dispacciamento potenza attivaattiva– Per determinare la funzione delle perdite p si può anche in

questo caso ricorrere a modelli di differente approssimazione.

– Modello di I approssimazione: le perdite vengono definite sulla base di una funzione semi-empirica del tipo

1

1

1

1,,2

1

1,10 2

1 n

i

n

jjiji

n

iii PPBPBBp

– Da tale espressione si ottiene:

1

1,,2,1

n

jjjii

i

PBBP

p

Dispacciamento potenza Dispacciamento potenza attivaattiva– Modello di II approssimazione: le perdite vengono definite sulla

base alla seguente formulazione (n = nodo di saldo)

pCPPP t

n

ini

n

ii

1

11

– Da cui derivando si ottiene:

ii

n

P

p

P

P

1

– Dalla esplicitazione delle variabili di stato in funzione di quelle di controllo si ottiene:

1

Jf

P

P n

T

i

n

Procedura di dispacciamento Procedura di dispacciamento potenza attivapotenza attiva

1. Determinazione delle potenze erogate da ciascun generatore con il criterio degli uguali costi incrementali;

2. Calcolo di load-flow e stima delle perdite di rete (nodo di saldo);

3. Calcolo dei con la voluta approssimazione;

4. Determinazione delle potenze erogate da ciascun generatore con il criterio degli uguali costi incrementali corretti;

5. Verifica di convergenza ed eventuale nuova iterazione dal punto 2.

iP

p

Dispacciamento potenza Dispacciamento potenza reattivareattiva– Il problema posto è quello di definire il dispacciamento della

potenza reattiva dei gruppi termoelettrici che sono “in giri” (a valle dello U.C.) avendo definito come funzione obiettivo:

– Un dispacciamento a minime perdite di potenza (attiva);

– Un dispacciamento a minima potenza reattiva prodotta;

– Un dispacciamento a massima tensione prodotta.

– La funzione obiettivo è naturalmente una sola tra quelle definite; quella che viene tipicamente utilizzata è la prima.

– Esaurita questa fase di dispacciamento si esegue una nuova iterazione di dispacciamento attivo/reattivo fino a raggiungere una soluzione stabile.

Analisi della sicurezzaAnalisi della sicurezza– Successivamente alla effettuazione del dispacciamento è

necessario verificare che la rete sia in grado di trasportare la potenza che viene prodotta dai singoli gruppi fino ai nodi di carico.

– Questa fase va sotto il nome di controllo della sicurezza del sistema; essa viene effettuata il giorno prima a valle del dispacciamento.

– La verifica di sicurezza viene generalmente effettuata mediante una procedura che va sotto il nome di verifica della sicurezza N-1 del sistema.

– In Italia tale verifica veniva effettuata dal CNC (Centro Nazionale di Controllo)

Il sistema di controllo.Il sistema di controllo.

Programmazione idraulico di breve termine –Obiettivo della programmazione di breve termine...

Documents

Transcript of Programmazione idraulico di breve termine –Obiettivo della programmazione di breve termine...