Optimal Component Sizing for Peak Shaving in Battery Energy Storage System for Industrial Applications

Rodrigo Martins 1,* , Holger C. Hesse 2 , Johanna Jungbauer 3, Thomas Vorbuchner 2 and Petr Musilek 1,4 1 Electrical and Computer Engineering, University of Alberta, Edmonton, AB T6G 1H9, Canada; pmusilek@ualberta.ca 2 Department of Electrical and Computer Engineering, Technical University of Munich (TUM), 80333 Munich, Germany; holger.hesse@tum.de (H.C.H.); thomas.vorbuchner@onlinehome.de (T.V.) 3 Smart Power GmbH & Co KG, 80333 Munich, Germany; jungbauer@smart-power.net 4 Electrical Engineering and Computer Science, VSB-Technical University Ostrava, 70800 Ostrava, Czech Republic * Correspondence: rcmartin@ualberta.ca; Tel.: +1-780-492-5368


Introduction
In power systems, the load profile can be characterized by the "peak load times" of the system-short periods of time when large amounts of power are required [1].The peak load periods can occur at different times during the day, depending on the season of the year and the load composition (residential, commercial, or industrial).Peaks of demand impact the network planning because the electrical infrastructure of transmission and distribution systems must be designed to support the maximal demand of the system [2].For this reason, the electrical power grid infrastructure may be underutilized most of the time, reaching its loading capacity limit at only a few moments of the year.Consequently, commercial and industrial customers are charged not only by their total energy consumption but also by their highest power demand that dominates the grid construction costs.The electricity charge can be discriminated in subcomponents like the generation cost, taxes, and fees which represent a small portion of the total electricity payment of the customers.Accordingly, commercial and industrial customers are interested in decreasing energy and power costs, which are the most significant part of the total charges, without lowering their energy consumption.In this context, energy storage systems (ESS) can be used to help customers flatten their demand profile by storing energy during off-peak periods and releasing it during peak load periods.
The deployment of ESS can achieve another benefit besides the reduction of demand charges for customers.For instance, system operators can reduce the need of network reinforcement by sizing the infrastructure for a more flat profile coupled with ESS, instead of designing it for the highest power demand [3].Depending on the market conditions, other benefits can be achieved.The customers can take advantage of time of use energy price [4] by discharging the ESS when the energy price at the peak load periods is more expensive than the price during the off-peak periods.This can lead to additional electricity bill reduction [5].
Energy storage system technologies are used for a variety of applications [6,7].They can be classified in many different ways, according to the application area [8], based on the energy conversion [9], or depending on the quantity of energy that the ESS can provide [10].For "power-type" applications like peak shaving, the ESS have to maintain a constant delivery of power [11].
Although the improvements of battery energy storage system (BESS) efficiency and life cycle are increasing the interest for this type of storage [12], the high investment costs necessary for BESS solutions still raise concerns about the economic viability of this technology in power system applications [13][14][15][16].Therefore, an important aspect of the deployment of any BESS project is their proper power and energy sizing [17].If a BESS is not sized properly, it can generate negative results from an economic perspective.While small BESS may result in excessive aging-related depreciation cost, over-sizing systems may not attain optimal cost-benefit ratio due to their relatively high initial investment cost.
In response to the need to properly size BEES, several studies aiming to find the optimal sizing of BESS have been conducted [6].Recent work by Merei et al. [18] concentrates on commercial applications of BESS.The authors use sensitivity analysis to study the maximization of energy self-consumption via storage integration.The techno-economic analysis reveals that, for most commercial applications, BESS is not favored economically when battery degradation is taken into account.Other previous work by Magnor and Sauer [19] and Merei et al. [20] analyzed the optimal sizing of storage in the context of island grids and home storage systems.The authors propose a genetic algorithm-based method to model a non-linear set of equations including battery-aging characteristics.However, the solver results may not find a globally optimal solution to the described problem, and the studies do not provide design rules for future storage systems.
A sophisticated optimization method applied to find the best-suited battery storage system located in a residential suburban area has been described by Tant et al. [21].A multi-objective function is used to find the balance between voltage regulation, peak power reduction, and annual cost.A grid operator can use this method to support the decision of temporarily installing a BESS in problematic feeders to postpone grid upgrades in the short term due to work planning issues.By comparing the cost of grid upgrades, the grid operator may conclude that the BESS is also a valuable alternative in the long term.Recent work by Rahmann et al. [1] proposed an approach to determine the break-even points for different BESS considering a wide range of life cycles, efficiency, energy price, and power price.The results presented in this work show that depending on the values of round trip efficiency, life cycles, and power price, there are BESS technologies that are already profitable when considering only peak shaving applications.Although the authors model an optimization algorithm used for the sizing of the storage system, only the distribution company perspective is considered.
In contrast to the important contributions mentioned above, this work proposes a linear optimization method to define a cost-optimal sizing of the battery and power electronics for peak shaving application in industrial settings.In addition, this paper also presents a case study conducted with real industrial profiles, a techno-economic analysis evaluating the return on investment (ROI) of the system and battery degradation, and a linear programming (LP) approach allowing exact solution determination for BESS sizing.At the same time, the power flow optimization reveals the best storage operation considering energy purchase, peaks of consumption, and battery aging.
The remainder of this paper is organized as follows: Section 2 introduces the system layout, parameters necessary as input for the subsequent optimization procedure, and the aging model.The concepts of LP methodology, including equations and constraints for optimization, and a case study where different industrial profiles verify effectiveness of the proposed model are presented in Section 3. Finally, Section 4 provides the conclusion and outlines possible directions for future research.

System Layout
The energy management system proposed in this study is derived from measured and simulated data for an exemplary BESS.The simulations involve a grid-connected system shown in a schematic diagram in Figure 1a.The arrows in this figure illustrate the power flow direction for all component links.Additionally, Figure 1b illustrates all price components for industrial customers: the total energy consumption E total = ∑ load i where i denotes an averaged time segment of 15 min, the maximum power peak in the billing period P max , and the maximum power peak after peak shave P PS .Other variables necessary for subsequent modeling are explained in more detail later, along with the optimization problem definition.

Economic and Legal Framework for Industrial Customers
BESS are very flexible devices that can be used for many different applications [6].Depending on the application, several factors influence the attractiveness for energy storage systems.Particularly in behind-the-meter-scenarios (BTM), the economic attractiveness of energy storage systems depends not so much on the electricity price itself, but on the pricing structure [22].
Peak shaving is a typical BTM application that concentrates on the reduction of the peak demand of consumers.Peak shaving systems are only attractive in markets where demand charges amount to a proportionally large part of the electricity price.Already at a very early stage of electricity system development, system operators introduced electricity tariffs that included a demand based part added to the usage based part of electricity costs [23].This scheme has been established to provide an incentive for efficient grid usage.This is a so-called cost reflective tariff, since the level of demand is the main driver of network costs (i.e., grid reinforcement and transformer overloading) [24].
Battery storage is still a new technology associated with high perceived investment risk.This is likely the reason why most storage projects are currently conducted in well-developed countries [25].According to a study by Azure International, the most attractive countries for demand charge management in the commercial and industrial (C&I) sector are Australia, France, the USA (California), Japan and Germany.
Electricity costs are paid via the utility company selected by the consumer.The utility company keeps a small percentage for itself to cover generation and retail costs.Also, it transfers taxes, fees and surcharges to the relevant authorities and transfers network costs to the system operator responsible for the corresponding system.Therefore, the location of the network connection point defines network costs.For instance, the two eastern transmission system operators (TSOs) in Germany charge significantly higher prices than the two western TSOs, but prices also differ from one distribution system operator to the next inside the same regulation zone.Specifically, commercial and industrial customers who (typically) exceed 100.000 kWh energy consumption per year or 500 kW of average power demand pay an additional power price per kW to the energy price per kWh.The electricity price in the C&I sector typically has the following components:

•
Electricity generation (wholesale prices and retail costs); prices depend on negotiations between customer and utility company.

•
The network costs (transmission and distribution) are subdivided into two categories.First, power price per kilowatt, based on the maximum power peak in the billing period; this is the only power-specific price component; prices vary with connected voltage level, billing period, distribution system operator and duration factor.Second, energy price per kilowatt hour, based on the total energy consumption.

•
The total for standard rates including taxes, fees, surcharges (including renewable energy surcharge, electricity tax, CHP surcharge etc.).
These prices are based on a load profile considering a duration factor calculated as where ∆ feh is equal to the full load equivalent hours, ∑ i P grid-load i is the total energy consumption per year, and P peak is the yearly peak power at network connection point.Table 1 summarizes the costs of electricity for industrial customers in Germany.For customers with ∆ feh ≤ 2500 hours per year, the energy price of 0.18 e/kWh and power price of 12.78 e/kW are assumed.Customers with ∆ feh ≥ 2500 hours per year are charged an energy price of 0.13 e/kWh and a power price of 139.12 e/kW.This pricing scheme produces a dependence of cost versus duration factor as shown in Figure 2. The total cost decreases as the duration factor increases.A typical lithium-ion battery available on the market can provide up to 3 C (i.e., a 50 kWh battery can be discharged with 150 kW or in 1/3 h).As specific capacity costs are higher than specific power costs, load profiles with peaks below 1 hour are ideal for peak shaving with BESS.Typical loads producing steep peaks are power intensive plants and machinery with short start-up times or heat-up periods, like furnaces in the steel industry.Another precondition for the feasibility of peak shaving is periodic, predictable behavior of the load.Forecasting algorithms ensure that the storage system will be able to discharge its maximum energy when needed [27].Although such prediction tasks are indispensable for achieving the best BESS operation, they are outside the scope of this work.
A nonrepresentative study of nearly 300 industrial load profiles, conducted by Smart Power (https://www.smart-power.net) in 2017, showed that about 10% of all load profiles result in a static ROI of five years or less, and thus can be directly considered for peak shaving application (cf. Figure 3).Under the assumption that storage system prices will decrease by about 30% and demand rate will rise by about 30%, the number of loads interesting for peak shaving will rise to about 33% in the next few years [14].Interestingly, Schmidt et al. [28] construct a comparative study for promising electrical energy storage technologies.The authors also investigate how the derived rates of future cost reduction influence when storage becomes economically competitive in transport and residential applications.In terms of price per energy capacity, the technology that brings the most energy density to market is likely to become the most cost-competitive.For instance, lithium-ion batteries can be used in multiple applications and secure high-capacity markets such as battery packs for electric vehicles.
For the sole battery storage investment without an inverter, the following price structure is considered.
where C batt represents the total battery investment cost, C fix corresponds to the fixed cost including the housing of storage and all the peripheries, C var,batt denotes the energy specific cost of a storage system, and C opex,batt is the storage operation and maintenance (OPEX) cost within the battery lifetime.
As such, the overall cost C storage for the energy storage system can be expressed as: which includes battery storage with energy content E nom batt , and inverter with nominal power P nom inv .As container storage systems predominantly have battery racks and inverter units assembled to the same casing, no separate fixed costs for inverters are assumed, but are given as part of the overall storage fixed cost C fix .

Battery Aging Model
As described in our previous work [17], storage deterioration is a significant cost driver during energy storage operation.As a result, the aging of storage devices must be taken into consideration when simulating BESS operations.Lithium-ion batteries [7] suffer from continuous aging.For most batteries of this type, it is possible to separate the degradation into a pure time-dependent irreversible loss of battery capacity called calendric aging, and an energy throughput dependent cyclic aging [29][30][31][32].For calendric aging, the growth of the solid electrolyte interphase (SEI) is considered to be the major aging factor [33].The SEI protects the negative electrolyte from decomposition and corrosion and it is mainly formed during the first charging process [34].With time and usage of the battery, the SEI undergoes a structural conversion, reformation, and a slow growth of the SEI takes place [33].The cyclic aging can be attributed to either Lithium plating (particularly at low temperatures and high current rates), to exfoliation/particle cracking (particularly at higher currents and often at high SoC levels), or to irreversible structure changes (induced by frequent intercalation/de-intercalation of lithium ions) [31].All of the above effects could be summarized as loss of lithium inventory and/or capacity reduction at the anode/cathode side.The battery cyclic and calendric lifetime define the remaining state of health (SoH) until a certain capacity fade for a battery cell becomes evident.In this paper, we assume that the BESS must be replaced when SoH drops to 80% of the nominal capacity.The overall aging can be estimated using the superposition principle [29] aging tot ≈ aging cyc + aging cal . ( Value aging tot = 0 represents a new, unused battery, while aging tot = 1 corresponds to a situation when the remaining capacity of the battery is 80% of its original capacity.However, it is important to note that additional use of the storage system with aging tot > 1 may be allowed if the replacement of storage is set to below 80% of the SoH.A detailed analysis and validation of battery performance and aging models is provided in [35]. Accelerated aging tests performed at the Technical University Munich were used to build an equivalent circuit based aging model coupled with a thermic model [36] for a cell with graphite anode and nickel manganese cobalt (NMC) cathode.The key factors making this model appropriate are the analytic equations for estimating the cell degradation and the superposition of calendric and cyclic aging.Both are essential for an implementation of a linear model.The aging model and its linearization can be described as follows [37].
The analytical term for the calendar capacity fade C fade,cal as a function of battery state of charge SoC (%), temperature T ( • C), and time ∆ t , takes the following form C fade,cal (SoC, T, ∆ t ). ( Considering that a BESS enclosure can maintain the temperature constant, it is possible to fit a linear capacity fade for each time step.This is shown in Figure 4 where the piecewise line represents the calendric aging per time step dependent of the SoC [36], and the straight line shows its linearization.The accelerated aging test reveals that the aging rate increases faster at very low and very high values of SoC.This resulted in the kink at 90% SoC and small bends at the other test points.The linearization of the calendric aging for each timestep i can be expressed as follows: The NMC model also holds a detailed cycle aging model, but this can be neglected as very few peaks occur [37].For this reason, in this work, the cyclic aging is represented by the number of full equivalent cycles (FEC) that provides the overall energy throughput (counting either only charge or discharge direction) with any Depth of Discharge (DoD) per cycle divided by the available capacity battery storage [38].The number of FEC can be defined as The factor of 0.5 results from the conversion of charge throughput to full cycle counting.SoC denotes the state of charge, P batt the power flow via the battery, and E nom batt the nominal energy capacity of the battery.The maximum charge/discharge throughput is achieved at 80% SoH if there is no calendric aging.

Case Study
This section presents the application of the introduced model for dimensioning BESS for industrial peak shaving application.The industrial customer is responsible for buying, installing, maintaining, and operating the storage system.In this model, the energy used to charge the battery and the energy used for immediate consumption have the same cost, and both are considered in the industrial customer peak power calculation.As a result, the usage of a storage system is transparent from the point of view of the utility company.

Linear Optimization of BESSs
The economically optimal battery storage component sizing for an industrial customer equipped with storage system is obtained using LP.The load demand profiles considered in this study cover one full year, to capture all seasons with their characteristic.As the intent is to minimize the overall electricity cost, three types of costs are considered: the energy cost C energy_tot , the power cost C power_max , and the battery degradation cost C storage_deg .The energy cost is composed of the base energy price, fees, taxes, and stock exchange price.The power cost is charged by the network operator on the basis of the duration factor.The battery degradation cost, also called aging cost, is the major cost driver during storage operation, caused by cyclic and calendric aging.
The annual cost flow analysis presented here takes into account the discounted storage cost caused by degradation.As such, this simulation allows to estimate the profitability of a BESS for a full life of the battery.
All variables and parameters considered in this study are described in Table 2. To meet the electrical demand, P load i , the system attempts to use power from the battery, P batt-load i , or it draws power from the grid, P grid-load i , i.e., In the same way, the power imported from the grid (P grid−load i + P grid−batt i ) in each time step i is restricted to the maximum power for the period.The two constraints can be represented as where P peak−shave j represents the highest point of demand in the billing period j.For instance, considering only the highest load of the year, all data points i should be limited to the same maximum annual limit P peak−shave j .However, if we consider the seasonal billing period where there are two independent thresholds, each season is limited to its own limit.The peak power is used to calculate the optimal solution power cost.The bidirectional power flow from the storage inverter to the battery is stored in an auxiliary variable, P batt i , and correlated with the inverter efficiency, η inv , as follows where η inv is the average one-way efficiency of the inverter.The reciprocal efficiencies are the battery charge power P grid−batt i and the discharge power, P batt−load i , both limited by the nominal power flow from the inverter to the battery 0 ≤P grid−batt i ≤ P nom inv , 0 ≤P batt−load i ≤ P nom inv , (11) where P nom inv corresponds to the inverter size.The battery energy content at time step i (E batt i ) satisfies the recurrence relation where SD batt represents the self-discharge factor of the battery and d = 96 the conversion factor of time steps per day.The energy content of the storage system is furthermore confined by an upper boundary, that decreases upon usage and aging according to the SoH.The SoH is defined as the irreversible capacity fade over time, related to the nominal battery capacity, and E batt i is a fraction of the total energy content of the battery installed: The state of health of the storage system at time step i also satisfies the recurrence relation Using Equations ( 6) and ( 7), the calendric and cyclic aging can be estimated as and The calendric aging is affected by the storage temperature and its SoC level according to Swierczynski et al. [32].Despite the fact that the charge/discharge process leads to dissipative heat generation and unavoidable temperature changes within the battery, the very low utilization ratio of the storage system and the restriction to a maximum C-rate of 3 limits the effects of temperature variations significantly.
As a result, the additional cyclic aging degradation of time step i is estimated by the energy throughput in time step i (P batt j • ∆t res ) divided by the energy content of the system E batt i and is normalized with the factor of 0.5 and the technology specific cycle life indicator Life 80% Cyc .Similarly, the SoC can be expressed as The inverter nominal power is limited to three times the battery nominal capacity The optimal solution must satisfy all constraints described above.It aims to reduce the overall cost by minimizing the expenses for energy purchase and implicit cost caused by battery degradation.This cost model is divided into three components, i.e., minimize The first component C energy_tot comprises the cost of energy purchased from the grid, while the second component C power_max is the peak induced cost based on the highest point of demand (or peak) within billing period (monthly or annually).These two components are evaluated as follows: where C buy and C power are the retail electricity price and the peak-power tariff, respectively.The third component estimates the storage system degradation cost that can be represented as where ∆t denotes the time span covered with the simulation (here one year) and ∆SoH the total battery aging.The full battery related cost is then calculated in consideration of the initial installation investment cost.

Case Description
Four industrial load profiles (A-D) shown in Figure 5 are used to verify the effectiveness of the proposed model.Data used for simulations was adapted from real measurements and averaged with a resolution of ∆t res = 15 min [39].This time discretization results from the fact that, in the model region, the 15 minutes demand average is registered and its maximum value is used for tariff calculation over a period of one month or one year [40].It is assumed that temperature is kept stable at approximately 25 • C. Parameters and price information for the BESS/inverter system used in the simulations are listed in Table 3.As described in Section 2, the electricity cost has two main components: the total energy consumption, and the power peak cost in the billing period.According to German law StromNEV §19I, every grid operator is obligated to offer a monthly billing scheme, i.e., instead of the highest load of the year, the highest load of the month is the basis for the price per kW.Hence, this study considers both yearly and monthly billing schemes.The results presented in the next section describe not only the optimal BESS/inverter component sizes, but also the optimal billing scheme.

Effect of Sizing Considering BESS Degradation Cost
The objective function and the constraints structured in this have linear relationships.This means that the effect of changing a decision variable is proportional to its magnitude.For this reason, the economically optimal battery storage component sizing for peak shaving is obtained using LP.The linear optimization was implemented in MATLAB (MathWorks, Natick, MA, USA) code using a dual-simplex algorithm, which is based on a conventional simplex algorithm on the dual problem [41].Each one-year simulation considered 15-min time resolution, co-optimized the storage and inverter size, and took on average 700 s on a workstation with Intel Core i5 processor at 3.5 Ghz and 16 GB of memory.
The optimal storage and inverter size for each profile (A-D), as well as a number relevant technical and economical indicators, are presented in Table 4.The investment comprises the overall cost C storage described in Equation ( 3).The operation cost, C opex,batt , reflects the German market and is calculated as 0.6% of the investment plus 6 e/kW.Table 5 shows the OPEX components considered in this paper.The total return equal to the internal return rate (IRR) [42], is calculated without inflation or price changes based on the total savings of the first year.Likewise, total savings and amortization time are static calculations 'Profile A' has an annual load of 9350 MWh and features weekday peaks and small load during weekends.As shown in Table 4, this profile exhibits similar results for yearly and monthly billing scheme.Figures 6 and 7 illustrate the results obtained for the two billing schemes.The yearly billing scheme requires an initial investment of almost e20.000less compared to the cost of the system optimized for monthly billing, and it can generate additional 5% of total return in a shorter time.Although monthly billing scheme with a 51 kWh battery and 152 kW inverter (Figure 7) can increase the peak load capping, it also shortens the battery end of life by seven years (Figure 8).Therefore, all things considered, the yearly billing scheme is more suitable for 'Profile A' because it delays the battery replacement and provides several extra years of saving grid charges before it becomes necessary to invest in a new battery system.Further considering the monthly billing scheme, the optimizer determines the optimal storage size of 51 kWh and inverter nominal power of 152 kW. Figure 9 illustrates the power flows for a three-day period during the first week of May.The left panel shows the load consumption, the power flow imported from the grid for direct use or to charge the battery, as well as the maximum power peak after shave.The right panel shows the periodically changing charge level of the storage system (SoC), and the evolution of battery degradation (SoH).It clearly shows that the capacity fade is stronger when the energy throughput is high.Figure 10 shows the battery power profile and the same capacity fade in the terms of C-rate.In contrast to 'Profile A', results for 'Profile B' show that the yearly billing scheme is not suitable for profiles with low average load and relatively high peaks.Although it is possible to reduce 30% of the peak load using a 57 kWh BESS with 171 kW inverter, this system configuration provides no savings to support the initial investment.In this case, the monthly scheme is more profitable, resulting in a peak reduction of 13% with a 21 kWh battery and 63 kW inverter.
Similarly, 'Profile C' has a negative IRR when considering the yearly billing scheme.This is caused by the seasonal nature of the load.As can be seen in Figure 5, the consumption in the last three months of the year is very high compared to the rest of the year.Analyzing the total return value in Table 4, the monthly billing scheme appears to be the right choice for this profile.However, a peak load reduction of only 1% is too slow to justify the installation of a BESS.
Finally, 'Profile D' presents the most extreme case.Considering the exposed investment cost for BESS and price schemes, there is no advantage to installing a BESS for peak shaving purpose for this profile.Figure 11 illustrates the optimal annual peak shaving limit for the profiles B, C, and D, as well as the state of charge and state of health for the storage system used in each case.Similarly, Figure 12 shows the optimal monthly peak shaving limit, the SoC, and SoH for the same three profiles.The Appendix A provides the battery power profiles for all investigated scenarios.
It can be seen that the optimization process minimizes expenses using the capacity of the storage system to decrease the peak power.The optimal power flow shows that the battery cycles are short, meaning that the battery is charged to the maximum necessary level just before being drained.This occurs due to the presence of SoC-dependent calendric degradation as one of the optimization criteria.At the same, cyclic aging is not a determinant in peak shaving applications because the BESS has only a low number of charging/discharging cycles and energy is never stored in the battery for a long time.For this reason, calendric degradation is the most important cost driver in storage systems for peak shaving applications.
To analyze the relation between load size and return of investment, consider Table 6.Optimization runs were performed scaling the load size of Profile A from 10 to 40,000 MWh/a.The relation between the peaks and the loads were kept the same as in the original profile, resulting in exactly the same shape of battery SoC profile.As an overall trend, customers with large loads require BESS with large storage size and large nominal power of the inverter.Loads smaller than 1000 MWh/a have a negative IRR and an extensive payback period, rendering them unsuitable for BESS-based peak shaving applications.On the other hand, larger load profiles have a substantial improvement in the payback period.The results show that the BESS can be used for almost 18 years before reaching end of life at 80% of SoH.Although the peak capping is the same in all simulations, the battery usage differs for each load size because of the assumed discrete sizing of BESS and inverters in steps of 10 kWh and 10 kW respectively.It is clear that larger load sizes can benefit from BESS-based peak shaving with better economical results.In addition, it is interesting to analyze the impact of the peak capping variation on the payback period as well as the battery life time.Refer to Figure 13 and Table 7 for a detailed comparison.To generate these results, an additional constraint was added to the linear model described in Section 3.1 P peak−shave = max(P load ) • (1 − PC); PC = [0.01,0.025, 0.05, 0.10, 0.15, 0.20, 0.25, 0.30], (25) where PC is the fixed amount of the peak that must be shaved.As shown in Table 7, all scenarios are profitable.However, the best IRR is obtained for peak capping equal to 5% of the total load.Smaller peak capping values extend the battery lifetime, but also extend the payback time as less savings of peak power reduction may be attained.In contrast, larger peak capping values increase the payback, but shorten the battery life.As an overall trend, the increase of inverter size has a direct relation to the peak shaved load, i.e., 2.5% of load shaving needs a 60 kW inverter, and 25% load shaving requires ten times more.It is a straightforward relation because the inverter is sized according to the load power peak.Interestingly, the battery sizing does not follow the same trend because it is related to the number of peaks that must be shaved.In contrast to the optimal result where the battery is charged closer to the load peaks, larger peak load capping results in a smaller C-rate because the battery is charged slowly to avoid violating the maximum power peak allowed, and the battery must keep energy content for a longer period.

Conclusions and Future Work
This article describes a linear optimization model to size the most cost-effective BESS for a variety of industrial load profiles and multiple billing schemes.The optimization approach formulated in this work minimizes the storage degradation cost and the maximum power peak in the billing period.
The optimal BESS size and the number of capped peaks are directly related to the load profile.As an overall trend, for the exemplary load profiles under investigation, the monthly billing scheme is more attractive for industrial customers because of the number of peaks that can be capped with acceptable BESS sizes.For instance, a 51 kWh/152 kW BESS can shave 243 peaks which represents 6% of the maximum load peak and results in a 15,223 e of annual savings.As a general remark, considering the current cost of storage and retail energy tariff valid in Germany for 2016, most scenarios favor storage system installation.The expected increases of electricity prices and the reduction of BESS costs are likely to accelerate this trend.Although this work uses parameters corresponding to German market conditions and regulations, the described methodology can be easily adapted to other jurisdictions that use or consider peak power penalty, albeit with a different billing period scheme and retail electricity tariff models.
This work is limited to the optimization of storage systems using historical data on specific industrial load profile.Load forecast is the key to developing online energy management controllers for BESS.Such a forecasting task is outside the scope of this paper, but it will be considered in the future.Future work will also examine the possibility of decreasing the idle time of the battery system by sharing the same BESS between multiple industrial customers.

Figure 1 .
Figure 1.System configuration showing all (a) considered power flows and (b) customer load curve with price components.

Figure 3 .
Figure 3. Static return on invest (ROI) of peak shaving storage systems in years based on 288 industrial load profiles analyzed by Smart Power in 2017 (blue), and the static ROI projection where the investment is reduced by 30% and the energy rate is raised by 30% (yellow).

Figure 6 .
Figure 6.Industrial load profile A with yearly billing scheme (left), and battery state of charge and state of health (right).

Figure 7 .
Figure 7. Industrial load profile A with monthly billing scheme (left), and battery state of charge and state of health (right).

Figure 8 .
Figure 8. Extrapolation of the one-year results to represent the degradation after 10 years of usage for Profile A.

Figure 9 .
Figure 9. Power flow analysis for a three-day period: load and power flows within the system (left); time correlated evolution of battery state of charge (SoC) and resulting state of health (SoH) decline (right).

Figure 10 .
Figure 10.Battery power profile analysis for a three-day period and resulting state of health decline.

Figure 13 .
Figure 13.Profile A-Impacts of peak capping variation.Initial investment Amortization time (ROI), the number of years before the battery end of life (EOL), and the number of years the BESS will keep being used and generating savings through peak shaving after ROI being achieved.

Full Load Equivalent Hours (∆ feh ) <2500 h/a >= 2500 h/a
Network cost vs.duration factor at constant total energy consumption of 100 MWh in a given network area, medium voltage tariff, 2017.

Table 2 .
Variables and parameters used for the battery modeling and optimization routines.

Table 4 .
Economical and technical comparison of system optimization results.

Table 6 .
Profile A with yearly billing scheme and duration factor of 4431 h/a.Peak load capping of 5% which represents 20 capped peaks per annum.
Figure 12.Industrial load profile B, C, and D with monthly billing scheme (left), and battery state of charge and state of health (right).

214.27 MWh/a ; Storage system: 57 kWh 171 kW
Battery power profile with yearly billing scheme and resulting state of health decline.