Linear Formulation for Short-Term Operational Scheduling of Energy Storage Systems in Power Grids

This paper presents linear programming (LP) formulations for short-term energy time-shift operational scheduling with energy storage systems (ESSs) in power grids. In particular, it is shown that the conventional nonlinear formulations for electric bill minimization, peak shaving, and load leveling can be formulated in the LP framework. New variables for the peak and off-peak values are introduced in peak shaving and load leveling model, and the historical peak value for demand charge are considered in the electric bill minimization model. The LP formulations simplify computation while maintaining the accuracy for including linear technical constraints of ESSs, such as the state-of-charge, charging/discharging efficiency, output power range, and energy limit considering the life cycle of ESS. Proposed LP formulations have been implemented and verified in practical power systems and a large-scale industrial customer using historical data.


Introduction
Energy storage systems (ESSs), both mechanical types, such as pumped-hydro energy storage (PHES), compressed-air energy storage (CAES) and flywheels, and battery types, such as Li-ion, NaS and lead-acid, are given opportunities to enter into power systems to supplement conventional generators as they show technological advances and price competitiveness.The ESSs can be applied for frequency regulation, energy arbitrage, peak-shaving, renewable resource integration, transmission congestion relief, demand response, and the energy management of end-user or micro grid [1,2].This paper is focused on developing optimal operational scheduling models of an ESS for short-term (minutes or hour basis) energy time-shift applications based on various objectives in a power grid as described in Table 1.
Strategies of ESS for energy arbitrage are determined in response to energy market prices, either in an investor or owner perspective [3].Moreover, in [4][5][6][7][8], the energy arbitrage of ESS is also tested for power systems with renewable resources.These linear formulated energy arbitrage models of ESSs can be flexibly applied in a transaction model with renewables considering uncertainty [4], operation problems of microgrid with other DERs [5][6][7], and the optimal size determination models of ESS considering the installed renewable resources [8].The production cost minimization model of ESS can be expressed in a form that the energy arbitrage model is included in the unit commitment problem by a system (or market) operator, and it can be effectively formulated as mixed integer linear programming (MILP) [9,10].Other applications of ESSs, except the energy arbitrage, are often solved by heuristic algorithms.Peak-shaving and load-leveling problems of ESSs are also tested in power grids from the interest of a system operator using dynamic programming in [11][12][13].The electric bill minimization problem in the customer side is formulated with nonlinear model considering both demand charge and energy charge, where dynamic programming, Markov decision processes, and particle swarm optimization are used in [14][15][16][17], respectively.Assuming that demand or renewable resource production is perfectly forecasted, the objective functions are defined with respect to the applications of an ESS, while technical constrains of ESS are commonly contained in the optimization model.The solution technique for an optimization problem is determined according to the type of the formulated problem.Often optimization problems of ESS have been formulated as non-linear problems and solved with heuristic and evolutionary algorithms, except the energy arbitrage for short-term optimal operation scheduling model of an ESS.The main disadvantage of these algorithms is to end up with suboptimal solutions.In contrast, if the same problem can be formulated with linear model without approximation, we can find the optimum solution efficiently and accurately using linear programming (LP) [18].
The objective of this paper is to develop linear expressions of optimal operational scheduling models of an ESS for all time-shifting applications in power grids.It was demonstrated that the energy arbitrage model from the investor or owner perspective as a price-taker in electricity market can be simply expressed in a linear formulation, as mentioned in the above literatures.Therefore, in this paper, we will first develop linear optimization models, which will ensure simplification as well as accuracy for peak-shaving and load-leveling in the system operator perspective.Moreover, we will develop a novel linear model for the electric bill minimization of end-users.Thus, we can find all optimum solutions for the time-shifting applications of ESSs listed in Table 1 using LP.
This paper is composed as follows: Section 2 presents the linear expression of general technical constraints of ESS.Section 3 presents linearly expressed objective functions and additional necessary conditions with respect to the energy-shifting applications of ESS.Proposed linear models are tested in Section 4 and conclusions are drawn in Section 5.

Linear Expression of ESS Technical Constraints
This section presents linear expression for constraints in the optimization problem by considering general technical characteristics of an ESS.The major constraints of an ESS are charging/discharging efficiency, limits of power output and the state-of-charge (SOC).Additionally technical constraints on the cumulated energy used during the operation period are included to consider the aging effect of an ESS.

Charging/Discharging Efficiency
where η j , ep j,t and EP j,t respectively are efficiency, AC power and DC power in MW; t and j represent the t-th time interval and the j-th EES unit; and superscripts c and d respectively represent charging and discharging.Energy loss occurs when an ESS are under charging or discharging.Constraints (1) and ( 2) ensure charging and discharging power of ESS-j considering AC/DC and DC/AC conversion efficiencies, respectively.

Output Power Limits
where EP c j and EP d j respectively are the maximum charging and discharging power outputs of ESS-j.Constraints (3) and ( 4) respectively represent the bounds on charging and discharging ranges of ESS-j at time interval-t.If an ESS has the minimum power outputs EP c j and EP d j , constraints (3) and ( 4) have to be changed as mixed integer linear expressions by introducing binary variables as follows: where u c j,t and u d j,t respectively represent the binary charging and discharging status of ESS-j at the t-th time interval.Constraint (7) ensures that the status of ESS-j is either under charge, discharge, or idle.
where ∆t is the time interval in hours, T is the number of time intervals over a decision horizon, SOC j,t is the SOC of ESS-j in MWh at the t-th time interval, SOC max respectively are ESS-j's lower and upper bounds of SOC at the final time interval.Constraint (8) determines the SOC of ESS-j at each time interval-t.Constraint (9) bounds ESS-j's SOC at each time interval-t.It is possible to adjust the range of the final SOC by constraint (10).

Total Energy Limits
One of the expressions for the lifetime of a battery ESS (BESS) is the cycle-life, which is the expected number of cycles for charging and discharging depending on the depth of discharge (DOD) of a BESS [19], which can be represented by the limited SOC range.When the possible SOC ranges of ESSs are fixed, the excessive number of round-trip operating cycles can be prevented by limiting total charged and discharged energies over an operational period with the following constraints: where α c j and α d j respectively are positive constants to limit the number of charging and discharging cycles of ESS-j over a decision horizon.
These two constraints are helpful not only for the BESS but also for mechanical ESS (MESS) operational scheduling, because, with them, we can find the least amount of the energy to meet the objectives, and prevent excessive wear and tear by choosing suitable values of α c j and α d j in the optimization models.

Linear Objective Functions for Applications in Power Grid
This section presents linear objective functions for the purpose of short-term ESS applications: energy arbitrage, peak-shaving, load-leveling and minimization of end-user's electricity bill.For energy arbitrage of ESSs from the investor or owner perspective, the existing linear energy arbitrage model as a price-taker in electricity market is already introduced in many existing literatures such as [3][4][5][6][7][8].Therefore, in this paper, we develop novel LP expressions, which will ensure simplification and accuracy for the electric bill minimization in an end-user perspective as well as the peak-shaving and load-leveling in the system operator perspective.

Peak-Shaving
The objective of peak-shaving is focused on reducing maximum load in a decision horizon as shown below: where N is the number of ESS, T is a decision horizon, D t is forecasted load demand at the t-th time interval, D t is a modified demand by ESSs at the t-th time interval, ep c j,t is a charging power of ESS-j from the grid at the t-th time interval, and ep d j,t is a discharging power of ESS-j injected to the grid at the t-th time interval.
The objective function ( 13) is a nonlinear function.For a linear formulation, we can modify the objective function (13) as following: where D max is a new variable which ensures that the peak value among the modified demands by ESSs is below this level in a decision horizon.Finally we can solve the peak-shaving problem using LP with the objective function (16) subject to ESS technical constraints (1)-( 12) and additional constraints ( 14), (15), and (17).

Load-Leveling
Strictly speaking, the peak-shaving and load-leveling should be distinguished because the goal of the load-leveling is to make the fluctuating demand flat by increasing the off-peak as well as decreasing the on-peak, while peak-shaving is only focused on the reduction of the peaks [20].Therefore to formulate as a linear model we can define an objective function of the load-leveling as following: where D min is a new variable which ensures that the off-peak value among the modified demands by ESSs is above this level in a decision horizon.The linear objective function ( 18) is to minimize the gap between the on-peak and off-peak demands subject to ESS technical constraints ( 1)-( 12) and additional constraints ( 14), ( 15), (17), and (19).

End-User Electricity Bill Minimization
The customer electric bill is composed of demand charge and energy charge.The demand charge is paid to the utility which operates and plans transmission and distribution network, and is based on the highest average load for 15 min during the billing period.To determine the demand charge, some electricity retail markets consider the historical peak loads within a specified period.For example, non-residential customers are served at the regulated time-of-use (TOU) rate by Korea Electric Power Corporation (KEPCO) in Korea, i.e., they are required to pay the demand charge for the highest load between the peak load during current billing month and the peak load of July-September and December-February during the prior 11 months [21].Energy charge is calculated for the amount of electricity consumption.Generally a customer have a right to choose an appropriate time-varying pricing among TOU, critical peak price (CPP), real-time-pricing (RTP), etc.
Therefore we can formulate the optimization model for electric bill minimization with ESSs in the end-user perspective as follows: where R D is the demand charge rate imposed on peak D max of the end-user, and R E t is the energy charge rate at the t-th time interval.The linear objective function ( 20) is subject to constraints ( 14), ( 15), (17), and ESS technical constraints (1)- (12).We can solve this model by LP because all of objective functions and constraints are expressed in linear forms.
If historical peak load D max,o is necessarily considered to calculate the demand charge, D max in the first term of the objective function ( 20) can be replaced with max(D max,o , D max ), then the objective function can be changed as following: By defining D max,g = max(D max,o , D max ), the optimization model for end-user's electric bill minimization with ESS can be generalized as the following linear formulation: where the linear objective function ( 22) is subject to constraints ( 14), ( 15), ( 17), ( 23), (24), and ESS technical constraints (1)-( 12).As mentioned above, the constant D max,o is the historical peak load during a billing period.In this generalized model, we can consider peak-load reduction or the energy-cost reduction, as priority, by setting D max,o = 0 or D max,o = ∞ (a large positive constant), respectively.In the peak-load reduction priority mode, the ESS charges when energy prices are low and discharges first when the loads are on-peak, then discharges extra charged energy at the other times when energy prices are high.On the other hand, ESS charges when energy prices are low and discharges when the prices are high in the energy-cost reduction priority mode.

Case Studies
Assuming that all estimation data are perfectly known, the proposed linear short-term optimal scheduling problems for the use of an ESS have been tested in this section.Peak-shaving and load-leveling in a system operator perspective are tested on a 500 MW/4000 MWh PHES.Since the round-trip efficiency of PHES is 0.75 in [2], we assume that the each of pumping and generating efficiencies is set to √ 0.75.The electric-bill minimization is implemented in a large industrial customer with a 4 MW/8 MWh Li-ion BESS, for which conversion efficiencies on charging and discharging are equally set to 0.95 [2].The proposed linear models are implemented in IBM CPLEX 12.6, which gives robust solutions for optimization problems [22].

Case-I: Peak-Shaving and Load-Leveling with PHES
Figure 1 is the demand for a week provided from the Korea Power Exchange (KPX) during 2-8 August 2010 [23], which is described in Table A1 and it has 3707 MW of off-peak and 6273 MW of on-peak demand.The leveled demand pattern corresponding to the time of energy prices shown in Figure 2 is obtained from the historical data from Korea, forecasted and posted by KPX [23].The energy prices are described in Table A2.
As mentioned above, the constant o D max, is the historical peak load during a billing period.In this generalized model, we can consider peak-load reduction or the energy-cost reduction, as priority, by setting (a large positive constant), respectively.In the peak-load reduction priority mode, the ESS charges when energy prices are low and discharges first when the loads are on-peak, then discharges extra charged energy at the other times when energy prices are high.On the other hand, ESS charges when energy prices are low and discharges when the prices are high in the energy-cost reduction priority mode.

Case Studies
Assuming that all estimation data are perfectly known, the proposed linear short-term optimal scheduling problems for the use of an ESS have been tested in this section.Peak-shaving and loadleveling in a system operator perspective are tested on a 500 MW/4000 MWh PHES.Since the roundtrip efficiency of PHES is 0.75 in [2], we assume that the each of pumping and generating efficiencies is set to 75 .0 . The electric-bill minimization is implemented in a large industrial customer with a 4 MW/8 MWh Li-ion BESS, for which conversion efficiencies on charging and discharging are equally set to 0.95 [2].The proposed linear models are implemented in IBM CPLEX 12.6, which gives robust solutions for optimization problems [22].

Case-I: Peak-Shaving and Load-Leveling with PHES
Figure 1 is the demand for a week provided from the Korea Power Exchange (KPX) during 2-8 August 2010 [23], which is described in Table A1 and it has 3707 MW of off-peak and 6273 MW of on-peak demand.The leveled demand pattern corresponding to the time of energy prices shown in Figure 2 is obtained from the historical data from Korea, forecasted and posted by KPX [23].The energy prices are described in Table A2.When the ESS is scheduled for peak-shaving and/or load-leveling, the adjusted on-peak D max and off-peak D mim need to be determined, which would result from the allowed maximum charging and discharging energy over a given planning horizon as defined in ( 11) and (12).Figure 3  When the ESS is scheduled for peak-shaving and/or load-leveling, the adjusted on-peak D max and off-peak D mim need to be determined, which would result from the allowed maximum charging and discharging energy over a given planning horizon as defined in (11) and (12).Figure 3 shows an experimental result of on-peak D max , off-peak D mim , and the differences between them (D max − D mim ) according to α c j and α d j for the weekly operation mode, where both α c j and α d j in ( 11) and ( 12) are set to the same value α in this case study.The optimization models are solved by weekly operation mode where each of the initial and final state of reservoir is set to 12.5% of the PHES capacity.In the figure, the adjusted on-peak D max and the difference (D max − D min ) are becoming constant when the values of α are 1.02 and 2.62, respectively, and these values are used for the peak-shaving and load-leveling models because full peak-shaving and load-leveling with the PHES can be most efficiently solved by weekly operation mode when the values of α are 1.02 and 2.62, respectively.When the ESS is scheduled for peak-shaving and/or load-leveling, the adjusted on-peak D max and off-peak D mim need to be determined, which would result from the allowed maximum charging and discharging energy over a given planning horizon as defined in (11) and (12).Figure 3 shows an experimental result of on-peak D max , off-peak D mim , and the differences between them (D max − D mim ) according to c j α and d j α for the weekly operation mode, where both c j α and d j α in (11) and (12)   are set to the same value α in this case study.The optimization models are solved by weekly operation mode where each of the initial and final state of reservoir is set to 12.5% of the PHES capacity.In the figure, the adjusted on-peak D max and the difference (D max − D min ) are becoming constant when the values of α are 1.02 and 2.62, respectively, and these values are used for the peak-shaving and load-leveling models because full peak-shaving and load-leveling with the PHES can be most efficiently solved by weekly operation mode when the values of α are 1.02 and 2.62, respectively.Dashed line and solid line shown in Figure 4a are the results of peak-shaving and load-leveling with the PHES, where the values of α are set to 1.02 and 2.62, respectively, in the optimization models for peak-shaving and load-leveling solved by weekly operation mode.Both peak-shaving and load-leveling have reduced the on-peak to 5840 MW.Moreover the load-leveling has raised the off-peak to 4284 MW while peak-shaving still remains at 3707 MW.Therefore, the load-leveling can reduce the demand range more than the peak-shaving.On the other hand, the advantage of the peakshaving model is that the energy required to reduce the peak is smaller than the load-leveling energy, even though it cannot control the off-peak.As the a result, cumulated charged (pumping) energies including energy losses from the grid for a week are 4108 MWh and 10589 MWh for peak-shaving and load-leveling, respectively, and cumulated discharged (generated) energies delivered to the grid are 3081 MWh and 7942 MWh, respectively.Each value is 0.75 of cumulated charged energy which is the round-trip efficiency of PHES as assumed above.Dashed line and solid line shown in Figure 4a are the results of peak-shaving and load-leveling with the PHES, where the values of α are set to 1.02 and 2.62, respectively, in the optimization models for peak-shaving and load-leveling solved by weekly operation mode.Both peak-shaving and load-leveling have reduced the on-peak to 5840 MW.Moreover the load-leveling has raised the off-peak to 4284 MW while peak-shaving still remains at 3707 MW.Therefore, the load-leveling can reduce the demand range more than the peak-shaving.On the other hand, the advantage of the peak-shaving model is that the energy required to reduce the peak is smaller than the load-leveling energy, even though it cannot control the off-peak.As the a result, cumulated charged (pumping) energies including energy losses from the grid for a week are 4108 MWh and 10589 MWh for peak-shaving and load-leveling, respectively, and cumulated discharged (generated) energies delivered to the grid are 3081 MWh and 7942 MWh, respectively.Each value is 0.75 of cumulated charged energy which is the round-trip efficiency of PHES as assumed above.
Since load forecasting necessarily has error caused by uncertainties, the operator can decide an alternative strategy of scheduling period to take consideration of the forecasting error.To reduce forecasting error, operating period should be shorter and scheduling cycle should be more frequent.If the optimization is solved by daily operation mode, we can get the results as shown in Figure 4b.In this case, for a week, total charged energies as the results of peak-shaving and load-leveling respectively are 24,172 MWh and 24,421 MWh including the losses, and total discharged energies to the grid are 18,129 MWh and 18,316 MWh, where α is set to 6.00 for the peak-shaving and to 6.06 for the load-leveling.These are minimal cumulated energy requirements which respectively can contribute to full peak-shaving and full load-leveling with the PHES.Since load forecasting necessarily has error caused by uncertainties, the operator can decide an alternative strategy of scheduling period to take consideration of the forecasting error.To reduce forecasting error, operating period should be shorter and scheduling cycle should be more frequent.If the optimization is solved by daily operation mode, we can get the results as shown in Figure 4b.In this case, for a week, total charged energies as the results of peak-shaving and load-leveling respectively are 24,172 MWh and 24,421 MWh including the losses, and total discharged energies to the grid are 18,129 MWh and 18,316 MWh, where α is set to 6.00 for the peak-shaving and to 6.06 for the load-leveling.These are minimal cumulated energy requirements which respectively can contribute to full peak-shaving and full load-leveling with the PHES.

Case-II: Electric Bill Minimization of a Customer with Li-ion BESS
The proposed linear electric bill minimization model was tested for a large industrial customer being served with the time-of-use (TOU) rate from KEPCO.The customer is charged with the demand rate of 7380 KRW/kW and the energy rate (KRW/kWh) adjusted on seasonal and hourly basis as shown in Figure 5 and described in Table A3 [21].As mentioned above, the peak-load applied to demand charge by KEPCO is determined as the maximum value during prior 12 months including the current month.As described in Figure 6 and Table A4, the customer load data in summer, assuming that it is constant for each hour, has different daily pattern, but has a recursive tendency on a weekly basis.The minimum and maximum values of the load are 0.66 MW and 15.15 MW, respectively.In addition, daily on-peaks in the weekend are larger than those in weekdays.Since load forecasting necessarily has error caused by uncertainties, the operator can decide an alternative strategy of scheduling period to take consideration of the forecasting error.To reduce forecasting error, operating period should be shorter and scheduling cycle should be more frequent.If the optimization is solved by daily operation mode, we can get the results as shown in Figure 4b.In this case, for a week, total charged energies as the results of peak-shaving and load-leveling respectively are 24,172 MWh and 24,421 MWh including the losses, and total discharged energies to the grid are 18,129 MWh and 18,316 MWh, where α is set to 6.00 for the peak-shaving and to 6.06 for the load-leveling.These are minimal cumulated energy requirements which respectively can contribute to full peak-shaving and full load-leveling with the PHES.

Case-II: Electric Bill Minimization of a Customer with Li-ion BESS
The proposed linear electric bill minimization model was tested for a large industrial customer being served with the time-of-use (TOU) rate from KEPCO.The customer is charged with the demand rate of 7380 KRW/kW and the energy rate (KRW/kWh) adjusted on seasonal and hourly basis as shown in Figure 5 and described in Table A3 [21].As mentioned above, the peak-load applied to demand charge by KEPCO is determined as the maximum value during prior 12 months including the current month.As described in Figure 6 and Table A4, the customer load data in summer, assuming that it is constant for each hour, has different daily pattern, but has a recursive tendency on a weekly basis.The minimum and maximum values of the load are 0.66 MW and 15.15 MW, respectively.In addition, daily on-peaks in the weekend are larger than those in weekdays.

Case-II: Electric Bill Minimization of a Customer with Li-ion BESS
The proposed linear electric bill minimization model was tested for a large industrial customer being served with the time-of-use (TOU) rate from KEPCO.The customer is charged with the demand rate of 7380 KRW/kW and the energy rate (KRW/kWh) adjusted on seasonal and hourly basis as shown in Figure 5 and described in Table A3 [21].As mentioned above, the peak-load applied to demand charge by KEPCO is determined as the maximum value during prior 12 months including the current month.As described in Figure 6 and Table A4, the customer load data in summer, assuming that it is constant for each hour, has different daily pattern, but has a recursive tendency on a weekly basis.The minimum and maximum values of the load are 0.66 MW and 15.15 MW, respectively.In addition, daily on-peaks in the weekend are larger than those in weekdays.Optimal scheduling results of the Li-ion BESS installed in an industrial customer obtained by the linear electric bill minimization algorithm for weekly operation mode are shown in Figure 7, where both initial and final SOCs of the Li-ion BESS have been set to 5% capacity of the rated energy.
To concentrate on saving the bill, α has been set to a large constant.For the weekly operating mode, adjusted on-peak load is decreased to 11.98 MW with the full discharging power, when the optimization is solved to satisfy the objective function ( 22) without consideration of the historical peak value by setting D max,o = 0. Since on-peak loads in the first four days (weekdays) are relatively lower than in the last three days (weekends), charge and discharge of the ESS are focusing on saving the energy cost during weekdays while they are focusing on reducing on-peak loads during the weekends.If the demand charge considers the historical peak-value and it is set to D max,o = 13 MW, the customer's on-peak decreases only to 13 MW and the ESS is scheduled to save the energy cost elsewhere.If the demand charge considers the historical peak-value but it is set to D max,o = 16 MW, the optimization is solved only to concentrate on saving the energy cost (arbitrage) since D max,o is greater than the peak load during a week.Hourly operating schedules of the ESS according to D max,o are shown in Figure 8.
The results for the daily operation mode are shown in Figures 9 and 10.Compared with the weekly mode, if the demand charge does not consider the historical peak value by setting D max,o = 0 on the first day of billing month, the ESS is scheduled to concentrate on reducing daily peak and D max,o is increasingly updated every day during the first week as shown in Figure 9a.The ESS is scheduled to reduce daily peak to 11.98 MW and to save the energy cost elsewhere on the day during remaining three weeks as shown in Figure 9b since the values of D max,o are equally set to 11.98 MW until the end of the billing month.If the demand charge considers the historical peak-value for the daily operation mode and it is set to D max,o = 13 MW and D max,o = 16 MW, the schedules of ESS are similar with those of the weekly operation mode while the used energies of ESS are equal to 74.12 MWh and 81.05 MWh.Optimal scheduling results of the Li-ion BESS installed in an industrial customer obtained by the linear electric bill minimization algorithm for weekly operation mode are shown in Figure 7, where both initial and final SOCs of the Li-ion BESS have been set to 5% capacity of the rated energy.To concentrate on saving the bill, α has been set to a large constant.For the weekly operating mode, adjusted on-peak load is decreased to 11.98 MW with the full discharging power, when the optimization is solved to satisfy the objective function ( 22) without consideration of the historical peak value by setting D max,o = 0. Since on-peak loads in the first four days (weekdays) are relatively lower than in the last three days (weekends), charge and discharge of the ESS are focusing on saving the energy cost during weekdays while they are focusing on reducing on-peak loads during the weekends.If the demand charge considers the historical peak-value and it is set to D max,o = 13 MW, the customer's on-peak decreases only to 13 MW and the ESS is scheduled to save the energy cost elsewhere.If the demand charge considers the historical peak-value but it is set to D max,o = 16 MW, the optimization is solved only to concentrate on saving the energy cost (arbitrage) since D max,o is greater than the peak load during a week.Hourly operating schedules of the ESS according to D max,o are shown in Figure 8.(a) Load profiles   The operation results and customer's cost savings for the two operation modes are summarized in Table 2.In both modes, the modified peak values for demand charge are the same as 11.98 MW, 13 MW, and 16 MW when D max,o = 0, D max,o = 13 MW, and D max,o = 16 MW, respectively.Compared with the weekly mode, the used energy for the daily operation mode decreases 1.42%, but total cost savings is 1.38 million KRW smaller due to the lower energy savings when D max,o = 0. On the other The operation results and customer's cost savings for the two operation modes are summarized in Table 2.In both modes, the modified peak values for demand charge are the same as 11.98 MW, 13 MW, and 16 MW when D max,o = 0, D max,o = 13 MW, and D max,o = 16 MW, respectively.Compared with the weekly mode, the used energy for the daily operation mode decreases 1.42%, but total cost savings is 1.38 million KRW smaller due to the lower energy savings when D max,o = 0. On the other hand, the used energies and cost savings for the daily mode are equal to those for the weekly mode when D max,o = 13 MW and D max,o = 16 MW.Total cost savings from the ESS can be expected to be 45.27 million KRW, which is the sum of demand charge savings of 15.89 million KRW and energy charge savings of 29.38 million KRW, and 5.57% reduction from the total electric bill without ESS when D max,o = 13 MW.Only energy savings of 31.12 million KRW can be expected due to no reduction of the peak value when D max,o = 16 MW.As a result from the schedule of weekly operation mode, reduced cycle-life of ESS is calculated to 38.21 cycles, 40.24 cycles, and 44.00 cycles during four weeks when each D max,o is set to 0 MW, 13 MW, and 16 MW, respectively.Reduced cycle-life from the schedule of daily operation mode is calculated to 37.67 cycles, 40.24 cycles, and 44.00 cycles when each D max,o is set to 0 MW, 13 MW, and 16 MW, respectively.

Conclusions
This paper presented linear programming models for short-term optimal operational scheduling of ESSs as the energy time-shifting application in a power system.As a follow-up to the existing linear energy arbitrage model in an ESS owner perspective, we developed linear programming framework for peak-shaving and load-leveling from the aspect of system operator, and for electric bill minimization model from the perspective of end users.The proposed linear models have been implemented and verified in a practical power system and a large-scale industrial customer using historical demand and energy prices.Under the forecast data, they can effectively solve the short-term operational scheduling problem of ESSs, promising that they can also be applied in other operational problems.The proposed linear peak-shaving and load-leveling problems can give a system operator an effective tool to secure the system reliability by making use of the operational scheduling models of grid-scale ESSs while existing energy arbitrage model can be included in traditional unit commitment problem solved by ISO.Moreover, the electric bill minimization model can be also applied in the energy management system for a grid-connected microgrid consisting of ESSs as well as other distributed resources to minimize production cost for industrial or commercial customers.Developing short-term operational scheduling models of the ESS considering stochastic linear formulation will be our future works.
minimum SOC limits of ESS-j, SOC o j is the initial SOC of ESS-j, and SOC lower

Figure 1 .
Figure 1.Demand for a week.

Figure 1 .
Figure 1.Demand for a week.

Figure 3 .
Figure 3. Adjusted on-peak D max and off-peak D min while increasing α for the weekly operation mode.

Figure 4 .
Figure 4. Adjusted demands for a week by peak-shaving and load-leveling models of the ESS.

Figure 4 .
Figure 4. Adjusted demands for a week by peak-shaving and load-leveling models of the ESS.

Figure 4 .
Figure 4. Adjusted demands for a week by peak-shaving and load-leveling models of the ESS.

Figure 5 .
Figure 5. Customer's energy rate (TOU) for a week.Figure 5. Customer's energy rate (TOU) for a week.

Figure 6 .
Figure 6.Customer's load profile for a week.

Figure 6 .
Figure 6.Customer's load profile for a week.

Energies 2017, 10 , 207 10 of 17 Figure 7 .
Figure 7. Weekly operation results for the customer electric bill minimization.Figure 7. Weekly operation results for the customer electric bill minimization.

Figure 7 .
Figure 7. Weekly operation results for the customer electric bill minimization.Figure 7. Weekly operation results for the customer electric bill minimization.

Figure 7 .
Figure 7. Weekly operation results for the customer electric bill minimization.

Figure 8 .
Figure 8. Weekly schedules of ESS for the customer electric bill minimization.Figure 8. Weekly schedules of ESS for the customer electric bill minimization.

Figure 8 .
Figure 8. Weekly schedules of ESS for the customer electric bill minimization.Figure 8. Weekly schedules of ESS for the customer electric bill minimization.

Figure 9 .
Figure 9. Daily operation results for the customer electric bill minimization which does not consider the historical peak-value.

Figure 9 .
Figure 9. Daily operation results for the customer electric bill minimization which does not consider the historical peak-value.

Figure 10 .
Figure 10.Daily operation results for customer electric bill minimization which considers the historical peak-value.

Figure 10 .
Figure 10.Daily operation results for customer electric bill minimization which considers the historical peak-value.

Table 1 .
Short-term energy time-shifting applications of ESSs in power systems.

Table 2 .
Customer's cost savings for the electric bill minimization during four weeks.

Table A4 .
Forecasted load of an industrial customer (MW).