State feedback control aware stochastic transmitting latency in cyber-physical power system

Abstract

Distributed sensor networks with a great number of sensors realize information gathering, transmitting, and controlling, which will greatly improve the reliability and efficiency of industrial infrastructure systems, such as cyber-physical power system. Cyber-physical power system is one of the heuristic systems tightly coupled by a continuous-time electrical power system and a discrete-time information system. The transmitting latency, packets disorder, or loss will fail system convergence, even unstable under disturbance. This article studies the effect of long latency on cyber-physical power system frequency stability based on network control theory and establishes a heuristic model to express the continuous-time distributed generation system and the discrete-time distributed sensor network. A new state feedback controller is designed based on the stochastic optimal control theory to solve the micro-grid’s frequency stability problem with stochastic latency disturbance. Finally, we have conducted an extensive evaluation study using a real distributed generation micro-grid system. The simulation results show that the new controller reduces the influence of distributed sensor networks transmitting latency on the frequency stability, and the frequency dynamic process in the distributed generation micro-grid has smaller overshoot and obtains faster dynamic response speed.

Keywords

Distributed sensor networks cyber-physical power system network control theory transmitting latency distributed generation micro-grid system

Introduction

With the construction of smart grid, the traditional power system has gradually developed into the cyber-physical power system (CPPS) which is one of the heuristic systems tightly coupled by electrical power system and information system.^1,2 CPPS integrates ubiquitous sensors, communication networks, and high-performance information processing system.³ It combines real-time monitoring, data transmitting, dynamic controlling, and information service. Distributed sensor networks (DSNs), with their various sensors, self-organized communication network, and data-centric service, become one of the key techniques in CPPS.⁴ As we know, stochastic transmitting latency is one of the shortages in DSN, compared with wired communication network.⁵ And, the tens or even hundreds of milliseconds latency may cause CPPS control failure, which will make the electrical power system, the other side of tightly coupling, lose stability and fall into abnormal operation state. It results in significant economic losses and adverse social impact. Compared with traditional electrical power system with less monitoring or local measurement data and offline calculation, aware stochastic transmitting latency in CPPS control is a new problem.

In recent years, with the development of network control theory, scholars begin to research on the relationship between the transmission network and industrial control system. Network control system (NCS) is a type of closed-loop space-distributed control system. Each unit, such as sensors, controllers, and actuators, is distributed in different geographical locations, and they are connected by communication networks. Control methods with simplified network models with constant latency or packet loss rate are studied.^6,7 But, a CPPS controlling distributed generation and variable power load is one of the more complex networks and measuring data transmitting with stochastic transmitting latency. Therefore, there are two main difficulties; we proposed (1) to establish a frequency-controlled uncertain state space model of micro-grid, which is more consistent with the actual system considering the observation noise and (2) a controller to solve the micro-grid’s frequency stability problem with stochastic latency disturbance. This article deeply studies the effect of stochastic latency on CPPS frequency stability based on network control theory and established a heuristic model to express the continuous-time distributed generation system and the discrete-time communication system of micro-grid. A new the stochastic optimal state feedback controller is designed to stabilize power system frequency.

The remainder of the article is organized as follows: section “Background and related work” briefly reviews related work in the field. In section “Heuristic model of CPPS,” the heuristic model of CPPS is established, and the stochastic optimal state feedback controller in CPPS is proposed in section “Stochastic optimal state feedback controller designed in CPPS.” Section “Experiments, results, and algorithm performance evaluation analysis” presents the experimental results under different types of transmitting latency scenarios, and the frequency dynamic response control with different controllers are also evaluated. We conclude the article in section “Conclusion.”

Background and related work

In the study of NCS, Chen and Liu⁶ established an ideal communication model and then used the Lyapunov stability theory to reduce the system’s output error. Dong et al.⁸ used H_∞ norm constraints and linear matrix inequalities to obtain feedback control gain and system stability criterion. In Pang et al.,⁹ the output tracking error was converted to an additional state, and then, solving the output tracking problem is same as solving an augmented system stabilization problem. Based on the observer of the original state increment and the feedback of the output tracking error, a model-based networked predictive output tracking control scheme was proposed to compensate transmission latency. Shi and Yu¹⁰ formulated a closed-loop system as a special jump linear system and designed an output feedback controller based on mixed H₂ and H_∞ robust control theory. In Bai et al.,⁷ based on Lyapunov theory and linear matrix inequality methods, the sufficient and necessary condition for robust stability with transmission delay less than one sampling period and fixed delay is presented. The sufficient conditions for the stability of discrete switched systems with variable delays were studied in Nayeripour et al.,¹¹ and an algorithm with calculating a set of linear matrix inequalities to obtain state feedback controller was presented.

Besides the theoretical research, some scholars also tried to using NCS in industrial and energy fields. Bevrani and Hiyama¹² simplified the load frequency control problem as a static output feedback control problem and used proportional–integral (PI) control to design a robust decentralized PI controller. The latency in automatic generation control system is studied in Duan and He,¹³ and a time-delay system control algorithm was proposed.

From all the above analyses, the previous control methods always simplified the stochastic transmitting latency from communication network as constant, such as a fixed time slot less than one sampling period. But the topology of CPPS is more complex and there is serious electromagnetic noise, which reduce the DSNs’ transmission performance and result in stochastic transmitting latency. In this article, the micro-grid system with distributed generation system and storage is taken as the object to study the frequency stability when the long time-delay occurs in the control network. The continuous-time model of electrical power system and discrete-time model of DSNs with stochastic transmitting latency are established, and a new state feedback controller is proposed based on the stochastic optimal control theory.

Heuristic model of CPPS

Micro-grid model with distributed power generations

In many remote areas, micro-grid with diesel power generation is often used as the main standby power supply mode. With the increasing price of petroleum and the maturity of renewable energy technology, such as photovoltaic (PV) power and fuel cell power generation, cleanly distributed generations are more and more used in micro-grid as one type of power complementary.^14,15 Monitoring the operating status of distributed generations and controlling each part of power generation units are necessary to keep the stability of micro-grid system.¹⁶ A CPPS structure of micro-grid is shown in Figure 1.

Figure 1.

CPPS structure of micro-grid.

In Figure 1, dashed lines represent the flow of information. The distributed generation system consists of calculation control devices, data sensors, and electric devices. Information is exchanged by communication network, in which there are two parts of transmitting latency. One is, τ_ca, the process from data sensors to the control center, and the other is, τ_ac, from the control center to each actuator, such as distributed power generating units to adjust the operating status.

The latency τ = τ_ca + τ_ac makes not only the control center lose the ability to monitor the system but also the actuators unable to execute the correct issue instruction in time, which destroys the stability and leads to the frequency instability of the system.^17,18 A typical micro-grid system as an experiment example is employed to simulate the influence of communication transmission delay on system stability. There are one 150-kW diesel generator, one 120-kW micro gas generator, one 50-kW wind power generator, one 30-kW PV system, and one fuel cell energy storage system with the rated power as 50 kW. At 1.5 s, the power load reduces 10% suddenly. Figure 2 shows the frequency fluctuation curve with different transmission latencies, in which Δf = f − f_ref is the frequency deviation. f_ref is the reference frequency.

Figure 2.

Influence of communication delay on grid frequency.

From the experimental results, if there is no latency in the network, then the frequency of micro-grid can quickly smoothen during a disturbance processing; if there is 0.2 s transmitting latency, then the frequency of micro-grid can also be stabilized finally, but the regulation time increases significantly, from 3 s to 10 s, shown as a black curve; if the latency is 0.6 s and 1 s, then the micro-grid system’s frequency will not remain stable and the self-excited oscillation occurs. Based on the experiment results, the transmitting latency from communication network not only reduces the system convergence speed but also has the possibility to destabilize the system.

Continuous-time system mathematic model for micro-grid

A control mathematic model of the micro-grid with distributed generation system is presented in Figure 3. The distributed generation system consists of m diesel power generation systems, n PV power generation units, and one fuel cell generators. Each power generation unit receives frequency modulation task with participation factor α. For diesel power generation system, this article considers the speed governor and the turbine, which are described by first-order transfer function, and describes the PV and fuel cell power generation with simplified first-order transfer function.¹¹

Figure 3.

Model of distributed generation system.

Here, α_1p (p = 1, 2, …, m), α_2q (q = 1, 2, …, n), and α_F are the participation factors of distributed power generations and they meet $\sum_{p = 1}^{m} α_{1 p} + \sum_{q = 1}^{n}$ $α_{2 q} + α_{F} = 1$ ; G_G(s) and G_T(s) are the governor function and the turbine transfer function of diesel power generation system, respectively; G_PV(s) is the transfer function of PV power generation; G_FC(s) is the transfer function of fuel cell power generation; and ΔP_v and ΔP_m are the signal of prime mover increment and the prime mover output power increment, respectively. ΔP_PV is the output power increment of the PV power generation; ΔP_FC is the fuel cell power increment; G_{p_eq}(s) is the generator-load equivalent transfer function of micro-grid, and it can be calculated based on equipment equation (1); β is the frequency deviation factor of micro-grid; and ΔP_d is the load variation of micro-grid

G_{p_eq} (s) = \frac{1}{D_{eq} + s M_{eq}}

(1)

where D_eq and M_eq are the micro-grid system’s equivalent damping and inertia, respectively, which can be calculated by formula (2)

D_{e q} = \sum_{p = 1}^{m} D_{1 p} + \sum_{q = 1}^{n} D_{2 q} + D_{F C}, M_{e q} = \sum_{p = 1}^{m} M_{1 p} + \sum_{q = 1}^{n} M_{2 q} + M_{F C}

(2)

D _1p and M_1p are damping coefficients and inertia of diesel generator, respectively; D_2q, M_2q, D_FC, and M_FC are equivalent damping coefficients and equivalent inertia of PV power generation system and fuel cell power generation system, respectively. Each diesel power generation system meets the following relationships

{\begin{matrix} \frac{d Δ P_{mp}}{dt} = - \frac{Δ P_{mp}}{T_{Tp}} + \frac{Δ P_{vp}}{T_{Tp}} \\ \frac{d Δ P_{vp}}{dt} = - \frac{Δ P_{vp}}{T_{Gp}} + \frac{α_{1 p}}{T_{Gp}} u - \frac{1}{T_{Gp} R_{1 p}} Δ f \end{matrix}

(3)

Each PV unit meets the following relationships¹⁹

\frac{d Δ P_{PVq}}{dt} = - \frac{Δ P_{PVq}}{T_{PVq}} + \frac{α_{2 q}}{T_{PVq}} u - \frac{1}{T_{PV q} R_{2 q}} Δ f

(4)

The fuel cell power generation system meets the following relationships²⁰

\frac{d Δ P_{FC}}{dt} = - \frac{Δ P_{FC}}{T_{FC}} + \frac{α_{F}}{T_{FC}} u - \frac{1}{T_{FC} R_{FC}} Δ f

(5)

The following formula should be satisfied among the increment of power generation, the load variation of micro-grid, and the frequency shift

\sum_{p = 1}^{m} Δ P_{1 p} + \sum_{q = 1}^{n} Δ P_{2 q} \sum_{p = 1}^{m} Δ P_{m p} + Δ P_{F C} - Δ P_{d} = D_{e q} Δ f + M_{e q} \frac{d Δ f}{d t}

(6)

Besides, the frequency deviation factor β can be calculated by formula (7)

β = \sum_{p = 1}^{m} \frac{1}{R_{1 p}} + \sum_{q = 1}^{n} \frac{1}{R_{2 p}} + D_{eq}

(7)

In this article, the state variable is x (t) = [Δf, ΔP_m₁, ΔP_m₂, …, ΔP_mm, ΔP_v₁, ΔP_v₂, …, ΔP_vm, ΔP_PV₁, ΔP_PV₂, …, ΔP_PVn, ΔP_FC]^T. The continuous-time model of the system can be obtained from formulas (2)–(7), and the final formula is

{\begin{matrix} \frac{d x (t)}{dt} = Ax (t) + Bu (t) + H_{1} Δ P_{d} (t) \\ y (t) = Cx (t) \end{matrix}

(8)

where

\begin{matrix} A = [\begin{matrix} - \frac{D_{eq}}{M_{eq}} & A_{12} & O_{1 \times m} & A_{14} & \frac{1}{M_{eq}} \\ O_{m \times 1} & A_{22} & A_{23} & O_{m \times m} & 0 \\ A_{31} & O_{m \times m} & A_{33} & O_{m \times m} & 0 \\ A_{41} & O_{m \times m} & O_{m \times m} & A_{44} & 0 \\ - \frac{1}{T_{FC} R_{3}} & O_{m \times m} & O_{m \times m} & O_{m \times m} & - \frac{1}{T_{FC}} \end{matrix}] \\ A_{12} = A_{14} = [\begin{matrix} \frac{1}{M_{eq}} & \frac{1}{M_{eq}} & \dots & \frac{1}{M_{eq}} \end{matrix}] \\ A_{22} = - A_{23} = diag (\begin{matrix} - \frac{1}{T_{T 1}}, & - \frac{1}{T_{T 2}}, & \dots & - \frac{1}{T_{Tm}} \end{matrix}) \\ A_{31} = {[\begin{matrix} - \frac{1}{T_{G 1} R_{11}} & - \frac{1}{T_{G 2} R_{12}} & \dots & - \frac{1}{T_{Gm} R_{1 m}} \end{matrix}]}^{T} \\ A_{33} = diag (\begin{matrix} - \frac{1}{T_{G 1}}, & - \frac{1}{T_{G 2}}, & \dots & - \frac{1}{T_{Gm}} \end{matrix}) \\ A_{41} = {[\begin{matrix} - \frac{1}{T_{PV 1} R_{21}} & - \frac{1}{T_{PV 2} R_{22}} & \dots & - \frac{1}{T_{PVn} R_{2 n}} \end{matrix}]}^{T} \\ A_{44} = diag (\begin{matrix} - \frac{1}{T_{PV 1}}, & - \frac{1}{T_{PV 2}}, & \dots & - \frac{1}{T_{PVn}} \end{matrix}) \\ B = {[\begin{matrix} 0 & O_{m \times 1} & B_{3} & B_{4} & \frac{α_{F}}{T_{FC}} \end{matrix}]}^{T} \\ B_{3} = {[\begin{matrix} \frac{α_{11}}{T_{G 1}} & \frac{α_{12}}{T_{G 2}} & \dots & \frac{α_{1 m}}{T_{Gm}} \end{matrix}]}^{T} \\ B_{4} = {[\begin{matrix} \frac{α_{21}}{T_{PV 1}} & \frac{α_{22}}{T_{PV 2}} & \dots & \frac{α_{2 n}}{T_{PVn}} \end{matrix}]}^{T} \\ H_{1} = {[\begin{matrix} - \frac{1}{M_{eq}} & O_{m \times 1} & O_{m \times 1} & O_{n \times 1} & 0 \end{matrix}]}^{T} \\ C = [\begin{matrix} β & O_{1 \times m} & O_{1 \times m} & O_{1 \times n} & 0 \end{matrix}] \end{matrix}

In the model, the power load fluctuation of micro-grid ΔP_d(t) is the power loads’ fluctuation, which is a random distribution with zero-means.²¹ To be consistent with the actual system operating status, observation noise ε (t) is taken into account, which is also an independent normal random distribution with zero-means. The frequency-controlled uncertain state space model of distributed generation system under continuous-time can be obtained as formula (9)

{\begin{matrix} \frac{d x (t)}{dt} = Ax (t) + Bu (t) + ξ (t) \\ y (t) = Cx (t) + ε (t) \end{matrix}

where ξ (t) = H ₁ΔP_d is the load fluctuation of micro-grid and obeys normal distribution; ξ (t) also accords with normal distribution.

Discrete-time system mathematic model of DSNs

To collaborative control continuous-time power system and discrete-time communication network, this article discretized the continuous-time model, formula (9), into a discrete-time model with the sampling period T, in which the upper bound of transmitting latency is τ. The electrical values sampling period T is defined with “IEC 61850-Part 5: Communication Requirements for Functional and Device Models” or “Q/GDW 131-2006: Real-time Dynamic Monitoring System Technical Standard of Electric Power System.” In this article, T is selected 10 ms, the fastest one among three grades of 40, 20, or 10 ms defined in the Q/GDW 131-2006.

In order to study the problem conveniently, this article makes the following assumptions: (1) sensors are all time driven, and the controllers and actuators are event driven. (2) Latency in network control is bound and subject to a certain probability distribution. (3) Random sequences { ξ _k}, { ε _k}, and ${t_{i}^{k}}$ are all independent and identically distributed sequences; furthermore, they are also independent from each other. Obviously, the above three assumptions are reasonable in the actual DSNs.

Considering the bounded latency τ, two issues can be easily represented: one is τ > μT, where μ is a positive integer. It can be treated as packet lost with enough long transmission latency, which is one type of system failures because of control instructions lost. Fault tolerant control approach can solve this problem, which is not in this article’s research scope. The other one is 0 ≤ τ < μT, which means the control instruction arrives with a stochastic latency. This article uses the stochastic optimal control theory to solve this issue.

Because of the event driven, the control signal u(t) arrived at controlled object is piecewise continuous within a sampling period [kT, (k + 1)T), and it has at most h + 1 number of different values. It is assumed that the piecewise continuous values of u(t) occur at time points $kT + t_{i}^{k}$ (i = 0, 1, 2, …, h) during one sampling period [kT, (k + 1)T). Because the controller is also the type of event driven, u_k can be used to represent the output when the sensing signal X_k arrives at the controller, as shown in Figure 4. Formula (9) is integrated in the sampling period [kT, (k + 1)T), and the discrete-time state equations can be obtained as formula (10)

{\begin{matrix} x_{k + 1} = A_{s} x_{k} + \sum_{i = 0}^{h} B_{i}^{k} u_{k - i} + ξ_{k} \\ y_{k} = C x_{k} + ε_{k} \end{matrix}

(10)

where x _k= x (kT), y _k=y(kT), u _k= u (kT), A _s=e^AT, $B_{i}^{k} = \int_{t_{i}^{k}}^{t_{i - 1}^{k}} e^{A (T - s)} ds B$ , $ξ_{k} = \int_{kT}^{(k + 1) T} e^{A [(k + 1) T - s]} ξ (s) ds$ , and ε _k = ε(kT).

Figure 4.

Time sequence diagram of signal transmission in NCS.

Stochastic optimal state feedback controller designed in CPPS

In CPPS with stochastic transmitting latency, the obvious characteristics are as follows: (1) because of the stochastic long transmitting latency, the intermittency and fluctuation of renewable energy, and the randomness of the power load in micro-grid, it is hard to obtain the state feedback controller with zero-pole placement algorithm; (2) it is difficult to directly measure all of the state variables in an actual system because of system noise. So, it is necessary to design the observer to obtain full state variables. The Kalman observer and the state feedback controller will be designed, respectively, in this article.

Design of Kalman observer

Kalman observer can obtain the optimal state estimation signal $\hat{X} (t)$ with the minimum covariance matrix $E = {[x (t) - \hat{X} (t)] [x (t) - \hat{X} (t)]^{T}}$ . And then, this set of estimated signals can be used to replace the actual full state variables. The gain matrix K_f of the Kalman observer can be obtained by formula (11), and the structure of the Kalman observer is shown in Figure 5

K_{f} = P C^{T} Q_{2}^{- 1}

(11)

where P satisfies the following the Riccati algebraic equation

P A^{T} + AP - P C^{T} Q_{2}^{- 1} CP + Q_{1} = 0

(12)

Figure 5.

Structure of Kalman observer.

Design of stochastic optimal state feedback controller

From the discrete stochastic optimal control theory, the control performance index in the case of the selected state variable x can be calculated by formula (13)

J = E {x_{N}^{T} Q_{0} x_{N} + \sum_{k = 0}^{N - 1} [{x^{T}}_{k} Q_{1} x_{k} + {u^{T}}_{k} Q_{2} u_{k}]}

(13)

where Q ₀ is a symmetric nonnegative definite matrix, and Q ₁ and Q ₂ are symmetric positive definite matrices.

Since the discrete-time state equation in formula (10) contains { u _-1, u _-1, …, u _-h}, the augmented vector $z_{k} = [\begin{matrix} {\hat{X}}_{k}^{T} & u_{k - 1}^{T} & u_{k - 2}^{T} & \dots & u_{k - h}^{T} \end{matrix}]$ is used, and the original formula (10) can be represented as formula (14)

{\begin{matrix} z_{k + 1} = Φ_{k} z_{k} + Γ_{k} u_{k} + H ξ_{k} \\ y_{k} = C_{0} z_{k} + ε_{k} \end{matrix}

(14)

where

Φ_{k} = [\begin{matrix} A_{s} & B_{1}^{k} & B_{2}^{k} & … & B_{h - 1}^{k} & B_{h}^{k} \\ 0 & 0 & 0 & … & 0 & 0 \\ 0 & I & 0 & … & 0 & 0 \\ … & … & … & … & … \\ 0 & 0 & 0 & … & I & 0 \end{matrix}] \begin{matrix} Γ_{k} = {[\begin{matrix} B_{0}^{k} & Ι & 0 & … & 0 \end{matrix}]}^{T} \\ H = {[\begin{matrix} I & 0 & 0 & … & 0 \end{matrix}]}^{T} \\ C_{0} = {[\begin{matrix} C & 0 & 0 & … & 0 \end{matrix}]}^{T} \end{matrix}

The control system performance index J can be calculated as formula (15)

J = E {z_{N}^{T} {Q^{'}}_{0} z_{N} + \sum_{k = 0}^{N - 1} [z_{k}^{T} {Q^{'}}_{1} z_{k} + u_{k}^{T} {Q^{'}}_{2} u_{k}]}

(15)

where

{Q^{'}}_{0} = [\begin{matrix} Q_{0} & 0 & \dots & 0 \\ 0 & (h / h + 1) Q_{2} & \dots & 0 \\ \dots & \dots & \dots \\ 0 & 0 & \dots & (1 / h + 1) Q_{2} \end{matrix}], {Q^{'}}_{1} = [\begin{matrix} Q_{1} & 0 & \dots & 0 \\ 0 & (h / h + 1) Q_{2} & \dots & 0 \\ \dots & \dots & \dots \\ 0 & 0 & \dots & (1 / h + 1) Q_{2} \end{matrix}], and {Q^{'}}_{2} = \frac{1}{h + 1} Q_{2}

The control strategy is u _k=− L _k z _k, where $L_{k} = {(Q_{2}^{'} + Γ_{k}^{T} S_{k + 1} Γ_{k})}^{- 1} Γ_{k}^{T} S_{k + 1} Φ_{k}$ and $S_{k} = {(Φ_{k} - Γ_{k} L_{k})}^{T}$ $S_{k + 1} Φ_{k} Q_{1}, S (N) = Q'_{0}$ . The control system block diagram based on the stochastic optimal control theory is shown in Figure 6, in which r(t) is the reference value of frequency deviation.

Figure 6.

Control system block diagram based on the stochastic optimal control theory.

Experiments, results, and algorithm performance evaluation analysis

Experiment settings and parameters

In the experiments, the micro-grid with distributed generation system consists of two identical diesel generators, two identical PV power generation systems, and a fuel cell power generation system. The two diesel power generation systems are homogeneous, in which the time constant of governor is T_G = 0.1 s, time constant of steam turbine is T_T = 0.3 s, damping coefficient of generator is M = 0.15, inertia is D = 0.1, and the frequency modulation participation factor is α = 0.4. Two PV power generation systems are also homogeneous, in which the equivalent time constant is T_pv = 1.8 s, equivalent damping coefficients is M = 0.1, equivalent inertia is D = 0.15, and the frequency modulation participation factor α = 0.075. For the third type of power generation, fuel cell, the time constant T_FC = 0.25 s, equivalent damping coefficients M = 0.14, equivalent inertia D = 0.05, frequency modulation participation factor α = 0.05. The rated frequency of micro-grid is 50 Hz, the three weight matrices Q ₀ = Q ₁ = diag[0.1,10,10,1,1,10,10,10] and Q ₂ = 0.01, and the variance of measurement noise R₁ = 10^-8, R = 2.4 Hz/pukW. In the DSNs to exchange information among sensors, actuator, and controllers, IEEE 802.11 protocol is used, and the carrier-sense multiple access with collision avoidance (CSMA/CA) wireless channel visiting mechanism detects the link connectivity. Between each two communication terminals, the free space propagation channel model is employed with the 300 kbps bandwidth. Two types of packet lengths are set: one is 1500 bytes for data packets and the other is 128 bytes for control instructions.

Different sampling period for discrete-time systems

During the process of discretizing one continuous-time model into a discrete-time model, the sampling period T will directly affect the accuracy of dynamic process description. In Figure 7, the black curve is the continuous-time dynamic response. We set three different sampling periods T = 0.1 s, T = 0.01 s, and T = 0.001 s to evaluate the discrete-time systems’ frequency dynamic response tracking performance.

Figure 7.

Influence of sampling period on system.

The experiment results show that, when the sampling period is 0.1 s, as a long sampling period, sensors can’t acquire and transmit the system status in time. The controller cannot track the dynamic response completely, shown as the green curve. When the sampling period is too short, as T = 0.001 s, the maximum frequency deviation of the discrete system is too large to tolerate, Δf = 0.6 Hz. Only the sampling period is selected as 0.01 s, the dynamic response red curve of the discrete-time system coincides with the target curve, which can accurately describe the dynamic process of the system. The experiment results show that the overshoot of frequency deviation is reduce from 0.6 Hz (T = 0.001 s) to 0.14 Hz. In the following experiments, the sampling period is set 0.01 s to obtain the good control performance.

Evaluated controller performance in the constant transmitting latency

In the first experiment, we evaluated the proposed state feedback controller performance with the constant transmitting latency, compared with H₂ controller proposed in Matni.²² The frequency dynamic response curves are shown in Figure 8 under different transmitting latencies: 0.2, 0.6, and 1 s.

Figure 8.

Frequency dynamic response with constant transmitting latency: (a) frequency response with τ = 0.2 s, (b) frequency response with τ = 0.6 s, and (c) frequency response with τ = 1 s.

Figure 8 shows that, in different transmitting latencies, the proposed state feedback controller has good frequency dynamic response performance. The maximum frequency deviation kept to 0.06 Hz, compared with the 0.12 Hz by H₂ controller. The response speed is 5 s with the 0.2 s latency scenarios, 7 s with the 0.6 s latency scenarios, and 8 s with the 1 s latency scenarios. But in the same experiment scenarios, the response speed by H₂ controller is 8, 12, and 15 s, respectively. Less frequency deviation and faster response speed bring the better system control performance.

Evaluated controller performance in the stochastic transmitting latency

More general scenarios coinciding with the real DSN in CPPS are employed in the experiments, in which the transmitting latency is stochastic with three different intervals [0, 0.2 s], [0, 0.6 s], [0, 1 s]. The system frequency response with different controller is shown in Figure 9.

Figure 9.

Frequency response with stochastic transmitting latency: (a) frequency response when τ = 0 ~ 0.2 s, (b) frequency response when τ = 0 ~ 0.6 s, and (c) frequency response when τ = 0 ~ 1 s.

The results from Figure 9 show that, compared with the H₂ controller, the proposed state feedback controller also has better frequency dynamic response performance in the stochastic transmitting latency scenarios. Analyzing the results in detailed, the maximum frequency deviation kept to 0.05 Hz, compared with the 0.1 Hz by H₂ controller. The response speeds are 5, 9, and 12 s in different latency intervals. In the same experiment scenarios, the response speed by H₂ controller is 8 s in [0 0.2 s] latency interval, and great longer than 15 s with the [0 0.6 s] and [0 1 s] latency intervals. Long-term frequency oscillation will damage the electrical equipment and induce the cascade failures, which cannot be tolerated in the electrical power system.

Conclusion

CPPS with DSNs technique is widely used in electric power system. Cyber failure will induce the physical malfunction, even power system breakdown. Therefore, new network control theory and feedback controller should be studied to increase the capability of system robustness. In this article, state feedback controller aware stochastic transmitting latency is proposed. The heuristic model with continuous-time power system and discrete-time DSNs is established. A new state feedback controller co-operated Kalman observer is presented based on the stochastic optimal control theory. The experiment results show that the proposed controller effectively reduces the frequency deviation and obtain a faster response speed, which is beneficial to increase the capability of system’s frequency stability and improve the power supply quality.

Footnotes

Handling Editor: Christos Anagnostopoulos

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the National Key Research and Development Program of China (2017YF B0903000) and the National Natural Science Foundation of China (61571324).

References

Arghandeh

Meier

Mehrmanesh

et al . On the definition of cyber-physical resilience in power systems. Renew Sustain Energ Rev 2016; 58: 1060–1069.

Rahman

Mahmud

AMT

et al . Multi-agent approach for enhancing security of protection schemes in cyber-physical energy systems. IEEE T Ind Inform 2017; 13(2): 436–447.

Yang

Huang

Pen

et al . Optimal planning of communication system of CPS for distribution network. J Sensor 2017; 2: 1–10.

Ahmed

MA.

Wireless communication architectures based on data aggregation for internal monitoring of large-scale wind turbines. Int J Distrib Sensor Netw 2016; 12(8): 14.

Campobello

Segreto

Serrano

Data gathering techniques for wireless sensor networks: a comparison. Int J Distrib Sensor Netw 2016; 12(1): 1–17.

Chen

Liu

KY.

Network time-delay compensation based on Narandra stable adaptive control. Adv Mat Res 2014; 971–973: 1266–1271.

Bai

Sun

Gao

et al . Modeling and stabilization of a wireless network control system with packet loss and time delay. J Franklin Inst 2012; 349(7): 2420–2430.

Dong

Tian

Han

QLA

et al . A delay system method for designing event-triggered controllers of networked control systems. IEEE T Automat Contr 2013; 58(2): 475–481.

Pang

Liu

Zhou

et al . Output tracking control for networked systems: a model-based prediction approach. IEEE T Ind Electron 2014; 61(9): 4867–4877.

10.

Shi

Robust mixed H2/H∞control of networked control systems with random time delays in both forward and backward communication links. Automatica 2011; 47(4): 754–760.

11.

Nayeripour

Hoseintabar

Niknam

Frequency deviation control by coordination control of FC and double-layer capacitor in an autonomous hybrid renewable energy power generation system. Renew Energ 2011; 36(6): 1741–1746.

12.

Bevrani

Hiyama

On load–frequency regulation with time delays: design real-time implementation. IEEE T Energy Conver 2009; 24(1): 292–300.

13.

Duan

FY.

Networked AGC robust controller design in consideration of communication delay. Proc CSEE 2006; 26(22): 35–40.

14.

Teng

Gao

Zhang

et al . Key technologies and the implementation of wind, PV and storage co-generation monitoring system. J Mod Power Syst Cle 2014; 2(2): 104–113.

15.

Ding

Gao

Zhang

et al . Control of hybrid AC/DC microgrid under islanding operational conditions. J Mod Power Syst Cle 2014; 2(3): 223–232.

16.

Yunwei

Nejabatkhah

Overview of control, integration and energy management of microgrids. J Mod Power Syst Cle 2014; 2(3): 212–222.

17.

Oureilidis

Bakirtzis

Demoulias

CS.

Frequency-based control of islanded microgrid with renewable energy sources and energy storage. J Mod Power Syst Cle 2016; 4(1): 54–62.

18.

Yang

Chen

et al . Frequency hierarchical control for islanded micro-grid consisting of diesel generator and battery energy storage system. Control Theory Appl 2015; 32(8): 1098–1105.

19.

Yang

Wang

Frequency optimization control for isolated photovoltaic-diesel hybrid microgrid based on load estimation. Proc CSEE 2013; 33(34): 115–121.

20.

Kuroda

Magori

Ichimura

et al . A hybrid multi-objective optimization method considering optimization problems in power distribution systems. J Mod Power Syst Cle 2015; 3(1): 41–50.

21.

Meng

Yunfei

Jianzhong

et al . Dynamic frequency response from electric vehicles in the Great Britain power system. J Mod Power Syst Cle 2015; 3(2): 203–211.

22.

Matni

Communication delay co-design in H2 distributed control using atomic norm minimization. IEEE T Control Network Syst 2017; 4(2): 267–278.