Next Article in Journal
Methane Hydrate Formation in Marine Sediment from South China Sea with Different Water Saturations
Next Article in Special Issue
Switched Control Strategies of Aggregated Commercial HVAC Systems for Demand Response in Smart Grids
Previous Article in Journal
Synchronization of Low-Frequency Oscillation in Power Systems
Previous Article in Special Issue
Performance Analysis of Data-Driven and Model-Based Control Strategies Applied to a Thermal Unit Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Frequency Control Ancillary Service Provided by Efficient Power Plants Integrated in Queuing-Controlled Domestic Water Heaters

1
Key Laboratory of Smart Grid of Ministry of Education, Tianjin University, Tianjin 300072, China
2
State Power Economic Research Institute, Changping District, Beijing 102209, China
3
China Electric Power Research Institute, Haidian District, Beijing 100192, China
*
Author to whom correspondence should be addressed.
Energies 2017, 10(4), 559; https://doi.org/10.3390/en10040559
Submission received: 20 December 2016 / Revised: 9 April 2017 / Accepted: 13 April 2017 / Published: 19 April 2017
(This article belongs to the Special Issue Energy Management Control)

Abstract

:
Frequency is an important parameter of a power system. It is of great significance to maintain its stability, especially in the current development scenario of large-scale interconnected power systems. Thermostatically controlled appliances (TCAs) are good controllable resources for demand response owing to their rapid response capabilities and relatively wide controllable ranges. In this study, domestic water heaters, which have wider deadbands compared with other typical TCAs, such as heat pumps, are used as frequency regulation resources. The main contribution of this paper is that it proposes a queuing-controlled strategy with lock-on and off constraints for controlling an efficient power plant consisting of water heaters (EPP-WH). The queuing-controlled strategy enables TCAs to provide frequency regulation ancillary service for the normal operation of the power system. The thermal dynamic process of the water heater and the formation of the EPP-WH are first discussed. Based on the developed model, a series of strategies are proposed, including load shedding calculation, top layer optimization, and improved temperature priority list (TPL) strategy with lock-on and off constraints. Finally, typical case studies are discussed to illustrate the frequency regulation effects and the effects of two characteristic parameters—users’ willingness and lock time limits. Reasonable targets are generated based on various consideration from top layer optimization module. The results indicate that using the model and proposed strategies, the EPP-WH has good frequency regulation performance.

1. Introduction

Frequency is a critical parameter that indicates the balance between power generation and load consumption in power systems. When there is a severe disturbance or a drastic change in loading conditions, the operating state of the system may change from the normal state to the emergency state. Traditionally, the frequency is primarily controlled by adjusting resources on the generation side, including extra capacity from large generators and interconnection [1]. If the system’s frequency is still out of the normal range, effective measures must be taken to maintain the frequency stability, such as load shedding and modifying demand response. In such conditions, load shedding is regarded as one effective emergency control method, which has been an established practice to bring systems back to their normal operating state. It is a long-term measure often used after the primary frequency regulation fails. The main objective of load shedding is to bring the system’s frequency back to an acceptable range. This type of scheme is called under frequency load shedding (UFLS), which is an effective control measure for preventing frequency collapse after a large disturbance occurs or for addressing power shortages. Many studies have been carried out in this field. There are mainly three categories of control methods in the traditional UFLS schemes: the maximum available capacity methods, the linear methods, and techniques to enable system frequency recovery. In addition, fuzzy and binary control methods have also been implemented in frequency regulation [2]. The study in [3] proposed a frequency-based control strategy, taking into account renewable energy resources and energy storage. The impacts of emulated inertia from wind power were studied in [4]. In addition, some dispatch strategies were proposed using the aggregated electric vehicles to participate in the ancillary frequency regulation [5,6,7,8,9]. With the development of advanced metering infrastructure (AMI) and automatic meter reading (AMR), the control center of the power system can acquire more detailed consumption information about the electrical loads. For example, a simple low cost smart socket was proposed to switch off household appliances such as fans in order to provide a primary frequency response [10].
With respect to load control, studies have been carried out on providing energy market and ancillary services, such as load shifting [11], primary frequency response [12], and voltage stability enhancement [13]. Frequency control ancillary service provided by a wind farm dual battery energy storage system was studied in [14]. A control strategy for V2G aggregators participating in frequency regulation was studied in [15]. In [16], an automatic generation control (AGC) algorithm with EVs was proposed and the dynamic power system characteristics due to a given load disturbance were analyzed. In addition, a comprehensive central DR algorithm for frequency regulation in a smart microgrid was presented in [17]. In [18], use of the electricity demand as frequency controlled reserve (DFR) has been investigated to provide frequency reserve to power systems. A dynamic demand-side control algorithm is presented to improve the frequency stability of the power system, and to reduce the demand for system reserve capacity in [19]. In [20], a multi-stage demand response technology is introduced to participate in the primary frequency regulation of the system through direct load control. An adaptive scheme is studied in [21] for coordination of demand response and system spinning reserves for quick frequency recovery. In [22], the frequency regulation of demand response and wind power on power system is discussed in single area power system.
Due to the presence of considerable thermal inertia in thermostatically controlled appliances (TCAs), their short-term operation changes will not have too much impact on users’ comfort. As a consequence, TCAs are good controllable load resources for frequency regulation in power systems. In [23], a direct load control algorithm was developed to aggregate the water heater load for the purpose of frequency regulation. In [24], a supplementary load frequency control algorithm considering a number of EVs and the heat pump water heaters as controllable loads was proposed. A dynamic frequency regulation algorithm which used residential TCAs to alleviate frequency deviations caused by high penetration of renewable energy sources in the power system was proposed in [25]. In [26], a novel decentralized demand control strategy for family-friendly controllable refrigerators considering customer comfort level is proposed to regulate the frequency of autonomous microgrid in coordination with the energy storage system (ESS). Combining smart demand response technology with under frequency load shedding measures is involved in the third line of defense specified in a technical guide for electric power system security and stability control, whereby coordinated control could be carried out by using quick response resources and the third line of defense, and a new under frequency load shedding scheme considering the response of smart household appliances was proposed in [27].
Traditionally, UFLS involves two main steps: relatively more load shedding at the beginning of power deficit and less load shedding when the frequency is close to the lower limit of the normal operating range. However, after the second step, the frequency may rise over the upper limit of the normal range. As a result, a long-term method aimed at maintaining the system’s frequency should be continuously implemented. Considering the time scale of control, TCAs are suitable resources for frequency regulation. Owing to the unique features of thermostatic inertia, TCAs have less influence on the users’ comfortable levels when they participate in frequency regulation.
In this paper, a type of adaptive UFLS scheme is proposed, which is established based on a frequency differential and system frequency response model. The required load shedding to be implemented was acquired according to the rate of change of frequency. Water heaters, as a typical example of TCAs, were selected as the load shedding resource. The proposed strategy takes full advantage of the thermostatic energy storage characteristic of TCAs, demonstrating the frequency regulation capability of water heaters.
This paper is organized as follows: models of water heaters, efficient power plants with and without considering lock-on and off constraints are presented in Section 2. Related control strategies are introduced in Section 3. Case studies and simulation results are discussed in Section 4. Conclusions and future work are summarized in Section 5.

2. Modeling Methodologies

A detailed understanding of the thermal dynamic process of water heaters is necessary to evaluate the impacts of different factors on water heaters’ potential controllability. This section will discuss the thermal model of water heaters developed based on the ETP approach.

2.1. Single Water Heater Model

A single water heater’s thermal dynamic process can be described by Equation (1). Assuming the water temperature in the tank is θ t e m p t 1 at time t 1 . After Δ t , the water temperature in the tank at time t 1 + Δ t is calculated by:
{ θ t e m p t 1 + Δ t = θ a t 1 + Q R ( θ a t 1 + u t 1 Q R θ t 1 ) e Δ t / R C u t 1 { 0 , 1 } θ t e m p t 1 + Δ t = [ θ t e m p t 1 M θ s e t t 1 m n t 1 + θ a t 1 m n t 1 ] / M θ t e m p t 1 + Δ t > θ s e t t 1 θ t e m p t 1 + Δ t = [ θ t e m p t 1 ( M m n t 1 ) + θ a t 1 m n t 1 ] / M θ t e m p t 1 + Δ t θ s e t t 1
In this paper, it is assumed that users can use hot water when water heaters are in both modes. When there is hot water demand from a household, hot water is drawn from the tank, and cold water is added to its bottom; its water temperature can be expressed by the last line of Equation (1). M is the tank volume, m n t 1 stands for water consumption at time t 1 , θ a t 1 is the ambient temperature. By fitting the observed performance data, including user’s water demand and water temperature, thermal parameters, such as R, C, and Q can be determined.
As shown in Figure 1, the temperature in the tank develops between θ h i g h and θ l o w . θ s e t represents the setpoints of the water heater. When the water heater is on, the temperature increases along with power consumption, while when the water heater is off, it decreases. The curve in the ellipse simulates the temperature trend when the user demands a large quantity of hot water. The shaded rectangles below represent the power consumption. There is a unique correspondence between the temperature rising and the power consumption. The relationship of θ h i g h , θ l o w , θ s e t , and δ can be defined as follows:
{ θ h i g h = θ s e t + δ 2 θ l o w = θ s e t δ 2 δ = θ h i g h θ l o w

2.2. Efficient Power Plant Consisting of Water Heaters

Based on the single water heater’s thermal dynamic process, a group of water heaters can be regarded as an efficient power plant (EPP). Its output power limits P m i n t + Δ t , P m a x t + Δ t at t + Δ t , can be described as follows:
{ P m i n t + Δ t P E P P t + Δ t P m a x t + Δ t P m a x t + Δ t = P E P P t + P w h i   ( θ i , t + Δ t ( θ l o w i , t + Δ t , θ h i g h i , t + Δ t ) , u i , t = 0 ) P m i n t + Δ t = P E P P t P w h i   ( θ i , t + Δ t ( θ l o w i , t + Δ t , θ h i g h i , t + Δ t ) , u i , t = 1 )
In Figure 2, the shaded triangle represents the energy reduced as a result of adopting control strategies during the period Δ t , the gray dotted lines represent the EPP’s output limit, and the dark points are some possible operating status.
There are four main parameters: P m a x t and P m i n t , which are logically equivalent to the users’ comfortable level constraints; uncontrolled EPP state P E P P _ u t ; and controlled EPP state P E P P t . From Figure 2, it can be observed that if EPP-WH has a large feasibility range, and the operating point is near the upper limit, it has better conditions to become a good candidate for load shedding. Using these basic parameters, other main technical parameters can be defined as:
{ Δ P E P P t + Δ t = P E P P t + Δ t P E P P _ u t + Δ t Δ P u p t + Δ t = P m a x t + Δ t P E P P t + Δ t Δ P d o w n t + Δ t = P m i n t + Δ t P E P P t + Δ t
The controlled capacity Δ P E P P t + Δ t represents the capacity of the energy storage or the energy releasing capacity, that is to say that water heaters can absorb or release energy when the EPP-WH is under control; Δ P u p and Δ P d o w n represent the feasible regions of the EPP’s state, which are the feasible boundaries of the energy storage or the energy releasing capacity at this moment. In particular, the shedding capacity Δ P d o w n represents the EPP’s biggest power saving ability.

2.3. Lock On and Off Constraints

To prevent water heaters from being controlled too often, lock-on and off functions were included in the TPL [28] strategy. The basic logic of this constraint is to set minimum time limits t l i m i t o n or t l i m i t o f f . When a water heater is on or off, it should maintain its status for at least t l i m i t o n or t l i m i t o f f . The water heaters in lock-on and off status should not be controlled. Using parameter l o n i , t and l o f f i , t as controlling index, using t r o n i , t + 1 and t r o f f i , t + 1 as lock-on and off timer to record the lock time:
{ t r o n i , t + 1 = t r o n i , t + 1 ( l o n i , t = 1 ) t r o n i , t + 1 = 0 ( l o n i , t = 0 ) t r o f f i , t + 1 = t r o f f i , t + 1 ( l o f f i , t = 1 ) t r o f f i , t + 1 = 0 ( l o f f i , t = 0 )
l o n i , t = { 1 t r o n i , t < t l i m i t o n , t r o n i , t 0 0 t r o n i , t t l i m i t o n
l o f f i , t = { 1 t r o f f i , t < t l i m i t o f f , t r o f f i , t 0 0 t r o f f i , t t l i m i t o f f
If a water heater is being locked, the timer will keep counting, and during this period, EPP-WH cannot be controlled. After the timer reaches its time limits, the EPP-WH can be available again as shown in Equations (5)–(7). When calculating the limits of EPP-WH, Equation (3) needs to be updated:
{ P m i n _ l t + Δ t P E P P _ l t + Δ t P m a x _ l t + Δ t P m a x _ l t + Δ t = P E P P _ l t + P w h , i   ( θ i , t + Δ t ( θ l o w i , t + Δ t , θ h i g h i , t + Δ t ) , u i , t = 0 , l o n i , t = 0 ) P m i n _ l t + Δ t = P E P P _ l t P w h , i   ( θ i , t + Δ t ( θ l o w i , t + Δ t , θ h i g h i , t + Δ t ) , u i , t = 1 , l o n i , t = 0 )
The controlled capacity P E P P _ l t + Δ t represents the capacity of the energy storage or the energy releasing capacity, that is to say that water heaters can absorb or release energy when the EPP-WH is under control; Δ P u p and Δ P d o w n represent the feasible regions of the EPP’s state, which are the feasible boundaries of the energy storage or the energy releasing capacity at this moment. In particular, the shedding capacity Δ P d o w n represents the EPP’s biggest power saving ability:
{ Δ P E P P _ l t + Δ t = P E P P _ l t + Δ t P E P P _ u t + Δ t Δ P u p _ l t + Δ t = P m a x _ l t + Δ t P E P P _ l t + Δ t Δ P d o w n _ l t + Δ t = P m i n _ l t + Δ t P E P P _ l t + Δ t
Figure 3 shows the operating points of the EPP-WH, where the gray points represent the state without lock-on and off constraints. The upper and lower limits of the EPP-WH have shrunk. Compared with the uncontrolled point P E P P _ u t , the total energy restored can be seen as Section A plus section B. Taking lock-on and off constraints into account, some water heaters cannot be controlled even if they are within reasonable temperature ranges. Therefore, the target generation of the EPP-WH is becoming “conservative”, leading to a higher operating point Δ P E P P _ l t , and the energy restored decreased to only the area in section A. The load capacities of the rising, decreasing, and controlled were all reduced. But with this constraint, the users’ comfortable levels were maintained, which is a major advantage of this approach.

2.4. Users’ Willingness Model

Users’ willingness γ can be defined as the degree of initiative for water heater users to participate into demand response. It is affected by many factors, like educational level, family income, family members’ age [29,30] and so on. Generally, family members having a college or above education level have stronger willingness to participate in demand response; families with higher income increase their initiative; if the average age of family members is between 35–45, it is more likely to have a strong willingness to participate in demand response. In this part, three factors—educational level, income influence and family members’ age—are taken into account. Figure 4 shows the relationship of related factors and uses’ willingness. As the diversity of educational level is not as abundant as that of family income and family members’ age, three levels are set up to describe the influences of educational level on users’ willingness, E L stands for the education level.
γ E L = { b 1 0 < E L a 1 b 2 a 1 < E L a 2 b 3 a 2 < E L a 3
The low and middle income group account for the majority of the total group. This group of people have relatively weaker willingness to participate in demand response activities compared with the high income group, so an exponential model is set up to describe the influence of educational level on users’ willingness. Here F I stands for the family income:
γ F I = { e F I 0 < F I a 4 e a 4 F I > a 4
where a 4 stands for the maximum input, if F I is more than a 4 , the influence stays the same.
The influence of family members’ average age on users’ willingness can be described as follows:
γ F A = 1 2 π e ( ( F A a 5 ) 2 2 σ 2 )
The influences of family members’ average age obey a normal distribution with expectation a 5 and standard deviation σ , which means children and the aged are not the main participants in demand response programs. People aged between a 6 and a 7 (usually 35–45) will have more interest.
Based on the analysis above, an integrated users’ willingness model can be set up:
γ w h = γ E L γ F I γ F A m a x ( γ E L γ F I γ F A )
where γ w h has taken educational level, family income and family members’ average age into account. γ w h changes from 0 to 1 after normalization.

3. Optimization and Strategies

This section will discuss the main optimization methods and strategies. Figure 5 illustrates the different periods of frequency response and the corresponding resources. When the system is running normally, the frequency of the system is kept within a reasonable range, which is usually a ±0.2 Hz deviation of the system standard frequency (50 or 60 Hz). When there is a system disturbance, the frequency of the system may immediately drop out of the normal range. The first part of the UFLS will begin when the system frequency is below f e p p ; this part can be regarded as “emergency adjustment”, which undertakes the most load shedding task, with the EPP-WH, EV, and storage (such as batteries) as the frequency response resources.
When the frequency changes back to f n o r _ l , which is the lower limit of the system frequency, regulation adjustment begins in order to maintain the system frequency. A long-term approach will be implemented, which includes load shedding and rising, different from the emergency adjustment. Figure 6 illustrates the control structure and strategies.
By analyzing the status of the generators, the frequency deviation of the system Δ f and the required load to be shed Δ P can be calculated. With the evaluation of Δ f , the frequency response resource will be decided. Using information on the generator’s operating point, the total target Δ P can be calculated through a series of analyses which will be elaborated in more detail in the following section. The EPP-WH, EV, and storage receive their targets based on their percent power consumption in the system. In this study, we focused on the performance of the EPP-WH, assuming EV and storage has good and stable performance. To ensure that appropriate targets are dispatched, several elements were taken into account, including the EPP-WH “users’ willingness” and operating points. During the period of emergency adjustment, each type of resources decreases its power consumption to help pull the system frequency back to the normal range. While in the regulation adjustment period, response targets generated from the system includes both load increasing and load shedding.

3.1. Load Shedding Calculation

When a disturbance occurs in a multi-machine power system, with the measured dynamic system response, changes in rotor speed of synchronous generators can be detected first. In the process of dynamic frequency change, different generators have different influences on the system frequency. Therefore, the definition of the center of inertia (COI) frequency is put forward. In this paper, the rate of change of the COI frequency is regarded as the rate of change of the system frequency. Based on this, the rotor speed and the motion equation of the COI are expressed as follows:
{ ω C O I = 1 T J i = 1 N T J i ω i T J d w C O I d t = Δ P
where ω C O I is the rotor speed of the center of inertia (rad/s), T J = i = 1 N T J i , T J i is the input mechanical torque of the ith generator (kW·s·rad−1), f COI is the frequency of the center of inertia (Hz) based on the relationship ω C O I = 2 π f C O I , in per unit value, ω C O I = f C O I :
{ f C O I = 1 T J i = 1 N T J i f i T J d f C O I d t = Δ P
Using Equation (15), the target load shedding Δ P is achieved. Using information on water heater users’ willingness and the operating point, f e p p can be obtained as:
{ f e p p = f n o r γ × f P γ i = k = 1 N i γ w h i , k × P w h i , k P E P P i γ = i = 1 m γ i × P E P P i P E P P _ t o t a l i P E P P _ t o t a l i = i = 1 m P E P P i
where γ w h i , k stands for the kth water heater user’s willingness in the i-th EPP-WH, ranging from 0 to 1, P w h i , k is the power consumption of the k-th water heater. γ i can be calculated through a method of weighting considering users’ actual power consumption, and γ can be generated by weighting all the γ i in system. f P represents the maximum frequency range of the regulation frequency adjustment, which can be estimated by Equation (16). Using information from the system, appropriate targets can be identified to the execution units under different Δ f conditions. Equation (15) describes the relationship between the rate of frequency and the loads that need to be decreased or increased. Before the frequency returns back to f n o r _ l , Δ P is a negative parameter; during this emergency adjustment period, the system continues load shedding. While in the regulation adjustment process, Δ P may be a positive or negative parameter depending on the system operating status, the system has demands of shedding and increasing loads.

3.2. Target Assignment Method

Assuming there is a load shedding demand at time t, the EPP-WHs at higher operating points will be suitable to undertake more task, as they have larger Δ P d o w n to decrease their power consumption. Similarly, if one EPP-WH is operating at a lower point, it can be a good resource to response load rising target. To evaluate the influences of EPP-WH’s operating status, parameter λ is defined to calculate the degree of suitability for an EPP-WH to participate into demand response. λ can be obtained through analysis of the operating status of EPP-WH. Considering the lock-on effects and diversities of EPP-WHs, each EPP-WH has different operating status and power limits which are all time-varying, the degree of suitability for an EPP-WH to participate into demand response at time t can be defined as follows:
λ k , t = { Δ P u p _ l t / ( P m a x _ l t P m i n _ l t ) load shedding Δ P d o w n _ l t / ( P m a x _ l t P m i n _ l t ) load rising
Conclusion can be drawn that the larger the value of λ k , t is, the more targets should be distributed to this EPP-WH, and the more ith EPP-WH is likely to be a good response resource. Based on the analysis above, the objective function and constraints can be written as follows:
o b j min { k = 1 m ( 1 λ k , t ) | Δ P * k , t | } s t k = 1 m Δ P * k , t = Δ P t o t a l t    { Δ P d o w n _ l t + Δ t < Δ P * k , t < 0 load shedding 0 < Δ P * k , t < Δ P u p _ l t + Δ t load rising
where Δ P * k , t stands for the target distributed to ith EPP-WH, Δ P t o t a l t is the total target which can be generated by Equation (15). The equation constraint means all the EPP-WHs undertake the total target Δ P t o t a l t . Besides, every EPP-WH should satisfy its upper and lower limits considering the lock-on and off effects. By using 1 λ k , t as weighting factor, the target generated from objective function has considered the influences of users’ willingness λ k , t and the lock-on and off constraints. Using PSO method to solve this optimal problem, generating Δ P * k , t for each EPP-WH.

3.3. EPP-WH Optimization Strategies

In this study, a new TPL method with lock-on and off constraints was used to control the EPP-WH.
First, at each time step, taking into account the status of the water heaters, the controlled group were screened out, the upper and lower limits were calculated using Equation (3), the EPP-WH was developed, in which water heaters in “lock” status have been eliminated. Second, the water heaters were divided into two parts based on their switch’s status, making up two controlled groups as shown in Figure 7 and Equations (19)–(20):
O o p e n t = ( O 1 t , O 2 t , O 3 t , , O n 1 t ) ,
C c l o s e d t = ( C 1 t , C 2 t , C 3 t , , C n 2 t ) .
In O o p e n t , water heaters were sorted according to their temperature in descending order. While in C c l o s e d t , water heaters were sorted according to their temperature in ascending order.
In Figure 7, the temperature increases from left to right, the dark arrows represent water heaters in the O o p e n t at the time t, the gray arrows represent those in the C c l o s e d t , and the red arrows represent water heaters eliminated from the EPP-WH. All the arrows move in the clockwise direction, the water heaters at the end of each group will switch their status automatically and step into the other group at the next time step.
Finally, a group of water heaters will be turned off or turned on, that is to say, the EPP-WH decreases or increases the power consumption by a certain amount based on the type of targets sent to it. When the EPP-WH has to decrease its power consumption, the exact number of the water heaters to be turned off can be calculated by:
{ i = 1 k o + 1 P w h , i P E P P _ l t Δ P k , t i = 1 k o P w h , i < P E P P _ l t Δ P k , t
Then the top k o water heaters in the O o p e n t will be closed, and the EPP-WH decreases its power consumption for Δ P k , t . Similarly, when the EPP-WH has to increase its power consumption, the exact number of the water heaters to be turned on can be calculated by:
{ i = 1 k c + 1 P w h , i P E P P _ l t Δ P k , t i = 1 k c P w h , i < P E P P _ l t Δ P k , t
Then the top k c water heaters in the C c l o s e d t will be opened, and the EPP-WH increases its power consumption for Δ P k , t .

4. Simulation Results and Analysis

This section discusses the effects of several key parameters on the control methods. All the cases are based on the IEEE New England 10 generators 39 bus system. Figure 8 shows the topology of the system, while the active power of the load nodes are listed in Table 1. The total number of water heaters was 100,000 with an average power 6 kW. Assuming a smart controlling unit has been set, the control period of the EPP-WH was 10 s. A detailed description of the system is presented in the Appendix A. 3-order model was used to describe the status of the generators, the detailed parameters are presented in Table A1, the detailed slack generator parameters are listed below in Table A2.
In this study, to achieve a drastic frequency change, the load level increased by 28% at the start of the simulation, and the system frequency quickly decreased to 47.6 Hz as shown in Figure 9.

4.1. Case 1-Principle of Traditional UFLS Method

The traditional method requires four or more rounds to undertake load shedding. During this period, different percentages of Δ P load shedding will be carried out. As the load shedding percent has been set previously, it cannot achieve accurate frequency adjustment. Besides, the multi-rounds schedule comes up with time delay, which should not be ignored. A typical load shedding schedule is presented in Table 2.

4.2. Case 2-Comparision of Different UFLS Methods

In this section, a comparison of different frequency adjustment strategies is presented to illustrate the advantages of the proposed method. The two strategies are under the same simulation conditions, three frequency curves are shown in Figure 10, which represent the traditional method, the proposed method with sufficient resource, and the deficient resource, respectively. It can be observed that the proposed method had a shorter reaction time and less frequency deviations from the f n o r . This is because using frequency response resources, such as the EPP-WH, EV, and storage, targets and commands are generated in real-time; therefore, the entire response process can be more generalized. As the EPP-WH has a quick response speed, the influences of time delay can be decreased. However, if the frequency adjustment resource is not sufficient, the response curve may remain in the lower level as shown in Figure 10.
The EPP-WH power with deficient and sufficient resources are presented in Figure 11 and Figure 12. It was observed that only when the EPP-WH received a reasonable response target (within its response range), can the water heaters exhibit stable performance. If the target remains lower than the lower limit of the EPP-WH, a high number of water heaters will be switched on to maintain the users’ comfortable levels leading to a relatively large response error as shown in Figure 11. In this case, the critical point has divided the response process into two periods. In the beginning of a period of time, the EPP-WH had enough demand response resource to follow the target. While after the critical point, demand response resource decreased significantly. Since the number of water heaters in C c l o s e d t decreased over time, for a single water heater, being closed for a long time means that the water temperature will be continuously decreasing, more water heaters will switch on, and go into the O o p e n t again to maintain the basic temperature demand. Therefore, a reasonable response target is of great importance for the comparison of Figure 11 and Figure 12.

4.3. Case 3-Performance of EPP-WH Based on Different Water Heater Users’ Willingness

Figure 13 shows the EPP-WH’s performance of frequency regulation based on different users’ willingness: 10% and 80% degrees of participation of the EPP-WH were simulated. That is, 10% and 80% of all the water heater users agreed to be controlled during the period of demand response.
To study the effects of different users’ willingness independently, supplemental resources were forbidden. Considering the scale of the EPP-WH, appropriate targets were set. There is a time difference △t between the two curves first reaching a frequency of 50 Hz. It was observed that the more the users participate in the demand response, the less the time required for the recovery of system’s frequency. The normal range of frequency is 49.8–50.2 Hz in this case. At the beginning, a disturbance causes the frequency of the system to decrease from 50 Hz to approximately 49 Hz. After a time of 30 s, UFLS measures were implemented to prevent more severe frequency collapse. Targets were generated, and the EPP-WH responded to these targets by shedding a certain amount of loads. After a period of approximately 1 min, the trajectory of the frequency returned back to the normal range. Despite this, the willingness of the water heater users still had great effects on the targets sent to the EPP-WH as shown in Figure 14. From Figure 14, it can be observed that when the water heater users have stronger willingness to participate in load shedding, the top layer optimization module will send more targets to the EPP-WH. In the case of 80% degree of participation, the target increased after 6 min because most water heaters have already been turned off. To maintain the users’ comfortable levels, targets sent to the EPP-WH were adjusted, which also can be speculated, as the operating point was close to the lower limit of the EPP-WH.
If the system does not have abundant supplemental resources, and the EPP-WH has a lower heater users’ willingness to participate in the demand response at the same time, the system will have little controllable resources, resulting in poor frequency performance as shown in Figure 15.
Therefore, a certain amount of supplemental resources are of great significance as the UFLS is not the only method for maintaining the system’s frequency. Therefore, alternative schemes should be anticipated.

4.4. Case 4-Influence of Lock-On and Off Constrains

On the one hand, the UFLS is a type of long-term method to maintain the normal status of power systems. On the other hand, water heaters should not be controlled too often when there are no urgent frequency regulation demand. The control period of the EPP-WH will be extended to 1 min after the normal period of system’s frequency. Figure 16 shows the EPP-WH’s operating in the upper and lower limits without lock-on and off constraints and the status taking constrains into consideration with a 5 min time limit.
The limits of the EPP-WH were significantly reduced due to the lesser availability of resources. Though the feasible region was smaller, this constraint was beneficial because the controlled time of each water heater per minute has decreased from 3.2 to 1.3, which improved the users’ comfortable levels. Therefore, appropriate lock-on and off constraints should be encouraged on the premise of acceptable response errors.
Figure 17 shows the frequency performance with and without the lock-on and off constraints. As presented above, with the support of supplemental resources, the system’s frequency has been maintained in a reasonable region.

4.5. Case 5-Temperature Changes of EPP-WH When Controlled

Figure 18 shows the behavior of a single water heater with 16 °C dead-band and 50 °C setpoint, according to Equation (2), the upper and lower temperature limits can be calculated. The circle with dashed line on the left highlights the obvious curve of the temperature trajectory indicating the moment when hot water was used, according to Equation (1). The gray points indicate the moment when the water heater received the command to be turned off, while the black points indicate the time measures were implemented to turn on the water heater.
Figure 19 shows the results of the temperature trajectories of the EPP-WH with an average of 16 °C dead-band and 50 °C setpoint. Compared with Figure 1, it can be observed that a given number of water heaters have been controlled with some concaves left in the operating regions. As there is a unique correspondence between rising temperature and power consumption, the EPP-WH, which is integrated by TCAs, can produce good responses based on the targets sent from the control center.

5. Conclusions and Future Work

In this study, an efficient power plant was developed based on the thermostatic dynamic model of water heaters. A number of strategies, including load shedding calculation, top layer optimization, and improved TPL strategy with lock-on and off constraints, were proposed to establish the EPP-WH as a good candidate for frequency regulation. To evaluate the performance of the proposed method, case studies with different conditions were simulated under load control, and essential conclusions were drawn. The results showed that using the model and the proposed strategies, the EPP-WH had a good and stable frequency regulation performance. As a load-side solution, more types of TCAs, such as air-conditioners and heat pumps were recommended as frequency regulation resources. The optimization of lock-on and off time limits and price constraints will also be investigated in the future.

Acknowledgments

Acknowledgment: This research was supported by the National High Technology Research and Development Program (“863” Program) of China (No. 2015AA050403), the National Natural Science Foundation of China (No. 51407125). This study was carried out in cooperation with the APPLIED ENERGY UNiLAB-DEM: Distributed Energy & Microgrid. UNiLAB is an international virtual lab of collective intelligence in Applied Energy and Projects of State Grid Corporation of China (“Research on Key Technology of Intelligent Distribution Network Planning with Large Scale Flexible Loads”).

Author Contributions

Yebai Qi and Dan Wang conceived and designed the study; Yebai Qi and Dan Wang performed the study; Hongjie Jia, Tianjiao Pu, Naishi Chen and Kaixin Liu reviewed and edited the manuscript; Yebai Qi wrote the paper. All authors read and approved the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

In this paper, we used 3-order generator model; the detailed parameters are presented in Table A1 and the descriptions are as follows:
The effects of the leakage reactance and the armature resistance have been included. The differential equations are as follows:
δ ˙ = Ω b ( ω 1 ) ω ˙ = ( p m p e D ( ω 1 ) ) / M
The electrical power p e is defined as:
p e = ( υ q + r a i q ) i q + ( υ d + r a i d ) i d
In this model, all the axis electromagnetic circuits are neglected, whereas a lead lag transfer function is used for the d-axis inductance. The state variable e q is described by the following differential equation:
e ˙ q = ( f s ( e q ) ( x d x d ) i d + υ f * ) / T d 0
The voltage and current link are described by:
0 = υ q + r a i q e q + x d i d 0 = υ d + r a i d x q i q
υ d d-axis machine voltage
υ q q-axis machine voltage
Ω b base frequency in rad/s
υ f * the field voltage
r a armature resistance
e q q-axis transient voltage
Dequivalent damping winding
Mmechanical starting time
x d d-axis transient reactance
x q q-axis transient reactance
x d d-axis synchronous reactance
x q q-axis synchronous reactance
i d d-axis current
i q q-axis current
T d 0 d-axis open circuit transient time constant
T q 0 q-axis open circuit transient time constant
x l leakage reactance
Table A1. Generator parameters.
Table A1. Generator parameters.
No. x d (pu) x q (pu) x d (pu) x q (pu) T d 0 (s) T q 0 (s) x l (pu)
10.0060.0080.020.01970.70.003
20.06970.170.2950.2826.561.50.035
30.05310.08760.24950.2375.71.50.0304
40.04360.0660.2620.2585.691.50.0295
50.0320.1660.670.625.40.440.054
60.050.08140.2540.2417.30.40.0224
70.0490.0860.2950.2925.661.50.0322
80.0570.09110.290.286.70.410.028
90.0570.05870.21060.2054.791.960.0298
100.0910.0080.10.06910.200.0125
Slack generators are modeled as υ θ buses, i.e., constant voltage magnitude and phase generators, as follows:
υ = υ G 0 θ = θ G 0
Each network must include at least one slack generator. The angle θ G 0 is assumed to be the reference angle of the system. The detailed slack generator parameters are listed below in Table A2.
Table A2. Slack Generator data.
Table A2. Slack Generator data.
DataDescriptionUnit
31Bus numberint
100.0Power ratingMVA
100.00Voltage ratingkV
0.98200Voltage magnitudep.u.
0.00000Reference Anglep.u.
9999.00000Maximum reactive powerp.u.
−999.90000Minimum reactive powerp.u.
1.1Maximum voltagep.u.
0.9Minimum voltagep.u.
5.72835Active power guessp.u.
1Loss participation coefficient-
1Reference bus-
1Connection status-

References

  1. Rebours, Y.G.; Kirschen, D.S.; Trotignon, M.; Rossignol, S. A survey of frequency and voltage control ancil-lary services—Part I: Technical features. IEEE Trans. Power Syst. 2007, 22, 350–357. [Google Scholar] [CrossRef]
  2. Manjarres, P.; Malik, O. Frequency regulation by fuzzy and binary control in a hybrid islanded microgrid. J. Mod. Power Syst. Clean Energy 2015, 3, 429–439. [Google Scholar] [CrossRef]
  3. Oureilidis, K.O.; Bakirtzis, E.A.; Demoulias, C.S. Frequency-based control of islanded microgrid with renewable energy sources and energy storage. J. Mod. Power Syst. Clean Energy 2016, 4, 1–9. [Google Scholar] [CrossRef]
  4. Gonzalez-Longatt, F.M. Impact of emulated inertia from wind power on under-frequency protection schemes of future power systems. J. Mod. Power Syst. Clean Energy 2016, 4, 211–218. [Google Scholar] [CrossRef]
  5. Wu, C.; Mohsenian-Rad, H.; Huang, J. Vehicle-to-aggregator interaction game. IEEE Trans. Smart Grid 2012, 3, 434–442. [Google Scholar] [CrossRef]
  6. Liu, H.; Hu, Z.; Song, Y.; Wang, J.; Xie, X. Vehicle-to-Grid control for supplementary frequency regulation considering charging demands. IEEE Trans. Power Syst. 2015, 30, 3110–3119. [Google Scholar] [CrossRef]
  7. Wang, R.; Li, Y.; Wang, P.; Niyato, D. Design of a V2G aggregator to optimize PHEV charging and frequency regulation control. In Proceedings of the IEEE International Conference on Smart Grid Communications, Vancouver, BC, Canada, 21–24 October 2013; pp. 127–132. [Google Scholar]
  8. Karfopoulos, E.L.; Panourgias, K.A.; Hatziargyriou, N.D. Distributed coordination of electric vehicles providing V2G regulation services. IEEE Trans Power Syst. 2016, 31, 2834–2846. [Google Scholar] [CrossRef]
  9. Shi, W.; Wong, V.W.S. Real-time vehicle-to-grid control algorithm under price uncertainty. In Proceedings of the IEEE International Conference on Smart Grid Communication, Brussels, Belgium, 17–20 October 2011; pp. 261–266. [Google Scholar]
  10. Li, D. Studies on Some Key Problems of Wide-Area Measurements Based Power System on-Line Security Analysis; Huazhong University of Science & Technology: Wuhan, China, 2006. [Google Scholar]
  11. Heffner, G.; Goldman, C.; Kirby, B.J.; Meyer, M.K. Loads Providing Ancillary Services: Review of International Experience; Ernest Orlando Lawrence Berkeley National Laboratory: Berkeley, CA, USA, 2007. [Google Scholar]
  12. Samarakoon, K.; Ekanayake, J. Demand side primary frequency response support through smart meter control. In Proceedings of the 44th International Universities Power Engineering Conference (UPEC), Glasgow, UK, 1–4 September 2009. [Google Scholar]
  13. Wang, D.; Parkinson, S.; Miao, W.; Jia, H.; Crawford, C.; Djilali, N. Online voltage security assessment considering comfort-constrained demand response control of distributed heat pump systems. Appl. Energy 2012, 96, 104–114. [Google Scholar] [CrossRef]
  14. Yao, D.; Choi, S.S.; Tseng, K.J.; Lie, T.T. Frequency control ancillary service provided by a wind farm: Dual-BESS scheme. J. Mod. Power Syst. Clean Energy 2014, 2, 93–103. [Google Scholar] [CrossRef]
  15. Meng, J.; Mu, Y.; Wu, J.; Jia, H.; Dai, Q.; Yu, X. Dynamic frequency response from electric vehicles in the Great Britain power system. J. Mod. Power Syst. Clean Energy 2015, 3, 203–211. [Google Scholar] [CrossRef]
  16. Almeida, P.M.R.; Lopes, J.A.P.; Soares, F.J.; Vasconcelos, M.H. Automatic generation control operation with electric vehicles. In Proceedings of the Bulk Power System Dynamics and Control, Rio de Janeiro, Brazil, 1–6 August 2010. [Google Scholar]
  17. Pourmousavi, S.A.; Nehrir, M.H. Real-time central demand response for primary frequency regulation in microgrids. IEEE Trans. Smart Grid 2012, 3, 1988–1996. [Google Scholar] [CrossRef]
  18. Xu, Z.; Ostergaard, J.; Togeby, M. Demand as frequency controlled reserve. IEEE Trans. Power Syst. 2011, 26, 1062–1071. [Google Scholar] [CrossRef]
  19. Short, J.A.; Infield, D.G.; Freris, L.L. Stabilization of grid frequency through dynamic demand control. IEEE Trans. Power Syst. 2007, 22, 1284–1293. [Google Scholar] [CrossRef]
  20. Molina-Garcia, A.; Bouffard, F.; Kirschen, D.S. Decentralized demand-side contribution to primary frequency control. IEEE Trans. Power Syst. 2011, 26, 411–419. [Google Scholar] [CrossRef]
  21. Chang-Chien, L.R.; An, L.N.; Lin, T.W.; Lee, W.J. Incorporating demand response with spinning reserve to realize an adaptive frequency restoration plan for system contingencies. IEEE Trans. Smart Grid 2012, 3, 1145–1153. [Google Scholar] [CrossRef]
  22. Molina-García, A.; Muñoz-Benavente, I.; Hansen, A.D.; Gómez-Lázaro, E. Demand-side contribution to primary frequency control with wind farm auxiliary control. IEEE Trans. Power Syst. 2014, 29, 2391–2399. [Google Scholar] [CrossRef]
  23. Kondoh, J.; Lu, N.; Hammerstrom, D.J. An evaluation of the water heater load potential for providing regulation service. IEEE Trans. Power Syst. 2011, 26, 1309–1316. [Google Scholar] [CrossRef]
  24. Masuta, T.; Yokoyama, A. Supplementary load frequency control by use of a number of both electric vehicles and heat pump water heaters. IEEE Trans. Smart Grid 2012, 3, 1253–1262. [Google Scholar] [CrossRef]
  25. Rui, X.; Liu, X.; Meng, J. Dynamic frequency regulation method based on thermostatically controlled appliances in the power system. Energy Procedia 2016, 88, 382–388. [Google Scholar] [CrossRef]
  26. Jia, H.J.; Yan, Q.I.; Mu, Y. Frequency response of autonomous microgrid based on family-friendly controllable loads. Sci. China Technol. Sci. 2013, 56, 693–702. [Google Scholar] [CrossRef]
  27. Tang, Y.; Deng, K.; Sun, H.; Yi, J.; He, Q. Research on coordination scheme for smart household appliances participating under frequency load shedding. Power Syst. Technol. 2013, 37, 2861–2867. [Google Scholar]
  28. Katipamula, S.; Lu, N. Evaluation of residential HVAC control strategies for demand response programs. ASHRAE Trans. 2006, 112, 535–546. [Google Scholar]
  29. Faruqui, A.; Sergici, S. BGE’s Smart Energy Pricing Pilot Summer 2008 Impact Evaluation; The Brattle Group Inc.: Cambridge, UK, 2009; pp. 1–25. [Google Scholar]
  30. Faruqui, A.; Sergici, S. Pricing Programs: Time-of-Use and Real Time; Taylor & Francis Group: London, UK, 2007; pp. 1175–1183. [Google Scholar]
Figure 1. Thermal dynamic process of single water heater.
Figure 1. Thermal dynamic process of single water heater.
Energies 10 00559 g001
Figure 2. Efficient power plant model.
Figure 2. Efficient power plant model.
Energies 10 00559 g002
Figure 3. Changes in limits taking lock-on and off constraints into account.
Figure 3. Changes in limits taking lock-on and off constraints into account.
Energies 10 00559 g003
Figure 4. The relationships of related factors and uses’ willingness. (a) Educational level influence; (b) Family income influence; and (c) Family members’ age influence.
Figure 4. The relationships of related factors and uses’ willingness. (a) Educational level influence; (b) Family income influence; and (c) Family members’ age influence.
Energies 10 00559 g004
Figure 5. Different periods of frequency response and corresponding resources.
Figure 5. Different periods of frequency response and corresponding resources.
Energies 10 00559 g005
Figure 6. Flowchart of control structure and strategies.
Figure 6. Flowchart of control structure and strategies.
Energies 10 00559 g006
Figure 7. Structure of TPL method.
Figure 7. Structure of TPL method.
Energies 10 00559 g007
Figure 8. The 10 generator 39 bus system.
Figure 8. The 10 generator 39 bus system.
Energies 10 00559 g008
Figure 9. System frequency without adjustment strategies.
Figure 9. System frequency without adjustment strategies.
Energies 10 00559 g009
Figure 10. System frequency with different strategies.
Figure 10. System frequency with different strategies.
Energies 10 00559 g010
Figure 11. EPP-WH power with deficient resource.
Figure 11. EPP-WH power with deficient resource.
Energies 10 00559 g011
Figure 12. EPP-WH power with sufficient resource.
Figure 12. EPP-WH power with sufficient resource.
Energies 10 00559 g012
Figure 13. EPP-WH’s performance of frequency regulation with different users’ willingness.
Figure 13. EPP-WH’s performance of frequency regulation with different users’ willingness.
Energies 10 00559 g013
Figure 14. EPP-WH’s targets and responses based on different users’ willingness.
Figure 14. EPP-WH’s targets and responses based on different users’ willingness.
Energies 10 00559 g014
Figure 15. EPP-WH’s performance of frequency regulation without abundant supplemental resources.
Figure 15. EPP-WH’s performance of frequency regulation without abundant supplemental resources.
Energies 10 00559 g015
Figure 16. EPP-WH’s operating limits without lock-on and off constraints.
Figure 16. EPP-WH’s operating limits without lock-on and off constraints.
Energies 10 00559 g016
Figure 17. EPP-WH’s operating status taking lock-on and off constraints into account.
Figure 17. EPP-WH’s operating status taking lock-on and off constraints into account.
Energies 10 00559 g017
Figure 18. Response of a single water heater.
Figure 18. Response of a single water heater.
Energies 10 00559 g018
Figure 19. Temperature trajectories of EPP-WH.
Figure 19. Temperature trajectories of EPP-WH.
Energies 10 00559 g019
Table 1. Active power of load nodes.
Table 1. Active power of load nodes.
Bus-IDActive Power (MW)Water Heater NumberTotal Water Heater Power Capacity (MW)
3322524031.4
4500818049.1
7233.8380022.8
8522849050.9
128.51400.84
15320524031.4
16329535032.1
18158257015.4
2068011,06066.4
21274445026.7
23247.5402024.1
24308.6502030.1
25224364021.8
26139226013.6
27281456027.4
28206335020.1
29283.5461027.7
319.21500.9
39110417,870107.2
Total6150.1100,000600
Table 2. Load shedding schedule.
Table 2. Load shedding schedule.
RoundTrigger Frequency (Hz)Time Delay (s)Load Shedding Percent (%)
149.250.230
2490.230
348.750.220
448.50.220

Share and Cite

MDPI and ACS Style

Qi, Y.; Wang, D.; Wang, X.; Jia, H.; Pu, T.; Chen, N.; Liu, K. Frequency Control Ancillary Service Provided by Efficient Power Plants Integrated in Queuing-Controlled Domestic Water Heaters. Energies 2017, 10, 559. https://doi.org/10.3390/en10040559

AMA Style

Qi Y, Wang D, Wang X, Jia H, Pu T, Chen N, Liu K. Frequency Control Ancillary Service Provided by Efficient Power Plants Integrated in Queuing-Controlled Domestic Water Heaters. Energies. 2017; 10(4):559. https://doi.org/10.3390/en10040559

Chicago/Turabian Style

Qi, Yebai, Dan Wang, Xuyang Wang, Hongjie Jia, Tianjiao Pu, Naishi Chen, and Kaixin Liu. 2017. "Frequency Control Ancillary Service Provided by Efficient Power Plants Integrated in Queuing-Controlled Domestic Water Heaters" Energies 10, no. 4: 559. https://doi.org/10.3390/en10040559

APA Style

Qi, Y., Wang, D., Wang, X., Jia, H., Pu, T., Chen, N., & Liu, K. (2017). Frequency Control Ancillary Service Provided by Efficient Power Plants Integrated in Queuing-Controlled Domestic Water Heaters. Energies, 10(4), 559. https://doi.org/10.3390/en10040559

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop