**Distributed Energy Resources to Improve the Power Quality and to Reduce Energy Costs of a Hybrid AC/DC Microgrid**

Luisa Alfieri, Antonio Bracale, Pierluigi Caramia and Guido Carpinelli

Additional information is available at the end of the chapter

http://dx.doi.org/10.5772/intechopen.68766

#### Abstract

This chapter deals with microgrids (μGs), i.e., a group of interconnected loads and distributed energy resources that act as a single controllable entity with respect to the grid. The μGs can be classified into AC and DC μGs depending on the characteristics of the supply voltage, with both solutions characterized by advantages and challenges. Recently, hybrid AC/DC μGs have been developed with the aim to exploit the advantages of both AC and DC solutions. Hybrid μGs require being properly controlled to guarantee their optimal behavior, in both grid-connected and islanding mode. In this chapter, we propose an optimal control strategy for a hybrid μG to be realized in an actual Italian industrial facility. The control strategy operates with the aim to simultaneously minimize the energy costs and to compensate waveform distortions. The key result of the chapter consists in evidencing the technical and economic advantages of the proposed solution by means of real-time simulations of the hybrid μG performed through Matlab/Simulink development tool in the different conditions (grid-connected and islanding mode).

Keywords: hybrid microgrids, distributed energy resources, power quality, optimal control strategies

### 1. Introduction

Nowadays distributed generators, storage systems, and controllable loads much more actively and simultaneously contribute to an optimal operation of the electrical distribution systems in the frame of the new concepts of smart grids (SGs) and microgrids (μGs).

In this exciting context, new technologies and services are being introduced to make the electrical networks more reliable, efficient, secure, and environmentally friendly.

This chapter deals with μGs that CIGRÉ C6.22 Working Group defines as electricity distribution systems containing loads and distributed energy resources that can be operated in a controlled, coordinated way either while they are connected to the main power network (grid-connected mode), or when operate in islanding mode [1].

As well known, μGs can be AC and DC μGs on the basis of the nature of the supplied voltage [2–4]. AC μGs have the advantage of utilizing existing AC technologies as well as DC μGs seem particularly suitable for supporting the current needs because most distributed generation sources, such as photovoltaic plants, fuel cells, and storage systems, generate DC power directly. Moreover, DC grids can also guarantee high power quality levels to AC sensible loads, even though this provision requires the presence of additional DC/AC converters.

Recently, hybrid AC/DC μGs have been developed with the aim to exploit the advantages of both AC and DC solutions and, in particular, to accelerate the integration process of DC power technologies into the existing consolidated AC systems [5–10].

As well known, several static converters are installed in the hybrid μGs, either to connect renewable generation and storage systems to the AC or DC grids, or as an interface between AC and DC grids. In this context, optimization strategies are required to perform the converters optimal control with the aim of efficiently operating the whole hybrid μG. These control strategies should take into account all goals and operating constraints of both DC and AC sections.

Motivated by the above requirements, in the following sections of this chapter, we propose an optimal control strategy for a hybrid μG to be realized in an actual Italian industrial facility. In particular, the static converters of the hybrid μG are controlled with the aim to simultaneously minimize the energy costs and to compensate waveform distortions, thus ensuring an optimal technical and economical behavior of the whole μG.

The remainder of the chapter is organized as follows. The electrical system of the actual Italian industrial facility is briefly recalled in Section 2, where also the structure of the proposed hybrid microgrid is shown. Section 3 deals with the proposed control strategy. Numerical applications are shown in Section 4, Section 5 provides our conclusions, and some data of the test system under study are explained in Appendix.

### 2. The electrical industrial system under study

The electrical industrial system under study was a hybrid AC/DC microgrid obtained by the modification of the existing LV distribution system of an actual Italian industrial facility where transformers are assembled. The simplified electrical scheme before the modification is reported in Figure 1. The existing facility's electrical distribution system is connected to the MV grid through a 20/0.4–630 kVA transformer and includes four low-voltage feeders; each feeder is dedicated to a different manufacturing process, i.e.: (i) tanks and boxes manufactory, (ii) assembly, (iii) winding and coils and (iv) test. An automatically switched 120 kVAr capacitor bank with 10 kVAr step is installed at bus 3 to guarantee a power factor at PCC (bus 2) equal to 0.95.

Distributed Energy Resources to Improve the Power Quality and to Reduce Energy Costs of a Hybrid AC/DC… http://dx.doi.org/10.5772/intechopen.68766 23

Figure 1. Electrical simplified scheme of the existing electrical distribution system of the industrial plant under study.

The data of the lines, transformers and loads are reported in Appendix.

This chapter deals with μGs that CIGRÉ C6.22 Working Group defines as electricity distribution systems containing loads and distributed energy resources that can be operated in a controlled, coordinated way either while they are connected to the main power network

As well known, μGs can be AC and DC μGs on the basis of the nature of the supplied voltage [2–4]. AC μGs have the advantage of utilizing existing AC technologies as well as DC μGs seem particularly suitable for supporting the current needs because most distributed generation sources, such as photovoltaic plants, fuel cells, and storage systems, generate DC power directly. Moreover, DC grids can also guarantee high power quality levels to AC sensible loads, even though this provision requires the presence of additional DC/AC

Recently, hybrid AC/DC μGs have been developed with the aim to exploit the advantages of both AC and DC solutions and, in particular, to accelerate the integration process of DC power

As well known, several static converters are installed in the hybrid μGs, either to connect renewable generation and storage systems to the AC or DC grids, or as an interface between AC and DC grids. In this context, optimization strategies are required to perform the converters optimal control with the aim of efficiently operating the whole hybrid μG. These control strategies should take into account all goals and operating constraints of both DC and AC sections.

Motivated by the above requirements, in the following sections of this chapter, we propose an optimal control strategy for a hybrid μG to be realized in an actual Italian industrial facility. In particular, the static converters of the hybrid μG are controlled with the aim to simultaneously minimize the energy costs and to compensate waveform distortions, thus ensuring an optimal

The remainder of the chapter is organized as follows. The electrical system of the actual Italian industrial facility is briefly recalled in Section 2, where also the structure of the proposed hybrid microgrid is shown. Section 3 deals with the proposed control strategy. Numerical applications are shown in Section 4, Section 5 provides our conclusions, and some data of the

The electrical industrial system under study was a hybrid AC/DC microgrid obtained by the modification of the existing LV distribution system of an actual Italian industrial facility where transformers are assembled. The simplified electrical scheme before the modification is reported in Figure 1. The existing facility's electrical distribution system is connected to the MV grid through a 20/0.4–630 kVA transformer and includes four low-voltage feeders; each feeder is dedicated to a different manufacturing process, i.e.: (i) tanks and boxes manufactory, (ii) assembly, (iii) winding and coils and (iv) test. An automatically switched 120 kVAr capacitor bank with 10 kVAr step is installed at bus 3 to guarantee a power factor at PCC (bus 2)

(grid-connected mode), or when operate in islanding mode [1].

technologies into the existing consolidated AC systems [5–10].

technical and economical behavior of the whole μG.

test system under study are explained in Appendix.

2. The electrical industrial system under study

converters.

22 Development and Integration of Microgrids

equal to 0.95.

The modifications on the existing electrical system of Figure 1 and the applied control strategy implemented in the centralized control system (CCS) are finalized (i) to improve the continuity and quality of the energy delivered to the loads and (ii) to optimize the exchange of energy between industrial plant and MV distribution grid.

The following modifications are effected (Figure 2): (i) a DC micro grid (DCμG) is connected at AC bus 16 through an AC/DC grids interfacing converter; (ii) a Battery Energy Storage System (BESS) is installed at AC bus 3 through a DC/AC static converter; and (iii) a dispatchable micro-turbine is connected at AC bus 20 through an AC/AC static converter.<sup>1</sup>

The DC micro grid (DCμG) consists of a photovoltaic (PV) plant equipped with a maximum power point tracker control system and connected to the DC grid through a DC/DC static converter and three sensitive AC loads (Folding walls island robot, Sandblasting machine and PLCs used to control the automation of the whole electromechanical process<sup>2</sup> ) moved from the original electrical scheme of Figure 1.

The sensitive loads are controllable in terms of start-up times; in particular, their working time intervals are obtained by a day-ahead scheduling aimed at achieving the minimization of the costs for the electricity purchase of the hybrid AC/DC μG [12].<sup>3</sup>

<sup>1</sup> In the following we will use the term micro-turbine, BESS and PV generators to indicate the system constituted by the micro-turbine, BESS and PV panels together with their interfacing static converters.

<sup>2</sup> These loads were chosen in accordance with the technical operators of the industrial facility on the basis of a deep study on the industrial process [11]

