This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Decentralized Droop-based Finite-Control-Set Model Predictive Control of Inverter-based Resources in Islanded AC Microgrid

Ayobami Olajube, Koto Omiloli, Satish Vedula, and
Olugbenga Moses Anubi

Department of Electrical and Computer Engineering, the Center for Advanced Power Systems, Florida State University E-mail: {aolajube, kao23a , svedula , oanubi}@fsu.edu

1 Abstract

This paper presents an improved droop control method to ensure effective power sharing, voltage regulation, and frequency stabilization of inverter-based resources (IBRs) connected in parallel in an islanded AC microgrid. In the contemporary droop control algorithm, the distance between connected inverters affects the effectiveness of the active power-frequency and the reactive power-voltage droop characteristics which results in poor power sharing at the primary level of the microgrid. That is, high impedance emanating from long transmission lines results in instability, poor voltage tracking, and ineffective frequency regulation. Hence, in this work, we use a finite-control-set model predictive controller (FCS-MPC) in the inner loop, which gives efficient voltage tracking, good frequency regulation, and faster performance response. FCS-MPC is easy to implement in fast switching converters and does not suffer from computational burden unlike the continuous-set MPC and is also devoid of issues of multiple-loop, parameter variation, and slow response associated with conventional droop control methods. We derived the condition for bounded stability for the FCS-MPC and the proposed method is tested via a numerical simulation on three IBRs. The results show effective power sharing, capacitor voltage tracking, and efficient frequency regulation with reduced oscillations to changes in load.

Keywords: Droop control, finite-control-set model predictive control, active power sharing, reactive power sharing, voltage regulation, inverter-based resources, islanded microgrid.

2 Introduction

Traditional Microgrids (MGs) are a conglomeration of numerous fossil-fuel-driven power generation elements like synchronous machines (SMs). However, in recent times, an increase in penetration of inverter-based resources (IBRs) such as solar and wind generation has caused a shift in power generation trends in the existing infrastructure. The addition of energy storage elements (ESEs) further upped the choice of power allocation in the MGs in terms of the available degree of freedom. Based on their configuration and the mode of operation MGs are classified into two categories known as grid-connected and islanded modes [1]. Control of the traditional MG employs a hierarchical control structure [2], with droop-control being the lowest level followed by the power and the energy management (EM) layers. The impending integration of IBRs into the traditional MGs is giving rise to structural and operational changes in the form of transitioning away from the SMs which are ramp-limited due to their inertial properties. However, the low inertia and the fast switching properties of the IBRs pose a major challenge in the form of assessing their stability and robustness [3]. One such method that is widely embraced to restore the inertia to the grid is to introduce virtual inertia [4, 5].

Regardless of the inertial properties of the SMs and the IBRs, the control objective remains unequivocal in terms of balancing the active and reactive power in the MGs while simultaneously maintaining the grid frequency around the desired operating point. Droop control is the method widely employed in addressing the control of the power and the frequency namely voltage droop and frequency droop [6]. The IBRs are operated in parallel mode due to high reliability and flexibility in terms of operation. However, the application of droop control in such a parallel setting is not suitable while sharing a nonlinear load due to the presence of harmonic currents [7]. Moreover, the frequency and amplitude deviation in the IBRs droop settings are dependent on the load scenarios [8]. Addition of virtual control loops such as adding a virtual resistance into the droop methods and decentralizing the control (consensus-based) are proposed to address the complications arising from the presence of current harmonics [9]. Although Communication-based and consensus control algorithms offer efficient voltage regulation and effective power sharing among connected IBRs, the methods suffer from the huge cost and reliability of communication links. On the other hand, droop control does not require communication between IBRs to effect power sharing and voltage regulation [10]. Multiple loops of conventional droop techniques often result in poor power-sharing due to parameter variations, especially in the current loop that generates PWM signals for the inverters [11].

Thus, in this paper, we utilize FCS-MPC in the inner loop to directly generate optimal switching signals for the IBRs. Finite-control-set model predictive control (FCS-MPC) has emerged as a promising discrete-time model-based nonlinear predictive control algorithm applicable to fast-switching power electronics dominated systems. Unlike the continuous-set MPC that requires modulators and a longer horizon to generate the pulse-width modulated (PWM) signals, FCS-MPC eliminates the need for modulators and accounts for nonlinearities in the power converter dynamics [12]. It uses the horizon of h = 1, which significantly reduces the computational burden that may result from very low sampling time in high-frequency converters. To this end, we develop a decentralized droop-based FCS-MPC-based control algorithm to ensure adequate active and reactive power sharing between the parallel connected IBRs and the accurate tracking of the filter output filter capacitor voltage. The droop control method provides a reference current for the FCS-MPC algorithm which in turn generates the PWM signal to switch the parallel connected IBRs. We examine the behaviors of the IBRs with respect to their droop coefficients. Also, the effect of sudden changes in the load on the droop characteristics of the inverters is examined. Consequently, the key contributions in this paper are:

  1. 1.

    The unique stability analysis for finite-set model predictive control (FCS-MPC) algorithm.

  2. 2.

    The development of the decentralized droop-based finite-control-set model predictive control of several IBRs to regulate the AC grid voltage error to zero without using PWM modulators.

  3. 3.

    The evaluation of the performance of the droop-based FCS-MPC on both P-f and Q-V droop characteristics through simulations on three parallel connected IBRs.

The rest of the paper is organized as follows: section II shows the preliminaries and notations used throughout this paper, section III shows the model development of decentralized inverter-based resources in an islanded AC grid, section IV depicts the control design and algorithm development, section V enumerates the numerical simulation results and the conclusion follows in section VI.

