Impacts of Responsive Loads and Energy Storage System on Frequency Response of a Multi-Machine Power System

In recent decades, the power grid’s configuration is shifting towards a smart grid where responsive loads and energy storage systems (ESS) are finding an increased role in the power system operation. In the presented work, a mathematical formulation for frequency response analysis of a multi-machine power system is developed, considering the individual and combined roles of ESS and responsive loads. The validity of the developed model is demonstrated with the help of multiple case studies, which consider the various configurations of the power system. Moreover, different combinations of capacities of generation units, ESS and responsive loads are also simulated. With the help of mathematical model and simulation results, it is demonstrated that ESS and responsive loads may improve the economy and performance of the power system even during the failure of a certain portion of generation capacity. Though the case studies consider non-reheat turbines only, the mathematical model and conclusions are equally valid for other types of turbines as well.


Introduction
Amongst the primary requirements for a safe and reliable operation of a power system is a stable frequency response which depends on the balance between real power's supply and demand. Though power supply and demand continually change in a power system, the power system frequency response is maintained within a specific operating range through a process known as frequency regulation. Traditionally, the frequency regulation is realized using spinning and non-spinning reserves during normal conditions; and, under-frequency or over-frequency load-shedding methods are adopted during emergency situations [1][2][3]. However, other methods of frequency regulation are also gaining the popularity due to environmental and financial concerns. Utilization of renewable energy resources, energy storage systems (ESS) and responsive loads are some of the famous methods for the frequency regulation [4][5][6].
In a smart grid environment, the responsive loads are integrated through demand response (DR) programs which are advantageous to consumers and utility [7,8]. Several strategies of integration of responsive loads and their benefits in enhancing the power grid's performance and economy are reported in literature. The authors of Reference [9] developed a 2-layer strategy for microgrid frequency regulation with the help of responsive loads considering electric vehicles as a part of them and showed a better performance of the microgrid with the developed strategy. The authors of References [10,11] proposed a control mechanism based on manipulation of responsive loads for primary frequency regulation. In Reference [10], a modified load frequency control (LFC) model considering responsive loads for frequency regulation process is developed and verified by simulation results. This work presented in Reference [11] splits the microgrid operation into 3 modes depending upon the magnitude of frequency deviation. For each operating mode, different share of responsive loads takes part to minimize the frequency deviation and the overall scheme shows a promising role of DR in frequency regulation process. Similarly, the simulation results of work carried out in Reference [12] show the effectiveness of responsive loads for frequency regulation.
The use of ESS to enhance grid stability is also under increased consideration in recent decades. ESS makes an optimal fit to accommodate intermittent renewable energy resources in smart grids and can provide virtual inertia (VI) to overcome the low-inertia problem associated with renewable energy systems [13]. The literature shows the application of VI for frequency regulation services and it is expected that role of ESS will increase in future for power system's stability [14,15]. A control scheme is proposed in Reference [16] to suppress the deviation in system's voltage and frequency using VI. The study in Reference [17] proposed a virtual inertia control-based model predictive control for frequency regulation. In this work, the test system is exposed to a high share of intermittent renewable energy systems and the proposed scheme exhibits a robust performance. In References [18][19][20], detailed discussions on the mathematical modeling, control system development and experimental results involving hardware in the loop are presented.
Considering an escalated and simultaneous integration of responsive loads and ESS in smart grids, further research is needed to discuss their combined role. The authors of this work presented a modified LFC model considering responsive loads and ESS and discussed their individual and combined role in improving the frequency response [21]. In Reference [22], a droop controlled distributed generation framework integrated with DR is developed for frequency regulation with the objectives of emission reduction and cost minimization. The work presented in References [23,24] considers renewable energy resources and demand response in the context of a virtual power plant.
In comparison to our previous work and other literature, the contribution and highlights of this work are presented as follows: 1.

2.
A mathematical model is developed to include the role of responsive loads and ESS in multi-machine electric power system. For their integration, independent control loops for the DR and the VI powers are introduced, which provide a means of extra degrees of freedom for frequency regulation. Moreover, this approach provides a simpler structure to design the controllers.

3.
The developed model is verified with the help of simulation results carried out under a comprehensive set of case studies. The case studies range from consideration of an individual machine to combined role of responsive loads and ESS in a multi-machine system. The results show that the responsive loads and ESS can play a vital role in frequency regulation in case of partial or temporary unavailability of generation resources.
The paper is organized as follows: In Section 2, a model for a multi-machine power system considering responsive loads and ESS is developed. Case studies are discussed in Sections 3 and 4 summarizes and concludes the paper.  Figure 1 shows a general block diagram of the LFC model of a multi-machine power system. The system is equipped with primary or droop control (with droop control value represented by R in Figure 1) and secondary control. Typically, the secondary control is a part of automatic generation control and provides an integral control action to diminish the steady-state frequency deviation. For the objective of frequency response analysis and control under a load disturbance (∆P d (s)), Equation (1) models the frequency deviation as follows [1,2].

Development of the Model
where, D and H are the damping coefficient and inertia constant of the power system, respectively, ∆P i T (s) is the incremental power from i-th turbine out of a total N turbines. The system shown in Figure 1 considers general forms of turbines and governors. The transfer functions for governor and three commonly available turbine types (i.e. reheated, non-reheated and hydro turbines) are given in Equations (2) and (3), respectively [26].
where, D and H are the damping coefficient and inertia constant of the power system, respectively, ( ) i T P s Δ is the incremental power from i -th turbine out of a total N turbines. The system shown in Figure 1 considers general forms of turbines and governors. The transfer functions for governor and three commonly available turbine types (i.e. reheated, non-reheated and hydro turbines) are given in Equations (2) and (3), respectively [26]. where,

Responsive Loads
In the smart grids, responsive loads make an important part of the power system for long-term energy balance and to provide ancillary services. Integration of responsive loads is realized with the help of DR programs. There are two major categories of DR programs:

1.
Dispatchable programs: where the system operator can manipulate the responsive loads according to power grid's needs. For example, direct load control is a type of such programs and also considered in this work.

2.
Non-dispatchable programs: the programs in which the customer is financially encouraged to manipulate the loads but the system operator does not have direct control over the responsive loads. Its examples of are real-time pricing and peak time rebate [27].
The concept of manipulation of responsive loads for the frequency regulation is easy to understand. After the inception of a negative (positive) frequency deviation, a portion of responsive loads is turned off (on) to improve the real-power balance. Once the balance between real-power supply and demand is achieved, the frequency profile is maintained within its allowed range. In general, DR control results in quicker frequency regulation in comparison with conventional methods where change in valve or gate position takes considerable amount of time [1,13]. The communication delay in DR control mainly depends on the communication infrastructure and size of the power system. With the modern fast-paced communication network in smart grids, a quicker restoration of frequency profile with the help of responsive loads is achievable [9][10][11]21]. Figure 2 shows the model of a multi-machine power system augmented with battery energy storage system (BESS) and responsive loads connected through DR programs. The term ε represents the share of power to be adjusted with the help of responsive loads by DR program. The value of ε depends upon several factors including the price of electricity, network conditions and customer-side aspects. With reference to Equation (1) and Figure 2, the presence of responsive loads for frequency regulation results in Equation (5) where ∆P DR (s) represents the power contribution from DR program.

Responsive Loads
In the smart grids, responsive loads make an important part of the power system for long-term energy balance and to provide ancillary services. Integration of responsive loads is realized with the help of DR programs. There are two major categories of DR programs: 1. Dispatchable programs: where the system operator can manipulate the responsive loads according to power grid's needs. For example, direct load control is a type of such programs and also considered in this work. 2. Non-dispatchable programs: the programs in which the customer is financially encouraged to manipulate the loads but the system operator does not have direct control over the responsive loads. Its examples of are real-time pricing and peak time rebate [27].
The concept of manipulation of responsive loads for the frequency regulation is easy to understand. After the inception of a negative (positive) frequency deviation, a portion of responsive loads is turned off (on) to improve the real-power balance. Once the balance between real-power supply and demand is achieved, the frequency profile is maintained within its allowed range. In general, DR control results in quicker frequency regulation in comparison with conventional methods where change in valve or gate position takes considerable amount of time [1,13]. The communication delay in DR control mainly depends on the communication infrastructure and size of the power system. With the modern fast-paced communication network in smart grids, a quicker restoration of frequency profile with the help of responsive loads is achievable [9][10][11]21]. Figure 2 shows the model of a multi-machine power system augmented with battery energy storage system (BESS) and responsive loads connected through DR programs. The term ε represents the share of power to be adjusted with the help of responsive loads by DR program. The value of ε depends upon several factors including the price of electricity, network conditions and customer-side aspects. With reference to Equation (1) and Figure 2, the presence of responsive loads for frequency regulation results in Equation (5) Model of a multi-machine power system augmented with responsive loads and virtual inertia from BESS.

Energy Storage System
In the conventional power systems, massive machines of large power systems have inherent inertia which improves the system's stability. On the other hand, renewable energy systems in smart grids have no, or very small, rotating masses which cannot provide sufficient inertia as compared to the traditional power plants. So, the increased penetration of renewable energy resources poses a

Energy Storage System
In the conventional power systems, massive machines of large power systems have inherent inertia which improves the system's stability. On the other hand, renewable energy systems in smart grids have no, or very small, rotating masses which cannot provide sufficient inertia as compared to the traditional power plants. So, the increased penetration of renewable energy resources poses a challenge to the grid stability. A successful solution to the low inertia problem is proposed in the form of a control system called virtual synchronous machine that provides virtual inertia (VI) [13,28]. Such a control system utilizes BESS as a source of power and other main constituents are a power converter and a controller, as shown in Figure 3. The idea behind a virtual synchronous machine is to imitate the dynamic model of a real synchronous machine to improve the grid stability. From the operational point of view, the nominal state of charge of BESS is kept at 60~70% of rated capacity during normal conditions. In this way, BESS can absorb or inject power during contingency or emergency events. Thus, there are two main operating modes of a virtual synchronous machine as follows: 1.
Virtual generator mode: in this case, BESS provides power to the grid 2.
Virtual load mode: in this case, excess power from the grid is stored to BESS The transfer function (T VI (s)) between frequency deviation (∆ f (s)) and virtual synchronous machine's power change (denoted by ∆P VI (s)) is calculated according to Equation (6) as follows [25].
where, K I and K P emulate the inertia and damper winding characteristic of a synchronous generator. The value of K I is calculated such that the maximum defined rate of change of frequency makes the virtual synchronous machine to exchange its maximum power. In this study, maximum rate of change of frequency is considered as 1 Hz/s Hz. Similarly, the value of K P is calculated in a way that ∆P VI (s) is maximum when the frequency deviation reaches a predefined maximum value (which is 2.5 Hz in this study) [13,19]. Equation (6) signifies that though a virtual synchronous machine does not have rotating masses, its electrical appearance is similar to an electromechanical synchronous machine from the grid point of view. With reference to Figure 2, where the term ζ represents the share of power to be adjusted with the help of VI mechanism, the power balance equation is calculated as Equation (7).
Machines 2019, 7, x FOR PEER REVIEW 5 of 15 challenge to the grid stability. A successful solution to the low inertia problem is proposed in the form of a control system called virtual synchronous machine that provides virtual inertia (VI) [13,28]. Such a control system utilizes BESS as a source of power and other main constituents are a power converter and a controller, as shown in Figure 3. The idea behind a virtual synchronous machine is to imitate the dynamic model of a real synchronous machine to improve the grid stability. From the operational point of view, the nominal state of charge of BESS is kept at 60~70% of rated capacity during normal conditions. In this way, BESS can absorb or inject power during contingency or emergency events. Thus, there are two main operating modes of a virtual synchronous machine as follows: 1. Virtual generator mode: in this case, BESS provides power to the grid 2. Virtual load mode: in this case, excess power from the grid is stored to BESS  [13,19]. Equation (6) signifies that though a virtual synchronous machine does not have rotating masses, its electrical appearance is similar to an electromechanical synchronous machine from the grid point of view.
With reference to Figure 2, where the term ζ represents the share of power to be adjusted with the help of VI mechanism, the power balance equation is calculated as Equation (7).

Multi-Machine System with DR and ESS
The frequency deviation in a multi-machine power system augmented with DR and VI resources is calculated according to Equation (8) as follows: Referring to Figure 1 and Equation (4), ∆P i T (s) is calculated as: Combining Equations (8) and (9) and solving for ∆ f (s): where, The step load disturbance of magnitude ∆P d in time-domain, is mathematically written as Equation (12) by taking its Laplace transform as follows: Applying the final value theorem and solving for steady-state frequency deviation (∆ f SS ): where, Equation (13) indicates that, in addition to power from conventional resources, there are two other resources of power (i.e., responsive loads and ESS) which contribute to diminishing the effect of ∆P d (s). In other words, the system operator has extra degree of freedom to select the share between these three controls. The selection of each participant's share depends on several factors including electricity rate, network situation, capacity of ESS and customer-side factors. If K(s) is the total effort to counteract the effects of load disturbance and γ, ε and ζ are respectively the shares of conventional, DR and ESS controls, then the mathematical representation of this distribution is provided as follows [21]:

Case Studies
This section highlights the various features of the developed model with the help of numerical simulations performed on a 3-machine, 100MW-base test system whose parameters are shown in Table 1 [10,13,21].
Numerous simulations were performed to analyze the role of responsive loads and BESS on the system frequency response under a load disturbance of 0.05 p.u applied at t = 5 sec. A comprehensive description of 31 test cases is provided in Table 2.  The numerical study has been divided into 4 categories, each having 2 subgroups considering various configurations of the test system. For a broader perspective, these cases include the results for single-machine and multi-machine system with and without the consideration of responsive loads and energy storage systems. Machine-1 is utilized for the cases where a single machine operation is executed in the presence of responsive loads or ESS. The numerical simulations are performed in MATLAB/Simulink environment.

Conventional Model: Single-and Multi-Machine Cases
This case study considers a conventional multi-machine system and serves as a base case for comparison where responsive loads and ESS are not available for frequency regulation. A 3-machine test system equipped with non-reheat type turbines is shown in Figure 4. As detailed in Table 2, the response of individual and combined operation of machines is discussed in cases A1-A3 and B1-B4, respectively. Figure 5a shows the frequency response for individual machines when a load disturbance (∆P d = 0.05 p.u) is applied to the system at t = 5 sec. The response of Machine-2 (Case-A2) is slower due to higher values of governor and turbines time constants as compared to other machines. Response of Machine-1 is the fastest as its governor and turbine have lower time constants.
The frequency response for two or all machines when a load disturbance (∆P d (s)) is applied to the system is shown in Figure 5b. As expected, the lowest value for maximum frequency deviation is observed when all the three machines are operational as indicated by Case-B4. system frequency response under a load disturbance of 0.05 .
p u applied at 5sec t = . A comprehensive description of 31 test cases is provided in Table 2.
The numerical study has been divided into 4 categories, each having 2 subgroups considering various configurations of the test system. For a broader perspective, these cases include the results for single-machine and multi-machine system with and without the consideration of responsive loads and energy storage systems. Machine-1 is utilized for the cases where a single machine operation is executed in the presence of responsive loads or ESS. The numerical simulations are performed in MATLAB/Simulink environment.

Conventional Model: Single-and Multi-Machine Cases
This case study considers a conventional multi-machine system and serves as a base case for comparison where responsive loads and ESS are not available for frequency regulation. A 3-machine test system equipped with non-reheat type turbines is shown in Figure 4. As detailed in Table 2, the response of individual and combined operation of machines is discussed in cases A1-A3 and B1-B4, respectively. Figure 5a shows the frequency response for individual machines when a load disturbance ( 0.05 .   comprehensive description of 31 test cases is provided in Table 2.
The numerical study has been divided into 4 categories, each having 2 subgroups considering various configurations of the test system. For a broader perspective, these cases include the results for single-machine and multi-machine system with and without the consideration of responsive loads and energy storage systems. Machine-1 is utilized for the cases where a single machine operation is executed in the presence of responsive loads or ESS. The numerical simulations are performed in MATLAB/Simulink environment.

Conventional Model: Single-and Multi-Machine Cases
This case study considers a conventional multi-machine system and serves as a base case for comparison where responsive loads and ESS are not available for frequency regulation. A 3-machine test system equipped with non-reheat type turbines is shown in Figure 4. As detailed in Table 2, the response of individual and combined operation of machines is discussed in cases A1-A3 and B1-B4, respectively. Figure 5a shows the frequency response for individual machines when a load disturbance ( 0.05 .

Conventional Model Augmented with Responsive Loads
In this case study, responsive loads are also considered for frequency regulation in addition to a single/multi-machine system. A description of the test cases is provided in Table 2 (i.e. C1-C4 and D1-D4). Figure 6a shows the role of DR when only Machine-1 is available in the system. The impact of responsive loads in the presence of all three machines is shown in Figure 6b. The summary of important observations from this case study is presented as follows: 1. The presence of responsive loads helps in improving the frequency profile depending on the share ( ε ) of responsive loads participating for frequency regulation under DR program. 2. Depending on the availability and amount of responsive loads, the need for (spinning and nonspinning) reserves reduces as responsive loads can be manipulated according to grid's needs.

Conventional Model Augmented with Responsive Loads
In this case study, responsive loads are also considered for frequency regulation in addition to a single/multi-machine system. A description of the test cases is provided in Table 2 (i.e. C1-C4 and D1-D4). Figure 6a shows the role of DR when only Machine-1 is available in the system. The impact of responsive loads in the presence of all three machines is shown in Figure 6b. The summary of important observations from this case study is presented as follows: 1.
The presence of responsive loads helps in improving the frequency profile depending on the share (ε) of responsive loads participating for frequency regulation under DR program.

2.
Depending on the availability and amount of responsive loads, the need for (spinning and non-spinning) reserves reduces as responsive loads can be manipulated according to grid's needs.

Conventional Model Augmented with Responsive Loads
In this case study, responsive loads are also considered for frequency regulation in addition to a single/multi-machine system. A description of the test cases is provided in Table 2 (i.e. C1-C4 and D1-D4). Figure 6a shows the role of DR when only Machine-1 is available in the system. The impact of responsive loads in the presence of all three machines is shown in Figure 6b. The summary of important observations from this case study is presented as follows: 1. The presence of responsive loads helps in improving the frequency profile depending on the share ( ε ) of responsive loads participating for frequency regulation under DR program. 2. Depending on the availability and amount of responsive loads, the need for (spinning and nonspinning) reserves reduces as responsive loads can be manipulated according to grid's needs.

Conventional Model Augmented with ESS
A description of the test cases of this study is provided as E1-E4 and F1-F4 in Table 2. Here, ESS takes a part of frequency regulation in addition to a single/multi-machine system. The summary of important observations from this case study is presented as follows: 1. Figure 7a shows the role of ESS when only Machine-1 is available in the system. The presence of ESS improves the frequency profile in terms of maximum frequency deviation and settling time. Moreover, it can be deduced that the presence of virtual synchronous generator can reduce the size of conventional reserves.

2.
The impact of ESS in the presence of all three machines is shown in Figure 7b and depends on its share (ζ). Generally, higher the value of ζ, better the frequency response. Moreover, the response of ESS is quicker in comparison to the responsive loads and conventional schemes as well.

Conventional Model Augmented with ESS
A description of the test cases of this study is provided as E1-E4 and F1-F4 in Table 2. Here, ESS takes a part of frequency regulation in addition to a single/multi-machine system. The summary of important observations from this case study is presented as follows: 1. Figure 7a shows the role of ESS when only Machine-1 is available in the system. The presence of ESS improves the frequency profile in terms of maximum frequency deviation and settling time. Moreover, it can be deduced that the presence of virtual synchronous generator can reduce the size of conventional reserves. 2. The impact of ESS in the presence of all three machines is shown in Figure 7b and depends on its share ( ζ ). Generally, higher the value of ζ , better the frequency response. Moreover, the response of ESS is quicker in comparison to the responsive loads and conventional schemes as well.

Conventional Model Augmented with Responsive Loads and ESS
This test study involves both the responsive loads and the ESS with their various shares as detailed in Table 2. The observations for this study are in-line with the already discussed cases. The frequency response for these simulation cases are shown in Figure 8. In addition to frequency response, the changes in machine, DR and VI powers are also shown in Figure 9. The summary of important observations from this case study is presented as follows: 1. Due to presence of responsive loads and ESS, the response in this scenario is better than conventional control only and between the two already discussed cases. 2. Figure 9a shows the simulation results for Case-H1 where 3 machines contribute their powers to balance the impact of load disturbance. The slower response of Machine-2 as already discussed is again visible here. In this case, there is no contribution from responsive loads and ESS resources, hence their powers are zero. 3. Figure 9b shows the simulation results for Case-H4 where 50% contribution comes from responsive loads and ESS. The ESS has the fastest response as shown in Figure 9b. The

Conventional Model Augmented with Responsive Loads and ESS
This test study involves both the responsive loads and the ESS with their various shares as detailed in Table 2. The observations for this study are in-line with the already discussed cases. The frequency response for these simulation cases are shown in Figure 8. In addition to frequency response, the changes in machine, DR and VI powers are also shown in Figure 9. The summary of important observations from this case study is presented as follows: 1.
Due to presence of responsive loads and ESS, the response in this scenario is better than conventional control only and between the two already discussed cases. 2. Figure 9a shows the simulation results for Case-H1 where 3 machines contribute their powers to balance the impact of load disturbance. The slower response of Machine-2 as already discussed is again visible here. In this case, there is no contribution from responsive loads and ESS resources, hence their powers are zero.
3. Figure 9b shows the simulation results for Case-H4 where 50% contribution comes from responsive loads and ESS. The ESS has the fastest response as shown in Figure 9b. The contribution of powers from machines, DR and ESS at steady-state agrees with the model developed in Section 2. For example, the contribution of DR is calculated as: ε = 25% ⇒ ∆P SS DR = 0.25 × 0.05 = 0.0125p.u , which agrees with the results shown in Figure 9b. Similar calculations are valid for a power contribution from machines and ESS.   , which agrees with the results shown in Figure 9b.
Similar calculations are valid for a power contribution from machines and ESS.   Figure 10 highlights the contribution of DR and ESS in a comparative sense as explained below:

Comaprison
1. Figure 10a compares a 2-machine case (Case-B3) with a single machine system augmented with DR (Case-C4). The results show that the responsive loads can effectively handle the load disturbances even when there is partial shortage of conventional supply resources. It is notable that responsive loads cannot fully replace a generation unit if the magnitude of load disturbance is very high. Moreover, despite their benefits, DR programs affect the consumers' comfort level, so their contribution is limited. However, the presented comparison highlights the fact that responsive loads can play a vital role in frequency regulation in case of partial or temporary unavailability of generation resources. 2. Figure 10b compares a 3-machine case (Case-B4) with a single machine system augmented with ESS (Case-E4). The results show that the ESS can quickly restore the frequency deviation. Similar to the previous discussion, it is once again observed that a virtual synchronous machine can efficiently handle the load disturbances in the presence of limited generation resources.
Machines 2019, 7, x FOR PEER REVIEW 12 of 15 Figure 10 highlights the contribution of DR and ESS in a comparative sense as explained below:

Comaprison
1. Figure 10a compares a 2-machine case (Case-B3) with a single machine system augmented with DR (Case-C4). The results show that the responsive loads can effectively handle the load disturbances even when there is partial shortage of conventional supply resources. It is notable that responsive loads cannot fully replace a generation unit if the magnitude of load disturbance is very high. Moreover, despite their benefits, DR programs affect the consumers' comfort level, so their contribution is limited. However, the presented comparison highlights the fact that responsive loads can play a vital role in frequency regulation in case of partial or temporary unavailability of generation resources. 2. Figure 10b compares a 3-machine case (Case-B4) with a single machine system augmented with ESS (Case-E4). The results show that the ESS can quickly restore the frequency deviation. Similar to the previous discussion, it is once again observed that a virtual synchronous machine can efficiently handle the load disturbances in the presence of limited generation resources.

Conclusions
Modern power systems will have increased penetration of responsive loads under DR programs and VI provided by ESS. In the literature, there is limited work that focuses on the combined role of DR and VI in a multi-machine power system. The presented work adds the DR and VI loops to the basic model of a multi-machine system for the objective of frequency response analysis and control. The mathematical formulation for DR and VI powers and their role in steady-state frequency control is provided. It is observed that the DR and VI controls have the potential to provide the system operator with additional capabilities of frequency regulation with minimum time delays. Moreover, the results indicate that the ESS and responsive loads can possibly handle the load disturbances during a temporary or partial shortage of conventional power resources. It is, however, notable that due to socio-technical reasons, the ESS and responsive loads cannot replace a major part of the

Conclusions
Modern power systems will have increased penetration of responsive loads under DR programs and VI provided by ESS. In the literature, there is limited work that focuses on the combined role of DR and VI in a multi-machine power system. The presented work adds the DR and VI loops to the basic model of a multi-machine system for the objective of frequency response analysis and control. The mathematical formulation for DR and VI powers and their role in steady-state frequency control is provided. It is observed that the DR and VI controls have the potential to provide the system operator with additional capabilities of frequency regulation with minimum time delays. Moreover, the results indicate that the ESS and responsive loads can possibly handle the load disturbances during a temporary or partial shortage of conventional power resources. It is, however, notable that due to socio-technical reasons, the ESS and responsive loads cannot replace a major part of the generation. These reasons include stochastic nature of DR programs and customers' participation factor; expanding demand for electrical energy; efficiency and operational life of ESS and so on. For future research, the distribution of power shares among conventional, DR and VI controls can be optimized under various constraints.

Conflicts of Interest:
The authors declare no conflict of interest.