A Cyber Physical Model Based on a Hybrid System for Flexible Load Control in an Active Distribution Network

To strengthen the integration of the primary and secondary systems, a concept of Cyber Physical Systems (CPS) is introduced to construct a CPS in Power Systems (Power CPS). The most basic work of the Power CPS is to build an integration model which combines both a continuous process and a discrete process. The advanced form of smart grid, the Active Distribution Network (ADN) is a typical example of Power CPS. After designing the Power CPS model architecture and its application in ADN, a Hybrid System based model and control method of Power CPS is proposed in this paper. As an application example, ADN flexible load is modeled and controlled with ADN feeder power control by a control strategy which includes the normal condition and the underpowered condition. In this model and strategy, some factors like load power consumption and load functional demand are considered and optimized. In order to make up some of the deficiencies of centralized control, a distributed control method is presented to reduce model complexity and improve calculation speed. The effectiveness of all the models and methods are demonstrated in the case study.


Introduction
With the rapid development of network technology, and increasing requirements for accuracy and efficiency of industry control, the industry needs to strengthen information acquisition, application and coordination control with the hierarchy-distributed control system.Hence, a new concept of Cyber Physical Systems (CPS) is proposed to promote the tight integration of control, communication, information, industry process, and reveal the integration mechanism of the real and virtual world [1].
CPS stems from embedded systems.To satisfy the performance requirement of spatial layout and power consumption, the traditional embedded system focuses on the optimal resource utilization of hardware and software [2].However, the embedded system in CPS concentrates on the interaction with the real physical system (controlled object), and focuses on problems such as concurrent event which cannot be perceived if the controlled object is too complicated [3].
Several research results on CPS have been presented.The purpose of them is to explore principle of interaction between cyber and physical systems, promote the system ability of sensing, analysis, decision and control.For example, when operating cardiac surgery by robot arms, the arms should vibrate with the heart at the same frequency to maintain relative rest [4].CPS models are built in [5,6] for fuel management and airspace management of aviation systems to optimize flight operation.
Energies 2017, 10, 267; doi:10.3390/en10030267www.mdpi.com/journal/energies Energies 2017, 10, 267 2 of 20 A power system meets the basic form of CPS.Especially after the smart grid and micro grid are presented, intellectualization and informatization cover almost all parts of the power system from source side to load side, and then the relationship and coordination among these parts are tighter.So the power system has had basis and actual demand to construct a complete CPS.
There are a few researches of CPS in the power system at present.Most of them focus on the influence brought from the information system (including the secondary system and the communication network) to the primary system.For example, in the aspect of reliability of relay protection, cyber factors like fault rates of logical nodes and logical connections in IEC61850 based substation automation system [7] and the fault rate of a communication link [8] are all considered in a traditional physical device reliability assessment.In the aspect of communication network, information system elements including communication routing [9], network time delay [10] and system observability [11] have been integrated into the primary system to analyze stability margin.There are also some researches about information security in the power system, such as information attack occurring in SCADA [12] and substation automation [13].To test the interaction between the information and primary system, several simulation platforms have been built such as INSPIRE [14], GECO [15], and TASSCS [16].
Besides the influence of information system, some researches consider that continuity and discrete are the main attributes of Cyber and Physical system respectively.So they try to rebuild power grid model according to Hybrid System theory which describes both the continuous dynamic and discrete state, then the object-oriented hybrid power system model library is developed [17], and sliding mode control based on the hybrid models can be applied to DG integrated shipboard power grid [18].
The main target of the Active Distribution Network (ADN) is to promote Distributed Generation (DG) utilization and reduce energy from the external system [19].In ADN, the control system coordinate with the primary system, complex energy flow, and information flow coexist when operating.Hence, ADN is suitable to be a study object for Power CPS.But, most of the existing ADN researches is on DGs control, and load control has still not been fully considered.The hierarchical control system in [20] is not only used in DG power control, but also in load control.
The load control is always a hot point for research.For instance, the relationship of refrigerators and the frequency in the micro grid are explored in [21].A Model Predictive Control (MPC) based method is proposed in [22] to optimize the profit of power demand in the commercial office building.From the point of CPS, literature [23] introduces a method of a Real-Time Physical System to reduce system instantaneous power consumption by adjusting the load coincidence factor.
There are several types of model which can be used in CPS to solve different problems.The Hybrid System is a typical one which can reflect CPS connotations, and can be used in control.This paper builds a Hybrid System based modeling structure for Power CPS which is beneficial for analyzing the interactive mechanism between information and the primary system.Then according to the modeling structure, an integration model method and its application are studied to optimally and predictively control the power system by an example of a flexible load control in ADN.
The rest of the paper is organized as follows: Section 2 describes the composition and characteristics of the Power CPS model, and designs application schemes in ADN.Section 3 presents a Hybrid System based CPS model and the control strategy for an ADN flexible load.Aiming at some shortages in Section 3, a distributed ADN flexible load control method is proposed in Section 4. All the models and methods are applied to simulated examples in Sections 3 and 4 respectively.Section 5 provides concluding remarks.

Power CPS Model and Its Application in ADN
Power CPS technical features and key techniques have been summarized in [24].Among these key techniques, the Cyber Physical model is the most basic work.Figure 1 shows the relationship between technical features and model in Power CPS.According to the three features in Figure 1, the Power CPS model can be classified into two kinds: "Control Model" and "Information Model".The "Information Model" corresponds to Feature (3), and this type of model mainly studies a uniform information description mode which can be applied in information exchange after combining with an information transmission structure including the model mapping method and communication protocol.
The "Control Model" reflects Features ( 1) and (2), and this type of model focuses on promoting the interaction and eliminating the heterogeneity between the primary and the information system in the power grid.

Control Model Based on Hybrid System
Compared with the Information Model, the 'Control Model' undertakes the main tasks of Power CPS control.As is shown in Figure 1, the 'Control Model' strengthens the integration of cyber and physical in two aspects.
(1) Integrating physical process and control process This aspect reflects the interaction between physical characteristic and the control input of primary devices.During the operation, the state variable and output of the primary device change with time continuously, but the control input can be seen as a kind of disturbance which interrupts and changes the current's continuous dynamic.No matter how the input is transmitted to the primary device, the input affects the device's physical process in a certain time interval, thus the continuous dynamic switches to new states discretely.Therefore, when the control input is being generated to control a device, continuous process of this device and discrete state switching caused by input should both be considered.
(2) Integrating of primary system and information system This aspect describes the interaction between the information system, control process and primary system.Where, the information system not only refers to the nodes which may be traversed by information flow, but also the sequence and time of different information flow.So the influence caused by the information system can be represented as topology and sequential logic, and incorporated in the 'integrating of physical process and control process' introduced above.Figure 2 illustrates the relationship between two aspects of 'integrating' in the 'Control Model'.
It can be concluded that the controlled primary system contains continuous process and state switching; the information system works discretely and logically.Continuity and discrete are the major attributes of the physical and cyber system respectively.So, Power CPS needs a common model method which is suitable for the physical and the cyber system.Moreover, the common model needs to be used in optimal control.According to the three features in Figure 1, the Power CPS model can be classified into two kinds: "Control Model" and "Information Model".The "Information Model" corresponds to Feature (3), and this type of model mainly studies a uniform information description mode which can be applied in information exchange after combining with an information transmission structure including the model mapping method and communication protocol.
The "Control Model" reflects Features ( 1) and ( 2), and this type of model focuses on promoting the interaction and eliminating the heterogeneity between the primary and the information system in the power grid.

Control Model Based on Hybrid System
Compared with the Information Model, the 'Control Model' undertakes the main tasks of Power CPS control.As is shown in Figure 1, the 'Control Model' strengthens the integration of cyber and physical in two aspects.
(1) Integrating physical process and control process This aspect reflects the interaction between physical characteristic and the control input of primary devices.During the operation, the state variable and output of the primary device change with time continuously, but the control input can be seen as a kind of disturbance which interrupts and changes the current's continuous dynamic.No matter how the input is transmitted to the primary device, the input affects the device's physical process in a certain time interval, thus the continuous dynamic switches to new states discretely.Therefore, when the control input is being generated to control a device, continuous process of this device and discrete state switching caused by input should both be considered.
(2) Integrating of primary system and information system This aspect describes the interaction between the information system, control process and primary system.Where, the information system not only refers to the nodes which may be traversed by information flow, but also the sequence and time of different information flow.So the influence caused by the information system can be represented as topology and sequential logic, and incorporated in the 'integrating of physical process and control process' introduced above.Figure 2 illustrates the relationship between two aspects of 'integrating' in the 'Control Model'.
It can be concluded that the controlled primary system contains continuous process and state switching; the information system works discretely and logically.Continuity and discrete are the major attributes of the physical and cyber system respectively.So, Power CPS needs a common model method which is suitable for the physical and the cyber system.Moreover, the common model needs to be used in optimal control.Several researches have used the mixed continuous and discrete model in power system optimization.The Mixed Integer Quadratic Programming (MIQP) models are built for unit commitment problem in [25,26] and are solved by the cut-and-branch method.A scenario based dynamic economic dispatch model is presented by [27] to optimize the operation cost of uncertain the renewable energy integrated power system.The MIQP model is also used in [28] to diminish power losses in the microgrid laboratory system with several kinds of DGs.
The Hybrid System is such a kind of model that combines the continuous dynamic and discrete state together.It is also able to be transformed to the predictive model and used in receding horizon optimization, and then MPC listed in Figure 1 can be applied based on the Hybrid System model.The Hybrid System based method in Power CPS is not only to optimize the primary system, but also to incorporate the information process into the traditional power system.

CPS Model and Control in ADN
The ADN control system is shown in Figure 3 which contains three control levels which are the Global Energy Manage System (Global EMS), Area coordination controllers, and Device coordination controllers.
The Global EMS monitors and gathers information of the ADN, and computes a long-time operation mode according to the electrical constraints.ADN is divided into several coordination control areas, and the Area coordination controllers are allocated to every area and control the area in a short time interval.The Device coordination controllers are connected to each primary device, and control these devices directly.At present, this control system has been applied in ADN feeder power control [29,30], voltage control [31], and feeder automation [32].
The ADN control system meets the features of Power CPS: (1) The cyber system (controllers and communication network) and the physical system (primary devices) are heterogeneous systems; (2) Controllers implement closed-loop control to primary devices, and there are discrete state switching and continuous dynamic in this process; (3) Information flow exists in the information and control system.
However, ADN does not realize the integration of cyber and physical because of the lack of a suitable integration model and control method.Several researches have used the mixed continuous and discrete model in power system optimization.The Mixed Integer Quadratic Programming (MIQP) models are built for unit commitment problem in [25,26] and are solved by the cut-and-branch method.A scenario based dynamic economic dispatch model is presented by [27] to optimize the operation cost of uncertain the renewable energy integrated power system.The MIQP model is also used in [28] to diminish power losses in the microgrid laboratory system with several kinds of DGs.
The Hybrid System is such a kind of model that combines the continuous dynamic and discrete state together.It is also able to be transformed to the predictive model and used in receding horizon optimization, and then MPC listed in Figure 1 can be applied based on the Hybrid System model.The Hybrid System based method in Power CPS is not only to optimize the primary system, but also to incorporate the information process into the traditional power system.

CPS Model and Control in ADN
The ADN control system is shown in Figure 3 which contains three control levels which are the Global Energy Manage System (Global EMS), Area coordination controllers, and Device coordination controllers.
The Global EMS monitors and gathers information of the ADN, and computes a long-time operation mode according to the electrical constraints.ADN is divided into several coordination control areas, and the Area coordination controllers are allocated to every area and control the area in a short time interval.The Device coordination controllers are connected to each primary device, and control these devices directly.At present, this control system has been applied in ADN feeder power control [29,30], voltage control [31], and feeder automation [32].
The ADN control system meets the features of Power CPS: (1) The cyber system (controllers and communication network) and the physical system (primary devices) are heterogeneous systems; (2) Controllers implement closed-loop control to primary devices, and there are discrete state switching and continuous dynamic in this process; (3) Information flow exists in the information and control system.However, ADN does not realize the integration of cyber and physical because of the lack of a suitable integration model and control method.Problem (1) is related to the primary system modeling.The ADN primary system model may change if state events disturb continuous operation.Now, ADN computes global target by long time optimization, and the optimal model only describes the dynamic of a certain state, but the optimal switching of the multi-state is not considered.
Being one of the Hybrid System models, Mixed Logical Dynamical (MLD) describes objects by dynamic and logic rules, and puts the objects' discrete and continuous characteristics into a unified framework [33].If the ADN primary system is modeled by MLD, all the separated operation states can be optimized together.
Just like Problem (1), Problem (2) is also caused by the model, and can be solved by MLD.According to route and sequence rules, information flow transmits in the control system.If these rules are modeled by logical expression, and transformed to inequality constraints by MLD, then the information flow and also the structure of the control system will be involved in the control strategy.
Problem (3) is associated with the control method.The control period of ADN long time optimization is about 10-15 min.However, the system finds it hard to maintain an optimal target because of so many variable factors.So, it is necessary to predict the state of the controlled object and reduce the control period.Moreover, if some disturbances take place, the control model should be rectified in the next period.Hence, MLD based MPC can be used to solve Problem (3).
Obviously, MLD and the control method based on it are able to diminish the ADN's existing weakness.As a preliminary study of the Power CPS model, this paper pays attention to the ADN primary device control which corresponds to the "Integrating physical and control processes", and lays the foundation for future work.Problem (1) is related to the primary system modeling.The ADN primary system model may change if state events disturb continuous operation.Now, ADN computes global target by long time optimization, and the optimal model only describes the dynamic of a certain state, but the optimal switching of the multi-state is not considered.

