Hierarchical Two-Layer Distributed Control Architecture for Voltage Regulation in Multiple Microgrids in the Presence of Time-Varying Delays

: The Multiple Microgrids (MMGs) concept has been identiﬁed as a promising solution for the management of large-scale power grids in order to maximize the use of widespread renewable energies sources. However, its deployment in realistic operation scenarios is still an open issue due to the presence of non-ideal and unreliable communication systems that allow each component within the power network to share information about its state. Indeed, due to technological constraints, multiple time-varying communication delays consistently appear during data acquisition and the transmission process and their effects must be considered in the control design phase. To this aim, this paper addresses the voltage regulation control problem for MMGs systems in the presence of time-varying communication delays. To solve this problem, we propose a novel hierarchical two-layer distributed control architecture that accounts for the presence of communication latencies in the information exchange. More speciﬁcally, the upper control layer aims at guaranteeing a proper and economical reactive power dispatch among MMGs, while the lower control layer aims at ensuring voltage regulation of all electrical buses within each MG to the desired voltage set-point. By leveraging a proper Driver Generator Nodes Selection Algorithm, we ﬁrst provide the best choice of generator nodes which, considering the upper layer control goal, speeds up the voltage synchronization process of all the buses within each MG to the voltage set-point computed by the upper-control layer. Then, the lower control layer, on the basis of this desired voltage value, drives the reactive power capability of each smart device within each MG and compensates for possible voltage deviations. Simulation analysis is carried out on the realistic case study of an MMGs system consisting of two identical IEEE 14-bus test systems and the numerical results disclose the effectiveness of the proposed control strategy, as well as its robustness with respect to load ﬂuctuations. The proposed control architecture, based on cluster-oriented distributed cooperative control strategies, is able to guarantee an optimal energy management of the overall MMGs system despite the presence of heterogeneous time-varying communication delays. The numerical analysis was performed on the realistic case study of an MMGs system consisting of two IEEE 14-bus test systems. Simulation results, both in nominal and load ﬂuctuations scenarios, have disclosed the effectiveness and robustness of the proposed control architecture.


