# Edge Computing- andH∞-Switching-Based Networked Control for Frequency Control in Multi-Microgrids with Time Delays

Authors: Peng Yang; Wei Guo; Guanghua Wu; Cong Wang; Kai Zhang; Ran Zhang
Journal: Wireless Communications and Mobile Computing
 (2021)

---

## Abstract

The frequency stability of multi-microgrids is easily affected by random load fluctuations and intermittent renewable resources. Additionally, geographically distributed generation equipment usually cannot adopt the “point-to-point” dedicated communication scheme to realize the information exchange considering the construction and computation costs. Therefore, aH∞-switching frequency control strategy for multi-microgrids based on edge computing framework is proposed in this paper. Firstly, an edge computing device is set up in each microgrid to collect the operation statuses of local participating equipment and generate the control instructions to ensure the real-time local frequency stability. Secondly, the multihop data transmission process in edge computing environment is described as a cascade queuing model. Then, the frequency control system in each microgrid is described as a switching model dependent on the varying time delays. Finally, via constructing a Lyapunov function, the constraints of the controller gains ensuring the H∞-damping performance for external load demands and the renewable outputs are derived at the same time. Simulation results show that compared with the traditional centralized control schemes, the peak value of our proposed edge computing framework is reduced by 32.51% compared with the traditional centralized control scheme. Moreover, under the same edge computing framework, the integral of absolute error (IAE) of frequency with the proposed H∞ control strategy can be reduced by 37.19% at least. Therefore, a better transient performance can be obtained with our proposed method.


---

## Body

## 1. Introduction

The microgrid is usually used to supply power for the rural areas, islands, and so on. Integrating multi-microgrids into the power systems can effectively address the contradiction between the load growth and the power infrastructure expansion [1]. The main feature of the multi-microgrids is that each microgrid integrates the high-proportional renewable generators and exchanges active power with the neighbouring microgrids through tie lines [2, 3]. However, the uncertainties in renewable outputs and load demands usually cause the active power mismatch between the supply and demand in the power systems. Correspondingly, the frequency in each microgrid will deviate from the rated value [4]. To maintain the real-time balance of active power supply and demand at the rated frequency point, the global operation statuses of each microgrid are usually collected in traditional centralized control strategies while the output adjustment instructions are sent out to the participating equipment including synchronous generators and energy storage systems in each microgrid [5, 6]. However, because there is only one centralized control center in the whole power system, the computing burden in the control center is usually heavy. Besides, if the number of interconnected microgrids is too large, there will be significant disadvantages of the centralized control strategies in reliability and construction costs of communication facilities [7]. In addition, too many interconnected microgrids will cause the dimensional disaster in traditional centralized frequency control schemes, which may lead to no solution for the controller gains [8].With the rapid development of Internet of things and edge computing technologies, using decentralized control structure and sharing communication networks have attracted considerable attentions. Different from the traditional centralized control scheme, in edge computing framework, the edge computing device with data storing, processing, and analysing capabilities is set up in each microgrid which is considered as the local control center to realize the local power balance at the rated frequency point [9, 10]. Since the edge computing device is closer to the participating equipment, the corresponding transmission delay can be shorten compared with the traditional centralized control schemes. However, it still should be noted that the influences of stochastic time delays and packet losses cannot be ignored in the controller design process. The numerical results in [11, 12] show that even the millisecond-level time delay may cause the frequency deviation exceeding the allowable ranges or even instability risk. Hence, the literatures concerning the delay/packet loss-dependent frequency control in multi-microgrids can be divided into the following two categories.


### 1.1. Delay Margin Calculation

In this context, the controller parameters for frequency stabilization are firstly determined. Then, the maximum allowable time delay that guarantees frequency deviation within the permitted ranges is calculated. For example, the analytical relationship between delay margin and controller gain is derived by constructing a Lyapunov function in [13]. Similarly, in [14], a stability criterion of frequency control system concerning the controller parameters and delay margin is proposed on the basis of the regular polynomial method. In [15], an event-triggered communication mechanism is proposed for the frequency control in power systems. Only when the frequency deviation exceeds the preset threshold, the update of control instructions can be triggered. The authors in [15] also strictly proved the analytical relationship among controller parameters, delay margin, and the triggering threshold. However, the disadvantage of above studies is that the controller parameters must be given in advance. Therefore, the optimal dynamic performance of frequency control systems cannot be guaranteed when the power system suffers from the varying time delays. In addition, when the actual transmission delays of the packets exceed the allowable delay margin, the frequency control system will be instable.


### 1.2. Delay/Packet Loss-Dependent Controller Design

