METHOD FOR REGULATING PRIMARY FREQUENCY OF POWER GRID BASED ON AIR CONDITIONING LOAD CLUSTER IN LARGE BUILDING

Abstract

The present disclosure relates to a method for primary frequency regulation of an electric network based on large building air conditioning loads cluster. The method includes the use of a two layer control structure with a central coordinating layer and a local control layer. Each local controller performs a thermal model parameter identification and a local air conditioning autonomous control, and uploads local information to the central controller at the end of each communication interval t.sub.gap, the central controller broadcasts coordinating information to each local controller. Based on the coordinating information sent from the central controller, each local controller determines whether a power deviation is beyond an action dead zone at the beginning of each action period t.sub.act, if beyond, then perform a frequency regulation control action, else, perform no action and estimate operation states of all the air conditionings at the beginning of the next action period.

Claims

1. A method for regulating primary frequency of a power grid based on an air conditioning load cluster in a large building, wherein, a two-layer control structure comprising a central coordinating layer and a local control layer is used in the air conditioning load cluster, the central coordinating layer comprises a central controller, the local control layer comprises N local controllers, N air conditionings, and temperature sensors and frequency sensors provided in rooms the air conditionings located in; and the method comprises: 1) performing, by each local controller, a thermal model parameter identification and an air conditioning autonomous control to obtain local information corresponding to each of the air conditionings, and uploading the local information to the central controller at an end of each communication interval t.sub.gap, and broadcasting, by the central controller, coordinating information to each local controller; 2) when a communication between the central controller and each of the local controllers in step 1) is finished, based on the coordinating information sent from the central controller, determining, by each local controller, whether a power deviation in the air conditionings is beyond an action dead zone at a beginning of each action period t.sub.act, if yes, a frequency regulation control action is performed, else, no action is performed and operation states of all the air conditionings at a beginning of a next action period are estimated; if a current time reaches to a beginning of a next communication interval, step 1) is executed, else, step 2) is repeated.

2. The method according to claim 1, wherein, step 1) comprises: 1-1) performing, by each local controller i, i=1.2 . . . N, the room thermal model parameter identification according to air temperature data recorded at each temperature acquisition period to obtain thermal model parameters corresponding to each room; 1-2) identifying, by the local controller i, parameters of the thermal model corresponding to room i according to air temperature data recorded at each temperature acquisition period t.sub.temp in a communication interval t.sub.gap to obtain identified thermal model parameters of each room; 1-3) performing, by each local controller, the air conditioning autonomous control according to following equations: $\begin{array}{cc}{\mathrm{state}}_i=\{\begin{array}{cc}1,& {\mathrm{Ta}}_i\ge {\stackrel{\_}{T}}_i\&.Math.{\mathrm{state}}_i=0\\ 0,& {\mathrm{Ta}}_i\le {\underset{\_}{T}}_i\&.Math.{\mathrm{state}}_i=1\end{array}& \left(6\right)\\ {\stackrel{\_}{T}}_i={\mathrm{Ts}}_i+{\Delta }_i,{\underset{\_}{T}}_i={\mathrm{Ts}}_i-{\Delta }_i& \left(7\right)\end{array}$ where, i=1.2 . . . N, Ta.sub.i is an air temperature in the i th room, Δ.sub.i is a temperature control dead zone corresponding to an i.sup.th air conditioning, T.sub.i corresponds to an upper bound of the temperature control dead zone Δ.sub.i, T.sub.i corresponds to a lower bound of the temperature control dead zone Δ.sub.i, Ts.sub.i is the required temperature corresponding to the i th air conditioning set by the user, state.sub.i is the on-off state of the i th conditioning, wherein state.sub.i=1 corresponds to state ON, state.sub.i=0 corresponds to state OFF; 1-4) at an end moment of the communication interval t.sub.gap between the local controller and the central controller, uploading, by each local controller, the local information to the central controller, wherein the local information comprises the indoor air temperature acquired most recently Ta.sub.i of room i, the on-off state state.sub.i, the operation power P.sub.i, the required temperature Ts.sub.i, the temperature control dead zone Δ.sub.i, and the thermal model parameters α.sub.i.sup.ON,γ.sub.i.sup.ON, α.sub.i.sup.OFF,γ.sub.i.sup.OFF, Ta.sub.i.sup.tog, and t.sub.i.sup.tog; 1-5) collecting, by the central controller, all the local information from the local controllers and broadcasting all collected information to each local controller as the coordinating information, obtaining, by the central controller, a reference power P0.sub.i of each air conditioning after the thermal model parameters corresponding to each local controller are collected, obtaining a reference power P0 of all the air conditionings by summing reference powers of all the air conditionings, and broadcasting the reference power of all the air conditionings to each of the local controllers, wherein the coordinating information comprises the indoor air temperature Ta.sub.i, the on-off state state.sub.i, the operation power P.sub.i, the required temperature Ts.sub.i, the temperature control dead zone Δ.sub.i, and the thermal model parameters α.sub.i.sup.ON,γ.sub.i.sup.ON, α.sub.i.sup.OFF,γ.sub.i.sup.OFF, Ta.sub.i.sup.tog, and t.sub.i.sup.tog wherein the coordination information comprises indoor air temperatures, on-off states, operation powers, required temperatures, temperature control dead zones, and thermal model parameters uploaded by all the local controller.

3. The method according to claim 2, wherein a precision degree of the thermal model parameters is determined according to a hardware storage capability of the local controller and an error requirement between a thermal model identification curve and an actual temperature curve.

4. The method according to claim 3, wherein the thermal model corresponding to i.sup.th room comprises a zero-order thermal model, a first-order thermal model, or a second-order thermal model, represented by equations (1)-(3) respectively:
ΔT.sub.i=α.sub.iΔt.sub.i  (1)
ΔT.sub.i=α.sub.ie.sup.γ.sup.i.sup.Δt.sup.i−α.sub.i  (2)
ΔT.sub.i=α.sub.i1e.sup.γ.sup.i1.sup.Δt.sup.i+α.sub.i2e.sup.γ.sup.i2.sup.Δt.sup.i−α.sub.i1−α.sub.i2  (3) where, numbers of parameters to be identified in the three thermal models are 1, 2, and 4 respectively, α.sub.i in equation (1) is a thermal model parameter to be identified in the zero-order thermal model, α.sub.i,γ.sub.i in equation (2) are thermal model parameters to be identified in the first-order thermal model, α.sub.i1,γ.sub.i1,α.sub.i2,γ.sub.i2 in equation (3) are thermal model parameters to be identified in the second-order thermal model, ΔT.sub.i is a difference between a current temperature Ta.sub.i and an indoor temperature Ta.sub.i.sup.tog when an on-off state of the air conditioning is last switched, and Δt.sub.i is a difference between a current time and a time t.sub.i.sup.tog when an on-off state of the air conditioning is last switched, where,
ΔT.sub.i=Ta.sub.i−Ta.sub.i.sup.tog  (4)
Δt.sub.i=t.sub.i−t.sub.i.sup.tog  (5) wherein, Ta.sub.i.sup.tog and t.sub.i.sup.tog are thermal model parameters.

5. The method according to claim 4, if the thermal model corresponding to i.sup.th room is the first-order thermal model, the thermal model parameters comprises α.sub.i.sup.ON,γ.sub.i.sup.ON and α.sub.i.sup.OFF,γ.sub.i.sup.OFF.

6. The method according to claim 2, wherein, the reference power P0.sub.i corresponds to an average power of the i th air conditioning during an on-off period T.sub.i in a communication interval t.sub.gap, and step 1-5) comprises: obtaining, by the central controller, a first time t.sub.i.sup.(1), a second time t.sub.i.sup.(1), a third time t.sub.i.sup.(1), and a forth time t.sub.i.sup.(1) by solving the following equations respectively according to the upper T.sub.i, the lower temperature bound T.sub.i, the thermal model parameters α.sub.i.sup.ON,γ.sub.i.sup.ON, α.sub.i.sup.OFF,γ.sub.i.sup.OFF, Ta.sub.i.sup.tog, and t.sub.i.sup.tog, the required temperature Ts.sub.i, and the temperature control dead zone Δ.sub.i:
T.sub.i−Ta.sub.i.sup.tog=α.sub.i.sup.ONe.sup.γ.sup.i.sup.ON.sup.(t.sup.i.sup.(1).sup.-t.sup.i.sup.tog.sub.)−α.sub.i.sup.ON
T.sub.i−Ta.sub.i.sup.tog=α.sub.i.sup.ONe.sup.γ.sup.i.sup.ON.sup.(t.sup.i.sup.(2).sup.-t.sup.i.sup.tog.sub.)−α.sub.i.sup.ON
T.sub.i−Ta.sub.i.sup.tog=α.sub.i.sup.OFFe.sup.γ.sup.i.sup.OFF.sup.(t.sup.i.sup.(3).sup.-t.sup.i.sup.tog.sub.)−α.sub.i.sup.OFF
T.sub.i−Ta.sub.i.sup.tog=α.sub.i.sup.OFFe.sup.γ.sup.i.sup.OFF.sup.(t.sup.i.sup.(4).sup.-t.sup.i.sup.tog.sub.)−α.sub.i.sup.OFF wherein, t.sub.i.sup.(1) is a moment when the indoor temperature is equal to the upper bound temperature T.sub.i and the air conditioning is in an “ON” state; t.sub.i.sup.(2) is a moment when the indoor temperature is equal to the lower bound temperature T.sub.i and the air conditioning is in an “ON” state; t.sub.i.sup.(3) is a moment when the indoor temperature is equal to the upper bound temperature T.sub.i and the air conditioning is in an “OFF” state; t.sub.i.sup.(4) is a moment when the indoor temperature is equal to the lower bound temperature T.sub.i and the air conditioning is in an “OFF” state; obtaining a total time period Ton.sub.i when the i.sup.th air conditioning is in an “ON” state in an on-off period T.sub.i, and a total time period Toff.sub.i when the i.sup.th air conditioning is in an “OFF” state in an on-off period T.sub.i according to following equations:
Ton.sub.i=t.sub.i.sup.(1)−t.sub.i.sup.(2)  (8)
Toff.sub.i=t.sub.i.sup.(4)−t.sub.i.sup.(3)  (9) obtaining the reference power P0.sub.i of each air conditioning as: $\begin{array}{cc}P.Math..Math.0_i=\frac{{\mathrm{Ton}}_i}{T_i}.Math.P_i=\frac{{\mathrm{Ton}}_i}{{\mathrm{Ton}}_i+{\mathrm{Toff}}_i}.Math.P_i& \left(10\right)\end{array}$ wherein, P0.sub.i is an reference power of the i.sup.th air conditioning, Ton.sub.i is the total time period when the i.sup.th air conditioning is in an “ON” state in an on-off period T.sub.i, Toff.sub.i is the total time period when the i.sup.th air conditioning is in an “OFF” state in an on-off period T.sub.i, P.sub.i is an operation power of the i.sup.th air conditioning; obtaining the reference power of all the air conditionings P0 by summing all the reference powers P0.sub.i of the air conditionings according to following equation: $\begin{array}{cc}P.Math..Math.0=\underset{i}{.Math.}.Math..Math.P.Math..Math.0_i;& \left(11\right)\end{array}$ broadcasting, by the central controller, the reference power P0 of all the air conditionings to each local controller.

7. The method according to claim 1, wherein, step 2) comprises: 2-1) acquiring, by a frequency sensor, a frequency of the power grid every action period t.sub.act, and calculating, by each local controller, a power deviation δ of all the air conditionings according to the acquired frequency of the power grid at the beginning of each action period t.sub.act; 2-2) determining, by each local controller, whether the power deviation δ is in the action dead zone ξ, when the power deviation δ is in the action dead zone ξ, the air conditioning does not participate in the frequency regulation control; when the power deviation δ is not in the action dead zone ξ, the air conditioning participates in the frequency regulation control action in the present action period; 2-3) estimating, by each local controller, on-off states of all the air conditionings at a beginning of a next action period; 2-4) estimating, by each local controller, air temperatures in other rooms at the beginning of the next action period, and modifying the on-off state state.sub.i of the i th air conditioning at the beginning of the next action period t.sub.act according to the coordinating parameters transmitted from the central controller and estimated on-off states of all the air conditionings via the air conditioning autonomous control in step 1-2), and executing step 2-1) when the next action period comes or executing step 1) when a next communication interval begins.

8. The method according to claim 7, wherein, step 2-1) comprises: calculating, by each local controller, a real-time total power P(t) of all the air conditionings according to the received coordinating information broadcasted by the central controller via following equation: $\begin{array}{cc}P\left(t\right)=\underset{i}{.Math.}.Math..Math.P_i*{\mathrm{state}}_i& \left(13\right)\end{array}$ where. i=1.2 . . . N, P.sub.i is an operation power of the i.sup.th air conditioning, state.sub.i is an on-off state of the i.sup.th air conditioning; calculating the power deviation δ of all the air conditionings according to following equation:
δ=P(t)−P0−KΔf, where, P(t) is the real-time total power of all the air conditionings, P0 is the reference power of all the air conditionings, K is a power-frequency response coefficient set for all the local controllers, Δf is a real-time frequency deviation;

9. The method according to claim 7, wherein, step 2-2) comprises: 2-2-1) obtaining a temperature priority Tpri.sub.i of each local controller according to following equation: $\begin{array}{cc}{\mathrm{Tpri}}_i=\{\begin{array}{cc}\left({\mathrm{Ts}}_i-{\mathrm{Ta}}_i\right)/{\Delta }_i,& {\mathrm{state}}_i=1\\ \left({\mathrm{Ta}}_i-{\mathrm{Ts}}_i\right)/{\Delta }_i,& {\mathrm{state}}_i=0\end{array}& \left(14\right)\end{array}$ where, Tpri.sub.i is a temperature priority of i.sup.th local controller, Ta.sub.i is the indoor air temperature, Ts.sub.i is the required temperature corresponding to the i.sup.th air conditioning set by the user, Δ.sub.i is the temperature control dead zone, state.sub.i is the on-off state of the i th air conditioning; 2-2-2) when δ>ξ, selecting temperature priorities of air conditionings whose state.sub.i=1, and generating an array qu.sub.ON accordingly with its rows arranged according to values of the temperature priorities in a descending order, wherein, a first column of the array is the temperature priorities, a second column is operation powers corresponding to the temperature priorities, a third column is mark numbers of air conditionings corresponding to the temperature priorities, and a number of rows in the array qu.sub.ON is denoted as r; selecting a minimum regulation control set which can regulate the power deviation into the dead zone according to r*=min{r|Σ.sub.d=1.sup.rqu.sub.ON(d,2)≧δ−ξ}, extracting a set of mark numbers of air conditionings to be regulated in a present operation from the minimum regulation control set according to I.sub.ON=qu.sub.ON(j,3), j=1, 2, . . . , r*, and calculating I.sub.ON′={iεI.sub.ON|Ta.sub.i<Tgon.sub.i}; if a number of an air conditioning controller i.sub.localεI.sub.ON′, controlling an air conditioning corresponding to the an air conditioning controller i.sub.local to participate in the present frequency regulation control, i.e. switching a state of the air conditioning corresponding to the an air conditioning controller i.sub.local, else, performing no action: 2-2-3) when δ<−ξ, selecting temperature priorities of air conditionings whose state.sub.i=0, and generating an array qu.sub.OFF accordingly with its rows arranged according to values of the temperature priorities in a descending order, wherein, a first column of the array qu.sub.OFF is the temperature priorities, a second column is operation powers corresponding to the temperature priorities, a third column is mark numbers of air conditionings corresponding to the temperature priorities, and a number of rows in the array qu.sub.OFF is denoted as r; selecting a minimum regulation control set which can regulate the power deviation into the dead zone according to r*=min{r|Σ.sub.d=1.sup.rqu.sub.OFF(d,2)≧−δ−ξ}, extracting a set of mark numbers of air conditionings to be controlled in a present operation from the minimum regulation control set according to I.sub.OFF=qu.sub.OFF(j,3), j=1, 2, . . . , r*, and calculating I.sub.OFF′={iεI.sub.OFF|Ta.sub.i>Tgoff.sub.i}, if a number of an air conditioning controller i.sub.localεI.sub.OFF′, controlling the air conditioning corresponding to the an air conditioning controller i.sub.local to participate in the present frequency regulation control, i.e. switching a state of the air conditioning corresponding to the an air conditioning controller i.sub.local, else, performing no action.

10. The method according to claim 1, wherein the central controller and the local controllers communicate in both-way at every communication interval, the local controllers acquire data from the temperature sensors at each temperature sampling period.

11. The method according to claim 1, wherein the local controllers regulate and control the air conditionings once during each action period according to local information and coordination information transmitted from the central controller.

12. The method according to claim 1, wherein each of the air conditionings is constant power air conditioning with an operation power and having two states; ON and OFF.

13. The method according to claim 1, wherein each temperature sensor acquires an indoor air temperature of a corresponding room in real-time, and the local controller acquires temperature data from the temperature sensors every the temperature acquisition period.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0011] FIG. 1 is a schematic diagram showing a two-layer control structure used in embodiments of the present disclosure.

[0012] FIG. 2 is a flow chart showing a method regulating primary frequency of a power grid based on an air conditioning load cluster in a large building according to embodiments of the present disclosure; and