3 Notations and Preliminaries

In this paper, we use n\mathbb{R}^{n}, n×m\mathbb{R}^{n\times m} to denote the space of real numbers, real vectors of length nn, and real matrices of nn rows, and mm columns, respectively. +\mathbb{R}_{+} indicates a set of positive real numbers. XX^{\top} denotes the transpose of the quantity XX. Normal-face lower-case letters (xx\in\mathbb{R}) are used to represent real scalars, bold-face lower-case letters (𝐱n\mathbf{x}\in\mathbb{R}^{n}) represent vectors, while normal-face upper case (Xn×mX\in\mathbb{R}^{n\times m}) represents matrices. X0(0)X\succ 0\hskip 2.84526pt(\succeq 0) denotes a positive definite (semi-definite) matrix. Furthermore, J=[0110]J=\begin{bmatrix}0&1\\ -1&0\\ \end{bmatrix} is a skew-symmetric matrix, MN=[M11NM12NM21NM22N]M\otimes N=\begin{bmatrix}M_{11}N&M_{12}N\\ M_{21}N&M_{22}N\\ \end{bmatrix} denotes the Kronecker product of matrices MM and NN, 𝟏n\mathbf{1}_{n} is a column vector of ones, and the identity matrix of size nn is denoted by InI_{n}.

3.1 Stability Analysis of Finite-Set Predictive Control System

Consider a discrete-time system described by

x(k+1)=f(x(k),u(k))\textbf{x}(k+1)=f(\textbf{x}(k),\textbf{u}(k)) (1)

Where x(k)𝒳n\textbf{x}(k)\in\mathcal{X}\subset\mathbb{R}^{n} is the system state and u(k)𝒰m\textbf{u}(k)\in\mathcal{U}\subset\mathbb{R}^{m} is the control input. If we assume the equilibrium point exists at the origin, we have f(0,0)=0f(0,0)=0. Suppose there exists positive constants α1\alpha_{1}, α2\alpha_{2}, α3\alpha_{3}, α4\alpha_{4} and a positive definitive function V:nV:\mathbb{R}^{n}\rightarrow\mathbb{R} such that

α1x2V(x)α2x2V(f(x,u¯(x))V(x)α3x2|V(x)V(y)|α4(x+y)2α4z¯xy2α4z¯xy,\begin{array}[]{c}\alpha_{1}\left\lVert\textbf{x}\right\rVert^{2}\leq V(\textbf{x})\leq\alpha_{2}\left\lVert\textbf{x}\right\rVert^{2}\\ V(f(\textbf{x},\bar{u}(\textbf{x}))-V(\textbf{x})\leq-\alpha_{3}\left\lVert\textbf{x}\right\rVert^{2}\\ |V(\textbf{x})-V(\textbf{y})|\leq\underbrace{\alpha_{4}(\left\lVert\textbf{x}\right\rVert+\left\lVert\textbf{y}\right\rVert)}_{2\alpha_{4}\bar{\textbf{z}}}\left\lVert\textbf{x}-\textbf{y}\right\rVert\leq 2\alpha_{4}\bar{\textbf{z}}\left\lVert\textbf{x}-\textbf{y}\right\rVert,\end{array} (2)

where

z¯=supx𝒳x\bar{\textbf{z}}=\sup_{\textbf{x}\in\mathcal{X}}\left\lVert\textbf{x}\right\rVert

for all x𝒳n\textbf{x}\in\mathcal{X}\subset\mathbb{R}^{n}, where u¯:𝒰\bar{u}:\mathbb{R}\rightarrow\mathcal{U} is the associated feedback control and 𝒰m\mathcal{U}\subset\mathbb{R}^{m} and 𝒳n\mathcal{X}\subset\mathbb{R}^{n} are compact subsets. It follows that the origin of (1) is asymptotically stable under the control law u¯𝒰\bar{\textbf{u}}\in\mathcal{U} and 𝒳\mathcal{X} is positively invariant with respect to the closed loop dynamics.

Now let us={u1,u2,,uM}𝒰\textbf{u}_{s}=\left\{u_{1},u_{2},\ldots,u_{M}\right\}\subseteq\mathcal{U} be a finite sampling of the compact set 𝒰\mathcal{U} such that

supu𝒰min1iMuui2ϵ\sup_{\textbf{u}\in\mathcal{U}}\min_{1\ll i\leq M}\left\|\textbf{u}-\textbf{u}_{i}\right\|_{2}\leq\epsilon (3)

for some ϵ>0\epsilon>0. Here, ϵ\epsilon is the quantization error. We are interested in the stability of (1) using the control law

u(x)=argmin1iMu¯(x)ui2\textbf{u}(x)=\underset{1\leq i\leq M}{\operatorname{argmin}}\left\|\bar{\textbf{u}}(\textbf{x})-\textbf{u}_{i}\right\|_{2} (4)

Clearly, the projection error u(x)u¯(x)2ϵ\left\|{\textbf{u}}(\textbf{x})-\bar{\textbf{u}}(\textbf{x})\right\|_{2}\leq\epsilon follows from (3).

Theorem 1

Suppose there exists a positive definite function V:nV:\mathbb{R}^{n}\rightarrow\mathbb{R} satisfying the inequalities in (2) for some α1,α2,α3,α4+\alpha_{1},\alpha_{2},\alpha_{3},\alpha_{4}\in\mathbb{R}_{+} such that

α3α1<2,\displaystyle\frac{\alpha_{3}}{\alpha_{1}}<2, (5)

Then, the origin of the closed loop system x(k+1)=f(x(k),u(k))\textbf{x}(k+1)=f(\textbf{x}(k),\textbf{u}(k)) under the finite-set control law in (4) is stable and

limkx(k)2Lvα4z¯α3ϵ,\lim_{k\rightarrow\infty}\|\textbf{x}(k)\|\leq\sqrt{\frac{2L_{v}\alpha_{4}\bar{\textbf{z}}}{\alpha_{3}}\epsilon},

where LvL_{v} is a Lipschitz constant of the function f(x)f(x) in (1) and z¯\bar{\textbf{z}} is related to compact set 𝒳n\mathcal{X}\subset\mathbb{R}^{n} as given in (2).

Proof: Let us consider a lyapunov function V(x) and a Lipschitz continuous function f(x) such that:

V(f(x,u(x)))V(x)\displaystyle V(f(\textbf{x},\textbf{u}(\textbf{x})))-V(\textbf{x}) =V(f(x,u(x)))+\displaystyle=V(f(\textbf{x},\textbf{u}(\textbf{x})))+
V(f(x,u¯(x)))+V(f(x,u(x)))\displaystyle V(f(\textbf{x},\bar{\textbf{u}}(\textbf{x})))+V(f(\textbf{x},u(\textbf{x})))
V(f(x,u¯(x)))V(x)\displaystyle-V(f(\textbf{x},\bar{\textbf{u}}(\textbf{x})))-V(\textbf{x})
Lvf(x,u(x))f(x,u¯(x))\displaystyle\leq L_{v}\left\lVert f(\textbf{x},u(\textbf{x}))-f(\textbf{x},\bar{\textbf{u}}(\textbf{x}))\right\rVert
α3x2\displaystyle-\alpha_{3}\left\lVert\textbf{x}\right\rVert^{2}
2Lvα4z¯u(x)u¯(x)α3x2\displaystyle\leq 2L_{v}\alpha_{4}\bar{\textbf{z}}\left\lVert\textbf{u}(\textbf{x})-\bar{\textbf{u}}(\textbf{x})\right\rVert-\alpha_{3}\left\lVert\textbf{x}\right\rVert^{2}
2ϵLvα4z¯α3x2.\displaystyle\leq 2\epsilon L_{v}\alpha_{4}\bar{\textbf{z}}-\alpha_{3}\left\lVert\textbf{x}\right\rVert^{2}.

For all x𝒳\textbf{x}\in\mathcal{X}. Let VkV(x(k))V_{k}\triangleq V(\textbf{x}(k)), thus

Vk+1Vk\displaystyle V_{k+1}-V_{k} 2ϵLvα4z¯α3xk2\displaystyle\leq 2\epsilon L_{v}\alpha_{4}\bar{\textbf{z}}-\alpha_{3}\left\lVert\textbf{x}_{k}\right\rVert^{2}
2ϵLvα4z¯α3α1Vk.\displaystyle\leq 2\epsilon L_{v}\alpha_{4}\bar{\textbf{z}}-\frac{\alpha_{3}}{\alpha_{1}}V_{k}.

This implies that

Vk+1(1α3α1)Vk+2ϵLvα4z¯V_{k+1}\leq(1-\frac{\alpha_{3}}{\alpha_{1}})V_{k}+2\epsilon L_{v}\alpha_{4}\bar{\textbf{{z}}}
Vk(1α3α1)kVo+2ϵLvα4z¯i=0k1(1α3α1)i.V_{k}\leq(1-\frac{\alpha_{3}}{\alpha_{1}})^{k}V_{o}+2\epsilon L_{v}\alpha_{4}\bar{\textbf{z}}\sum_{i=0}^{k-1}(1-\frac{\alpha_{3}}{\alpha_{1}})^{i}.

Hence, V(xk)V(\textbf{x}_{k}) is bounded and, assuming bounded initial conditions

limkV(xk)2Lvα1α4z¯α3ϵ,\lim_{k\rightarrow\infty}V(\textbf{x}_{k})\leq\frac{2L_{v}\alpha_{1}\alpha_{4}\bar{\textbf{z}}}{\alpha_{3}}\epsilon,

which implies that

limkx(k)2Lvα4z¯α3ϵ\lim_{k\rightarrow\infty}\|\textbf{x}(k)\|\leq\sqrt{\frac{2L_{v}\alpha_{4}\bar{\textbf{z}}}{\alpha_{3}}\epsilon}

\square

3.2 Communicationless design of Decentralized IBRs

We consider an autonomous AC microgrid comprising inverter-based resources (IBRs) connected in parallel to supply power to constant power loads. It is represented as a graph 𝒢=(𝒱,,ik)\mathcal{G}=(\mathcal{V},\mathcal{E},\mathcal{B}_{ik}), where 𝒱={v|vin}\mathcal{V}=\{v|v_{i}\in\mathbb{R}^{n}\} represents a set of IBRs and 𝒱×𝒱={1,.,m}\mathcal{E}\in\mathcal{V}\times\mathcal{V}=\{1,....,m\} is a set of transmission lines connecting the IBRs. The incident matrix ikn×m\mathcal{B}_{ik}\in\mathbb{R}^{n\times m} showing the interconnections and current flow between the IBRs is shown below:

ik={+1, if IBR i starts i to k,1, if IBR i ends i to k,0, otherwise. \mathcal{B}_{ik}=\begin{cases}+1,&\text{ if IBR }i\text{ starts {i} to }k,\\ -1,&\text{ if IBR }i\text{ ends {i} to }k,\\ 0,&\text{ otherwise. }\end{cases}

When current flows from the i-th IBR to the k-th IBR, we assign a +1 sign and when current flows in the reverse direction, -1 is assigned. But, when there is no transmission line between any two inverters, we use zero.

4 Model Development

A hypothetical autonomous AC microgrid powered by IBRs is shown in Figure 1. Each IBR is connected to the load through an LC filter, which helps to dampen harmonics and ensure accurate extraction of the fundamental frequency. An IBR shares current with the neighboring IBR through the resistive-inductive transmission line.

Remark 1

The decentralized IBR dynamics in (LABEL:dq0_eqn) is a balanced and symmetrical system. Hence, the zero component of the dq0 rotating reference frame is zero. Then, we will stick to dq components only for the rest of this paper.

The following important assumptions are made to establish the generalized state-space representation of the system shown in Figure 1:

Assumption 1

For n homogeneous number of IBRs, switching losses are neglected during power transformation from DC to AC.

Assumption 2

The DC power supply to each IBR is stiff and the output voltages 𝐮dq\mathbf{u}_{dq} are local and unique to each inverter [13].

Assumption 3

The transmission line between each IBR is highly inductive (X\ggR) [14]. The power angle or phase shift angle under this condition becomes very small.

Assumption 4

The load currents iLdq\textbf{i}_{L{dq}} at the PCC are unknown, bounded, and obtainable from load voltage measurements at each bus.

Refer to caption
Figure 1: Decentralized Inverter-based Resources

The model of the decentralized scheme shown in Figure 1 in dq0 reference frame is given as:

Ld𝐢dqdt=𝐯dq(ωg(JL)+RIn)𝐢dq+𝐮dq,\displaystyle L\frac{d\mathbf{i}_{dq}}{dt}=-\mathbf{v}_{\mathrm{dq}}-(\omega_{g}(J\otimes L)+RI_{n})\mathbf{i}_{\mathrm{dq}}+\mathbf{u}_{\mathrm{dq}}, (6)
Csd𝐯dqdt=ωg(JCs)𝐯dq+𝐢dq+ik𝐢tdq𝐢Ldq,\displaystyle C_{s}\frac{d\mathbf{v}_{dq}}{dt}=-\omega_{g}(J\otimes C_{s})\mathbf{v}_{\mathrm{dq}}+\mathbf{i}_{\mathrm{dq}}+\mathcal{B}_{\mathrm{ik}}\mathbf{i}_{\mathrm{tdq}}-\mathbf{i}_{Ldq},
Ltd𝐢tdqdt=ik𝐯dq(ωg(JLt)+RIm)𝐢tdq,\displaystyle L_{t}\frac{d\mathbf{i}_{t{dq}}}{dt}=-\mathcal{B}_{\mathrm{ik}}^{\top}\mathbf{v}_{\mathrm{dq}}-(\omega_{g}(J\otimes L_{t})+RI_{m})\mathbf{i}_{\mathrm{tdq}},
𝐲dq=[𝐯dq;𝐢dq]\displaystyle\mathbf{y}_{\mathrm{dq}}=[\mathbf{v}_{\mathrm{dq}};\mathbf{i}_{\mathrm{dq}}]

where 𝐢dq=[𝐢d,𝐢q],𝐯dq=[𝐯d,𝐯q],𝐮dq=[𝐮d,𝐮q],𝐲dq=[𝐲d,𝐲q]2n\mathbf{i}_{dq}=[\mathbf{i}_{d}^{\top},\mathbf{i}_{q}^{\top}]^{\top},\mathbf{v}_{dq}=[\mathbf{v}_{d}^{\top},\mathbf{v}_{q}^{\top}]^{\top},\mathbf{u}_{dq}=[\mathbf{u}_{d}^{\top},\mathbf{u}_{q}^{\top}]^{\top},\newline \mathbf{y}_{dq}=[\mathbf{y}_{d},\mathbf{y}_{q}]^{\top}\in\mathbb{R}^{2n}, 𝐢Ldq=[𝐢Ld,𝐢Lq]2n\mathbf{i}_{L{dq}}=[\mathbf{i}_{Ld}^{\top},\mathbf{i}_{Lq}^{\top}]^{\top}\in\mathbb{R}^{2n} and 𝐢tdq=[𝐢td,𝐢tq]2m\mathbf{i}_{t{dq}}=[\mathbf{i}_{td}^{\top},\mathbf{i}_{tq}^{\top}]^{\top}\in\mathbb{R}^{2m}. Then, the filter resistance matrix becomes R=diag(r1,r2,,rn)R=\textsf{diag}(r_{1},r_{2},\ldots,r_{n}), filter inductance matrix L=diag(l1,l2,,ln)L=\textsf{diag}(l_{1},l_{2},\ldots,l_{n}), the filter capacitor matrix Cs=diag(C1,C2,,Cn)C_{s}=\textsf{diag}(C_{1},C_{2},\ldots,C_{n}), the transmission line resistance matrix Rt=diag(rt1,rt2,,rtm)R_{t}=\textsf{diag}(r_{t1},r_{t2},\ldots,r_{tm}) and transmission line inductance matrix Lt=diag(lt1,lt2,,ltm)L_{t}=\textsf{diag}(l_{t1},l_{t2},\ldots,l_{tm}).

Hence, the state space dynamics of the entire system in equation (LABEL:dq0_eqn) is given as:

x˙(t)=Ax(t)+Bu(t)+EiLdq(t),\displaystyle\dot{\textbf{x}}(t)=\mathrm{A\textbf{x}(t)}+\mathrm{B\textbf{u}(t)}+\mathrm{E\textbf{i}_{L{dq}}(t)}, (7)
y(t)=Cx(t),\displaystyle\mathrm{\textbf{y}(t)}=\mathrm{C\textbf{x}(t)},

where x(t)=[𝐯d,𝐯q,𝐢d,𝐢q,𝐢td,𝐢tq](4n+2m)\textbf{x}(t)=[\mathbf{v}_{d}^{\top},\mathbf{v}_{q}^{\top},\mathbf{i}_{d}^{\top},\mathbf{i}_{q}^{\top},\mathbf{i}_{td}^{\top},\mathbf{i}_{tq}^{\top}]^{\top}\in\mathbb{R}^{(4n+2m)} is the state vector, u(t)=[𝐮d,𝐮q]2n\textbf{u}(t)=[\mathbf{u}_{d}^{\top},\mathbf{u}_{q}^{\top}]^{\top}\in\mathbb{R}^{2n} is the input vector and 𝐲(t)=[𝐱d,𝐱q]2n\mathbf{y}(t)=[\mathbf{x}_{d}^{\top},\mathbf{x}_{q}^{\top}]^{\top}\in\mathbb{R}^{2n} is the output vector. A(4n+2m)×(4n+2m)A\in\mathbb{R}^{(4n+2m)\times(4n+2m)} is the state matrix of the system, B(4n+2m)×(2n)B\in\mathbb{R}^{(4n+2m)\times(2n)}, C(2n)×(4n)C\in\mathbb{R}^{(2n)\times(4n)} is the measurement or output matrix and E(4n+2m)×(2n)E\in\mathbb{R}^{(4n+2m)\times(2n)} is the disturbance matrix resulting from the bounded load current demands at the AC buses.

Refer to caption
Figure 2: Droop-based FS-MPC Control Diagram

5 Control Development

Figure 2 shows the control diagram of the i-th IBR supposedly connected to a j-th neighboring inverter through the transmission lines in the set 𝒱×𝒱={1,.,m}\mathcal{E}\in\mathcal{V}\times\mathcal{V}=\{1,....,m\}. The filter output voltage, which is the same as the capacitor voltage 𝐯dq\mathbf{v}_{dq} and the output current of the filter 𝐢fdq\mathbf{i}_{fdq} are used to compute the active and the reactive power references that are fed into the droop controller. The droop control controller provides the frequency and the voltage required by the PI controller in the voltage control loop to compute the current reference 𝐢ref\mathbf{i}_{ref} for the FCS-MPC, which in turn provides the optimum PWM signals SabcS_{abc} to switch the IBRs.

The control objectives in this paper are to ensure perfect tracking of the IBR filter output voltage, AC grid frequency regulation, active and reactive power sharing and the stability analysis for the FCS-MPC. The next subsections detail the droop control development for power sharing and the Finite-set model predictive control formulation for the filter output voltage and current tracking.

5.1 Droop Control development

According to [14], if δij\delta_{ij} is very small for a highly inductive transmission line. That is, sinδijδij\sin\delta_{ij}\approx\delta_{ij} and cosδij1\cos\delta_{ij}\approx 1, where δij\delta_{ij} is the phase shift angle, viv_{i} and vjv_{j} are the output voltages of each inverter and XjX_{j} is the inductive reactance of the transmission line. Then, equation for the output active power (P) and reactive power (Q) developed between an IBRi connected to a neighboring IBRj can be rewritten as

P\displaystyle P =vivjXjδij,\displaystyle=\frac{v_{i}v_{j}}{X_{j}}\delta_{ij}, (8)
Q\displaystyle Q =vj(vivjXj).\displaystyle=v_{j}\bigg{(}\frac{v_{i}-v_{j}}{X_{j}}\bigg{)}.

From (8), the control of the voltage difference 𝐯i𝐯j\mathbf{v}_{i}-\mathbf{v}_{j} depends on the reactive power Q and that of the phase angle which is dynamically equivalent to the grid frequency, depends on the active power P. Hence, droop control for several IBRs is given by

vi\displaystyle{v}_{i} =vnom,inqiQi,\displaystyle={v}_{nom,i}-n_{qi}Q_{i}, (9)
ωi\displaystyle\omega_{i} =ωgmpiPi,\displaystyle=\omega_{g}-m_{pi}P_{i},

where vnom,iv_{nom,i} is the nominal voltage rating of each IBR that must be met for perfect voltage regulation, ωg\omega_{g} is a set of constant dq rotating reference frame frequency, nqin_{qi} is the reactive power-voltage (QV)(Q-V) droop coefficient of each IBR, mpim_{pi} is the active power-frequency (Pω)(P-\omega) droop coefficient of each IBR and

Pi\displaystyle P_{i} =32𝐯dqiT𝐢dqi,\displaystyle=\frac{3}{2}\mathbf{v}_{dqi}^{T}\mathbf{i}_{dqi}, (10)
Qi\displaystyle Q_{i} =32𝐯dqiTJ𝐢dqi.\displaystyle=\frac{3}{2}\mathbf{v}_{dqi}^{T}J\mathbf{i}_{dqi}.

The droop coefficients must be carefully selected or computed because they affect stability and the way loads are being shared between a set of inverter-based resources (IBRs)[12]. Thus, for the total power P to be shared among n inverters, the load distribution between them must satisfy the following:

mp1P1=mp2P2=\displaystyle m_{p1}P_{1}=m_{p2}P_{2}=\ldots =mpnPn,\displaystyle=m_{pn}P_{n}, (11)
P1+P2+P3++Pn\displaystyle P_{1}+P_{2}+P_{3}+\ldots+P_{n} =P.\displaystyle=P.

When the grid synchronization is perfectly matched at a steady state, active power is efficiently shared between IBRi and IBRj as long as the condition in (11) is satisfied.

idiidj=mpjmpii,j𝒱.\frac{\textbf{i}_{di}}{\textbf{i}_{dj}}=\frac{m_{pj}}{m_{pi}}\quad\forall i,j\in\mathcal{V}. (12)

That is using PiP_{i} in equations (9), (10), and assuming perfect synchronization at steady state, ωi=ωg𝟏\omega_{i}=\omega_{g}\mathbf{1}. Thus,

PiPj=32𝐯dqiT𝐢dqi32𝐯dqjT𝐢dqj=𝐯di𝐢di+𝐯qi𝐢qi𝐯dj𝐢dj+𝐯qj𝐢qj=mpjmpii,j𝒱.\frac{P_{i}}{P_{j}}=\frac{\frac{3}{2}\mathbf{v}_{dqi}^{T}\mathbf{i}_{dqi}}{\frac{3}{2}\mathbf{v}_{dqj}^{T}\mathbf{i}_{dqj}}=\frac{\mathbf{v}_{di}\mathbf{i}_{di}+\mathbf{v}_{qi}\mathbf{i}_{qi}}{\mathbf{v}_{dj}\mathbf{i}_{dj}+\mathbf{v}_{qj}\mathbf{i}_{qj}}=\frac{m_{pj}}{m_{pi}}\quad\forall i,j\in\mathcal{V}. (13)

Since voltage does not vary much in a voltage-regulated power grid and the power factor is almost unity such that vd\textbf{v}_{d} component of the voltage is much higher than the vq\textbf{v}_{q} component. Hence,

𝐢di𝐢dj=mpjmpii,j𝒱.\frac{\mathbf{i}_{di}}{\mathbf{i}_{dj}}=\frac{m_{pj}}{m_{pi}}\quad\forall i,j\in\mathcal{V}. (14)
Remark 2

We assume that the AC microgrid is well-regulated, perfectly synchronized, balanced, and has a power factor approximately equal to 1.

5.2 Finite Control Set MPC Formulation

To develop the FCS-MPC algorithm for the state space dynamics of the decentralized system shown in (3). In discrete state space form, the dynamics in (LABEL:ss_cont) then becomes

x(k+1)\displaystyle\textbf{x}(k+1) =Adx(k)+Bdu(k)+EdiLdq(k),\displaystyle=\mathrm{A_{d}}\textbf{x}(k)+\mathrm{B_{d}}\textbf{u}(k)+\mathrm{E_{d}\textbf{i}_{L{dq}}(k)}, (15)
y(k)\displaystyle\mathrm{\textbf{y}(k)} =Cdx(k),\displaystyle=\mathrm{C_{d}}\textbf{x}(k),

where Ad=I+Aτs\mathrm{A_{d}}=I+\mathrm{A}\tau_{s}, Bd=Bτs\mathrm{B_{d}}=\mathrm{B}\tau_{s}, Cd=CC_{d}=C and Ed=Eτs\mathrm{E_{d}}=\mathrm{E}\tau_{s}. The dimensions of the state matrices and the state variables remain the same except for the discrete input matrix BdB_{d} and the control inputs u(k)\textbf{u}(k) in (15). The control input now belongs to the finite set of inputs 𝒰\mathcal{U}. That is:

u(k)𝒰={u1,u2,,uM}2n\textbf{u}(k)\in\mathcal{U}=\left\{u_{1},u_{2},\ldots,u_{M}\right\}\subset\mathbb{R}^{2n} (16)

Since a modulator is not required for PWM generation, the decision variables for FCS-MPC are the switching states (Sa,Sb,ScS_{a},S_{b},S_{c}) in each phase of the 3-phase, 2-level inverters used in this system and the relationship between the switching state vector SabcS_{abc} and the finite number of control variables u(k)\textbf{u}(k) is given as

Sabc\displaystyle S_{abc} =23(Sa+Sbej2π3+Scej4π3)\displaystyle=\frac{2}{3}(S_{a}+S_{b}e^{j\frac{2\pi}{3}}+S_{c}e^{j\frac{4\pi}{3}}) (17)
𝐮i(k)\displaystyle\mathbf{u}_{i}(k) =VdcSabci{1,2,,M}.\displaystyle=V_{dc}S_{abc}\quad\forall i\in\{1,2,\ldots,M\}.

The control objective here is to regulate the voltage vd\textbf{v}_{d} across the filter capacitor and track the injected current provided by the voltage loop using the following optimization problem:

Minimize:𝐯,𝐢\displaystyle\operatorname*{Minimize:}_{\mathbf{v},\mathbf{i}} 𝐲refC𝐱(k+1)2\displaystyle\left\|\mathbf{y}_{\text{ref}}-C\mathbf{x}(k+1)\right\|^{2} (18)
s.t. x(k+1)=Adx(k)+Bdu(k)+EdiLdq(k),\displaystyle\textbf{x}(k+1)=\mathrm{A_{d}}\textbf{x}(k)+\mathrm{B_{d}}\textbf{u}(k)+\mathrm{E_{d}\textbf{i}_{L{dq}}(k)},
𝐮(k)𝒰,\displaystyle\mathbf{u}(k)\in\mathcal{U},

where y(k)=[𝐯d,𝐯q,𝐢d,𝐢q]\textbf{y}(k)=[\mathbf{v}_{d}^{\top},\mathbf{v}_{q}^{\top},\mathbf{i}_{d}^{\top},\mathbf{i}_{q}^{\top}]^{\top} and C2n×4nC\in\mathbb{R}^{{2n}\times{4n}}. The reference currents idqref(k)=[idref(k),iqref(k)]\textbf{i}_{dqref}(k)=[\textbf{i}_{dref}(k)^{\top},\textbf{i}_{qref}(k)^{\top}]^{\top} for the FCS-MPC controller is generated from the filter dynamic by using PI controller in the voltage control loop as shown in (19) [15].

idqref=idqJωgCivdq+Kp(vdqrefvdq)+Ki0τ(vdqrefvdq)𝑑τ\begin{split}\textbf{i}_{dqref}=\textbf{i}_{dq}-J\omega_{g}C_{i}\textbf{v}_{dq}+K_{p}\left(\textbf{v}_{dqref}-\textbf{v}_{dq}\right)+\\ K_{i}\int_{0}^{\tau}\left(\textbf{v}_{dqref}-\textbf{v}_{dq}\right)d\tau\end{split} (19)

Reference voltages vdqref(k)=[vdref(k),vqref(k)]\textbf{v}_{dqref}(k)=[\textbf{v}_{dref}(k)^{\top},\textbf{v}_{qref}(k)^{\top}]^{\top} is the grid rated voltage of 220V220\mathrm{~{}V} with vqref=0\textbf{v}_{qref}=0, as mentioned in remark 2. And PI controller gains KpK_{p} and KiK_{i} are carefully selected. Then, since AC grid reference voltage doesn’t change much at every sampling time and the inverter switching frequency is much higher than the AC grid frequency, we can approximate both the voltage and the current references as vdref(k+1)vdref\textbf{v}_{dref}(k+1)\approx\textbf{v}_{dref} and idref(k+1)idref\textbf{i}_{dref}(k+1)\approx\textbf{i}_{dref}.

6 Numerical Simulation Results

We implement the proposed Droop-based FCS-MPC control algorithm on three IBRs (n=3n=3) with the current flow through two transmission lines (m=2m=2) connecting them as shown in Figure 3, where the incidence matrix ik3×2\mathcal{B}_{ik}\in\mathbb{R}^{3\times 2} is given as:

ik=[101101]\mathcal{B}_{ik}=\left[\begin{array}[]{cc}1&0\\ -1&1\\ 0&-1\end{array}\right]
Refer to caption
Figure 3: Current flow between connected IBRs

In our numerical simulation, we assume a homogeneous configuration of IBRs, and the parameters are given in Table 1. The transmission line impedances are Z1Z_{1} connecting IBR1 to IBR2 and Z2Z_{2} connecting IBR2 to IBR3 as shown in Figure 3. Z1Z_{1} and Z2Z_{2} have resistances rt1=rt2=0.2Ωr_{t1}=r_{t2}=0.2\Omega and inductances lt1=lt2=0.3mHl_{t1}=l_{t2}=0.3\mathrm{~{}mH}. Constant impedance loads are available at buses 1, 2 and 3 as SL1=SL2=1+j0.2kVAS_{L1}=S_{L2}=1+j0.2\mathrm{~{}kVA}, SL3=3+j0.2kVAS_{L3}=3+j0.2\mathrm{~{}kVA}. Additional load of SL4=2+j0.2kVAS_{L4}=2+j0.2\mathrm{~{}kVA} is suddenly added to bus 3 at t=1st=1s and later disconnected at t=3st=3s.

Table 1: parameters for the IBRs
Parameters Value
DC Source Voltage,VdcV_{dc} 600V600\mathrm{~{}V}
Rated Voltage (rms) 220V220\mathrm{~{}V}
Grid frequency, ωg/(2π)\omega_{g}/(2\pi) 50Hz50\mathrm{~{}Hz}
Switching frequency, fsf_{s} 10kHz10\mathrm{~{}kHz}
Filter inductor, lil_{i} 3.5mH3.5\mathrm{~{}mH}
Filter capacitor, CiC_{i} 50μF50\mathrm{~{}\mu F}
Filter resistor, rir_{i} 0.2Ω0.2\mathrm{~{}\Omega}
P-f droop coefficients, mp1=mp2=mp3m_{p1}=m_{p2}=m_{p3} 6e5rad/s/W6e^{-5}\mathrm{~{}rad/s/W}
Q-V droop coefficients, nq1=nq2=nq3n_{q1}=n_{q2}=n_{q3} 4e3V/Var4e^{-3}\mathrm{~{}V/Var}
Refer to caption
Figure 4: (a) Response of 3-phase currents to load changes and (b) Filter capacitor voltage tracking for IBR1, IBR2 and IBR3.
Refer to caption
Figure 5: FCS-MPC-based Droop for IBR1, IBR2 and IBR3: (a)Active-Power/Frequency (b) Reactive-power/Voltage

Figure 4(a) shows that 3-phase currents in each IBR change when the load changes. That is, the currents increase when load is added at bus 3, and return to its initial value when the load is disconnected. Figure 4(b) shows that the output voltage across the capacitor does not change when the load is disconnected and reconnected. The voltage remains at a peak value of 311V(220Vrms)311\mathrm{V}(220\mathrm{V_{rms}}) regardless of the change in load. This confirms that the FCS-MPC provides effective voltage regulation across the filter capacitor at the AC buses.

Figure 5(a) shows that efficient frequency synchronization and active power sharing with fewer oscillations are achieved at a steady state. Frequency deviation is minimal and reaches a steady state quickly during connection and disconnection of load at bus 3. Active power is equally shared proportional to the stiff frequency maintained in the AC grid. This shows the effectiveness of the decentralized droop-based FCS-MPC control for active power sharing. Slight frequency deviation is normal based on (9) but this can be further improved by a higher-level controller which is not considered here. Figure 5(b) shows the reactive power sharing for the output voltages of IBR1, IBR2, and IBR3. Reactive power changes with respect to changes in load and this deviation is more pronounced at bus 3 where additional load is added. Although the inverter voltages do not deviate significantly from each but the reactive is more pronounced at higher load changes. Droop-based FCS-MPC reduces chattering and oscillations in reactive power sharing and the slight deviation is due to the grid impedance mismatch [16]. And has been addressed in literature by using virtual impedance droop control [4, 17]. Unlike the P-f droop, the Q-V droop is local to each inverter depending on how long they are connected to each other with respect to line impedance variations. However, the reactive powers and the voltages of the IBRs reach a steady state quickly as the load connects and disconnects at bus 3 efficiently. This shows the effectiveness of FCS-MPC in the inner control loop.

7 Conclusion

A Droop-based finite set model predictive control technique is proposed to regulate the voltage across the filter capacitor at the load bus of an islanded AC microgrid. This method provides a PWM signal to control the inverters directly without the need for a modulator. Accurate frequency synchronization and effective active power sharing are achieved at a steady state with fewer oscillations. Filter capacitor voltage at the AC bus remains unchanged despite the load changes as depicted in the numerical simulations thereby confirming the affirmation in Figure 4(b). Reactive power sharing is also achieved with fewer oscillations and low inverter voltage variations. This method eliminates the need for tuning parameter gains in the inner current loop controller, making it faster and the slight frequency deviation at steady state can be further improved by introducing a secondary or supervisory controller in future work. Also, we developed the condition for bounded stability for FCS-MPC.

References

  • [1] N. Hatziargyriou, H. Asano, R. Iravani, and C. Marnay, “Microgrids,” IEEE Power and Energy Magazine, vol. 5, no. 4, pp. 78–94, 2007.
  • [2] J. M. Guerrero, J. C. Vasquez, J. Matas, L. G. de Vicuna, and M. Castilla, “Hierarchical control of droop-controlled ac and dc microgrids—a general approach toward standardization,” IEEE Transactions on Industrial Electronics, vol. 58, no. 1, pp. 158–172, 2011.
  • [3] F. Milano, F. Dörfler, G. Hug, D. J. Hill, and G. Verbiv̈c, “Foundations and challenges of low-inertia systems (invited paper),” 2018 Power Systems Computation Conference (PSCC), pp. 1–25, 2018.
  • [4] O. M. Anubi and S. Ameli, “Robust stabilization of inverter-based resources using virtual resistance-based control,” IEEE Control Systems Letters, vol. 6, pp. 3295–3300, 2022.
  • [5] S. T. Cady, M. Zholbaryssov, A. D. Domínguez-García, and C. N. Hadjicostis, “A distributed frequency regulation architecture for islanded inertialess ac microgrids,” IEEE Transactions on Control Systems Technology, vol. 25, no. 6, pp. 1961–1977, 2017.
  • [6] G. F. Silva, A. Donaire, M. M. Seron, A. McFadyen, and J. Ford, “String stability in microgrids using frequency controlled inverter chains,” IEEE Control Systems Letters, vol. 6, pp. 1484–1489, 2021.
  • [7] A. Khaledian and M. A. Golkar, “Analysis of droop control method in an autonomous microgrid,” Journal of applied research and technology, vol. 15, no. 4, pp. 371–377, 2017.
  • [8] J. M. Guerrero, M. Chandorkar, T.-L. Lee, and P. C. Loh, “Advanced control architectures for intelligent microgrids—part i: Decentralized and hierarchical control,” IEEE Transactions on Industrial Electronics, vol. 60, no. 4, pp. 1254–1262, 2013.
  • [9] X. Lu, J. Wang, J. M. Guerrero, and D. Zhao, “Virtual-impedance-based fault current limiters for inverter dominated ac microgrids,” IEEE Transactions on Smart Grid, vol. 9, no. 3, pp. 1599–1612, 2018.
  • [10] Y. Dong, B. Ren, and Q.-C. Zhong, “Bounded universal droop control to enable the operation of power inverters under some abnormal conditions and maintain voltage and frequency within predetermined ranges,” IEEE Transactions on Industrial Electronics, vol. 69, no. 11, pp. 11 633–11 643, 2022.
  • [11] J. He, Y. W. Li, J. M. Guerrero, J. C. Vasquez, and F. Blaabjerg, “An islanding microgrid reactive power sharing scheme enhanced by programmed virtual impedances,” in 2012 3rd IEEE International Symposium on Power Electronics for Distributed Generation Systems (PEDG), 2012, pp. 229–235.
  • [12] S. Bayhan, H. Abu-Rub, and R. S. Balog, “Model predictive control of quasi-z-source four-leg inverter,” IEEE Transactions on Industrial Electronics, vol. 63, no. 7, pp. 4506–4516, 2016.
  • [13] J. Watson, Y. Ojo, I. Lestas, and C. Spanias, “Stability of power networks with grid-forming converters,” in 2019 IEEE Milan PowerTech, 2019, pp. 1–6.
  • [14] H. Laaksonen, P. Saari, and R. Komulainen, “Voltage and frequency control of inverter based weak lv network microgrid,” in 2005 International Conference on Future Power Systems.   IEEE, 2005, pp. 6–pp.
  • [15] L. Subramanian and H. B Gooi, “Stabilizing droop variation of converter-connected generation in autonomous microgrids with virtual inertia control,” in 2018 IEEE Energy Conversion Congress and Exposition (ECCE), 2018, pp. 3746–3753.
  • [16] Q.-C. Zhong, “Robust droop controller for accurate proportional load sharing among inverters operated in parallel,” IEEE Transactions on industrial Electronics, vol. 60, no. 4, pp. 1281–1290, 2011.
  • [17] J. He and Y. W. Li, “Analysis, design, and implementation of virtual impedance for power electronics interfaced distributed generation,” IEEE Transactions on industry Applications, vol. 47, no. 6, pp. 2525–2538, 2011.