1. Introduction
In the last decade, electric distribution systems are undergoing deep changes to allow large penetration of Distributed Generation (DG) exploiting Renewable Energy Sources (RESs). Their configuration has changes from passive to Active Distribution Networks (ADNs) and new Distributed Energy Sources (DERs) are being connected, which include, in addition to DG, Battery Energy Storage Systems (BESSs), controllable loads and Electric Vehicles (EVs).
A key issue in ADNs is the regulation of the voltage profiles which are no longer monotonic along the feeders and can violate the quality limits of over and under voltage. In distribution systems voltage regulation is traditionally poorly automated and essentially centralized. Then, the new ADNs require all the DERs to support the voltage regulation. In particular, it is possible to act on the inverters interfacing the DERs to grid so as to vary their reactive power injections with no impact on the operation of the DC bus of the DER and, consequently, on the RES generator. However due to the high
ratio of distribution lines the control of the reactive power injection is not enough effective to assure an adequate voltage regulation. In this view, Active Power Curtailment (APC) has been proposed as an additional tool for voltage regulation. In [
1] the control system prioritizes the use of reactive power, while APC is performed only as a last resort; the controllers switch to active power curtailment only if all the reactive power capability has been exploited. In [
2] a method that includes APC in coordination with reactive power control for voltage control in distribution networks is presented: APC is applied if the inverter rated apparent power is reached.
However, the approaches exploiting APC generally present two drawbacks. The first one is the non linear behavior of the control systems which often adopt a switching strategy between reactive power control and APC. To overcome this issue, the proposed voltage control methods typically use different control layers, that are adequately coordinated to manage both reactive power and APC [
3,
4,
5]. The second important drawback is the economical impact of APC: it wastes energy produced from RESs and reduces revenues [
6]. To reduce the economical impact, some methods use the voltage-active sensitivity factors to define the optimal APC strategy with respect to some chosen objectives, such as acting with priority on the most voltage-effective PV system [
7] so as to limit the overall curtailment. To really overcome this problem, BESS should be installed to store excessive PV produced energy during peak generation periods and then use it to supply active power for voltage regulation [
8]. But the BESS equipment is usually expensive, making it difficult to be justified from the cost-benefit point of views [
9].
Recently, the asset of distribution systems with DERs is significantly changing because of the environmental requirements that push towards the complete substitution of conventional energy sources with RESs. To allow a further increase of RESs in distribution systems two main innovations are being introduced:
new Ancillary Services (ASs) provided by DERs at ADN level are being defined [
10] and commercialized in the deregulated markets [
9];
new market agents can aggregate [
11] to provide ASs at transmission system level, as well as new entities, such as energy communities [
12], are being introduced to promote the power balance and the AS provision at distribution system level.
Such innovations are favouring the installation of BESSs in conjunction with RES generators because ASs can provide additional incomes that justify the cost of BESSs. For this reason, the present paper focuses attention on DERs composed of DG exploiting RESs, equipped with BESS and interfaced to the ADN by an inverter. Such a type of DER is commonly equipped with a local controller that guarantees voltage regulation at node of connection with the network.
As well known, any local controller interacts with each other through the distribution lines of the ADN when they act on multiple DERs connected to the same feeder [
13]. This mutual interaction leads to system operating conflicts and also to voltage instability [
14,
15,
16]. The techniques employed to address this problem are categorized in two main approaches [
17,
18]. The first is based on centralized and/or distributed techniques which guarantee optimal solutions to the voltage control problem, for example, in term of power loss minimization [
19,
20,
21,
22]. A general issue with the aforementioned techniques is that they are based on field network’s measurements in a closed loop real time control system [
23,
24]. However, the implementation of the communication-based control strategies may not be practical in active distribution networks that lack of information infrastructures or the whole control scheme fails to work if a communication failure occurs [
25]. The main drawback of the centralized solutions is the computational cost to handle a big amount of data measured by the ADN. On the contrary, distributed control algorithms require less information, coming from the neighboring DERs, with respect to centralized solutions but may suffer from a slow convergence rate [
26]. The second approach relies on the use of decentralized techniques. Generally they do not give optimal solutions since use only local information and measurements. However, as previously noted, due to the requirement of a large amount of information exchanges and increased needed computing capacity, decentralized techniques are usually preferred since they can be implemented in the existing distribution networks easily. Furthermore, they are less complex and more reliable compared to centralized/distributed techniques [
27].
Keeping the simplicity of a completely decentralized architecture, our contribution to the voltage control problem is the design that we present in this paper applicable to local control of DERs composed of DG exploiting RESs and equipped with BESS. Each local controller acts on the reactive power of the inverter and, if necessary, provides an additional AS by varying the active power injected by the DER. In particular, the AS is defined by its droop characteristic and by the range of power variation allowed by the BESS. Moreover, the local control scheme should both avoid the introduction of any switching mechanism between reactive and active power control and mitigate the interaction’s level.
To pursue this objective, in the proposed approach each DER is equipped with a double IMC-based PI local control loop, which regulates the voltage and active power to the desired set-point, respectively. Moreover, the two control loops are coupled by a filtered droop control law which realizes the AS. Based on the Relative Gain Matrix (RGA) matrix of the adopted Multi-Input Multi-Output (MIMO) system dynamically coupled of the ADN, we propose a design method articulated in two sequential steps. In the first one, each PI-IMC voltage controller is designed for the corresponding Effective Open-Loop Transfer Function (EOTF) [
28], where the free parameters of all controllers are designed to mitigate the interactions among the voltage control loops in a frequency range of interest (external interaction). In the second step, due to the reduced external interaction, the free parameters of all PI-IMC active-power controllers and the cut off frequency of all filters are designed to mitigate the interaction between voltage and active power control loop of any DER (internal interaction). The robust stability of the IMC-PI controllers is also investigated in presence of model parameter uncertainty modeling unknown scenarios of the ADN.
Summarizing, the contributions of this paper can be synthesized into: i. the AS of voltage regulation provided by DERs locally varying, primarily, their reactive power injection according to an integral law, and, in case of reactive current saturation, their active power injection according to a given droop law; ii. decentralized approach without any need of communication among DERs; iii. a single parameter to design each IMC-PI controller; iv. reduction of the internal and external interaction level; v. robust stability without any need of communication among local controllers; v. simplicity of the controller’s structure easily implementable in real distribution systems.
The paper is organized as follows. In Section II we firstly recall the adopted model of the ADN; subsequently we illustrate the proposed technique and validate the robust stability. Finally Section III is devoted to the development of a detailed numerical case studies aimed at validating the performance of our technique also in comparison with the technique proposed in [
29].
2. The proposed control design
The scope of this section is to illustrate the design methodology for the decentralized voltage control with reduced loop interactions. The control scheme is shown in
Figure 1.
Any DER is equipped with two control loops; the former controls the active power
while the latter the voltage
(
). The
i-th decentralized control matrix is
. As it can be observed from
Figure 1, the two loops are coupled by the filtered voltage error
which realizes the AS of voltage regulation by active power. In practice, in absence of saturation of
, the integral action of
imposes
; hence the contribution of active power to voltage regulation is null at steady-state. Conversely, if
saturates,
, the signal
acts with the scope of varying the active power
in order to reduce the steady-state voltage error.
The transfer function
is given by:
where the value of the droop constant
is chosen according to the AS that is agreed to be provided by the active power of each DER according to economical/commercial criteria. The mission of the first-order low pass filter
is to reduce the effect of the filtered voltage error on the active power during the transient.
In the following, after briefly recalling the plant model of the ADN with DERs, the design of and is presented, starting from the nominal plant matrix obtained assuming that the ADN operates in a known operating condition.
2.1. ADN Model
Concerning the ADN, it is composed of the local transfer functions
which represent the generic
and the matrix
that represents the distribution grid. The model is briefly recalled in the following; refer to [
30] for its foundation.
Concerning the DERs, they are assumed to be composed of a renewable energy generator and a battery energy storage system, connected on a DC bus, which is interfaced to the distribution grid by an inverter, an AC filter and a transformer. According to [
31], the following TITO model is adopted to represent the closed-loop transfer function of the current control loops of the inverter, acting on the
components
where the transfer functions matrix
is expressed by
Concerning the distribution grid, it is modelled by the full algebraic matrix , which linearly represents the effect of the vector of the currents injected by the DERs on the output vector of the injected active powers and of the DERs nodal voltages.
The resulting model of the ADN, including both the DERs and the distribution grid, is
where
is the MIMO model with
inputs
, and outputs
y and
. Moreover
, see (
3), where
is the matrix that represents the steady-state (
) model.
Before illustrating the design technique, it is necessary to discuss about the mutual interaction among the control loops of DERs. For any DER there are two types of interaction. The former is the internal interaction, represented by the TITO model (
2)-(
3); the latter is an external interaction caused by other DERs through the distribution grid, represented by the matrix
.
One of the most popular tools used to quantify the degree of interaction is the Relative Gain Array (RGA) defined as [
32]
where the operator ⊗ denotes element-by-element multiplication or Schur product and
is the matrix that represents the steady-state (
) model.
From the analysis of the RGA matrix, paper [
30] has shown that the external interaction exists only between the voltage control loops. That is the regulation of
is affected by the regulation of all other
voltages and vice-versa. This means that the voltage control loops are fully coupled. Conversely, the interaction between the active power control loops is weak. Hence the regulation of
does not affect the regulation of all other
active powers and vice-versa. As concerns the internal interaction, the regulation of
does not influence the regulation of
; but, in the case of this paper, the converse is not true. In fact the regulation of
affects that of
by means of the filtered voltage error
.
Summarizing, external interaction exists only between the voltage control loops, whereas it has internal interaction from the voltage control loop to the active power one. Based on this result, we present a design procedure formed by two sequential steps. In both steps the Internal Model Control (IMC) technique is adopted since exhibits robustness to model parameter uncertainty. In particular, in the first step, control matrix is designed for achieving voltage regulation and reduction of the external interactions. Conversely, in the subsequently second step, control matrix is designed for achieving active power regulation and reduction of the internal interaction.
2.2. First Step of the Design
Let now address the problem of designing
. Since the active power control loops are decoupled, it is possible to consider only the transfer functions corresponding to the pairings
, see
Figure 1. To this scope the matrix plant
is written by evidencing the
inputs and outputs as follows
The square matrix plant with inputs
and outputs
used in this step is then given by.
Employing the method of the Effective Open-Loop Transfer Function (EOTF), , the design of the IMC controller is developed for model rather than , thus avoiding a more complex MIMO design.
The expression of
is given by [
28]
where
and
denote, respectively, the
i-th diagonal element of
and of the Dynamic Relative Gain Array (
DRGA) calculated by
The IMC voltage control matrix is
where
are positive integers chosen to make
realizable.
From
it is obtains the feedback control matrix
The adjustable control parameters
are selected with the objective of reducing the external interactions. To this scope we formulate the controller design problem into an optimization framework. In details,
are obtained by solving the following problem
subject to
where
is a a frequency range of interest and
is the interaction matrix obtained from the following closed voltage control matrix
by setting
. In this way the response
due to a change of
(
) is effectively attenuated.
Remark 1.
As above described, the design of uses model in (6). The EOTF relates to when the i-th loop is open while all other loops are closed under the condition of perfect control. To explain this concept, let consider the actual relationship between to given by [28]
where , denotes a transfer functions matrix where both the i-th row and column are removed from , and are the i-th row and column vector of matrix where is discarded, respectively.
If the following approximation
holds in the frequencies range smaller than the cross-over frequency , then eq. (9) can be reasonably simplified as follows
Condition (10) guarantees the perfect control and it is fulfilled if , (, ) has a pole in and its gain is greater than that of , respectively.
2.3. Second Step of the Design
In the second step of the procedure, any controller
is independently designed from all the others since the active control loos are decoupled. The expression of the
i-th IMC active power controller is
with
a positive integer chosen to make
realizable. The corresponding feedback controller is
Since the external interactions have been minimized in the first step of the design procedure, the relationship between
and
(internal interaction), can be modeled by the element
of the following TITO closed transfer matrix
with
The internal interaction can be reduced by finding parameters
and
which solve the following problem
subject to
with
(
) positive quantities.
The described design technique assumes that
is a constant matrix; however, its value depends on factors such as the power injected through the DERs, the power consumed by the loads at each node, the topology of the distributions system, etc... Hence in presence of uncertainty in the elements of matrix
it is necessary to investigate the robust stability of the control system. To this aim, we have determined a matrix plant
for each of the considered operating conditions. In particular, the following ones have been considered: the reactive power ranges in the interval
pu; the voltage at the slack bus varies in the interval
pu; loads in the range
. At the end of the procedure, a set formed by 180 different matrices
has been created. Subsequently, for any
the corresponding MIMO closed-loop matrix
is built and the largest pole computed. The obtained set formed by the 180 poles is shown in
Figure 2. One can see that the requirement of robust stability is fulfilled.