In this context, the controller parameters are obtained according to the actual transmission delays. Therefore, the control performance of the frequency control system can be improved effectively. For example, in [16], queuing theory is adopted to calculate the delay ranges in power systems. Besides, a decentralized control strategy based on linear matrix inequality-linear quadratic regulator (LMI-LQR) is proposed to guarantee the closed-loop asymptotic stability in the maximum delay case. The effects of packet losses on the control performance have not been discussed in [16]. In [17], a robust model predictive control (MPC) method is proposed to tackle the frequency stability problem in power systems with stochastic time delays. However, the MPC method needs to store continuous historical data to generate the control instructions. Moreover, the instruction calculation in MPC also introduces additional delays. In [18], a decentralized robust sliding control strategy is proposed to damping the frequency deviations in power systems with time delays. However, the external active power distributions from load demands and renewable outputs are assumed to be known in advance during the sliding surface construction process. In fact, it is the uncertainties of load demand and renewable outputs that lead to the frequency deviation from the rated value, so the practicability of the method in [18] is indeed open to debate. Furthermore, the current studies mainly pay close attention to the closed-loop asymptotic stability in the maximum transmission delay case. The dynamic performance during frequency restoration process is sacrificed to some extent.Based on the above analysis, aH∞-switching control strategy for frequency stability in multi-microgrids based on edge computing framework is proposed in this paper. The main contributions are given as follows:
(1)
An edge computing framework for frequency stability in multi-microgrids is proposed. Rather than the traditional centralized control schemes which require global operation statuses of all the microgrids, the edge computing device is set up in each microgrid to be responsible for maintaining the local power balance between supply and demand sides at the rated frequency point. Correspondingly, the lighter communication burden and lower computation cost can be realized(2)
Based on the queuing theory, the analytical relationship between the transmission delay and the network parameters such as packet size, transmission rate, and transmission hops in the process of multihop data transmission process under the edge computing framework is calculated. Then, the dynamics of frequency control system in each microgrid are further described as a switching model which depends on the time-varying delay. Hence, the mapping relationship between time-varying delay and dynamic frequency deviation response can be revealed more clearly(3)
By constructing a Lyapunov function, a stability criterion of the closed-loop delay-dependent frequency control system withH∞-damping performance for external power disturbances is strictly derived. Furthermore, by aiming at minimizing the integral of absolute error (IAE) of the frequency in each microgrid, a constrained controller optimization algorithm is proposed. Then, the dynamic performance during frequency restoration process can be improvedThe rest is organized as follows. State-space model for frequency control system in each microgrid is established in Section2. In Section 3, the relationships among transmission delay and key network parameters are investigated. The dynamics of closed-loop frequency control system with time-varying delays are discussed in Section 4. In Section 5, an optimization method combined with iterative linear matrix inequality and heuristic algorithm is proposed. Simulations are discussed in Section 6. Finally, conclusions are given in Section 7.


## 1.1. Delay Margin Calculation

In this context, the controller parameters for frequency stabilization are firstly determined. Then, the maximum allowable time delay that guarantees frequency deviation within the permitted ranges is calculated. For example, the analytical relationship between delay margin and controller gain is derived by constructing a Lyapunov function in [13]. Similarly, in [14], a stability criterion of frequency control system concerning the controller parameters and delay margin is proposed on the basis of the regular polynomial method. In [15], an event-triggered communication mechanism is proposed for the frequency control in power systems. Only when the frequency deviation exceeds the preset threshold, the update of control instructions can be triggered. The authors in [15] also strictly proved the analytical relationship among controller parameters, delay margin, and the triggering threshold. However, the disadvantage of above studies is that the controller parameters must be given in advance. Therefore, the optimal dynamic performance of frequency control systems cannot be guaranteed when the power system suffers from the varying time delays. In addition, when the actual transmission delays of the packets exceed the allowable delay margin, the frequency control system will be instable.


## 1.2. Delay/Packet Loss-Dependent Controller Design

In this context, the controller parameters are obtained according to the actual transmission delays. Therefore, the control performance of the frequency control system can be improved effectively. For example, in [16], queuing theory is adopted to calculate the delay ranges in power systems. Besides, a decentralized control strategy based on linear matrix inequality-linear quadratic regulator (LMI-LQR) is proposed to guarantee the closed-loop asymptotic stability in the maximum delay case. The effects of packet losses on the control performance have not been discussed in [16]. In [17], a robust model predictive control (MPC) method is proposed to tackle the frequency stability problem in power systems with stochastic time delays. However, the MPC method needs to store continuous historical data to generate the control instructions. Moreover, the instruction calculation in MPC also introduces additional delays. In [18], a decentralized robust sliding control strategy is proposed to damping the frequency deviations in power systems with time delays. However, the external active power distributions from load demands and renewable outputs are assumed to be known in advance during the sliding surface construction process. In fact, it is the uncertainties of load demand and renewable outputs that lead to the frequency deviation from the rated value, so the practicability of the method in [18] is indeed open to debate. Furthermore, the current studies mainly pay close attention to the closed-loop asymptotic stability in the maximum transmission delay case. The dynamic performance during frequency restoration process is sacrificed to some extent.Based on the above analysis, aH∞-switching control strategy for frequency stability in multi-microgrids based on edge computing framework is proposed in this paper. The main contributions are given as follows:
(1)
An edge computing framework for frequency stability in multi-microgrids is proposed. Rather than the traditional centralized control schemes which require global operation statuses of all the microgrids, the edge computing device is set up in each microgrid to be responsible for maintaining the local power balance between supply and demand sides at the rated frequency point. Correspondingly, the lighter communication burden and lower computation cost can be realized(2)
Based on the queuing theory, the analytical relationship between the transmission delay and the network parameters such as packet size, transmission rate, and transmission hops in the process of multihop data transmission process under the edge computing framework is calculated. Then, the dynamics of frequency control system in each microgrid are further described as a switching model which depends on the time-varying delay. Hence, the mapping relationship between time-varying delay and dynamic frequency deviation response can be revealed more clearly(3)
By constructing a Lyapunov function, a stability criterion of the closed-loop delay-dependent frequency control system withH∞-damping performance for external power disturbances is strictly derived. Furthermore, by aiming at minimizing the integral of absolute error (IAE) of the frequency in each microgrid, a constrained controller optimization algorithm is proposed. Then, the dynamic performance during frequency restoration process can be improvedThe rest is organized as follows. State-space model for frequency control system in each microgrid is established in Section2. In Section 3, the relationships among transmission delay and key network parameters are investigated. The dynamics of closed-loop frequency control system with time-varying delays are discussed in Section 4. In Section 5, an optimization method combined with iterative linear matrix inequality and heuristic algorithm is proposed. Simulations are discussed in Section 6. Finally, conclusions are given in Section 7.


## 2. State-Space Model of Frequency Control in Multi-Microgrids Based on Edge Computing Framework

The power resources in multi-microgrids mainly contain synchronous generators, renewable generators (such as wind turbines and photovoltaic cells), and energy storage systems. Due to the uncertain fluctuations of renewable outputs and load demands, the frequency may deviate from the rated value. Hence, it is required to timely adjust the outputs of synchronous generators and energy storage systems to maintain real-time power balance between supply and demand at the rated frequency point [19, 20]. Considering that there exist the heavy communication and computation burdens in the traditional centralized control scheme which requires global operation statuses of all the microgrids, in this paper, a frequency control strategy based on edge computing framework is proposed. As shown in Figure 1, the local operating statuses are transmitted to the edge computing device installed in each microgrid. In other words, in our proposed edge computing framework, an edge computing device considered as the local control center is set up in each microgrid. Since the local controller is closer to the participating equipment and the corresponding control structure is simpler, lighter communication burden and lower computation cost can be realized. Similar to the centralized scheme, the function of the edge computing device is sending the active power output commands of synchronous generators and energy storage systems and then restoring the frequency to the rated value. Note that the edge computing devices and devices participating in frequency control are geographically dispersed; the uploading of operation statuses usually relays on the sharing communication network. Hence, there exist time delay and packet loss problems during the data transmission process. In this section, the state-space model of frequency control system in each microgrid is firstly discussed. The effects of time delay and packet loss on the dynamics of frequency control system will be analysed in the next sections.Figure 1
Framework of multi-microgrids with edge computing framework.Let the number of microgrids in a multi-microgrid beN. For the i-th (i=1,2,⋯,N) microgrid, the equivalent block diagram of the frequency control system is shown in Figure 2. Moreover, let the number of synchronous generators and energy storage systems in i-th microgrid be MSG_i and MESS_i, respectively. Therefore, the dynamics of the frequency control system in i-th microgrid satisfy [4]
(1)ddtΔPvij=−1TgijΔPvij+1Tgijαijui−1RijΔfi,ddtΔPmij=−1TtijΔPmij+1TtijΔPvij,ddtΔPtie_i=2π∑i=1,i≠lNTilΔfi−Δfl,ddtΔfi=−Deq_iMeq_iΔfi+1Meq_i∑j=1MSG_iΔPmij+ΔPres_i+∑k=1MESS_iΔPESS_ik−ΔPdi−ΔPtie_i,ΔPESS_ik=βikui−1RESS_ikΔfi,where ΔPvij is the valve opening of the j-th synchronous generator; ΔPmij is the mechanical output power of the j-th synchronous generator; ΔPtie_i is the tie line power fluctuation in i-th microgrid; Δfi and Δfl are the frequency deviations in the i-th and l-th microgrids, respectively; Til is the synchronization coefficient of the tie line between the i-th and l-th microgrids; Tgij and Ttij are the time constants of governor and turbine in the j-th synchronous generator; Rij and RESS_ik are the drooping coefficients of the j-th synchronous generator and k-th energy storage system, respectively; ΔPESS_ik is the active output power of the k-th energy storage system; ΔPdi and ΔPres_i are the fluctuations of load demand and renewable output in the i-th microgrids, respectively; Meq_i and Deq_i are the equivalent inertia and damping coefficient of the i-th microgrids, respectively; αij and βik are the participation factors of the j-th synchronous generator and k-th energy storage system, respectively, and satisfy
(2)∑j=1MSG_iαij+∑k=1MESS_iβik=1.Figure 2
Control block of networked frequency control system.By defining the state vector asXi=Δfi,ΔPtie_i,ΔPmi1,⋯,ΔPmiMSG_i,ΔPvi1,⋯,ΔPviMSG_iT, the state-space model of the frequency control system in the i-th microgrid is given by
(3)ddtXit=AiXit+Biuit+Hiwit,Yit=CiXit,where
(4)Ai=−Deq_iMeq_i−1Meq_i∑k=1MESS_i1RESS_ik−1Meq_i1Meq_i⋯1Meq_i0⋯02π∑i=1,i≠lNTilΔfi00⋯00⋯000−1Tti1⋯01Tti1⋯0⋮⋮⋮⋱⋮⋮⋱⋮000⋯−1TtiMSG_i0⋯1TtiMSG_i1Tgi1Ri100⋯0−1Tgi1⋯0⋮⋮⋮⋱⋮⋮⋱⋮1TgiMSG_iRiMSG_i00⋯00⋯−1TgiMSG_i,Bi=∑k=1MESS_iβikMeq_i00⋮0αi1Tgi1Ri1⋮αiMSG_iTgiMSG_i,Hi=−1Meq_i00−2π00⋮⋮0000⋮⋮00,Ci=100⋮00⋮0T,wi=ΔPdi−ΔPres_i∑i=1,i≠lNTilΔfl.Moreover, considering that the frequency control system is a sampling control system under the edge computing framework essentially, the respective discrete-time model with the sampling timeTs is given by
(5)Xik+1=EiXik+Fiuik+Giwik,Yik=CiXik,where Ei=eAiTs, Fi=∫0TseAitdtBi, and Gi=∫0TseAitdtHi.


## 3. Transmission Delay Analysis in Edge Computing Environment

Generally, the total time delay during the transmission process from the underlying participating equipment to the edge computing device includes the following three parts: serial delay, propagation delay, and routing delay [16]. Hence, the following equation holds.
(6)τΣ=τserial+τpropagation+τrouting,where τΣ is the total time delay. Besides, the serial delay (τserial) is proportional to the size of packets (denoted as Psize) and inversely proportional to the transmission rate (denoted as Dtransmission). The propagation delay is proportional to the transmission distance (denoted as LengthΣ) and inversely proportional to the propagation velocity (denoted as V) in a certain physical media. Obviously, both of the serial delay and propagation delay are constants and can be calculated by the following equations.
(7)τserial=PsizeDtransmission,τpropagation=LengthΣV.In particular, the uncertain routing delay is caused by the queuing and forwarding of multiple data flows in the switches. Specifically, as shown in Figure3, the operating statuses transmitted from the underlying equipment to the edge computing device need to pass through multiple switches with limited forward capacities. Therefore, at the entrance of a switch, the data flows from different receivers form a queue. Such queuing process can be described as a cascade M/M/1 model based on the queuing theory [21].Figure 3
Multihop transmission model in edge computing environment.Letμr denote the forwarding rate of the r-th (r=1,2,⋯,R) switch provided for the data flow and τr denote the queuing delay of the data flow transmitted by the r-th switch. Assume that a data flow passes through the switches numbered as 1,2,⋯,R in turn and the data arrival process at the r-th switch obeys the Poisson process with parameter λr. Based on the queuing theory, the probability density function (PDF) of the queuing delay at the r-th switch satisfies
(8)gτr=μr−λre−μr−λrτr,τr>0.The corresponding cumulative distribution function (CDF) is given by(9)g¯τr=∫0τrgtdt=1−e−μr−λrτr,τr>0.It can be seen that the queuing delay inr-th switch satisfies an exponential distribution with the parameter μr−λr. Furthermore, when the switches 1−R are independent to each other, the total routing delay (τrouting) satisfies
(10)τrouting=∑r=1Rτr.The corresponding PDF and CDF are given by Equations (11) and (12), respectively.
(11)gτrouting=∏r=1Rμr−λre−∑r=1Rμr−λrτr,τr>0,(12)g¯τrouting=∏r=1R1−e−μr−λrτr,τr>0.Finally, considering that the serial delay and propagation delay are constants, the PDF and CDF of the total transmission delay (τΣ) in the proposed edge computing environment equate to those of the total routing delay (τrouting).


## 4. Dynamic Characteristic Analysis of Closed-Loop Frequency Control System with Time Delay and Packet Loss

In this paper, to economize the memory space of edge computing devices, each microgrid adopts the memoryless state-feedback control mode, as shown in(13)uik=KiXinewestk,where Ki is the controller gain and Xinewestk is the newest packet arrived at the edge computing device.Noting that the total transmission delay (τΣ) always satisfies τΣ≤LTsL∈Z+, in this paper, LTs is defined as the preset maximum transmission delay. In other words, if the transmission delay of a packet is less than LTs, then this packet is called an effective packet and will be used to generate the control instructions uik. Otherwise, if the transmission delay of a packet exceeds LTs, this packet is viewed as a dropped one and will not be used to generate the control instruction uik. Obviously, the different transmission delay may cause that the edge computing device use different packets to generate the control instructions. Hence, the corresponding dynamic characteristics of closed-loop frequency control system will switch with the varying time delays.Without loss of generality, letpv and pv+1 be the serial number of the two neighbouring effective packets sampled at t=pvTs and t=pv+1Ts. Both the two packets are used to generate the control instructions due to the fact that their transmission delays are within 0,LTs. Obviously, the packets sampled within pv+1Ts,pvTs are dropped due to the transmission delays exceeding LTs. The maximum number of consecutive dropped packets can be denoted as Dmax=maxpv+1−pv. Then, the dynamics of closed-loop frequency control system during any two neighbouring effective packets (i.e., t∈pvTs,pv+1Ts) can be described as a switching model dependent on the time-varying delays. In order to clearly illustrate the modelling process for the closed-loop frequency control system with time delay and packet loss, in this section, we take L=1 as an example. As shown in Figure 4, there are two possible scenarios of the closed-loop frequency control system with varying time delays.Figure 4
Effects of time delays on the closed-loop frequency control system.Scenario 1:τΣ=0, i.e., the packets sampled at t=pvTs arrive at the edge computing device without any delay. Therefore, the newest packet satisfies Xinewestk=Xipv, t∈pvTs,pv+1Ts. According to the discrete-time model of frequency control system (i.e., Equation (5)), the dynamics of close-loop frequency control systems during the time interval pvTs,pv+1Ts can be derived as follows:
(14)Xipv+1=Ei+FiKiXipv+Giwipv,Xipv+2=Ei+FiKiXipv+1+Giwipv+1=Ei2+EiFiKi+FiKiXipv+EiGiwipv+Giwipv+1,⋮Xipv+1=EiDmax+1+∑ξ=0DmaxEiξFiKiXipv+∑ξ=0DmaxEiDmax−ξGiwipv+ξ.Scenario 2:τΣ≤1×Ts, i.e., the packet sampled at t=pvTs is transmitted to the edge computing device with a delay less than 1×Ts. In this scenario, the control instruction during the time interval pvTs,pv+1Ts satisfies the following piecewise function
(15)Xinewestk=Xipv−1,pvTs≤t<pv+1Ts,Xipv,pv+1Ts≤t<pv+1Ts.Similar to Scenario 1, the dynamics of close-loop frequency control systems during the time intervalpvTs,pv+1Ts can be derived as follows:
(16)Xpv+1=EiDmax+1+∑ξ=0Dmax−1EiξFiKiXipv+EiDmaxFiKiXipv−1+∑ξ=0DmaxEiDmax−ξGiwipv+ξ.More generally, the above derivation process in this section can be extended to case ofL∈N+. Define the following two augmented vectors:
(17)X¯ipv=XiTpv,XiTpv−1,⋯,XiTpv−LT,Wipv=wiTpv,wiTpv+1,⋯,wiTpv+LT.Hence, there exist2L possible scenarios of the closed-loop frequency control system in the i-th microgrid, uniformly described by
(18)X¯ipv+1=EiDmax+1+∑ξ=0Dmax−ζEiξFiKi∑ξ=Dmax−ζ+1ψ1EiξFiKi⋯∑ξ=Dmax−ζ+l−1ψL−1EiξFiKi∑ξ=Dmax−ζ+LψLEiξFiKiI0⋯000I⋯00⋮⋮⋮⋮⋮00⋯I0⏟ΦiζX¯ipv+EiDmaxGiEiDmax−1Gi⋯Gi00⋯0⋮⋮⋮⋮00⋯0⏟ΛiWipv,Ypv=Ci,0,0,⋯,0⏟L⏟C¯iX¯ipv,where ζ=0,1,⋯,L; the values of Ψ1,Ψ2,⋯,ΨL are shown in Equation (19) where the arrow “⟶” represents Ψ1,Ψ2,⋯,ΨL take values from the same row.
(19)


## 5.H∞-Controller Optimization considering Dynamic Performance Improvement

The controller gainKi not only needs to ensure that the closed-loop frequency control systems are asymptotically stable during any two neighbouring effective packets but also require to have a robust attenuation ability against the external power disturbances from the load demands and renewable outputs. In other words, there exists the following inequality between the output Yipv and the external power disturbances Wipv at any time t=pvTs.
(20)Yipv2≤γ2Wipv2,where γ is the attenuation factor and ·2 is the 2-norm operator. Let kpv and kpv+1∈1,2,⋯,2L denote the serial numbers of the possible operation scenarios at time t=pvTs and t=pv+1Ts. A candidate Lyapunov function is constructed by
(21)Vipv=X¯iTpvΩkpvX¯ipv,where Ωkpv is a symmetric positive definite matrix. When the state vector of closed-loop frequency control system changes from X¯ipv to X¯ipv+1, the increment of Equation (21) is given by
(22)ΔVipv=Vipv+1−Vipv=X¯iTpv+1Ωkpv+1X¯ipv+1−X¯iTpvΩkpvX¯ipv.Substituting Equation (18) and inequality (20) in Equation (22) results in
(23)ΔVipv≤X¯ipvWipvTΦiζTΩkpvΦiζ−Ωkpv+1ΦiζTΩkpvΛi∗ΛiTΩkpvΛiX¯ipvWipv−YiTpvYipv+γ2WiTpvWipv=X¯ipvWipvTΦiζTΩkpvΦiζ−Ωkpv+1+C¯iTC¯iΦiζTΩkpvΛi∗ΛiTΩkpvΛi−γ2IX¯ipvWipv,where I is an identity matrix with suitable dimensions. When the following matrix inequality
(24)ΦiζTΩkpvΦiζ−Ωkpv+1+C¯iTC¯iΦiζTΩkpvΛi∗ΛiTΩkpvΛi−γ2I<0holds, the integral of Lyapunov function on [0, +∞) with zero initial conditions satisfies
(25)0≤Vi+∞−Vi0≤∑pv=0+∞−YiTpvYipv+γ2WiTpvWipv=−Yipv2+γ2Wipv2.Overall, if the matrix inequality (24) holds, the control gain Ki can guarantee the closed-loop frequency control system has the asymptotic stability in the case of stochastic communication changes and H∞-robust attenuation performance against external power disturbances simultaneously.Besides, in order to improve the dynamic performance during frequency restoration process, the integral of absolute error (IAE) of the frequency is chosen as the objective function, as shown in Equation (26). The IAE of frequency reflects the transient response performance of the frequency control system to external power disturbances
(26)Ji=∑k=0∞Δfik.In conclusion, a constrained optimization model for the controller gainKi is established, as follows.
(27)objective:minJis.t.:ΦiζTΩkpvΦiζ−Ωkpv+1+C¯iTC¯iΦiζTΩkpvΛi∗ΛiTΩkpvΛi−γ2I<0Ωkpr=ΩkprT>0,Ωkpr+1=Ωkpr+1T>0.According to Equation (18), the matrix Φiζ contains the controller gain Ki which needs to be optimized. Besides, there exists a multiplicative relation between matrix Φiζ and unknown matrix Ωkpv. Therefore, the constraints (27) do not satisfy the linear matrix inequality (LMI) forms and cannot be solved by directly using the robust control toolbox of MATLAB directly. In this paper, an iterative relaxation technology combined with heuristic search method is proposed to optimize the controller gain Ki. The flowchart is shown in Figure 5. Firstly, an initial solution set containing considerable matrixes with the same dimensions of Ki is generated randomly. Then, the initial feasible solution set is selected from the initial solution set according to the constraints in (27). Secondly, the initial solution set is updated via the operations such as mutation and pheromone update. The constraints (27) are used again to select the feasible solution set for the next iteration. Since the candidates in initial solution set for Ki are fixed during the whole optimization process, the matrix nonlinear matrix inequality in (27) is relaxed into the LMIs.Figure 5
Flowchart for controller gain optimization.


## 6. Simulations and Discussions

In this section, the feasibility of the proposed cascade M/M/1 model for communication network in edge computing framework is verified by using OPNET Modeler 14.5. In addition, assume that the number of microgrids in the multi-microgrid is 2 and each microgrid contains two synchronous generators and one energy storage system. The frequency control in multi-microgrid is simulated in MATLAB/Simulink environment, as shown in Figure6. The time delays during the simulation process adopt the simulation results obtained from the OPNET. The controller gain is obtained according to the optimization algorithm proposed in Section 5 where the particle swarm optimization method is used. The synchronization coefficient between the two microgrids satisfies T12=T21=0.52p.u./Hz. Besides, the communication system configurations in both two microgrids are assumed to be the same. To evaluate the dynamic performance of the proposed control strategy, the step changes of renewable output and load demand with 0.1 p.u. and -0.1 p.u. occur at t=0s in each microgrid. Other simulation parameters are demonstrated in Table 1.Figure 6
Simulation diagram of multi-microgrids.Table 1
Simulation parameters.
ParametersValuesIllustrationsTg11 (s)0.31Microgrid 1Tt11 (s)0.05Tg12 (s)0.28Tt12 (s)0.07R11 (Hz/p.u.)3R12 (Hz/p.u.)2.8RESS_11 (Hz/p.u.)5Meq_1 (p.u. s)0.2308Deq_1 (p.u./Hz)0.016α11,α12,β11(0.4, 0.4, 0.2)Tg21 (s)0.35Microgrid 2Tt21 (s)0.06Tg22 (s)0.33Tt22 (s)0.08R21 (Hz/p.u.)2.87R22 (Hz/p.u.)2.5RESS_21 (Hz/p.u.)4.5Meq_2 (p.u. s)0.2408Deq_2 (p.u./Hz)0.018α21,α22,β21(0.3, 0.6, 0.1)Ts (ms)10Communication networksPsize (bits)200Dtransmission (Mbps)10LengthΣ (km)50V (km/s)1.8×105λr (packet/s)50μr (packet/s)500R5


### 6.1. Transmission Delays in Edge Computing Framework

As shown in Figure7, a multihop data transmission network is established in OPNET environment to verify the proposed cascade M/M/1 model for time delay calculation under the edge computing framework. It should be noted that since the number of transmission hops is equal to the number of forwarding nodes, a chain topology is adopted in the simulation. Figure 8 shows the statistical results of transmission delays of 100 consecutive packets and the theoretical mean value of time delays calculated by the proposed M/M/1 model.Figure 7
Schematic diagram of communication network in edge computing framework.Figure 8
Statistical results of transmission delays of 100 consecutive packets in edge computing environment.According to Figure8, the simulation mean value of the transmission delays is 10.15 ms. Meanwhile, the theoretical mean value calculated by the proposed cascade M/M/1 model is 11.41 ms. The relative error between theoretical and simulation values is 11.04%. Therefore, the simulation results illustrate that the proposed cascade M/M/1 model can effectively reflect the transmission characteristics in edge computing environment. In addition, it can be found that the number of continuous packets with transmission delays exceeding 20 ms is less than 3. Therefore, in controller design process, the preset maximum transmission delay is assumed to be 20 ms (i.e., 2×10ms) while the maximum number of consecutive dropped packets (i.e., Dmax) is assumed to be 3.


### 6.2. Control Performance Analysis with the Proposed Method

#### 6.2.1. Comparisons with the Traditional Centralized Control Scheme

In the traditional centralized control schemes, all the operation statuses of the microgrids should be transmitted to the only control center in the whole power system. Without loss of generality, the simplest case that there is only one switch between the edge computing device and the control center is considered in this subsection, as shown in Figure9. Figure 10 shows the statistical results of transmission delays of 100 consecutive packets and the theoretical mean value of time delays calculated by the proposed M/M/1 model. According to Figure 10, the number of continuous packets with transmission delays which exceeds 20 ms is less than 6. Therefore, the preset maximum transmission delay is assumed to be 20 ms (i.e., 2×10ms) while the maximum number of consecutive dropped packets (i.e., Dmax) is assumed to be 6. In addition, in order to quantitatively evaluate the control performance, the following indexes are selected:
(1)
Peak value of frequency deviation: the maximum value of the absolute value of frequency deviation when step changes of external power disturbances occur(2)
Restoration time of frequency deviation: the starting time when the frequency deviation restores and maintains within ±5% of the peak value(3)
IAE of the frequency: the result calculated according to Equation (26)Figure 9
Schematic diagram of communication network in traditional centralized framework.Figure 10
Statistical results of transmission delays of 100 consecutive packets in traditional centralized framework.In this subsection, the controller gains of both centralized scheme and edge computing framework are obtained with the proposedH∞-switching control method in Section 5. Figure 10 shows the frequency deviation responses with the traditional centralized and the proposed edge computing framework. Table 2 demonstrates the corresponding performance indexes. For the first microgrid, according to Figures 10 and 11 and Table 2, although the restoration times of the two control schemes are close, the peak value of our proposed edge computing framework is reduced by 32.51% compared with the traditional centralized control scheme. Meanwhile, for the second microgrid, not only the peak value of our proposed edge computing framework is reduced by 43.51% but also the corresponding restoration time can be shortened by 52.93% compared with the centralized control scheme. In addition, the IAEs in both two microgrids using the proposed edge computing framework are less than the results using the centralized scheme. The simulation results show that the proposed edge computing framework can provide a better frequency deviation damping performance. This is because the proposed control strategy based on edge computing framework only requires the local operation statuses and the edge computing device which is closer to the underlying participating equipment. Hence, compared with the traditional centralized control scheme, the communication network with edge computing framework can provide a higher quality of service (QoS) for the data flows and has the less consecutive dropped packets.Table 2
Performance indexes with traditional centralized framework and proposed edge computing framework.
MicrogridsControl strategiesPeak value (Hz)Restoration time (s)IAE (Hz)1This paper0.01641.360.4182Centralized framework0.02431.140.60922This paper0.02092.010.8920Centralized framework0.03704.273.5346Figure 11
Frequency deviations with traditional centralized framework and proposed edge computing framework.


#### 6.2.2. Comparisons with Different Controller Design Methods in Edge Computing Framework

In this subsection, the superiority of dynamic performance with our proposedH∞-switching control strategy compared with other design methods is discussed. The comparison methods are as follows: (1) H∞ control strategy in [9]; (2) MPC strategy in [17]; (3) sliding control strategy in [18]. All the control strategies in this subsection adopt the decentralized scheme and can be applied into the edge computing framework.Figure12 shows the frequency deviation responses with different controller design methods. Table 3 demonstrates corresponding dynamic performance indexes. For the first microgrid, according to Figure 12 and Table 3, compared with the sliding control strategy, the maximum frequency deviation with our proposed H∞-switching control is 32.3% larger while the restoration time is shortened by 58.28%. Compared with H∞ control, the maximum frequency deviation and restoration time are reduced by 33.87% and 17.58%, respectively. Compared with the MPC method, the maximum frequency deviation and restoration time are reduced by 1.80% and 72.53%, respectively. Note that the IAE of frequency with the proposed H∞-switching control is reduced by 37.19%, 42.40%, and 68.54%, respectively, compared with the sliding control, H∞ control, and MPC. Similar results can be obtained for the second microgrid. The reason is that our proposed controller design method optimizes the dynamic performance during frequency restoration process on the premise of ensuring the asymptotical stability of the closed-loop frequency control systems. Therefore, the shorter restoration time and better transition process can be obtained.Figure 12
Frequency deviations with different control strategies.Table 3
Performance indexes with different control strategies.
MicrogridsControl strategiesPeak value (Hz)Restoration time (s)IAE (Hz)1This paper0.01641.360.4182Sidling control0.01243.260.6658H∞ control0.02481.650.7259MPC0.01674.951.32942This paper0.02092.010.8920Sidling control0.02661.731.0061H∞ control0.03802.932.4598MPC0.02703.311.6974


## 6.1. Transmission Delays in Edge Computing Framework

As shown in Figure7, a multihop data transmission network is established in OPNET environment to verify the proposed cascade M/M/1 model for time delay calculation under the edge computing framework. It should be noted that since the number of transmission hops is equal to the number of forwarding nodes, a chain topology is adopted in the simulation. Figure 8 shows the statistical results of transmission delays of 100 consecutive packets and the theoretical mean value of time delays calculated by the proposed M/M/1 model.Figure 7
Schematic diagram of communication network in edge computing framework.Figure 8
Statistical results of transmission delays of 100 consecutive packets in edge computing environment.According to Figure8, the simulation mean value of the transmission delays is 10.15 ms. Meanwhile, the theoretical mean value calculated by the proposed cascade M/M/1 model is 11.41 ms. The relative error between theoretical and simulation values is 11.04%. Therefore, the simulation results illustrate that the proposed cascade M/M/1 model can effectively reflect the transmission characteristics in edge computing environment. In addition, it can be found that the number of continuous packets with transmission delays exceeding 20 ms is less than 3. Therefore, in controller design process, the preset maximum transmission delay is assumed to be 20 ms (i.e., 2×10ms) while the maximum number of consecutive dropped packets (i.e., Dmax) is assumed to be 3.


## 6.2. Control Performance Analysis with the Proposed Method

### 6.2.1. Comparisons with the Traditional Centralized Control Scheme

In the traditional centralized control schemes, all the operation statuses of the microgrids should be transmitted to the only control center in the whole power system. Without loss of generality, the simplest case that there is only one switch between the edge computing device and the control center is considered in this subsection, as shown in Figure9. Figure 10 shows the statistical results of transmission delays of 100 consecutive packets and the theoretical mean value of time delays calculated by the proposed M/M/1 model. According to Figure 10, the number of continuous packets with transmission delays which exceeds 20 ms is less than 6. Therefore, the preset maximum transmission delay is assumed to be 20 ms (i.e., 2×10ms) while the maximum number of consecutive dropped packets (i.e., Dmax) is assumed to be 6. In addition, in order to quantitatively evaluate the control performance, the following indexes are selected:
(1)
Peak value of frequency deviation: the maximum value of the absolute value of frequency deviation when step changes of external power disturbances occur(2)
Restoration time of frequency deviation: the starting time when the frequency deviation restores and maintains within ±5% of the peak value(3)
IAE of the frequency: the result calculated according to Equation (26)Figure 9
Schematic diagram of communication network in traditional centralized framework.Figure 10
Statistical results of transmission delays of 100 consecutive packets in traditional centralized framework.In this subsection, the controller gains of both centralized scheme and edge computing framework are obtained with the proposedH∞-switching control method in Section 5. Figure 10 shows the frequency deviation responses with the traditional centralized and the proposed edge computing framework. Table 2 demonstrates the corresponding performance indexes. For the first microgrid, according to Figures 10 and 11 and Table 2, although the restoration times of the two control schemes are close, the peak value of our proposed edge computing framework is reduced by 32.51% compared with the traditional centralized control scheme. Meanwhile, for the second microgrid, not only the peak value of our proposed edge computing framework is reduced by 43.51% but also the corresponding restoration time can be shortened by 52.93% compared with the centralized control scheme. In addition, the IAEs in both two microgrids using the proposed edge computing framework are less than the results using the centralized scheme. The simulation results show that the proposed edge computing framework can provide a better frequency deviation damping performance. This is because the proposed control strategy based on edge computing framework only requires the local operation statuses and the edge computing device which is closer to the underlying participating equipment. Hence, compared with the traditional centralized control scheme, the communication network with edge computing framework can provide a higher quality of service (QoS) for the data flows and has the less consecutive dropped packets.Table 2
Performance indexes with traditional centralized framework and proposed edge computing framework.
MicrogridsControl strategiesPeak value (Hz)Restoration time (s)IAE (Hz)1This paper0.01641.360.4182Centralized framework0.02431.140.60922This paper0.02092.010.8920Centralized framework0.03704.273.5346Figure 11
Frequency deviations with traditional centralized framework and proposed edge computing framework.


### 6.2.2. Comparisons with Different Controller Design Methods in Edge Computing Framework

In this subsection, the superiority of dynamic performance with our proposedH∞-switching control strategy compared with other design methods is discussed. The comparison methods are as follows: (1) H∞ control strategy in [9]; (2) MPC strategy in [17]; (3) sliding control strategy in [18]. All the control strategies in this subsection adopt the decentralized scheme and can be applied into the edge computing framework.Figure12 shows the frequency deviation responses with different controller design methods. Table 3 demonstrates corresponding dynamic performance indexes. For the first microgrid, according to Figure 12 and Table 3, compared with the sliding control strategy, the maximum frequency deviation with our proposed H∞-switching control is 32.3% larger while the restoration time is shortened by 58.28%. Compared with H∞ control, the maximum frequency deviation and restoration time are reduced by 33.87% and 17.58%, respectively. Compared with the MPC method, the maximum frequency deviation and restoration time are reduced by 1.80% and 72.53%, respectively. Note that the IAE of frequency with the proposed H∞-switching control is reduced by 37.19%, 42.40%, and 68.54%, respectively, compared with the sliding control, H∞ control, and MPC. Similar results can be obtained for the second microgrid. The reason is that our proposed controller design method optimizes the dynamic performance during frequency restoration process on the premise of ensuring the asymptotical stability of the closed-loop frequency control systems. Therefore, the shorter restoration time and better transition process can be obtained.Figure 12
Frequency deviations with different control strategies.Table 3
Performance indexes with different control strategies.
MicrogridsControl strategiesPeak value (Hz)Restoration time (s)IAE (Hz)1This paper0.01641.360.4182Sidling control0.01243.260.6658H∞ control0.02481.650.7259MPC0.01674.951.32942This paper0.02092.010.8920Sidling control0.02661.731.0061H∞ control0.03802.932.4598MPC0.02703.311.6974


## 6.2.1. Comparisons with the Traditional Centralized Control Scheme

In the traditional centralized control schemes, all the operation statuses of the microgrids should be transmitted to the only control center in the whole power system. Without loss of generality, the simplest case that there is only one switch between the edge computing device and the control center is considered in this subsection, as shown in Figure9. Figure 10 shows the statistical results of transmission delays of 100 consecutive packets and the theoretical mean value of time delays calculated by the proposed M/M/1 model. According to Figure 10, the number of continuous packets with transmission delays which exceeds 20 ms is less than 6. Therefore, the preset maximum transmission delay is assumed to be 20 ms (i.e., 2×10ms) while the maximum number of consecutive dropped packets (i.e., Dmax) is assumed to be 6. In addition, in order to quantitatively evaluate the control performance, the following indexes are selected:
(1)
Peak value of frequency deviation: the maximum value of the absolute value of frequency deviation when step changes of external power disturbances occur(2)
Restoration time of frequency deviation: the starting time when the frequency deviation restores and maintains within ±5% of the peak value(3)
IAE of the frequency: the result calculated according to Equation (26)Figure 9
Schematic diagram of communication network in traditional centralized framework.Figure 10
Statistical results of transmission delays of 100 consecutive packets in traditional centralized framework.In this subsection, the controller gains of both centralized scheme and edge computing framework are obtained with the proposedH∞-switching control method in Section 5. Figure 10 shows the frequency deviation responses with the traditional centralized and the proposed edge computing framework. Table 2 demonstrates the corresponding performance indexes. For the first microgrid, according to Figures 10 and 11 and Table 2, although the restoration times of the two control schemes are close, the peak value of our proposed edge computing framework is reduced by 32.51% compared with the traditional centralized control scheme. Meanwhile, for the second microgrid, not only the peak value of our proposed edge computing framework is reduced by 43.51% but also the corresponding restoration time can be shortened by 52.93% compared with the centralized control scheme. In addition, the IAEs in both two microgrids using the proposed edge computing framework are less than the results using the centralized scheme. The simulation results show that the proposed edge computing framework can provide a better frequency deviation damping performance. This is because the proposed control strategy based on edge computing framework only requires the local operation statuses and the edge computing device which is closer to the underlying participating equipment. Hence, compared with the traditional centralized control scheme, the communication network with edge computing framework can provide a higher quality of service (QoS) for the data flows and has the less consecutive dropped packets.Table 2
Performance indexes with traditional centralized framework and proposed edge computing framework.
MicrogridsControl strategiesPeak value (Hz)Restoration time (s)IAE (Hz)1This paper0.01641.360.4182Centralized framework0.02431.140.60922This paper0.02092.010.8920Centralized framework0.03704.273.5346Figure 11
Frequency deviations with traditional centralized framework and proposed edge computing framework.


## 6.2.2. Comparisons with Different Controller Design Methods in Edge Computing Framework

In this subsection, the superiority of dynamic performance with our proposedH∞-switching control strategy compared with other design methods is discussed. The comparison methods are as follows: (1) H∞ control strategy in [9]; (2) MPC strategy in [17]; (3) sliding control strategy in [18]. All the control strategies in this subsection adopt the decentralized scheme and can be applied into the edge computing framework.Figure12 shows the frequency deviation responses with different controller design methods. Table 3 demonstrates corresponding dynamic performance indexes. For the first microgrid, according to Figure 12 and Table 3, compared with the sliding control strategy, the maximum frequency deviation with our proposed H∞-switching control is 32.3% larger while the restoration time is shortened by 58.28%. Compared with H∞ control, the maximum frequency deviation and restoration time are reduced by 33.87% and 17.58%, respectively. Compared with the MPC method, the maximum frequency deviation and restoration time are reduced by 1.80% and 72.53%, respectively. Note that the IAE of frequency with the proposed H∞-switching control is reduced by 37.19%, 42.40%, and 68.54%, respectively, compared with the sliding control, H∞ control, and MPC. Similar results can be obtained for the second microgrid. The reason is that our proposed controller design method optimizes the dynamic performance during frequency restoration process on the premise of ensuring the asymptotical stability of the closed-loop frequency control systems. Therefore, the shorter restoration time and better transition process can be obtained.Figure 12
Frequency deviations with different control strategies.Table 3
Performance indexes with different control strategies.
MicrogridsControl strategiesPeak value (Hz)Restoration time (s)IAE (Hz)1This paper0.01641.360.4182Sidling control0.01243.260.6658H∞ control0.02481.650.7259MPC0.01674.951.32942This paper0.02092.010.8920Sidling control0.02661.731.0061H∞ control0.03802.932.4598MPC0.02703.311.6974


## 7. Conclusions

In this paper, an edge computing-based control scheme is proposed to deal with the frequency stability problem in multi-microgrids. Different from the traditional centralized scheme, the edge computing device is set up in each microgrid to realize the local frequency stability with higher QoS for data transmission and the lower computation burden for controller design. Via describing the closed-loop frequency control system as a delay-dependent switching model, the constraints of controller gain guaranteeing asymptotic stability andH∞-attenuation performance are strictly derived. Moreover, an optimization algorithm is proposed aiming at improving the transition performance during frequency restoration process. Simulation results show that our proposed method has shorter restoration time in stochastic time delay and packet loss cases. Further researches will focus on the frequency stability problem in edge computing framework suffering from malicious network attacks.



---

*Source: 6670591-2021-01-06.xml*