Sintesi Tecnica: Formulazioni Compatte per l'Unit Commitment con Diversi Livelli di Dettaglio – Parte I

1. Definizione del Problema

Il problema dell'Unit Commitment (UC) è una sfida critica di Programmazione Lineare Intera Mista (MILP) nella gestione dei sistemi elettrici, che mira a determinare il programma operativo ottimale per i generatori (termici, nucleari, rinnovabili, ecc.) su un orizzonte temporale a breve termine (tipicamente 24–48 ore). Il problema consiste nel minimizzare i costi operativi dell'intero sistema soddisfacendo la domanda e rispettando i vincoli specifici delle unità, come i limiti di rampa, i tempi minimi di accensione/spegnimento (up/down time) e le capacità di avvio/arresto.

Nonostante la sua importanza, il problema UC affronta gravi limitazioni computazionali nelle applicazioni su larga scala, in particolare quando viene integrato in modelli di investimento o framework stocastici/robusti. La presenza di variabili binarie (che rappresentano lo stato di impegno dell'unità) causa una crescita esponenziale della complessità del problema. Gli approcci attuali per mitigare ciò spesso prevedono:

• Riduzione della Fedeltà del Modello: Raggruppamento (clustering) di unità simili o uso di giorni rappresentativi, il che può introdurre errori sovrastimando la flessibilità o sottostimando i costi.
• Rilassamento dell'Interezza: Risoluzione della rilassazione Lineare (LP) invece della MILP, il che può portare a grandi gap di ottimalità e soluzioni inammissibili nella pratica.

Esiste una mancanza di una chiara panoramica a priori riguardante quale formulazione offra il miglior equilibrio tra compattezza (lo spazio tra la regione ammissibile MILP e la sua rilassazione LP) e dimensione del modello (numero di variabili e vincoli) per diversi tipi di generatori e livelli di dettaglio.

2. Metodologia

Questo articolo (Parte I) si concentra sullo sviluppo teorico e l'analisi delle formulazioni UC. Gli autori adottano il seguente approccio metodologico:

• Classificazione per Livello di Dettaglio: Gli autori categorizzano i modelli UC in base al numero di variabili binarie per periodo temporale:
  • Modelli 1bin: Contengono solo la variabile di impegno ($u_{gt}$).
  • Modelli 2bin/3bin: Includono variabili aggiuntive per l'avvio ($v_{gt}$) e lo spegnimento ($w_{gt}$).
• Analisi dell'Inviluppo Convesso: L'obiettivo primario è identificare formulazioni basate sull'inviluppo convesso della regione ammissibile. La formulazione più compatta corrisponde all'inviluppo convesso. Tuttavia, poiché trovare l'intero inviluppo convesso richiede spesso un numero esponenziale di vincoli, gli autori si concentrano su sottoinsiemi di facce che:
  1. Descrivano l'inviluppo convesso (o una sua approssimazione stretta) per sottoproblemi specifici (es. rampe, capacità di avvio/arresto).
  2. Mantengano una dimensione del modello simile alla più piccola formulazione valida.
• Dimostrazioni Teoriche: L'articolo fornisce nuove dimostrazioni e intuizioni teoriche riguardanti la compattezza di formulazioni note. Utilizza tecniche come l'eliminazione di Fourier-Motzkin e la programmazione disgiuntiva (metodo di Balas) per derivare e verificare le descrizioni dell'inviluppo convesso per problemi di singola unità su due periodi temporali.
• Limitazione dell'Ambito: Tutte le intuizioni teoriche riguardanti la compattezza sono derivate per un problema operativo di singola unità, escludendo il vincolo di bilancio energetico dell'intero sistema. La funzione obiettivo è generalmente considerata irrilevante per la determinazione dell'inviluppo convesso, a meno di specifiche diverse.

3. Contributi Chiave

L'articolo presenta quattro contributi specifici, di cui i contributi 1 e 2 sono dettagliati nella Parte I:

1. Panoramica Strutturata dei Modelli UC: Gli autori presentano una tassonomia strutturata di modelli UC con diversi livelli di dettaglio (dai limiti di generazione di base alle complesse traiettorie di avvio). Per ogni livello, propongono formulazioni specifiche che soddisfano un compromesso tra compattezza e dimensione.
2. Intuizioni Teoriche e Dimostrazioni:
   • Rampa (Ramping): Vengono fornite nuove dimostrazioni per la compattezza dei vincoli di rampa. L'articolo stabilisce che specifiche formulazioni (Modello I-E-T) costituiscono l'inviluppo convesso a due periodi per le rampe in salita e in discesa.
   • Costi di Avvio/Spegnimento: L'articolo analizza due formulazioni comuni per modellare i costi di avvio e spegnimento (Modello II e Modello II-S). Dimostra che, sebbene il Modello II-S abbia meno vincoli, entrambe le formulazioni definiscono lo stesso spazio di soluzione nella direzione della funzione obiettivo, rendendole ugualmente compatte ai fini dell'ottimizzazione.
   • Capacità: L'articolo presenta le formulazioni dell'inviluppo convesso 1bin e 3bin per le capacità di avvio e spegnimento (output massimo durante l'avvio/spegnimento), provando che combinazioni specifiche di vincoli (es. Modello I-E2 e Modello II-E2) definiscono l'inviluppo convesso per problemi a due periodi.
3. Confronto tra Formulazioni: L'articolo confronta direttamente diverse formulazioni per lo stesso problema (es. 1bin vs 3bin, o diversi set di facce) per evidenziare le differenze teoriche in termini di compattezza e dimensione.
4. Fondamenta per la Parte II: Questo lavoro pone le basi teoriche per la Parte II, che confronta numericamente questi modelli in contesti di investimento e operativi su larga scala.

4. Risultati e Formulazioni

L'articolo definisce diversi modelli specifici, che vanno dai basilari agli estesi:

• Livello I (Modelli 1bin):

  • Modello I: Limiti di generazione di base. Dimostrato essere la formulazione più compatta per il suo ambito.
  • Modello I-E: Aggiunge vincoli di rampa.
  • Modello I-E-T: Una versione più compatta del Modello I-E che include ulteriori facce per la rampa in salita/discesa. Dimostrato essere l'inviluppo convesso a due periodi per le rampe.
  • Modello I-E2: Include rampe e capacità di avvio/spegnimento.
  • Modello I-E2-T: La versione più compatta, dimostrata essere l'inviluppo convesso a due periodi per rampe e capacità.

• Livello II (Modelli 3bin):

  • Modello II: Include variabili e costi di avvio/spegnimento. Dimostrato essere la formulazione più compatta per questo ambito.
  • Modello II-S: Una formulazione alternativa che utilizza variabili di costo continue invece di variabili binarie di avvio. Il Teorema III.4 dimostra che è ugualmente compatto del Modello II nella direzione della funzione obiettivo, sebbene definisca un inviluppo convesso illimitato senza il vincolo superiore $u_t \leq 1$.
  • Modello II-E: Aggiunge vincoli di tempo minimo di accensione/spegnimento (up/down time). Dimostrato essere la formulazione più compatta per questo ambito (dimensione lineare).
  • Modello II-E2: Il modello 3bin più dettagliato, che incorpora rampe, capacità di avvio/spegnimento e tempi minimi di accensione/spegnimento. Il Teorema III.5 dimostra che definisce l'inviluppo convesso per due periodi.
  • Modello II-E3: Estende il Modello II-E2 per includere le traiettorie di avvio e spegnimento (rampe multi-periodo per raggiungere l'output minimo), dimostrato definire l'inviluppo convesso a due periodi.

Risultati Teorici Chiave:

• Per problemi di singola unità, le rilassazioni LP dei Modelli I, II e II-E descrivono l'inviluppo convesso dei rispettivi spazi ammissibili.
• Per problemi a due periodi, i Modelli I-E-T, I-E2-T e II-E2 descrivono l'inviluppo convesso.
• Il rilassamento dell'interezza delle variabili di avvio/spegnimento (da binarie a continue) non cambia l'insieme ammissibile del MIP, ma i solver moderni possono beneficiare del mantenimento delle stesse come intere per le strategie di branching.
• I modelli 3bin possono essere riscritti come modelli 2bin (eliminando la variabile di spegnimento) senza perdere compattezza, sebbene il numero di vincoli possa variare.

5. Significato e Rivendicazioni

L'articolo sostiene che la sua importanza primaria risieda nel fornire un quadro strutturato e teoricamente fondato per la selezione delle formulazioni UC. Definendo modelli con diversi livelli di dettaglio e dimostrando le loro proprietà di compattezza, gli autori permettono ai modellatori di:

• Selezionare il livello di dettaglio più appropriato per generatori specifici o scale di problemi senza sacrificare inutilmente l'efficienza computazionale.
• Comprendere i compromessi tra dimensione del modello e compattezza a priori, anziché affidarsi esclusivamente a test computazionali a posteriori che possono dipendere dal solver e non essere generalizzabili.
• Incorporare dettagli operativi ad alta fedeltà in modelli di investimento su larga scala (come dimostrato nella Parte II) rimanendo entro i limiti computazionali.

Gli autori sottolineano che, sebbene trovare l'intero inviluppo convesso sia spesso intrattabile a causa della crescita esponenziale, le loro formulazioni proposte offrono il "miglior" equilibrio utilizzando sottoinsiemi di facce che sono simili in dimensione alle più piccole formulazioni valide, ma significativamente più compatte rispetto alle rilassazioni standard. Questo lavoro mira a guidare la selezione di modelli UC che siano computazionalmente efficienti ma sufficientemente accurati da evitare la sovrastima della flessibilità del sistema e la sottostima dei costi, fenomeni spesso riscontrati nei modelli di investimento semplificati.