Next Article in Journal
Optimizing Water-Light Complementary Systems for the Complex Terrain of the Southwestern China Plateau Region: A Two-Layer Model Approach
Previous Article in Journal
Dynamic Shear Responses of Combined Contaminated Soil Treated with Nano Zero-Valent Iron (nZVI) under Controlled Moisture
Previous Article in Special Issue
Electric Vehicle Load Estimation at Home and Workplace in Saudi Arabia for Grid Planners and Policy Makers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Probabilistic Planning for an Energy Storage System Considering the Uncertainties in Smart Distribution Networks

by
Ahmed A. Alguhi
1,2,*,
Majed A. Alotaibi
1,2 and
Essam A. Al-Ammar
1,2
1
Electrical Engineering Department, Faculty of Engineering, King Saud University, Riyadh 11421, Saudi Arabia
2
K.A.CARE Energy Research and Innovation Center in Riyadh, King Saud University, Riyadh 11421, Saudi Arabia
*
Author to whom correspondence should be addressed.
Sustainability 2024, 16(1), 290; https://doi.org/10.3390/su16010290
Submission received: 18 October 2023 / Revised: 27 November 2023 / Accepted: 11 December 2023 / Published: 28 December 2023
(This article belongs to the Special Issue Smart Grid Technologies and Renewable Energy Applications)

Abstract

:
Today, many countries are focused on smart grids due to their positive effects on all sectors of a power system, including those of operators, utilities, and consumers. Furthermore, the usage of renewable energy sources for power production is quickly expanding due to the depletion of fossil fuels and the emissions caused by their use. Additionally, intermittent power generation from renewable energy sources, such as wind and solar, necessitates the use of energy storage devices with which to ensure a continuous power supply to meet demand. This can be accomplished by employing an appropriate storage device with a sufficient storage capacity, thus enabling a grid-connected solar PV and wind system to have enhanced performance and to reduce adverse effects on the power quality of the grid. In this study, a probabilistic planning model that takes the intermittent natures of solar irradiances, wind speeds, and system demands into account is introduced. A novel criterion is also adopted to map the three-dimensional spaces of intermittency with the proposed model for optimizing BESS charging/discharging decisions. This planning model is intended to minimize the economic costs of investment and operation of a battery energy storage system (BESS) for a planning period. Moreover, the substation and feeder upgrade costs, as well as the overall system loss costs, are included in the proposed model. Particle swarm optimization (PSO) is utilized to find the optimal sizing, location, and operation of energy storage systems. The proposed methodology was validated using a 69-bus distribution system.

1. Introduction

Electricity distribution networks hold a pivotal position within power delivery systems, which is primarily because of their proximity to end users. Over the last few years, distribution networks have witnessed the integration of numerous innovative technologies in response to various technical and economic factors. Among these advancements, energy storage systems (ESSs) have emerged as a technology poised to assume a vital role in the energy sector’s imminent future. ESSs are anticipated to deliver electricity with the utmost cost efficiency and the requisite level of quality [1].
ESSs come in various forms, from traditional batteries to cutting-edge technologies, such as pumped hydro energy storage and thermal energy storage. The planning, deployment, and management of these systems are at the heart of the energy transition, enabling us to harness the full potential of renewable resources and enhance grid resilience.
The challenges of ESS planning are formidable but surmountable. Energy storage planning requires interdisciplinary collaboration among engineers, policymakers, economists, and environmentalists. It demands a holistic approach that considers not only technical feasibility but also the socioeconomic and environmental impacts of energy storage projects. It necessitates innovation in policy frameworks, market structures, and regulatory mechanisms to incentivize investments in energy storage infrastructure.
Saudi Arabia has taken significant measures to incorporate renewable energy sources—alongside its traditional reliance on oil and gas—within the national energy blend. Saudi Vision 2030 lays out a strategic plan for the nation to reduce its dependence on oil in the energy mix, acknowledging that Saudi Arabia is currently progressing in the development of a competitive renewable energy sector. One of the main forces behind sustainability is the development of renewable energy projects, which will help reduce emissions and replace high-value fuel in power production. Indeed, energy storage can help address the problems with renewable energy, such as intermittent solar and wind output power levels.
ESSs possess the ability to swiftly address substantial shifts in energy demand, enhancing grid responsiveness, minimizing the necessity for constructing backup power facilities, and effectively mitigating the intermittency challenges associated with solar and wind energy generation. However, incorporating ESSs in a distribution system without optimizing their size, location, and operation mechanisms will affect the stability and dependability of the power system. Therefore, optimizing the allocation and capacity of ESSs within a network is essential for enhancing their performance and improving the overall performance and quality of a power system.
The selection of an ESS is dependent on the application of its use. For example, suppose that it is used to supply power during transitions between power sources. In that case, one should select either capacitor storage or an SMES; however, when it is used for power quality applications, one should choose an ESS such as a PHS, flywheel, SMES, CAES, or capacitor. Further, if energy management is the desired application, a battery energy storage system (BESS) is the most proper choice of ESS; therefore, in this study, a BESS was chosen as the candidate ESS in the application under study since its rated power and discharge time were suitable for the application of interest [2].
Based on a survey of the literature, the following referenced articles contributed to the understanding and advancement of BESSs from various perspectives. Ref. [3] presented a comprehensive review encompassing BESS technologies, optimization objectives, constraints, approaches, and outstanding issues, providing a broad overview of the current landscape. Ref. [4] focused on the specific context of New York State while investigating the impacts of BESS technologies and the integration of renewable energy on the energy transition, offering insights into regional dynamics. Ref. [5] presented a modeling approach for a large-scale BESS and emphasized its application to power grid analysis. Ref. [6] contributed to the integration of BESSs into multi-megawatt grid-connected photovoltaic systems and addressed the complexities of large-scale renewable integration. Ref. [7] offered a methodological framework for site selection and capacity setting for BESSs in distribution networks with renewable sources, contributing practical strategies for optimal deployment. Additionally, there has been much research providing comprehensive overviews of ESSs in distribution networks with a focus on placement, sizing, operation, and power quality to enhance the performance and reliability of distribution networks. Ref. [8] offered an extensive overview that covered various aspects of ESSs in distribution networks. Ref. [9] contributed insights into the optimal location, selection, and operation of BESSs and renewable distributed generation in medium–low-voltage distribution networks. Ref. [10] presented an approach based on genetic algorithms for the integration of energy storage systems in AC distribution networks that addressed their optimal location, selection, and operation. Ref. [11] focused on low-voltage residential networks in the UK while addressing the optimal placement, sizing, and dispatch of multiple BESSs. Ref. [12] introduced a PSO algorithm for optimizing energy storage capacities in distribution networks while considering probability correlations between wind farms. Ref. [13] investigated the installation of battery energy storage systems with renewable energy resources in distribution systems while considering various load models. Ref. [14] discussed energy storage optimization in the configuration of active distribution networks using distributed approaches. Ref. [15] explored the optimal allocation and operation of an energy storage system with high-penetration grid-connected photovoltaic systems, contributing to the sustainable integration of renewable energy in distribution networks. Ref. [16] introduced the Whale Optimization Algorithm for optimizing the placement and sizing of BESSs to reduce losses, contributing to efficiency improvements in distribution networks. Ref. [17] focused on energy storage system scheduling for peak demand reduction by employing evolutionary combinatorial optimization techniques to enhance sustainable energy technologies. Ref. [18] addressed distributed generation and energy storage system planning for a distribution system operator and highlighted the importance of integrated planning for improving system efficiency. Refs. [19,20] contributed to optimal ESS allocation for benefit maximization and load shedding to improve the reliability of distribution systems and employed advanced optimization techniques, respectively. Ref. [21] provided an economic analysis model for an ESS applied to a distribution substation and offered insights into the economic feasibility of ESS integration. Ref. [22] conducted an analysis of the adequacy and economics of distribution systems integrated with electric energy storage and renewable energy resources, and they emphasized the importance of considering both adequacy and economic aspects in system planning.
Even though numerous studies have explored various optimization methods for energy storage technologies (ESTs), there is a noticeable scarcity of comprehensive information and up-to-date data pertaining to the comprehensive planning and application potential of ESSs for the integration of renewable energy sources (RESs) into distribution planning. In addition, in response to the escalating energy demand, the heightened integration of renewable sources, and the recent evolution of grid demand, there is a critical need for comprehensive energy storage planning strategies. Additionally, there is a pressing need for further research that assesses the impacts of ESSs on distribution systems while considering both technical and economic constraints. Thus, this study aims to bridge this knowledge gap by investigating the optimization of energy storage within distribution systems that incorporate renewable energy sources, such as DGs, and a novel criterion for mapping the three-dimensional spaces of intermittency with the proposed model is adopted to optimize BESS charging/discharging decisions. The objective is to minimize the total planning and operational costs while considering technical and economic constraints and, ultimately, identifying the most cost-effective solutions for the placement and sizing of energy storage units within a smart distribution network.
The main contributions of this study are the following.
  • A robust probabilistic planning model for BESSs in distribution networks is developed in order to optimize the location, sizing, and operation of BESSs and to determine the most economical BESS with the lowest overall planning cost while taking uncertainties in wind speed, solar irradiance, and system demand into account and considering technical and economic factors.
  • The proposed model brings a significant advantage for distribution companies, as it strategically addresses challenges such as energy losses, deferral of system upgrades, and effective energy management during off-peak and on-peak periods. By maximizing benefits through these considerations, the model contributes to the overall efficiency and sustainability of distribution networks.
  • In contrast to previous approaches that focused solely on the demand state for charging and discharging BESSs, this study introduces a novel paradigm. It employs three-dimensional spaces encompassing the wind state, PV state, and demand state to optimize charging and discharge decisions. By evaluating the interplay of these dimensions, the optimization model makes informed decisions regarding energy production or absorption, ultimately minimizing the planning and operational costs for the entire system.
  • This research makes a valuable contribution to the electricity sector by providing an opportunity for electricity utilities to leverage the study’s findings to advance the transition toward renewable energy sources.
  • This study’s findings can be leveraged to advance the transition toward sustainability in the energy sector.

