Next Article in Journal
Short-Term Load Forecasting Algorithm Using a Similar Day Selection Method Based on Reinforcement Learning
Next Article in Special Issue
Smart Control Strategies for Primary Frequency Regulation through Electric Vehicles: A Battery Degradation Perspective
Previous Article in Journal
Wave Power Assessment in the Middle Part of the Southern Coast of Java Island
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Model Predictive Control Strategy for Distribution Grids: Voltage and Frequency Regulation for Islanded Mode Operation

DIBRIS—Department of Informatics, Bioengineering, Robotics and Systems Engineering, University of Genoa, 16145 Genoa, Italy
*
Author to whom correspondence should be addressed.
Energies 2020, 13(10), 2637; https://doi.org/10.3390/en13102637
Submission received: 23 March 2020 / Revised: 12 May 2020 / Accepted: 19 May 2020 / Published: 21 May 2020
(This article belongs to the Special Issue Optimal Control of Smart Distributed Power and Energy Systems)

Abstract

:
In the last few years, one of the most important challenges of power technologies has been the integration of traditional energy production systems and distributed energy resources. Large-scale photovoltaic systems and wind farms may decrease the quality of the electrical grid service, mainly due to voltage and frequency peaks and fluctuations. Besides, new functionalities, such as the operation in islanded mode of some portions of the medium-voltage grid, are more and more required. In this respect, a model predictive control for voltage and frequency regulation in interconnected local distribution systems is presented. In the proposed model, each local system represents a collection of intelligent buildings and microgrids with a large capacity in active and reactive power regulation. The related model formalization includes a linear approximation of the power flow equations, based on stochastic variables related to the electrical load and to the production from renewable sources. A model predictive control problem is formalized, and a closed-loop linear control law has been obtained. In the results section, the proposed approach has been tested on the Institute of Electrical and Electronics Engineers(IEEE) 5 bus system, considering multiple loads and renewable sources variations on each local system.

Graphical Abstract

1. Introduction

The increase of intermittent renewable energy sources (RES) has created instability issues, requiring new controllers for modern smart microgrids [1,2,3]. Actually, RES may cause voltage/frequency peaks and fluctuations that negatively affect the quality of the electrical grid service. On the other hand, RES represents an important alternative to reduce global fossil fuel consumption. The need for new control strategies is also motivated by the increase of distributed small generators that cause millions of new producers potentially involved in the energy market. In such a scenario, in the balancing market, new actors, such as aggregators [4], are present and new functionalities are more and more required. For example, in order to preserve the electrical grid from faults and emergency situations, it may be necessary to reduce the power demand and to operate some portions of the medium-voltage grid in islanded mode. These new set-ups require new fast controllers, which can take into account voltage and frequency models.
In addition, the grid operators, i.e., the transmission system operator (TSO) and distribution system operator (DSO), are forced to deeply change their roles. As an example, the DSO must become active in its activities as a grid manager [5]. The DSO, in particular, is responsible for the normal operation of the distribution grid and its subsystems (i.e., managing power quality and energy losses). Among others, the new roles for the DSO are: operational security management, dynamic voltage and frequency control, outages response, performing restoration actions, and market operations. In this framework, it becomes hard to define an effective optimal control, which derives from the structure of the power grid, and, specifically, from the presence of several issues: renewable and traditional power production, bidirectional power flows, and stochastic modelling issues (due to uncertainties in forecasting power generated from renewable resources). The main challenge is to define and to solve control problems embedding prediction models for the regulation of frequency, voltage, and active and reactive power flows. Indeed, the complete power flow equations are nonlinear and nonconvex [6]. Therefore, optimization problems can hardly apply them to compute the control in real time. On the other hand, there is the need to represent in detail the physical system to provide reliable control strategies [7].
This work focuses on the optimal control of interconnected local systems in a distribution grid. A local system can represent a group of renewable generators, passive and active loads, or microgrids with a large capacity for active and reactive power regulation. The proposed approach can be used:
  • To reduce frequency oscillations and to improve inertia response due to a sudden change in the grid (i.e., a variation on a renewable power plant production), and
  • To manage a portion of the distribution grid that is operated in islanded mode because of demand response programs and/or emergency requirements.
In fact, an increase in the RES penetration can decrease the number of generation units that provide reserve power for primary and secondary control. So, the frequency deviation increases [8]. New equations for inverters’ control to regulate frequency and voltage are required when a change in the grid occurs (e.g., a loss of load).
The modelling approach proposed in this paper is based on the concept of virtual inertia [9,10,11]. An electrical model is defined, including voltage (magnitude and phase), frequency, active and reactive power, and stochastic components. So, the definition of an optimal control problem for frequency and voltage regulation in multi-local systems is introduced. The aim is to define a closed-loop controller for the mitigation of frequency and voltage variations.
The paper is organized as follows. Section 2 reports the state-of-the-art and discusses the innovation of the proposed work. In Section 3 and Section 4, the system and the optimal control problem are formalized. Finally, the application to a case study is described in Section 5, while conclusions and future developments are reported in Section 6. The Appendix A reports proofs related to demonstrations not included in the text for sake of readability.

2. Literature Review

Nowadays, the control of interconnected local distribution systems is a challenging problem. In this respect, several papers focus on the definition and solution of optimization and control problems. A detailed survey [12] on multi-microgrids systems has been recently published, where operational management and control are quoted as challenging problems to be faced. Zheng et al. [13] propose a distributed model predictive control (MPC) approach using discrete-time Laguerre functions for load frequency control in multi-area interconnected power systems, in order to overcome the problems of computational burden and online optimization. The advantage of an MPC approach is that it allows for implementing a multivariable controller that controls the outputs simultaneously by taking into account all the interactions among system variables. Another strength of MPC (not used in our work, where a closed fast feedback is adopted) is that it can handle constraints. Khooban and Niknam [14] propose a heuristic algorithm and a fuzzy logic approach to tune parameters of classical proportional integral (PI) controllers for multi-area load frequency control. However, electrical models are not detailed. In the literature of microgrids, there are several interesting approaches based on MPC but related to the single microgrid, and thus not interconnected with other local systems. In Reference [15], a novel fractional-order model predictive control method is proposed to achieve the optimal frequency control of an islanded microgrid by introducing a fractional order integral cost function into the MPC approach. As in the present work, frequency regulation is attained, even if hereinafter we propose a more detailed system model, where active and reactive powers are present and voltage is regulated.
In Reference [16], a feedback linearization is applied to design a predictive controller for the secondary control of the voltage and frequency in an islanded microgrid. The goal is to improve the performance of a closed-loop nonlinear system. The control scheme is fully distributed, and necessary and sufficient conditions for convergence and stability of the whole system are described. As a similar approach, we model both active and reactive power, but providing a control strategy for the whole optimization problem that can be used under a MPC scheme. In Reference [17], a MPC-based controller uses a simplified voltage prediction model to predict the voltage behavior in an islanded microgrid. Apart from the fact that we consider multiple microgrids, in our model, we use a more detailed electrical model and we derive an optimal control law. The same considerations can be done in comparison with References [18,19].
The approach considered in the present paper gives a strong emphasis on the model of the electrical grid. Specifically, the approach is based on the concept of virtual inertia [20,21]. The problem of virtual inertia control is mainly felt in microgrids and grids where RES are present, as it is necessary to cope with small inertia and uncertainties. In Reference [22], the authors propose an approach to increase the inertia of the photovoltaic (PV) system through inertia emulation, which can be obtained by controlling the charging/discharging of the DC-link capacitor and adjusting the PV generation. Wu et al. [23] consider an islanded microgrid with renewables and storage systems, in which virtual inertia allows power curtailment of RES units and demand response (loads’ shifting). In fact, these approaches are also useful for islanded microgrids, and several works are present in the recent literature for frequency and voltage regulation, and for stability [24,25,26,27,28].
In Reference [24], the authors assess that droop control methods are not suitable for renewable energy-based microgrids due to their limited capability for power delivery. They propose a supercapacitor-based DC-link voltage stabilizer and a battery energy storage system-based frequency stabilizer. In particular, as in Reference [25], a fuzzy-based controller is proposed. In Reference [26], the authors highlight the drawbacks of PI controller gains and propose the Grasshopper Optimization Algorithm to optimize the PI controller parameters.
An additional critical issue is represented by the demand of electric vehicles in a smart grid. To face such problem, in Reference [27], distributed controllers are proposed. Other papers [28] focus on improvements of droop controllers and on new approaches for the optimal dispatch in hierarchical microgrids [25]. Specifically, two control modes are considered for the generation units: a voltage control mode, with an active droop control loop, and a power control mode, which allows setting the output power in advance. Some recent works that study virtual inertia are References [29,30,31,32,33]. In Reference [29], the authors study a novel control scheme for virtual inertia converters to operate in unbalanced settings, operating test on real hardware. In Reference [30], an adaptive controller method for a PV-based virtual inertia system is discussed. The approach is studied in the case of the small-signal stability and low-frequency oscillation, which are characteristics of the PV generators. In Reference [31], the authors aim at improving the inertia of a microgrid by reducing the time delay of the communication system to overcome frequency oscillations, finding a compromise with the increase of operational costs.
In Reference [32], the authors propose an approach to find the proper value of inertia constant, together with the frequency droop coefficient of DERs (distributed energy resources) and tuning load frequency controllers’ parameters to improve the frequency stability. A multi-objective decision problem is formulated in which the following goals are considered: minimizing the maximum frequency deviation while promptly bringing the frequency back to the desirable region, and minimizing a higher inertia constant that results in a higher cost of energy storages. With a significant difference, here, we propose an optimal control problem to minimize frequency and voltage deviations for multiple microgrids and the optimization problem is analytically solved taking into account both active and reactive power.
In Reference [33], virtual inertia in islanded microgrids is studied. An H-infinite (H∞) robust control approach to the virtual inertial control loop is implemented, taking into account the high penetration of renewables, thus enhancing the robust performance and the stability of the microgrid during contingencies. In the present paper, we propose a control law to minimize frequency and voltage deviations taking into account both active and reactive power under an MPC scheme.
Hereinafter, the electrical grid connecting different local systems is modelled. This paper aims at studying the regulation of voltage and frequency within local systems interconnected through a distribution network. Since these systems have a strong renewable, uncontrolled, and intermittent component, they intrinsically represent a disturbance that, if not controlled, may cause a decrease in the system performance and quality. It is worthwhile to observe that the problem faced is not to be confused with the classic problem of inter-area oscillations between systems connected by a high-voltage grid.
In this work, an MPC approach is proposed, taking into account uncertainties, and deriving an optimal solution that can be used in closed-loop and in real time. In the literature of power systems, similar control approaches can be found [34,35]. Delfino et al. [36] propose a multilevel architecture for interconnected microgrids in which the DC power flow equations are used at a higher level, an analytical solution emulates each microgrid, and the AC power flow equations are used at the lower level. Minciardi and Robba [37] formalize a decision model in which a stochastic optimal control problem is defined and solved. An optimal control strategy can be also found in Reference [38]. The authors consider the problem of load-frequency control in a multi-area power system proposing a hierarchical optimal robust controller. They decompose the overall problem in subsystems, solve the lower level through a control strategy, and finally, iteratively coordinate local systems. However, a detailed representation of the electrical network is not present. Annaswamy and Kiani, in Reference [39], include a discrete-time decentralized linear quadratic regulator (LQR) controller in a transactive architecture for hierarchical smart grids. In the present paper, with respect to References [34,35], a detailed representation of the electrical network is included in order to regulate voltage, and control active and reactive power. With respect to References [36,37], the decision problem has been solved analytically. In addition, References [34,35,36,37,38,39] do not consider frequency variations and equivalent inertia. Other optimal control strategies for the smart grid are present in the literature [40,41] but not in the area of interconnected local systems.
In summary, the main contributions of the proposed paper are:
  • The statement of the frequency and voltage control problem of distribution grids in islanded mode expressed as an MPC of an affine system affected by additive noise.
  • The inclusion of models for the control of frequency, voltage, and active and reactive power in distribution grids.
  • The derivation of an analytical solution for the stochastic MPC problem based on a two-stage decomposition.
  • An optimal control strategy for the inter-local systems variation problem that can be used either for frequency and voltage regulation or the operation in the islanded mode of portions of the distribution grid for demand response purposes.

