#### Abstract

A distributed cooperative control scheme is proposed in order to implement a distributed secondary control for hybrid lossy microgrids. The designed distributed control is able to synchronize the frequency of inverse-based distributed generators (DGs) and minisynchronous generators (MSGs/SGs) to the desired state with a virtual leader DG/SG (reference value) in a distribution switching network under the existence of time-varying communication delays. The secondary control stage selects suitable frequencies of each DG/SG such that they can be synchronized at the desired set point. Using the proposed algorithm, each DG/SG only needs to communicate with its neighboring DGs/SGs intermittently even if the communication networks are local, the topology is time-varying, and the communication delays may exist. Therefore, the failure of a single DG/SG will not produce the failing down of the whole system. Sufficient conditions on the requirements for the network connectivity and the delays boundedness which guarantees the stability and synchronization of the controlled hybrid lossy microgrid power systems are presented. The feasibility of the proposed control methodology is verified by the simulation of a given lossy microgrid test system.

#### 1. Introduction

The recent changes in the structure of power generation towards a distributed generation have motivated the increasing interest in the control of so-called microgrids. As the main building blocks of smart grids, microgrids are small-scale power systems that facilitate the effective integration of highly hybrid and heterogeneous DGs and storage devices including solar (photovoltaic array), wind, microturbines, supercapacitor, and batteries [1, 2]. The structure of hybrid microgrid is as shown in Figure 1. Many of these energy sources and storage devices generate or reserve variable frequency AC/DC power and are interfaced with a synchronous AC grid via power electronic DC/AC inverters [3]. Meanwhile, with the rapid growth of power electronics techniques on MSGs or virtual SGs, the increasing application of clear energy sources, including diesel plants, wind plants, and geothermal plants, also makes them possible to be integrated into microgrids. Therefore, new control strategies on how to preserve the synchronization and proportional power sharing in the connected or isolated manner are demanded for the network [4].

Recently, hierarchical control for microgrids has been proposed in order to standardize their operation and functionalities, that is, primary control [5] (droop control maintaining voltage and frequency stability of the microgrid subsequent to the islanding process), secondary control [6] (compensating the voltage and frequency derivations caused by primary control), and tertiary control [7] (optimal operation in both operating modes and power flow control in grid-connected mode).

However, the conventional secondary control of microgrids assumes a centralized control structure which requires a complex communication network and a central controller. A central controller reduces the system reliability [8]. Additionally, some literatures point that physical and communication structures of microgrid can be time-varying because of the desired plug-and-play capability of microgrids [9]. From this perspective, the distributed control structure provides a robust secondary control framework [10, 11]. It could largely reduce the system reliability [12, 13]. In contrast, decentralised control structures are adopted in [14, 15] to avoid the above drawbacks. Nevertheless, for the large-scale microgrids, it is impossible to determine and maintain a feasible operating condition if all of the DGs are running independently. Alternatively, affected by the idea of cooperation control for multiagent systems [16–21], there are increasing researchers who begin to focus on the distributed control structures which allow all the DGs to communicate with their neighbors and share information among neighboring units via local sparse communication networks [22–26].

Particularly, in paper [27], the authors proposed a secondary voltage and frequency control scheme based on the distributed cooperative control of multiagent systems and thus improved the system reliability. In order to regulate the power output of a microgrid dynamics consisting of a large amount of photovoltaic generators, [28] presented a self-organizing strategy in a distributed time-varying network.

The main common features of the above literatures lie in that a distributed communication network has been taken as a basis for the control strategies designed. In this case, the network variability and communication delays are usually unavoidable. They cannot be neglected for the system stability analysis [29, 30]. In fact, in normal circumstances, the microgrids frequency fluctuation and the amplitudes of each DG unit response to primary frequency control are very small. For the purpose of efficient signal transmission, it is inevitable to induce some disturbance and time-delays during the signal processing. Hence, in [29], decentralized inverter control algorithms for power sharing with nonnegligible constant communication delays are investigated. Also, by using the small-signal analysis method, a distributed secondary voltage and frequency control for reactive power sharing are discussed in [30], in which the constant communication delays are analyzed only in an experimental simulation way.

In view of the above analysis, most research on synchronization and power sharing of microgrids has focused on purely inverter-based systems. However, from a practical consideration, most present and near-future applications involve networks of mixed generation mode including SGs and inverter-based DGs and/or storage devices [31]. This paper aims to develop a distributed cooperative control scheme for hybrid microgrids that is able to synchronize the frequency of DGs/SGs to the desired state with a virtual leader DG/SG (reference value) in a distribution switching network under the existence of time-varying communication delays. Combining with the conventional droop control technique, the distributed cooperative control algorithm is designed and deployed into the secondary control stage, in which manner the control derivation produced during the primary control stage can be well restored. And the secondary control stage selects suitable frequencies of each DG/SG such that they can synchronize to the desired state. With the proposed algorithms, each DG/SG only needs to communicate with its neighboring DGs/SGs intermittently even if the communication networks are local, the topology is time-varying, and the communication delays may exist. Sufficient conditions on the requirements for the network connectivity and the delays boundedness that guarantees the stability and reliability of the controlled microgrid power systems are presented. The effectiveness of the proposed control methodology is verified by the simulation of a lossy microgrid test system.

The main contributions of this paper lie in the following aspects. Firstly, in order to restore the frequency in hybrid lossy microgrids, one first designs a distributed secondary control algorithm to synchronize the frequency of each DG/SG to the desired state (reference value 50 Hz). Secondly, the considered distributed communication network is local and time-varying and with time-varying communication delays. By Lyapunov-Krasovskii stable analysis technique, sufficient conditions on the requirements for the network connectivity and the delays boundedness that guarantees the stability and synchronization of the controlled microgrid power systems are presented.

The rest of this paper is organized as follows. In Section 2, the droop-based distributed cooperative control for microgrids problem is formulated and several necessary lemmas are then given. The distributed cooperative algorithm is designed for realizing the frequency restoration in Section 3. Section 4 gives some typical examples to verify the effectiveness of the proposed algorithm. Finally, some concluding remarks are drawn in Section 5.

#### 2. Problem Formulation

In normal situation, the microgrid is operating with the nominal frequency and voltage magnitudes. The objective of primary control is to maintain the frequency and/or voltage magnitudes near around its nominal value once user demands and/or power supply changes. Herein our attentions are restricted to the frequency regulation, so the inverters are considered as voltage source with constant voltage amplitude. And we adopt the physically reasonable and widely adopted assumption of an inductive microgrid [26].

The procedure of primary and secondary control is as shown in Figure 2. In the steady state, the load power curve (after linearization processing) and the generator power curve (after linearization processing) have a balance point at , whose coordinate is . The frequency will decrease by in order to make the output power of the DG/SG increase to the point . Due to the inherent defects of primary control, the increase of the generator power in the above case is smaller than the actual increase in the load power . In order to compensate the frequency deviation produced by the primary control, secondary frequency control (the green lines in Figure 2) is deployed to regulate the frequency such that the intersection point moves to .

##### 2.1. Droop Characters of SGs in Electrical Networks

For the dynamics of SGs, we consider the conventional expression of a regulated SG by a constant voltage behind transient reactance [24]. The so-called regulated machine is always equipped with a speed governor connected to a governing system. The speed governor is controlled to adjust the output of turbine mechanical power and its speed. Droop control is the intrinsic characteristic to maintain system stability. Here, we assume that there is an approximate linear relationship between the valve position and the mechanical power. It is worth noting that the mechanical speed is different from the electrical speed , but their relationship is , where is the number of machine poles [31]. Droop control eventually can be expressed as follows:where is an input signal, is an output signal, and denotes the droop coefficient. The is a constant and the nominal set point for the mechanical power.

##### 2.2. Droop Control for Inverter-Based DGs Networks

From the above theoretical analysis, the droop technique deployed in primary controller for DGs during the multiple loop control process can be given bywhere is chosen from the nominal set point of frequency of the DG unit, is the measured active powers at the DG terminal, is the associated droop coefficient that is usually selected based on the active power rating, and is chosen from the nominal set point of active power of the DG unit.

The control process of each DG generally consists of three control loops (as shown in Figure 3), that is, the power controller, voltage magnitudes controller, and current controller.

The active power and reactive power of the th DG are calculated based on the low-pass filter (with cutoff frequency ), measured output voltage magnitude, and output current. They can be described by the following equations, respectively:

The nominal value of output frequency of the power controller is used by the sinusoidal pulse width modulation (SPWM) inverter as frequency reference, while the voltage reference of the SPWM inverter needs to be regulated by the following voltage and current controllers. Based on the references provided by the power controller, and , the output of the voltage magnitude controller is given bywhere and are the proportional and integral gains of the voltage magnitudes controller of th DG, respectively, while and are the capacitance of the LC filter and feedforward gain, respectively. On the other hand, based on the references provided by the voltage magnitudes controller, and , the output of the current controller is given by

The dynamics of a frequency droop-controlled inverter (2) and a regulated SG (1) are equivalent (see [9] for details on relevant proof). Therefore, next on the analysis of the droop control, DGs and SGs will be considered to be of the same features and deploy the same dynamic model.

Since inverters may connect pure storage devices, for example, batteries, to the network, can also take negative values. In that case, the storage device is charged in dependency of the excess power available in the network and thus functions as a frequency dependent load.

##### 2.3. Secondary Control for Hybrid Microgrid

The conventional secondary control is a kind of centralized control methods, the nominal value is used in the primary control stage, and it will be regulated gradually in each cycle of secondary control process. Note that the final desired frequency value can be measured from the main grid in the connected mode or obtained by command generator or many forms of DG/SG such as fuel-cells and microturbines in the islanded mode.

Different from the conventional secondary control with centralized mode, this paper in the secondary frequency control stage proposes one distributed cooperative control algorithm for each DG/SG. The basic idea of the proposed distributed control is by incorporating local communication networks that share information among neighboring units.

From this perspective, a microgrid can be considered as a multiagent system, where each DG/SG is an agent. The secondary control design for frequency resembles a tracking synchronization problem while that design for proportional power sharing is carried out. For this purpose, the important network digraph theory is introduced in the following.

##### 2.4. Graph Theory

In this paper, the required communication network for microgrids can be modeled by a digraph , where the node set denotes the set of DGs/SGs, the set of edges denotes the communication links among DGs/SGs, and is a weighted adjacency matrix defined as and . if and only if the edge . The set of neighbors of the th DG/SG is given by . The degree matrix is defined as with , and then is the Laplacian matrix, which has all row sums equal to zero; that is, , with .

In the following, the digraph is used to describe the interconnection topology of a microgrid consisting of one virtual leader DG/SG, denoted by , and followers DGs/SGs, denoted by . Diagonal matrix is called the leader adjacency matrix, where if follower DG/SG is connected to the leader across the communication link ; otherwise . The following lemma plays an important role in the proof of our main results.

Lemma 1 (see [32]). *If the digraph corresponding to the Laplacian matrix is strongly connected, then (i) if is a left eigenvector of corresponding to eigenvalue , then holds for all , and (ii) all the eigenvalues of have positive real parts, where , , and for all .*

#### 3. Main Results

In this section, one will discuss the distributed cooperative control schemes for frequency and power sharing in Sections 3.1 and 3.2, respectively.

It should be noted that the dynamics of the voltage and current control loops are much faster than the dynamics of the power control loop [30]. Therefore, it is necessary and intentional to design a typical inverter system by applying the time separation operation [6]. In order to facilitate our analysis, this paper will focus on the stability of power control loop, while the investigation of current controllers follows traditional proportional-integral (PI) controllers.

##### 3.1. Distributed Cooperative Control for Frequency

In this paper, the - reference frame transformation is considered, where the -axis and -axis of the reference frame of each inverter are rotating at the common reference frequency. Based on the traditional droop control strategy, the reference of frequency of the th inverter-based DG or SG is determined by the power control loop and can be abstracted aswhere is the nominal set point of output frequency and is the nominal set point of output active, respectively, while is the frequency droop coefficient.

In the following, a distributed pinning control algorithm is designed to pin the frequency of DGs/SGs to synchronize the reference frequency , respectively.

According to the thought of so-called input-output feedback linearisation [27], differentiating the equations in (6) yieldswhere is the pinning controllers which will be designed later such that for .

Then the nominal control input used in the droop control procedure can be computed by the following integrator:

Since it is assumed that each DG/SG only needs to communicate with its neighboring DGs/SGs through a communication network, the pinning controllers are designed based on the own information of each DG/SG and the information of its neighbors:for , where the time-varying delay is a continuously differentiable function and the factors of adjacency matrix are also allowed to be dynamically changing.

In order to facilitate analysis, some denotations are given. Let , with and ; then combining (7) and (9) and the fact that , one can deduce the following error system:where with the Laplacian matrix of digraph , denotes the Kronect product, denotes the two-dimensional identity matrix, and ( denotes the total number of all possible directed graphs) is a switching signal that determines the communication topology . If is a constant function, then the corresponding topology is fixed.

Theorem 2. *Suppose that the digraph is always strongly connected and there exists at least one such that for any switching signal . If the time-varying delays satisfy and , then the distributed controllers (9) can pin the DG/SG output frequency to the reference states asymptotically.*

*Proof. *Define a common Lyapunov-Krasovkii function for system (10) as follows:where with the positive left eigenvector corresponding to the zero eigenvalue of the Laplacian matrix . Note that since is strongly connected, one can always choose such that and for all (by Lemma 1).

Along the trajectory of system (10), we haveBy Newton-Leibniz formula,and noting that hold for any appropriate positive definite matrix , we haveConsequently,Then, a sufficient condition for iswhich leads toThen, one can deduce thatfor given . This completes the proof.

*Remark 3. *As we all know, the communication latency is a key factor that cannot be ignored in the smart microgrid control problem since many cases have proven that neglect of the communication latency will lead to serious and unpredictable consequences. When analyzing the nonlinear systems containing time-varying communication latency, most existing literatures presented some delay independent stability conditions or just require that the delay is sufficiently small. However, in order to investigate the upper bound of the amount of time-varying delay, this paper constructs a Lyapunov-Krasovskii functional (i.e., ) with delay-related parameters. It is worthy pointing out that, for the presented nonlinear system with directed and switching topology, the asymmetry of the adjacency matrix and the nonlinearity of the controlled systems will increase the difficulty to construct and analyze its stability. By using the special matrix theory, inequality theory, and theory of delay differential equations, this paper solves the above problems.

##### 3.2. Active Power Sharing

It is [26] that has proposed a method on how to choose the frequency droop coefficients and nominal setpoints to keep DGs/SGs realizing the active power sharing according to their power ratings. This control goal is analyzed in many literatures and spread to actual applications [5]. The proposed proportional power sharing method also can guarantee that storage devices in charging mode, that is, for some , are charged proportionally.

In the following, the DGs/SGs proportionally share their power in terms of their power ratings if the gains and setpoints are chosen such that, for all , the condition is standard in the microgrid research:where and are the instantaneous maximum capacity of and and all the and are assumed to have the same sign. The active powers of DGs/SGs are able to be controlled according to condition (19) and droop control (2) the active power sharing proportionally. Condition (19) that also has a guiding significance achieves the reactive power sharing.

#### 4. Simulation Results and Discussion

In this section, the effectiveness of the proposed distributed control algorithms will be verified by simulating an islanded lossy microgrid in MATLAB. Based on the above analysis for the identity property on dynamics of DGs and SGs, without loss of generality, five DGs have been chosen to the test system. Figure 4 shows the basic diagram of the test system which consists of five DGs. The specifications of the DGs, lines, and loads are summarized in Table 1. Note that its communication topology can be abstracted as digraph shown in Figure 5. The associated Laplacian matrix of isand the corresponding leader adjacency matrix can be chosen as .

In the following simulation, the time-varying delay is always taken as ; then and . For the fixed digraph , by simple calculation, one can choose and then yield , , , and with . Therefore, the delay constraint condition in Theorem 2 can be verified as

For the switching digraph shown in Figure 5, one can verify the above conditions similarly.

*(1) Case 1 (fixed topology).* Consider a directed communication network with five DGs as shown in . Figures 6 and 7 show the state evolution processes of frequency, corresponding to active power with and without delay, respectively. It is assumed that the demands of loads increase suddenly the islanded operation microgrid at s. As seen in Figure 6(a), after islanding, each DG’s frequency goes to different values less than reference value rad/s. Then the primary control and secondary control are applied at s and s, respectively. Due to the droop control technique, each DG’s frequency converges to stability value still less than reference value. However, the secondary control restores the operating frequency to synchronize their reference value after 1.6 s.

**(a)**

**(b)**

**(a)**

**(b)**

Comparing Figures 6 and 7, respectively, one finds that the existence of communication delay postpones the convergence speed of frequency.

*(2) Case 2 (switching topology).* Consider a directed network with switching topology as shown in Figure 5. The controlled system begins at the state and switches at to the next stage with the following switching order: , where the switching signal . All the weights are set to be , and the corresponding leader adjacency matrix for to is chosen as , , , and . The associated system parameters are similar to the fixed topology case.

Figures 8 and 9 show the state evolution processes of frequency, corresponding to active power with and without delay, respectively. Comparing Figures 6-7 and Figures 8-9, it is not difficult to find that the evolution curves in Figures 8-9 contain some slight oscillation phases. It is because that the controlled system in this case contains frequent switching signals. Similarly, the communication delay also postpones the convergence speed of the whole switching system.

**(a)**

**(b)**

**(a)**

**(b)**

However, if we take the large delay which does not satisfy the constraint conditions in Theorem 2, the proposed distributed control algorithms may be ineffective. As shown in Figure 10, the state evolution curves do not converge to the desired values in switching case.

**(a)**

**(b)**

#### 5. Conclusion

In this paper, a distributed cooperative control scheme for hybrid lossy microgrids has been proposed. In this method, a distributed secondary control encompasses each DG/SG local controller and the communication system; that is, each DG/SG only needs to communicate with its neighboring DGs/SGs intermittently even though the communication delays may be time-varying. The designed distributed control is able to synchronize the frequency of hybrid lossy microgrids to the desired state (50 Hz) in a distribution switching network under the existence of time-varying communication delays. All the frequency of DGs/SGs can be synchronized to the desired value in both fixed and switching distributed networks. Meanwhile, numerical simulations show higher robustness in front large communication latency.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

This work was supported by National Key Technology Support Program of China under Grant no. 2013BAA01B01; by the National Natural Science Foundation of China under Grant nos. 61403133, 61374064, and 61304152; by the Young Teachers Growth Plan of Hunan University under Grant no. 531107040651; by the Postdoctoral Foundation of China under Grant no. 2013M540627; by the Postdoctoral Special Foundation of Hunan Province under Grant no. 2013RS4042; by the Natural Science Foundation of Hunan Province under Grant no. 14JJ3051; by the Doctoral Fund of Ministry of Education of China under Grant no. 20130161120016; by the Fundamental Research Funds for the Central Universities no. 2014208020201.