[0013] FIG. 3 is a flow chart showing a method regulating primary frequency of a power grid based on an air conditioning load cluster in a large building according to embodiments of the present disclosure.

DETAILED DESCRIPTION

[0014] A method for regulating primary frequency of a power grid based on an air conditioning load cluster in a large building provided in the present disclosure will be described in combination with embodiments and with reference to drawings as follows.

[0015] The present embodiment is shown as FIG. 1, a two-layer, i.e. a central coordinating layer and a local control layer, control structure is used in an air conditioning loads cluster. The central coordinating layer includes a central controller. The local control layer includes N local controllers. N air conditionings and temperature sensors and frequency sensors provided in rooms the air conditionings located in.

[0016] In some embodiments, the central controller and the local controllers may communicate in both-way at every communication interval. The local controllers acquire data from the temperature sensors at each temperature sampling period. A communication among the central controller and each of the local controllers is in a way of wireless communication. A communication between each of the local controllers and each of the air conditionings, each of the temperature sensors or the frequency sensors may be in a way of wireless communication or wire communication. The local controllers regulate and control the local air conditionings once during each action period according to local information and coordinating information transmitted from the central controller.

[0017] As shown in FIG. 2 and FIG. 3, a method for regulating primary frequency of a power grid based on an air conditioning load cluster in a large building provided in the present disclosure includes following acts.

[0018] In block 1, each local controller performs a thermal model parameter identification and an air conditioning autonomous control to obtain local information corresponding to each of the air conditionings, and uploads the local information to the central controller at an end of each communication interval t.sub.gap, the central controller broadcasts coordinating information to each local controller.

[0019] In an embodiment of the present disclosure, the communication interval t.sub.gap may be a period between 15 seconds to 1 minute.

[0020] In block 2, when a communication between the central controller and each of the local controllers in block 1 is finished, based on the coordinating information sent from the central controller, each local controller begins to determine whether a power deviation in the air conditionings is beyond an action dead zone at a beginning of each action period t.sub.act, if yes, a frequency regulation control action is performed, else, no action is performed and operation states of all the air conditionings at a beginning of a next action period are estimated.

[0021] If a current time reaches to a beginning of a next communication interval, block 1 is executed, else, block 2 is repeated.

[0022] In an embodiment of the present disclosure, the action period t.sub.act may be 1 second or other preset times.

[0023] In an embodiment of the present disclosure, the next action period is an action period t.sub.act interval from current moment, i.e. there is one action period t.sub.act between the beginning of the next action period and the current moment.

[0024] In some embodiments, block 1 includes following sub-acts.

[0025] In block 11, each local controller i, i=1.2 . . . N, performs the room thermal model parameter identification according to air temperature data recorded at each temperature acquisition period to obtain thermal model parameters corresponding to each room. N is a number of the local controllers.

[0026] In some embodiments, a precision degree of the thermal model parameters is determined according to a hardware storage capability of the local controller and an error requirement between a thermal model identification curve and an actual temperature curve.

[0027] In some embodiments, for each air conditioning room i, i=1.2 . . . N, three precision degrees of thermal model may be determined. The three precision degrees of thermal model include a zero-order thermal model, a first-order thermal model, and a second-order thermal model represented by equations (1)-(3) respectively.

ΔT.sub.i=α.sub.iΔt.sub.i  (1)

ΔT.sub.i=α.sub.ie.sup.γ.sup.i.sup.Δt.sup.i−α.sub.i  (2)

ΔT.sub.i=α.sub.i1e.sup.γ.sup.i1.sup.Δt.sup.i+α.sub.i2e.sup.γ.sup.i2.sup.Δt.sup.i−α.sub.i1−α.sub.i2  (3)

[0028] where, numbers of parameters to be identified in the three precision degrees of thermal models are 1, 2, and 4 respectively, i.e. α.sub.i in equation (1) is a thermal model parameter to be identified in the zero-order thermal model, α.sub.i,γ.sub.i in equation (2) are thermal model parameters to be identified in the first-order thermal model, α.sub.i1,γ.sub.i1,α.sub.i2,γ.sub.i2 in equation (3) are thermal model parameters to be identified in the second-order thermal model, ΔT.sub.i (an initial value of switching temperature for short) is a difference between a current temperature Ta.sub.i and an indoor temperature Ta.sub.i.sup.tog when an on-off state of the air conditioning is last switched, and Δt.sub.i (an initial value of switching time for short) is a difference between a current time and a time t.sub.i.sup.tog when an on-off state of the air conditioning is last switched, and,

ΔT.sub.i=Ta.sub.i−Ta.sub.i.sup.tog  (4)

Δt.sub.i=t.sub.i−t.sub.−.sup.tog  (5)

[0029] The initial value of switching temperature Ta.sub.i.sup.tog and the initial value of switching time t.sub.i.sup.tog are taken as parameters of the thermal model as well as α.sub.i in equation (1), α.sub.i,γ.sub.i in equation (2), or α.sub.i1,γ.sub.i1,α.sub.i2,γ.sub.i2 in equation (3).

[0030] In block 12, the local controller i identifies parameters of the thermal model corresponding to room i according to air temperature data recorded at each temperature acquisition period t.sub.temp (a period between 1 to 4 seconds) in a communication interval t.sub.gap to obtain thermal model parameters corresponding to each room.

[0031] In some embodiments, k is denoted as a number of times that the temperature is recorded, and the on-off state state.sub.i (where state.sub.i=1 corresponding to state ON, state.sub.i=0 corresponding to state OFF) of the air conditioning i is recorded at each time a temperature is recorded. When a zero-order (linear) thermal model is selected, a corresponding parameter identification model is

[00001]$\min .Math.\left\{\underset{k}{.Math.}.Math..Math.{\left[\left({\mathrm{Ta}}_i^k-{\mathrm{Ta}}_i^{\mathrm{tog}}\right)-\alpha \left(t_i^k-t_i^{\mathrm{tog}}\right)\right]}^2\right\},$

when an on-off state of the air conditioning i is state.sub.i=1, state parameter α.sub.i.sup.ON corresponding to the ON state is identified according to currently recorded k sets of switching temperature and switching time, when the on-off state of the air conditioning i is state.sub.i=0, state parameter α.sub.i.sup.OFF corresponding to the OFF state is identified according to currently recorded k sets of switching temperature and switching time. When a first-order model is selected, a corresponding parameter identification model is

[00002]$\min .Math.\left\{\underset{k}{.Math.}.Math..Math.{\left[\left({\mathrm{Ta}}_i^k-{\mathrm{Ta}}_i^{\mathrm{tog}}\right)-\left(\alpha .Math..Math.e^{\gamma \left(t_i^k-t_i^{\mathrm{tog}}\right)}-\alpha \right)\right]}^2\right\},$

similarly, two sets of parameters α.sub.i.sup.ON,γ.sub.i.sup.ON and α.sub.i.sup.OFF,γ.sub.i.sup.OFF are identified in different states ON and OFF respectively. When a two-order model is selected, a corresponding parameter identification model is

[00003]$\min .Math.\left\{\underset{k}{.Math.}.Math..Math.{\left[\left({\mathrm{Ta}}_i^k-{\mathrm{Ta}}_i^{\mathrm{tog}}\right)-\left({\alpha }_1.Math..Math.e^{{\gamma }_1\left(t_i^k-t_i^{\mathrm{tog}}\right)}+{\alpha }_2.Math.e^{{\gamma }_2\left(t_i^k-t_i^{\mathrm{tog}}\right)}-{\alpha }_1-{\alpha }_2\right)\right]}^2\right\},$

similarly, two sets of parameters α.sub.i1.sup.ON,γ.sub.i1.sup.ON,α.sub.i2.sup.ON,γ.sub.i2.sup.ON and α.sub.i1.sup.OFF,γ.sub.i1.sup.OFF,α.sub.i2.sup.OFF,γ.sub.i2.sup.OFF are identified in different states ON and OFF respectively.

[0032] In some embodiments of the present disclosure, a same thermal model is selected to use for all the local controllers. In an embodiment, a first-order model is selected to use for all the local controllers, the identified parameters of the i th room are α.sub.i.sup.ON,γ.sub.i.sup.ON and α.sub.i.sup.OFF,γ.sub.i.sup.OFF via common algorithms.

[0033] In block 13, each local controller performs the air conditioning autonomous control, according to following equations.

[00004]$\begin{array}{cc}{\mathrm{state}}_i=\{\begin{array}{cc}1,& {\mathrm{Ta}}_i\ge {\stackrel{\_}{T}}_i\&.Math.{\mathrm{state}}_i=0\\ 0,& {\mathrm{Ta}}_i\le {\underset{\_}{T}}_i\&.Math.{\mathrm{state}}_i=1\end{array}& \left(6\right)\\ {\stackrel{\_}{T}}_i={\mathrm{Ts}}_i+{\Delta }_i,{\underset{\_}{T}}_i={\mathrm{Ts}}_i-{\Delta }_i& \left(7\right)\end{array}$

[0034] In the above equations, i=1.2 . . . N, Ta.sub.i is an air temperature in the i th room, Δ.sub.i is a temperature control dead zone corresponding to the i th air conditioning, T.sub.i corresponds to an upper bound of the temperature control dead zone Δ.sub.i, T.sub.i corresponds to a lower bound of the temperature control dead zone Δ.sub.i, Ts.sub.i is the required temperature set by the user, state.sub.i is the on-off state of the i th conditioning (state.sub.i=1 corresponds to state ON, state.sub.i=0 correspond to state OFF).

[0035] Denoting an i th local controller, an i th air conditioning and an i th room with mark i, i=1.2 . . . N. An air temperature in the i th room is Ta.sub.i. An on-off state of the i th air conditioning is state.sub.i (where state.sub.i=1 corresponding to state ON, state.sub.i=0 corresponding to state OFF). In some embodiments, it is assumed that the i th air conditioning is a constant power air conditioning with an operation power P.sub.i and with on-off state controlled only.

[0036] Each temperature sensor acquires an indoor air temperature of a corresponding room in real-time. Each local controller acquires the temperature data from a corresponding temperature sensor every temperature acquisition period t.sub.temp (for example, a period between 1 to 4 seconds).

[0037] A required temperature Ts.sub.i of each air conditioning i is set directly by the user, and each air conditioning i has a temperature control dead zone Δ.sub.i, which is a factory setting attribute, and in an embodiment, Δ.sub.i is assumed to be 1° C. Equation (6) shows that if state.sub.i=0, i.e. the air conditioning is in an OFF state, when the room air temperature Ta.sub.i rises to the upper bound T.sub.i of the temperature control dead zone of the i th local air conditioning, the i th local controller turns on the i th air conditioning autonomously; if state.sub.i=1, i.e. the air conditioning is in an ON state, when the room air temperature Ta.sub.i drops to the lower bound T.sub.i of the temperature control dead zone of the i th local air conditioning, the i th local controller turns off the i th air conditioning autonomously.

[0038] In block 14, at an end moment (i.e. a communication moment) of the communication interval t.sub.gap between the local controller and the central controller, each local controller uploads the local information to the central controller.

[0039] The local information includes the indoor air temperature acquired most recently Ta.sub.i of the room, the on-off state state.sub.i of the air conditioning, the operation power P.sub.i, the required temperature Ts.sub.i, the temperature control dead zone Δ.sub.i, and the thermal model parameters α.sub.i.sup.ON,γ.sub.i.sup.ON, α.sub.i.sup.OFF,γ.sub.i.sup.OFF, Ta.sub.itog, and t.sub.i.sup.tog.

[0040] In block 15, the central controller collects all the local information from the local controllers and broadcasts all collected information to each local controller as the coordinating information, and the central controller obtains a reference power P0.sub.i of each air conditioning after the thermal model parameters corresponding to each local controller are collected, a sum of reference powers of all the air conditionings is obtained as a reference power P0 of all the air conditionings, and the reference power of all the air conditionings P0 is broadcasted to each local controller.

[0041] The coordinating information includes the indoor air temperature Ta.sub.i, the on-off state state.sub.i, the operation power P.sub.i, the required temperature Ts.sub.i, the temperature control dead zone Δ.sub.i, and the thermal model parameters α.sub.i.sup.ON,γ.sub.i.sup.ON, α.sub.i.sup.OFF,γ.sub.i.sup.OFF, Ta.sub.i.sup.tog, and t.sub.i.sup.tog.

[0042] In some embodiments, the reference power P0.sub.i corresponds to an average power of the i th air conditioning during an on-off period T.sub.i (referring to a time period during which the i th air conditioning switches its on-off state in one cycle according to a local autonomous control logic) in a communication interval t.sub.gap.

[0043] Block 15 may include following acts.

[0044] The central controller calculates a first time t.sub.i.sup.(1), a second time t.sub.1.sup.(1), a third time t.sub.i.sup.(1), and a forth time t.sub.i.sup.(1) by solving the following equations respectively according to the upper T.sub.i, the lower temperature bound T.sub.i, the thermal model parameters α.sub.i.sup.ON,γ.sub.i.sup.ON, α.sub.i.sup.OFF,γ.sub.i.sup.OFF, Ta.sub.i.sup.tog, and t.sub.i.sup.tog, the required temperature Ts.sub.i, and the temperature control dead zone Δ.sub.i.

T.sub.i−Ta.sub.i.sup.tog=α.sub.i.sup.ONe.sup.γ.sup.i.sup.ON.sup.(t.sup.i.sup.(1).sup.-t.sup.i.sup.tog.sub.)−α.sub.i.sup.ON

T.sub.i−Ta.sub.i.sup.tog=α.sub.i.sup.ONe.sup.γ.sup.i.sup.ON.sup.(t.sup.i.sup.(2).sup.-t.sup.i.sup.tog.sub.)−α.sub.i.sup.ON

T.sub.i−Ta.sub.i.sup.tog=α.sub.i.sup.OFFe.sup.γ.sup.i.sup.OFF.sup.(t.sup.i.sup.(3).sup.-t.sup.i.sup.tog.sub.)−α.sub.i.sup.OFF

T.sub.i−Ta.sub.i.sup.tog=α.sub.i.sup.OFFe.sup.γ.sup.i.sup.OFF.sup.(t.sup.i.sup.(4).sup.-t.sup.i.sup.tog.sub.)−α.sub.i.sup.OFF

[0045] In above equations, t.sub.i(1) is a moment when the indoor temperature is equal to the upper bound temperature T.sub.i and the air conditioning is in an “ON” state; t.sub.i.sup.(2) is a moment when the indoor temperature is equal to the lower bound temperature T.sub.i and the air conditioning is in an “ON” state; t.sub.i.sup.(3) is a moment when the indoor temperature is equal to the upper bound temperature T.sub.i and the air conditioning is in an “OFF” state; t.sub.i.sup.(4) is a moment when the indoor temperature is equal to the lower bound temperature T.sub.i and the air conditioning is in an “OFF” state.

[0046] A total time period Ton.sub.i when the i.sup.th air conditioning is in an “ON” state in an on-off period T.sub.i, and a total time period Toff.sub.i when the i.sup.th air conditioning is in an “OFF” state in an on-off period T.sub.i are obtained according to following equations.

Ton.sub.i=t.sub.i.sup.(1)−t.sub.i.sup.(2)  (8)

Toff.sub.i=t.sub.i.sup.(4)−t.sub.i.sup.(3)  (9)

[0047] The reference power P0.sub.i of each air conditioning is calculated as:

[00005]$\begin{array}{cc}P.Math..Math.0_i=\frac{{\mathrm{Ton}}_i}{T_i}.Math.P_i=\frac{{\mathrm{Ton}}_i}{{\mathrm{Ton}}_i+{\mathrm{Toff}}_i}.Math.P_i& \left(10\right)\end{array}$

[0048] In which, P0.sub.i is an reference power of the i.sup.th air conditioning, Ton.sub.i is the total time period when the i.sup.th air conditioning is in an “ON” state in an on-off period T.sub.i, Toff.sub.i is the total time period when the i.sup.th air conditioning is in an “OFF” state in an on-off period T.sub.i, P.sub.i is an operation power of the i.sup.th air conditioning.

[0049] The reference power of all the air conditionings P0 is obtained by summing all the reference powers P0.sub.i of the air conditionings according to following equation.

[00006]$\begin{array}{cc}P.Math..Math.0=\underset{i}{.Math.}.Math..Math.P.Math..Math.0_i& \left(11\right)\end{array}$

[0050] The central controller broadcasts the reference power P0 of all the air conditionings to each local controller.

[0051] In some embodiments, a control objective of primary frequency regulation response of the air conditioning cluster is set as making a difference ΔP between a real-time total power P(t) of all the air conditionings and the reference power P0 of all the air conditionings to be directly proportional to a real-time frequency deviation Δf, satisfying following equation.

ΔP=P(t)−P0=K(f(t)−f.sub.0)=KΔf  (12)

[0052] In which, f.sub.0 is a reference frequency, being 50 Hz for Chinese mainland, f(t) is a real-time frequency obtained by the frequency sensor, K is a power-frequency response coefficient and set to a same value for all the local controllers. K may be determined according to a ratio of a total power of the air conditioning cluster to a maximum frequency fluctuation in history. The greater K is, the more the air conditioning involves in the frequency regulation, and the smaller K is, the less the air conditioning involves in the frequency regulation.

[0053] Block 2) includes the following actions.

[0054] In block 21, a frequency of the power grid is acquired by a frequency sensor every action period t.sub.act, and each local controller calculates a power deviation δ of all the air conditionings according to the acquired frequency of the power grid at the beginning of each action period t.sub.act.

[0055] In some embodiments, each local controller calculates the real-time total power P(t) of all the air conditionings according to the received coordinating information (air temperature Ta.sub.i, on-off state state.sub.i and power Pi of the air conditionings in all the rooms) broadcasted by the central controller via following equation:

[00007]$\begin{array}{cc}P\left(t\right)=\underset{i}{.Math.}.Math..Math.P_i*{\mathrm{state}}_i& \left(13\right)\end{array}$

[0056] where. i=1.2 . . . N, P.sub.i is the operation power of the i.sub.th air conditioning, state.sub.i is an on-off state of the i.sup.th air conditioning.

[0057] Then, the power deviation δ of all the air conditionings is calculated according to following equation.

δ=P(t)−P0−KΔf,

[0058] where, P(t) is the real-time total power of all the air conditionings, P0 is the reference power of all the air conditionings, K is a power-frequency response coefficient set for all the local controllers, Δf is a real-time frequency deviation.

[0059] In block 22, each local controller determines whether the power deviation δ is in the action dead zone ξ, when the power deviation δ is in the action dead zone ξ, the air conditioning does not participate in the frequency regulation control, when the power deviation δ is not in the action dead zone ξ, the air conditioning participates in the frequency regulation control action in the present action period.

[0060] In some embodiments, ξ may be set according to accuracy requirement, for example, in an embodiment, ξ is 1 KW.

[0061] In some embodiments, if|ξ|≧ξ, the power deviation δ is determined to be in the action dead zone ξ.

[0062] In some embodiments, block 22 includes following actions.

[0063] In block 221, a temperature priority Tpri.sub.i of each local controller is obtained according to following equation.

[00008]$\begin{array}{cc}{\mathrm{Tpri}}_i=\{\begin{array}{cc}\left({\mathrm{Ts}}_i-{\mathrm{Ta}}_i\right)/{\Delta }_i,& {\mathrm{state}}_i=1\\ \left({\mathrm{Ta}}_i-{\mathrm{Ts}}_i\right)/{\Delta }_i,& {\mathrm{state}}_i=0\end{array}& \left(14\right)\end{array}$

[0064] In which, Tpri.sub.i is a temperature priority of i.sup.th local controller, Ta.sub.i is the indoor air temperature, Zs.sub.i the required temperature corresponding to the i.sup.th air conditioning set by a user, Δ.sub.i s the temperature control dead zone, state.sub.i is the on-off state of the i th air conditioning (ON corresponds 1, OFF corresponds to 0).

[0065] Equation (14) means that the air conditioning in a room where the air temperature is closer to a boundary of the temperature control dead zone corresponds a higher priority. When state.sub.i=1, i.e. the air conditioning is in an “ON” state, the lower the air temperature Ta.sub.i is, the higher the priority Tpri.sub.i is, and the air conditioning will be turned off more preferentially in a local frequency regulation process; When state.sub.i=0, i.e. the air conditioning is in an “OFF” state, the higher the air temperature Ta.sub.i is, the higher the priority Tpri.sub.i is, and the air conditioning will be turned on more preferentially in the local frequency regulation process.

[0066] In block 222, when δ>ξ, temperature priorities Tpri.sub.i of air conditionings whose state.sub.i=1 are selected, and an array qu.sub.ON is generated accordingly with its rows arranged according to values of the temperature priorities Tpri.sub.i in a descending order, the first column of the array is Tpri.sub.i, the second column is P.sub.i, the third column is i, and the number of rows in the array qu.sub.ON is denoted as r, a minimum control set which can regulate the power deviation into the dead zone is selected according to r*=min{r|Σ.sub.d=1.sup.rqu.sub.ON(d,2)≧δ−ξ}, a set I.sub.ON of numbers of the air conditionings to be regulated in the present operation is extracted from the minimum regulation control set according to I.sub.ON=qu.sub.ON(j,3), j=1, 2, . . . , r*, and I.sub.ON′={iεI.sub.ON|Ta.sub.i<Tgon.sub.i} is calculated (in which, as a parameter represents a participating degree of the air conditioning in the frequency regulation, Tgon.sub.i may be preset by users of the air conditionings and people who controls the frequency regulation system, for example, Tgon.sub.i may be set as T.sub.i+0.8Δ.sub.i), if a number of an air conditioning controller i.sub.localεI.sub.ON′, an air conditioning corresponding to the an air conditioning controller i.sub.local is controlled to participate in the present frequency regulation control, i.e. a state of the air conditioning corresponding to the an air conditioning controller i.sub.local is switched (turn off the local air conditioning), else, no action is performed.

[0067] In block 223, when δ<−ξ, temperature priorities Tpri.sub.i of air conditionings whose state.sub.i0 are selected, and an array qu.sub.OFF is generated accordingly with its rows arranged according to values of the temperature priorities Tpri.sub.i in a descending order, the first column of the array is Tpri.sub.i, the second column is P.sub.i, the third column is i, and the number of rows in the array qu.sub.OFF is denoted as r, a minimum control set which can regulate the power deviation into the dead zone is selected according to r*=min{r|Σ.sub.d=1.sup.rqu.sub.OFF(d,2)≧−δ−ξ}, a set I.sub.OFF of numbers of the air conditionings to be regulated in the present operation is extracted from the minimum regulation control set according to qu.sub.OFF(j,3), j=1, 2, . . . , r*, and I.sub.OFF′={iεI.sub.OFF|Ta.sub.i>Tgoff.sub.i} is calculated (in which, as a parameter represents a participating degree of the air conditioning in the frequency regulation, Tgoff.sub.i may be preset by users of the air conditionings and people who controls the frequency regulation system, for example, Tgoff.sub.i may be set as T.sub.i−0.4Δ.sub.i), if a number an air conditioning controller i.sub.localεI.sub.OFF′, the air conditioning corresponding to the an air conditioning controller i.sub.local is controlled to participate in the present frequency regulation control, i.e. a state of the local air conditioning corresponding to the an air conditioning controller i.sub.local is switched (turn on the local air conditioning), else, no action is performed.

[0068] In block 23, after the frequency regulation in block 22 is finished, each local controller estimates the on-off states of all the air conditionings at a beginning of a next action period.

[0069] In some embodiments, the next action period is one action period t.sub.act ahead from the present moment.

[0070] Because on-off states of some of the air conditionings have been changed in the present action period, each local controller estimates the on-off states of all the air conditionings at the beginning of the next action period after each frequency regulation in block 22 is finished.

[0071] The process (block 23) of estimating the on-off states of all the air conditionings includes following actions.

[0072] In block 231, the set I.sub.ON′ or I.sub.OFF′ in block 22) are obtained.

[0073] In block 232, it is set that i=1.

[0074] In block 233, it is determined whether iεI.sub.ON to determine whether the i th air conditioning participates in the frequency regulation action, if iεI.sub.ON, the present state of the i th air conditioning is “OFF”, i.e. state.sub.i0, and the air temperature Ta.sub.i.sup.tog before the switch is flipped and the moment t.sub.i.sup.tog when the switch is flipped are recorded. if iεI.sub.OFF, the present state of the i th air conditioning is “ON”, i.e. state.sub.i=1, and the air temperature Ta.sub.i.sup.tog before the switch is shifted and the moment t.sub.i.sup.tog when the switch is shifted are recorded;

[0075] In block 234, let i=i+1, if i≦N block 233 is executed, else, block 24 is executed.

[0076] In block 24, each local controller estimates air temperatures in other rooms at the beginning of the next action period, and modifies on-off state state.sub.i of the i th air conditioning at the beginning of the next action period t.sub.act according to the coordinating parameters transmitted from the central controller and the on-off states of all the air conditionings estimated in block 23 via the autonomous control method.

[0077] Block 21 is executed when the next action period comes, or block 1 is executed when the next communication interval begins.

[0078] The autonomous control method can refer to equations (6) and (7) illustrated in block 13.

[0079] Block 24 may include following acts.

[0080] In block 241, a first-order thermal model is used for estimating the temperature in a present embodiment, for the i th air conditioning, a time variance relative to t.sub.i.sup.tog at moment t is Δt.sub.i, and i=1.

[0081] In block 242, if the on-off state of the i th air conditioning stored locally is “ON”, i.e. state.sub.i=1, the room air temperature stored locally is Ta.sub.i(t)=α.sub.i.sup.ONe.sup.γ.sup.i.sup.ON.sup.Δt.sup.i−α.sub.i.sup.ON+Ta.sub.i.sup.tog; if the on-off state of the i th air conditioning stored locally is “OFF”, i.e. state.sub.i=0, the room air temperature stored locally is Ta.sub.i(t)=α.sub.i.sup.OFFe.sup.γ.sup.i.sup.OFF.sup.Δt.sup.i−α.sub.i.sup.OFF+Ta.sub.i.sup.tog.

[0082] In block 243, if the room air temperature meets a condition Ta.sub.i(t)≦T.sub.i, the state of the i th air conditioning stored locally is state.sub.i=0, and Ta.sub.i.sup.tog and t.sub.i.sup.tog are recorded, if the room air temperature meets a condition Ta.sub.i(t)≧T.sub.i, the state of the i th air conditioning stored locally is state.sub.i=1, and Ta.sub.i.sup.tog and t.sub.i.sup.tog are recorded, else, the state of the i th air conditioning state.sub.i remains unchanged.

[0083] In block 244, let i=i+1, if i≦N block 242 is executed, else, block 21 is executed when the next action period t.sub.act comes, or, block 1 is executed if next communication moment comes.

[0084] The method for regulating primary frequency of a power grid based on an air conditioning load cluster in a large building according to embodiments of the present disclosure has following characteristics.

[0085] With the method according to embodiments of the present disclosure, by taking advantage of heat capacity of large buildings, the two layers control structure including the central coordinating layer and the local control layer is formed, via rapid control of the air conditioning cluster, it is now possible to involve the air conditioning cluster in primary frequency regulation with a linear power-frequency characteristic similar to that of an electric generator. Each air conditioning performs a primary frequency regulation response locally to increase the speed of entire response, and communicate with the central controller every a certain time interval to upload local information and obtain overall coordinating information so as to ensure accuracy of the entire power-frequency linear response. Meanwhile, a temperature monitor threshold guarantees comfort of users and useful life of the equipment will not be significantly affected.

[0086] With the method according to embodiments of the present disclosure, a two-layer control structure including a slow centralized coordination and a rapid distributed local control is provided. In the centralized coordinating layer (i.e. the central coordinating layer), information of each room, such as operation power and state of corresponding air conditioning, the temperature in the room, thermal model parameters of the room, etc., is collected according to coordinating control period, and then is broadcasted to each of the local controllers. The problems of long time delay in centralized control and lack of coordination in distributed control are solved, thus improving the control accuracy and solving the problem of rebound effect.

[0087] The controllers in the local control layer estimate operation states of all the air conditionings and temperatures in all rooms based on dynamic thermal models respectively, and sequence the air conditioning cluster accordingly. Whether a local frequency regulation will be triggered is determined by an order of the air conditioning in the sequence and a real-time frequency deviation. The method is a local algorithm, thus improving speed of response to the power-frequency deviation.

[0088] With the method according to embodiments of the present disclosure, loads of the air conditioning cluster is a linear response to a frequency deviation and the comfort of users is not effected, the contradiction of slow response in central control and lack of coordinating information in local decentralized control when an air conditioning cluster participating in primary frequency regulation is eliminated. The control accuracy is increased and the problem of rebound effect is solved.