3. The System Model and the Optimization Problem

Consider a distribution power grid (Figure 1) modeled as an undirected graph G = (A, L), where the set of nodes A = { 1 , , N } characterizes the set of local systems and L = { 1 , , b } is the set of links that connect the nodes. Specifically, each local system is interfaced with the distribution grid by a power converter which can perform a virtual inertia control (allowing the inverter to behave as a fictitious synchronous generator) with equivalent inertia time constant H i ( i A ) and damping factor K D i ( i A ). The set of arcs, L, represents the set of lines that connect the local systems—each line has a proper impedance Z i , j = R i , j + j X i , j     ( i , j A ) , with R i , j ( i , j A ) being the line resistance and X i , j ( i , j A ) being the line inductance. All local systems are composed of different devices: distributed controllable resources with high regulation capacity (denoted by G), reactive power-regulating sources (denoted by NG), renewable non-controllable resources (denoted by RES), and passive loads (denoted by L). Each local system can be modeled as an equivalent synchronous generator with equivalent inertia to consider its contribution to the frequency regulation of the whole grid.
In the following paragraphs, the state and the control variables are defined, and the optimization problem is formalized. All the quantities (but the phases) and the equations are expressed at discrete time per unit (p.u.) values. As suggested by the literature, this method offers computational simplicity by eliminating units and expressing system quantities as dimensionless ratios.
The decision variables are:
  • P i , t i A , t = 0 , , T 1 : active power injected at node i.
  • Q i , t i A , t = 0 , , T 1 : reactive power injected at node i.
  • P G i , t i A , t = 0 , , T 1 : distributed controllable source active power at node i.
  • Q N G i , t i A , t = 0 , , T 1 : reactive power regulating source reactive power at node i.
  • P m i , t : regulation active power injected by local system i.
  • V i , t i A , t = 0 , , T 1 : voltage amplitude at node i
The parameters are:
  • P L i , t i A , t = 0 , , T 1 : load active power at node i.
  • Q L i , t i A , t = 0 , , T 1 : load reactive power at node i.
  • Q G i , t i A , t = 0 , , T 1 : distributed controllable source reactive power at node i.
  • P N G i , t i A , t = 0 , , T 1 : reactive power regulating source active power at node i.
  • P R E S i , t i A , t = 0 , , T 1 : renewable source active power at node i.
  • Q R E S i , t i A , t = 0 , , T 1 : renewable source reactive power at node i.
The state variables are:
  • δ i , t i A , t = 0 , , T 1 : voltage phase at node i (rad).
  • ω i , t i A , t = 0 , , T 1 : frequency at node i.
The decision variables detail for a local system is reported in Figure 2.
The electrical grid inside each local system i is not considered, and thus only the power balance must be satisfied:
P G i , t + P N G i , t + P R E S i , t P L i , t = P i , t i A , t = 0 , , T 1
Q G i , t + Q N G i , t + Q R E S i , t Q L i , t = Q i , t i A , t = 0 , , T 1
Each local system is associated with a time-varying voltage phasor of magnitude V i , t and phase δ i , t . The discrete-time system dynamics ( Δ t is the time interval length) are governed by the equations of the virtual inertia controller, linking the node angle, δ i , t , to the frequency difference to the nominal value ( ω i , t ω 0 ):
δ i , t + 1 = δ i , t + ( ω i , t ω 0 ) Δ t i A t = 0 , , T 1
where the frequency is related to the active power supply–demand imbalance:
ω i , t + 1 = ω i , t + P m i , t P G i , t K D i ( ω i , t ω 0 ) 2 H i Δ t i A , t = 0 , , T 1
where K D i is the local system damping factor. The distribution grid model includes the linearized power flow equations for a medium voltage AC grid, with G i , j = 1 / R i , j ( i , j A ) as the line conductance, and B i , j = 1 / X i , j ( i , j A ) as the line susceptance. That is:
P i , t = j A G i , j V j , t + B i , j ( δ i , t δ j , t ) i , j A t = 0 , , T 1
Q i , t = j A B i , j V j , t + G i , j ( δ i , t δ j , t ) i , j A t = 0 , , T 1
The model described by Equations (5) and (6) is the well-known linearized version of power flow equations, obtained considering the first term of Taylor expansion of the sinusoidal terms of the equations and just considering the voltage difference between nodes [42,43,44].
In the following paragraphs, the linear model described by the previous equations is reformulated in the state space canonical form.
Proposition I: The system model, described by Equations (1)–(6), can be expressed as an affine system affected by additive noise:
x t + 1 = A x t + B u t + C h t + C ^ h ^ t t = 0 , , T 1
with
  • u t = [ P m i , t , i = 1 , , N 1 , Q T N i , t , i = 1 , , N 1 ] ,
  • x t = [ δ i , t , i = 1 , , N 1 ; ω i , t , i = 1 , , N 1 ] ,
where A, B, C, and C ^ are known matrices, h t is a vector of parameters, and h ^ t represents stochastic noise with E { h ^ t } = 0 .
Proof of Proposition I: The Proof is shown in the Appendix A.1.
In the day ahead, the DSO plans the production and energy flow based on different objectives, such as the minimization of losses and the capacity of the lines. It is reasonable to consider variations concerning an already predetermined working point. However, given the strong variability of part of the energy sources involved in our problem, control policies are required to mitigate deviations from an initial working point.
The proposed objective function aims at minimizing four different contributions: the square of the difference between the phases of the grid nodes, the quadratic deviation of the frequency from the steady-state value, the voltage quadratic deviation from the desired value, V 0 , and the quadratic deviation of the active regulation power to the set point P m :
min   [ ( t = 0 T i , j A i j ( δ i , t δ j , t ) 2 ) + ( t = 0 T i A ( ω i , t ω 0 ) 2 ) + ( t = 0 T 1 i A ( V i , t V 0 ) 2 ) + ( t = 0 T 1 i A ( P m i , t P m i ) 2 ) ]
Proposition II: The objective function (8) can be expressed as a stochastic quadratic function with mixed terms:
min J ( x t , u t ) = 1 2 E { t = 0 T 1 [ ( x t x ^ ) Q t ( x t x ^ ) + ( u t u ^ t ) R t ( u t u ^ t ) + 2 ( x t x ^ ) F t ( u t u ^ t ) ] + ( x T x ^ ) Q T ( x T x ^ ) }
where x t and u t are the state and control vectors respectively, x ^ , u ^ are the tracking parameters, and Qt, Rt, Ft , QT are known matrices.
Proof of Proposition II: The Proof is shown in the Appendix A.2.

4. The Model Predictive Control Problem

4.1. The MPC Scheme

The MPC scheme is defined in Figure 3. The optimization horizon (T) is the length over which the optimization problem runs and (S) is the simulation horizon. Both horizons (T) and (S) can be chosen according to the available information about uncertain parameters/forecasts (current states, forecasted demands, and renewables). In the MPC approach, in the first run, the solution related to the first time interval is considered. Then, uncertain parameters (new information from sensors in the field and new forecasts) are inserted in the procedure and a new optimization problem is run (second run), and so on for the whole duration of the simulation horizon (S). The MPC approach allows for updating new information when available and reducing uncertainties related to renewables’ availability and load forecasting.

4.2. Explicit Solution of the MPC Problem

The optimization problem is in the form:
min J ( x t , u t ) = 1 2 E { t = 0 T 1 [ ( x t x ^ ) Q t ( x t x ^ ) + ( u t u ^ t ) R t ( u t u ^ t ) + 2 ( x t x ^ ) F t ( u t u ^ t ) ] + ( x T x ^ ) Q T ( x T x ^ ) }
subject to (s.t.) Equation (7).
By a change of variables z t = x t x ^ and e t = u t u ^ t , the optimization problem is given by:
min J ( z , e ) = 1 2 E { t = 0 T 1 [ z t Q t z t + e t R t e t + 2 z t F t e t ] + z T Q T z T }
s.t.
z t + 1 = A z t + B e t + μ t + C ^ h ^ t t = 0 , , T 1
where:
μ t = C h t + ( A I ) x ^ + B u ^ t t = 0 , , T 1
The Equations (11)–(13) are a stochastic MPC problem, which is “non-standard,” due to the presence of the known input sequence, μ t , in the state equation. In the following claims, it is shown how to solve the MPC problem in a closed-form through a two-stage decomposition.
Claim I: Due to its linearity, the decision problem in Equations (11)–(13) can be decomposed separating the stochastic ( z t + 1 s ) and deterministic ( z t + 1 d ) components [45]:
J ( z , e ) = J ( z s , e s ) + J ( z d , e d )
s.t.
z t + 1 s = A z t s + B e t s + C ^ h ^ t t = 0 , , T 1
z 0 s = 0
z t + 1 d = A z t d + B e t d + μ t t = 0 , , T 1
z 0 d = z 0
The optimal control strategy is given by:
z t = z t s , + z t d , t = 0 , , T 1
e t = e t s , + e t d , t = 0 , , T 1
Proof of Claim I: The Proof is shown in the Appendix A.3.
Claim II: The optimal control problem to determine e t d , , due to its linearity, can be written as:
min J ( z d , e d ) = 1 2 E { t = 0 T 1 ( z t d 1 r t ) Q t ( z t d 1 r t ) + e t d R t e t d + 2 ( z t d 1 r t ) F t e t d + ( z T d 1 r T ) Q T ( z T d 1 r T ) }
where:
z t d = z t d 1 + z t d 2 t = 0 , , T 1
z t + 1 d 1 = A z t d 1 + B e t d t = 0 , , T 1
z 0 d 1 = 0
z t + 1 d 2 = A z t d 2 + μ t t = 0 , , T 1
z 0 d 2 = z 0
z t d 2 = r t t = 0 , , T 1
Proof of Claim II: The Proof is shown in the Appendix A.4.
Proposition III: The optimal control law for the problem in Equations (11)–(13), and taking into account Claim I and Claim II, is given by:
e t = K t ( z t s + z t d 1 ) + K t g ( B g t + 1 + F t r t )   t = 0 , , T 1
with
g t = [ ( A M t ) P t + 1 ( I + E t P t + 1 ) 1 ( E t g t + 1 + M t r t ) ( A M t ) g t + 1 ] + ( Q t D t ) r t t = 0 , , T 1
K t = ( R + B P t + 1 B ) 1 ( B P t + 1 A + F t )   t = 0 , , T 1
P t = Q t D t + ( A M t ) P t + 1 ( I + E t P t + 1 ) 1 ( A t M t ) t = 0 , , T 1
K t g = ( R t + B P t + 1 B ) 1   t = 0 , , T 1
r t + 1 = A r t μ t
P T = Q T   t = 0 , , T 1
g T = Q T r T   t = 0 , , T 1
being: M t = B R t 1 F t , D t = F t R t 1 F t , D t = F t R t 1 F t and E t = B R t 1 B .
Proof of Proposition III: The Proof is shown in the Appendix A.5.
As reported in Claim I and Claim II, the original problem can be so decomposed into a stochastic problem and a deterministic discrete tracking problem. The closed-loop optimal control strategy is obtained as follows [46]:
Sub-problem 1: Stochastic problem.
The optimal control problem to determine e t s , is given by:
min J ( z , e ) = 1 2 E { t = 0 T 1 [ z t s Q t z t s + e t s R t e t s + 2 z t s F t e t s ] + z T s Q T z T s }
s.t.
z t + 1 s = A z t s + B e t s + C ^ h ^ t t = 0 , , T 1
z 0 s = 0 t = 0 , , T 1
Proposition IV: The solution of the optimal control problem, Equations (36)–(38), is given by the following recursive equations:
e t s = K t z t s   t = 0 , , T 1
with
K t = ( R + B P t + 1 B ) 1 ( B P t + 1 A + F t ) t = 0 , , T 1
P t = Q t D t + ( A M t ) P t + 1 ( I + E t P t + 1 ) 1 ( A t M t )   t = 0 , , T 1
P T = Q T
being M t = B R t 1 F t , D t = F t R t 1 F t and E t = B R t 1 B .
Proof of Proposition IV: The Proof is shown in the Appendix A.6.
Sub-problem 2: Deterministic problem.
As reported in Claim II, the optimal control problem to determine e t d , can be written as:
min J ( z d , e d ) = 1 2 E { t = 0 T 1 [ ( z t d 1 r t ) Q t ( z t d 1 r t ) + e t d R t e t d + 2 ( z t d 1 r t ) F t e t d ] + ( z T d 1 r T ) Q T ( z T d 1 r T ) }
s.t.
z t + 1 d 1 = A z t d 1 + B e t d t = 0 , , T 1
z 0 d 1 = 0
being z t d 2 = r t , z 0 d 2 = z 0 , z t + 1 d 2 = A z t d 2 + μ t , μ t = C h t + ( A I ) x ^ + B u ^ t .
Proposition V: The optimal control strategy for the problem, Equations (43)–(45), is given by
e k d , = K t z t d 1 + K t g ( B g t + 1 + F t r t ) t = 0 , , T 1
with
g t = [ ( A M t ) P t + 1 ( I + E t P t + 1 ) 1 ( E t g t + 1 + M t r t ) ( A M t ) g t + 1 ] + ( Q t D t ) r t t = 0 , , T 1
K t = ( R t + B P t + 1 B ) 1 ( B P t + 1 A + F t ) t = 0 , , T 1
K t g = ( R t + B P t + 1 B ) 1 t = 0 , , T 1
P t = Q t D t + ( A M t ) P t + 1 ( I + E t P t + 1 ) 1 ( A M t ) t = 0 , , T 1
r t + 1 = A r t μ t t = 0 , , T 1
P T = Q T
g T = Q T r T

5. Application to a Case Study

5.1. Case Study Description

