Power System Stability Analysis of the Sicilian Network in the 2050 OSMOSE Project Scenario †
Abstract
:1. Introduction
- SI provided by wind generators;
- FFR provided by BESSs;
- SI and FFR provided by controlled loads (DSR);
- Voltage regulation support provided by RESs (PV and wind generators).
- Large-perturbation angle stability (electromechanical stability): To analyze the response of the power system to the variations in frequency and voltage on a time scale from tens of milliseconds up to tens of seconds. In particular, the electromechanical stability is assessed, taking into account the dynamics of generators and loads and the triggering of protection schemes. The objective is to state if synchronous machines can keep synchronism after a severe transient disturbance [32].
- Small-perturbation angle stability: To evaluate the power system stability when small variations in loads and generation occur, as continuously happens in a real scenario, not necessarily related to a transient disturbance.
- Voltage security: To assess the capability of a power system to keep an acceptable steady voltage at all busses, either under normal conditions or after disturbances.
- Is based on validated realistic scenarios for the year 2050;
- Takes into a complete set of flexibility options, as listed before (SI provided by wind generators; FFR provided by BESSs; SI and FFR provided by controlled loads (DSR); voltage regulation support provided by RESs);
- Uses a detailed network model composed by more than 600 busses, provided by the TSO that owns and manages the modeled system.
2. Network Model and Scenarios
- More than 600 busbars at 400, 230, 150, and 132 kV;
- 441 lines;
- 516 substations;
- 30 large (nominal power higher than 10 MW) static generators, representing the wind and solar plants connected to the HV grid;
- 379 loads, mimic both equivalents and HV loads;
- 72 synchronous machines representing both thermal and hydro plants with their associated controllers.
- The frequency overtakes the threshold of 49.7 Hz and the frequency derivative is lower than −0.2 Hz/s;
- The frequency is lower than 49.5 Hz.
2.1. Scenarios
2.1.1. Data Scenario Application
- Power rating of the cables with the mainland were increased to 2300 MW to meet the Available Transfer Capacity (ATC) data provided.
- An equivalent synchronous machine was installed close to the continental terminal of the link Sicily-continental Italy to mimic the dynamic behavior of the rest of the Italian system. Its rated power equals the sum of the active power of the synchronous machines of the other Italian market zones at the time frame considered, divided by a power factor (0.8), and starting time constant equal to 10 s. This approach is conservative in the sense that the equivalent machine represents only a portion of the generating units of continental Italy able to provide primary frequency control. This machine is equipped with a governor, an automatic voltage regulator, and a power system stabilizer properly tuned.
- 15 Full-Converter Wind Turbines (FCWTs) and 20 Doubly-Fed Induction Generators wind turbines (DFIG) were equipped with a controller able to provide SI (detailed in Section 2.3).
- Already existing HV PV plants ere modelled according to [40], including its LV/MV transformer. The same PV model was installed at the MV busbars of the 285 primary substations to mimic the contribution of the dispersed generation connected to the distribution grids. All PV models are equipped with the overfrequency protections, set according to the Italian standards [38,41].
- A BESS model able to provide FFR (detailed in Section 2.4) was considered. BESSs are installed either near each new wind plant or at the sites of the few synchronous generators not decommissioned. The total amount of storage power was installed according to Task 1.1. These BESS were sized equal to the 20% of the capacity of the already installed generation plant.
- Controlled loads able to provide both SI and FFR (detailed in Section 2.5) were installed near 16 already existing loads considered large enough (at least 13 MW in the Terna base model).
- The 230 kV ring circuit of the Island was doubled, even if no indications regarding the zonal grid reinforcements are given by the other Tasks. This solution was carried out to accommodate the demand expected for 2050.
2.1.2. Typical Generation and Demand Conditions—Dispatching Profiles
- Very low load/very low rotating generation in operation;
- High load/low rotating generation in operation;
- Maximum export/import of areas;
- Operational conditions with a weak network (lines out of service);
- Islanding conditions.
- High Export: represented by 24th of May of MC 4 at 2:00 a.m., characterized by quite high wind production, no photovoltaic, and low load demand.
- High Import: represented by 11th of January of MC 1 at 1:00 a.m., characterized by almost no renewable production and medium/high demand.
- High Load: represented by 27th of May of MC 1 at 4:00 p.m., characterized by high load demand, high wind, and photovoltaic production.
- Island: represented by the 23rd of June of MC 5 at 3:00 p.m.; with high photovoltaic and wind production and the link with the mainland out of service.
- Low Load: represented by the 4th of June of MC2 at 2:00 a.m., characterized by low load and almost zero wind production.
- Lines out of service: it is the Low Load profile with the Favara/Chiaramonte and Caracoli/Sorgente 230 kV lines out of service.
2.2. Wind Generators Models
2.2.1. Full-Converter with Synthetic Inertia
- -
- Two linear Proportional-Integral (PI) regulators (suitably equipped with saturation) are introduced. These two independent controllers are associated with the DC bus stabilization and reactive injection. They provide the direct and quadrature references for the internal current loop.
- -
- An additional signal , introduced on the direct-axis reference and proportional (by the inertia coefficient ) to the frequency approximate derivative, enables the inertial control. Its negative sign is introduced to increase active power injection in case of underfrequency events, exploiting the discharge of the DC capacitor. Indeed, the DC voltage can temporarily vary in response to the system frequency disturbance and allow the inertial response of the converter.
- -
- is the nominal power of the equivalent wind generator exploited;
- -
- is the maximum time-length of the frequency transient supported;
- -
- is the nominal DC voltage;
- -
- is the minimum allowed DC voltage.
- -
- When DC voltage is lower than the minimum allowed value (), discharge is inhibited, and only charging operations are allowed;
- -
- Discharging signals can be sent to the control whenever the voltage is close to its maximum and charging operations are inhibited.
2.2.2. DFIG with Synthetic Inertia
2.3. Battery Energy Storage Systems Model
2.4. DSR
- FCR is not artificially delayed and begins as soon as possible after a frequency deviation;
- If the frequency deviation is equal to or larger than 200 mHz, at least 50% of the full FCR capacity must be delivered at the latest after 15 s; and 100% of the full FCR capacity must be delivered at the latest after 30 s. Moreover, the activation of the full FCR capacity must rise at least linearly from 15 to 30 s;
- If the frequency deviation is smaller than 200 mHz, the activated FCR capacity must be at least proportional with the same time behavior referred in points (A) and (B).
2.5. Methodology
3. Large Perturbation Stability Analysis
- -
- The outage of one or two of the three connections with the mainland (except for the Island DP);
- -
- A 500 MW load-step occurring in continental Italy (load-step in continental Italy is not considered in the Island DP);
- -
- The outage of groups of large generating units.
3.1. Frequency Stability Definitions
- -
- In Normal Operating Conditions (NOC), the frequency must remain within the range 49.9–50.1 Hz;
- -
- For the special case of Sicily disconnected from the mainland, the NOC range is 49.5–50.5 Hz;
- -
- In Emergency Operating Conditions (EOC), the frequency must be kept in the range 47.5–51.5 Hz.
- Strongly stable, if stability is guaranteed and frequency remains within the NOC limits;
- Stable, if stability is guaranteed, but NOC are violated;
- Unstable, if stability is not guaranteed.
3.2. Network Configurations
3.2.1. Configurations with BESS FFR
3.2.2. Configurations with DSR
3.2.3. Simulations Procedure
3.3. Results
4. Small Signal Stability Analysis
4.1. Small Signal Stability Definitions
4.2. Network Configurations
4.2.1. Base Configuration
4.2.2. Configuration with FCWT SI
4.2.3. Combined Configuration
5. Voltage Stability Analysis
5.1. Procedure
5.2. Results
5.2.1. High Import and Island Dispatching Profiles
5.2.2. RES Contribution
5.2.3. Critical Dispatching Profiles
6. Conclusions
Author Contributions
Funding
Institutional Review Board Statement
Informed Consent Statement
Data Availability Statement
Conflicts of Interest
References
- OSMOSE—Optimal System-Mix of Flexibility Solutions for European Flexibility. Available online: https://www.osmose-h2020.eu (accessed on 26 April 2022).
- Liu, H.; Yang, S.; Yuan, X. Inertia Control Strategy of DFIG-Based Wind Turbines Considering Low-Frequency Oscillation Suppression. Energies 2022, 15, 29. [Google Scholar] [CrossRef]
- Wang, J.; Xu, Y.; Wu, X.; Huang, J.; Zhang, X.; Yuan, H. Enhanced Inertial Response Capability of a Variable Wind Energy Conversion System. Energies 2021, 14, 8132. [Google Scholar] [CrossRef]
- Li, Y.; Xu, Z.; Wong, K.P. Advanced Control Strategies of PMSG-Based Wind Turbines for System Inertia Support. IEEE Trans. Power Syst. 2017, 32, 3027–3037. [Google Scholar] [CrossRef]
- Xu, Y.; Chen, P.; Zhang, X.; Yang, D. An Improved Droop Control Scheme of a Doubly-Fed Induction Generator for Various Disturbances. Energies 2021, 14, 7980. [Google Scholar] [CrossRef]
- Bolzoni, A.; Terlizzi, C.; Perini, R. Analytical Design and Modelling of Power Converters Equipped with Synthetic Inertia Control. In Proceedings of the 2018 20th European Conference on Power Electronics and Applications (EPE’18 ECCE Europe), Riga, Latvia, 17–21 September 2018; pp. 1–10. [Google Scholar]
- Fang, J.; Zhang, R.; Li, H.; Tang, Y. Frequency Derivative-based Inertia Enhancement by Grid-Connected Power Converters with a Frequency-Locked-Loop. IEEE Trans. Smart Grid 2018, 10, 4918–4927. [Google Scholar] [CrossRef]
- Conte, F.; Massucco, S.; Silvestro, F. Frequency control services by a building cooling system aggregate. Electr. Power Syst. Res. 2016, 141, 137–146. [Google Scholar] [CrossRef]
- Tindemans, S.H.; Strbac, G. Low-Complexity Decentralized Algorithm for Aggregate Load Control of Thermostatic Loads. IEEE Trans. Ind. Appl. 2021, 57, 987–998. [Google Scholar] [CrossRef]
- Trovato, V.; Sanz, I.M.; Chaudhuri, B.; Strbac, G. Advanced Control of Thermostatic Loads for Rapid Frequency Response in Great Britain. IEEE Trans. Power Syst. 2017, 32, 2106–2117. [Google Scholar] [CrossRef] [Green Version]
- Lin, Y.; Barooah, P.; Meyn, S.; Middelkoop, T. Experimental evaluation of frequency regulation from commercial building HVAC systems. IEEE Trans. Smart Grid 2015, 6, 776–783. [Google Scholar] [CrossRef]
- Cheng, B.; Powell, W.B. Co-Optimizing Battery Storage for the Frequency Regulation and Energy Arbitrage Using Multi-Scale Dynamic Programming. IEEE Trans. Smart Grid 2018, 9, 1997–2005. [Google Scholar] [CrossRef]
- Zhao, H.; Wu, Q.; Hu, S.; Xu, H.; Rasmussen, C.N. Review of energy storage system for wind power integration support. Appl. Energy 2015, 137, 545–553. [Google Scholar] [CrossRef]
- Conte, F.; Massucco, S.; Schiapparelli, G.P.; Silvestro, F. Day-ahead and intra-day planning of integrated BESS-PV systems providing frequency regulation. IEEE Trans. Sustain. Energy 2020, 11, 1797–1806. [Google Scholar] [CrossRef]
- Berizzi, A.; Bovo, C.; Ilea, V.; Merlo, M.; Miotti, A.; Zanellini, F. Decentralized reactive power control of wind power plants. In Proceedings of the 2012 IEEE International Energy Conference and Exhibition, ENERGYCON 2012, Florence, Italy, 9–12 September 2012. [Google Scholar]
- Christakou, K.; Tomozei, D.C.; Bahramipanah, M.; Le Boudec, J.Y.; Paolone, M. Primary voltage control in active distribution networks via broadcast signals: The case of distributed storage. IEEE Trans. Smart Grid 2014, 5, 2314–2325. [Google Scholar] [CrossRef]
- Ilea, V.; Bovo, C.; Falabretti, D.; Merlo, M.; Arrigoni, C.; Bonera, R.; Rodolfi, M. Voltage control methodologies in active distribution networks. Energies 2020, 13, 3293. [Google Scholar] [CrossRef]
- Ma, K.; Yao, T.; Yang, J.; Guan, X. Residential power scheduling for demand response in smart grid. Int. J. Electr. Power Energy Syst. 2016, 78, 320–325. [Google Scholar] [CrossRef]
- Bañales, S.; Dormido, R.; Duro, N. Smart Meters Time Series Clustering for Demand Response Applications in the Context of High Penetration of Renewable Energy Resources. Energies 2021, 14, 3458. [Google Scholar] [CrossRef]
- Wang, J.; Padullaparti, H.; Ding, F.; Baggu, M.; Symko-Davies, M. Voltage Regulation Performance Evaluation of Distributed Energy Resource Management via Advanced Hardware-in-the-Loop Simulation. Energies 2021, 14, 6734. [Google Scholar] [CrossRef]
- Stephen, A.A.; Musasa, K.; Davidson, I.E. Voltage Rise Regulation with a Grid Connected Solar Photovoltaic System. Energies 2021, 14, 7510. [Google Scholar] [CrossRef]
- Szultka, A.; Szultka, S.; Czapp, S.; Karolak, R.; Andrzejewski, M.; Kapitaniak, J.; Kulling, M.; Bonk, J. Voltage Profiles Improvement in a Power Network with PV Energy Sources—Results of a Voltage Regulator Implementation. Energies 2022, 15, 723. [Google Scholar] [CrossRef]
- Ding, Z.; Huang, X.; Liu, Z. Active Exploration by Chance-Constrained Optimization for Voltage Regulation with Reinforcement Learning. Energies 2022, 15, 614. [Google Scholar] [CrossRef]
- Li, W.; Tang, M.; Zhang, X.; Gao, D.; Wang, J. Operation of Distributed Battery Considering Demand Response Using Deep Reinforcement Learning in Grid Edge Control. Energies 2021, 14, 7749. [Google Scholar] [CrossRef]
- Weibezahn, J.; Göke, L.; Orrù, L. OSMOSE D1.1: European Long-Term Scenarios Description. 2019. Available online: https://www.osmose-h2020.eu/download/d1-1-european-long-term-scenarios-description (accessed on 29 March 2022).
- Weibezahn, J.; Göke, L.; Poli, D.; Protard, V.; Bourien, Y.M.; Ažman, G. OSMOSE D1.2: Flexibility Cost and Operational Data Outlook. 2020. Available online: https://www.osmose-h2020.eu/download/d1-2-flexibility-cost-and-operational-data-outlook (accessed on 29 March 2022).
- Bourmaud, J.-Y.; Lhuillier, N.; Orlic, D.; Kostic, M.; Heggarty, T.; Grisey, N. OSMOSE D1.3: Optimal Mix of Flexibility. 2022. Available online: https://www.osmose-h2020.eu/wp-content/uploads/2022/04/D1.3-Optimal-Mix-of-Flexibility-1.pdf (accessed on 29 March 2022).
- ENSIEL. Available online: https://www.consorzioensiel.it (accessed on 28 February 2022).
- Heggarty, T.; Bourmaud, J.Y.; Girard, R.; Kariniotakis, G. Quantifying power system flexibility provision. Appl. Energy 2020, 279, 115852. [Google Scholar] [CrossRef]
- Berizzi, A.; Ilea, V.; Vicario, A.; Conte, F.; Massucco, S.; Adu, J.A.; Nucci, C.A.; Pontecorvo, T. Stability analysis of the OSMOSE scenarios: Main findings, problems, and solutions adopted. In Proceedings of the 2021 AEIT International Annual Conference (AEIT), Milan, Italy, 4–8 October 2021; pp. 1–6. [Google Scholar]
- DIgSILENT GmbH. DIgSILENT PowerFactory, User Manual; Version 20; DIgSILENT GmbH: Gomaringen, Germany, 2019. [Google Scholar]
- Kundur, P. Power System Stability And Control; McGraw-Hill: New York, NY, USA, 1994. [Google Scholar]
- Adrees, A.; Milanović, J.V.; Mancarella, P. Effect of inertia heterogeneity on frequency dynamics of low-inertia power systems. IET Gener. Transm. Distrib. 2019, 13, 2951–2958. [Google Scholar] [CrossRef]
- IEEE/NERC. Task Force on Short-Circuit and System Performance Impact of Inverter Based Generation. In Impact of Inverter Based Generation on Bulk Power System Dynamics and Short-Circuit Performance; IEEE: Piscataway, NJ, USA, 2018. [Google Scholar]
- Eftekharnejad, S.; Vittal, V.; Heydt, G.T.; Keel, B.; Loehr, J. Small signal stability assessment of power systems with increased penetration of photovoltaic generation: A case study. IEEE Trans. Sustain. Energy 2013, 4, 960–967. [Google Scholar] [CrossRef]
- Map of the Sicilian 380–220 kV System. Available online: http://www.regione.sicilia.it (accessed on 28 February 2022).
- EC Directive (EU). 2018/2001 of the European Parliament and of the Council of 11 December 2018 on the promotion of the use of energy from renewable sources. Off. J. Eur. Union 2018, 2001, 1–128. [Google Scholar]
- TERNA S.p.A. Grid Code Annexes 11, “Criteri Generali per la Taratura delle Protezione delle Reti a Tensione Uguale o Superiore a 110 kV”. 2018. Available online: https://download.terna.it/terna/0000/0105/28.pdf (accessed on 29 March 2022). (In Italian).
- TERNA S.p.A. GAUDI’ Portal. Available online: https://www.terna.it/it/sistema-elettrico/gaudi (accessed on 28 February 2022).
- Western Electricity Coordinating Council Modeling and Validation Work Group. WECC Renewable Energy Modeling Task Force WECC Solar Plant Dynamic Modeling Guidelines; World Environment Center: Washington, DC, USA, 2014. [Google Scholar]
- TERNA S.p.A. Grid Code Annexes 12, “Criteri di Taratura dei Relè di Frequenza del Sistema Elettrico e Piano di Alleggerimento”. 2018. Available online: https://download.terna.it/terna/0000/1127/88.pdf (accessed on 29 March 2022). (In Italian).
- Berizzi, A.; Bosisio, A.; Ilea, V.; Marchesini, D.; Perini, R.; Vicario, A. Analysis of Synthetic Inertia Strategies from Wind Turbines for Large System Stability. In IEEE Transactions on Industry Applications; IEEE: New York, NY, USA, 2022. [Google Scholar] [CrossRef]
- Wang, B.; Ni, J.; Geng, J.; Lu, Y.; Dong, X. Arc flash fault detection in wind farm collection feeders based on current waveform analysis. J. Mod. Power Syst. Clean Energy 2017, 5, 211–219. [Google Scholar] [CrossRef] [Green Version]
- Bolzoni, A.; Perini, R. Experimental validation of a novel angular estimator for synthetic inertia support under disturbed network conditions. In Proceedings of the 2019 21st European Conference on Power Electronics and Applications (EPE’19 ECCE Europe), Genova, Italy, 3–5 September 2019; pp. 1–10. [Google Scholar] [CrossRef] [Green Version]
- Hu, J.; Huang, Y.; Wang, D.; Yuan, H.; Yuan, X. Modeling of grid-connected DFIG-based wind turbines for dc-link voltage stability analysis. IEEE Trans. Sustain. Energy 2015, 6, 1325–1336. [Google Scholar] [CrossRef]
- Adu, J.A.; Napolitano, F.; Nucci, C.A.; Diego Rios Penaloza, J.; Tossani, F. A DC-Link Voltage Control Strategy for Fast Frequency Response Support. In Proceedings of the 2020 IEEE 20th Mediterranean Electrotechnical Conference ( MELECON), Palermo, Italy, 16–18 June 2020; pp. 470–475. [Google Scholar] [CrossRef]
- Adu, J.A.; Rios Penaloza, J.D.; Napolitano, F.; Tossani, F. Virtual Inertia in a Microgrid with Renewable Generation and a Battery Energy Storage System in Islanding Transition. In Proceedings of the SyNERGY MED 2019—1st International Conference on Energy Transition in the Mediterranean Area, Cagliari, Italy, 28–30 May 2019. [Google Scholar]
- ENTSO-E. Demand Connection Code. Available online: https://www.entsoe.eu/major-projects/network-code-development/demand-connection (accessed on 20 June 2020).
- TERNA S.p.A. Codice di trasmissione, Dispacciamento, Sviluppo e Sicurezza Della Rete. 2018. Available online: https://download.terna.it/terna/0000/0107/31.pdf (accessed on 29 March 2022). (In Italian).
- ENTSO-E. Frequency Stability Evaluation Criteria for the Synchronous Zone of Continental Europe: Requirements and Impacting Factors; ENTSO-E: Brussels, Belgium, 2016. [Google Scholar]
- Machowski, J.; Bialek, J.W.; Bumby, J.R. Power System Dynamics: Stability and Control; John Wiley & Sons: Hoboken, NJ, USA, 2008; ISBN 9780470725580. [Google Scholar]
- Gonzalez-Longatt, F.; Rueda, J. PowerFactory Applications for Power System Analysis; Springer: Berlin/Heidelberg, Germany, 2014; ISBN 978-3-319-12957-0. [Google Scholar]
- TERNA S.p.A. Grid Code Annexes 17, “Centrali Eoliche—Condizioni Generali di Connessione Alle Reti AT. Sistemi di Protezione, Regolazione e Controllo”. 2019. Available online: https://download.terna.it/terna/0000/0105/34.pdf (accessed on 29 March 2022). (In Italian).
- TERNA S.p.A. Grid Code Annexes 68, “Centrali Fotovoltaiche—Condizioni Generali di Connessione Alle Reti AT. Sistemi di Protezione, Regolazione e Controllo”. 2019. Available online: https://download.terna.it/terna/0000/0105/84.pdf (accessed on 29 March 2022). (In Italian).
- CEI 0–16. Regola Tecnica di Riferimento per la Connessione di Utenti Attivi e Passivi Alle Reti AT ed MT Delle Imprese Distributrici di Energia Elettrica; Comitato Elettrico Italiano: Milano, Italy, 2019. (In Italian) [Google Scholar]
Plant | Owner | Type | Rating (MW) |
---|---|---|---|
Anapo (SR) | Enel | Storage Hydro | 500 |
Augusta (SR) | Enel | Fossil Fuel | 210 |
Priolo Gargallo, Nuce Nord (SR) | ERG | CCGT * | 480 |
Priolo Gargallo, ISAB energy (SR) | ERG | IGCC * | 528 |
Priolo Gargallo Archimede (SR) | Enel | CCGT | 750 |
Trapani (TP) | E.On | OCGT * | 169 |
Gela (CL) | Eni | Fossil Fuel | 260 |
San Filippo del Mela (ME) | A2A Energie Future | Fossil Fuel | 1280 |
Milazzo (ME) | Termica Milazzo | CCGT | 365 |
Termini Imerese (PA) | Enel | CCGT | 1340 |
Threshold | Starting Frequency (Hz) | Frequency Derivative Threshold (Hz/s) | Pure Frequency Threshold (Hz) | Percentage of Shed Load |
---|---|---|---|---|
1 | 49.3 | −0.3 | 49.0 | 9% |
2 | 49.2 | −0.6 | 48.9 | 8% |
3 | 49.1 | −0.9 | 48.8 | 7% |
4 | 49.1 | −1.2 | 48.7 | 7% |
Technology | Total Installed (MW) | MV Connected (%) | HV Connected (%) |
---|---|---|---|
Photovoltaic | 1422 | 97.00 | 3.00 |
Wind | 1887 | 6.00 | 94.00 |
Hydro | 274 | 4.00 | 96.00 |
Zone | Battery | PV | Hydro | Wind | Waste | Gas | P2G |
---|---|---|---|---|---|---|---|
56IT | 1572 | 6075 | 313 | 6360 | 566 | 162 | 1947 |
Location | Identification Name | Rating (MVA) |
---|---|---|
Contrasto | CNTP | 24 |
Dittaino | DITP | 23 |
Priolo Gargallo | EGNP | 576 |
Città Giardino | EGSP | 90 |
Augusta | ESSP | 80 |
Priolo Gargallo | ISBP | 344 |
Paternò | PATP | 9 |
Priolo Gargallo | PRGP | 658 |
Milazzo | TEMP | 185 |
Termini Imerese | TIMP | 946 |
Troina | TROP | 14 |
Dispatching Profile | High Export | High Import | High Load | Island | Low Load | Lines Out of Service |
---|---|---|---|---|---|---|
Time | 24 May 2050 02:00 | 11 January 2050 01:00 | 27 May 2050 16:00 | 23 June 2030 15:00 | 4 June 2050 02:00 | 4 June 2050 02:00 |
MC Year | 4 | 1 | 1 | 5 | 2 | 2 |
Zonal export (MW) | 725 | −1188 | 0 | 0 | −578 | −578 |
Loads | ||||||
Total Load (MW) | 1837 | 3218 | 4099 | 4294 | 1514 | 1514 |
Traditional load (MW) | 1837 | 3218 | 2318 | 2500 | 1514 | 1514 |
DSR Electrical Vehicle (MW) | 0 | 0 | 0 | 0 | 0 | 0 |
Pumping (MW) | 0 | 0 | 0 | 0 | 0 | 0 |
Electrolyser (MW) | 0 | 0 | 1781 | 1794 | 0 | 0 |
Battery Storage (MW) | 0 | 0 | 0 | 0 | 0 | 0 |
DSR Heat Pump (MW) | 0 | 0 | 0 | 0 | 0 | 0 |
Generation | ||||||
Total Generation (MW) | 2562 | 2030 | 4099 | 4294 | 936 | 936 |
ROR (MW) | 30 | 11 | 30 | 25 | 25 | 25 |
WIND (MW) | 1480 | 253 | 2059 | 1124 | 195 | 195 |
SOLAR (MW) | 0 | 0 | 1211 | 2429 | 0 | 0 |
NUCLEAR (MW) | 0 | 0 | 0 | 0 | 0 | 0 |
COAL (MW) | 0 | 0 | 0 | 0 | 0 | 0 |
GAS (MW) | 0 | 0 | 0 | 0 | 0 | 0 |
BATTERY (MW) | 786 | 0 | 0 | 0 | 0 | 0 |
PSP (MW) | 0 | 0 | 0 | 0 | 0 | 0 |
P2G (MW) | 0 | 1500 | 533 | 450 | 450 | 450 |
CHP (MW) | 0 | 0 | 0 | 0 | 0 | 0 |
BIOENERGY (MW) | 266 | 266 | 266 | 266 | 266 | 266 |
H. STOR (MW) | 0 | 0 | 0 | 0 | 0 | 0 |
SPIL. ENRG (MW) | 0 | 0 | 0 | 0 | 0 | 0 |
UNSP. ENRG (MW) | 0 | 0 | 0 | 0 | 0 | 0 |
Continental Italy Equivalent Generator | ||||||
NOMINAL POWER (MVA) | 10,411 | 1493 | 23,053 | 9607 | 13,533 | 13,533 |
Active Power (MW) | ||||||
---|---|---|---|---|---|---|
Load | High Export | High Import | High Load | Island | Low Load | Lines Out of Service |
Load 1 | 42.8 | 42.8 | 59.92 | 59.92 | 21.4 | 21.4 |
Load 2 | 13.8 | 27.6 | 38.64 | 38.64 | 13.8 | 13.8 |
Load 3 | 14.5 | 29 | 40.6 | 40.6 | 14.5 | 14.5 |
Load 4 | 13.9 | 27.8 | 38.92 | 38.92 | 13.9 | 13.9 |
Load 5 | 13.3 | 26.6 | 37.24 | 37.24 | 13.3 | 13.3 |
Load 6 | 15 | 30 | 42 | 42 | 15 | 15 |
Load 7 | 126.6 | 189.9 | 177.24 | 177.24 | 0 | 0 |
Load 8 | 13.5 | 27 | 37.8 | 37.8 | 13.5 | 13.5 |
Load 9 | 16.9 | 33.8 | 47.32 | 47.32 | 16.9 | 16.9 |
Load 10 | 15 | 30 | 42 | 42 | 15 | 15 |
Load 11 | 13.4 | 26.8 | 37.52 | 37.52 | 13.4 | 13.4 |
Load 12 | 52.4 | 52.4 | 73.36 | 73.36 | 26.2 | 26.2 |
Load 13 | 13.8 | 27.6 | 38.64 | 38.64 | 13.8 | 13.8 |
Load 14 | 13.6 | 27.2 | 38.08 | 38.08 | 13.6 | 13.6 |
Load 15 | 14.8 | 29.6 | 41.44 | 41.44 | 14.8 | 14.8 |
Load 16 | 13.3 | 26.6 | 37.24 | 37.24 | 13.3 | 13.3 |
TOTALS | 406.6 | 654.7 | 827.96 | 827.96 | 232.4 | 232.4 |
Parameter | Symbol | Value |
---|---|---|
Single cell capacity | 50 Ah | |
Voltage of full cells | 15 V | |
Voltage of empty cells | 12 V | |
Number of series cells | 60 | |
Battery nominal voltage | 0.9 kV | |
Converter AC side nominal voltage | 0.4 kV | |
Cells internal resistance | 0.001 Ω |
Event | Base Configuration | Best BESSs—DSR Configuration | ||||
---|---|---|---|---|---|---|
Stable? | Strongly Stable? | Configuration | FCR | Stable? | Strongly Stable? | |
High Export | ||||||
Load-step in continental Italy (500 MW) | YES | NO | BESS-1 DSR-25 | 655 MW | YES | YES |
Outage of 1 link with continental Italy | YES | YES | Not required | - | - | - |
Outage of 2 links with continental Italy | YES | YES | Not required | - | - | - |
Generation outage (220 MW) | YES | YES | Not required | - | - | - |
High Import | ||||||
Load-step in continental Italy (500 MW) | YES | NO | BESS-5 DSR-15 | 203 MW | YES | YES |
Outage of 1 link with continental Italy | YES | YES | Not required | - | - | - |
Outage of 2 links with continental Italy | NO | NO | BESS-1 | 522 MW | YES | YES |
Generation outage (321 MW) | YES | YES | Not required | - | - | - |
Generation outage (569 MW) | YES | NO | BESS-3 DSR-25 | 338 MW | YES | YES |
Generation outage (817 MW) | YES | NO | BESS-1 DSR-35 | 751 MW | YES | NO |
High Load | ||||||
Load-step in continental Italy (500 MW) | YES | YES | Not required | - | - | - |
Outage of 1 link with continental Italy | YES | YES | Not required | - | - | - |
Outage of 2 links with continental Italy | YES | YES | Not required | - | - | - |
Generation outage (337 MW) | YES | YES | Not required | - | - | - |
Island | ||||||
Generation outage (337 MW) | YES | NO | BESS-3 | 174 MW | YES | YES |
Low Load | ||||||
Load-step in continental Italy (500 MW) | YES | NO | BESS-1 DSR-30 | 362 MW | YES | YES |
Outage of 1 link with continental Italy | YES | YES | Not required | - | - | - |
Outage of 2 links with continental Italy | YES | YES | Not required | - | - | - |
Generation outage (337 MW) | YES | NO | BESS-1 | 292 MW | YES | YES |
Lines out of service | ||||||
Load-step in continental Italy (500 MW) | YES | NO | BESS-1 DSR-30 | 362 MW | YES | NO |
Outage of 1 link with continental Italy | YES | YES | Not required | - | - | - |
Outage of 2 links with continental Italy | YES | YES | Not required | - | - | - |
Generation outage (337 MW) | YES | NO | BESS-3 DSR-30 | 362 MW | YES | YES |
Mode | Eigenvalues | Frequency (Hz) | Damping (%) | Involved Plants |
---|---|---|---|---|
M1 | −0.27 ± 4.81j | 0.77 | 5.66 | TIMP, PRGP, SLACK |
D1 | −0.67 ± 11.29j | 1.80 | 5.90 | DFIGs in ANPP |
M2 | −0.59 ± 8.18j | 1.30 | 7.24 | EGSP, PRGP |
M3 | −0.88 ± 9.98j | 1.59 | 8.83 | EGSP |
M4 | −0.89 ± 9.97j | 1.59 | 8.90 | EGSP |
M5 | −1.04 ± 10.44j | 1.66 | 9.91 | DITP, ESSP, PRGP |
DFIG Active Power Setpoint (MW) | Frequency (Hz) | Damping (%) |
---|---|---|
0 | 1.80 | 6.1 |
0.1 | 1.80 | 5.9 |
0.2 | ✓ | ✓ |
0.3 | ✓ | ✓ |
Mode | Eigenvalues | Frequency (Hz) | Damping (%) | Involved Plants |
---|---|---|---|---|
M1 | −0.19 ± 4.95j | 0.79 | 3.88 | TIMP, PRGP, SLACK |
M6 | −0.80 ± 10.86j | 1.73 | 7.38 | CNTP, PRGP |
M7 | −0.85 ± 11.09j | 1.77 | 7.67 | PATP, CNTP |
Mode | Eigenvalues | Frequency (Hz) | Damping (%) | Involved Plants |
---|---|---|---|---|
M1 | −0.05 ± 4.73j | 0.75 | 1.08 | TIMP, PRGP, SLACK |
M6 | −0.86 ± 10.82j | 1.72 | 7.95 | CNTP, PRGP |
M7 | −0.92 ± 11.05j | 1.76 | 8.26 | PATP, CNTP |
M1 | D1 | M2 | M3 | M4 | M5 | |||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
f (Hz) | ζ (%) | f (Hz) | ζ (%) | f (Hz) | ζ (%) | f (Hz) | ζ (%) | f (Hz) | ζ (%) | f (Hz) | ζ (%) | |
0 | 0.77 | 5.66 | 1.80 | 5.9 | 1.3 | 7.24 | 1.59 | 8.83 | 1.59 | 8.90 | 1.66 | 9.91 |
10 | 0.77 | 5.67 | 1.80 | 5.9 | 1.3 | 7.23 | 1.59 | 8.83 | 1.59 | 8.90 | 1.66 | 9.91 |
30 | 0.77 | 7.39 | 1.80 | 5.9 | 1.3 | 7.07 | 1.59 | 8.83 | 1.59 | 8.90 | 1.66 | 9.69 |
50 | 0.78 | 7.35 | 1.80 | 5.9 | 1.3 | 7.03 | 1.59 | 8.83 | 1.59 | 8.90 | 1.67 | 9.66 |
M1 | M6 | M7 | ||||
---|---|---|---|---|---|---|
f (Hz) | ζ (%) | f (Hz) | ζ (%) | f (Hz) | ζ (%) | |
0 | 0.79 | 3.88 | 1.73 | 7.38 | 1.77 | 7.67 |
10 | 0.79 | 4.37 | 1.73 | 7.34 | 1.77 | 7.67 |
30 | 0.79 | 5.18 | 1.73 | 7.30 | 1.77 | 7.67 |
50 | 0.79 | 5.37 | 1.73 | 7.28 | 1.77 | 7.67 |
M1 | M6 | M7 | ||||
---|---|---|---|---|---|---|
f (Hz) | ζ (%) | f (Hz) | ζ (%) | f (Hz) | ζ (%) | |
0 | 0.75 | 1.08 | 1.72 | 7.95 | 1.76 | 8.26 |
10 | 0.75 | 3.00 | 1.72 | 7.92 | 1.76 | 8.27 |
30 | 0.77 | 5.17 | 1.73 | 7.88 | 1.76 | 8.27 |
50 | 0.78 | 4.85 | 1.73 | 7.87 | 1.76 | 8.27 |
Case | M1 | D1 | M2 | M3 | M4 | M5 | ||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
f (Hz) | ζ (%) | f (Hz) | ζ (%) | f (Hz) | ζ (%) | f (Hz) | ζ (%) | f (Hz) | ζ (%) | f (Hz) | ζ (%) | |
Base | 0.77 | 5.66 | 1.80 | 5.9 | 1.3 | 7.24 | 1.59 | 8.83 | 1.59 | 8.90 | 1.66 | 9.91 |
Comb. | ✓ | ✓ | 1.80 | 5.9 | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ |
Case | M1 | M6 | M7 | |||
---|---|---|---|---|---|---|
f (Hz) | ζ (%) | f (Hz) | ζ (%) | f (Hz) | ζ (%) | |
Base | 0.79 | 3.88 | 1.73 | 7.38 | 1.77 | 7.67 |
Combined | ✓ | ✓ | 1.73 | 7.30 | 1.77 | 7.68 |
Case | M1 | M6 | M7 | |||
---|---|---|---|---|---|---|
f (Hz) | ζ (%) | f (Hz) | ζ (%) | f (Hz) | ζ (%) | |
Base | 0.75 | 1.08 | 1.72 | 7.95 | 1.76 | 8.26 |
Combined | ✓ | ✓ | 1.73 | 7.88 | 1.76 | 8.27 |
Dispatching Profile | Initial Load (MW) | Final Load (MW) | Loadability (MW) | Loadability (%) |
---|---|---|---|---|
High Export | 1835 | 3717 | 1882 | 102.5 |
High Import | 3220 | 3596 | 376 | 11.7 |
High Load | 4088 | 5498 | 1410 | 34.5 |
Low Load | 1510 | 3346 | 1836 | 121.5 |
Island | 4289 | 4857 | 568 | 13.2 |
Lines out of service | 1510 | 3225 | 1715 | 113.5 |
Dispatching Profile | Initial Load (MW) | Loadability (MW) | Loadability (%) | Loadability Base Case (%) |
---|---|---|---|---|
High Export | 4592 | 2757 | 150.2 | 102.5 |
High Import | 3993 | 773 | 24.0 | 11.7 |
High Load | 5702 | 1614 | 40.0 | 34.5 |
Low Load | 4887 | 3376 | 223.5 | 121.5 |
Island | 5791 | 1502 | 35.0 | 13.2 |
Lines out of service | 4615 | 3105 | 205.5 | 113.5 |
Dispatching Profile | Final Load (MW) | Loadability (MW) | Loadability (%) | Loadability without Priolo (%) |
---|---|---|---|---|
High Import | 5144 | 1924 | 60.0 | 24.0 |
Island | 6007 | 1718 | 40.0 | 35.0 |
Dispatching Profile | Final Load (MW) | Loadability (MW) | Loadability (%) | Loadability without Priolo (%) |
---|---|---|---|---|
High Import | 5625 | 2405 | 74.6 | 24.0 |
Island | 6036 | 1747 | 41.0 | 35.0 |
Type of Stability | Is Stability Guaranteed? | Does the Operational Condition Need Specific Additional Flexibilitiy Tools? | Solutions Identified/Notes |
---|---|---|---|
High Export DP | |||
Large perturbation angle stability | YES | YES | DSR 1 and BESS FCR 2 are required to guarantee the normal operating conditions |
Small perturbation angle stability | YES | NO | - |
Voltage stability | YES | YES | Extend the reactive capability of RES units |
High Import DP | |||
Large perturbation angle stability | NO | YES | DSR and BESS FCR are required |
Small perturbation angle stability | YES | YES | SI 3 by FCWTs 4 and PSS tuning are required |
Voltage stability | NO | YES | Extend the reactive capability of RES units + Local compensation in Priolo substation |
High Load DP | |||
Large perturbation angle stability | YES | NO | - |
Small perturbation angle stability | YES | NO | - |
Voltage stability | NO | YES | Extend the reactive capability of RES units |
Island DP | |||
Large perturbation angle stability | NO | YES | DSR and BESS FCR are required |
Small perturbation angle stability | YES | NO | - |
Voltage stability | NO | YES | Extend the reactive capability of RES units + Local compensation in Priolo substation |
Low Load DP | |||
Large perturbation angle stability | NO | YES | DSR and BESS FCR are required |
Small perturbation angle stability | NO | YES | SI by FCWTs and PSS tuning are required |
Voltage stability | YES | YES | Extend the reactive capability of RES units |
Lines out of service DP | |||
Large perturbation angle stability | NO | YES | DSR and BESS FCR are required |
Small perturbation angle stability | NO | YES | SI by FCWTs and PSS tuning are required |
Voltage stability | YES | YES | Extend the reactive capability of RES units |
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Adu, J.A.; Berizzi, A.; Conte, F.; D’Agostino, F.; Ilea, V.; Napolitano, F.; Pontecorvo, T.; Vicario, A. Power System Stability Analysis of the Sicilian Network in the 2050 OSMOSE Project Scenario. Energies 2022, 15, 3517. https://doi.org/10.3390/en15103517
Adu JA, Berizzi A, Conte F, D’Agostino F, Ilea V, Napolitano F, Pontecorvo T, Vicario A. Power System Stability Analysis of the Sicilian Network in the 2050 OSMOSE Project Scenario. Energies. 2022; 15(10):3517. https://doi.org/10.3390/en15103517
Chicago/Turabian StyleAdu, James Amankwah, Alberto Berizzi, Francesco Conte, Fabio D’Agostino, Valentin Ilea, Fabio Napolitano, Tadeo Pontecorvo, and Andrea Vicario. 2022. "Power System Stability Analysis of the Sicilian Network in the 2050 OSMOSE Project Scenario" Energies 15, no. 10: 3517. https://doi.org/10.3390/en15103517
APA StyleAdu, J. A., Berizzi, A., Conte, F., D’Agostino, F., Ilea, V., Napolitano, F., Pontecorvo, T., & Vicario, A. (2022). Power System Stability Analysis of the Sicilian Network in the 2050 OSMOSE Project Scenario. Energies, 15(10), 3517. https://doi.org/10.3390/en15103517