2. Modeling of Uncertainties

The uncertainty in smart distribution networks is a critical area of research and development in the field of electrical power systems. As the integration of renewable energy sources, advanced monitoring and control technologies, and the increasing complexity of grid operation become more prevalent, it is essential to have a robust theoretical foundation for dealing with uncertainty. In this context, uncertainty refers to the variations and unpredictability in generation, consumption, and grid conditions that can affect the safe and efficient operation of an electrical grid.
The most significant theoretical enhancements and considerations related to the description of uncertainty in a smart distribution network are as follows.
  • Stochastic models: Traditional deterministic models of power system operation are no longer sufficient to capture the complexity of modern distribution networks. Stochastic models, which incorporate probability distributions for various parameters, are gaining importance. These models allow for a more realistic representation of uncertainties, such as the intermittent nature of renewable energy generation or the variability in electricity demand.
  • Incorporating renewable energy sources: With the increasing penetration of renewable energy sources, such as wind and solar, uncertainty in generation becomes a significant concern. Theoretical improvements involve better modeling of the variability and intermittency of renewable generation, as well as their spatial and temporal correlations. This can be achieved through advanced statistical methods, such as time-series analysis and probabilistic forecasting techniques.
  • Demand uncertainty: Electricity demand can also be uncertain due to factors such as weather conditions, economic fluctuations, and changes in consumer behavior. Theoretical improvements should focus on capturing this demand uncertainty, including short-term load forecasting and the development of probabilistic demand models.
Therefore, the location of renewable resources in a distribution system requires effective modeling of demand and DGs. In some research, such as that in Ref. [13], the authors used deterministic methods that represented the demand and DGs as constant power based on their average or maximum power; although these methods can be easily executed, they may lead to unrealistic and inaccurate results. Therefore, much research considered the stochastic nature of demand and renewable resources to find more accurate and realistic results. This stochastic nature can be taken into account by using either time-series modeling or probabilistic modeling, which are the two main methods for modeling the nature of such systems [18,19,20,23,24,25]. However, in a recent lecture [26], the methodologies employed for modeling the nature of systems were thoroughly examined in order to move beyond conventional time-series and probabilistic approaches. The lecture spotlighted alternative techniques, such as robust optimization, information gap decision theory (IGDT), and interval approaches, among others, offering a mosaic of tools for researchers. The selection of a specific method depends on the nature of the system, the type of data available, and the goals of the analysis.
Much research, such as that in Ref. [18], has used the time series to model the stochastic nature of system components; this approach proposes a forecasting horizon of one year ahead for both demand and renewable resources. Then, it uses the forecasted shape to determine an ESS’s charging and discharging power at each hour. This is hard to implement for a year-long planning model, and it requires a very long time for the simulation for vast sets of historical data. This may lead to unrealistic and inaccurate results. Therefore, probabilistic methods are effective for year-long planning models, as they can handle the stochastic nature of demand, wind speed, and solar irradiance, as described in [19,20,23,24,25]. This is because they can be easily implemented and require less time for simulations with vast sets of historical data in comparison with time-series models. A probabilistic model was applied in this study to model the stochastic nature of wind speed, solar irradiance, and demand fluctuation. In this method, historical data on the wind speed, solar irradiance, and demand fluctuation were described by using a specific probability distribution function (PDF) (1). Then, this PDF was divided into several states, which were used to generate a probabilistic model for the load, wind speed, and solar irradiance.
The wind speed, solar irradiance, and load are assumed to be independent in this work.
p ( x ) = P [ a < x b ] = a b f ( x ) . d x
These were modeled as described in the following.

2.1. Modeling of the Wind Speed and Wind Turbine Output Power

In this study, we opted for the Weibull distribution to capture the intermittent nature of wind speed data, a choice that is commonly employed in numerous studies [23,24]. We developed a probabilistic wind speed model with a step size of 1.1 m/s. For any given wind speed dataset, it is crucial to ascertain the mean ( μ ) and the standard deviation ( σ ) of the data. Additionally, the wind speed characteristics at any location are characterized by the Weibull scale parameter (c) and the Weibull shape parameter (k). The key parameters of a wind turbine include the cut-in speed ( v i n ), rated speed ( v r ), cut-out speed ( v o u t ), and rated or maximum output power ( P w r ). Equation (2) defines the Weibull distribution function, while Equation (3) is used to compute the Weibull distribution parameters [24].
f ( v ) = k c v c k 1 e x p v c k
k = σ μ 1.086 , c = μ Γ 1 + 1 k
As can be seen in Figure 1, the characteristic curve of wind power can be categorized into three distinct regions: no wind power, de-rated power, and rated power. The area in which there is no wind power is defined as the area in which the wind speed is less than the turbine’s cut-off speed. The turbine’s blades are unable to overcome the friction generated prior to a speed reduction because there is not enough torque. In the de-rated speed region, the output power of the turbine dramatically increases to achieve the rated power of the turbine as the wind speed rises over the cut-off speed. In the rated power region, the wind speed exceeds the turbine’s cut-out speed, and the wind turbine will shut down as a preventative measure to protect the rotor from the strong forces acting on the turbine’s structure. A wind turbine’s output power can be obtained from (4) [23,24]. Table 1 provides the data parameters for the wind turbine and required wind speed in this study, while Table 2 presents the multistate probability model of a wind-based DG.
P w v = 0 ,           0 < v < v i n    a n d   v > v o u t P w r v v i n v r v i n ,              v i n < v < v r P w r                                  v r < v < v r  

2.2. Solar Irradiance and PV Output Power Models

Similar to when modeling the wind speed, there are various distributions available for studying the uncertainty associated with solar irradiance based on its statistical properties. The probability density function that mirrors the probabilistic behavior of solar irradiance (s) can be represented by using the Beta probability density function f s , as outlined in Equation (5). The calculation of the Beta distribution parameters β and  α is demonstrated in Equation (6) [23,24]. Table 3 provides the data parameters of the solar system needed in this study.
f s = Γ α + β Γ ( α ) Γ ( β ) s α 1 1 s β 1   ,   0 s 1   0    o . w
β = ( 1 μ ) σ 2 μ ( 1 μ ) σ 2     and     α = β × μ ( 1 μ )
The level of PV generation power is intricately linked to the solar irradiance. Consequently, the PV output power can be obtained by using Equation (7) [27]. Table 4 presents the multi-state probability model of a PV-based DG.
P P V s = P P V r s 2 s S T D × r C   f o r   0 < s < r C P P V r s s S T D   f o r   r C s < s S T D P P V r   f o r   s > s S T D

2.3. Demand Models

The normal distribution function ( f ( d ) ), as outlined in Equation (8), was employed to serve as a probabilistic model for representing the system demand behavior, a choice aligned with prior studies [23,24]. The analysis encompassed 8760 hours of load data points—specifically, those of the IEEE’s Reliability Test System (RTS)—with a mean of μ = 0.6142 and a standard deviation of σ = 0.1448. Table 5 illustrates the multi-state probability model for demand.
f ( d ) = 1 σ 2 π e 1 2 d μ σ 2

2.4. Generating Operating Scenarios for the Overall System

After establishing the probability distribution functions for the wind speed, solar irradiance, and system demand, the next step involved partitioning these PDFs into numerous discrete states for their integration into the calculations. The process of selecting these states carried significant importance, as it entailed striking a delicate balance between the precision of the outcomes and the intricacy of the analysis. The division of the PDFs into multiple equidistant intervals depended on factors such as the maximum value and the number of intervals necessary, as shown in Ref. [25].
Once all of the states for the wind power, solar power, and system load were defined, a three-column matrix encompassing every conceivable combination of these states—or scenario—was constructed. In this matrix, the first column signified the various levels of load states, the second column represented the states of the solar DG output power, and the third column signified the wind-based DG output, as shown in Table 6. This multi-scenario matrix was structured with a number of rows equal to the total number of scenarios, which was determined by multiplying the numbers of wind states, solar states, and load states. The probability associated with each scenario was calculated as the product of the probabilities of the wind state, solar state, and load state for that specific scenario, assuming that the wind speed, solar irradiance, and load were independent events.

3. Problem Formulation and Planning Model

This section outlines the formulation of the problem of the proposed probabilistic planning model for a BESS and explains how BESSs can be integrated into power distribution networks alongside renewable energy resources. The probabilistic planning model encompassed the modeling of probabilistic PV and wind power outputs, as well as the probabilistic representation of system demand. Within this model, various factors were considered, including the costs associated with investment in and operation of BESSs, substation and transmission line upgrades, and power losses. The primary objective of this planning was to minimize the cumulative investment and operational costs over the planning horizon, and this optimization was achieved by determining the optimal placement, sizing, and operation of BESSs.

3.1. Objective Function

The objective function is composed of BESS investment costs ( B C i ) and economic costs of operation ( B C e o p ) , along with substation and feeder costs ( S C s u ) and power loss costs ( L C L ). Equation (9) presents the mathematical representation of this objective function.
Min   c o s t = m i n ( N P V ( B C i n + B C e o p + S C s u + L C L i n e ) )
The mathematical formulations for each term in the objective function are displayed in Equations (10) through (13). The net present value (NPV) can be computed either from the future value (FV), as expressed in Equation (14), or from the annual value (AV), as described in Equation (15).
B C i n = i = 1 N ( C p × P B E S S , N r a t e d + C E × E B E S S , N r a t e d )
B C e o p y r = ξ ( ( 8760 × p ( ξ c h ) × C p r i c e _ o f f × P y r , ξ c h c h )                                             ( 8760 × p ( ξ d i s ) × C p r i c e _ o n × P y r , ξ d i s d i s ) ) , y r
L C L i n e y r = ξ ( ( C p r i c e _ o f f × 8760 × p ( ξ c h ) × i j = 1 N l i n e P i j ,    y r , ξ c h , l o s s )                                             + ( C p r i c e _ o n × 8760 × p ( ξ d i s ) × i j = 1 N l i n e P i j ,    y r , ξ d i s , l o s s ) ) , y r
S C s u y r = T r a n s m i s s i o n   c o s t + s u b s t a t i o n   c o s t
N P V = F V ( 1 + I R ) n
N P V = A V ( 1 + I R ) n 1 I R ( 1 + I R ) n
where C p is the BESS’s power cost in $/kW, C E is the BESS’s energy cost in $/kWh, P B E S S , N r a t e d is the BESS’s rated power at bus number N, E B E S S , N r a t e d is the BESS’s rated energy at bus number N, ξ c h   a n d   ξ d i s , are the charging and discharging states, p ( ξ ) is the probability of each state ξ , C p r i c e _ o f f is the energy price at off-peak hours, C p r i c e _ o n is the energy price during peak hours, and P y r , ξ c h c h and P y r , ξ d i s d i s are the charging and discharging of power at state ξ in year y, respectively. P i j , y r , ξ c h , l o s s is the power loss of line ij at state ξ , P i j , y r , ξ d i s , l o s s is the power loss of line ij at state ξ in year yr, IR is the interest rate, and n represents the number of years.

3.2. Model Constraints

3.2.1. Network Constraints

Equations (16) and (17) define the power-flow equation during the charging and discharging of the BESS at each bus i and state ξ . Equation (18) is used to calculate the reactive power at each bus i and state ξ .
P G i , ξ , y r + C w ξ × P w r i + C P V ξ × P P V r i P D   i , ξ , y r η c h × P B E S S ξ , y r c h = j = 1 N V i , ξ , y r × V j , ξ , y r × Y i j × c o s ( θ i j + δ j , ξ , y r δ i , ξ , y r ) , i , ξ ξ c h ,   y r
P G i , ξ , y r + C w ξ × P w r i + C P V ξ × P P V r i P D   i , ξ , y r η d i ξ × P B E S S ξ , y r d i s = j = 1 N V i , ξ , y r × V j , ξ , y r × Y i j × c o s ( θ i j + δ j , ξ , y r δ i , ξ , y r ) , i , ξ ξ d i s ,   y r
Q G i , ξ , y r Q D   i , ξ , y r = j = 1 N V i , ξ , y r × V j , ξ , y r × Y i j × s i n θ i j + δ j , ξ , y r δ i , ξ , y r ,   i , ξ , y r
where P G i , , ξ , y r is the generated power at bus i in state ξ during year yr, P D i , ξ , y r is the demand power at bus i in state ξ during year yr, P w r is the rated power of the wind turbine, C w ξ is the percentage of output power of wind in each state, P P V r is the rated power of the PV unit, C P V ξ is the percentage of the PV output power in each state, η c h and η d i ξ are the charging and discharging efficiency of the BESS, ( V i , ξ , y r ) is the defined voltage at bus i in state ξ in year yr, ( Y i j ) defines the magnitude of admittances of the line between buses i and j, θ i j is the angle of admittances, and δ is the voltage angle.
Equation (19) is used to calculate the total power loss in the line in state ξ during year yr, and Equation (20) defines the line current between bus i and bus j in state ξ during year yr. Equations (21) and (22) define the per-unit limit of voltage ( V i , ξ y r ) at bus i in state ξ during year yr and the slake voltage ( V i , ξ , y r ) in state ξ during year yr, respectively. Equation (23) constrains the feeder capacity limit, where I i j , m a x is the maximum permissible current in the feeder between i and j.
P ξ , y r l o s s = j = 1 N I i j , ξ , y r 2 R ( i j ) , i , ξ , y r
I i j , ξ , y r = P j , ξ , y r + j Q j , ξ , y r V j , ξ , y r * , i , ξ , y r
V m i n   V i , ξ , y r V m a x   , i , ξ , y r
V ξ , 1 , y r = 1   ,   δ ξ , 1 , y r = 0 ,   ξ , y r
0 I i j , ξ , y r I i j m a x ξ , y r

3.2.2. DG Constraints

Equations (24)–(26) constrain the output power of the DG between zero and the rated power.
( P _ D G ) 2 + ( Q _ D G ) 2 ( S _ D G _ r a t e d ) 2
0 P _ D G P _ D G _ r a t e d
0 Q _ D G Q _ D G _ r a t e d

3.2.3. BESS Constraints

Equations (27) and (28) define the charging and discharging power of the BESS and constrain its rated power.
P B E S S c h = P B E S S r a t e d ,                    P B E S S c h P B E S S _ c h r a t e d P B E S S c h   ,            P B E S S _ c h r a t e d   P B E S S c h 0
P B E S S d i s = P B E S S d i s     ,   0   P B E S S d i s P B E S S _ d i s r a t e d   P B E S S r a t e d    , P B E S S _ d i s   _ M a x   P B E S S d i s
where P B E S S c h is the charging power, P B E S S r a t e d is the rated charging power, P B E S S d i s is the discharging power of the BESS, and P B E S S r a t e d is the rated discharging power of the BESS.

3.3. Algorithm for Solving the Proposed Planning Model

In this study, the particle swarm optimization (PSO) algorithm was utilized to solve the optimization model, as it is gaining recognition as an efficient technique for solving this type of optimization problem. Probabilistic models of the demand, wind-based DG units, PV-based DG units, load growth, system topology, energy price, system data, planning period, discount rate, BESS data, and BESS costs were the main methodology’s inputs, and the main outputs of the optimization were the optimal size and location of the BESS units and the BESS units’ optimized operation. The proposed methodology’s justification was the optimization of the planning costs for distribution utilities through the deployment and management of BESSs. In the process, as shown in Figure 2, before the simulation started, a matrix consisting of the magnitude and the corresponding probability of each state of the load, wind turbines, and PV units was defined (this represented all of the operating scenarios); then, all required data of the system under study and the PSO algorithm were input to start the simulation. First, the simulation began with the first year to optimize the location, operation, and sizing of the BESSs, and this was controlled by all states of the load, wind, and PV. Secondly, a load flow analysis was performed to calculate the overall system losses and determine the required system upgrades. Finally, the total NPV of the planning was computed, and the most economical solution for PSO was obtained. Then, the algorithm moved to the next year; the load was increased by 5%, and the simulation was performed again until the final planning year was computed. After that, the particle velocity and position in the PSO algorithm were updated, and all steps were repeated until the simulation’s termination criteria were fulfilled; then, the final output results were obtained. Table 7 shows a comparison of our proposed model featuring a three-dimensional ESS dispatch with the probabilistic model outlined in Refs. [19,20], which was designed for a load-based ESS dispatch. This comparison took place for the scenario of a lithium-ion PV-based DG, and key performance metrics were assessed. Our model exhibited superior accuracy, resulting in a reduced total cost of 2.451 M$ in comparison with the load-based ESS dispatch probabilistic model’s 3.08 M$. Despite the high complexity introduced by our three-dimensional ESS dispatch strategy, our model demonstrated competitive computational efficiency, as it completed the analysis in about 2.5 h, whereas the previous probabilistic model required approximately 2 h. These findings underscore the advancements of our proposed model, which offers improved accuracy and complexity without compromising computational efficiency.

4. Case Studies and Results

4.1. The Network System under Study

The proposed planning model was evaluated by using the IEEE69 radial-distribution system, as shown in Figure 3, and its data can be found in [28]. The total system load was 3.802 MW and 2.7 MVAR, with base values of 12.66 KV and 100 MVA [29]. The planning period spanned 15 years, as proposed in this study, with an assumed annual load growth rate of 5%. The costs for off-peak and on-peak power were considered to be $23.6/MWh and $32.5/MWh, respectively [24]. The fixed and variable transmission line costs were assumed to be $150,000/km and $1000/MW, respectively. Furthermore, the fixed and variable substation costs were estimated to be $200,000 and $50,000/MW, respectively [30]. An interest rate of 10% was assumed, and the system operated with a power factor of 0.9.
Two renewable resources were used in this study as DGs: PV and wind-based DGs. The choice of DG power capacity was contingent on the capacity factor of each DG, and DGs with higher capacity factors possessed greater power ratings; in addition, the rating of DG power and the location of DGs in the system in this study depended on the fluctuations in the voltage profile and power loss of the system. Therefore, the sizing and placement of DGs were determined after studying the impacts of different locations and sizes on the voltage profiles and power losses of the system under study. The optimal size for the wind-based DG was determined to be 1 MW when connected to bus #61, while the optimal size for the PV-based DG was found to be 0.5 MW when connected to bus #17. It is noteworthy that the wind-based DG possessed a higher capacity than that of the PV-based DG, which was primarily due to its large capacity factor derived from the wind speed and solar irradiance data utilized in this study. Figure 4, Figure 5 and Figure 6 depict the voltage profiles, line currents, and active power losses of the system when different types of DGs operated at their rated capacities within the system.
Lead–acid batteries, lithium-ion batteries, redox flow batteries, and sodium sulfur batteries were chosen as candidate BESSs because they are the most suitable for use in power distribution and storage systems [18,19]. The power and energy of these BESSs were considered to be available on the market with step sizes of 100 KW/100 kWh, respectively. The candidate buses for the installation of BESSs in the system under study were assumed to include the following: 12, 16, 21, 29, 40, 48, 52, 59, 64, 67, and 69. The comprehensive cost of a BESS encompasses the capital expenditure, power conversion expenses, the balance of the plant cost, construction and commissioning expenditures, fixed operation and maintenance costs, and variable operational costs. Data on the costs and characteristics of BESSs are listed in Table 8.

4.2. Results and Discussion

Four cases are presented in this study, and each case had five scenarios in which the most economical BESS that has the lowest operation and investment costs and lowest power loss was found. These cases included those with no DGs, a PV-based DG, wind-based DG, and PV and wind-based DGs together. Each case had four scenarios: no BESS, lithium-ion batteries, lead–acid batteries, redox flow batteries, and sodium sulfur batteries, as shown in Table 9. All results will be discussed in detail in the following.

4.2.1. No-DG Case

In this case, it was assumed that the system did not have any DGs, and the impacts of the optimal placement, operation, and sizing of BESSs on the minimization of the investment costs and economic costs of operation of the BESSs were studied. Additionally, we examined the impacts of BESS optimization on system losses and upgrade costs.
In the first scenario, there were no BESSs in the system. The only cost components that were taken into consideration were the system upgrade costs and energy loss costs. In the following four scenarios, BESSs were added to the system, with each having investment and operation costs, as well as system upgrade and energy loss costs. As can be seen in Table 10, it was found that the optimal placement and sizing of the BESSs were at bus #64 and 100 kW, respectively. Each BESS has a different investment cost because BESSs have varying lifespans and different costs on the market. For instance, a lead–acid BESS has the highest NPV investment cost due to its lower lifespan, so it should be replaced many times during a planning period, while a lithium-ion BESS has the lowest investment cost. The economic cost of operation of each BESS is different due to the differences in their efficiency. For instance, lead–acid BESSs boast the highest cost, whereas lithium-ion BESSs have a lower cost than that of lead–acid and sodium sulfur BESSs, and they share the same operational expenses as those of redox flow BESSs due to their identical charging and discharging efficiency. The system upgrade cost was the same in all BESS scenarios, since they had the same placement and sizing. However, it was reduced by 10.45% in comparison with the first scenario due to the deferment of the system upgrade. For instance, feeder #1 was deferred from the third year to the fourth year in this case due to the installation of BESSs. The energy loss cost was slightly reduced by about 0.9% in comparison with that in scenario #1 due to the charging of BESSs during off-peak hours and discharging during on-peak hours. In this case, it was concluded that the lithium-ion BESS had the lowest cost, resulting in savings of 7.04% in comparison with the first scenario. Figure 7 shows the convergence of particle swarm optimization (PSO) in the case in which there were no DGs and in the scenario with lithium-ion batteries.

4.2.2. PV-Based DG Case

In this case, it was previously determined that the PV-based DG was located at bus #17 in order to study its impact on the placement and sizing of BESSs and the impacts of BESS optimization on the economic costs of investment and operation cost, the system cost, and the energy loss cost.
As can be noticed in Table 11, in the ‘No BESS’ scenario, the system cost remained the same as that in Case #1. This was because the calculation of the system upgrade cost was obtained in the worst-case scenario for the power flow in the network (which was at its maximum in the peak load state and the minimum PV-based state). Nevertheless, there was a 4.36% reduction in energy loss costs due to the installation of the PV-based DG. In the BESS scenario, it was found that the optimal placement and sizing of these BESSs were at bus #64 and 600 kW, respectively. These increases in BESS sizing are attributed to the charging power of the BESSs in the minimum and maximum off-peak load states of the PV-based DG. Each BESS had a different investment cost due to varying lifespans and different costs on the market. In this case, the investment cost significantly increased because of the increase in BESS sizing. For instance, lead–acid batteries incurred an investment cost of $1.011 million, and lithium-ion batteries boasted the lowest investment cost, amounting to just $0.57 million. The economic cost of operation of each BESS varied due to the efficiency of each BESS. In this case, the operation cost significantly increased due to an increase in the power discharge during the peak load state. It was observed that lithium-ion batteries had the lowest operation cost, while lead–acid batteries had the highest cost, since they had the lowest efficiency. The system upgrade cost in the BESS scenarios was the same for all BESS technologies because they shared the same placement and sizing, and their efficiencies were very close to each other. However, it was reduced by 38.3% in comparison with the first scenario due to the deferral of the system upgrade and by 31% in comparison with the first case in the same scenario—for instance, the upgrade decision for feeder #1 was deferred from the third year to the seventh year in the case of BESS inclusion. The energy loss cost experienced significant reductions from the previous case; for example, in the lithium-ion scenario, it was reduced by around 12% in comparison with the first scenario and by about 15.3% in comparison with the first case in the same scenario. It was concluded that the lithium-ion BESS was the most economical choice in this case, with savings of 26.2% in comparison with the first scenario and a 20.85% reduction in costs in comparison with the first case. In contrast, the lead–acid BESS was the least economical BESS option due to its high investment cost, lower efficiency, and shorter lifespan.

4.2.3. Wind-Based DG Cases

A.
Wind-based DG at 1 MW
In this case, it was previously determined that a 1 MW wind-based DG was located at bus #61 in order to study the impact of a wind-based DG on the placement and sizing of BESSs and the impact of BESS optimization on BESS planning, including the investment cost, operation cost, system cost, and energy loss costs.
Table 12 showcases the NPV of the total planning cost associated with the wind-based DG and BESSs. In the first scenario (no BESSs), the system cost remained the same as that in Case #1 and Case #2, and this was attributed to the calculation of the system upgrade that was obtained in the worst-case scenarios for the power flow in the network, (these occurred in the maximum peak load state and the minimum wind-based state). Nevertheless, there was a 16.1% reduction in the energy loss costs in comparison with Case #1 and a 5.6% reduction in comparison with Case #2. This was due to the installation of the wind-based DG. For the BESS scenarios, it was observed that the optimal placement and sizing for the BESSs were at bus #64 and 1100 kW, respectively. These increases in BESS sizing were attributed to the charging power of the BESSs in the minimum off-peak load state and the maximum state of the wind-based DG. In the following four scenarios, each BESS had a different investment cost because each of them had a different lifespan and different costs on the market. In this case, the investment cost significantly increased due to the increase in BESS sizing. For instance, lead–acid batteries had an investment cost of $1.58 million, and lithium-ion batteries boasted the lowest investment cost, amounting to just $1.05 million.
The economic cost of operation of each BESS varied due to the efficiency of each. In this case, the operation cost significantly increased due to an increase in power discharging during on-peak hours. For instance, the lead–acid battery boasted the highest cost, whereas the lithium-ion battery had a lower cost than that of the lead–acid battery and shared the same operational expenses as those of the redox flow battery due to their identical charging and discharging efficiency. The system upgrade cost in the BESS scenarios was the same for all BESS technologies because they shared the same placement and sizing, and their efficiencies were very close to each other, but it was reduced by 58% in comparison with the first scenario due to the deferral of the system upgrade decision. For instance, feeder #1 was deferred from the third to the tenth year, resulting in a 53.2% reduction in its system upgrade cost in comparison with the first case and a 32.2% reduction in comparison with Case #2 in the same scenarios.
The energy loss cost experienced significant reductions from those in the previous cases; for instance, in the lithium-ion scenario, it decreased by around 27.4% in comparison with the first scenario and by about 33.8% in comparison with the first case, and it decreased by 21.8% in comparison with Case #2 in the same scenarios. It was concluded that the lithium-ion BESS was the most economical choice in this case, with cost savings of 43.4% in comparison with the first scenario, a cost reduction of about 39.8% in comparison with the first case, and a reduction of approximately 23.85% in comparison with Case #2.
Table 13 shows the optimal operation and charging and discharging power of the BESSs in the case of a 1 MW wind-based DG (i.e., in each wind-based DG state and load state); this could help system operators control the operation of BESSs to maximize their benefits for the distribution system. It can be observed in Table 13 that in the initial load state, the charging of the BESSs increased with the increase in the wind state. Conversely, in the final load state, the discharging power of the BESSs decreased with the increase in the wind state. Therefore, the lowest state (i.e., load state = 1 and wind state = 1) had the lowest charging power. On the other hand, the highest state (i.e., load state = 12 and wind state = 12) had the lowest discharging power.
B.
Wind-based DG at 4 MW
In this case, the penetration level of the wind-based DG was increased to 4 MW, which was greater than the rated demand of the system under study. This was to demonstrate how the DG’s state 19 controlled the BESSs’ status (charging and discharging). In this case, the generated power was greater than the demand, which resulted in excess energy in some scenarios. In order to solve this issue, the system’s operator should curtail the excess power from the wind-based DG to reduce its impact on the system’s operation (i.e., maintaining the voltages and line currents within their operation limits), and this is not efficient from an economic perspective. Therefore, BESSs can solve this issue due to their ability to absorb excess power. Table 14 shows the optimal operation and charging and discharging power of BESSs in the 4 MW wind-based case (i.e., in each wind-based DG state and load state). Figure 8 illustrates the operation of BESSs and the amounts of charging and discharging power in different selected states of the load and wind (i.e., states #1, 4, 9, and 12 for the load state and wind state) in the 1 MW wind-based DG case and in the 4 MW wind-based DG case. Unlike in the previous case (i.e., the 1 MW wind-based DG), where the charging decisions of the BESSs took place in an off-peak load state and discharging decisions happened in an on-peak load state, the 4 MW wind-based DG case showcased that the BESSs were charged in the maximum wind state (i.e., state #12) regardless of the load state. This was because of the increase in the penetration level of the wind-based DG, which produced excess energy that could be stored by BESSs and used when needed (i.e., the on-load state and off-wind state). Figure 9 shows the voltage profile at each bus in the no-DG, 4 MW wind-based DG, and 4 MW-wind based DG cases with BESSs in the maximum load state and maximum wind state.

4.2.4. PV and Wind-Based DG Case

In this case, it was determined previously that a 500 kW PV-based DG and a 1 MW wind-based DG were connected to bus #17 and bus #61, respectively, in order to study the common impact of the PV-based DG and wind-based DG on the placement and sizing of BESSs and the impact of BESS optimization on BESS planning, including the BESS investment cost, economic cost of operation, system cost, and energy loss cost.
Table 15 showcases the NPV of the total planning cost associated with the PV and wind-based DGs and BESSs. In the first scenario (no BESSs), the system cost remained the same as that in Cases #1–#3, and this was attributed to the calculation of the system upgrades that was obtained in the worst-case scenarios for the power flow in the network (which occurred in the maximum on-peak load state and the minimum PV and wind-based state). Nevertheless, there was a significant reduction in the energy loss costs in comparison with Case #1, as well as a 14.1% reduction and 10% reduction in comparison with Case #2 and a 4.8% reduction in comparison with Case #3. This reduction was attributed to the installation of PV-based and wind-based DGs together in the system.
For the BESS scenario, it became essential to increase the number of BESS units to store a larger amount of energy due to the installation of two DGs with large capacities in the system. Therefore, two BESS units were installed in these cases. It was observed that the optimal placement and sizing of BESS #1 and BESS #2 were at bus #12 and bus #64 with capacities of 600 kW and 1100 kW, respectively. These increases in BESS sizing are attributed to the charging power of the BESSs during the minimum off-peak load state and the maximum output power state of the wind-based and PV-based DGs.
In the following four scenarios, each BESS unit had a different investment cost due to varying lifespans and different costs on the market. In this case, the investment cost significantly increased due to the increase in BESS sizing. Lead–acid batteries had an investment cost of $2.86 million, and lithium-ion batteries boasted the lowest investment cost, amounting to just $1.62 million. The operational cost of each BESS unit varied due to the efficiency of each. In this case, the economic cost of operation significantly increased due to the increase in power charging during the peak load hours. For instance, lead–acid batteries boasted the highest cost, whereas lithium-ion batteries had a lower cost than that of lead–acid batteries and shared the same operational expenses as those of redox flow batteries due to their identical charging and discharging efficiencies.
The system upgrade cost in the BESS scenarios remained the same for all BESS technologies because they shared the same placement and sizing, and their efficiencies were very close to each other. However, the system upgrade cost in the BESS scenario was reduced by 77.9% in comparison with that in the first scenario. This reduction was attributed to the deferral of the system upgrade, as shown in Table 16. For example, feeder #1 was deferred from the 3rd to the 13th year of the planning horizon, and feeders #6, 7, 8, and 9 did not require any upgrades during the planning period. This was due to the installation of BESS units with a large capacity in this case, resulting in a 75% reduction in comparison with the first case, a 64% reduction in comparison with Case #2, and a 47% reduction in comparison with Case #3 in the same scenarios.
The energy loss cost experienced a significant reduction from those in the previous cases. For example, in the lithium-ion scenario, it decreased by around 29.3% in comparison with the first scenario; it decreased by about 38.7% in comparison with Case #1, by 27.6% in comparison with Case #2, and by 7.4% in comparison with Case #3 in the same scenarios.
In conclusion, the lithium-ion BESS proved to be the most economical choice in this case, with total cost savings of 51.9% in comparison with the first scenario and total cost reductions of about 48.9% in comparison with the first case, about 35.5% in comparison with Case #2, and about 15.3% in comparison with Case #3.

4.3. Sensitivity Analysis

The sensitivity analysis of our proposed model involved varying the key parameters in PSO, such as by exploring different inertia weights (e.g., 0.1, 0.5, 0.9), adjusting the cognitive and social coefficients (e.g., 1.5, 2.0, 2.5), and introducing variations in the initial particle positions. Surprisingly, the results did not exhibit substantial differences when compared with our main PSO configuration with an inertia weight of 0.99 and cognitive and social coefficients of 1.5 and 2, respectively. This suggested that our primary model was robust and displayed a notable insensitivity to parameter variations for the specific optimization task. While these findings affirm the stability of our chosen parameters, future adjustments may be considered based on the requirements of different optimization scenarios.

5. Conclusions

In recent years, numerous innovative technologies such as ESSs and distributed energy resources (DERs) have been seamlessly integrated into distribution networks in response to various technical and economic considerations in order to deliver electricity at the most cost-effective rates while maintaining the necessary level of quality. However, the integration of these ESSs into power systems without thorough research, design, and optimized planning and operational processes could potentially compromise the overall quality and security of electrical grids. Consequently, the optimization of the planning and operation of ESSs is indispensable for bolstering their performance and enhancing the overall quality and reliability of power systems.
The primary objectives of this research were to integrate energy storage systems into distribution systems alongside renewable energy sources and to optimize their placement, operation, and sizing to maximize their benefits while minimizing the total planning and economic costs operation and considering both technical and economic aspects.
This study proposed a probabilistic planning model for BESSs in distribution networks to optimize the location, sizing, and operation of BESSs and to determine the most economically efficient BESS technologies that minimized the overall planning costs. Four different BESS technologies were examined and compared against the base scenario (i.e., without BESSs) and each other. It was concluded that the integration of these BESSs into a distribution system in conjunction with a wind-based DG and PV-based DG had a significant impact on reducing energy losses, deferring required system upgrades, and maximizing benefits through energy purchases during off-peak hours and energy selling during on-peak hours. The proposed model demonstrated that lithium-ion battery energy storage was the most cost-effective BESS option for the case under study.
The proposed model proved its effectiveness in accommodating all possible operating scenarios of the system, as well as the intermittent nature of renewable-based DGs. Unlike in previous research work, the proposed model allowed the optimizer to freely make charging and discharging decisions without imposing any unrealistic constraints, leading to more efficient outcomes.

Author Contributions

Conceptualization, A.A.A., M.A.A. and E.A.A.-A.; Funding acquisition, M.A.A.; Investigation, A.A.A., M.A.A. and E.A.A.-A.; Methodology, A.A.A. and M.A.A.; Resources, A.A.A. and M.A.A.; Software, A.A.A.; Supervision, M.A.A. and E.A.A.-A.; Writing—original draft, A.A.A.; Writing—review and editing, M.A.A. and E.A.A.-A. All authors have read and agreed to the published version of the manuscript.

Funding

This research and the APC were funded by K.A.CARE Energy Research and Innovation Center in Riyadh, King Saud University, Riyadh, Saudi Arabia.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

The authors extend their appreciation to K.A.CARE Energy Research and Innovation Center in Riyadh, King Saud University, Riyadh, Saudi Arabia for funding this work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Farhangi, H. The path of the smart grid. IEEE Power Energy Mag. 2010, 8, 18–28. [Google Scholar] [CrossRef]
  2. Hossain, E.; Faruque, H.; Sunny, M.; Mohammad, N.; Nawar, N. A Comprehensive Review on Energy Storage Systems: Types, Comparison, Current Scenario, Applications, Barriers, and Potential Solutions, Policies, and Future Prospects. Energies 2020, 13, 3651. [Google Scholar] [CrossRef]
  3. Hannan, M.; Wali, S.; Ker, P.; Abd Rahman, M.; Mansor, M.; Ramachandaramurthy, V.; Muttaqi, K.; Mahlia, T.; Dong, Z. Battery energy-storage system: A review of technologies, optimization objectives, constraints, approaches, and outstanding issues. J. Energy Storage 2021, 42, 103023. [Google Scholar] [CrossRef]
  4. Huang, W.-C.; Zhang, Q.; You, F. Impacts of battery energy storage technologies and renewable integration on the energy transition in New York State. Adv. Appl. Energy 2023, 9, 100126. [Google Scholar] [CrossRef]
  5. Rancilio, G.; Lucas, A.; Kotsakis, E.; Fulli, G.; Merlo, M.; Delfanti, M.; Masera, M. Modeling a large-scale battery energy storage system for power grid application analysis. Energies 2019, 12, 3312. [Google Scholar] [CrossRef]
  6. Rallabandi, V.; Akeyo, O.M.; Jewell, N.; Ionel, D.M. Incorporating battery energy storage systems into multi-MW grid connected PV systems. IEEE Trans. Ind. Appl. 2019, 55, 638–647. [Google Scholar] [CrossRef]
  7. Peng, S.; Zhu, L.; Dou, Z.; Liu, D.; Yang, R.; Pecht, M. Method of Site Selection and Capacity Setting for Battery Energy Storage System in Distribution Networks with Renewable Energy Sources. Energies 2023, 16, 3899. [Google Scholar] [CrossRef]
  8. Das, C.K.; Bass, O.; Kothapalli, G.; Mahmoud, T.S.; Habibi, D. Overview of energy storage systems in distribution networks: Placement, sizing, operation, and power quality. Renew. Sustain. Energy Rev. 2018, 91, 1205–1230. [Google Scholar] [CrossRef]
  9. Valencia, A.; Hincapie, R.A.; Gallego, R.A. Optimal location, selection, and operation of battery energy storage systems and renewable distributed generation in medium–low voltage distribution networks. J. Energy Storage 2021, 34, 102158. [Google Scholar] [CrossRef]
  10. Grisales-Noreña, L.F.; Montoya, O.D.; Gil-González, W. Integration of energy storage systems in AC distribution networks: Optimal location, selecting, and operation approach based on genetic algorithms. J. Energy Storage 2019, 25, 100891. [Google Scholar] [CrossRef]
  11. Johnson, R.C.; Mayfield, M.; Beck, S.B.M. Optimal placement, sizing, and dispatch of multiple BES systems on UK low voltage residential networks. J. Energy Storage 2018, 17, 272–286. [Google Scholar] [CrossRef]
  12. Wang, B.; Li, H.; Cao, X.; Zhu, X.; Gan, Z. The optimization of energy storage capacity for distribution networks with the consideration of probability correlation between wind farms based on PSO algorithm. IOP Conf. Ser. Earth Environ. Sci. 2017, 61, 012064. [Google Scholar] [CrossRef]
  13. Saini, P.; Gidwani, L. An investigation for battery energy storage system installation with renewable energy resources in distribution system by considering residential, commercial, and industrial load models. J. Energy Storage 2022, 45, 103493. [Google Scholar] [CrossRef]
  14. Zhe, J.; Liang, W.; Xuxiang, H.; Xueshan, H.; Dong, Y.; Ning, Z. Distributed Z.: Energy Storage Optimization Configuration of Active Distribution Network. IOP Conf. Ser. Mater. Sci. Eng. 2020, 729, 012010. [Google Scholar] [CrossRef]
  15. Wang, H.; Wang, J.; Piao, Z.; Meng, X.; Sun, C.; Yuan, G.; Zhu, S. The optimal allocation and operation of an energy storage system with high penetration grid-connected photovoltaic systems. Sustainability 2020, 12, 6154. [Google Scholar] [CrossRef]
  16. Wong, L.A.; Ramachandaramurthy, V.K.; Walker, S.L.; Taylor, P.; Sanjari, M.J. Optimal placement and sizing of battery energy storage system for losses reduction using whale optimization algorithm. J. Energy Storage 2019, 26, 100892. [Google Scholar] [CrossRef]
  17. Agamah, S.U.; Ekonomou, L. Energy storage system scheduling for peak demand reduction using evolutionary combinatorial optimization. Sustain. Energy Technol. Assess. 2017, 23, 73–82. [Google Scholar]
  18. Qiu, J.; Xu, Z.; Zheng, Y.; Wang, D.; Dong, Z. Distributed generation and energy storage system planning for a distribution system operator. IET Renew. Power Gener. 2018, 12, 1345–1353. [Google Scholar] [CrossRef]
  19. Awad, A.; EL-Fouly, T.; Salama, M. Optimal ESS Allocation for Benefit Maximization in Distribution Networks. IEEE Trans. Smart Grid 2017, 8, 1668–1678. [Google Scholar] [CrossRef]
  20. Awad, A.; EL-Fouly, T.; Salama, M. Optimal ESS Allocation and Load Shedding for Improving Distribution System Reliability. IEEE Trans. Smart Grid 2014, 5, 2339–2349. [Google Scholar] [CrossRef]
  21. Leou, R.C. An economic analysis model for the energy storage system applied to a distribution substation. Int. J. Electr. Power Energy Syst. 2012, 34, 132–137. [Google Scholar] [CrossRef]
  22. Xu, Y.; Singh, C. Adequacy and economy analysis of distribution systems integrated with electric energy storage and renewable energy resources. IEEE Trans. Power Syst. 2012, 27, 2332–2341. [Google Scholar] [CrossRef]
  23. Alotaibi, A.M.; Salama, M. A multi-state model for renewable resources in distribution systems planning. In Proceedings of the 2014 IEEE Electrical Power and Energy Conference (EPEC), Calgary, AB, Canada, 12–14 November 2014; pp. 48–53. [Google Scholar]
  24. Alotaibi, A.M.; Salama, M. An incentive-based multistage expansion planning model for smart distribution systems. IEEE Trans. Power Syst. 2018, 33, 5469–5485. [Google Scholar] [CrossRef]
  25. Alotaibi, A.M.; Salama, M. An efficient probabilistic-chronological matching modeling for DG planning and reliability assessment in power distribution systems. Renew. Energy 2016, 99, 158–169. [Google Scholar] [CrossRef]
  26. Boroumandfar, G.; Khajehzadeh, A.; Eslami, M.; Syah, R.B. Information gap decision theory with risk aversion strategy for robust planning of hybrid photovoltaic/wind/battery storage system in distribution networks considering uncertainty. Energy 2023, 278, 127778. [Google Scholar] [CrossRef]
  27. Park, J.; Liang, W.; El-Keib, A.A.; Shahidehpour, M.; Billinton, R. A probablistic reliability evaluation of a power system including solar/photovoltaic cell generator. In Proceedings of the 2009 IEEE Power and Energy Society General Meeting, Calgary, AB, Canada, 26–30 July 2009; pp. 1–6. [Google Scholar]
  28. Wong, S.; Bhattacharya, K.; Fuller, J.D. Electric power distribution system design and planning in a deregulated environment. IET Gener. Transm. Distrib. 2009, 3, 1061–1078. [Google Scholar] [CrossRef]
  29. El-Ela, A.A.A.; El-Sehiemy, R.A.A.; Kinawy, M.; Ali, E.S. Optimal placement and sizing of distributed generation units using different cat swarm optimization algorithms. In Proceedings of the 2016 Eighteenth International Middle East Power Systems Conference (MEPCON), Cairo, Egypt, 7–29 December 2017; pp. 975–981. [Google Scholar]
  30. Pisica, I.; Bulac, C.; Eremia, M. Optimal distributed generation location and sizing using genetic algorithms. In Proceedings of the 15th International Conference on Intelligent System Applications to Power Systems, ISAP, Curitiba, Brazil, 8–12 November 2009; pp. 1–6. [Google Scholar]
  31. Mongird, K.; Viswanathan, V.; Balducci, P.; Alam, J.; Fotedar, V.; Koritarov, V.; Hadjerioua, B. An evaluation of energy storage cost and performance characteristics. Energies 2020, 13, 3307. [Google Scholar] [CrossRef]
  32. Aquino, T.; Zuelch, C.; Koss, C. Energy Storage Technology Assessment; HDR Report No. 10060535-0ZP-C1001; Prepared for the Public Service Company of New Mexico by HDR; HDR: Omaha, NE, USA, 2017. [Google Scholar]
Figure 1. Power characteristic curve of a wind turbine.
Figure 1. Power characteristic curve of a wind turbine.
Sustainability 16 00290 g001
Figure 2. Flowchart of the proposed method.
Figure 2. Flowchart of the proposed method.
Sustainability 16 00290 g002
Figure 3. The IEEE69 distribution system.
Figure 3. The IEEE69 distribution system.
Sustainability 16 00290 g003
Figure 4. The voltage profiles for different DGs.
Figure 4. The voltage profiles for different DGs.
Sustainability 16 00290 g004
Figure 5. Line currents for different DGs.
Figure 5. Line currents for different DGs.
Sustainability 16 00290 g005
Figure 6. System power losses in different cases.
Figure 6. System power losses in different cases.
Sustainability 16 00290 g006
Figure 7. Convergence of PSO in the no-DG case and the lithium-ion scenario.
Figure 7. Convergence of PSO in the no-DG case and the lithium-ion scenario.
Sustainability 16 00290 g007
Figure 8. Optimal operation of BESSs in the 1 MW and 4 MW wind-based DG cases in different states. + indicates the charging mode and − indicates the discharging mode.
Figure 8. Optimal operation of BESSs in the 1 MW and 4 MW wind-based DG cases in different states. + indicates the charging mode and − indicates the discharging mode.
Sustainability 16 00290 g008
Figure 9. The voltage profile at each bus in the 4 MW wind-based DG case with and without BESSs.
Figure 9. The voltage profile at each bus in the 4 MW wind-based DG case with and without BESSs.
Sustainability 16 00290 g009
Table 1. Wind speed and wind turbine parameters.
Table 1. Wind speed and wind turbine parameters.
ParametersValue
Rated power ( P r )1000 kW
Cut-in speed ( v i n ) m/s3
Rated speed ( v r ) m/s14
Cut-out speed ( v o u t ) m/s25
Weibull scale parameter c4.2483
Weibull shape parameter k1.6515
Table 2. Multi-state probability model of a wind-based DG.
Table 2. Multi-state probability model of a wind-based DG.
StateLower (m/s)Upper (m/s)% Output PowerProbability
103 and >2500.4305
234.150.18007
34.15.2150.14195
45.26.3250.10046
56.37.4350.06501
67.48.5450.0389
78.59.6550.0217
89.610.7650.01134
910.711.8750.00558
1011.812.9850.00259
1112.914950.00114
1214251000.000772
Table 3. Solar irradiance and PV parameters.
Table 3. Solar irradiance and PV parameters.
Parameters Value
Rated PV power ( P P V r )500 kW
Adjustable value of irradiance r C 200 W/m2
Standard irradiance s S T D 1000 W/m2
Beta parameter β 0.45
Beta parameter α 1.438
Table 4. Multi-state probability model of a PV-based DG.
Table 4. Multi-state probability model of a PV-based DG.
StateLower (kW/m2)Upper (kW/m2)% Output PowerProbability
100.08400.395786
20.0840.1687.940.138345
30.1680.252210.098823
40.2520.33429.30.076266
50.3340.41837.60.064414
60.4180.502460.054077
70.5020.58654.40.045772
80.5860.6762.80.03867
90.670.75471.20.032253
100.7540.83879.60.026061
110.8380.922880.019489
120.922196.10.010005
Table 5. Multi-state probability model of demand.
Table 5. Multi-state probability model of demand.
State Lower Load (pu)Upper Load (pu)Probability
100.350.03402
20.350.410.045205
30.410.470.08042
40.470.530.1208
50.530.590.1532
60.590.650.164
70.650.710.14825
80.710.770.1131
90.770.830.0729
100.830.890.0397
110.890.950.01821
120.9510.00634
Table 6. The column matrix for all possible scenarios for all states.
Table 6. The column matrix for all possible scenarios for all states.
Scenarios Demand   Power   ( P L ) PV   Power   ( P P V ) Wind   Power   ( P w )
1LLL
2LLH
3LHL
----
----
----
----
----
----
----
----
NHHH
L represents a lower-power state and H represents a higher-power state; N = 12 × 12 × 12.
Table 7. Comparison of our proposed model with an existing model from the literature.
Table 7. Comparison of our proposed model with an existing model from the literature.
Proposed ModelObjective Function
For Lithium-Ion Scenario
AccuracyComplexityComputation Time
Probabilistic [19,20]
(load-based ESS dispatch)
3.08 M$HighModeratedAbout 2 h
Our proposed model
(three-dimensional ESS dispatch)
2.451 M$Very HighRelatively HighAbout 2.5 h
Table 8. Data on the cost and characteristics of BESSs [31,32].
Table 8. Data on the cost and characteristics of BESSs [31,32].
Lead–AcidLithium IonRedox Flow Sodium Sulfur
Capital Cost—Energy
Capacity ($/kWh)
260271555661
Power Conversion
System ($/kW)
350288350350
Balance of Plant ($/kW)100100100100
Construction and
Commissioning Cost
($/kWh)
176101190133
Fixed O&M ($/kW-yr)10101010
Efficiency 85%95%95%90%
Life Span (yr)3101515
Table 9. The cases and scenarios.
Table 9. The cases and scenarios.
CasesDG TypeDG LocationScenarios
1No DG--1. No BESSs
2. Lead–acid
3. Lithium ion
4. Redox flow
5. Sodium sulfur
2PV171. No BESSs
2. Lead–acid
3. Lithium ion
4. Redox flow
5. Sodium sulfur
3wind611. No BESSs
2. Lead–acid
3. Lithium ion
4. Redox flow
5. Sodium sulfur
4PV–wind61–171. No BESSs
2. Lead–acid
3. Lithium ion
4. Redox flow
5. Sodium sulfur
Table 10. Results for the first case (no DGs).
Table 10. Results for the first case (no DGs).
No DGsNo BESSsLead–AcidLithium IonRedox Flow Sodium Sulfur
NPV of Investment Cost (M$)00.1690.0950.1270.130
NPV of Economic Cost of Operation (M$)0−0.008−0.0094−0.0094−0.0089
NPV of System Upgrade Cost (M$)3.0332.7162.7162.7162.716
NPV of Energy Loss (M$)0.2980.29540.29520.29520.2953
NPV of Total Cost (M$)3.33143.17123.09683.12873.1323
Saving %04.87.046.15.98
Execution Time (h)≈1 h
Table 11. Results for the second case (PV).
Table 11. Results for the second case (PV).
PV-Based DGNo BESSLead–AcidLithium IonRedox Flow Sodium Sulfur
NPV of Investment Cost (M$)01.0110.57230.7640.781
NPV of the Economic Cost of Operation (M$)0−0.271−0.243−0.243−0.23
NPV of System Upgrade Cost (M$)3.0331.8721.8721.8721.872
NPV of Energy Loss (M$)0.2850.2530.24990.24990.251
NPV of Total Cost (M$)3.31872.91802.45122.64282.6743
Saving %012.126.220.419.4
Execution Time (h)≈2.5 h
Table 12. Results for the third case (wind).
Table 12. Results for the third case (wind).
Wind-Based DGNo BESSLead–AcidLithium IonRedox Flow Sodium Sulfur
NPV of Investment Cost (M$)01.851.051.41.43
NPV of Economic Cost of Operation (M$)0−0.577−0.644−0.644−0.611
NPV of System Upgrade Cost (M$)3.0331.271.271.271.27
NPV of Energy Loss (M$)0.2690.19970.19540.19540.1975
NPV of Total Cost (M$)3.302.74261.86662.21782.2857
Saving %016.943.432.830.7
Execution Time (h) ≈3 h
Table 13. Optimal operation and capacity of BESSs (kW) in the 1 MW wind-based case.
Table 13. Optimal operation and capacity of BESSs (kW) in the 1 MW wind-based case.
StateWind State
1 2 3 4 5 6 7 8 9 10 11 12
168118218318418.0004518.0018618.002718.0119818.0231918.041018.0641068.078
254.38101.88196.88291.88386.8805481.882576.886671.8937766.9067861.9263956.95371004.471
340.7681.16160.76240.7601323.7605404.1622484.767565.7754646.3901726.8124807.8436848.3628
427.1465.54142.14219.1401295.1405372.1424449.143525.757602.1735679.1984756.2333794.3749
513.5249.52122.52195.5201282.0206341.5226413.528487.5386560.1567633.5842703.6227742.6466
6−0.134.9103.9169.9001254.9006310.9028379.906449.9201518.9398589.9698658.012690.0382
7−703.72−671.72−603.72−532.72−463.719−393.717−324.711−268.698−183.677−117.645−48.1989−13.5705
8−756.34−717.34−647.34−573.94−501.339−427.337−355.33−295.817−209.294−136.259−63.7101−27.1795
9−807.96−769.96−692.96−615.96−539.559−462.957−385.95−308.936−232.911−155.874−79.1215−40.7887
10−904.58−821.58−740.58−660.18−579.579−498.577−417.569−337.554−254.529−174.489−94.433−54.3981
11−1018.2−970.7−918.9−780.7−685.699−2158.2−495.689−400.673−305.646−210.604−115.545−68.0077
12−1079.5−1029.5−929.55−829.55−729.549−629.546−529.539−429.522−329.494−229.45−129.388−79.3492
− indicates the discharging mode.
Table 14. Optimal operation and capacity of BESSs (kW) in the 4 MW wind-based case.
Table 14. Optimal operation and capacity of BESSs (kW) in the 4 MW wind-based case.
StateWind State
1 2 3 4 5 6 7 8 9 10 11 12
1−1262.7−862.7−62.7737.31537.31868.0022268.052668.0122068.0233468.043868.0644068.078
2−1466.42−1066.42−266.42533.581333.581754.3822184.382554.3942964.4073350.7263738.5544004.471
3−2050.34−1656.44−868.64−80.8399706.96051640.7622070.762440.7752808.793177.3123550.8443918.863
4−2252.9−1869.86−1100.56−331.56437.44051206.4421975.452317.1572665.1733016.6983367.2333717.255
5−2456.5−2081.93−1330.23−578.53173.1706924.87261676.582193.5392526.0572861.0843196.1233530.647
6−2660.1−2294−1559−824−88.9994646.00281381.012115.922392.442711.473030.5123350.038
7−2865.7−2505.72−1786.72−1067.72−348.719370.2831089.291808.3022153.8232568.6552872.5013176.429
8−3068.9−2717.04−2013.34−1309.44−605.63998.16312801.971505.6831905.8062429.2412718.293010.821
9−3272.7−2928.16−2239.16−1550.16−861.159−172.157516.8511205.8641605.8892301.7262577.3792849.211
10−3476.4−3138.58−2462.98−1787.38−1111.78−436.177239.431915.04551315.0712266.3112583.0672701.602
11−3680.1−3348.85−2686.4−2023.85−1361.4−698.847−36.3389625.82671025.8541951.2462421.9552556.992
12−3881.6−3556.55−2906.55−2256.55−1606.55−956.546−306.539343.4777743.5061643.552295.6122420.651
− indicates the discharging mode.
Table 15. Results for fourth case (PV and wind).
Table 15. Results for fourth case (PV and wind).
Wind-Based DGNo BESSLead–AcidLithium IonRedox Flow Sodium Sulfur
NPV of Investment Cost (M$)02.8641.6222.1652.214
NPV of the Economic Cost of Operation (M$)0−0.798−0.892−0.892−0.845
NPV of System Upgrade Cost (M$)3.0330.670.670.670.67
NPV of Energy Loss (M$)0.2560.1850.1810.1810.183
NPV of Total Cost (M$)3.28932.92121.58122.12412.2221
Saving %011.251.935.432.4
Execution Time (h)≈8 h
Table 16. Years in which distribution feeders were upgraded with and without BESSs in different cases.
Table 16. Years in which distribution feeders were upgraded with and without BESSs in different cases.
StateNO DGsPV-Based DGWind-Based DG PV- and Wind-Based DGs
Upgrade Year without BESSsUpgrade Year with BESSsUpgrade Year without BESSsUpgrade Year with BESSsUpgrade Year without BESSsUpgrade Year with BESSsUpgrade Year without BESSsUpgrade Year with BESSs
13437310313
23437310313
33437310313
45559512515
534383123No upgrade required during planning period
634383123No upgrade required during planning period
734393133No upgrade required during planning period
845494134No upgrade required during planning period
9565105155No upgrade required during planning period
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Alguhi, A.A.; Alotaibi, M.A.; Al-Ammar, E.A. Probabilistic Planning for an Energy Storage System Considering the Uncertainties in Smart Distribution Networks. Sustainability 2024, 16, 290. https://doi.org/10.3390/su16010290

AMA Style

Alguhi AA, Alotaibi MA, Al-Ammar EA. Probabilistic Planning for an Energy Storage System Considering the Uncertainties in Smart Distribution Networks. Sustainability. 2024; 16(1):290. https://doi.org/10.3390/su16010290

Chicago/Turabian Style

Alguhi, Ahmed A., Majed A. Alotaibi, and Essam A. Al-Ammar. 2024. "Probabilistic Planning for an Energy Storage System Considering the Uncertainties in Smart Distribution Networks" Sustainability 16, no. 1: 290. https://doi.org/10.3390/su16010290

APA Style

Alguhi, A. A., Alotaibi, M. A., & Al-Ammar, E. A. (2024). Probabilistic Planning for an Energy Storage System Considering the Uncertainties in Smart Distribution Networks. Sustainability, 16(1), 290. https://doi.org/10.3390/su16010290

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

Article Metrics

Back to TopTop