<sup>3</sup> An Uninterruptible Power Supply (UPS) system is dedicated for PLC power backup to provide a soft shut-down of the load control devices in case AC μG out of service.

Figure 2. Scheme of the hybrid AC/DC µG.

The above modifications lead to a hybrid μG, which includes AC buses (ACμG) and DC buses (DCμG).

The data of the newly added components are reported in Appendix.

### 3. The proposed control strategy

The proposed control strategy is based on the solution of optimization problems performed at fundamental and at harmonic frequencies, where the efficient operation of the DC and AC grids is guaranteed by a proper control of the static converters installed in the hybrid μG, either to connect micro-turbine and storage system to the AC grid, or as interface between AC and DC grids.

The control strategy allows the hybrid AC/DC μG to operate both in grid-connected and in islanding mode as described in details in the following subsections.

### 3.1. Grid-connected mode

When the hybrid μG is in grid-connected mode, the converters are controlled with the aim to simultaneously minimize: (i) the energy costs sustained for the electrical energy imported from the main distribution network and the production cost of the micro-turbine and (ii) the waveform distortions of the AC bus voltages.

The converter sizes have to be appropriately chosen to correctly perform both the required services.

Two different optimization models are formulated and solved by CCS to provide the reference signals for the energy cost minimization (reference signals at fundamental frequency) and for the waveform distortion compensation (reference signals at harmonic frequencies), as shown in Figure 3.

Figure 3. Optimal control strategy for grid-connected mode.

The above modifications lead to a hybrid μG, which includes AC buses (ACμG) and DC buses

The data of the newly added components are reported in Appendix.

(DCμG).

Figure 2. Scheme of the hybrid AC/DC µG.

24 Development and Integration of Microgrids

#### 3.1.1. Energy cost minimization

At the fundamental frequency, the real-time operation of the hybrid AC/DC µG is optimized by controlling the active and reactive power of the BESS, the active and reactive power of the micro-turbine and the reactive power provided by the AC/DC grids interfacing converter. The reference signals for the converters, which minimize the total energy cost, are obtained solving a non-linear constrained optimization problem, such as:

$$\min f\_{obj}(\mathbf{x})\tag{1}$$

s.t.

$$\beta\_h(\mathbf{x}) = 0 \; h = 1, \ldots, N\_{eq} \tag{2}$$

$$\text{type}\_k(\mathbf{x}) \le 0 \; k = 1, \dots, N\_{\text{ineq}} \tag{3}$$

where fobj is the objective function and β<sup>h</sup> and γ<sup>k</sup> are the hth equality and the kth inequality constraints to be met, respectively. The vector x includes the state and control variables.

Dividing the day into Nt time slots of the same duration Δt, the optimization problem in Eqs. (1)–(3) is repeatedly solved at each ith time slot and its solution furnishes the reference signals needed at the successive time slot (i + 1)th for the controllers of micro-turbine, BESS and AC/DC grids interfacing converters.

Input and output data, the objective function and the constraints structure are listed in the following.<sup>4</sup>

The input and output data of the optimization problem are:

Input data:


<sup>4</sup> It should be noted that a fixed value is assigned to the energy stored in the battery at the last time slot of the day (Nt), thus influencing the battery power profile over the whole day; then, the optimization problem becomes multi-period with the consequence that, at the i th time slot, objective function and constraints must be formulated along all the rest of the day (i.e., time slots included in {i+1,…, Nt}).

<sup>5</sup> The real time control procedure follows a day ahead scheduling procedure (performed off-line), as the one proposed in [12], finalized to schedule the powers of DCµG controllable loads.

• the charge and discharge periods of the BESS in ACµG which are fixed a-priori in function of the energy hourly tariff (typically, the discharge period is during the hours of peak energy tariff).

Output data:

3.1.1. Energy cost minimization

26 Development and Integration of Microgrids

AC/DC grids interfacing converters.

the successive day;

• the energy charge;

the consequence that, at the i

(i.e., time slots included in {i+1,…, Nt}).

4

5

s.t.

following.<sup>4</sup>

Input data:

a non-linear constrained optimization problem, such as:

The input and output data of the optimization problem are:

• the power required by the controllable loads5

in [12], finalized to schedule the powers of DCµG controllable loads.

At the fundamental frequency, the real-time operation of the hybrid AC/DC µG is optimized by controlling the active and reactive power of the BESS, the active and reactive power of the micro-turbine and the reactive power provided by the AC/DC grids interfacing converter. The reference signals for the converters, which minimize the total energy cost, are obtained solving

where fobj is the objective function and β<sup>h</sup> and γ<sup>k</sup> are the hth equality and the kth inequality constraints to be met, respectively. The vector x includes the state and control variables.

Dividing the day into Nt time slots of the same duration Δt, the optimization problem in Eqs. (1)–(3) is repeatedly solved at each ith time slot and its solution furnishes the reference signals needed at the successive time slot (i + 1)th for the controllers of micro-turbine, BESS and

Input and output data, the objective function and the constraints structure are listed in the

• the state of charge of the BESS at the end of the ith time slot and at the end of the day (i.e., at the end of time slot Nt). In particular, the value at the end of the ith time slot is obtained as result of the optimization problem in Eqs. (1)–(3) solved at the previous (i-1)th time slot; the value at the end of the day is fixed on the basis of the desired value at the beginning of

• the forecasted powers of the PV system and of the non-controllable loads, all provided by

It should be noted that a fixed value is assigned to the energy stored in the battery at the last time slot of the day (Nt), thus influencing the battery power profile over the whole day; then, the optimization problem becomes multi-period with

The real time control procedure follows a day ahead scheduling procedure (performed off-line), as the one proposed

;

th time slot, objective function and constraints must be formulated along all the rest of the day

proper forecasting tools for all the time slots from (i + 1)th to Nt [13–15];

min f objðxÞ ð1Þ

βhðxÞ ¼ 0 h ¼ 1, …, Neq ð2Þ

γkðxÞ ≤ 0 k ¼ 1, …, Nineq ð3Þ


All the output data are calculated from the (i + 1)th time slot to the end of the day, but, obviously, only the values of the (i + 1)th time slot are sent to the controllers of the ACµG BESS, micro-turbine and AC/DC grids interfacing converter.

The objective function in Eq. (1) to be minimized is the daily total costs sustained by the industrial facility for the electrical energy. This cost is the sum of the cost for the energy imported from the MV electrical distribution grid and of the production cost of the microturbine. The objective function fobj(x) results in:

$$f\_{\text{obj}}(\mathbf{x}) = \sum\_{k=i+1}^{N\_l} (\mathbf{Pr}\_{\text{E},k} \, P\_k^{grid} \Delta t + \mathbf{C}\_k^{WTG} \Delta t) \tag{4}$$

In Eq. (4), Pgrid <sup>k</sup> is the active power furnished by the distribution grid at the kth time interval, PrE,k and CmTG <sup>k</sup> are the energy charge and the production cost of the micro-turbine generator at the kth time interval, expressed as [16]:

$$\mathbf{C}\_{k}^{\rm mTG} = a\_{\rm mTG} \cdot (P\_{20,k}^{\rm mTG})^2 + b\_{\rm mTG,n} P\_{20,k}^{\rm mTG} + c\_{\rm mTG} \,\mu(k), k = \, i+1, \ldots, \, \text{Nt} \tag{5}$$

In Eq. (5), amTG, bmTG and cmTG are the cost coefficients of the micro-turbine, PmTG <sup>20</sup>, <sup>k</sup> is the power injected by the micro-turbine connected at busbar 20, u(k) is equal to 1 (if the micro-turbine is switched ON) or 0 (if the micro-turbine is switched OFF). The start-up cost is not considered in Eq. (5) [17].

The equality and inequality constraints related to the DCµG, the ACµG, and the AC/DC grids interfacing converter are illustrated in Tables 1–3, respectively.

With reference to the DCµG constituted by MDC buses (Table 1), the equality constraints are the DC load flow Eq. (6), while the inequality constraints refer to limits on bus voltage magnitudes Eq. (7) and line currents Eq. (8).

In Eq. (6), PDC s, <sup>k</sup> is the active power injected in the sth DC bus at the kth time slot, VDC s, <sup>k</sup> is the sth bus voltage magnitude in the kth time slot, Gsd is the s-d element of the conductance matrix.6

<sup>6</sup> The terms of the conductance matrix don't depend on the time interval, since no time variation of the DCµG topology is assumed.

#### Load flow equations

$$P\_{s,k}^{\rm DC} = \sum\_{d=1}^{M\_{\rm DC}} V\_{s,k}^{\rm DC} V\_{d,k}^{\rm DC} \mathbf{G}\_{sd} \quad , \quad k = 1 + 1, \ldots, N\_{\rm t}, \mathbf{s} = 1, \ldots, M\_{\rm DC} \tag{6}$$