Introduction
The increasing use of Renewable Energy Sources (RESs), Distributed Generators (DGs) and storage systems is imperative for the introduction of Microgrids (MGs) in actual power distribution networks because of the many benefits it could lead to in terms of reductions in power losses and electrical system performance improvements. In general, an MG is a small distribution power system have been proposed for the coordination of the MMGs system [25], namely master-slave control [26], peer-to-peer control [27] and hierarchical control [28]. Leveraging this latter framework, a hierarchical decentralized SoS architecture for the energy management of an MMGs system has been proposed in [15] by formulating a bi-level optimization problem and considering the RESs uncertainty, while [29] exploits a smart transformer-based approach. Again, a hierarchical two-layer control architecture for a cluster of islanded MGs sharing information via intermittent communications has been proposed in [22,30].
Although the solutions proposed herein allow for the coordination control problem for MMGs to be effectively and reliably solved, their design leverages the restrictive assumption of an ideal and always reliable communication among the different electrical components constituting the whole power grid system. However, when considering a wireless communication network, based on, for example the IEEE 802.11 protocol, each communication link that connects a pair of agents is affected by a different variable time-delay whose value depends on actual conditions, or possible impairments, of the communication channel [11]. It follows that the hypothesis commonly made in the technical literature of an ideal communication network in the MMGs may be unrealistic. Hence, the control architecture deployment in realistic operation scenarios is still an open problem, which warrants further investigation [17,22].
To overcome this limitation and solve the problem of determining the optimal power dispatch in practical MMGs systems, we propose a novel cooperative cluster-oriented hierarchical control architecture that accounts for the presence of communication time-varying delays. Exploiting the Multi-Agent System (MAS) mathematical framework, the proposed control architecture is based on a network of cooperative smart controllers (called nodes), each one regulating the voltage magnitude of a specific bus within each MG cluster constituting the whole MMGs system. All nodes/controllers are able to share information about their states with their neighbors (within their communication range) so that the control actions can be cooperatively computed by embedding within the online decision-making process not only information coming from local sensing, but the delayed network information about the surroundings. More specifically, the proposed control solution consists of two cyber layers, namely the upper and lower control layers. The upper control layer aims at guaranteeing a suitable and economical reactive power dispatch among multiple MGs. This can be done by controlling only a specific subset of the DGs within each MG, i.e., the driver generator nodes, via a cooperative control action that, based on the knowledge of the rated voltage value to be imposed to the whole MMGs and the neighboring delayed information shared with other drivers belonging to other MGs through the inter-cluster communication network, determines the voltage reference value for each MG. The driver generator nodes set, as well as the inter-cluster communication network topology, is computed by leveraging the Master Stability Function (MSF) approach [31,32] which determines the best choice of these nodes for speeding up the voltage synchronization process of all buses within each MG to the safe voltage set-point computed by the upper-control layer. Conversely, based on this latter point, the lower control layer aims at ensuring voltage regulation of all PV/PQ buses within each MG to the desired voltage set-point. This is done via a cooperative control action that, by exploiting delayed state information shared among the smart devices within the single MG via the intra-cluster communication network, drives the reactive power generation capability of each smart device and compensates for possible voltage deviations.
To prove the effectiveness of the proposed cluster-oriented cooperative control strategy, we consider an MMGs system consisting of two identical MGs, i.e., the IEEE 14-bus test system. The results confirm the ability of the proposed solution to guarantee the two-fold voltage regulation of the MMGs system, despite the presence of communication time-varying delays and/or possible load variations.
The paper is organized as follows. Section 2 describes the voltage control problem for MMGs systems, as well as the dynamics of the cooperative smart controllers network that was exploited for its control. The Cluster-Oriented Cooperative control strategy, together with Driver Generator Nodes Selection Algorithm based on MSF formalism, is detailed in Section 3. Section 4 discloses the effectiveness and robustness of the proposed approach for the exemplar MMGs system consisting of two IEEE 14-bus test system. Finally, conclusions are drawn in Section 5.

Multiple MicroGrids Modeling
Consider an MMGs system consisting of M MGs, labeled as MG 1 , . . . , MG k , . . . , MG M as in Figure 1. Each power grid MG k (∀k ∈ I M ) consists of N k c capacitor banks, N k g distributed generators and N k = N k g + N k c buses managed by N k cooperative smart controllers. Leveraging the Cluster-Oriented Cooperative Control strategy [28] and Distributed Hierarchical Cooperative (DHC) framework [22], each power grid MG k is able to share information via wireless inter-cluster communication networks with its neighbors MG k , being k ∈ I M with k = k . Moreover, each smart device j within the single MG k (being j ∈ I N k ) can communicate via an intra-cluster wireless communication network with its neighboring smart devices q, being q ∈ I N k with q = j. The intra-cluster communication network structure corresponds to the information exchange topology among the N k smart controllers within the MG k . The inter-cluster communication, instead, is only enabled for n k smart devices within each MG k , associated to some specific generators buses and named drivers/pinners. This inter-cluster communication network structure, indicated withG, defines the information flow shared among the different l MMGs pinners, being l ∈ I pin .
This decomposition of the MMGs into two cyber layers, i.e., the upper and the lower ones, enables a twofold voltage regulation control, as well as an optimal power management of the entire MMGs system [22]. In particular, supposing that the pinners of each MG can access the voltage-rated value V rated to be imposed to the whole MMGs system, the upper inter-cluster layer aims at guaranteeing a correct and economical reactive power dispatch among multiple MGs by computing the reference voltage set value, i.e., V nom k , to be sent to each pinner s ∈ I pin n k (being k ∈ I M ). Conversely, the lower intra-cluster layer aims at ensuring a faster voltage synchronization process of all PV/PQ buses within the MG k to the voltage set-point V nom k by controlling the reactive power generation capability of each MG k bus to compensate for voltage deviations. Note that, for each MG k , the optimal identification of the pinner generator nodes, i.e., n k , as well as their best positions in the grid are identified according to the MSF formalism within Pinning Control Theory, which proves that the most influential nodes are highly dependent on the network structure [31,33,34].

