1 Introduction
The share of renewable energy sources (RES) is growing steadily in power
systems. It is essential to facilitate the growth of RES penetration to reduce
the carbon emission from fossil fuel based generators. There are however some
obstacles that limit the applicability of RES. RES are volatile in nature and
forecasting them is subject to uncertainties. Hence, integrating them in the
power system is challenging. Moreover RES are usually decoupled from the
system, and therefore they do not add any inertia to the system. This is par-
ticularly important in small power systems like islands, as they typically suffer
from inertia scarcity, and are therefore more prone to frequency volatility. For
that reason, when integrating RES in such systems, it can be very challenging
to maintain frequency stability in case of contingencies.
To address this issue, researchers have included frequency dynamics in short
term scheduling processes like Unit Commitment (UC) to form a frequency
constrained UC (FCUC), in [1], [2], [3], and etc. The standard (non-frequency
constrained) UC problem can be formulated as a mixed integer linear program-
ming (MILP) problem, which can be solved efficiently using standard solvers.
Unfortunately, the frequency dynamics of a power system is highly nonlinear
and non-convex, complicating how the UC problem can still be formulated as
a MILP problem. There is valuable research work in the literature, addressing
this very issue ([4], [5], [6], and [7]). Frequency dynamics after outages are usu-
ally described by the rate of change of frequency (RoCoF), frequency nadir, and
steady-state frequency. RoCoF and steady-state frequency can be formulated
linearly, but frequency nadir cannot. In previously mentioned studies, the non-
linear constraint on the frequency nadir (derived from the well-known swing
equation) has been simplified or approximated so that it still can be used in
the MILP formulation of UC problem. These formulations are based on simpli-
fying assumptions and usually are computationally demanding. More recently,
data-driven approaches are being introduced to more accurately model the fre-
quency dynamics in the UC problem, instead of relying on analytical methods
([8], [9], [10], [11]). These methods try to estimate the dynamics of the system
accurately, while keeping the solution time of UC reasonably low.
Among the analytical methods, in [1], a linear formulation of inertial re-
sponse and the frequency response of the system is added to the UC problem,
which makes sure that in case of the largest outage, there is enough ancillary
service to prevent under frequency load shedding (UFLS). To linearize frequency
nadir constraint, first-order partial derivatives of its equation with respect to
higher-order non-linear variables are calculated. Then the frequency nadir is
presented by a set of piecewise linear constraints. In [2], different frequency ser-
vices are optimized simultaneously with a stochastic unit commitment (SUC)
approach, targeting low inertia systems that have high levels of RES penetra-
tion. The stochastic model uses scenario trees, generated by the quantile-based
scenario generation method. To linearize frequency nadir, an inner approxima-
tion method is used for one side of the constraint, and for the other side, a
binary expansion is employed to approximate the constraint as a MILP using
the big-M technique. In [3], a stochastic unit commitment approach is intro-
duced for low inertia systems, that includes frequency-related constraints. The
problem considers both the probability of failure events and wind power uncer-
tainty to compute scenario trees for the two-stage SUC problem. An alternative
4