Voltage magnitudes

$$V\_{s\_{\rm LB}}^{D\complement} \le V\_{s,k}^{D\complement} \le V\_{s\_{\rm LB}}^{D\complement} \quad k = i + 1, \ldots, N\_t,\\ \mathbf{s} = \mathbf{2}, \ldots, M\_{\rm DC} \tag{7}$$

Line currents

$$I\_{l,k}^{\rm DC} \le I\_{l,\max'}^{\rm DC} \quad k = l+1,\ldots,N\_{t\prime} \text{ } l = 1,\ldots,L\_{\rm DC} \tag{8}$$

Table 1. DCµG's constraints.

In particular, PDC s, <sup>k</sup> with s = 2,…, MDC is the scheduled active power of the controllable loads at load busbars and the active power injected by PV generator. In particular, at busbars 2, 3 and 4 only loads are connected and their active powers are the ones reported in Table A. VI of Appendix (with minus sign) corresponding to the busbars 17, 9 and 22 of the actual industrial system shown in Figure 1; at bus 5, only PV generator is connected and the active power is equal to the value obtained by forecasted PV solar irradiance; at bus 1, the voltage VDC <sup>1</sup>, <sup>k</sup> <sup>¼</sup> <sup>V</sup>DC sp , for each kth time slot, since AC/DC interfacing converter controls the DC voltage at the constant value VDC sp .

In Eq. (7), VDC sLB and <sup>V</sup>DC sUB are the lower and upper bounds of voltage amplitude of the sth DC bus. In Eq. (8), I DC l, <sup>k</sup> is the value of current flowing in the lth line during the kth time slot, LDC is the number of the DC grid lines and I DC l,max is the ampacity of the lth line.

With reference to the ACµG constituted by MAC buses (Table 2), the equality constraints are the AC load flow in Eq. (9) and the balance equation of the energy stored in the BESS at the end of the day (Eq. (16)), while the inequality constraints refer to limits on bus voltage magnitudes (Eq. (10)), line currents (Eq. (11)), active and reactive power at PCC (Eq. (12)), MV/LV transformers apparent power (Eq. (13)), micro-turbine active and apparent power (Eqs. (14) and (15)), energy stored, charging/discharging powers and apparent power of BESS in the generic time slot (Eqs. (17)–(19)).

In Eq. (9), PAC s, <sup>k</sup> <sup>ð</sup>QAC s, <sup>k</sup> Þ is the active (reactive) power injected in busbar s at the kth time slot; Vs,k ϑs,k are the magnitude and argument of the voltage ofsth AC bus at the kth time slot, respectively, Gs,j (Bs,j) is the s-j term of the conductance (susceptance) matrix.7 The first busbar of the grid in

<sup>7</sup> The terms of the conductance and susceptance matrix don't depend on the time interval, since no time variation of the AC microgrid topology is assumed.

#### AC grid

#### Load flow equations

$$\begin{aligned} P\_{s,k}^{AC} &= \left. V\_{s,k}^{AC} \sum\_{j=1}^{M\_{AC}} V\_{j,k}^{AC} [\mathbf{G}\_{s,j} \cos \left( \mathfrak{G}\_{s,k} - \mathfrak{G}\_{j,k} \right) + B\_{s,j} \sin(\mathfrak{G}\_{s,k} - \mathfrak{G}\_{j,k})] \right| \\\\ Q\_{s,k}^{AC} &= \left. V\_{s,k}^{AC} \sum\_{j=1}^{M\_{AC}} V\_{j,k}^{AC} [\mathbf{G}\_{s,j} \sin \left( \mathfrak{G}\_{s,k} - \mathfrak{G}\_{j,k} \right) - B\_{s,j} \cos(\mathfrak{G}\_{s,k} - \mathfrak{G}\_{j,k})] \right| \\ k &= i + 1, \dots, N\_t \, \mathbf{s} = 1, \dots, M\_{AC} \end{aligned} \tag{9}$$

Voltage magnitudes

$$V\_{s\_{\rm LB}}^{\rm AC} \le V\_{s,k}^{\rm AC} \le V\_{s\_{\rm LB}}^{\rm AC}, \quad k = \text{ i} + 1, \ldots, N\_t \text{ s} = 2, \text{ ...} \quad M\_{\rm AC} \tag{10}$$

Line currents

$$I\_{l,k}^{AC} \le I\_{l,\max}^{AC}, \quad k = \text{ $i$ } + 1, \dots, \text{N}\_t \text{ } l = 1, \dots, L\_{AC} \tag{11}$$

PPCC powers

In particular, PDC

Table 1. DCµG's constraints.

Load flow equations

28 Development and Integration of Microgrids

Voltage magnitudes

Line currents

PDC s, <sup>k</sup> <sup>¼</sup> <sup>X</sup> MDC

d¼1

VDC sLB <sup>≤</sup>VDC

> I DC l, <sup>k</sup> ≤ I DC

VDC s, <sup>k</sup> VDC

s, <sup>k</sup> ≤VDC

constant value VDC

In Eq. (7), VDC

bus. In Eq. (8), I

time slot (Eqs. (17)–(19)).

s, <sup>k</sup> <sup>ð</sup>QAC

AC microgrid topology is assumed.

In Eq. (9), PAC

7

sp .

sLB and <sup>V</sup>DC

the number of the DC grid lines and I

DC

s, <sup>k</sup> with s = 2,…, MDC is the scheduled active power of the controllable loads at

d, <sup>k</sup>Gsd , k ¼ i þ 1, …, Nt, s ¼ 1, …, MDC ð6Þ

sUB , k ¼ i þ 1, …, Nt, s ¼ 2, …, MDC ð7Þ

l,max, k ¼ i þ 1, …, Nt, l ¼ 1, …, LDC ð8Þ

sUB are the lower and upper bounds of voltage amplitude of the sth DC

l, <sup>k</sup> is the value of current flowing in the lth line during the kth time slot, LDC is

l,max is the ampacity of the lth line.

s, <sup>k</sup> Þ is the active (reactive) power injected in busbar s at the kth time slot; Vs,k ϑs,k

<sup>1</sup>, <sup>k</sup> <sup>¼</sup> <sup>V</sup>DC sp ,

load busbars and the active power injected by PV generator. In particular, at busbars 2, 3 and 4 only loads are connected and their active powers are the ones reported in Table A. VI of Appendix (with minus sign) corresponding to the busbars 17, 9 and 22 of the actual industrial system shown in Figure 1; at bus 5, only PV generator is connected and the active power is

for each kth time slot, since AC/DC interfacing converter controls the DC voltage at the

With reference to the ACµG constituted by MAC buses (Table 2), the equality constraints are the AC load flow in Eq. (9) and the balance equation of the energy stored in the BESS at the end of the day (Eq. (16)), while the inequality constraints refer to limits on bus voltage magnitudes (Eq. (10)), line currents (Eq. (11)), active and reactive power at PCC (Eq. (12)), MV/LV transformers apparent power (Eq. (13)), micro-turbine active and apparent power (Eqs. (14) and (15)), energy stored, charging/discharging powers and apparent power of BESS in the generic

are the magnitude and argument of the voltage ofsth AC bus at the kth time slot, respectively, Gs,j (Bs,j) is the s-j term of the conductance (susceptance) matrix.7 The first busbar of the grid in

The terms of the conductance and susceptance matrix don't depend on the time interval, since no time variation of the

equal to the value obtained by forecasted PV solar irradiance; at bus 1, the voltage VDC

DC

$$P\_{2-3,k}^{AC} \leq P\_{PCC}^{Max} \tag{12}$$

$$Q\_{2-3,k}^{AC} \ge 0 \quad \quad k = i + 1, \dots, N\_t$$

MV/LV transformer apparent power

$$\sqrt{P\_{\text{tr},k}\,^2 + Q\_{\text{tr},k}} \le \text{S}\_{\text{TR}}\,\, \, \, \, k = \, i+1, \dots, N\_t \tag{13}$$

Micro-turbine

Active and apparent power

$$P\_{\min}^{mTG} \text{ } \mu(k) \le P\_{20,\ k}^{mTG} \le P\_{\max}^{mTG} \text{ } \mu(k), \quad k = 1, \pm 1, \ldots, N\_t \tag{14}$$

$$\sqrt{P\_{20,\,k}^{mTG2} + Q\_{20,\,k}^{mTG2}} \le S\_{nom}^{mTG}, \quad k = 1, \ldots, N\_t \tag{15}$$

BESS

Energy stored at the end of the day

$$E\_i - \Delta t \sum\_{k=i+1}^{N\_l} \gamma\_k \, P\_{3,k}^{\text{BESS}} = E\_{sp}^{in}, \quad \gamma\_k = \begin{cases} \frac{1}{\eta\_{\text{lens}}} & \text{if } P\_{3,k}^{\text{BESS}} \ge 0 \\\\ \eta\_{\text{lens}} & \text{if } P\_{3,k}^{\text{BESS}} < 0 \end{cases} \tag{16}$$

#### AC grid

Energy stored in the generic time slot

$$E\_{sp}^{\rm lb} \le E\_{k-1} - \Delta t \ P\_{3,k}^{\rm BESS} \le E\_{sp}^{\rm ub}, \quad k = \text{ i} + 1, \ldots, \left(\text{N}\_{\rm t} - 1\right). \tag{17}$$

Charging and discharging power in the generic time slot

$$\begin{aligned} 0 - P\_{\text{max}}^{ch} \le P\_{3,k}^{BESS} \le 0 \quad k \in \varOmega\_{ch,k} \\ 0 \le P\_{3,k}^{BESS} \le P\_{\text{max}}^{dch} \qquad k \in \varOmega\_{dch,k} \end{aligned} \tag{18}$$

Apparent power

$$\sqrt{P\_{3,k}^{BESS} + Q\_{3,k}^{BESS}} \le S\_{\text{nom}}^{BESS} \quad k = \text{ i} + 1, \ldots, N\_t \tag{19}$$

Table 2. ACµG's constraints.

Figure 1 is the slack busbar. Then, VAC <sup>1</sup>, <sup>k</sup> <sup>¼</sup> <sup>V</sup>AC sp and ϑ1,<sup>k</sup> = 0, for each kth time slot, where VAC sp is the fixed slack voltage magnitude. In particular, with reference to each kth time slot:


In Eq. (10), VAC sLB and <sup>V</sup>AC sUB are the values of lower and upper bounds of voltage at the sth AC bus, respectively.

In Eq. (11), I AC l, <sup>k</sup> is the value of current flowing in the lth line at the kth time interval, LAC is the number of lines of the AC grid and I AC l,max is the ampacity of the lth line.

In Eq. (12), the active power PAC <sup>2</sup>�3, <sup>k</sup> at the PCC busbar during the <sup>k</sup>th time interval (coincident with Pgrid <sup>k</sup> in Eq. (4)) is limited by contractual active power PMax PCC and the reactive power QAC 2�3, k

at the PCC busbar during the kth time interval has to be positive since reactive power cannot be injected into distribution network. The power flows PAC <sup>2</sup>�3, <sup>k</sup> and <sup>Q</sup>AC <sup>2</sup>�3, <sup>k</sup> can be expressed in a simple way using state variables.

In Eq. (13), Ptr,<sup>k</sup> and Qtr,k are the values of active and reactive power flowing through the transformer at the kth time interval which can be expressed in function of state variables; STR is the size of transformer (630 kVA).

In Eq. (14), the active power produced by micro-turbine <sup>ð</sup>PmTG <sup>20</sup>, <sup>k</sup> Þ is within a normal operating interval [PmTG min , PmTG max ].

In Eq. (15), QmTG <sup>20</sup>, <sup>k</sup> is the reactive power injected by converter of the micro-turbine and SmTG nom is its nominal apparent power.

In Eq. (16), Ei is the energy stored in the BESS at the ith time interval, PBESS <sup>3</sup>, <sup>k</sup> is the BESS power at the kth time slot, Ein sp is the desired energy stored in the BESS at the end of the day that is equal to the desired energy at the beginning of the successive day, ηbess is the battery efficiency in charging/discharging mode.

Figure 1 is the slack busbar. Then, VAC

• at the load bus bars, PAC

Table 2. ACµG's constraints.

Energy stored in the generic time slot

30 Development and Integration of Microgrids

Elb

Charging and discharging power in the generic time slot

sp <sup>≤</sup> Ek�<sup>1</sup> � <sup>Δ</sup>t PBESS

�Pch

ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi

<sup>2</sup> <sup>þ</sup> <sup>Q</sup>BESS 3, k 2

PBESS 3, k

q

0 ≤ PBESS

power of the BESS <sup>ð</sup>QBESS

of the micro-turbine;

interfacing converter;

AC

<sup>20</sup>, <sup>k</sup> and <sup>Q</sup>AC

<sup>16</sup>, <sup>k</sup> and <sup>Q</sup>AC

sLB and <sup>V</sup>AC

number of lines of the AC grid and I

In Eq. (12), the active power PAC

• at bus 3, PAC

AC grid

Apparent power

• at bus 20, PAC

• at bus 16, PAC

In Eq. (10), VAC

bus, respectively.

In Eq. (11), I

with Pgrid

<sup>1</sup>, <sup>k</sup> <sup>¼</sup> <sup>V</sup>AC

<sup>3</sup>, <sup>k</sup> <sup>≤</sup> <sup>E</sup>ub

max ≤ PBESS

<sup>3</sup>, <sup>k</sup> ≤ Pdch

≤ SBESS

<sup>3</sup>, <sup>k</sup> ≤ 0 k ∈ Ωch, <sup>k</sup>

max k ∈ Ωdch, <sup>k</sup>

the fixed slack voltage magnitude. In particular, with reference to each kth time slot:

non-controllable loads, whose rated powers are reported in Table A. VI;

<sup>3</sup>, <sup>k</sup> is the sum of the active power of the BESS <sup>ð</sup>PBESS

connected at bus 3 and the reactive power injected by capacitor bank;

<sup>20</sup>, <sup>k</sup> are, respectively, the active <sup>ð</sup>PmTG

power of the non-controllable load connected at bus 3 while QAC

s, <sup>k</sup> and <sup>Q</sup>AC

• at buses, where no generation or loads are connected, the PAC

AC

<sup>k</sup> in Eq. (4)) is limited by contractual active power PMax

sp and ϑ1,<sup>k</sup> = 0, for each kth time slot, where VAC

nom , k ¼ i þ 1, …, Nt ð19Þ

sp , k ¼ i þ 1, …,ðNt � 1Þ : ð17Þ

<sup>3</sup>, <sup>k</sup> Þ and the forecasted active

<sup>20</sup>, <sup>k</sup> <sup>Þ</sup> and reactive <sup>ð</sup>QmTG

s, <sup>k</sup> and <sup>Q</sup>AC

<sup>3</sup>, <sup>k</sup> is the sum of the reactive

s, <sup>k</sup> values are zero.

PCC and the reactive power QAC

<sup>20</sup>, <sup>k</sup> Þ powers

2�3, k

s, <sup>k</sup> are the forecasted active and reactive powers of the

<sup>3</sup>, <sup>k</sup> Þ, the forecasted reactive power of the non-controllable load

<sup>16</sup>, <sup>k</sup> are the active and reactive power injected by the AC/DC grids

sUB are the values of lower and upper bounds of voltage at the sth AC

<sup>2</sup>�3, <sup>k</sup> at the PCC busbar during the <sup>k</sup>th time interval (coincident

l, <sup>k</sup> is the value of current flowing in the lth line at the kth time interval, LAC is the

l,max is the ampacity of the lth line.

sp is

ð18Þ

In Eq. (17), the upper value Eub sp is the battery size while the lower value <sup>E</sup>lb sp must be specified on the basis of the maximum allowable depth of discharge, thus preserving an adequate value of the battery lifetime [18]. Ek–<sup>1</sup> is BESS energy stored at the end of (k-1)th time slot.

In Eq. (18), Pch max and Pdch max are the maximum power rates in charging and discharging modes; Ωdch,k and Ωch,k are the set of time slots from (i + 1)th to Ntth in which the BESS is allowed to discharge and charge, respectively.

In Eq. (19), QBESS <sup>3</sup>, <sup>k</sup> is the reactive power injected by converter of the BESS and SBESS nom is its nominal apparent power.

With reference to AC/DC grids interfacing converter (Table 3), the equality constraints are the DC and AC voltage balance equation (Eq. (20)) and the DC and AC active power balance equation (Eq. (21)), while an inequality constraint (Eq. (22)) refers to a limit on apparent power flowing through the AC/DC grids interfacing converter.

In Eq. (20), ma,k is the amplitude modulation ratio of the PWM interfacing converter, VAC <sup>16</sup>, <sup>k</sup> is the RMS voltage of the bus at AC side of converter and VDC <sup>1</sup>, <sup>k</sup> is the voltage amplitude of the bus at DC side of converter (the first node of the DCµG) at the kth time slot. As before mentioned, VDC <sup>1</sup>, <sup>k</sup> , for each kth time slot, is equal to VDC sp .

In Eq. (21), ηrect (ηinv) is the converter efficiency in rectifier (inverter) operation of the power converter.

In Eq. (22), the active (PAC <sup>16</sup>, <sup>k</sup><sup>Þ</sup> and reactive (QAC <sup>16</sup>, <sup>k</sup>) powers injected by the converter from DCµG to AC µG are limited by its size (Sconv nom ) at the kth time slot.

#### Interfacing converter

#### DC and AC voltage

$$V\_{1,k}^{\rm DC} = \frac{2\sqrt{2}}{m\_{a,k}\sqrt{3}} \ V\_{16,k}^{\rm AC}{}\_{\prime} \quad k = \, i+1, \ldots, N\_t \tag{20}$$

DC and AC active power

$$P\_{16,k}^{\rm AC} = -\xi\_k P\_{1,k}^{\rm DC} \qquad \xi\_k = \begin{cases} \frac{1}{\eta\_{\rm rect}} & \text{if } P\_{16,k}^{\rm AC} \le 0 \\\ \eta\_{\rm inv} & \text{if } P\_{16,k}^{\rm AC} > 0 \end{cases}, k = 1, \ldots, N\_t \tag{21}$$

Apparent power

$$\sqrt{P\_{16,k}^{\rm AC} + Q\_{16,k}^{\rm AC}} \le S\_{nom}^{\rm conv}, \quad k = 1, \ldots, N\_t \tag{22}$$

Table 3. Constraints on AC/DC grids interfacing converter.

#### 3.1.2. Waveform distortion compensation

At harmonic frequencies, the centralized control system furnishes reference current signals that ACµG power converters have to be injected in order to compensate waveform distortions at AC busbars. These reference currents are obtained as the solution of an optimization problem in the frequency domain, based on the optimal control theory [19–21]. The optimization problem consists of minimizing the voltage harmonics at AC busbars without significantly oversizing of the converters employed for the compensation.

The Fourier transform U of the converters reference currents to be injected is obtained by minimizing the following system index:

$$J = \int\_{-\infty}^{+\infty} (\overline{Y}^\* \cdot \mathcal{Q} \cdot \overline{Y} + \overline{U}^\* \cdot \mathcal{R} \cdot \overline{U}) d\omega \tag{23}$$

In Eq. (23), Y is the Fourier transform of the complex output vector y, which includes the harmonic voltages at all network buses and at all considered harmonic frequencies; Q and R are weight matrices and the symbol \* indicates the complex conjugate.

The analytical Fourier transform of the vector of the injected currents U corresponding to the minimization of index of Eq. (23) is given by [20]:

$$
\overline{U} = -\dot{\mathbf{G}}^{-1} \cdot \dot{\mathbf{P}} \cdot \overline{\mathbf{D}} \tag{24}
$$

where <sup>G</sup>\_ ¼ ð<sup>L</sup> � <sup>C</sup> � <sup>Φ</sup>\_ � <sup>B</sup><sup>Þ</sup> � � <sup>L</sup> � <sup>C</sup> � <sup>Φ</sup>\_ � <sup>B</sup> <sup>þ</sup> <sup>R</sup> and <sup>P</sup>\_ ¼ ð<sup>L</sup> � <sup>C</sup> � <sup>Φ</sup>\_ � <sup>B</sup><sup>Þ</sup> � � <sup>L</sup> � <sup>C</sup> � <sup>Φ</sup>\_ � <sup>F</sup>. D is the vector of the disturbances (harmonic currents injected at buses where non-linear loads are present at all considered harmonic frequencies) and Φ ¼ ðjωI � AÞ �<sup>1</sup> being I the identity matrix and Q factorized in Q = L\*L.

The matrices A, B, C, and F are the system state equation matrices obtained from the following linear state space model, derived from the application of the Kirchoff's laws in the time domain to the ACµG [21]:

$$\begin{cases} \frac{d\mathbf{x}}{dt} = A\mathbf{x} + B\mathbf{u} + F\mathbf{d} \\ \mathbf{y} = \mathbf{C}\mathbf{x} \end{cases} \tag{25}$$

where x is the state vector (inductor currents and capacitor voltages), d is the vector of the disturbances, u is the vector of the compensating currents injected by the converters for voltage harmonic compensation, and y is the output vector. A, B, C, and F depend on the parameters of the distribution network (matrices A and C) and on the busbars where the disturbing loads and the converters participating to system harmonic compensation are connected (matrices F and B).

It should be noted that the values of the element of vector D elements are estimated thanks to a disturbance-estimation algorithm based on a Kalman Filter [22, 23].

It should be also noted that the values of the weight matrix Q elements depend on targets assigned to the waveform distortion compensation levels at the nodes, with the rule that the values of the weights are assigned greater in the busbars where particularly sensitive loads are connected (area compensation). The values of the weights of matrix R elements are assigned on the basis of targets on the rate of converter sizes used for the harmonic compensation: greater values of the R elements correspond to the smaller rate of the converter and vice versa [19]. With reference to the considered industrial study system, the converters installed at buses 3, 16 and 20 participate, using a part of their size, to harmonic compensation.

#### 3.2. Islanding mode

3.1.2. Waveform distortion compensation

PAC

<sup>16</sup>, <sup>k</sup> ¼ �ξkPDC

Interfacing converter DC and AC voltage

32 Development and Integration of Microgrids

DC and AC active power

Apparent power

minimizing the following system index:

minimization of index of Eq. (23) is given by [20]:

where <sup>G</sup>\_ ¼ ð<sup>L</sup> � <sup>C</sup> � <sup>Φ</sup>\_ � <sup>B</sup><sup>Þ</sup>

At harmonic frequencies, the centralized control system furnishes reference current signals that ACµG power converters have to be injected in order to compensate waveform distortions at AC busbars. These reference currents are obtained as the solution of an optimization problem in the frequency domain, based on the optimal control theory [19–21]. The optimization problem consists of minimizing the voltage harmonics at AC busbars without signifi-

The Fourier transform U of the converters reference currents to be injected is obtained by

In Eq. (23), Y is the Fourier transform of the complex output vector y, which includes the harmonic voltages at all network buses and at all considered harmonic frequencies; Q and R

The analytical Fourier transform of the vector of the injected currents U corresponding to the

� � <sup>L</sup> � <sup>C</sup> � <sup>Φ</sup>\_ � <sup>B</sup> <sup>þ</sup> <sup>R</sup> and <sup>P</sup>\_ ¼ ð<sup>L</sup> � <sup>C</sup> � <sup>Φ</sup>\_ � <sup>B</sup><sup>Þ</sup>

<sup>U</sup> ¼ �G\_ �<sup>1</sup>

<sup>ð</sup>Y� � <sup>Q</sup> � <sup>Y</sup> <sup>þ</sup> <sup>U</sup>� � <sup>R</sup> � <sup>U</sup>Þd<sup>ω</sup> <sup>ð</sup>23<sup>Þ</sup>

<sup>16</sup>, <sup>k</sup> , k ¼ i þ 1, …, Nt ð20Þ

, k ¼ i þ 1, …, Nt

nom , k ¼ i þ 1, …, Nt ð22Þ

ð21Þ

� <sup>P</sup>\_ � <sup>D</sup> <sup>ð</sup>24<sup>Þ</sup>

� � <sup>L</sup> � <sup>C</sup> � <sup>Φ</sup>\_ � <sup>F</sup>.

cantly oversizing of the converters employed for the compensation.

VDC

<sup>1</sup>, <sup>k</sup> <sup>¼</sup> <sup>2</sup> ffiffiffi

<sup>1</sup>, <sup>k</sup> ξ<sup>k</sup> ¼

ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi

<sup>2</sup> <sup>þ</sup> <sup>Q</sup>AC 16, k 2

PAC 16, k

q

Table 3. Constraints on AC/DC grids interfacing converter.

ma, <sup>k</sup> ffiffiffi 3 <sup>p</sup> <sup>V</sup>AC

2 p

> 1 ηrect

≤ Sconv

8 < :

if PAC <sup>16</sup>, <sup>k</sup> ≤ 0

<sup>16</sup>, <sup>k</sup> > 0

<sup>η</sup>inv if PAC

J ¼ ðþ<sup>∞</sup> �∞

are weight matrices and the symbol \* indicates the complex conjugate.

In addition, the control strategy applied by CCS is able to disconnect the hybrid μG from the distribution network when abnormal conditions occur. We will consider the behavior of the AC/DC µG in the case of supply interruptions. In such conditions, local energy resources and loads are optimally controlled in order to guarantee high power quality levels to priority loads. As priority loads, we refer to electrical loads requiring high levels of continuity of supply.

In particular, in islanding mode, the so-called "Single Mater Operation" control strategy is adopted [24]. In fact, while in grid-connected mode all inverters installed into the hybrid μG operate in PQ mode (inverter is used to supply a given active P and reactive Q power set point), in islanding mode the BESS interface converter operates feeding the μG with predefined values of voltage and frequency (V-f control mode), as well as the remaining inverters, on the contrary, continue to operate in PQ mode.

As soon as the hybrid μG is in the islanding mode, all loads are disconnected except the priority loads whose power demand can be provided by the BESS and micro-turbine (the PV system is not considered in the power balance because its availability depends on the uncertain presence of the primary source). As an example, in the considered case study, the priority loads are assumed to be the sensitive loads connected to the DCµG. In this case, the maximum power required by the priority loads is 72 kW (corresponding to the power absorbed when all sensitive loads are working). Also in the worst conditions (long interruption arising when BESS charge is Elb sp = 102 kWh), the energy that should be generated by BESS and micro-turbine in 1 hour<sup>8</sup> (equal to 72 kWh) is significantly lower than the energy they can provide (equal to 132 kWh) .9

Eventually, in islanding mode, the proposed control strategy starts to supply the priority loads. Anyway, according to a predefined list of load priorities, some of the disconnected loads can be reconnected taking into account of: (i) the available power from PV generator; (ii) the scheduled power required by sensitive loads; (iii) the maximum time fixed for the autonomous operation of the hybrid µG; and (iv) the value of the SOC of battery. The available energy for the islanding operation can be easily calculated from the knowledge of the SOC of battery, the forecasted PV power generation in the time intervals successive to the disconnection and the rated power of micro-turbine.

### 4. Numerical applications

A simulation of the hybrid micro grid of Section 2 using Matlab® Simulink® toolbox was implemented in order to show the feasibility of the proposed control strategy both in gridconnected and in islanding mode. Several case studies were performed, and, for sake of conciseness, the results relative to a typical working day are shown. In particular, in the following subsections, results are presented for the two considered operation modes, separately.

Remind that, in the grid-connected mode, the converters participating in the energy cost minimization problem are: (i) the BESS converter at bus 3; (ii) the AC/DC grids interfacing converter at bus 16; and (iii) the micro-turbine converter at bus 20. In particular, the AC/DC grids interfacing converter (micro-turbine) participates optimally exchanging only reactive (active) power with the ACµG.

<sup>8</sup> The maximum time interval in islanding mode operation of hybrid µG was fixed to 60 minutes considering the mean value of the time duration of long interruptions in the area where the industrial facility is located.

<sup>9</sup> We assume that in the worst conditions the PV production is not available and the SOC of battery is initially equal to its minimum value (20%). In these conditions, in order to supply the sensitive loads (72kWh), the micro-turbine provides 30kWh and the BESS has to furnish only 42 of the 102kWh available. This, however, determines a further discharge of the BESS and a final SOC equal to 11.76%, causing a reduction of the expected battery life. However, since 60 minutes long interruptions occur in distribution networks only few times during the battery life and taking into account that the worst conditions are infrequent, the actual life reduction of the BESS can be considered largely tolerable.

The converters participating in the harmonic compensation are: (i) the AC/DC grids interfacing converter at bus 16 and (ii) the micro-turbine converter at bus 20. Note also that sizes of the converters, shown in Table A. VIII, are chosen according to the aforesaid operating conditions. Specifically, while the 120 kVA of the BESS converter is fully utilized for the energy cost minimization, only 80 kVA of the AC/DC grids interfacing converter and 30 kVA of the microturbine converter are addressed for the same operation. The remaining 20 kVA of the AC/DC grids interfacing converter and 10 kVA of the micro-turbine converter are used for the harmonic compensation.

#### 4.1. Grid-connected mode

As soon as the hybrid μG is in the islanding mode, all loads are disconnected except the priority loads whose power demand can be provided by the BESS and micro-turbine (the PV system is not considered in the power balance because its availability depends on the uncertain presence of the primary source). As an example, in the considered case study, the priority loads are assumed to be the sensitive loads connected to the DCµG. In this case, the maximum power required by the priority loads is 72 kW (corresponding to the power absorbed when all sensitive loads are working). Also in the worst conditions (long interruption arising when

in 1 hour<sup>8</sup> (equal to 72 kWh) is significantly lower than the energy they can provide (equal to

Eventually, in islanding mode, the proposed control strategy starts to supply the priority loads. Anyway, according to a predefined list of load priorities, some of the disconnected loads can be reconnected taking into account of: (i) the available power from PV generator; (ii) the scheduled power required by sensitive loads; (iii) the maximum time fixed for the autonomous operation of the hybrid µG; and (iv) the value of the SOC of battery. The available energy for the islanding operation can be easily calculated from the knowledge of the SOC of battery, the forecasted PV power generation in the time intervals successive to the disconnection and the

A simulation of the hybrid micro grid of Section 2 using Matlab® Simulink® toolbox was implemented in order to show the feasibility of the proposed control strategy both in gridconnected and in islanding mode. Several case studies were performed, and, for sake of conciseness, the results relative to a typical working day are shown. In particular, in the following subsections, results are presented for the two considered operation modes, sepa-

Remind that, in the grid-connected mode, the converters participating in the energy cost minimization problem are: (i) the BESS converter at bus 3; (ii) the AC/DC grids interfacing converter at bus 16; and (iii) the micro-turbine converter at bus 20. In particular, the AC/DC grids interfacing converter (micro-turbine) participates optimally exchanging only reactive

The maximum time interval in islanding mode operation of hybrid µG was fixed to 60 minutes considering the mean

We assume that in the worst conditions the PV production is not available and the SOC of battery is initially equal to its minimum value (20%). In these conditions, in order to supply the sensitive loads (72kWh), the micro-turbine provides 30kWh and the BESS has to furnish only 42 of the 102kWh available. This, however, determines a further discharge of the BESS and a final SOC equal to 11.76%, causing a reduction of the expected battery life. However, since 60 minutes long interruptions occur in distribution networks only few times during the battery life and taking into account that the worst

value of the time duration of long interruptions in the area where the industrial facility is located.

conditions are infrequent, the actual life reduction of the BESS can be considered largely tolerable.

sp = 102 kWh), the energy that should be generated by BESS and micro-turbine

BESS charge is Elb

34 Development and Integration of Microgrids

rated power of micro-turbine.

4. Numerical applications

(active) power with the ACµG.

132 kWh) .9

rately.

8

9

The considered day was divided into 72-time slot of equal duration (20 min), and the inputs of the proposed control strategy for the energy cost minimization are:


Figure 4. Forecasted power profile for the first time slot: (a) power production of the PV system; (b) active (solid line) and reactive (dotted line) power of the non-controllable load at bus 33.

Figure 5. Scheduled profile of the power of the controllable loads.

Figure 6. Energy charge.

The time behavior of the main quantities of interest is illustrated in the following relatively to the considered working day.

Figure 7 shows the active (solid line) and reactive (dotted line) power profiles of the BESS converter. In Figure 8, the active (solid line) and reactive (dotted line) power profiles of the AC/DC grids interfacing converter are plotted. Figure 9 shows the reactive power supplied by the controllable capacitor bank. Eventually, in Figure 10, the active (solid line) and reactive (dotted line) power profiles at the PCC are reported. The active power of the micro-turbine is not reported since it always operates at a constant value of 30 kW.

From the analysis of Figures 7–10, the following considerations can arise:

Figure 5. Scheduled profile of the power of the controllable loads.

36 Development and Integration of Microgrids

Figure 6. Energy charge.


Figure 7. Active (solid line) and reactive (dotted line) power profiles related to the BESS converter at bus 3.

Figure 8. Active (solid line) and reactive (dotted line) power profiles related to the AC/DC grids interfacing converter at bus 16.

Figure 9. Reactive power profile of the capacitor bank at bus 3.

Distributed Energy Resources to Improve the Power Quality and to Reduce Energy Costs of a Hybrid AC/DC… http://dx.doi.org/10.5772/intechopen.68766 39

Figure 10. Active (solid line) and reactive (dotted line) power profiles at the PCC.

Figure 8. Active (solid line) and reactive (dotted line) power profiles related to the AC/DC grids interfacing converter at

bus 16.

38 Development and Integration of Microgrids

Figure 9. Reactive power profile of the capacitor bank at bus 3.

Anyway, the BESS totally discharges (i.e. the SOC reaches 20%) in the interval where the tariff is higher, maximizing in this way the economic benefit.

Note also that the reactive power production inside the µG allows not only the reduction of the power losses but also the maximization of the internal active power generation, determining the minimization of the objective function. In fact, an increasing of the mere internal active power production determines a corresponding decreasing of power factor at the PCC; therefore, a production of reactive power inside the hybrid grid is needed to guarantee an adequate power factor at the PCC.

In Figure 11, a comparison among the energy charge profile (dotted line) and the active power (solid line) of BESS converter (Figure 11a) and the active power of AC/DC grids interfacing converter (Figure 11b) is reported. Figure 12 shows a comparison among the energy charge profile (dotted line), the active power profile at the PCC of the hybrid µG (solid line) and of the actual industrial facility (dashed line).

The analysis of these last figures clearly evidences that:


Figure 11. Comparison among the energy charge profile (dotted line), the active power profiles (solid line) related to: (a) the BESS converter at bus 3 and (b) to the AC/DC grids interfacing converter at bus 16.

Figure 12. Comparison among the energy charge profile (dotted line), the active power profile at the PCC of the industrial system with µG (solid line) and of the actual industrial system (dashed line).

• the active power at the PCC of hybrid µG increases respect to the value of the existing industrial facility when the tariff price is minimum, due to: (i) the BESS charge profile and (ii) the shifting of the controllable loads; on the other hand, the active power at the PCC of hybrid µG decreases when the tariff price assumes higher values, due to: (i) the shifting of the controllable loads; (ii) the PV system production; and (iii) the BESS discharge profile.

The aforesaid observations justify the reduction, shown in Figure 13, of the hourly costs to be sustained by the hybrid µG (solid line) respect to the ones sustained by the existing industrial facility (dotted line). In particular:


Globally, the proposed procedure allows a daily reduction of the total costs for the energy purchase of about 26%.

Figure 13. Hourly costs to be sustained by the industrial system with hybrid µG (solid line) and by the actual industrial system (dotted line).

• the active power at the PCC of hybrid µG increases respect to the value of the existing industrial facility when the tariff price is minimum, due to: (i) the BESS charge profile and

Figure 12. Comparison among the energy charge profile (dotted line), the active power profile at the PCC of the

industrial system with µG (solid line) and of the actual industrial system (dashed line).

Figure 11. Comparison among the energy charge profile (dotted line), the active power profiles (solid line) related to: (a)

the BESS converter at bus 3 and (b) to the AC/DC grids interfacing converter at bus 16.

40 Development and Integration of Microgrids

With reference to the waveform distortion compensation, an analysis of the load currents in the existing industrial system revealed that main contributions to voltage distortions were due to the harmonic of order h = 5, 7, 11 and 13. So, the procedure described in Section 3.1.2 was applied only for the compensation of the aforesaid harmonic orders. The compensation action is applied in a working day using the AC/DC grids interfacing converter and the micro-turbine converter.

As illustrated in Section 3, the inputs needed to solve the minimization problem and to obtain the reference currents of converters (see Eq. (23)) are the estimated harmonic disturbances D and the Q and R matrices.

The loads producing waveform distortion are marked with an asterisk in Table A. VI. As an example, Figures 14 and 15 report the harmonic currents amplitudes injected by disturbing loads at the buses 23 and 33 during the time slot from 16:20 to 16:40.

The terms of the weight matrix R are obtained from the knowledge of the size of the converters dedicated to harmonic compensation, and giving priority to the compensation of the harmonic with the lowest orders. As an example, the value of R term related to the compensation of the 5th harmonic order is 0.59, while the value of R term related to the compensation of the 13th harmonic order is 0.23.

Figure 14. DFT estimation of harmonic currents amplitudes injected by disturbing loads at the bus 23 during the time slot from 16:20 to 16:40.

Distributed Energy Resources to Improve the Power Quality and to Reduce Energy Costs of a Hybrid AC/DC… http://dx.doi.org/10.5772/intechopen.68766 43

Figure 15. DFT estimation of harmonic currents amplitudes injected by disturbing loads at the bus 33 during the time slot from 16:20 to 16:40.

The terms of weight matrix Q are set giving priority to the harmonic compensation for the feeder "tanks and boxes manufactory"; thus the elements of Q corresponding to nodes of this feeder are 50 times greater than those of the other feeders.

Figures 16–18 show the comparison of total harmonic distortion (THD) with (solid bars) and without (empty bars) the proposed compensation procedure. Specifically, Figure 16 shows the THD comparison over the buses in the time slot corresponding to the time interval from 16:20 to 16:40, while Figures 17 and 18 show the THD comparison, in the whole day, for the buses 23 and 35, respectively. The benefits in term of THD reduction with the application of the proposed procedure are evident and significant, since, as observed in Figures 16–18, the THD values in all of the buses of the hybrid µG in time stayed always under the 0.6% with the waveform distortion compensation. Without compensation, the THD overcame the 4.5% in some buses for some time slots.

#### 4.2. Islanding mode

With reference to the waveform distortion compensation, an analysis of the load currents in the existing industrial system revealed that main contributions to voltage distortions were due to the harmonic of order h = 5, 7, 11 and 13. So, the procedure described in Section 3.1.2 was applied only for the compensation of the aforesaid harmonic orders. The compensation action is applied in a working day using the AC/DC grids interfacing converter and the micro-turbine

As illustrated in Section 3, the inputs needed to solve the minimization problem and to obtain the reference currents of converters (see Eq. (23)) are the estimated harmonic disturbances D

The loads producing waveform distortion are marked with an asterisk in Table A. VI. As an example, Figures 14 and 15 report the harmonic currents amplitudes injected by disturbing

The terms of the weight matrix R are obtained from the knowledge of the size of the converters dedicated to harmonic compensation, and giving priority to the compensation of the harmonic with the lowest orders. As an example, the value of R term related to the compensation of the 5th harmonic order is 0.59, while the value of R term related to the compensation of the 13th

Figure 14. DFT estimation of harmonic currents amplitudes injected by disturbing loads at the bus 23 during the time slot

loads at the buses 23 and 33 during the time slot from 16:20 to 16:40.

converter.

and the Q and R matrices.

42 Development and Integration of Microgrids

harmonic order is 0.23.

from 16:20 to 16:40.

To simulate islanding mode of operation of the hybrid μG, an interruption in the MV distribution network supplying the industrial facility was simulated. The interruption determines the switching from grid-connected mode to islanding mode. The control strategy described in Section 3.2 is applied; the BESS converter switches from P-Q control mode to V-f control mode

Figure 16. THD comparison with (solid bars) and without (empty bars) the proposed compensation procedure corresponding to the time interval from 16:20 to 16:40.

Figure 17. THD comparison with (solid bars) and without (empty bars) the proposed compensation procedure corresponding to the bus 23.

Distributed Energy Resources to Improve the Power Quality and to Reduce Energy Costs of a Hybrid AC/DC… http://dx.doi.org/10.5772/intechopen.68766 45

Figure 18. THD comparison with (solid bars) and without (empty bars) the proposed compensation procedure corresponding to the bus 35.

Figure 16. THD comparison with (solid bars) and without (empty bars) the proposed compensation procedure corresp-

Figure 17. THD comparison with (solid bars) and without (empty bars) the proposed compensation procedure corresp-

onding to the time interval from 16:20 to 16:40.

44 Development and Integration of Microgrids

onding to the bus 23.

since the RMS voltage at PCC becomes lower than 5% of nominal voltage [26], that is the threshold for the switching of the BESS converter control strategy. The remaining inverters, on the contrary, continue to operate in P-Q mode. The interruption arises at 6 h 25<sup>0</sup> 1<sup>00</sup> of the working day analyzed in Section 4.1. As soon as the islanding condition arises, all loads of the hybrid μG are disconnected except the sensitive loads present in the DC μG. As illustrated in Figure 5, at 6 h 25<sup>0</sup> all sensitive loads are working and absorbing 72 kW, micro-turbine and photovoltaic generators produce 30 and 31 kW, respectively. BESS is characterized by a SOC of about 100%.

In Figure 19, the voltage profile at bus 3 versus time is reported. In Figure 20, the active power injected by BESS and AC/DC grids interfacing converters are shown. Figures 19 show how the proposed approach allows to supply the hybrid μG during the interruption with obvious advantages in terms of continuity of the supplied energy.

Figure 20a shows the increasing of BESS active power from the value at 6 h 25<sup>0</sup> in grid-connected mode (see also Figures 7 and 11) to the value needed to supply the hybrid μG in islanded mode. In particular, this value, in the considered application, is related to the grid active power losses and active power required by the sensitive loads connected to DCμG reduced by the active power furnished by micro-turbine and PV system. As expected, the active power flowing in the AC/DC grids interfacing converter remains constant to the value of grid-connected mode (see also Figures 8 and 11).

In Figure 20b, reactive power variations from the values at 6 h 25<sup>0</sup> in grid-connected mode (see also Figures 7 and 8) to the values needed to provide the reactive power required by the hybrid

Figure 19. Phase-A voltage at bus 3 versus time starting from 6 h 25<sup>0</sup> .

μG in islanding mode is shown. In particular, since only the sensitive loads connected to DC μG are working, only reactive power losses are required and provided by BESS converter.

Note that the active and reactive power transients shown in Figure 20 are due to the switching of the BESS converter regulator from P-Q control mode to voltage control mode. Moreover, although the other converters in the hybrid µG continue to be controlled in P-Q mode, also their responses are influenced by the transient at bus 3.

Figure 20. Active (a) and reactive (b) powers injected by BESS, micro-turbine and AC/DC grids interfacing converter starting from 6 h 25<sup>0</sup> .

Accordingly to the aim of the proposed procedure, all the distributed resources in the hybrid μG guarantee the continuity of supplying of sensitive loads.

### 5. Conclusions

A hybrid µG including two interconnected sections operating in DC and AC was considered starting from an existing Italian industrial facility. The DC section included non-dispatchable generation systems and controllable loads; the AC section included a BESS, a micro-turbine, non-controllable loads and non-linear loads.

A control strategy based on a real-time control, at both the fundamental frequency and harmonic frequencies, was applied. The control strategy required the solution of optimization problems and was aimed at minimizing the operation costs and improving the power quality levels.

The results of a numerical application relative to a typical working day are shown and demonstrated the effectiveness of the procedure in reducing the daily operating costs while guaranteeing higher values of the quality of the power supply.

### Acknowledgements

This work was supported by Italian Ministero dell'Istruzione dell'Università e della Ricerca (MIUR) in the frame of PON03PE\_00178\_1: Microgrid Ibride in Corrente Continua e Corrente Alternata (MICCA) and by Università of Napoli Parthenope in the framework of "Bando per il Sostegno alla Ricerca individuale triennio 2015–2017".

### Appendix

μG in islanding mode is shown. In particular, since only the sensitive loads connected to DC μG

.

Note that the active and reactive power transients shown in Figure 20 are due to the switching of the BESS converter regulator from P-Q control mode to voltage control mode. Moreover, although the other converters in the hybrid µG continue to be controlled in P-Q mode, also

Figure 20. Active (a) and reactive (b) powers injected by BESS, micro-turbine and AC/DC grids interfacing converter

are working, only reactive power losses are required and provided by BESS converter.

their responses are influenced by the transient at bus 3.

starting from 6 h 25<sup>0</sup>

.

Figure 19. Phase-A voltage at bus 3 versus time starting from 6 h 25<sup>0</sup>

46 Development and Integration of Microgrids

In the following, the data of the lines are provided. In particular, for each line, Tables A. I, A. II, A. III, and A. IV show the starting and ending buses, the length, resistance, and reactance. Tables A. V and A. VI report the transformers and loads data. Tables A. VII and A. VIII report the main data of the new components installed in the proposed hybrid μG.



Table A.I. Line parameters for the feeder "Tanks and boxes manufactory".


Distributed Energy Resources to Improve the Power Quality and to Reduce Energy Costs of a Hybrid AC/DC… http://dx.doi.org/10.5772/intechopen.68766 


Table A. II. Line parameters for the feeder "Assembly".

Buses ℓ [m] R [mΩ/m] X [mΩ/m]

 10 9.3 0.163 0.130 11 3 0.163 0.130 12 2.8 0.163 0.130 13 3.5 0.163 0.130 14 3.5 0.163 0.130 15 11 0.163 0.130 16 19.1 0.163 0.130 17 4 1.410 0.112 18 1.9 0.163 0.130 19 4 1.410 0.112 20 3 0.163 0.130 21 10 0.236 0.097 22 42 1.410 0.112 23 61 2.240 0.119 24 61 1.410 0.112

Buses ℓ [m] R [mΩ/m] X [mΩ/m]

 25 31 0.041 0.014 26 3.5 0.163 0.130 27 16 2.240 0.119 28 9 0.163 0.130 29 12.5 0.163 0.130 30 10.5 0.163 0.130 31 8 0.641 0.101 32 1.5 0.163 0.130 33 10 0.641 0.101 34 12.5 0.163 0.130 35 8 0.328 0.096 36 13 0.163 0.130 37 1 0.163 0.130 38 35 0.665 0.260

Table A.I. Line parameters for the feeder "Tanks and boxes manufactory".

From To

Development and Integration of Microgrids

From To


Table A. III Line parameters for the feeder "Winding and coils".

#### 50 Development and Integration of Microgrids


Table A. IV Line parameters for the feeder "Test".


#### Table A. V. Transformer data.


Distributed Energy Resources to Improve the Power Quality and to Reduce Energy Costs of a Hybrid AC/DC… http://dx.doi.org/10.5772/intechopen.68766 51


Table A. VI Load data.

Buses ℓ [m] R [mΩ/m] X [mΩ/m]

 66 30 0.094 0.090 67 87 0.163 0.130 68 7 1.410 0.112 69 0.5 0.163 0.130 70 7 0.473 0.101

Nominal voltage 20 kV/400 V Nominal power 630 kVA Connection Dyn Group 11 Short circuit voltage 6.0% Load losses 5721 W No load losses 925 W No load current 5%

Bus Type Rated power [kVA] Power factor

3 Heating and cooling system 200 0.8 6 Painting machine 75 0.8 7 Box overturning machine 4 0.99 9\* Sandblasting machine 55 0.75 11 Welder aspirators 11 0.99 15 Manual bender 8 0.99 17\* Folding walls island robot 24 0.99 19\* Wave welding machine 30 0.65 21 Automated bending robot for corrugated panels 122 0.65 22\* PLC + computer 3 0.62 23\* Plasma cutting machine 15 0.8 27 Crane 5.5 0.8 31\*, 33\* Core cutting machine n. 1 60 0.8 35 Autoclaves 86 0.8 44 Furnace 5 0.99

From To

50 Development and Integration of Microgrids

Table A. V. Transformer data.

Table A. IV Line parameters for the feeder "Test".


1 5 35 0.660 0.00029

Table A. VII DCμG components: PV generator (a) and lines data (b).


Table A. VIII ACμG components data.

### Author details

Luisa Alfieri<sup>1</sup> , Antonio Bracale1 \*, Pierluigi Caramia<sup>1</sup> and Guido Carpinelli<sup>2</sup>


### References


of the International Conference on Renewable Energies and Power Quality (ICREPQ15); 25- 27 March 2015; La Coruña, Spain


Author details

52 Development and Integration of Microgrids

, Antonio Bracale1

1 University of Napoli Parthenope, Naples, Italy

2 University of Napoli Federico II, Naples, Italy

\*Address all correspondence to: antonio.bracale@uniparthenope.it

tems and Technologies (EDST15), September 2015

IEEE Transactions on Smart Grid. 2014;5:1494-1505

IEEE Transactions on Smart Grid. 2011;2:278-286

Transmission & Distribution. 2014;8:1700-1711

IEEE Transactions on Industrial Electronics. 2013;60:1263-1270

\*, Pierluigi Caramia<sup>1</sup> and Guido Carpinelli<sup>2</sup>

[1] CIGRÉ. Working Group C6.22 Microgrids Evolution Roadmap, Microgrids 1: Engineer, Economics, & Experience. International Symposium on Smart Electric Distribution Sys-

[2] Hatziargyriou N. Microgrids: Architectures and Control. Wiley-IEEE Press; United King-

[3] Farhangi H. The path of the smart grid. IEEE Power and Energy Magazine. 2010;8:18-28 [4] Hatziargyriou N, Asano H, Iravani R, Marnay C. Microgrids. IEEE Power and Energy

[5] Eghtedarpour N, Farjah E. Power control and management in a hybrid AC/DC microgrid.

[6] Guerrero JM, Loh PC, Tzung-Lin L, Chandorkar M. Advanced control architectures for intelligent microgrids—Part II: Power quality, energy storage, and AC/DC microgrids.

[7] Xiong L, Peng Wang P, Loh C. A hybrid AC/DC microgrid and its coordination control.

[8] Loh PC, Ding L, Yi Kang C, Blaabjerg F. Autonomous operation of hybrid microgrid with AC and DC subgrids. IEEE Transactions on Power Electronics. 2013;28:2214-2223

[9] Long B, Jeong TW, Lee JD, Jung YC, Chong KT. Energy management of a hybrid AC–DC

[10] Baboli PT, Shahparasti M, Moghaddam MP, Haghifam MR, Mohamadian M. Energy management and operation modeling of hybrid AC-DC microgrid. IET Generation,

[11] Bracale A, Caramia P, Carpinelli G, Mancini E, Mottola F. Optimal control strategy of a DC micro grid. International Journal of Electrical Power & Energy Systems. 2015;67:25-38

[12] Alfieri L, Carpinelli G, Proto D, Russo G. Day-ahead optimal scheduling of loads and dispatchable resources in a hybrid AC/DC microgrid of an industrial system. In: Proceedings

micro-grid based on a battery testing system. Energies. 2015;8:1181-1194

Luisa Alfieri<sup>1</sup>

References

dom, 2014

Magazine. 2007;5:78-94

[26] Standard CEI EN 50160: Voltage characteristics of electricity supplied by public distribution systems. 2015.
