1 Introduction
A rising interest to investigate the impact of high penetration of renewable generation on frequency control of power system and the capability of delivering frequency support by full Converter Control-Based Generators (CCBGs) is emerged during recent years [1]. Decrease of system inertia and conventional spinning reserve and increase of unpredictable uncertainties imposed by integrating renewables in the generation mix deteriorates the frequency support. By implementation of the virtual synchronous machine behavior [2, 3] through the new grid-supporting convertors of storages and renewables, the frequency support can be improved by means of the provided fast synthetic inertia and timely generation control in severe conditions [4].
In severe frequency decline situations, Under-Frequency Load Shedding (UFLS) programs (generally including event-based and response-based UFLSs) as the last automated measure of power system prevention from collapsing are designed and implemented by Transmission System Operator (TSO) through dynamic off-line simulations in operational planning time-scale. In response-based UFLS programs, the relay settings are usually pre-defined, fixed, and non-robust which may not be a comprehensive solution for the wide range of combinational/cascading events and operational variability and uncertainty imposed by high penetration of renewables [5]. Grid-interactive BESS can support the system operation during inertial and primary response through the inertia and droop parameters tuning of the SPC (Synchronous Power Controller).
Several TSOs have announced new grid codes requiring Electronically Interfaced Resources (EIRs) like energy storages, PVs, HVDC links, and wind power plants to provide frequency response [6]. The mutual interaction of UFLS as a traditional emergency control action has not been investigated enough under high share of renewable energy sources with BESS frequency support. The influence of the fast acting power controller of BESS on UFLS scheme and also mal-operation and readjustment of UFLS settings in high penetration of wind generation are presented in [7] and [8]. In [9], a systematic method is presented for controlling an energy storage system output power for preventing transient load shedding. Existing UFLS schemes suggested for conventional power systems has not been yet incorporated in grid-supporting convertors tuning to provide better inertial and primary frequency response.
A SPC-based Battery Energy Storage System (BESS) [10–12] as a linear state feedback controller could be tuned to control the performance of system frequency response measured by conflicting multi-dimensional performance measures such as minimum RoCoF and frequency, steady state error, and settling time. Good performance of UFLS schemes as another state feedback controller in conjunction with BESS frequency response is of great importance considering the fact that frequency support of the BESS may be harmful to some frequency performance measures or UFLS scheme performance. On the other hand, UFLS actions as a step change in RHS of the swing equation declines the BESS frequency response. Although, UFLS scheme has similar operation horizon as BESS frequency support, they used to be adjusted in different operational planning horizons, which makes the conflicts more complex and inevitable. In this paper, it is shown that the gains of the SPC and UFLS relays parameters need be adjusted in a coordinated manner to compromise between diverse technical and economic performance criteria of the frequency response.
The rest of the paper is constructed as follows: Sects. 2 and 3 provide an investigation of the impact of UFLS scheme and BESS inertia and droop characteristic on frequency response based on the unified discretized system frequency response using UFLS and BESS frequency responses formulation. Simulation results and discussion is presented in Sect. 4, respectively. The conclusion is drawn in Sect. 5.
2 Impact of BESS Inertia and Droop Characteristic on Frequency Response
2.1 Unified Discretized System Frequency Response
the mechanical torque (Tm) (power input),
the electrical torque (Te) (power output), and
the damping torque (Td) [Nm]
![$$ \frac{2H}{{f_{0} }}\frac{df\left( t \right)}{dt} = \sum\nolimits_{i = 1}^{{N_{g} }} {P_{mi} - P_{ei} } $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_Equ1.png)
![$$ \frac{d\Delta f\left( t \right)}{dt} = \frac{{f_{0} }}{2H}\Delta P^{im} \left( t \right) $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_Equ2.png)
![$$ \Delta P^{im} \left( t \right) = \Delta P^{gov} \left( t \right) - \Delta P^{c} \left( t \right) + \Delta P^{sh} \left( t \right) + \Delta P^{wind} \left( t \right) + \Delta P^{PV} \left( t \right) + \Delta P^{ESS} \left( t \right) + \Delta P^{DG} \left( t \right) + \Delta P^{Tie} \left( t \right) - D\Delta f\left( t \right) $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_Equa.png)
![$$ \Delta P^{im} \left( t \right) $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_IEq1.png)
![$$ \frac{d\Delta f\left( t \right)}{dt} $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_IEq2.png)
![$$ \Delta f\left( t \right) $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_IEq3.png)
![$$ (\Delta P^{c} ) $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_IEq4.png)
![$$ (\Delta P^{gov} ) $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_IEq5.png)
![$$ (\Delta P^{wind} , \Delta P^{PV} \left( t \right),\Delta P^{ESS} \left( t \right) ) $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_IEq6.png)
![$$ (\Delta P^{sh} \left( t \right)) $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_IEq7.png)
![$$ (D\Delta f\left( t \right)) $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_IEq8.png)
![$$ \Delta P^{im} $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_IEq9.png)
![$$ \Delta X\left( {n\Delta t} \right) = \Delta X_{n} $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_IEq10.png)
![$$ \Delta f_{n + 1} = \Delta f_{n} + \Delta t\frac{{f_{0} }}{2H}(\Delta P_{n}^{gov} - \Delta P^{c} + \Delta P_{n}^{sh} + \Delta P_{n}^{wind} + \Delta P_{n}^{PV} + \Delta P_{n}^{ESS} - D\Delta f_{n} ) $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_Equ3.png)
![$$ \Delta P_{n + 1}^{gov} = \Delta P_{n}^{gov} + \frac{ - \Delta t}{T}\left( {\Delta P_{n}^{gov} + \frac{{\Delta f_{n} }}{R}} \right) $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_Equ4.png)
The main focus of the rest of the paper is on the BESS influence on system frequency response and UFLS scheme. According to, other terms of (3) except than and
are ignored in the following equations. However, in simulation step, all of the resources contributing in the system frequency response are modelled and reflected.
2.2 System Acceleration Behaviour
![../images/489280_1_En_6_Chapter/489280_1_En_6_Fig1_HTML.png](../images/489280_1_En_6_Chapter/489280_1_En_6_Fig1_HTML.png)
Frequency performance measures illustration.
For a step-wise decrease in from 0 to
at the disturbance time t0, system frequency starts to drop fast with a high RoCoF. Higher RoCoF may activate the RoCoF relays of DGs and loads. The inertial and primary frequency response of the generations and loads try to limit the frequency deviation by increasing the system kinetic energy. At minimum frequency (Nadir frequency) time tm, RoCoF reaches zero and the frequency deviation
starts to decrease. Subsequently,
oscillates and finally stabilizes at the new steady-state frequency
.
![$$ \Delta f $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_IEq18.png)
![$$ \Delta f_{ss} $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_IEq19.png)
![$$ \Delta f - \Delta f_{ss} $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_IEq20.png)
![../images/489280_1_En_6_Chapter/489280_1_En_6_Fig2_HTML.png](../images/489280_1_En_6_Chapter/489280_1_En_6_Fig2_HTML.png)
Acceleration periods of power system frequency response.
![$$ {\text{Accelerating}}:\left( {\Delta f\left( t \right) - \Delta f_{ss} } \right).\frac{d\Delta f\left( t \right)}{dt} > 0 $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_Equ5.png)
![$$ {\text{Decelerating}}:\left( {\Delta f\left( t \right) - \Delta f_{ss} } \right).\frac{d\Delta f\left( t \right)}{dt} < 0 $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_Equb.png)
2.3 BESS as a State Feedback Controller
![$$ \Delta f\left( t \right) = \frac{{\Delta P^{ESS} \left( t \right)}}{{K_{e} }}\quad \;K_{e} = \frac{{ - \left( {2H_{e} s + \frac{1}{{R_{e} }}} \right)}}{{1 + T_{e} s}} $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_Equ6.png)
![$$ H_{e} $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_IEq21.png)
![$$ R_{e} $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_IEq22.png)
![$$ T_{e} $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_IEq23.png)
![$$ - 2H_{e} \frac{d\varDelta f\left( t \right)}{dt} - \frac{1}{{R_{e} }}\varDelta f\left( t \right) = \Delta P^{ESS} \left( t \right) $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_Equ7.png)
![$$ - 2H_{e} \frac{{\Delta f_{n + 1} - \Delta f_{n} }}{\Delta t} - \frac{1}{{R_{e} }}\Delta f_{n} = \Delta P^{ESS} \left( n \right) $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_Equ8.png)
![$$ \Delta P^{ESS} $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_IEq24.png)
![$$ H_{e} $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_IEq25.png)
![$$ R_{e} $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_IEq26.png)
![$$ \Delta f_{n + 1} = \Delta f_{n} + \Delta t\frac{{f_{0} }}{{2\left( {H + H_{e} } \right)}}(\Delta P_{n}^{gov} - \Delta P^{c} + \Delta P_{n}^{sh} - \frac{1}{{R_{e} }}\Delta f_{n} - D\Delta f_{n} ) $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_Equ9.png)
![$$ f,\dot{f} $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_IEq27.png)
![../images/489280_1_En_6_Chapter/489280_1_En_6_Fig3_HTML.png](../images/489280_1_En_6_Chapter/489280_1_En_6_Fig3_HTML.png)
SPC-based BESS as a State feedback control.
2.4 Discussion on BESS Frequency Response
We may state that feedback can cause a system that is originally stable to become unstable. Certainly, feedback is a double-edged sword; when it is improperly used, it can be harmful [14]. Based on the (5) and (6), during accelerating time-intervals, to arrest the rate of change of frequency the system needs inertia emulated by the BESS. As soon as the system starts decelerating, inertia should reduce again to prevent overshoot in frequency response. The unused inertia imposes the frequency to return to the steady-state frequency slower, which deteriorates the overall frequency response.
Additionally, damping should be increased during decelerating phases to help decrease the settling time of the frequency deviation. The additional inertia required during accelerating phases is proportional to the RoCoF and the required additional damping is proportional to the frequency deviation from the steady state. Without adaptive damping and inertia control, the BESS response is not canceled out when the system has sufficiently stabilized which may result in unstable oscillations.
Due to the conflicting influence of the battery inertia and droop responses in accelerating and decelerating periods of frequency response, conflict of interests appears between different performance measures of frequency response with respect to the provided inertia and droop responses of BESS.
3 Response-Driven UFLS as a Non-linear State Feedback Controller
The response driven based load shedding, widely uses Under Frequency-based Load Shedding solutions. The inherent closed loop/ feedback based control scheme in decentralized response driven UFLS system makes them efficient in acting against the disturbances.
As frequency is a very good indicator of Power Mismatch, the operation based on under frequency relays has high control precision and robustness with respect to uncertainties. However, it can have many difficulties and challenges. The setting of the UFLS is highly complex, which involves many parameters including the number of LS stages, percentage of load allowed to be shed, the time delay for each stage, real time topology etc. Together these variables make the UFLS settings a non-linear and multi-dimensional problem. Secondly, UFLS is triggered after the frequency has already declined to certain low values causing a time delay, which makes the solution more reactive in nature and the system stability requires more time.
![$$ \Delta P^{sh} \left( t \right) $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_IEq28.png)
![$$ \Delta P^{sh} \left( n \right) = \sum\nolimits_{k = 1}^{L} {\Delta P_{k} .u\left( {n - n_{k} } \right)} $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_Equ10.png)
![$$ u\left( {n - n_{k} } \right) = 1\;{\text{if}}\;f_{0} + \Delta f_{{n - n_{k} + i}} < f_{th}^{k} \quad i = 0, \ldots ,n_{k} $$](../images/489280_1_En_6_Chapter/489280_1_En_6_Chapter_TeX_Equ11.png)
(10) and (11) could be linearized using binary and big variables [16]. As stated above, the relay logic dictates that the block of load be shed when the corresponding timer exceeds
. When a contingency occurs, for each load-shedding stage k with
, the relay disconnects a block of load after
time step when the frequency trajectory [computed using (3)–(4)] violates the frequency set point for a predetermined time delay.
3.1 Discussion on UFLS and BESS Mutual Effect
![../images/489280_1_En_6_Chapter/489280_1_En_6_Fig4_HTML.png](../images/489280_1_En_6_Chapter/489280_1_En_6_Fig4_HTML.png)
SPC-based BESS and UFLS scheme as State feedback controls control.
UFLS action at is a positive step change in generation-load imbalance which is reflected in system DAEs as a disturbance. Subsequent change in the rocof and frequency values alter the BESS inertial and primary response and reduce the BESS frequency support.
In some situations, however, frequency support of EIRs may be harmful to frequency control as they generate an extra energy in a short period of time, which reduces frequency decay and derivative but could not be maintained over time. Reduction of frequency derivative through emulating inertia via EIRs looks positive at first sight as the reduced frequency derivative triggers less under-frequency relays, but in some conditions, the shedded load is smaller than the amount what the event requires. As the EIRs cannot maintain extra generation over time, the frequency may continue to decay until the shedding of the next load step. It is also possible to activate other system and components protection by resulting in lower steady state frequency.
Based on (10) and (11), load shedding steps are activated, if the frequency reduces temporary or recurrently below a certain UFLS threshold. If the steady state frequency settles at a frequency above the same threshold, load shedding may be considered as nonvital. BESS can provide fast injection of power with limited energy capacity to prevent temporary frequency decreases and unnecessary load shedding. On the other hand, BESS should withdraw from support when the frequency decline is under arrest and additional load shedding would not be activated.
In order to avoid mal operation of UFLS scheme in the presence of the BESS frequency support, UFLS schemes should be incorporated in the controller tuning of the BESS as an economic objective to be optimized in conjunction with other frequency performance measures. In comparison with BESS controller tuning, most power systems operate with predetermined load shedding schemes and use fixed relay settings. Hence, resetting UFLS thresholds or changing load shedding amounts are not often viable options.
Moreover, with the existing supervisory and control capabilities, TSOs are not tendentious to re-design UFLS scheme for a short while. Then, the frequency thresholds, time delays and shed load percentage could be assumed fixed during the optimal gain-tuning of BESS. Using synchrophasor measurements for control of grid interactive energy storage system and UFLS relays [17] provide the possibility of coordinated online control of these resources of inertia and primary frequency response.
4 Studied System
4.1 System Modelling and Characteristics
There are 13 load shedding stages (~68% of the total load) spread over all Cyprus, which are presented in the following Table 1. This UFLS scheme is implemented in the system by definition and setting of the under frequency relays in DIgSILENT PowerFactory. A 50 MW, 100MWh BESS is added to the network equipped by SPC with tunable inertia and droop parameters of 15 and 60 as the maximum gains, as it is presented in [18]. The initial values of 0.5 and 0 are assumed for the SOC and power output states of the BESS, respectively. SOC limits are supposed as 0.2 and 0.8.
Other RHS terms of the swing equation in (2) are enclosed in the system modelling and control like turbine and governors, and Wind turbine standard controllers accorded to Grid Code requirements. Load damping effect are not considered in this study. The coincident disturbance of two generation units with 90 and 60 MW dispatch power is considered to tune the BESS parameters.
4.2 Conflict of Interests
![../images/489280_1_En_6_Chapter/489280_1_En_6_Fig5_HTML.png](../images/489280_1_En_6_Chapter/489280_1_En_6_Fig5_HTML.png)
Effect of UFLS plan and BESS frequency support on frequency response: red (b) and green (c).
- 1.
In (a), the system frequency decreases continuously until the system becomes unstable due to the cascading outage of conventional and renewable resources due to their under frequency relay activation. Primary and inertia support of conventional generation are not enough to arrest the frequency decline.
- 2.By activation of UFLS in (b), the first three blocks of the UFLS plan are activated with 0.2 s delay after the frequency reaches the thresholds (taction = [2.89, 3.15, 3.95]). The non-vital shedding stage is avoided by BESS frequency support. Re-producing the frequency responses for three different combination of activated shedding stages in Fig. 6 (red: stage one, green: stage one and two, blue: all stages) analyses the necessity of the activated stage with respect to steady state frequency. While two first stages are categorized as vital, third stage is non-vital as the frequency settles above 48.8 without the activation of this stage.Fig. 6.
Vitality analysis of activated load shedding stages.
- 3.
Addition of BESS in (c), the first two blocks of the UFLS plan are activated with 0.2 s delay after the frequency reaches the thresholds (taction = [3.24, 3.74]). Although the non-vital shedding stage is avoided by BESS frequency support, other performance measures like RoCoF and minimum frequency are improved and response time and steady state frequency affected harmfully.
Load shedding stages spread over all Cyprus.
Stages | Frequency thresholds | Time delay | Shed Percentage | Shed MW |
---|---|---|---|---|
1 | 49 | 0.2 | 4 | 37.69 |
2 | 48.9 | 0.2 | 4 | 37.69 |
3 | 48.8 | 0.2 | 3 | 28.27 |
4 | 48.7 | 0.2 | 8 | 75.38 |
5 | 48.6 | 0.2 | 4 | 37.69 |
6 | 48.5 | 0.2 | 4 | 37.69 |
7 | 48.4 | 0.2 | 4 | 37.69 |
8 | 48.3 | 0.2 | 1 | 9.42 |
9 | 48.2 | 0.2 | 7 | 65.95 |
10 | 48.1 | 0.2 | 6 | 56.53 |
11 | 48 | 0.2 | 5 | 47.11 |
12 | 47.75 | 0.2 | 9 | 84.80 |
13 | 47.5 | 0.2 | 9 | 84.80 |
Frequency performance measures.
Performance measures | With BESS | Without BESS |
---|---|---|
fNadir | 48.87 | 48.78 |
RoCoFmin | −0.0845 | −0.0859 |
Fss | 49.53 | 49.67 |
Response time | 8.57 | 8.44 |
tNadir | 3.74 | 3.96 |
Pshed | 75.34 | 103.48 |
5 Conclusion
First, it was shown that by considering BESS droop and inertia gains along with the multistage UFLS plan, the overall performance of the load shedding plan and BESS frequency support are mutually influenced by each other in the both favorable and unfair directions. It was also revealed that the BESS droop and inertia gains are affecting frequency performance measures as a double-edged sword regarding system acceleration behaviour. Accordingly, the conflict of interests raised with respect to the UFLS and BESS state feedback controls interference (with similar operation and different planning horizons) should be handled through a gain tuning approach. Enabling coordinated online control of BESS parameters UFLS scheme using synchrophasor measurements provides better inertia and primary frequency response considering demonstrated conflicts of interest.
Acknowledgments
This work was partially supported by the European Commission under project FLEXITRANSTORE—H2020-LCE-2016-2017-SGS-774407 and by the Spanish Ministry of Science under project ENE2017-88889-C2-1-R. Any opinions, findings and conclusions or recommendations expressed in this material are those of the authors and do not necessarily reflect those of the host institutions or funders.
![Creative Commons](../css/cc-by.png)
Open Access This chapter is licensed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license and indicate if changes were made.
The images or other third party material in this chapter are included in the chapter's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the chapter's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.