Double-Layer Communication Network
The double-layer communication network describing the information exchange among the smart controllers within the MMGs in the cyber-physical space can be modeled according to graph theory.
Regarding the intra-cluster network for each MG k (k ∈ I M ), the communication topology of the N k c cooperative smart device that controls the capacitor bank buses can be described as a directed graph G N k The associated adjacency matrix with non-negative elements is being p, ρ ∈ I N k c with p = ρ. We assume α k pρ = 1 in the presence of a communication link from the device p to device ρ, otherwise α k pρ = 0. Moreover, α k pp = 0, i.e., self-edges (p, p) k are not allowed. The presence/absence of connections among the N k c cooperative smart controller and the ν k = N k g − n k smart controller for the non-pinner generation buses is instead described by the overall with the n k pinners, we introduce the leader-adjacency matrix B k = diag{b k js } ∈ R k × k , whose value is b k js = 1 if the j-th bus controller receives information from s-th driver generator bus (being s ∈ I pin n k ), otherwise b k js = 0. Similarly, the upper inter-cluster cyber-network is described via the digraphG = {Ṽ,Ẽ ,Ã} with virtual node setṼ = {V pin 1 , . . . , V pin M } being V pin k the set of driver generator nodes for the k-th MG with cardinality card(V pin k ) = n k < N k g .Ẽ ⊆Ṽ ×Ṽ represents the set of cyber communication links, while the adjacency matrixÃ = [a sl ] ∈ R n×n (with n = ∑ M k=1 n k ) is such that a sl = 1 if there is a link from pinner generator s ∈ I pin to pinner generator l ∈ I pin , with s = l.
Finally note that, since, in practice, communication networks are commonly affected by latency in the information delivery due to the current conditions of the communication infrastructure, an unknown heterogeneous time-varying delay τ qj (t) (being q and j generic electrical nodes) can be associated to each direct edge. Although delays are time-varying, they are usually bounded during the normal operating conditions of technological communication networks, hence τ qj (t)

Cooperative Smart Agents Dynamics
The dynamics of each smart device p within the k-th MG and associated to the p-th capacitor bank bus is described ∀p ∈ I N k c and ∀k ∈ I M by the following dynamical system [11]: where Q k,p [p.u.] is the reactive power of the p-th capacitor bank within the k-th MG; u k,p (t, τ k pj (t), τ k ps (t)) is the distributed cooperative control protocol that, by exploiting the local measurement and the ones shared via the intra-cluster communication network, need to be designed in order to drive the reactive power of the capacitor bank bus, as well as its voltage magnitude; , p = j) and τ k ps (t) (s ∈ I pin n k ) are the heterogeneous time-varying communication delays affecting the communication among each pair of electrical nodes. Their actual values, at a given time-instant, depend on the current availability of the wireless channels themselves.
Note that, since the capacitor bank is composed of PQ buses and its voltage magnitude can be regulated by imposing a proper reactive power variation, (1) well represents this kind of phenomena and allows for the reactive power of the capacitor bank to be adapted to guarantee voltage regulation (see [11,36] and references therein).
The dynamics of the i-th smart device controlling the non-pinner generators nodes within MG k (∀k ∈ I M ) is described ∀i ∈ I non−pin k as [11]: where V k,i [p.u.] represents the voltage magnitude of the i-th non-pinner generator bus within the k-th is the control action to be designed, exploiting both the local measurement and the ones shared via the intra-cluster communication network in order to synchronize the voltage of the i-th non pinner generator bus to the desired nominal value V nom k,s ; τ k ij (t) is the heterogeneous time-varying delay affecting the information shared between the non-pinner DG i (being i ∈ I non−pin k ) and the smart device j (being j ∈ I lower k ), while τ k is (t) is the communication delay affecting the exchanging data between non-pinner i-DG and s-th DG pinner (s ∈ I pin n k ) within the k-th MG. The behaviour of each smart device controlling the driver generator bus s (∀s ∈ I pin n k ) within MG k (∀k ∈ I M ) is described by the following dynamical system: where V nom k,s [p.u.] is the voltage magnitude of the s-th driver generator within the k-th MG; V rated [p.u.] is the known voltage magnitude set-point to be imposed to the whole MMGs system; u k,s (t, V rated , τ ls (t)) is the cooperative control action which, leveraging both local measurements and the information shared with the other pinners l (∀l ∈ I pin ) via the inter-cluster communication networkG, computes the voltage nominal reference value V nom k,s = V nom k to be imposed to the k-th MG based on both the knowledge of V rated and the reactive power mismatch among the MGs; τ ls (t) is the heterogeneous time-varying delay affecting the data exchange between the pinner generator node l within MG k and the pinner s within MG k .
We remark that the s-th pinner generator node, ∀s ∈ I pin n k , acts as a leader for the whole MG k in the lower control layer by imposing the reference voltage magnitude V nom k for both the capacitor bank and non-pinner generator nodes.

Design of Cluster-Oriented Cooperative Control Strategy
In what follows, we address the two-layer voltage regulation control problem and the related power management of the MMGs system. The design of the cooperative controllers starts with the appropriate selection of the n k driver generator nodes within each MG k , ∀k ∈ I M . The pinner selection algorithm leverages the MSF approach [31] and guarantees the choice of the optimal nodes for speeding up the voltage synchronization process of all buses within each MG to a safe voltage set-point via the double-layer cooperative control strategy. Once the driver nodes set I pin n k is defined for each MG k , the upper cooperative control layer, leveraging the inter-cluster cyber communication topology and the knowledge of V rated , computes the nominal reference value V nom k to be sent to each MG k . On the basis of this latter point, the lower cooperative control layer, exploiting the intra-cluster cyber communication topology, guarantees that each electrical node within the MG k reaches this reference voltage value despite the presence of communication delays.

Driver Generator Nodes Selection
Here, we present the selection algorithm based on MSF that we leverage to find the best driver generator nodes [31,32]. It provides a suitable criteria to find the most influential nodes set in each MG k , k ∈ I M , i.e., I pin n k , and hence, the overall set I pin = M k=1 I pin n k involved in the upper inter-cluster control layer, needed to perform an accurate economical power allocation within the whole MMGs system. Specifically, for each MG k , the algorithm selects within the set I N k g , the n k < N k g driver generator nodes. This is performed through the evaluation of the matrix C structure. The latter is obtained by considering all the possible configurations for the Laplacian matrix of the intra-cluster communication network, i.e., L = [l ij ] = D − A, with A being the corresponding adjacency matrix and D the in-dregree matrix, augmented with the driver generator nodes on the main diagonal as where β i = 1 if the generic node i is a driver generator node, 0 otherwise, ∀i = 1, · · · , N k . Iteratively, for each possible matrix C, the algorithm computes the eigenratio R, defined as with λ N and λ 1 being the biggest and smallest eigenvalues of the matrix C, respectively. Then, among all the possible driver node configurations, the algorithm selects the ones such that R <R, with R being a positive threshold that guarantees that the Master Stability Function is negative [32]. Finally, the algorithm solves an optimization problem whose solution returns the best driver node configuration, indicated as C with the minimum eigenratio R . This optimization problem can, hence, be formulated as where β = [β 1 , β 2 , . . . , β N k ] is the optimization variable representing the number and the position of the driver nodes as well as n K the total number of them. The selection algorithm for the driver generator nodes set I pin n k for each MG k is summarized in Algorithm 1.

Inter-Cluster Cooperative Control Strategy Design
The upper inter-cluster cooperative control layer aims at adjusting the power flow among the M MGs to achieve the following output reactive power balance: where Q k (t) [p.u.] and Q k (t) [p.u.] are the total reactive power outputs of MG k and MG k , respectively, ∀k, k ∈ I M with k = k . Considering that the inter-cluster communication is enabled only for the n driver generator nodes within the MMGs system, the control goal as in (7) can be recast as the following cooperative consensus problem: where Q k ,l (t) [p.u.] and Q k,s (t) [p.u.] are the reactive power outputs of pinner generator nodes l and s, respectively, with l ∈ I pin n k , s ∈ I pin n k being ∀k, k ∈ I M k = k . Since each pinner generator node s within MG k (i.e., s ∈ I pin n k , k ∈ I M ) knows the voltage rated value to be imposed on the whole MMGs system, i.e., V rated ), in order to achieve (8), we design the cooperative control action for the pinner s within each MG k in (3), as u k,s (t, V rated , τ ls (t)) = V rated + k Q ∑ l∈I pin a sl Q k ,l (t − τ ls (t)) − Q k,s (t − τ ls (t)) (9) where k Q is the control gain; a sl models the upper inter-cluster cyber-network emerging from the presence/absence of communication link between the s-th and the l-th pinner generator of MG k and MG k , respectively; τ ls (t) is the heterogeneous time-varying delay depending on the specific communication link, with s ∈ I n k and l ∈ I n k . Indeed, since smart devices share information through a wireless communication network, there can be a time delay in the communication among devices as a result of the non-ideality of communication channels and the limited bandwidth. The above-mentioned considerations lead to the need of running the controller based on outdated information and of using the time stamp inserted into messages to correctly correlate the information. Finally, note that the cooperative control strategy in (9) also generates the voltage nominal reference value V nom k,s = V nom k to be imposed on all the PV/PQ buses within each k-MG via the lower intra-cluster control layer.
is the optimal configuration of the n k driver generator nodes.

Intra-Cluster Cooperative Control Strategy Design
The intra-cluster voltage regulation control problem for each MG k requires the achievement of the following control goals ∀k ∈ I M : 1. To design a fully distributed cooperative controller u k,i (t, τ k ij (t), τ k is (t)) in (2), which, based on local measurements and outdated networked information, regulates the voltage magnitude V k,i [p.u.] of the i-th non-pinner generator bus within MG k (i ∈ I non−pin k ) in order to reach and maintain the desired reference voltage value V nom k,s [p.u.] (∀s ∈ I n k ) computed by the upper control level, i.e.,: where V k,j [p.u.] is the voltage magnitude of the generic bus j ∈ I lower k with i = j and V nom k,s = V nom k [p.u.] is the nominal voltage magnitude of the s-th driver generator node, ∀s ∈ I pin n k . 2. To design a fully distributed cooperative control protocol u k,p (t, τ k pj (t), τ k ps (t)) in (1), which opportunely drives the reactive power of the capacitor bank p (p ∈ I N k c ) by updating its voltage magnitude V k,p [p.u.] until it reaches the desired reference behavior as imposed by the N k g generators within the MG k ; i.e: where V k,p [p.u.] is the voltage magnitude of the generic bus p ∈ I lower k with p = j and V nom k,s = V nom k [p.u.] is the nominal voltage magnitude of the s-th driver generator node, ∀s ∈ I pin n k .
To achieve control objectives (10) and (11), for each non-pinner generator node i ∈ I non−pin k , we design the following delayed consensus-based control strategy which weights both the voltage error with respect to the desired nominal value V nom k,s of the s-th driver generator node and the local voltage networked information as where a k ij model the presence/absence of a communication link between the bus j and the bus i; b k is is the leader-adjacency matrix element (i.e., B k , see Section 2.1); τ k ij (t) and τ k is (t) are the communication latencies arising from the information exhange among the smart controllers; γ 1 and γ 2 ∈ R are the control gains to be tuned.
Conversely, to fulfill the control objectives in (12) and (13), we propose for each capacitor bank p ∈ I N k c , the following delayed consensus-based controller: where a k pj models the communication topology emerging from the presence/absence of a communication link between the p-th and j-th bus, while b k ps models the presence/absence of communication among bus p and pinner bus s; τ k pj (t) and τ k ps (t) are the heterogeneous time-varying delays; δ 1 and δ 2 ∈ R are the control gains to be tuned in order to maintain the reactive power of the bus p within its operating range [Q min k,p , Q max k,p ].

Case Study
In this section, we prove the effectiveness of the proposed cluster-oriented cooperative control strategy (based on the Driver Generator nodes Selection Algorithm in Section 3.1) in guaranteeing the voltage regulation of an MMGs system, despite the presence of communication time-varying delays and/or possible load variations. More specifically, we consider an MMGs system consisting of k = 2 identical IEEE 14-bus test systems MGs sharing information via a double-layer communication network, as in Figure 2. Each MG consists of N g 1 = N g 2 = 5 distributed generators, i.e., buses 1, 2, 3, 6, 8, and N c 1 = N c 2 = 9 capacitor banks with twenty power transmission lines, respectively. Details about line impedances, load and reactive power limits were provided in [37], while initial conditions for the N k electrical bus within each MG k (k = 1, 2) were reported in Table 1, where the rated voltage value V rated to be imposed to the whole MMGs system was also reported. Note that, the initial conditions for each electrical bus within MG 1 and MG 2 are randomly chosen according to their acceptable operative ranges.
Without loss of generality, we assume that, within each MG k , the intra-cluster communication topology perfectly matches the physical one described by the power transmission lines. Indeed, since the electrical topology satisfies the global reachability property of generator nodes, the above assumption can be made. However, we remark that any other communication topology meeting the global reachability property of generator nodes can be adopted to prove the effectiveness and the robustness of our control strategies [11]. The inter-cluster communication topologyG is instead determined via Algorithm 1 to guarantee a faster synchronization process of all the electrical buses to the voltage-rated value V rated . Namely, by setting n k = 3 ∀k = 1, 2, the algorithm returns, for each MG k , the optimal configuration for driver generator nodes corresponding to generator nodes 2, 6 and 8 with a communication topology described by the following Laplacian matrix: The multiple time-varying delays, affecting both the inter-cluster and the intra-cluster communication networks, have been emulated as random variables with uniform distribution within the range [0, τ ], where τ = 0.1 [s] is the upper bound setting above the typical average end-to-end communication delay in a wireless network [11,38,39].
Moreover, the control gains parameters for the cooperative control strategies (9) (14) (15) are selected to avoid the reactive power exceeding its maximum/minimum allowable values. Their values are listed in Table 1.

Initial Conditions for the Multiple MG Cluster System
Voltage magnitude of generator bus ∀k = 1, 2 [p.u.]

Nominal Scenario
In this scenario, where no load variations are considered, the two MGs are electrically disconnected from each other at t = 0 [s] and, then, connected at t = 5 [s]. Therefore, the lower controller and the upper controllers (as well as its corresponding pinning links) are disabled at t = 0 [s] and switch-on at The results in Figure 3 confirm the effectiveness of the proposed control approach in guaranteeing the cluster-oriented double-layer cooperative control despite the presence of heterogeneous time-varying communication delays. Specifically, since upper-layer controllers u k,s ∀s ∈ I pin n k , ∀k ∈ I M in (9) are inactive during the time interval [0, 5] [s], generator voltage magnitudes do not reach the nominal value V nom k = 1.03 [p.u.], k = 1, 2 (Figure 3a,b). When, at t = 5 [s] the cluster-oriented cooperative control strategy is enabled, the voltage magnitudes of all the generator buses (including the driver generator nodes) within both the MG 1 and the MG 2 converge to the desired nominal values V nom k = 1.03 [p.u.], k = 1, 2 (see Figure 3 (a) and (b), respectively). As a result, the control objectives in (7) and (8) are fulfilled. As a consequence, the lower cooperative control protocol u k,p ∀p ∈ N c k , k = 1, 2 in (15) drives the reactive power Q k,p [p.u.] of each capacitor bank p to guarantee that the corresponding voltage magnitudes V k,p reach the nominal value V nom k = 1.03, [p.u.], as shown in Figure 3c-e. Figure 4 highlights the time histories of reactive powers of capacitor banks within MG 1 and MG 2 . Thus, the control goals in (10) and (11), as well as the control objectives in (12) and (13)

Load Fluctuations Scenario
In this section, we present the results of the performance analysis of our cooperative controllers, confirming voltage regulation of the overall MMGs when load fluctuations occur. Indeed, since in a real practical situation, the load demand is sensitive to frequent changes, it is particularly desirable to evaluate the robustness of the proposed control strategies in this crucial scenario [1].
The appraised load fluctuations profile L(t) is the one depicted in Figure 5a, where a maximum load variation of ±50% can be observed.  Figure 5b,c). Accordingly, the voltage magnitudes of all buses within the MMGs system achieve the desired set point (see Figure 5d), thus proving the ability of the designed controllers to react to load fluctuations by inducing a reactive power variation ( see Figure 6a-c) and promptly restore the voltage to the required level. Finally, Figure 5e,f show the time histories of the voltages of each bus p ∈ N c k k = 1, 2 in the MG 1 and the MG 2 , respectively. In conclusion, all the control goals (10)- (13) are also fulfilled in this simulation scenario.

Comparison Analysis
Here, we compare the performances of our fully delayed control strategies in (9), (14) and (15), which were designed according to delayed information processing with respect to the ones achievable with the following control strategy: which only involves the transmission delays of the information received from neighbors. Therefore, each controller in (17a)-(17c) consists of two different terms-a local action depending on the state variables of the generic node itself (measured onboard) without any delay, and an action depending on the information received from neighboring nodes (distributed generators and/or capacitor banks) through the communication network that was affected by a time-varying communication delay. Exemplar results, reported in Figure 7, refer to the nominal scenario detailed in Section 4.1. As it is possible to observe, the proposed fully delayed Cluster-Oriented Double-Layer Cooperative Control approach achieves better leader tracking performance in terms of voltages and reactive power with respect to the results obtained via controllers in (17a)-(17c) (depicted in Figure 7). This is due to the benefit introduced by networked-induced delays on the overall MMGs system [40,41] which improves the closed-loop performances. Hence, the advantages obtained by considering random delays and their effects on the closed-loop network in the control design phase have been shown. Therefore, by implementing our distributed control strategies via fully outdated information, any possible instability sources can be counteracted [42,43].

Conclusions
In this paper, the voltage regulation problem, as well as the optimal economical power sharing, for a multiple MGs cluster system in the presence of heterogeneous time-varying communication latency has been investigated and solved via a two-layer distributed control approach, which also leverages the MSF method. The proposed control architecture, based on cluster-oriented distributed cooperative control strategies, is able to guarantee an optimal energy management of the overall MMGs system despite the presence of heterogeneous time-varying communication delays. The numerical analysis was performed on the realistic case study of an MMGs system consisting of two IEEE 14-bus test systems. Simulation results, both in nominal and load fluctuations scenarios, have disclosed the effectiveness and robustness of the proposed control architecture. The set of k nodes within k-MG, involved into the lower cyber layer n The number of total driver generator nodes within MMGs I pin = M k=1 I n k The set of all n driver generator nodes within the MMGs G k The intra-cluster cyber network of the k-th MG (except for driver generator nodes) G The inter-cluster cyber network Q k Output Reactive power of the k-th MG [Q min k,j , Q max k,j ] Reactive Power Operating Range of the j-th electrical bus within k-th MG Q k,p Reactive power of the p-th capacitor bank within k-th MG Q k,s Reactive power of the s-th pinner generator within k-th MG V k,p Voltage magnitude of the p-th capacitor bank within k-th MG V k,i Voltage magnitude of the i-th non pinner generator node within k-th MG V nom k,s Voltage magnitude of the s-th pinner generator node within k-th MG V rate Rated voltage magnitude imposed on the MMGs V nom k Nominal Voltage magnitude imposed on the k-th MG