To show the effectiveness of the theoretical results presented in the previous section, the proposed control strategy has been applied to a modified IEEE 5-bus system (N = 5) [47], shown in Figure 4. The test network has been adapted to our problem through the following assumptions. Each bus corresponds to a local system, and, keeping the network topology unchanged, the values of cables’ resistance, Ri,j, and reactance, Xi,j, have conveniently rescaled to a distribution grid system (assuming base quantities of 10 MVA for apparent power, V 0 = 15 KV for voltage and ω 0 = 50 Hz for frequency). Table 1 shows the detailed values of the line parameters Ri,j, and Xi,j for the case study.
The parameters that characterize each local system, i.e., the inertia constant, Hi, and the equivalent damping factor, KDi, are reported in Table 2. The local system 5 is the slack node, in which the voltage, V 5 , t , is fixed to 1 (p.u.) and the phase, δ 5 , t , is equal to 0 (rad). The case study can be representative of small local systems connected by a distribution grid (such as interconnected microgrids), in which there is a great presence of generation from renewable resources.
The simulation horizon is 60 s, with a discretization step Δ t of 0.5 s. This horizon length has been chosen as a trade-off between the following two aspects: (a) to avoid delays due to the centralized communication from the controller to the local systems, and (b) to enhance the reliability of the short-term forecasting of RES and loads.
To make the approach more effective, the MPC strategy is applied, in which only a part of the optimal control law is applied (e.g., 10 seconds over a T = 60 second horizon), and there is a continuous recalculation of the control law based on new measurements and forecasting. The length of the MPC control horizon is justified by recent developments on the nowcasting of renewable power sources that predict 1 min future PV values [48,49].
At this point, by the means of the local systems’ parameters and the distribution network, it is possible to quantify the matrices A and B, defined in the Appendix A.1, described by Equations (A19) and (A20), and needed by Equation (7):
A = ( 1 0 0 0 0.5 0 0 0 0 1 0 0 0 0.5 0 0 0 0 1 0 0 0 0.5 0 0 0 0 1 0 0 0 0.5 0.22 0.21 3.8 e 19 0.035 0.76 0 0 0 0.16 0.47 0.062 0.046 0 0.66 0 0 2.08 e 18 0.07 0.54 0.070 0 0 0.58 0 0.0421 0.07 0.093 0.20 0 0 0 0.80 )
B = ( 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.05 0 0 0 0.33 0 6.93 e 18 2.77 e 17 0 0.062 0 0 0 0.33 4.51 e 17 6.93 e 17 0 0 0.055 0 4.6 e 17 5.55 e 17 0.33 2.77 e 17 0 0 0 0.041 5.55 e 17 2.77 e 17 1.38 e 17 0.33 )
In each MPC step, the optimal control law is the sum of the solution of two different control problems: a stochastic problem defined by Equations (36)–(38), with the optimal solution, Equations (39)–(42), and a deterministic problem, Equations (43)–(45), with the optimal solution, Equations (46)–(53). Both control laws have been recursively defined, i.e., all the parameters that define the control can be calculated offline, the parameters are known, and the forecasts of renewable generation and load.
The variations on the system are caused by the combined variation of load and RES Psc,i,t on the local systems defined by Equation (A4), represented into the model in the terms h t and h ^ t defined by Equations (A21) and (A22) respectively, deterministic and stochastic parts of the noise. Figure 5 reports the values, P s c i , t d , for the whole simulation horizon. Data have been obtained by adapting real power measurements from photovoltaic, residential, and industrial loads, and rescaled to obtain congruent measurements for the distribution networks. The stochastic part of the noise, P s c i , t s ~ N(0,1), affects the state at each instant t.
The power profiles on each local system could be critical for the whole system. Without any upper-level controller but only with the damping given by each local system, the generators could not be able to efficiently compensate (i.e., in a reasonable time, according to the dynamics of the system) the given power variation, causing frequency and voltage variations.
The objective of the proposed controller is to reduce the voltage and frequency variations (during the transient period) and, after the perturbation, bringing them back to their reference value. In a practical context, numerous frequency variations and voltage peaks on the grid could cause the intervention of protection systems, sensitive to variations in voltage or frequency, disconnecting some portions of the network itself, causing a blackout for many users.

5.2. Optimal Control Results

Optimal solutions for the state variables in the case with control or without control are reported in Figure 6 and Figure 7. It can be seen how the presence of the proposed control can limit frequency variations (i.e., the power quality remains within acceptable limits, the voltage deviation from the reference value is less than 5%). The dotted line represents the system with no control (only the damping coefficient is considered). In this case, each local system operates without any central controller (like the one developed here), but the system only reacts based on the physics of the system with constant inputs equal to those of the equilibrium point. The bold line represents the controlled system using the developed approach. It is important to note that without the proposed control strategy, the system does not reach a steady state for the whole simulation.
Figure 8 shows the voltage pattern for each local system in the presence or absence of the developed control. The control strategy helps in containing the variations affecting this variable, improving the system’s power quality. Finally, Figure 9 and Figure 10 show the optimal active and reactive power at each local system.
When there is no control, the objective function value is equal to 0.1098, while when the proposed control strategy is applied, it is 0.0151.
Table 3 shows the objective function values and contributions.
The test has been developed in MATLAB-SIMULINK (2018b, Mathworks, Natick, Massachusetts, USA), running on a PC Intel i7, 16 GB RAM. The MATLAB tool has been used to calculate a priori the control laws defined recursively in Section 4.

5.3. Closed-Loop Stability Analysis

This subsection aims to perform a sensitivity analysis concerning the asymptotic stability properties of the proposed system, with respect to the damping parameter, KDi. The choice of this parameter is fundamental with regards to the design of the control system based on virtual inertia. It is well known from power systems theory that the asymptotic stability of a synchronous machine is directly linked to its damping factor, similarly for the virtual inertia that mimics this behavior.
The first part of Figure 11 shows how different choices of damping factors at each node affect the norm of the eigenvalues of the open-loop system for each element of the state vector x . The values are chosen in an interval between ± 50 % of the ones used for the use case, 20 different damping factors are considered.
It is noteworthy that for lower values of KDi, some eigenvalues have a norm greater than 1, so the system is unstable. This is consistent with the power system theory. To this end, the asymptotical closed-loop stability of the proposed controller for an infinite time control law defined is analyzed, as in Reference . As can be seen from the second part of Figure 11, the proposed optimal controller also guarantees the asymptotical stability in the cases in which the system was designed as unstable. This aspect is crucial for the optimal control of this kind of problem, as it also fixes design instabilities.

6. Conclusions and Future Developments

The proposed work introduced an approach based on optimal control for the regulation of voltage and frequency in multi-local systems (such as interconnected microgrids) in a distribution grid. This represents one of the most challenging problems regarding the dynamic aspects of a smart grid.
If not considered and properly mitigated, these aspects can drastically damage the entire system. To reduce these effects, the paper considers accurate grid models and the presence of innovative devices that allow for better regulating the power flows on the grid and, consequently, voltage and frequency at each node. It is worthwhile to observe that the main contribution of this paper is not to assess the controller performance (that is optimal by definition to the linearized model) but to propose a model that can be linearized and controlled by a closed-form control law.
To cope with the variability represented by renewable sources and load, this approach was used in an MPC framework. In this case, just the control law corresponding to the first time interval would be applied, while an update of the law would be performed according to the new forecasts made in the short term. This would allow a further increase in the quality of the system. The IEEE 5 bus system was used to test the control strategy, considering multiple loads and renewable sources variations on each local system. The results showed a strong decrease in the number of variations as well as strong damping of the frequency and voltage peaks in each local system. This greatly increases the power quality of the electrical system. A sensitivity analysis about the asymptotical stability was also given, showing that the proposed control can also prevent instabilities caused by an incorrect design of the damping parameters.
Future developments will concern several aspects. First, delays in the communication between the central controller and each local system could be inserted in the optimal control problem.
The particular model of the local systems will be the focus of future investigations. Furthermore, distributed and team robust approaches to overcome the poor robustness of totally centralized control approaches in case of emergencies will be considered.

Author Contributions

“Conceptualization, G.F., M.R. and R.S.; Methodology, G.F., M.R. and R.S.; Software, G.F., M.R. and R.S.; Validation, G.F., M.R. and R.S.; Writing-Review & Editing, G.F., M.R. and R.S.; Visualization, G.F., M.R. and R.S.; Supervision, G.F., M.R. and R.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Appendix A.1. Proof of Proposition I

Equations (5) and (6) in the absence of shunt admittances can be written in the compact form:
P ¯ t = G V ¯ t B p δ ¯ t + K p   t = 0 , , T 1
Q ¯ t = B p V ¯ t G δ ¯ t + K q   t = 0 , , T 1
with
V ¯ t = [ V i , t , i = 1 , , N 1 ] δ ¯ t = [ δ i , t , i = 1 , , N 1 ] Q ¯ t = [ Q i , t , i = 1 , , N 1 ] P ¯ t = [ P i , t , i = 1 , , N 1 ]
where K p and K q are known parameters coming from considering the voltage slack bus equal to 1. Equation (1) can be inserted in the frequency state Equation (4):
ω i , t + 1 = ω i , t + P m i , t + P N G i , t + P s c i , t P i , t K D i ( ω i , t 1 ) 2 H i Δ t i A t = 0 , , T 1
being
P s c i , t = P R E S i , t P L i , t = P s c i , t d + P s c i , t s i A t = 0 , , T 1
where P s c i , t d and P s c i , t s are the deterministic and stochastic parts of the net load, respectively.
The frequency state equation is given by:
ω ¯ t + 1 = A ω ω ¯ t + B ω P ¯ t + B ω , m P ¯ m , t + C ω h ¯ ω , t + C ω , s h ¯ ω , s , t t = 0 , , T 1
being:
ω ¯ t = [ ω i , t , i = 1 , , N 1 ]
P ¯ m , t = [ P m i , t , i = 1 , , N 1 ]
h ¯ ω , s , t = [ P s c i , t s , i = 1 , , N 1 ]
h ¯ ω , t = [ P N G i , t + K D i + P s c i , t d , i = 1 , , N 1 ]
where
A ω = [ 1 K D 1 2 H 1 Δ t 0 0 0 0 0 0 0 0 1 K D N 1 2 H N 1 Δ t ]
B ω = [ Δ t 2 H 1 0 0 0 0 0 0 0 0 Δ t 2 H N 1 ]
B ω , m = C ω = C ω , s = B ω
Substituting Equation (A1) in the frequency state equation, we obtain:
ω ¯ t + 1 = A ω ω ¯ t + B ω ( G V ¯ t B p δ ¯ t + K p ) + B ω , m P ¯ m , t + C ω h ¯ ω , t + C ω , s h ¯ ω , s , t t = 0 , , T 1
then, the voltage can be obtained from the Equation of reactive power (A1). That is:
V ¯ t = B p 1 Q ¯ t B p 1 G δ ¯ t + B p 1 K q t = 0 , , T 1
and substituting Equation (A7) in the frequency state Equation (A6) we have:
ω ¯ t + 1 = A ω ω ¯ t + B ω , m P ¯ m , t + C ω h ¯ ω , t + C ω , s h ¯ ω , s , t B ω G B p 1 Q ¯ t + ( B ω G B p 1 G + B ω B p 1 ) δ ¯ t + B ω G B p 1 K q + B ω K p t = 0 , , T 1
substituting Equation (2) (where the control variable is Q N G i , t ) into Equation (A8), we get:
ω ¯ t + 1 = A ω ω ¯ t + B ω , m P ¯ m , t + C ω h ¯ ω , t + C ω , s h ¯ ω , s , t B ω G B p 1 ( Q ¯ G , t + Q ¯ N G , t + Q ¯ R E S , t Q ¯ L , t ) ( B ω G B p 1 G + B ω B p ) δ ¯ t + B ω G B p 1 K q + B ω K p t = 0 , , T 1
with
Q ¯ G , t = [ Q G i , t , i = 1 , , N 1 ] Q ¯ N G , t = [ Q N G i , t , i = 1 , , N 1 ] Q ¯ R E S , t = [ Q R E S i , t , i = 1 , , N 1 ] Q ¯ L , t = [ Q L i , t , i = 1 , , N 1 ]
Let us pose:
A ω , δ = B ω G B p 1 G + B ω B p
h ¯ t = [ h ¯ ω , t ; Q G i , t + Q R E S i , t Q L i , t K q B p K p , i = 1 , , N 1 ] D ω = ( B ω G B p 1 )
C ω = C ω , s = d i a g ( Δ t )
C d = [ C ω D ω ] T
Q ¯ T N , t = B p 1 G 1 B p Q ¯ N G , t
then, we obtain:
ω ¯ t + 1 = A ω ω ¯ t + B ω , m P ¯ m , t + C d h ¯ t + C ω , s h ¯ ω , s , t Q ¯ T N , t + A ω , δ δ ¯ t t = 0 , , T 1
It is important to note that Q ¯ T N , t = [ Q T N i , t , i = 1 , , N 1 ] is an auxiliary variable that is crucial to define the objective function accurately (as it will be explained in the following). Thus, in the state equation, we have to satisfy the following condition: C T N Q ¯ T N , t = C T N B p 1 G 1 B p Q ¯ N G , t = B ω G B p 1 Q ¯ N G , t . Thus, one should write C T N = [ B p 1 G 1 B p ] 1 B ω G B p 1 .
Then, we obtain:
ω ¯ t + 1 = A ω ω ¯ t + B ω , m P ¯ m , t + C d h ¯ t + C ω , s h ¯ ω , s , t C T N Q ¯ T N , t + A ω , δ δ ¯ t t = 0 , , T 1
as regards the other state equation, one can write:
δ ¯ t + 1 = A δ δ ¯ t + A δ , ω ω ¯ t A δ , ω t = 0 , , T 1
where A δ ( N 1 x N 1 ) , A δ , ω ( N 1 x N 1 ) are known matrices.
Finally, posing x t = [ δ i , t , i = 1 , , N 1 ; ω i , t , i = 1 , , N 1 ] and u t = [ P m i , t , i = 1 , , N 1 ; Q T N i , t , i = 1 , , N 1 ] , it is possible to write:
x t + 1 = A x t + B u t + C h t + C ^ h ^ t t = 0 , , T 1
where:
A = [ A δ A δ ω A ω δ A ω ]
B = [ 0 0 B ω , m C T N ]
h t = [ I h ¯ t ]
h ^ t = [ 0 h ¯ ω , s , t ]
C = [ A δ , ω 0 0 C d ]
C ^ = [ 0 0 0 C ω , s ]

Appendix A.2. Proof of Proposition II

As regards the objective function, we consider the four terms in Equation (7):
  • The first term can be written as δ ¯ t T Q δ δ ¯ t .
  • The second term is given by ( ω ¯ t ω 0 ) T Q ω ( ω ¯ t ω 0 ) .
  • The fourth term is given by ( P ¯ m , t P ¯ m ) T R p ( P ¯ m , t P ¯ m ) .
With
P ¯ m = [ P m i , i = 1 , , N 1 ]
where Q δ Q ω and R p are the identity matrix of appropriate dimensions.
The third objective is the most complicated and it is given by (substituting voltage with Equation (A7)):
( B p 1 ( Q ¯ G , t + Q ¯ N G , t + Q ¯ R E S , t Q ¯ L , t K q ) B p 1 G δ ¯ t I ) T R v ( B p 1 ( Q ¯ G , t + Q ¯ N G , t + Q ¯ R E S , t Q ¯ L , t K q ) B p 1 G δ ¯ t I ) t = 0 , , T 1
Equation (A25) can be written as posing Y t = B p 1 ( Q ¯ G , t + Q ¯ R E S , t Q ¯ L , t + K q ) I :
( B p 1 Q ¯ N G , t B p 1 G δ ¯ t Y t ) T R v ( B p 1 Q ¯ N G , t B p 1 G δ ¯ t Y t ) t = 0 , , T 1
to reformulate Equation (A26), we can write, posing:
( B p 1 Q ¯ N G , t B p 1 G δ ¯ t Y t ) 2 = ( B p 1 G ) 2 ( B p 1 G 1 B p Q ¯ N G , t + δ ¯ t + G 1 B p Y t ) 2 t = 0 , , T 1
and posing Q ¯ T N , t = B p 1 G 1 B p Q ¯ N G , t , γ t = G 1 B p Y
( B p 1 G ) 2 ( Q ¯ T N , t + δ ¯ t γ t ) 2 = ( B p 1 G ) 2 ( δ ¯ t 2 + ( Q ¯ T N , t γ t ) 2 + 2 ( Q ¯ T N , t γ t ) δ ¯ t ) t = 0 , , T 1
and thus, posing M = ( B p 1 G ) 2 , the overall objective function can be expressed as:
J = t [ δ ¯ t T Q δ δ ¯ t + ( ω ¯ t ω 0 ) T Q ω ( ω ¯ t ω 0 ) + ( P ¯ m , t P ¯ m ) T R p ( P ¯ m , t P ¯ m ) + δ ¯ t T M δ ¯ t + ( Q ¯ T N , t γ t ) T M ( Q ¯ T N , t γ t ) + 2 M δ ¯ t ( Q ¯ T N , t γ t ) ]
Equation (A29) can be expressed in a vector form by using the following notation: x t = [ δ ¯ t ; ω ¯ t ] , u t = [ P ¯ m , t ; Q ¯ T N , t ] , x ^ = [ 0 ¯ ; ω ¯ 0 ] , u ^ t = [ P ¯ m ; γ t ]
J = t = 0 T 1 ( x t x ^ ) Q t ( x t x ^ ) + ( u t u ^ t ) R t ( u t u ^ t ) + 2 ( x t x ^ ) F t ( u t u ^ t )
with
Q = [ M + Q δ 0 0 Q ω ]
R = [ R p 0 0 M ]
F = [ 0 M 0 0 ]

Appendix A.3. Proof of Claim I