ADN Flexible Load Control Based on the Hybrid System Model
Being one of the Hybrid System models, Mixed Logical Dynamical (MLD) describes objects by dynamic and logic rules, and puts the objects' discrete and continuous characteristics into a unified framework [33].If the ADN primary system is modeled by MLD, all the separated operation states can be optimized together.
Just like Problem (1), Problem (2) is also caused by the model, and can be solved by MLD.According to route and sequence rules, information flow transmits in the control system.If these rules are modeled by logical expression, and transformed to inequality constraints by MLD, then the information flow and also the structure of the control system will be involved in the control strategy.
Problem (3) is associated with the control method.The control period of ADN long time optimization is about 10-15 min.However, the system finds it hard to maintain an optimal target because of so many variable factors.So, it is necessary to predict the state of the controlled object and reduce the control period.Moreover, if some disturbances take place, the control model should be rectified in the next period.Hence, MLD based MPC can be used to solve Problem (3).
Obviously, MLD and the control method based on it are able to diminish the ADN's existing weakness.As a preliminary study of the Power CPS model, this paper pays attention to the ADN primary device control which corresponds to the "Integrating physical and control processes", and lays the foundation for future work.

Flexible Load Control Strategy
There are usually some high-power loads in the distribution network such as conditioning, cold storage, and the pump station.These loads need to achieve a function such as controlling temperature or water level.Two features can be concluded from them: (1) Operate segmentally in accordance with a different state; (2) Load functional demand (temperature, water level, etc.) is accumulative during the operation period.Loads meet these two features and are suitable to be modeled by MLD and applied MPC.
The optimal target of flexible load is to reduce power consumption on the premise of satisfying load functional demand.But in the perspective of the system, load consumption needs to be controlled near a stable target.
ADN feeder power control is affected by loads to a great extent.The control purpose of ADN feeder power is to sustain exchange power between control area and feeder (named: Area Power), also the exchange power between the feeder and the external system (named: Feeder Power).Uncontrolled load may switch state optionally in a short time and make exchange power fluctuate frequently and sharply, so the ADN feeder is apt to deviate from the operation target.
If flexible loads are incorporated into ADN control, and feeder power distribution adjusted with DGs, the feeder and loads will both share the benefit.
As shown in Table 2, ADN flexible control is classified into two conditions.When the feeder works normally, in every area, the DGs power satisfies the loads' demand.Based on global optimal computation, Feeder Power and Area Power should be maintained near the optimal point.Under this condition: (1) Each area coordination controller controls the flexible loads of this area to run at a minimum consumption mode; (2) All the DGs in this area sustain the area's Area Power.When a bigger power shortage occurs in the feeder (for example when some DGs quit because of the weather), Area Power of this feeder severely deviates from the target, so the main purpose at this time is to reduce Feeder Power deviation by flexible load coordination control.Under this condition: (1) The global EMS calculates the operation mode for all the flexible loads in the feeder, and controls the Feeder Power near its target value as well as being more stable; (2) The rest of the DGs are controlled to stabilize Feeder Power fluctuation.

Hybrid System Model of Flexible Load
The flexible load which meets the two features in 3.1 is modeled by MLD in this section. (

1) State model
The state model describes the load's functional demand, control input and the relation with outside disturbance.After discretization, n flexible loads can be expressed as a state function as Equation (1).
where x(t) = [x 1 (t), x 2 (t) . . .x n (t)] is the functional index which reflects the loads' functional demand at time t.u(t) = [u 1 (t), u 2 (t) . . .u n (t)] T is the control input, and the element u n (t) ∈ U = {u n1 , u n2 . . .u ni } denotes that each load contains i kinds of input and can only chose one of these inputs to control the load at time t.A, B 2 are n × n coefficient matrixes of the x(t) and u(t), they are determined by load type.B 1 (t) = [B 11 (t), B 12 (t) . . .B 1n (t)] T is the outside disturbance which will affect the function index at time t + 1.
(2) Power consumption model This model describes the loads' power under different control inputs.Equation ( 2) is the discrete power consumption model.
where y(t) denotes the total power of n flexible loads at time t.D = [p 1 (t), p 2 (t) . . .p n (t)] consists of the power of each load at time t, and p n (t) ∈ P = {p n1 , p n2 . . .p ni }.The elements of P have a one-to-one correlation with U and denote the power of i-th control inputs of n flexible loads.
(3) MLD model If the flexible load of Equation ( 1) contains a finite control input selection, then it can be considered as states switching.By setting logical variables, the separate states can be integrated.
During a time step ∆t, n loads work under different states, so there are j = i n kinds of operation modes that u(t) = u m = [u 1 (t), u 2 (t) . . .u n (t)] T and m ∈ {1, 2 . . .j}.According to Equation (3), set a logical variable for each of the j operation modes respectively, and make up a column vector δ(t) = [δ 1 (t), δ 2 (t) . . .δ m (t) . . .δ j (t)] T , where δ(t) ∈ {0, 1}, m ∈ {1, 2 . . .j}.Each element of δ(t) means a state combination of the loads' operation mode.When the loads work under this state combination, set the corresponding element of δ(t) to 1, that is Based on Equation (3), the Equation ( 1) is denoted by δ(t), and then a n × j matrix B 2t = {B 2t1, B 2t2 . . .B 2ti } is deduced.The elements (they are all n-dimension column vectors) of B 2t is calculated by B 1 (t) + B 2 •u m which represents the influence factors brought to a functional index which are caused by load disturbance or load control under the operation mode δ m (t).
Replace the control input u(t) of Equation ( 2) with δ(t), then the meaning of D changes from the description of each load's power to the description of n loads' total power under j state combinations.So set a j-dimension row vector D 2t = [P 1 , P 2 . . .P m . . .P j ] which consists of j possible total powers corresponding to j kinds of state combinations, where At any time t, all the flexible loads can work under only one combination, so there is only one element of δ(t) can be set to 1.This constrain is written as Equation ( 4), where 1 is the n-dimension vector of only 1.
With Equations ( 1)-(4), and considering the limit of functional index x(t) ∈ [a, b], a MLD based flexible load model Equation ( 5) is acquired, where 0 n is a n-dimension zero matrix and I n is a n-dimension identity matrix.

