Coordinated Control of Distributed and Bulk Energy Storage for Alleviation of Post-contingency Overloads

This paper presents a novel corrective control strategy that can effectively coordinate distributed and bulk energy storage to relieve post-contingency overloads. Immediately following a contingency, distributed batteries are implemented to provide fast corrective actions to reduce power flows below their short-term emergency ratings. During the long-term period, Pumped Hydro Storage units work in pumping or generation mode to aid conventional generating units keep line flows below the normal ratings. This problem is formulated as a multi-stage Corrective Security-constrained OPF (CSCOPF). An algorithm based on Benders decomposition was proposed to find the optimal base case solution and seek feasible corrective actions to handle all contingencies. Case studies based on a modified RTS-96 system demonstrate the performance and effectiveness of the proposed control strategy.


Introduction
Power systems have been operating in the way that requires the base case operating point can withstand an unexpected loss of components [1].However, the traditional preventive dispatch is conservative, very costly, and even infeasible to implement for potentially dangerous contingencies [2].Changing the dispatch paradigm from preventive to corrective would be useful to manage this situation,

OPEN ACCESS
as it can take into account the operator's ability to take corrective actions after an outage occurs [3].Corrective actions can eliminate post-contingency overloads on affected transmission lines, allowing the system operates with relaxed security but lower costs in the normal state.The primary tools for achieving this goal are programs for solving Corrective Security-constrained Optimal Power Flow (CSCOPF) problems [4][5][6][7][8].
Generation redispatch [3][4][5][6][7]9] is a commonly used corrective action incorporated in the CSCOPF to remedy post-contingency overloads, in which operators send orders to generators to increase/decrease their output in the event that any single generating unit or transmission line suddenly trips.However, this form of corrective action is only effective in the case of post-contingency power flows on lines are lower than their short-term emergency (STE) ratings [10], because most generators cannot adjust their output fast enough to relieve the STE overloads within several min, due to their slow ramp up/down rates [11].In the literature, many researchers [4][5][6][7] have not considered the STE limits in their CSCOPF models, which eases the real-time computational burden, but ignoring such short-term emergency overloads may lead to cascading line tripping before the operator has redispatched generators [12].In current practice, if there are no fast-response corrective actions available in the system, the generating units should be dispatched out-of-merit in the pre-contingency state to make sure that no post-contingency power flows on the affected lines would surpass their STE ratings [2,[12][13][14][15].This preventive/corrective combined measure can guarantee security during the post-contingency short-term state, whereas the preventive actions raise the operating costs.
Load shedding [16,17] is considered as the last resort for operators to restore system security when other corrective actions are not sufficient.If this desperate action is executed, significant amounts of loads would be shed in critical areas.To mitigate the amount of lost loads, new techniques and devices are necessarily introduced into the network to extract excess generation or provide back-up power support for relieving N-1congestion.
Energy storage (ES) can deliver multiple benefits that enhance grid performance, such as load following, peak shaving, spinning reserve, stability enhancement and power quality improvements [18,19].With the recent rapid development of energy storage technologies and their flexibility during operation, interest in introducing energy storage to remove transmission congestion has been growing [20].The investigation by U.S. EPRI [21] indicates the technical feasibility and potential benefits of energy storage to remedy N-1 emergency overloads in thermally constrained networks.In [22], an enhanced SCOPF formulation was studied, which incorporates the battery energy storage systems for relieving transmission emergency congestion.The combination of energy storage and FACTS devices was proposed to allow higher power transfer levels [23] and minimize power not served after the severe fault conditions [24].
In this paper, we focus on the use of different types of energy storage as corrective control resources to remove post-contingency overloads.Since the effectiveness of different energy storage technologies for corrective control is constrained not only by their ramp speed, but also the power and energy capacity that they are able to inject or store to cope with a contingency [25], the control of various forms of energy storage to relieve overloads has to be coordinated.The main contribution of this work is two-fold: (1) A novel control strategy that can effectively coordinate distributed and bulk energy storage to relieve post-contingency overloads is proposed.Under the normal operating condition, this paradigm would reduce generation costs and increase the transmission capability.Following a contingency, control on a short-term and long-term timescale would reallocate the distributed and bulk energy storage to alleviate the violations.(2) This problem is formulated as a multi-stage CSCOPF, and an algorithm based on Benders decomposition was proposed for the solution of the optimal base case dispatch and feasible corrective actions.
The remainder of this paper is organized as follows: Section 2 explains the coordinate control strategy of distributed and bulk energy storage.Section 3 details the mathematic model formulation.The solution methodology is given in Section 4. Section 5 presents the test results of case studies, and finally, Section 6 states the conclusions.

Post-Contingency Timeline
Every conductor used in transmission networks has associated thermal ratings [10], such as normal (continuous), STE, and long-term emergency (LTE) ratings.As in CSCOPF, the LTE rating is usually set the same as the normal rating, given the STE and normal rating for a circuit, the post-contingency timeline can be divided into two periods: (1) Short-term period: Time period from the end of the transient stability region until the power flows on overloaded lines are brought back within their long-term ratings.This time period can only be last for a very short time period (typically 15 min).(2) Long-term period: Time period during which the power flows on overloaded lines are limited below the long-term (normal) ratings.The time period is typically 30 min to several hours.
An understanding of the post-contingency timeline will help in choosing which storage technologies are best for a specific corrective control application.

Available Storage Technologies for Corrective Control
In some severe contingency scenarios, the power flow on an affected line may exceed its normal and STE ratings.Distributed ES (i.e., batteries) and bulk ES, such as Pumped Hydro Storage (PHS), Compressed Air Energy Storage (CAES), and some types of large batteries, can be controlled to correct the short-term and long-term violations.In this paper, two different energy storage systems, namely batteries and PHS are illustrated adapting to the proposed control strategy.Table 1 makes comparison of the key characteristics between current battery storage and PHS [25].Batteries are one of the most cost-effective energy storage options available, which normally have limited power and energy capacity but with very fast ramp rate, capable of transitioning from zero to full output in microseconds to several seconds [26,27].In addition, batteries can be containerized and easily installed in distributed applications across the transmission network [28].Contribution to these characteristics, distributed batteries are suitable for providing fast corrective actions in the post-contingency short-term period.
PHS is the most widely implemented grid-scale energy storage, which typically provides hundreds to thousands of megawatts of capacity in a single facility [29,30].Originally, these units were used for daily energy management [31,32], while the newer adjustable speed system design allows PHS with a flexible ability to adjust their input and output power in the pumping mode as well, PHS can thus provide ancillary services (frequency regulation, contingency reserves, etc.) or support the integration of renewable energy sources, which greatly increases overall plant efficiency [18,29,33,34].We only consider the PHS units with adjustable speed system for enhancing system corrective security.Compared to the battery storage, PHS units have relatively slower response speed but much larger power and energy capacity, could be utilized to alleviate post-contingency long-term congestion.

Coordinated Control Strategy
The operator must deploy the distributed and bulk energy storage in a timely manner to clear the post-contingency short-and long-term violations.To coordinated control these storage units for relief of transmission N-1 congestion, two types of corrective control are proposed: (1) Short-term Corrective Control: During the short-term period, distributed batteries are implemented to provide fast corrective actions.In this way, the immediately post-contingency flows could be larger than the STE rating, however the violations can be alleviated instantly by batteries.This control method ensures that the post-contingency system be stable to endure until the generators and PHS units have completely adjusted their output.(2) Long-term Corrective Control: During the long-term period, PHS units work in pumping mode to absorb excess generation, or in generation mode to provide back-up power support.These corrective actions aid generating units keep line flows below the normal ratings.In this way, load shedding could be mitigated or avoided under some severe contingency conditions (i.e., failures of a large generator or a major transmission line).
The detailed process of the proposed coordinated control strategy of distributed batteries and PHS units to relieve the post-contingency overloads is given as follows.To help understand, Figure 1 illustrates the evolution of power flow on an overloaded line (F 0 in this figure represents the pre-contingency flow).
Step 1: Immediately after the fault occurs at time t 0 , distributed batteries are controlled instantly to discharge or charge power to bring the post-contingency line flows (F') back down with their short-term emergency ratings (F STE ).Then, the power injections from distributed batteries will remain constant until time t 1 when generators and PHS units start ramping.
Step 2: During the ramping period from time t 1 to t 2 , generators and PHS units adjust their output continuously.On the other hand, distributed batteries reduce continuously their injections and extractions, and consequently stop the charging/discharging process at time t 2 .During this time interval, the power flows on the overloaded lines decrease linearly until they reach their normal ratings (F max ).
Step 3: During the long-term period (after t 2 ), the output of generators and PHS units keep constant to prevent the power flows exceeding the normal ratings.By implementing the proposed coordinated control, the post-contingency overloads on a transmission line could be corrected within the perspective of short-and long-term timeframes.To fulfill this control strategy, two key questions are required to be resolved: (1) How much power should be discharged or charged instantly by batteries after a contingency occurs?(2) How much power the PHS units should produce or absorb to aid generators removing long-term violations at time t 2 (the operating point at the end of the short-term period)?
Given corrective actions at these two operating points, the post-contingency flows then could be controlled below the STE ratings and decreasing linearly during the ramping period (t 1 to t 2 ).Since security during this period is guaranteed, constraints associated with the ramping process are not necessarily considered in the CSCOPF model.

Problem Formulation
This problem is formulated as a three-stage CSCOPF formulation, explicitly modeling the timeframe to dispatch storage units to comply with the post-contingency security limits.The structure of this problem is shown in Figure 2. The first-stage problem determines the cheapest pre-contingency generation schedule with no continuous contribution from the energy storage units.For a given base case solution, the latter two-stage problems handle the post-contingency conditions, relieving the short-and long-term violations, respectively.Corrective actions in the post-contingency conditions involve those from distributed and bulk storage units, as well as those from conventional generators.For simplicity sake, the pre-and post-contingency power flows are calculated using the DC power flow method, and generator outages are not considered.It is assumed that the energy stored and spare energy capacity available in each storage unit is large enough to cope with the contingencies.The best location for siting a battery to remove post-contingency short-term overloads is determined through the sensitivity analysis method given in [21].

Objective Function
The objective function is to minimize the base case generation cost of the generators, which is calculated using the quadratic cost functions: where PG i 0 is the power produced by generator i. a i , b i , c i are the coefficients of cost function of generator i. N G is the set of generators.Superscript 0 represents the base case.
The associated operating costs of storage units under contingency conditions are not included in the objective function [25,35,36].The main reasons are given as follows: (1) Energy storage normally have a high investment while relatively low operating cost [28].
(2) The probability of a contingency would occur is very small (close to zero), therefore, the expected operating costs for storage units to relieve overloads are not necessarily considered.(3) Since the actual implementation of corrective actions occurs in real time, the proposed model is concerned about the feasibility of distributed and bulk storage units to comply with the post-contingency security, but not the optimal operating costs of such storage units to relieve overloads.
The objective function is subject to the pre-contingency, as well as the post-contingency short-and long-term security constraints.

Power Balance Equation
The total generation should meet the load demand: where PG i 0 is the real power output of generator i; PL l is the load at bus l; N L is the set of load buses.

Base Case Transmission Limits
At base case, the power flows must not exceed the normal limits: where T 0 is the power transfer distribution factor (PTDF) matrices for the base case.PG 0 , PL are respectively the vector of generator power output and load injections for the base case.

Power Output Limits
The power output of a generator should be within its limits: where PG i min and PG i max are the minimum and maximum real power output of generator i.

Power Balance Equation
The power output of generators remains the same immediately after the transmission contingency occurs.Thus, the power injections and extractions of the distributed batteries should be balanced to maintain the power balance in the system: where PBD m k , PBC m k are the discharging power and charging power of battery m to deal with the short-term violations following contingency k. k is an index of the contingency set (N C ). N B is the set of distributed batteries.

Short-Term Transmission Limits
Immediately following a contingency, distributed batteries respond instantly to bring line flows back within their short-term emergency ratings: where γ is the vector of factors relating the short-and long-term ratings of the branches.T k is the PTDF matrices for the k-th contingency.PBD k , PBC k are the vector of power that batteries must discharge and charge to deal with contingency k.

Power Limits of the Distributed Batteries
The battery is capable of transitioning from zero to full output continuously: where PBD m max , PBC m max are the maximum discharging and charging power limits of battery m.

Power Balance Equation
In the long-term period, the system power balance equation is as follows: where PG i k is the real power output of generator i after redispatching.PSD n k , PSC n k are the dispatch power of the n-th PHS unit in generation and pumping mode to deal with the long-term violations following contingency k.N S is the set of PHS units.

Long-Term Transmission Limits
The corrective actions from generators and PHS units should reset the power flows on each line be within its normal limit: max ( ) where PSD k , PSC k are the vector of power that PHS units must inject to or extract from the grid to deal with contingency k.

Generator Ramping Constraints
The redispatching amount of a generator should be within its ramping limits: where is the maximum possible redispatch of generator i during the ramping period.

Generation Limits
The redispatched generation of a generator should be within its limits:

Online Status of the PHS Units
The online status of a PHS unit is modeled as follows: where UD n k , UC n k are binary variables.UD n k is 1 if the PHS unit works in generation mode for contingency k, and 0 if it is either idle or pumping.UC n k is 1 indicating the PHS unit working in pumping mode, and 0 if it is either idle or generating.

Power Limits of the PHS Units
The power a PHS unit releases or absorbs should be within its limit: where PSD n min , PSD n max are the minimum and maximum generation power limits of PHS unit n, and PSC n min , PSC n max are the minimum and maximum pumping power limits of the n-th PHS unit.

Solution Approach
The proposed model is a large-scale mixed-integer programming (MIP) problem.It is difficult to solve the problem directly due to its high computing dimensionality, especially when the system is large and many N-1 contingencies are considered.Alternatively, Benders decomposition [37][38][39] can be used to solve this problem.The primal problem is decomposed into a master problem linked with two sets of Sub-problems.The master problem finds an optimal solution under normal state, which is an ordinary DC-OPF formulation augmented by the feasibility Benders cuts generated by the two sets of sub-problems.Sub-problems of type 1 and type 2 correspond to the second-and third-stage, seeking feasible corrective actions to remove the post-contingency short-and long-term overloads.It should be noted that non-negative slack variables are introduced to sub-problems 1 and 2 to ensure the problems are feasible.See the Appendix section for detailed models of the master and two types of sub-problems.
The master problem is solved using standard quadratic programming.Sub-problems 1 are linear, and can be solved using linear programming.Sub-problems 2 involve binary variables, the Branch-and-Cut method [40] is used to obtain feasible solution within a small duality gap.If the duality gap is set to zero, sub-problems 2 will obtain the optimal solution with a longer computation time.
The two sets of sub-problems are solved iteratively with the master problem until all short-and long-term post-contingency violations are removed using corrective actions (zero slack variables achieved).If, after solving a sub-problem, the corresponding slack variables are not equal to 0, then it is labeled an uncontrollable contingency, a feasibility Benders cut must be generated and added to the master problem for mitigating the violations in the next iteration: where f k is the objective function value of the sub-problem, PG i 0* is the base case trial operating point obtained from solving the master problem, λ i is the multiplier associated with the i-th generator output in sub-problem 1 or sub-problem 2, and can be determined as follows: The controllable contingencies are handled by corrective actions by storage and generating units without requiring any revisions to the existing base case solution.
The flowchart of the Benders decomposition based algorithm for solving the proposed multi-stage CSCOPF formulation is shown in Figure 3.
Step 1: Solve the master problem to determine an initial operating point; set the contingency index K = 1.
Step 2: For contingency K, calculate the post-contingency power flows.
Step 3: Check if the post-contingency flows are within their short-term emergency limits, if not, solve sub-problem 1 and generate a Benders cut for each sub-problem that is not feasible.
Step 4: Check if all the post-contingency flows are within their long-term limits; for those that are not, solve a sub-problem of type 2 and generate a Benders cut when the problem is not feasible.
Step 6: When K > NC, if any sub-problems are not feasible, add the Benders cuts generated for each infeasible sub-problem to the master problem and repeat steps 1 to 5. Otherwise, stop and print solutions.

Case Studies
The proposed coordinated control strategy and algorithm have been tested using a modified RTS-96 system as shown in Figure 4, which consists of three interconnected areas (RTS-79 network), with 96 generators, 73 buses, and 120 lines.The modified line parameters of the RTS-96 system are given in Table A1 in the Appendix section.The STE rating of all lines is assumed to be 1.2 times larger than the normal rating.We consider all N-1 line outage contingencies as constraints of the CSCOPF model.All the experiments were performed on a personal computer with 4 Intel (R) Core (TM) i7-4700MQ CPU (2.4 GHz) and 8 Gb of memory.The programs are implemented using Matlab and To analyze the performance of the proposed coordinated control strategy, the following three cases were considered:

Case 1:
Traditional CSCOPF without distributed or bulk energy storage installed in the system.Generators are preventively controlled to prevent the post-contingency flows exceeding the STE ratings.Case 2: CSCOPF with distributed batteries.12 batteries each rated at 10 MW/10 MWh are distributedly installed in the system to provide short-term corrective actions.Only conventional generators are available to provide long-term corrective actions.The buses location of these batteries are: (106,108,113,114), (206,208,213,214), (306,308,313,314).Case 3: CSCOPF with distributed batteries and PHS units.The short-term corrective actions come from 12 batteries as described in Case 2. Both generators and 2 PHS units could provide long-term corrective actions.The location and power limits of the PHS units are given in Table 2.

Results Obtained in Different Load Level
The simulation is performed for several load levels, in order to estimate the impact of distributed and bulk energy storage at different load levels.The base case load is 6122.2MW, and the load level is increased in steps of 2.5%.
Tables 3 to 5 show the results obtained by the three cases at each load level.NSC and NLC are the numbers of contingencies that result in a violation of the short-term and long-term limits, respectively.As increasing the load level, all the tables indicate an increase in the generation cost, and the CSCOPF models have to deal with more contingencies that violate the shot-term and long-term limits.The last column of Tables 3 to 5 shows that the computation time for the three types of CSCOPF increases with the load level.The CSCOPF without storage units is relatively costly.As no batteries were installed in the system, preventive actions have to be applied in the pre-contingency state to make sure that no short-term violations would occur, which increases the base case generation cost.On the other hand, since the model relies on generators to guard all the long-term security, under stress conditions (load level > 1.05), the overloads cannot be completely removed, load shedding have to be implemented to handle the contingencies, otherwise the program would be infeasible.
As shown in Table 4, the use of distributed batteries to prevent the power flows exceeding the STE ratings decreases the generation cost.For example, if the load level is 1.05, the generation cost obtained by the CSCOPF with batteries is $141,378, saving 0.3% compared to that without batteries ($141,812).It allows the system operates in a cheaper pre-contingency state, however generators are required to deal with more long-term violations.Note that distributed batteries in this model would rather aid generators to alleviate the long-term violations, nor further improve the system's maximum load ability.
Table 5 gives the results obtained by the CSCOPF when both distributed batteries and PHS units are available in the grid.Similar to Case 2, the post-contingency short-term overloads can be removed using corrective actions from the batteries.Moreover, the CSCOPF is able to withstand a much larger range of loads, owing to the PHS units.No loads are required to be shed or curtailed if the load level is not larger than 1.175.The maximum load level of the system increased 11.9% compared to those of Case 1 and Case 2. Therefore, the implementing of PHS units for long-term corrective control would greatly improve the system load ability and reduce the need for transmission expansion.

F
As can be seen, if no storage units are placed, the pre-contingency loading level on line L 118 is 0.91.After the outage of line L 107 , the immediately flow becomes 1.15, which is lower than the STE rating (1.2).In the long-term period, the flow (0.98) is adjusted to below the continuous rating (1.0).
If both distributed and bulk storage are installed, the pre-contingency loading level on line L 118 is increased to 0.96.After the outage of line L 107 , the immediately flow (1.28) on line L 118 would violate the STE rating.However, this flow is first reduced to the STE rating using the distributed batteries and then the continuous rating through power adjustment of the conventional generating units and PHS units.Hence, the coordinated control of batteries and PHS units provides means for efficient operation in the post-contingency state while still maintaining security in a robust way.It can thus enhance the transmission capability of certain transmission corridors, allowing the lines operating with higher power flows during the pre-and post-contingency state.As can be seen from Table 6, 71 lines in the RTS-96 network show an increase in their pre-contingency transmission capacity, including 27 lines in Area 1, 20 lines in Area 2, 21 lines in Area 3, and three tie-lines (L 24 , L 118 , L 120 ) connect those two areas.

Corrective Actions of CSCOPF with Distributed and Bulk Storage
The short-term and long-term corrective actions in Case 3 (the load level is 1.05) are given in this part to further explain how the post-contingency overloads were cleared using coordinated control of distributed and bulk storage units.
Figure 6 shows the corrective actions provided by the 12 batteries immediately after the outage of line L 107 .It can be seen that, batteries located in area 3 only need to provide corrective actions in the form of discharge, while batteries placed in area 1 and area 2 need charge or discharge.This is because that, to remove the post-contingency emergency overloads, batteries located downstream from the overloaded lines will inject power in the network, and other batteries located upstream from these lines might extract power from the network.Table 7 illustrates that 12 contingencies, out of those that result in a violation of long-term limits (NLC = 34), require PHS units to provide corrective actions.In the case of the outage of line L 27 , L 28 , or L 41 , the PHS unit (PHS 2 ) located at bus 317 in area 3 has to work in generation mode to provide back-up power to aid generators, while under other contingency conditions (L 66 , L 67 , L 72 , L 86 , L 103 , L 104 , L 105 , L 106 , L 117 ), the PHS units located in area 2 or area 3 would work in pumping mode to absorb excess power from the network.If the long-term period is set to 1 hour, it would require that 50 × 1 = 50 MWh spare energy capacity be available in PHS 1 , while 230 MWh energy should be stored and 225 MWh energy margin in PHS 2.

Effect of the Power Capacity of Distributed Batteries
The effect of the power capacity of distributed batteries on the generation cost is studied here (Case 3 with load level equals to 1.05).Figure 7 shows how the generation cost varies as the power capacity of the batteries.The power capacity of each battery was varied from 0 to 50 MW.As one would expect, lower generation cost could be achieved as increasing the power capacity of batteries.The choice of how much power capacity of a battery is adopted for correcting short-term overloads would depend on its investment costs.

Effect of the Number and Power Capacity of PHS Units
To analyze the effect of CSCOPF with different number of PHS units, three cases were considered: • Case 4: 1 PHS unit is installed at bus 117 in area 1.
• Case 5: 2 PHS units located at buses 117 and 217 (Area 2).• Case 6: 3 PHS units located at buses 117, 217, and 317 (Area 3).For all the cases, each PHS unit has the same power capacity with PHS 2 as given in Table 2, the system base case load is 6122.2MW.It can be seen from Table 8, if more PHS units in the network are used for post-contingency corrective control, the system could have a larger load level, and the CSCOPF allows more long-term violations to be removed by corrective actions using storage units.The power capacity of each PHS unit in case 6 was increased in steps of 40 MW.To help understand the impact of the power rating of PHS units on the generation dispatch, the expected amount of generation redispatch (E G ), and the expected amount of power adjustment of PHS units (E S ) are calculated: C S ( ) where ρ k is the probability of contingency k, which is calculated from historical failure rates [41].Table 9 indicates that, at base case, the expected amount of power adjustment of PHS units increases as the power capacity of PHS units increases, while the expected amount of generation redispatch decreases.This demonstrates that if larger PHS units are available in the system, less post-contingency generation redispatch is required, as more corrective actions can be provided by the PHS units to move back the power flows on affected lines below their normal ratings.Table 9 also shows that the increasing of power capacity of the PHS units impacts the system's overall load ability.

Conclusions
A novel corrective control strategy that can coordinate distributed and bulk energy storage for relief of post-contingency overloads has been presented in this paper.This problem is formulated as a multi-stage CSCOPF incorporating the distributed batteries and PHS units.Immediately after a contingency, batteries are used to provide fast-response corrective actions to prevent the power flows on affected lines exceeding their short-term emergency ratings.In the long-term period, PHS units work in generation or pumping mode to aid generators bring the flows down within the normal limits.
An algorithm based on Benders decomposition was proposed to solve the proposed CSCOPF model.The primal problem was decomposed into a pre-contingency master problem linked with two sets of post-contingency sub-problems.The master problem determines the optimal base case solution, while the two types of sub-problems seek feasible corrective actions to handle all contingencies.
Test results on a modified RTS-96 system demonstrate that the proposed control strategy offers the following advantages: (1) Coordinated control of distributed batteries and PHS units following an outage could effectively remove post-contingency overloads and guarantee system operational reliability.(2) It lightens the requirement of preventive/corrective actions from conventional generators, thus decreases the generation costs.(3) It allows the system operates with higher pre-and post-contingency power flows, therefore reinforces the available transfer capabilities.It would thus reduce the need for investments in additional or upgraded transmission lines.
Although this paper is focused on the intraday single operating point, the proposed coordinated control method can also be extended to the multi-period CSCOPF or day-ahead Security-constrained Unit Commitment problems.These topics are left for future work.

A2. Sub-Problem 1
Slack variables s k are added to the post-contingency short-term transmission limits to ensure the sub-problem 1 is feasible.The objective of sub-problem 1 is to minimize the sum of slack variables, subject to the post-contingency shot-term security constraints: Minimize: where N T is the set of transmission lines.

A3. Sub-Problem 2
Slack variables r k are added to the generator ramping constraints to ensure the sub-problem 2 is feasible.The objective of sub-problem 2 is to minimize the sum of slack variables, subject to the post-contingency long-term security constraints: Minimize:

Figure 1 .
Figure 1.Power flow on an overloaded line.

Figure 3 .
Figure 3. Flowchart of the proposed algorithm.

Figure 5
Figure 5 shows the evolution of the power flow on line L 118 (the line connects buses 121 and 325) in the case of the outage of line L 107 (the line connects buses 316 and 317).The dashed line in the figure represents the pre-and post-contingency flows obtained by the CSCOPF without storage (Case 1), the solid line represents the flows obtained by the CSCOPF with batteries and PHS units (Case 3).In both cases, the load level is 1.05.

Figure 5 .
Figure 5. Evolution of the power flow on line L 118 .

Figure 6 .
Figure 6.Short-term corrective actions of distributed batteries to cope with the outage of line L 107 .

Figure 7 .
Figure 7. Minimum cost achieved by the CSCOPF as a function of the power capacity of batteries.

Table 1 .
Key technical characteristics of batteries and PHS.

Table 2 .
Location and power limits of the PHS plants.
ID BUS PSD

Table 3 .
Results obtained in Case 1.

Table 4 .
Results obtained in Case 2.

Table 5 .
Results obtained in Case 3.

Table 6 .
List of lines those with an increase in their loading level.

Table 7 .
Long-term corrective actions of PHS units to cope with all contingencies.

Table 8 .
Results obtained in three cases with different number of PHS units.

Table 9 .
Results obtained in case 6 with different power capacity of PHS units.

Table A1 .
The modified line parameters of the RTS-96 system are given in the following table: Cont.