Due to its linearity, the system state equation can be written separating the stochastic ( z t + 1 s ) and deterministic ( z t + 1 d ) components [45]:
z t + 1 = z t + 1 s + z t + 1 d = A ( z t s + z t d ) + B ( e t s + e t d ) + μ t + C ^ h ^ t t = 0 , , T 1
thus, the system can be so decomposed into a stochastic and deterministic subsystem. That is:
z t + 1 s = A z t s + B e t s + C ^ h ^ t t = 0 , , T 1 z 0 s = 0
z t + 1 d = A z t d + B e t d + μ t t = 0 , , T 1 z 0 d = z 0
the objective function is:
min J ( z , e ) = 1 2 E { t = 0 T 1 [ ( z t + 1 s + z t + 1 d ) Q t ( z t + 1 s + z t + 1 d ) + ( e t s + e t d ) R t ( e t s + e t d ) + 2 ( z t + 1 s + z t + 1 d ) F t ( e t s + e t d ) ] + ( z t + 1 s + z t + 1 d ) Q T ( z t + 1 s + z t + 1 d ) }
Also, the cost function (Equation (A36)) can be decomposed into its stochastic and deterministic components represented by the functions J ( z s , e s ) and J ( z d , e d ) , defined straightforwardly. This means that for each admissible solution ( z , e ) , and the corresponding stochastic ( z s , e s ) and deterministic ( z d , e d ) components, the following identity holds:
J ( z , e ) = J ( z s , e s ) + J ( z d , e d ) .
This immediately follows from the assumptions on the stochastic process, C ^ h ^ t , which implies E { z t s } = 0 , and E { e t s } = 0 t . In fact, in this way, the expected values of terms z t s Q t z t d , z t d Q t z t s , e t d R t e t s , e t s R t e t d , z t d F t e t s , z t s F t z t d are equal to zero.
Thus, the overall optimal solution is given by:
z t = z t s , + z t d , t = 0 , , T 1
e t = e t s , + e t d , t = 0 , , T 1
being, of course ( z , d , e , d ) , ( z , s , e , s ) , the optimal solutions of the problem corresponding to the optimization of the deterministic and stochastic costs, respectively. In fact, for every admissible solution ( z , e ) , the following relationship can be written:
J ( z , e ) = J ( z s , e s ) + J ( z d , e d ) J ( z , s , e , s ) + J ( z , d , e , d ) = J ( z , e ) .

Appendix A.4. Proof of Claim II

The proof follows the results obtained in the continuous domain and reported in Reference [50].
The optimal control problem to determine e t d , is given by:
min J ( z , e ) = 1 2 E { t = 0 T 1 [ z t d Q t z t + e t d R t e t d + 2 z t d F t e t d ] + z T d Q T z T d }
s.t.
z t + 1 d = A z t d + B e t d + μ t t = 0 , , T 1
z 0 d = z 0
The deterministic problem can be introduced observing that z t d can be decomposed in two parts:
z t d = z t d 1 + z t d 2 t = 0 , , T 1
where:
z t + 1 d 1 = A z t d 1 + B e t d t = 0 , , T 1
z 0 d 1 = 0
z t + 1 d 2 = A z t d 2 + μ t t = 0 , , T 1
z 0 d 2 = z 0
The vector z t d 2 may be exactly determined for each instant t. This means that the deterministic component of the cost function can only be influenced by the values of z t d 1 . Then, putting z t d 2 = r t , the optimization of J ( z d 1 , e d ) is equivalent to the optimization of
min J ( z , e ) = 1 2 E { t = 0 T 1 ( z t d 1 r t ) Q t ( z t d 1 r t ) + e t d R t e t d + 2 ( z t d 1 r t ) F t e t d + ( z T d 1 r T ) Q T ( z T d 1 r T ) }

Appendix A.5. Proof of Proposition III

As reported in Reference [47], the following steps have been considered for the solution of discrete-time optimal control (Equations (35)–(37)):
  • Calculate the Hamiltonian function:
H ( z t s , e t s , λ t + 1 ) = 1 2 t = 0 T 1 [ z t s Q t z t s + e t s R t e t s + 2 z t s F t e t s ] + λ t + 1 ( A z t s + B e t s )
2.
Obtain state, co-state, and control from the Hamiltonian:
H ( z t s , e t s , λ t + 1 ) λ t + 1 = z t + 1 s , = A z t s , + B e t s , t = 0 , , T 1
H ( z t s , e t s , λ t + 1 ) z t = λ t = A λ t + 1 + Q t z t s , + F t e t s , t = 0 , , T 1
H ( z t s , e t s , λ t + 1 ) e t = 0 = B λ t + 1 + R   t e t s , + F t z t s , t = 0 , , T 1
3.
Obtain open-loop optimal control:
Solving Equation (A53) and replacing it into Equations (A51) and (A52):
e t s , = R t 1 B λ t + 1 R t 1 F t z t s , t = 0 , , T 1
z t + 1 s , = A z t s , E t λ t + 1 M t z t s , t = 0 , , T 1
λ t = A λ t + 1 + Q t z t s , M t λ t + 1 D t z t s , t = 0 , , T 1
For the sake of simplicity, we have defined M t = B R t 1 F t , D t = F t R t 1 F t and E t = B R t 1 B .
The final (boundary) condition is given by:
λ T = 1 2 ( z T s , Q T z T s , ) z T s , = Q T z T s , t = 0 , , T 1
4.
Obtain the Riccati equation:
We make the assumption:
λ t = P t z t s , t = 0 , , T 1
and replacing Equation (A62) into Equation (A59), we obtain:
z t + 1 s , = A z t s , E t P t + 1 z t + 1 s , M t z t s , t = 0 , , T 1
z t + 1 s , = ( I + E t P t + 1 ) 1 [ A z t s , M t z t s , ] t = 0 , , T 1
We then replace Equations (A62) and (A64) into Equation (A60):
P t z t s , = A t P t + 1 { ( I + E t P t + 1 ) 1 [ A z t s , M t z t s , ] } + Q t z t s , D t z t s , M t P t + 1 { ( I + E t P t + 1 ) 1 [ A z t s , M t z t s , ] } t = 0 , , T 1
this equation must hold for all values of the state z t , which in turn leads to the fact that the coefficient of z t in Equation (A61) must individually vanish. That is,
P t = Q t D t + ( A M t ) P t + 1 ( I + E t P t + 1 ) 1 ( A t M t ) t = 0 , , T 1
5.
Closed-loop optimal control:
Obtaining these solutions off-line, use Equation (A58) in the control relation Equation (A54) to get the closed-loop optimal control as
e t s , = R t 1 B P t + 1 ( A z t s , + B e t s , ) R t 1 F t z t s , t = 0 , , T 1
Now, pre-multiplying by R and solving for the optimal control, we have:
e t s = K t z t s t = 0 , , T 1
K t = ( R + B P t + 1 B ) 1 ( B P t + 1 A + F t ) t = 0 , , T 1
P t = Q t D t + ( A M t ) P t + 1 ( I + E t P t + 1 ) 1 ( A t M t ) t = 0 , , T 1
P T = Q T

Appendix A.6. Proof of Proposition IV

The proof follows the same steps of Proposition III. Specifically, the Hamiltonian function and related conditions are given by:
H ( z t d 1 , e t d , λ t + 1 ) = 1 2 t = 0 T 1 { ( z t d 1 r t ) Q t ( z t d 1 r t ) + e t d R t e t + 2 ( z t d 1 r t ) F t e t d } + λ t + 1 ( A z t d 1 + B e t d ) t = 0 , , T 1
H ( z t d 1 , e t d , λ t + 1 ) λ t + 1 = z t + 1 d 1 , = A z t d 1 , + B e t d , t = 0 , , T 1
H ( z t d 1 , e t d , λ t + 1 ) z t = λ t = A λ t + 1 + Q t z t d , Q t r t + F t e t t = 0 , , T 1
H ( z t d 1 , e t d , λ t + 1 ) e t = 0 = B λ t + 1 + R t e t d , + F t ( z t d 1 , r t ) t = 0 , , T 1
e t d , = R t 1 B λ t + 1 R t 1 F t ( z t d 1 , r t ) t = 0 , , T 1
z t + 1 d 1 , = A z t d 1 , E t λ t + 1 M t ( z t d 1 , r t ) t = 0 , , T 1
λ t = A λ t + 1 + Q t z t d 1 , Q t r t M t λ t + 1 D t ( z t d 1 , r t ) t = 0 , , T 1
λ T = 1 2 ( ( z T d 1 r T ) Q T ( z T d 1 r T ) ) z T d 1 = Q T z T d 1 Q T r T t = 0 , , T 1
being M = B R 1 F , D = F R 1 F and E = B R 1 B .
We may assume that λ t = P t z t d 1 , g t , and we obtain
z t + 1 d 1 , = A z t d 1 , E t P t + 1 z t + 1 d 1 , + E t g t + 1 M t ( z t d 1 , r t ) t = 0 , , T 1
z t + 1 d 1 , = ( I + E t P t + 1 ) 1 [ A z t d 1 , + E t g t + 1 M t ( z t d 1 , r t ) ] t = 0 , , T 1
P t z t d 1 , g t = A { P k + 1 { ( I + E t P t + 1 ) 1 [ A z t d 1 , + E t g t + 1 M t ( z t d 1 , r t ) ] } g t + 1 } + Q t z t d 1 , Q t r t D t ( z t d 1 , r t ) M t { P k + 1 { ( I + E t P t + 1 ) 1 [ A z t d 1 , + E t g t + 1 M t ( z t d 1 , r t ) ] } g t + 1 } t = 0 , , T 1
This equation must hold for all values of the state z t , which in turn leads to the fact that the coefficient of z t and the rest of the terms in Equation (A77) must individually vanish. That is
P t = Q t D t + ( A M t ) P t + 1 ( I + E t P t + 1 ) 1 ( A M t ) t = 0 , , T 1
g t = [ ( A M t ) P t + 1 ( I + E t P t + 1 ) 1 ( E t g t + 1 + M t r t ) ( A M t ) g t + 1 ] + ( Q t D t ) r t t = 0 , , T 1
Finally, the closed-loop optimal control is
e k d , = K t z t d 1 + K t g ( B g t + 1 + F t r t ) t = 0 , , T 1
with
K t = ( R t + B P t + 1 B ) 1 ( B P t + 1 A + F t ) t = 0 , , T 1
K t g = ( R t + B P t + 1 B ) 1
P t = Q t D t + ( A M t ) P t + 1 ( I + E t P t + 1 ) 1 ( A M t ) t = 0 , , T 1
r t + 1 = r t μ t t = 0 , , T 1
P T = Q T
g T = Q T r T

References

  1. Arya, Y.; Kumar, N. Optimal control strategy-based AGC of electrical power systems: A comparative performance analysis. Optim. Control Appl. Methods 2017, 38, 982–992. [Google Scholar] [CrossRef]
  2. Dagdougui, H.; Ouammi, A.; Sacile, R. Optimal control of a networkof power microgrids using the Pontryagin’s minimum principle. IEEE Trans. Control Syst. Technol. 1942, 22, 1942–1948. [Google Scholar] [CrossRef]
  3. Bersani, C.; Dagdougui, H.; Ouammi, A.; Sacile, R. Distributed Robust Control of the Power Flows in a Team of Cooperating Microgrids. IEEE Trans. Control Syst. Technol. 2016, 25, 1473–1479. [Google Scholar] [CrossRef]
  4. Ferro, G.; Minciardi, R.; Parodi, L.; Robba, M.; Rossi, M. Optimal Control of Multiple Microgrids and Buildings by an Aggregator. Energies 2020, 13, 1058. [Google Scholar] [CrossRef] [Green Version]
  5. Pandey, R.K.; Gupta, D.K. Intelligent Multi-Area Power Control: Dynamic Knowledge Domain Inference Concept. IEEE Trans. Power Syst. 2017, 32, 4310–4318. [Google Scholar] [CrossRef]
  6. Trip, S.; Cucuzzella, M.; De Persis, C.; Ferrara, A.; Scherpen, J.M. Robust load frequency control of nonlinear power networks. Int. J. Control 2018, 93, 346–359. [Google Scholar] [CrossRef]
  7. Datta, S. Small Signal Stability Criteria for Descriptor Form Power Network Model. Int. J. Control 2018, 2018, 1–16. [Google Scholar] [CrossRef]
  8. Dreidy, M.; Mokhlis, H.; Mekhilef, S. Inertia response and frequency control techniques for renewable energy sources: A review. Renew. Sustain. Energy Rev. 2017, 69, 144–155. [Google Scholar] [CrossRef]
  9. Rezaei, N.; Kalantar, M. Smart microgrid hierarchical frequency control ancillary service provision based on virtual inertia concept: An integrated demand response and droop controlled distributed generation framework. Energy Convers. Manag. 2015, 92, 287–301. [Google Scholar] [CrossRef]
  10. Hammad, E.; Farraj, A.; Kundur, D. On Effective Virtual Inertia of Storage-Based Distributed Control for Transient Stability. IEEE Trans. Smart Grid 2017, 10, 327–336. [Google Scholar] [CrossRef]
  11. Markovic, U.; Chu, Z.; Aristidou, P.; Hug, G. Fast Frequency Control Scheme through Adaptive Virtual Inertia Emulation. IEEE Trans. Smart Grid 2015, 6, 2627–2638. [Google Scholar] [CrossRef]
  12. Xu, Z.-R.; Yang, P.; Zheng, C.; Zhang, Y.; Peng, J.; Zeng, Z. Analysis on the organization and Development of multi-microgrids. Renew. Sustain. Energy Rev. 2018, 81, 2204–2216. [Google Scholar] [CrossRef]
  13. Zheng, Y.; Zhou, J.; Xu, Y.; Zhang, Y.; Qian, Z. A distributed model predictive control based load frequency control scheme for multi-area interconnected power system using discrete-time Laguerre functions. ISA Trans. 2017, 68, 127–140. [Google Scholar] [CrossRef]
  14. Khooban, M.-H.; Niknam, T. A new intelligent online fuzzy tuning approach for multi-area load frequency control: Self Adaptive Modified Bat Algorithm. Int. J. Electr. Power Energy Syst. 2015, 71, 254–261. [Google Scholar] [CrossRef]
  15. Chen, M.-R.; Zeng, G.; Dai, Y.-X.; Lu, K.-D.; Bi, D.-Q. Fractional-Order Model Predictive Frequency Control of an Islanded Microgrid. Energies 2018, 12, 84. [Google Scholar] [CrossRef] [Green Version]
  16. Guo, Z.; Li, S.; Zheng, Y. Feedback Linearization Based Distributed Model Predictive Control for Secondary Control of Islanded Microgrid. Asian J. Control 2018, 22, 460–473. [Google Scholar] [CrossRef]
  17. Falahi, M.; Butler-Purry, K.L.; Ehsani, M. Dynamic Reactive Power Control of Islanded Microgrids. IEEE Trans. Power Syst. 2013, 28, 3649–3657. [Google Scholar] [CrossRef]
  18. Yuan, C.; Xie, P.; Yang, D.; Xiao, X. Transient Stability Analysis of Islanded AC Microgrids with a Significant Share of Virtual Synchronous Generators. Energies 2018, 11, 44. [Google Scholar] [CrossRef] [Green Version]
  19. Marinovici, L.D.; Lian, J.; Kalsi, K.; Du, P.; Elizondo, M. Distributed Hierarchical Control Architecture for Transient Dynamics Improvement in Power Systems. IEEE Trans. Power Syst. 2013, 28, 3065–3074. [Google Scholar] [CrossRef]
  20. Lu, W.; Liu, M.; Lin, S.; Li, L. Fully Decentralized Optimal Power Flow of Multi-Area Interconnected Power Systems Based on Distributed Interior Point Method. IEEE Trans. Power Syst. 2018, 33, 901–910. [Google Scholar] [CrossRef]
  21. Im, W.-S.; Wang, C.; Liu, W.; Liu, L.; Kim, J.-M. Distributed virtual inertia based control of multiple photovoltaic systems in autonomous microgrid. IEEE/CAA J. Autom. Sin. 2016, 4, 512–519. [Google Scholar] [CrossRef]
  22. Wu, D.; Tang, F.; Dragicevic, T.; Guerrero, J.M.; Vasquez, J. Coordinated Control Based on Bus-Signaling and Virtual Inertia for Islanded DC Microgrids. IEEE Trans. Smart Grid 2015, 6, 1. [Google Scholar] [CrossRef] [Green Version]
  23. Asghar, F.; Talha, M.; Kim, S.H. Robust Frequency and Voltage Stability Control Strategy for Standalone AC/DC Hybrid Microgrid. Energies 2017, 10, 760. [Google Scholar] [CrossRef] [Green Version]
  24. Asghar, F.; Talha, M.; Kim, S.H. Fuzzy logic-based intelligent frequency and voltage stability control system for standalone microgrid. Int. Trans. Electr. Energy Syst. 2017, 28, e2510. [Google Scholar] [CrossRef]
  25. Jumani, T.A.; Mustafa, M.W.; Rasid, M.; Mirjat, N.; Leghari, Z.H.; Saeed, M.S. Optimal Voltage and Frequency Control of an Islanded Microgrid using Grasshopper Optimization Algorithm. Energies 2018, 11, 3191. [Google Scholar] [CrossRef] [Green Version]
  26. Qi, X.; Bai, Y.; Luo, H.; Zhang, Y.; Zhou, G.; Wei, Z. Fully-distributed Load Frequency Control Strategy in an Islanded Microgrid Considering Plug-In Electric Vehicles. Energies 2018, 11, 1613. [Google Scholar] [CrossRef] [Green Version]
  27. Bouzid, A.E.M.; Sicard, P.; Chaoui, H.; Cheriti, A.; Sechilariu, M.; Guerrero, J. A novel Decoupled Trigonometric Saturated droop controller for power sharing in islanded low-voltage microgrids. Electr. Power Syst. Res. 2019, 168, 146–161. [Google Scholar] [CrossRef]
  28. Wu, J.; Liu, T.; Qiu, T.; Xu, D. Current and power quality multi-objective control of virtual synchronous generators under unbalanced grid conditions. J. Power Electron. 2020, 20, 580–589. [Google Scholar] [CrossRef]
  29. Deng, J.; Xia, N.; Yin, J.; Jin, J.; Peng, S.; Wang, T. Small-Signal Modeling and Parameter Optimization Design for Photovoltaic Virtual Synchronous Generator. Energies 2020, 13, 398. [Google Scholar] [CrossRef] [Green Version]
  30. Alizadeh, G.A.; Rahimi, T.; Nozadian, M.H.B.; Padmanaban, S.; Leonowicz, Z. Improving Microgrid Frequency Regulation Based on the Virtual Inertia Concept while Considering Communication System Delay. Energies 2019, 12, 2016. [Google Scholar] [CrossRef] [Green Version]
  31. Fini, M.H.; Golshan, M.E.H. Determining optimal virtual inertia and frequency control parameters to preserve the frequency stability in islanded microgrids with high penetration of renewables. Electr. Power Syst. Res. 2018, 154, 13–22. [Google Scholar] [CrossRef]
  32. Kerdphol, T.; Rahman, F.S.; Mitani, Y.; Watanabe, M.; Kufeoglu, S. Robust Virtual Inertia Control of an Islanded Microgrid Considering High Penetration of Renewable Energy. IEEE Access 2018, 6, 625–636. [Google Scholar] [CrossRef]
  33. Vergara, P.P.; Rey, J.M.; Shaker, H.R.; Guerrero, J.M.; Jørgensen, B.N.; Da Silva, L.C.P. Distributed Strategy for Optimal Dispatch of Unbalanced Three-Phase Islanded Microgrids. IEEE Trans. Smart Grid 2019, 10, 3210–3225. [Google Scholar] [CrossRef]
  34. Minciardi, R.; Sacile, R. Optimal Control in a Cooperative Network of Smart Power Grids. IEEE Syst. J. 2011, 6, 126–133. [Google Scholar] [CrossRef]
  35. Dagdougui, H.; Sacile, R. Decentralized Control of the Power Flows in a Network of Smart Microgrids Modeled as a Team of Cooperative Agents. IEEE Trans. Control Syst. Technol. 2013, 22, 510–519. [Google Scholar] [CrossRef]
  36. Delfino, F.; Minciardi, R.; Pampararo, F.; Robba, M. A Multilevel Approach for the Optimal Control of Distributed Energy Resources and Storage. IEEE Trans. Smart Grid 2014, 5, 2155–2162. [Google Scholar] [CrossRef]
  37. Minciardi, R.; Robba, M. A Bilevel Approach for the Stochastic Optimal Operation of Interconnected Microgrids. IEEE Trans. Autom. Sci. Eng. 2016, 14, 1–12. [Google Scholar] [CrossRef]
  38. Rahmani, M.; Sadati, N. Hierarchical optimal robust load-frequency control for power systems. IET Gener. Transm. Distrib. 2012, 6, 303. [Google Scholar] [CrossRef]
  39. Bejestani, A.K.; Annaswamy, A.; Samad, T. A Hierarchical Transactive Control Architecture for Renewables Integration in Smart Grids: Analytical Modeling and Stability. IEEE Trans. Smart Grid 2014, 5, 2054–2065. [Google Scholar] [CrossRef]
  40. Ghosh, S.; Kamalasadan, S. An Energy Function-Based Optimal Control Strategy for Output Stabilization of Integrated DFIG-Flywheel Energy Storage System. IEEE Trans. Smart Grid 2017, 8, 1922–1931. [Google Scholar] [CrossRef]
  41. Nguyen, D.B.; Scherpen, J.M.; Bliek, F. Distributed Optimal Control of Smart Electricity Grids With Congestion Management. IEEE Trans. Autom. Sci. Eng. 2017, 14, 494–504. [Google Scholar] [CrossRef] [Green Version]
  42. Cai, N. Linearized and Distributed Methods for Power Flow Analysis and Control in Smart Grids and Microgrids; Michigan State University: East Lansing, MI, USA, 2014. [Google Scholar]
  43. Delfino, F.; Ferro, G.; Robba, M.; Rossi, M. An Energy Management Platform for the Optimal Control of Active and Reactive Powers in Sustainable Microgrids. IEEE Trans. Ind. Appl. 2019, 55, 7146–7156. [Google Scholar] [CrossRef]
  44. Yang, Z.; Zhong, H.; Xia, Q.; Bose, A.; Kang, C. Optimal power flow based on successive linear approximation of power flow equations. IET Gener. Transm. Distrib. 2016, 10, 3654–3662. [Google Scholar] [CrossRef]
  45. Bruni, C.; Iacoviello, D. Some results about the optimal LQG tracking problem. Int. J. Control 2001, 74, 977–987. [Google Scholar] [CrossRef]
  46. Pindyck, R. An application of the linear quadratic tracking problem to economic stabilization policy. IEEE Trans. Autom. Control 1972, 17, 287–300. [Google Scholar] [CrossRef]
  47. Kumar R., S.; Jomoah, I.M.; Bafail, A.O. Optimal Placement of Unified Power Flow Controller for Minimization of Power Transmission Line Losses. Int. J. Comput. Theory Eng. 2014, 6, 377–381. [Google Scholar] [CrossRef] [Green Version]
  48. Zhang, J.; Verschae, R.; Nobuhara, S.; LaLonde, J.-F. Deep photovoltaic nowcasting. Sol. Energy 2018, 176, 267–276. [Google Scholar] [CrossRef] [Green Version]
  49. Oneto, L.; Laureri, F.; Robba, M.; Delfino, F.; Anguita, D. Data-Driven Photovoltaic Power Production Nowcasting and Forecasting for Polygeneration Microgrids. IEEE Syst. J. 2017, 12, 2842–2853. [Google Scholar] [CrossRef]
  50. Naidu, D. Optimal Control Systems; CRC Press: Boca Raton, FL, USA, 2003. [Google Scholar]