MPC of Flexible Load
The MPC problem is to obtain a series of control inputs which satisfy optimal function Equation (6) and meet the constraint conditions Equation ( 5).The Q δ , Q x , Q y in Equation ( 6) is weight factors of logical variable, functional index, and total power consumption.
ADN load control has been classified into Normal Condition and Underpowered Condition in 3.1, and the conditions are distinguished in accordance with the feeder power shortage ∆P(t) of Equation (7).
where (t|t  Figure 4 shows the control procedure of ADN flexible load.The ∆P(t) will be computed firstly to judge whether its value exceeds a preset threshold P threshold .If ∆P(t) ≤ P threshold , the feeder is in the Normal Condition.All the flexible loads run independently to reduce power consumption in a control period T, and Equation ( 5) should set target points as δ f = δ min , δ min → (P min = 0) or y f = P min = 0.

MPC of Flexible Load
The MPC problem is to obtain a series of control inputs which satisfy optimal function Equation (6) and meet the constraint conditions Equation (5).The Qδ, Qx, Qy in Equation ( 6) is weight factors of logical variable, functional index, and total power consumption.
ADN load control has been classified into Normal Condition and Underpowered Condition in 3.1, and the conditions are distinguished in accordance with the feeder power shortage ΔP(t) of Equation (7).  Figure 4 shows the control procedure of ADN flexible load.The ΔP(t) will be computed firstly to judge whether its value exceeds a preset threshold Pthreshold.If ΔP(t) ≤ Pthreshold, the feeder is in the Normal Condition.All the flexible loads run independently to reduce power consumption in a control period T, and Equation ( 5) should set target points as δf = δmin, δmin → (Pmin = 0) or yf = Pmin = 0.After condition judgment and target setting, the optimal control models of Equations ( 5) and ( 6) are transformed to a Mixed Integer Quadratic Programming (MIQP) problem.Then the problem is solved and only the first sample of the optimal sequence is applied to the flexible loads by terminal controllers before the next control period arrives.
The flexible load model is able to be updated before the next period.Updating requests may be proposed from feeder side or load side.For instance, the feeder or load both can adjust the limit of the functional index x(t), and this will lead to a negotiation mechanism between the two sides.

Case Study
The MLD based CPS model of ADN flexible load and its MPC method is tested using the 10 kV feeder system shown in Figure 5.
If ΔP(t) ≥ Pthreshold, the area coordination control will be closed at once.All the flexible loads, DGs in the feeder are controlled to mitigate fluctuation and diminish deviation of the Feeder Power.Then choose a consumption target Pm from the flexible loads' achievable power set D2t.The chosen Pm  n m max ctrlLoad m m P P P t to be minimal, and it means that the flexible loads try their best to counteract the power shortage.
After condition judgment and target setting, the optimal control models of Equations ( 5) and ( 6) are transformed to a Mixed Integer Quadratic Programming (MIQP) problem.Then the problem is solved and only the first sample of the optimal sequence is applied to the flexible loads by terminal controllers before the next control period arrives.
The flexible load model is able to be updated before the next period.Updating requests may be proposed from feeder side or load side.For instance, the feeder or load both can adjust the limit of the functional index x(t), and this will lead to a negotiation mechanism between the two sides.

Case Study
The MLD based CPS model of ADN flexible load and its MPC method is tested using the 10 kV feeder system shown in Figure 5.The feeder contains a main line and three branches, and it is divided into three coordination control areas (dashed boxes).There are five DGs (including three green energies and two energy storages) and three flexible loads distributed in the three branches.There are also five common loads locate in Node 1, 4, 5, 9, 12.The global optimization determines the feeder's operation targets for a long period.Table 3 shows the main parameters and targets of the feeder system.The feeder contains a main line and three branches, and it is divided into three coordination control areas (dashed boxes).There are five DGs (including three green energies and two energy storages) and three flexible loads distributed in the three branches.There are also five common loads locate in Node 1, 4, 5, 9, 12.The global optimization determines the feeder's operation targets for a long period.Table 3 shows the main parameters and targets of the feeder system.The three flexible loads are all middle-type cold storages which keep inner temperatures within the limits by electric refrigeration equipment.The inner temperatures of the three cold storages can be modeled as Equation (1), where x(t) is inner temperature at time t; B 1 (t) indicates temperature disturbance caused by outside heat, and it is invariant in this example; B 2 denotes the refrigeration effect coefficient matrix; the control input u(t) has two states including Stop and Start, the two states represented by 0 and 1 respectively, that is u n (t) ∈ U = {0, 1}, n ∈ {1, 2, 3}; set the control step ∆t = 2 s.Table 4 lists the typical parameters of the three types of cold storages.Now the three flexible loads' control model as Equation ( 5) under the Normal and Underpowered condition can be obtained. (

1) Normal condition control
The control purpose is to minimize the load's consumption of each area.Set control period T = 10∆t (20 s), and simulation period is 0-200 s, so there are 100 times MPC computations during this period.
From Figure 6, it can be found that the flexible load temperature controlled by the model and method of this paper (flexible group) is better than the loads controlled by the traditional method (traditional group).In Figure 6a, at any time, three temperatures are all restricted in the low and up boundaries and the curves seems more stable.This is because the load models predict the future trend and avoid the possibility of out-of-limit.In addition, Equation (6) contains reduction of the deviation of x(t), and this contributes to the stability.
The traditional group is controlled by simple if-then modes.Because of the lack of prediction, the controllers do not act only when the x(t) exceeds the limit.Moreover, by this mode, x(t) always occupies the whole range of the limit and leads to sharp fluctuation.The three flexible loads are all middle-type cold storages which keep inner temperatures within the limits by electric refrigeration equipment.The inner temperatures of the three cold storages can be modeled as Equation (1), where x(t) is inner temperature at time t; B1(t) indicates temperature disturbance caused by outside heat, and it is invariant in this example; B2 denotes the refrigeration effect coefficient matrix; the control input u(t) has two states including Stop and Start, the two states represented by 0 and 1 respectively, that is un(t) ∈ U = {0, 1}, n ∈ {1, 2, 3}; set the control step Δt = 2 s.Table 4 lists the typical parameters of the three types of cold storages.5) under the Normal and Underpowered condition can be obtained. (

1) Normal condition control
The control purpose is to minimize the load's consumption of each area.Set control period T = 10Δt (20 s), and simulation period is 0-200 s, so there are 100 times MPC computations during this period.
From Figure 6, it can be found that the flexible load temperature controlled by the model and method of this paper (flexible group) is better than the loads controlled by the traditional method (traditional group).In Figure 6a, at any time, three temperatures are all restricted in the low and up boundaries and the curves seems more stable.This is because the load models predict the future trend and avoid the possibility of out-of-limit.In addition, Equation (6) contains reduction of the deviation of x(t), and this contributes to the stability.
The traditional group is controlled by simple if-then modes.Because of the lack of prediction, the controllers do not act only when the x(t) exceeds the limit.Moreover, by this mode, x(t) always occupies the whole range of the limit and leads to sharp fluctuation.In Figure 7, loads of the flexible group act more times than the traditional group.This phenomenon implies that load devices in the flexible group need higher reliability, and it can be adjusted by weight factors for different demand.In Figure 7, loads of the flexible group act more times than the traditional group.This phenomenon implies that load devices in the flexible group need higher reliability, and it can be adjusted by weight factors for different demand.Table 5 indicates the power consumption of two groups.It is obvious that each load's energy consumption and mean power of flexible group is lower than that of the traditional group.However, because the flexible loads are not optimized together under the normal condition, loads in different control areas do not improve the load coincidence factor.8 reflects the Feeder Power and 3 Area power.Every curve involves power pulse which is due to the loads' frequent switching.The power pulses also imply that when load disturbance occurs, the ADN feeder power control will change the DGs' output for the corresponding area and the external input to counteract the disturbance and maintain the optimal target.(2) Underpowered condition control If at time t = 1, because of the weather, DG1 reduces its output 215 kW, and DG2 reduce 250 kW, there is a power shortage of ΔP = 365 kW in the feeder according to Equation (7).Under this condition, all the flexible loads are controlled, unified, and undertake some of the power shortage.In this case, Table 5 indicates the power consumption of two groups.It is obvious that each load's energy consumption and mean power of flexible group is lower than that of the traditional group.However, because the flexible loads are not optimized together under the normal condition, loads in different control areas do not improve the load coincidence factor.Figure 8 reflects the Feeder Power and 3 Area power.Every curve involves power pulse which is due to the loads' frequent switching.The power pulses also imply that when load disturbance occurs, the ADN feeder power control will change the DGs' output for the corresponding area and the external input to counteract the disturbance and maintain the optimal target.Table 5 indicates the power consumption of two groups.It is obvious that each load's energy consumption and mean power of flexible group is lower than that of the traditional group.However, because the flexible loads are not optimized together under the normal condition, loads in different control areas do not improve the load coincidence factor.8 reflects the Feeder Power and 3 Area power.Every curve involves power pulse which is due to the loads' frequent switching.The power pulses also imply that when load disturbance occurs, the ADN feeder power control will change the DGs' output for the corresponding area and the external input to counteract the disturbance and maintain the optimal target.(2) Underpowered condition control If at time t = 1, because of the weather, DG1 reduces its output 215 kW, and DG2 reduce 250 kW, there is a power shortage of ΔP = 365 kW in the feeder according to Equation (7).Under this condition, all the flexible loads are controlled, unified, and undertake some of the power shortage.In this case, (2) Underpowered condition control If at time t = 1, because of the weather, DG1 reduces its output 215 kW, and DG2 reduce 250 kW, there is a power shortage of ∆P = 365 kW in the feeder according to Equation (7).Under this condition, all the flexible loads are controlled, unified, and undertake some of the power shortage.In this case, the power target of the three loads is 356 kW, set control step ∆t = 2 s, control period T = 10∆t, and simulation period is 0-200 s.
The same as the normal condition, because of the strict constraints and prediction in Equation ( 5), the temperatures of the three cold storages of the flexible group in the underpowered condition are all in the feasible range, and more stable than the traditional group.(Shown in Figure 14a).
Figure 9 and Table 6 compare the total power of the two groups.From Figure 9, the flexible group undertakes the power shortage effectively, and operates very close to the target of 356 kW and the RMS value deviates from 356 kW only about 30% which is lower than the traditional group.Additionally, the flexible group avoids a sharp power change, and the maximal load difference between the two steps is only 370 kW.However, in the traditional group, the power changes sharply and irregularly.Obviously, the phenomenon in the traditional group is negative for ADN feeder power.From the point of view of energy conservation, the flexible group never works at the peak power of 726 kW, and the total energy is 4.5% lower than the traditional group.
Energies 2017, 10, 267 12 of 20 the power target of the three loads is 356 kW, set control step Δt = 2 s, control period T = 10Δt, and simulation period is 0-200 s.
The same as the normal condition, because of the strict constraints and prediction in Equation ( 5), the temperatures of the three cold storages of the flexible group in the underpowered condition are all in the feasible range, and more stable than the traditional group.(Shown in Figure 14a).Figure 9 and Table 6 compare the total power of the two groups.From Figure 9, the flexible group undertakes the power shortage effectively, and operates very close to the target of 356 kW and the RMS value deviates from 356 kW only about 30% which is lower than the traditional group.
Additionally, the flexible group avoids a sharp power change, and the maximal load difference between the two steps is only 370 kW.However, in the traditional group, the power changes sharply and irregularly.Obviously, the phenomenon in the traditional group is negative for ADN feeder power.From the point of view of energy conservation, the flexible group never works at the peak power of 726 kW, and the total energy is 4.5% lower than the traditional group.The Feeder Power curve is shown in Figure 10.It is easy to see that the flexible group keeps near the target.However, the traditional group has loss of stability, the frequent and large power change makes the Feeder Power oscillate severely, and the target value cannot be held.The Feeder Power curve is shown in Figure 10.It is easy to see that the flexible group keeps near the target.However, the traditional group has loss of stability, the frequent and large power change makes the Feeder Power oscillate severely, and the target value cannot be held.
Energies 2017, 10, 267 12 of 20 the power target of the three loads is 356 kW, set control step Δt = 2 s, control period T = 10Δt, and simulation period is 0-200 s.
The same as the normal condition, because of the strict constraints and prediction in Equation ( 5), the temperatures of the three cold storages of the flexible group in the underpowered condition are all in the feasible range, and more stable than the traditional group.(Shown in Figure 14a).Figure 9 and Table 6 compare the total power of the two groups.From Figure 9, the flexible group undertakes the power shortage effectively, and operates very close to the target of 356 kW and the RMS value deviates from 356 kW only about 30% which is lower than the traditional group.
Additionally, the flexible group avoids a sharp power change, and the maximal load difference between the two steps is only 370 kW.However, in the traditional group, the power changes sharply and irregularly.Obviously, the phenomenon in the traditional group is negative for ADN feeder power.From the point of view of energy conservation, the flexible group never works at the peak power of 726 kW, and the total energy is 4.5% lower than the traditional group.The Feeder Power curve is shown in Figure 10.It is easy to see that the flexible group keeps near the target.However, the traditional group has loss of stability, the frequent and large power change makes the Feeder Power oscillate severely, and the target value cannot be held.

Distributed Load Control in the Underpowered Condition
Section 3 introduces a Hybrid System based CPS model and control method for the flexible load.The model and method promote load participation in ADN control, improve load consumption and ADN operation, and form positive interaction of loads and system.
However, there are also three aspects to the problems in the unified load control of the underpowered condition.
(1) Complex control model of the large scale system.When all loads are controlled together, there are too many possible state switching and combinations, this leads to a huge model for Equation (5).
After the model is transformed to MIQP, the scale and complexity of the calculation model expand with the increase of the control period T.
(2) Non-convergence caused by a distinction of model objects.There are always some differences between different kinds of objects like the model parameter and the control period.If these objects are modeled together, the control problem may last too long or even be unsolvable.(3) The interaction between area coordination controllers is insufficient.The modeling and controlling mostly rely on global EMS, and do not take full advantage of the ability of lower-level controllers.
In this section, the distributed model and the MPC control method for flexible load, based on Section 3, are proposed to help with solving the three problems above.

Distributed MPC in ADN
The distributed control usually has two types.The first one; the control problem is split into several parts at the boundary point.These sub-problems are solved by alternating iteration.The second one; the control problem is separated into several coupled sub-problems according to the arrangement of controllers and category of the controlled object.These coupled sub-problems are solved by transmitting coupling information [34].
Considering the calculation speed and the existing control system, the second type is more suitable for ADN. Figure 11 shows the differences between distributed control and centralized control in Section 3 of the ADN flexible load when power shortage happens in the feeder.There are three main differences.

Distributed Load Control in the Underpowered Condition
Section 3 introduces a Hybrid System based CPS model and control method for the flexible load.The model and method promote load participation in ADN control, improve load consumption and ADN operation, and form positive interaction of loads and system.
However, there are also three aspects to the problems in the unified load control of the underpowered condition.
(1) Complex control model of the large scale system.When all loads are controlled together, there are too many possible state switching and combinations, this leads to a huge model for Equation (5).After the model is transformed to MIQP, the scale and complexity of the calculation model expand with the increase of the control period T. (2) Non-convergence caused by a distinction of model objects.There are always some differences between different kinds of objects like the model parameter and the control period.If these objects are modeled together, the control problem may last too long or even be unsolvable.(3) The interaction between area coordination controllers is insufficient.The modeling and controlling mostly rely on global EMS, and do not take full advantage of the ability of lowerlevel controllers.
In this section, the distributed model and the MPC control method for flexible load, based on Section 3, are proposed to help with solving the three problems above.

Distributed MPC in ADN
The distributed control usually has two types.The first one; the control problem is split into several parts at the boundary point.These sub-problems are solved by alternating iteration.The second one; the control problem is separated into several coupled sub-problems according to the arrangement of controllers and category of the controlled object.These coupled sub-problems are solved by transmitting coupling information [34].
Considering the calculation speed and the existing control system, the second type is more suitable for ADN. Figure 11 shows the differences between distributed control and centralized control in Section 3 of the ADN flexible load when power shortage happens in the feeder.There are three main differences.Difference (1): Control structure.The centralized control only has two control levels, global EMS and terminal controller, and all the information is exchanged between these two levels.The distributed mode contains three levels, and the communication not only appears between levels, but also appears between controllers in the same level.
Difference (2): Control task.The hollow circles in Figure 11 indicate the main tasks of each control level.In centralized mode, the global EMS undertakes almost all the modeling and computation but in distributed mode, some tasks are assigned to the area coordination control level.Difference (1): Control structure.The centralized control only has two control levels, global EMS and terminal controller, and all the information is exchanged between these two levels.The distributed mode contains three levels, and the communication not only appears between levels, but also appears between controllers in the same level.Difference (2): Control task.The hollow circles in Figure 11 indicate the main tasks of each control level.In centralized mode, the global EMS undertakes almost all the modeling and computation but in distributed mode, some tasks are assigned to the area coordination control level.
Difference (3): Control information.The solid circles in Figure 11 show the main information exchanged in the control process.The distributed model adds coupling information transmission of the same level, and because the coupling information is not generated in real time, there is therefore some information of a different time.
Even though, in distributed control, the optimal performance may be lost and the complexity of communication is increased, the control model is simplified significantly and computation time is reduced, also the flexibility of the control is strengthened.

Distributed Control Model and Method of Flexible Load
When the centralized control solves the optimal model, the coupling relationship between controlled objects exists naturally during the iteration.Only if all the variables including coupling variables satisfy the constraints and optimal purpose, is then the feasible solution acquired.However, the coupling is impossible to automatically appear in the distributed control process, because the sub-problems are separated; additionally, independent control makes it that every controller cannot obtain accurate and real-time state information of the coupled controllers.These two problems are solved by reconstructing the models of Section 3. (

1) Classifying Load Controlled Object
The Load Controlled Object is a flexible loads set, and all the flexible loads of the set are controlled with the same controller.The classification relates to the load type and area, in general, loads of the same type and with similar parameter or the same time step are placed into a category.Let (0, 1) s×n express the set of n × n matrices which only consists of elements 0 and 1. Introduce an incidence matrix A I ∈ (0, 1) s×n which defines relation of n loads of s classes in ADN.Only when the matrix's element A I (i,j) = 1, i = 1, . . ., s, j = 1, . . .,n, the j-th load belongs to the i-th class LCO i .
(2) Coupling relationship of LCOs The coupling relationship means that several LCOs affect the optimal target, and this relationship determines if information exchange exists between controllers of the same level.Introduce an incidence vector A I I−q ∈ (0, 1) 1×s which defines the relation between the q-th LCO q and other s − 1 LCOs.If the q-th LCO q couples with the p-th LCO p , the vector's elements A II-q (1,p) = 1, p = 1, . . .,s, p = q, and A II-q (1,q) = 0.
(3) Coupling variable In Equation (6), every flexible load contributes to the optimization of total power y(t), and LCOs are the same after the n loads are classified into s LCOs.Hence, the s LCOs has a coupling relation to However, each LCO is controlled by the corresponding controller, and the q-th controller Ctrl q knows the model of q-th LCO q including the state model (Equation ( 1)) and the power consumption model (Equation ( 2)).Moreover, the Ctrl q does not consider the rest of s − 1 LCOs.Because of this, other controllers such as the p-th controller Ctrl p must offer power information at time t to the Ctrl q , that is y(t) = y 1|t (t) + • • • + y q|t (t) + • • • + y s|t (t), where y p|t is the knowable coupling variable, p = 1, . . ., s, and p = q.
Unfortunately, at time t, due to the same reason of lack of coupling information, y p|t is unable to be obtained.As shown in Figure 12, the coupling information can be replaced by the predictive control results of [t 1 − 1, t 1 − 1 + T] which are calculated at time t = t 1 − 1, and T is the control period.However, for time t = t1 + T − 1, there is no yp(t1 + T − 1|t1 − 1) that can be assigned to yp|t1+T-1.Define constants yf-q and yf-p, p = 1, …, s, and p ≠ q. yf-q and yf-p consist of a power combination of all the LCOs.The power combination corresponds to a kind of load operation mode δm and satisfies No matter how the coupling variable gets its value, the Ctrlp may not follow the yp|t which has been sent to Ctrlq.So, there should be some restrictions in optimal function or constraint to keep the actual trajectory of LCO consistent with the coupling variables of the coupled controllers.
In optimal function, during a control period t = [t1, t1 + T], minimizing the deviation between yq(t) and yq|t, gives Equation (8).  y tT yq q q t tt J y t y Q (8) In control model, add inequality constraints, and make the deviation between yq(t) and yq|t of any step during the control period t = [t1, t1 + T] not larger than the maximal deviation in the last control period t = [t1 − 1, t1 − 1 + T].That is Equation (9).
q q t t t t T q q t t t t T y t y y t y (9) (5) Optimal function The optimal function of LCOq includes the functional index deviation Jx, load consumption deviation Jy, and coupling variable deviation Jy.
With the action of u q (t 1 |t 1 − 1), Ctrl q moves to time t = t 1 , and a new control loop begins.At this time, set y p|t = y p (t|t 1 − 1) for every step of t = [t 1 , However, for time t = t 1 + T − 1, there is no y p (t 1 + T − 1|t 1 − 1) that can be assigned to y p|t1+T-1 .Define constants y f-q and y f-p , p = 1, . . ., s, and p = q.y f-q and y f-p consist of a power combination of all the LCOs.The power combination corresponds to a kind of load operation mode δ m and satisfies No matter how the coupling variable gets its value, the Ctrl p may not follow the y p|t which has been sent to Ctrl q .So, there should be some restrictions in optimal function or constraint to keep the actual trajectory of LCO consistent with the coupling variables of the coupled controllers.
In optimal function, during a control period t = [t 1 , t 1 + T], minimizing the deviation between y q (t) and y q|t , gives Equation (8).
In control model, add inequality constraints, and make the deviation between y q (t) and y q|t of any step during the control period t = [t 1 , t 1 + T] not larger than the maximal deviation in the last control period (5) Optimal function The optimal function of LCO q includes the functional index deviation J x , load consumption deviation J y , and coupling variable deviation J y .
Like Equation (6), J x keeps x q (t) tracing the target x f-q in the control period t = [t 1 , t 1 + T], that is Equation (10).
Energies 2017, 10, 267 16 of 20 J y contains coupling relationship of LCO q and LCO p , and reflects the difference of total consumption and its target.Define matrix Y ∈ R s×1 , Y = [y 1|t , y 2|t , • • • , y s|t ] is the coupling information set of time t which should have been sent by s controllers to each other.Also the load consumption deviation of LCO q is shown as Equation (11).
J yq has been presented in Equation (8).So the optimal function is Equation (12).
(6) Control model and method Model the s LCOs respectively based on 3.2, and MLD models as Equation ( 5) are acquired.Then add Equations ( 9) and (12) to each LCO, and the control model is completed.
The optimization problem is that s LCOs are controlled by s controllers on the basis of coupling information to get control inputs.The solution must satisfy Equations ( 5) and ( 9), and make Equation (12) optimal.
The control procedure is shown in Figure 13.By the prediction of DGs and common loads, the ADN feeder is in the underpowered condition; every controller accepts the coupling information of the other controllers and the power consumption target from the global EMS; all the controllers generate control models of their own LCO and solve them, then act on the first sample of the optimal sequence; all the controllers send coupling information to other controllers and prepare for the next period of control.
Energies 2017, 10, 267 16 of 20 Jy contains coupling relationship of LCOq and LCOp, and reflects the difference of total consumption and its target.Define matrix y y y Y is the coupling information set of time t which should have been sent by s controllers to each other.Also the load consumption deviation of LCOq is shown as Equation (11).II q m tt J y t P Q AY (11) Jyq has been presented in Equation (8).So the optimal function is Equation (12).
(6) Control model and method Model the s LCOs respectively based on 3.2, and MLD models as Equation ( 5) are acquired.Then add Equations ( 9) and (12) to each LCO, and the control model is completed.
The optimization problem is that s LCOs are controlled by s controllers on the basis of coupling information to get control inputs.The solution must satisfy Equations ( 5) and ( 9), and make Equation (12) optimal.
The control procedure is shown in Figure 13.By the prediction of DGs and common loads, the ADN feeder is in the underpowered condition; every controller accepts the coupling information of the other controllers and the power consumption target from the global EMS; all the controllers generate control models of their own LCO and solve them, then act on the first sample of the optimal sequence; all the controllers send coupling information to other controllers and prepare for the next period of control.

Case Study
The distribution ADN flexible load control model and method of underpowered condition is tested by the case in 3.4.
According to 4.2, the three flexible loads of Load 1 , Load 2 , Load 3 are classified to three LOCs of LOC 1 , LOC 2 , LOC 3 .Based on Table 3 and Equations ( 5) and ( 9), three MLD models of LOCs are built.Set the total power consumption target P m = 356 kW, then the three optimal functions can be obtained by Equation (12).The models and optimal function are controlled 100 times according to the procedure in Figure 13.
Area coordination controllers are simulated by three IBM Flex System x240 Compute Nodes, and configure another Compute Node to simulate the physical system of the ADN feeder.Each Compute Node contains two Intel Xeon E5-2600 Core6 2.1 GHz processors and 16G memory.An inbuilt Gigabit Switch is used for communication.
Matlab/Simulink is the simulating and computing tool for every Compute Node, and a Socket interface is used for software to communicate with each other.There is also a synchronous program to prevent possible asynchronous problems occurring between different Compute Nodes.So the simulation of ADN feeder including three loads does not begin, only when the feeder has received all the load control inputs from other Compute Nodes.
Figure 14 shows the temperature of the distributed and centralized mode in the underpowered condition.The curves of the two modes all locate in the feasible range.The three loads of the distributed mode change their states about 56, 48, and 52 times respectively, and the loads in the centralized mode change only 48, 45, 46 times.The more frequent state switching makes the loads of the distributed mode have a lower fluctuation.Especially the Load 3 , it not only operates steadily, but also works at a lower temperature.
This phenomenon is caused by the method difference of how to balance the optimal target of the two modes.Load 3 is the largest load, and the centralized mode considers the optimization from the overall model, so the action of Load 3 should be diminished and the lower loads work more to counteract disturbance.However, in the distributed mode, the controller manages only one load, and does not need to consider other problems.
Figure 15 shows the state switching and total power of distributed control.From Figure 15b, the three loads' total consumption is able to be kept near 356 kW.But compared with Figure 9a, a sharp power change occurs in the distributed mode, and the curve is more unstable.This is caused by inaccurate coupling information between controllers.So, the optimal performance of the distributed mode is worse.The distribution ADN flexible load control model and method of underpowered condition is tested by the case in 3.4.
According to 4.2, the three flexible loads of Load1, Load2, Load3 are classified to three LOCs of LOC1, LOC2, LOC3.Based on Table 3 and Equations ( 5) and ( 9), three MLD models of LOCs are built.Set the total power consumption target Pm = 356 kW, then the three optimal functions can be obtained by Equation (12).The models and optimal function are controlled 100 times according to the procedure in Figure 13.
Area coordination controllers are simulated by three IBM Flex System x240 Compute Nodes, and configure another Compute Node to simulate the physical system of the ADN feeder.Each Compute Node contains two Intel Xeon E5-2600 Core6 2.1 GHz processors and 16G memory.An inbuilt Gigabit Switch is used for communication.
Matlab/Simulink is the simulating and computing tool for every Compute Node, and a Socket interface is used for software to communicate with each other.There is also a synchronous program to prevent possible asynchronous problems occurring between different Compute Nodes.So the simulation of ADN feeder including three loads does not begin, only when the feeder has received all the load control inputs from other Compute Nodes.
Figure 14 shows the temperature of the distributed and centralized mode in the underpowered condition.The curves of the two modes all locate in the feasible range.The three loads of the distributed mode change their states about 56, 48, and 52 times respectively, and the loads in the centralized mode change only 48, 45, 46 times.The more frequent state switching makes the loads of the distributed mode have a lower fluctuation.Especially the Load3, it not only operates steadily, but also works at a lower temperature.
This phenomenon is caused by the method difference of how to balance the optimal target of the two modes.Load3 is the largest load, and the centralized mode considers the optimization from the overall model, so the action of Load3 should be diminished and the lower loads work more to counteract disturbance.However, in the distributed mode, the controller manages only one load, and does not need to consider other problems.Figure 15 shows the state switching and total power of distributed control.From Figure 15b, the three loads' total consumption is able to be kept near 356 kW.But compared with Figure 9a, a sharp power change occurs in the distributed mode, and the curve is more unstable.This is caused by inaccurate coupling information between So, the optimal performance of the distributed mode is worse.The detailed comparison of total consumption and calculation performance between the two modes is listed in Table 7.Though there is a loss of optimal performance, the speed of distributed control is promoted.The detailed comparison of total consumption and calculation performance between the two modes is listed in Table 7.Though there is a loss of optimal performance, the speed of distributed control is promoted.

Conclusions
The characteristics of discrete and logical of the control system determine that in order to realize integrated analyzing and control of Power CPS, the first thing is to build a Hybrid System model which is able to describe the continuous dynamic process and discrete process for the primary system of the power grid.
This paper applies the MLD method to build a Hybrid System based Power CPS control model for ADN flexible load control.Combining with the ADN feeder power control, a flexible load control strategy is designed for normal and underpowered conditions of ADN feeder.To solve the problems of centralized control mode when the feeder is in the underpowered condition, a model and method of distributed control for flexible load is proposed.
From the study case, it can be concluded that compared with the traditional load control, the proposed flexible load control contains two advantages: (1) With the help of the hybrid model the predictability of the control is better, by which the flexible loads' operation constraints are satisfied strictly when the traditional loads violate the limit; (2) The hybrid model has the ability to combine the discrete process into the optimal control of continuous variables, in that the ADN flexible load state switching can be optimized to meet continuous variables, so that ADN flexible load state switching can be optimized to meet consumption demand.In the case study simulation of this paper, the flexible load power consumption is lower, and the most load's energy conservation reaches 11% of the traditional load; when the coordinates with DGs in the ADN, flexible load power consumption can be sustained near a target value and the deviation of RMS is 30% which is lower than 46% of the traditional load control.
Simulations of centralized and distributed control models of flexible load show that both of them have their advantages and disadvantages.The centralized control has higher efficiency and the target deviation is 7% lower than the distributed control in the study case.However, the distributed control is suitable for some complex situations, and computes faster (0.794 s in study case) than the centralized mode.Besides, the distributed control makes full use of the ADN control system, and has some practical value.
In conclusion, the Hybrid System based model proposed in this paper combines continuous dynamic, discrete state and even complex logical description in a unified expression, therefore physical process, control process and information flow is able to be described by mathematics directly.According to the model, the control method optimized continuous variables and state switching together in a control step.The ADN flexible load control is an example to illustrate the "integrating physical process and the control process" of the "Control Model", and the characteristics of the information system will be considered in the future to show the "integration of the primary and information systems".

Figure 1 .
Figure 1.Relationship of technical features and model in Power Cyber Physical Systems (CPS).

Figure 1 .
Figure 1.Relationship of technical features and model in Power Cyber Physical Systems (CPS).

Figure 3 .
Figure 3.Control system of Active Distribution Network (ADN).

Figure 3 .
Figure 3.Control system of Active Distribution Network (ADN).
ctrlLoad−m indicates the maximum consumption of n flexible loads.
of n flexible loads.

StartFigure 4 .
Figure 4. Block diagram of control of the ADN flexible loads.Figure 4. Block diagram of control of the ADN flexible loads.

Figure 4 .
Figure 4. Block diagram of control of the ADN flexible loads.Figure 4. Block diagram of control of the ADN flexible loads.

Figure 5 .
Figure 5. Topological connection of ADN feeder containing flexible load.

Figure 5 .
Figure 5. Topological connection of ADN feeder containing flexible load.

Figure 6 .
Figure 6.The function indexes comparison of refrigerator loads under normal condition; (a) flexible group; (b) traditional group.

Figure 6 .
Figure 6.The function indexes comparison of refrigerator loads under normal condition; (a) flexible group; (b) traditional group.

Figure 7 .
Figure 7.The power consumption comparison of refrigerator loads under normal condition.

Figure 8 .
Figure 8.The exchange power of ADN feeder and control areas under normal condition; (a) feeder and external system; (b) feeder and areas.

Figure 7 .
Figure 7.The power consumption comparison of refrigerator loads under normal condition.

Figure 8 .
Figure 8.The exchange power of ADN feeder and control areas under normal condition; (a) feeder and external system; (b) feeder and areas.

Figure 8 .
Figure 8.The exchange power of ADN feeder and control areas under normal condition; (a) feeder and external system; (b) feeder and areas.

Figure 9 .
Figure 9.The total power consumption comparison of ADN flexible loads in the underpowered condition; (a) flexible group; (b) traditional group.

Figure 10 .
Figure 10.The comparison of the ADN Feeder Power in underpowered condition; (a) state switching; (b) power consumption.

Figure 9 .
Figure 9.The total power consumption comparison of ADN flexible loads in the underpowered condition; (a) flexible group; (b) traditional group.

Figure 9 .
Figure 9.The total power consumption comparison of ADN flexible loads in the underpowered condition; (a) flexible group; (b) traditional group.

Figure 10 .
Figure 10.The comparison of the ADN Feeder Power in underpowered condition; (a) state switching; (b) power consumption.

Figure 10 .
Figure 10.The comparison of the ADN Feeder Power in underpowered condition; (a) state switching; (b) power consumption.

Figure 11 .
Figure 11.Comparison of centralized and distributed ADN load control in underpowered condition; (a) centralized control; (b) distributed control.

Figure 12 .
Figure 12.The relationship of coupling information among controllers.

P( 4 )
, where the meaning of δm and Pm can be found in 3.3.Let yp|t1+T -1 = yf-p, in other words, the unknown coupling variable yp|t1+T-1 is set artificially based on Pm.Besides, in the first loop of control, for time t = [0, T − 1], set yp|t = yf-p.Keep the consistency of coupling variables

Figure 12 .
Figure 12.The relationship of coupling information among controllers.
y f −p = P m , where the meaning of δ m and P m can be found in 3.3.Let y p|t1+T-1 = y f-p , in other words, the unknown coupling variable y p|t1+T -1 is set artificially based on P m .Besides, in the first loop of control, for time t = [0, T − 1], set y p|t = y f-p .(4) Keep the consistency of coupling variables

Figure 13 .
Figure 13.Block diagram of distributed predictive control of ADN flexible load.

Figure 14 .
Figure 14.The function indexes comparison of centralized and distributed control; (a) centralized control; (b) distributed control.

Figure 15 .
Figure 15.The state switching and power consumption of distributed control; (a) state switching; (b) power consumption.
Table 1 lists existing problems in the ADN control process, solutions are also given in the table.
Table 1 lists existing problems in the ADN control process, solutions are also given in the table.

Table 1 .
The problems of Active Distribution Network and solution based on Cyber Physical Systems (CPS).

Table 1 .
The problems of Active Distribution Network and solution based on Cyber Physical Systems (CPS).

Table 2 .
ADN flexible load control strategy.

Table 3 .
Main electrical specifications of feeder system and operation target.

Table 3 .
Main electrical specifications of feeder system and operation target.

Table 4 .
Parameters of flexible loads in three refrigerators.

Table 4 .
Parameters of flexible loads in three refrigerators.

Table 5 .
The comparison of power consumption under normal condition.

Table 5 .
The comparison of power consumption under normal condition.

Table 5 .
The comparison of power consumption under normal condition.

Table 6 .
The comparison of power consumption under power shortage condition.

Table 6 .
The comparison of power consumption under power shortage condition.

Table 6 .
The comparison of power consumption under power shortage condition.

Table 7 .
The comparison of centralized and distributed control.

Table 7 .
The comparison of centralized and distributed control.