Figure 1. An interconnected local systems power system.
Figure 1. An interconnected local systems power system.
Energies 13 02637 g001
Figure 2. Local systems’ variables detail.
Figure 2. Local systems’ variables detail.
Energies 13 02637 g002
Figure 3. Model predictive control (MPC) scheme.
Figure 3. Model predictive control (MPC) scheme.
Energies 13 02637 g003
Figure 4. Modified IEEE 5 bus test system.
Figure 4. Modified IEEE 5 bus test system.
Energies 13 02637 g004
Figure 5. Deterministic part of the net load, P s c i , t d (difference between renewable generation and load), data on each local system.
Figure 5. Deterministic part of the net load, P s c i , t d (difference between renewable generation and load), data on each local system.
Energies 13 02637 g005
Figure 6. Phase optimal results at each local system.
Figure 6. Phase optimal results at each local system.
Energies 13 02637 g006
Figure 7. Frequency optimal results at each local system.
Figure 7. Frequency optimal results at each local system.
Energies 13 02637 g007
Figure 8. Voltage optimal results at each local system.
Figure 8. Voltage optimal results at each local system.
Energies 13 02637 g008
Figure 9. Optimal active power at each local system.
Figure 9. Optimal active power at each local system.
Energies 13 02637 g009
Figure 10. Optimal reactive power at each local system.
Figure 10. Optimal reactive power at each local system.
Energies 13 02637 g010
Figure 11. Sensitivity analysis of the stability property for the proposed system.
Figure 11. Sensitivity analysis of the stability property for the proposed system.
Energies 13 02637 g011
Table 1. Lines data.
Table 1. Lines data.
LineResistance Ri,j (p.u.)Reactance Xi,j (p.u.)
1–20.0880.266
1–30.3551.066
2–30.2660.800
2–40.2660.800
2–50.1770.533
3–40.1770.533
4–50.3551.066
Table 2. Local systems’ parameters.
Table 2. Local systems’ parameters.
ParameterLocal System 1Local System 2Local System 3Local System 4
KDi (p.u)343.55
Hi (s)545.56
Table 3. Objective function values and contributions.
Table 3. Objective function values and contributions.
ContributionPhaseFrequencyVoltageActive Power
No Control0.008548.71 × −40.00950
With Control0.008538.64 × −40.00420.0015

Share and Cite

MDPI and ACS Style

Ferro, G.; Robba, M.; Sacile, R. A Model Predictive Control Strategy for Distribution Grids: Voltage and Frequency Regulation for Islanded Mode Operation. Energies 2020, 13, 2637. https://doi.org/10.3390/en13102637

AMA Style

Ferro G, Robba M, Sacile R. A Model Predictive Control Strategy for Distribution Grids: Voltage and Frequency Regulation for Islanded Mode Operation. Energies. 2020; 13(10):2637. https://doi.org/10.3390/en13102637

Chicago/Turabian Style

Ferro, Giulio, Michela Robba, and Roberto Sacile. 2020. "A Model Predictive Control Strategy for Distribution Grids: Voltage and Frequency Regulation for Islanded Mode Operation" Energies 13, no. 10: 2637. https://doi.org/10.3390/en13102637

APA Style

Ferro, G., Robba, M., & Sacile, R. (2020). A Model Predictive Control Strategy for Distribution Grids: Voltage and Frequency Regulation for Islanded Mode Operation. Energies, 13(10), 2637. https://doi.org/10.3390/en13102637

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop