Next Article in Journal
Topographical Optimization of a Battery Module Case That Equips an Electric Vehicle
Next Article in Special Issue
End-to-End Direct-Current-Based Extreme Fast Electric Vehicle Charging Infrastructure Using Lithium-Ion Battery Storage
Previous Article in Journal
High-Performance Anodes Made of Metallic Lithium Layers and Lithiated Silicon Layers Prepared by Vacuum Technologies
Previous Article in Special Issue
Data-Driven Thermal Anomaly Detection in Large Battery Packs
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Capacity and Cost Analysis of Battery Energy Storage System in Standalone Microgrid Considering Battery Lifetime

by
Pinit Wongdet
1,
Terapong Boonraksa
2,
Promphak Boonraksa
3,
Watcharakorn Pinthurat
4,
Boonruang Marungsri
1,* and
Branislav Hredzak
4,*
1
School of Electrical Engineering, Suranaree University of Technology, Nakhon Ratchasima 30000, Thailand
2
School of Electrical Engineering, Rajamangala University of Technology Rattanakosin, Nakhon Pathom 73170, Thailand
3
School of Electrical Engineering, Faculty of Engineering and Architecture, Rajamangala University of Technology Suvarnabhumi, Nonthaburi 11000, Thailand
4
School of Electrical Engineering and Telecommunications, The University of New South Wales, Sydney 2052, Australia
*
Authors to whom correspondence should be addressed.
Batteries 2023, 9(2), 76; https://doi.org/10.3390/batteries9020076
Submission received: 30 November 2022 / Revised: 18 January 2023 / Accepted: 22 January 2023 / Published: 23 January 2023
(This article belongs to the Collection Recent Advances in Battery Management Systems)

Abstract

:
In standalone microgrids, the Battery Energy Storage System (BESS) is a popular energy storage technology. Because of renewable energy generation sources such as PV and Wind Turbine (WT), the output power of a microgrid varies greatly, which can reduce the BESS lifetime. Because the BESS has a limited lifespan and is the most expensive component in a microgrid, frequent replacement significantly increases a project’s operating costs. This paper proposes a capacity optimization method as well as a cost analysis that takes the BESS lifetime into account. The weighted Wh throughput method is used in this paper to estimate the BESS lifetime. Furthermore, the well-known Particle Swarm Optimization (PSO) algorithm is employed to maximize battery capacity while minimizing the total net present value. According to simulation results, the optimal adjusting factor of 1.761 yields the lowest total net present value of US$200,653. The optimal capacity of the BESS can significantly reduce the net present value of total operation costs throughout the project by extending its lifetime. When applied to larger power systems, the proposed strategy can further reduce total costs.

1. Introduction

Recently, Renewable Energy Resources (RESs) have become attractive sources in the electrical power system due to various advantages, such as being sustainable sources and being environmentally friendly [1,2]. The primary RESs are solar and wind energy generation sources. However, power generation from the RESs is intermittent because of its behavior. Therefore, renewable energy generation sources cannot be directly connected to the existing system or used without the appropriate management system [3,4]. A microgrid is a small power system constructed to manage Distributed Generators (DGs) from renewable energy and load clusters. The microgrid that connects to the bulk power system is called to be in “on-grid mode”, and when it disconnects from the bulk power system in an emergency, it is called to be in “islanded mode”. The microgrid that operates in the islanded mode can be called a standalone microgrid [5,6]. The standalone microgrid is common in remote or rural areas. Thus, the standalone microgrid’s stability is low due to the lack of assistance from the bulk power system. The primary goal in the microgrid is to balance the power load and power generation without interrupting the load. Hence, it is necessary to have manageable power generation sources, such as diesel engine generators or the BESS, to control the microgrid’s power [7,8]. A standalone microgrid consisting of renewable energy generation sources, load clusters, and BESS is widely used [9] and avoids using diesel generators for environmental reasons. One of the significant problems for BESS applications is finding optimal capacity that considers the lifetime of BESS. Because of the high cost of the BESS, BESSs with a short life have been widely used in a microgrid. The performance assessment of the grid is primarily based on cost and reliability related to the system’s lower greenhouse gas emissions [10]. The optimal size of PV and BESS in a microgrid for maintaining the longest BESS lifetime using the Loss of Power Supply Probability (LPSP) and only the initial cost as the indicators was presented in [11]. The cost reduction analysis was estimated for the battery State of Health (SOH) and artificial intelligence was applied to help with the battery management system. A proper battery health assessment can extend the battery’s lifetime and reduce the cost of the battery [12]. An accurate SOH is essential to optimizing battery lifespan. Battery deterioration tracking and SOH estimates are based on electrochemical impedance spectroscopy and the distribution of relaxation time measurements. They were used to look for indicators linked to deterioration [13]. The BESS lifetime estimation in a PV system using a practical model was presented by [14]. The authors also considered the impact of BESS sizes and types on its lifetime. The BESS lifetime increases with increased BESS size, and upfront costs also increase. The authors also introduced a strategy to optimise the total cost, including upfront costs and the replacement cost of BESS. The lifetime prediction method and sizing of lead-acid BESS in microgrids were applied by varying the BESS’s size and the weighted Wh throughput method to estimate the BESS lifetime. However, the authors did not consider the cost analysis to obtain the optimal size [15]. This paper proposes an optimization of BESS on a microgrid with PV and WT as the power generation source. The objective function of this study is to reduce total costs. According to the literature review, the articles in the references [11,14] used a microgrid system with only PV and BESS. References [12,13] considered the SOH estimation method to reduce costs, but this paper uses the optimization method. The work in [14] only studied the BESS lifetime estimation in a PV system. Moreover, the authors in [15] did not consider the cost analysis to obtain the optimal size. Therefore, in this paper, both PV and WT generation systems are analyzed. The sizing and optimal NPV of the BESS were analyzed with the optimization method. Implementing a microgrid system requires forecasting investment costs and profits, as well as maintenance over its lifecycle because the batteries are expensive components of the microgrid system. If the battery is replaced prematurely, the cost of the system will increase. Forecasting and estimation methods are generally used for the life cycle and the replacement of the battery. However, this paper proposes optimization to get the best results and reduce the total cost of the BESS system. The weighted Wh method and the PSO algorithm are applied for optimizing the cost of BESS. In a standalone microgrid system, prolonging the life of the equipment is necessary to reduce the cost of its replacement. However, the size and installation costs of the storage systems must be appropriate. Therefore, this paper provides an appropriate weighting to minimize the cost of the microgrid system. The PSO algorithm is applied to determine the optimum weighting parameter, which results in faster and more accurate analysis results than statistical analysis.
The major contributions of the paper are summarized as follows: The capacity optimization method and cost analysis of the BESS by taking into account the BESS lifetime are proposed in this paper. The BESS lifetime is estimated using the weighted Wh throughput method, which is fully described in [16]. The total cost of the microgrid and the notion of Net Present Value ( N P V ) constructs serve as the objective function. The objective function is solved using the PSO algorithm to determine the BESS’s ideal capacity. Due to its dependability and affordability, the lead-acid battery is used in this study. According to an economic analysis, it is a better choice than a modern battery like a lithium-ion battery [17,18].
The rest of this paper is structured as follows: Section 2 explains the mathematical modeling of DGs and BESS. Section 3 introduces the optimization model which consists of the objective function and the PSO algorithm. Section 4 provides the simulation results and result discussions. Finally, Section 5 concludes the paper.

2. DGs and BESS Models

In this section, the mathematical models of PV, WT and BESS used in the proposed optimization problem are briefly explained. A small industrial load is used for the case study in which PV and WT power generation systems are installed. The battery energy storage systems are used for power demand periods where the DGs are unable to supply the load for only some periods. Hence, BESS is small in size, and costs are reduced accordingly. However, the proper size of a BESS affects its longevity and maintenance or replacement costs. This model can be applied to the larger systems.

2.1. PV Model

The amount of PV power depends on the solar intensity and temperature and is given by [16],
P P V = P S T C · G c G S T C · [ 1 + k ( T c ) T S T C ] ,
where P P V is the output power; P S T C is the rated output power at standard conditions (the solar irradiance G S T C of 1000 W/m 2 and the temperature T S T C of 25 C; G C is the solar irradiance during operation; k is the discrete-time index and T c is the temperature during operation. The installed capacity of the PV system in this paper is 65 kW p .

2.2. Wind Turbine Model

The amount of generated WT power depends on the wind speed [19]. The output power of the WT system is expressed by,
P w t = 0 , v v c i OR v v c o , P r a t e d w t ( v v c i ) ( v r v c i ) , v c i v v r , P r a t e d w t , v r v v c o ,
where P W T is the output power of the wind turbine; v is the wind speed and P r a t e d w t is the rated output power of WT. The capacity of the WT system is selected as 40 kW p .

2.3. BESS Model

In this paper, a lead-acid battery is used for the calculation of the BESS cost because it is more cost-effective and safer compared to Li-ion battery. Although price of the Li-ion battery is continuing to decrease, it is still expensive in Thailand. In Thailand, the batteries widely used for energy storage in PV power generation systems are lead-acid batteries. In order to simulate the operation of the BESS, mathematical models for calculating the charge and discharge parameters and State of Charge (SOC) of the BESS are required. Batteries are bidirectional energy sources, meaning they can store energy in a manner comparable to an electrical load. In contrast, they provide energy as an energy source. The battery energy balance equation and its constraint are,
E B ( t + Δ t ) = E B ( t ) + ( P B , c h ( t ) · η B , c h P B , d ( t ) / η B , d ) Δ t ,
E B m i n E B ( t ) E B m a x ,
where E B is the energy of the BESS at time t; P B , c h , P B , d are the charging and discharging powers of BESS; η B , c h , η B , d are the charging and discharging efficiency of BESS respectively.
The maximum permissible output power of the BESS is limited as,
P B , c h m a x P b ( t ) B B , d m a x ,
where P B , c h m a x , P B , d m a x are the maximum charging and discharging powers of BESS.
The battery SOC and DOD can be obtained as,
S O C ( t ) = E B ( t ) E B , r a t e d , 0 S O C ( t ) 1 ,
D O D = 1 S O C ,
where SOC is the battery state of charge; E B is the energy of BESS at time t; E B , r e t e d is the rated energy and D O D is the battery depth of discharge.
The S O C = 1 or 0 means that the battery is fully charged or completely discharged respectively. The SOC at time t can be expressed as,
S O C ( t + Δ t ) = S O C ( t ) · ( 1 β ) + E B ( t + Δ t ) E B , r a t e d ,
where β is the battery’s self-discharge rate.
The battery SOC should be maintained within the specific ranges that depend on the type of the battery,
S O C m i n S O C ( t ) S O C m a x ,
where S O C m i n , S O C m a x are the minimum and maximum limits of battery SOC.

2.4. BESS Capacity Model

BESS is used to maintain the microgrid’s power balance. Thus, the required power of a BESS is the power mismatch between the load and the power generation. The balance between load and generation can be achieved through,
P B , r e q u i r e d ( t ) = P L ( t ) P W T ( t ) P P V ( t ) ,
where P B , r e q u i r e d is the BESS’ required power; P L is the load demand; P W T is the output power of WT and P P V is the output power of PV.
Taking into account the charging and discharging efficiencies, the BESS power and its rated power can be rewritten as follows: [20],
P B ( t ) = P B , r e q u i r e d ( t ) / η B , d , P B , r e q u i r e d ( t ) > 0 , P B , r e q u i r e d ( t ) × η B , c h , P B , r e q u i r e d ( t ) < 0 , 0 , P B , r e q u i r e d ( t ) = 0 ,
P B , r a t e d = max ( | P B ( t ) | ) ,
where P B , r a t e d is the BESS’s rated power.
The required energy from BESS can be calculated as,
E B , r e q u i r e d = max ( E B ( t ) ) min ( E B ( t ) ) .
The SOC of BESS should be kept within a specific range. Therefore, the rated energy considering DOD limits is expressed by
E B , r a t e d = E B , r e q u i r e d D O D m a x D O D m i n ,
where D O D m a x , D O D m i n are the DOD limits of the BESS.
In order to reduce the SOC level of the BESS, oversizing of the BESS is considered in this paper, as discussed in the next section. The oversized BESS does not have the same lifetime as an optimal size BESS [21,22]. BESS’s non-optimal size also increases investment costs. An adjusting factor is used to obtain the BESS’ oversized energy as,
E B , O r a t e d = E B , r a t e d × q B ,
where E B , O r a t e d is the BESS’s oversized energy, E B , r a t e d is the BESS’ rated energy and q B is the adjusting factor.
The adjusting factor q B is equal to or greater than 1. As the size of the battery can affect the DOD range, it follows that a battery with a higher capacity will last longer under identical conditions. However, the cost of installation will rise. Consequently, the parameter that must be optimized in this paper is q B . The optimal values will result in longer battery life and a reduction in total system costs.

2.5. BESS Lifetime Estimation

In this paper, a weighted Wh throughput method is used to calculate the lifetime of a lead-acid battery. This method has been employed by many researchers [15,23], as it can roughly estimate the battery lifetime by using the SOC as the most significant factor for battery lifetime. The battery has an overall throughput calculated from the manufacturer’s datasheet, which usually provides the relationship between the Cycles to Failure (CTF) and the different DODs as shown in Figure 1. In this study, the employed battery is the Sun Xtender PVX-2580L model (12 V, 3.096 kWh) [24].
The energy throughput corresponding to a specific DOD is expressed as,
E t h r o u g h p u t ( D O D ) = 2 × E B , r a t e d × D O D × C T F ,
where E t h r o u g h p u t is the energy throughput corresponding to a specific DOD and C T F is the cycles to failure.
Because the DOD varies within the range of DOD limits, the energy throughput varies with the DOD, as shown in Figure 2. By considering Equation (16) and the data shown in the Figure 1, the relationship between energy throughput and DOD in the Figure 2 can be obtained. The nonlinear curve in the Figure 2 is caused by the nonlinear function in the Figure which is the relationship between the Cycles to Failure (CTF) and different DODs. The data from the Figure 1 was obtained from the datasheet of the lead-acid battery in [24]. The overall energy throughput is the average energy throughput between the DOD limits and can be determined as,
E t h r o u g h p u t , a v g = a v g { E t h r o u g h p u t ( D O D ) } D O D ( m i n ) D O D ( m a x ) .
The cumulative energy throughput E W s o c multiplied by the SOC weighting factor W S O C during period Δ t is given as,
E W s o c = W S O C × | P B ( t ) | × Δ t .
The relationship between SOC and the weighting factor W S O C is shown in Figure 3. For example, using 1 Wh of battery at 0.2 p.u. of SOC is equivalent to 1.3 Wh of the cumulative energy throughput. On the other hand, using 1 Wh at 1 p.u. of SOC is equal to 0.55 Wh of the cumulative throughput. Based on the Figure 3 that shows the relationship between SOC and SOC weighting factor and on Equations (18) and (19), battery life loss mainly depends on SOC and SOC weighting factor ( S O C and W S O C ). If the SOC of the battery is kept lower, the battery life loss will be affected (see Equation (19)). Thus, the battery should operate at a high SOC or within its upper SOC limits to extend its lifetime. It has to be noted that this method ignores many influencing factors for simplicity, such as temperature, because temperature can be controlled. Therefore, the battery life loss can be calculated as,
L l i f e = E W s o c E t h r o u g h p u t , a v g ,
where L l i f e is the battery life loss and E t h r o u g h p u t , a v g is the average energy throughput.
The battery is considered to reach the end of life when the total cumulative energy throughput is equal to the overall energy throughput.
The following exponential function can be used to express the relationship between the number of cycles and the depth of discharge,
N C T F = f ( D O D ) = a × e b × D O D B + c × e d × D O D B ,
where a , b , c , d are constants.
The depreciation ( D B ) of BESS can then be calculated as,
D B = i = 1 m N B , c y c l e ( D O D i ) N C T F ( D O D i ) ,
where N B , c y c l e is the number of cycles according to the depth of discharge.
Therefore, the BESS lifetime L B can be obtained as,
L B = 1 D B , 1 y e a r × y e a r s ,
where D B , 1 y e a r is the depreciation of the battery in one year.
Hence, the total number of BESS replacements N R B throughout the project is given by,
N R B = L p L B 1 .
where L p is the project life.

3. Optimization Model

The optimization problem is introduced in this section, and the PSO algorithm used to solve it is briefly described.

3.1. Objective Function

The cost of the microgrid is represented by the concept of the Net Present Value ( N P V ). The total cost includes the initial cost, fixed operation and maintenance costs, and replacement costs. Therefore, the objective function used to minimize the N P V of the total cost can be defined as:
Min ( N P V C t o t a l ) = m i n ( C i n i t i a l + N P V C o m + N P V C r e ) ,
where N P V C t o t a l represents the net present value of the total cost; N P V C o m is the net present value of the fixed operation and maintenance costs; N P V C r e is the net present value of the replacement cost and C i n i t i a l is the initial cost.
The initial cost, the net present values of the replacement cost, and the fixed operation and maintenance costs respectively are given as,
C i n i t i a l = ( k B E , i n i t i a l × E B , O r a t e d ) + ( k B P , i n i t i a l × P B , r a t e d ) ,
N P V C r e = n = 1 N R B C i n i t i a l ( 1 + d ) N R B × L B ,
N P V C o m = k B , o m × E B , O r a t e d × U S W ,
where k B E , i n i t i a l is the BESS initial cost per kWh ($/kWh); E B , O r a t e d is the BESS’ oversized energy; k B P , i n i t i a l is the BESS initial cost per kW ($/kW); k B , o m is the operation and maintenance costs set to 5% of the initial cost ($/kWh/year); d is the discount rate (5%); L B is the BESS lifetime and U S W is the uniform series present worth factor.
The uniform series present worth factor is given as,
U S W = ( 1 + d ) L p 1 L p ( 1 + d ) L p ,
where d is the discount rate.
Figure 4 depicts the procedure flowchart for the BESS capacity determination. It starts by obtaining the input power of WT, PV, and load, and then calculating the rated power and energy capacity of the battery. Then, it estimates the BESS lifetime using the BESS model and obtains the objective function’s value. If N P V C t o t a l is minimal, the calculation ends. If not, q B is adjusted by incremental increases of 0.001 between 1 and 5 and the procedure is repeated.

3.2. PSO Algorithm

Particle swarm optimization (PSO) is an intelligent optimization algorithm introduced by J. Kennedy and R. Eberhart [25]. It has been widely used in many works owing to its effectiveness and simplicity. The basics of PSO can be found in [26,27,28]. The reason for choosing the PSO algorithm is that it has a simple structure with few parameters and provides a quick solution. It creates particle swarms, which simulate a system or structure of a group of organisms that group together for some specific purpose in each organism. The PSO chooses each i-th particle initial position and velocity at random. Then the i-th particle position and velocity are updated until the termination conditions are met. The particle’s position and velocity of the particle i are expressed by,
x i ( k + 1 ) = x i ( k ) + v i ( k + 1 ) ,
v i ( k + 1 ) = w × v i ( k ) + c 1 r 1 i ( p i , b e s t x i ( k ) + c 2 i r 2 i ( g i , b e s t x i k ) ) ,
where x i is the position of particle i; v i is the velocity of the particle i; w is the inertia factor; c 1 , c 2 are the acceleration factors; r 1 , r 2 are the random number [0–1]; p b e s t is the best particle and g b e s t is the best global solution.
Both acceleration factors ( c 1 , c 2 ) are commonly set to 2 [29]. The inertia factor which linearly decreases in each iteration i t e is calculated as,
w i t e = w m a x i t e × w m a x w m i n i t e m a x ,
where w m a x , w m i n are the weights of the maximum and minimum speeds and i t e m a x is the maximum iteration.
This paper applies the PSO algorithm to determine the optimal installation size of the BESS based on the battery size multiplier and the cost of the energy storage system throughout the project as a measure of fitness or sufficiency. Each particle in the herd has the battery size adjuster ( q B ) as a variable and is assigned to represent the answer point. The initialization parameters for the PSO are the number of particles ( n = 50), the maximum number of iterations ( i t e m a x = 100), the range of velocity weighted values ( w m a x / w m i n = 0.9/0.4), and the learning weighting value ( c 1 , c 2 = 2). The optimization toolbox in MATLAB is used to solve the PSO algorithm [30]. The effectiveness of the proposed strategy was verified using MATLAB running on MacBook Air (Early 2015), 1.6 GHz Dual-Core Intel Core i5, RAM 8 GB 1600 MHz DDR3. The verification findings for the suggested optimization approach are presented in the following section.

4. Results

In this study, a standalone microgrid with an industrial load, a WT, a PV and a BESS illustrated in Figure 5 is considered. The daily load demand, PV and WT generations with a 5-min resolution shown in Figure 6, tasken from [31], are used to verify the effectiveness of the proposed method. The industrial load has a peak power of 100 kW at 12:30 p.m. The PV generates power between 7:00 a.m. and 5:00 p.m. with a maximum value of 63 kW at noon. The power generated by wind turbines fluctuates around 35 kW throughout the day. The required power from the BESS is obtained based on the mismatch between the load demand and the total power generation. In Figure 7, the BESS’ discharging power is positive and the charging power is negative. The initial SOC is set at 50%. The standalone microgrid project life is set to 20 years. The case parameters used in the optimization are given in Table 1.
The rated energy of BESS E B , r a t e d   is 82 kWh ( q B = 1). For q B = 1, the BESS lifetime is 1.2 years, and the N P V C t o t a l throughout the project is US$204,436. When q B = 1, the SOC value is in the range of 20% to 80%. As the q B is increases, the battery size is increasing, and the deviation of the battery SOC from 50% is decreasing. However, the increased battery capacity results in higher total costs. Thus, the PSO optimization method is applied to find the optimal value of q B . The BESS SOCs for q B = 1, 1.761, and 5 are illustrated in Figure 8.
Figure 9 shows a stepped pattern of the relationship between N R B and q B . As the q B increases, the number of battery replacements will decrease. As the size of the battery increases, the SOH and the battery lifetime increase as well. The relationship between N R B and L B is shown as the purple line (saw-tooth pattern) and is used to obtain the N P V C r e using Equation (26).
Using the proposed method, the obtained value of US$200,653 was the lowest N P V C t o t a l at the q B of 1.761. This corresponds to the optimal capacity of the BESS of 144.4 kWh. The effects of q B on the C i n i t i a l , the N P V C o m and the N P V C r e are depicted in Figure 10. As the q B increases, the size of the battery and the cost increase, and the number of replacements decreases. If the battery size is small, the cost of battery replacement is low, but the cost of replacing the battery is high. As the q B varied from 1 to 5, in steps of 0.001, the C i n i t i a l increased from US$18,419 to US$78,724, and the N P V C o m increased from US$23.48 to US$117.43. In contrast, the C r e has downward trends because the q B can extend the BESS lifetime and the BESS replacement was reduced.
The effects of q B on the BESS lifetime and the number of BESS replacements are shown in Figure 11. The number of BESS replacements will decrease step by step with respect to the DOD and the size of the battery in relation to the q B factor. However, when calculating the BESS replacement cost in the low q B range, the replacement cost is low, and when q B is high, the cost of replacing the battery increases. The minimum saw-tooth pattern range is the point where the number of BESS replacements is reduced. The q B has a positive effect on the BESS lifetime. The BESS lifetime also increased from 1.2 years to 5.44 years when the q B increased from 1 to 5. The number of BESS replacements decreased because the number of replacements remains equivalent even when q B is adjusted slightly, such as when N R B = 8 and q B = 2.0 and 2.2.
The q B has an effect on the N P V C t o t a l of the microgrid, as shown in Figure 12. This curve is similar to a saw-tooth line due to the nonlinear characteristic of the replacement cost. When q B = 1.761 and q B = 2.0 are compared, the N P V C t o t a l of the q B = 1.761 is the lowest. Even though the C i n i t i a l and C o m for q B = 1.761 are greater than those for q B = 2.0, the C r e is lower because it can reduce the N R B from 9 to 8. This shows how the q B affects the N P V C t o t a l . The BESS lifetime and cost analysis versus some of the q B values are given in Table 2. Battery lifetime estimates show that increasing the size of the battery by the q B factor increases battery life because it reduces the number of cycles at a high DOD value.

5. Conclusions

This paper proposed a capacity optimization method for a BESS in a standalone microgrid while taking the BESS’ lifetime into account. The BESS’ capacity influenced the initial cost, operation and maintenance costs, and replacement cost. The case study demonstrated the efficacy of the proposed method. According to the PSO algorithm results, the optimal capacity of the BESS ( q B = 1.761, E B , r a t e d = 144.4 kWh, and N P V C t o t a l = US$200,653) has the lowest NPV of the total cost. According to the simulation results, the capacity and lifetime of the BESS were the major factors influencing the total cost of the system. Based on the results, the overall cost was reduced slightly. This is due to the small size of the test system, which caused the battery system to be small as well. However, if the proposed method is applied to a larger system, the total cost can be greatly reduced. Future research will concentrate on scaling up the proposed method and comparing it to other current battery technologies.

Author Contributions

Conceptualization, P.W., W.P., B.M. and B.H.; methodology, P.W. and T.B.; software P.W. and T.B.; validation, P.W., T.B., W.P., B.M. and B.H.; formal analysis, P.W., T.B., P.B., W.P., B.M. and B.H.; investigate T.B. and B.M.; data curation, T.B. and P.B.; writing—original draft preparation, P.W. and T.B.; writing—review and editing, T.B., P.B., W.P., B.M. and B.H.; supervision, B.M.; funding acquisition, B.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Suranaree University of Technology, Thailand (kittibundit-2559-3).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
BESSBattery Energy Storage System
c 1 , c 2 Acceleration factors
C i n i t i a l Initial cost
C o m Fixed operation and maintenance cost
C r e Replacement cost
C T F Cycles to failure
C t o t a l Total cost
dDiscount rate
DGsDistributed generation system
D O D Depth of discharge
E B Battery energy
E B , r a t e d Battery capacity
E B , O r a t e d Battery’s oversize energy
E t h r o u g h p u t Battery throughput corresponding to a specified DOD
E t h r o u g h p u t , a v g Average battery throughput
g b e s t The best global solution
G c The solar irradiance on the operating time
G S T C The solar irradiance on the standard test condition (STC) (1000w/ m 2 )
iParticle index
i t e m a x Maximum iteration
k c Temperature coefficient
kDiscrete time index
k B , o m The operation and maintenance cost set to 5% of the initial cost ($/kWh/year)
k B E , i n i t i a l The BESS initial cost per energy ($/kWh)
k B P , i n i t i a l The BESS initial cost per power ($/kW)
L B BESS lifetime
L l i f e Life loss
LPSPLoss of power supply probability
nNumber of particles in the swarm
N P V C Net present value of cost
N P V Net present value
N R B Number of BESS replacement throughput the project
p b e s t Personal pest solution
P B Output power limit of BESS
P B , c h , P B , d Charging and discharging power of BESS
P B , r e q u i r e d BESS’s required power
P B , r a t e d BESS’s rated power
P B , c h m a x , P B , d m a x Upper limits of charging and discharging power of BESS
P P V Output power of PV system
P r a t e d w t Rated output power of WT
P S T C Rated output power at standard test condition
PVPhotovoltaic
P W T Output power of WT
q B Adjusting factor
r 1 , r 2 Random number [0, 1]
S O C State of charge
T C The PV temperature on the operating time
U S W Uniform series presents the worth factor
vWind speed (m/s)
v i Velocity of the particle i
v c i , v r , v c o Cut-in speed, rated speed and cut-off speed of WT
wInertia factor
W S O C Weight factor
xPosition of particle i
β Battery’s self discharge rate
η B , c h , η B , d Charging and discharging efficiency of BESS

References

  1. Momoh, J.A. Smart Grid: Fundamentals of Design and Analysis; John Wiley & Sons: Hoboken, NJ, USA, 2012; Volume 63. [Google Scholar]
  2. Pinthurat, W.; Hredzak, B. Distributed Control Strategy of Single-Phase Battery Systems for Compensation of Unbalanced Active Powers in a Three-Phase Four-Wire Microgrid. Energies 2021, 14, 8287. [Google Scholar] [CrossRef]
  3. Annathurai, V.; Gan, C.K.; Ibrahim, K.A.; Baharin, K.A.; Ghani, M. A review on the impact of distributed energy resources uncertainty on distribution networks. Int. Rev. Electr. Eng. 2016, 11, 420. [Google Scholar] [CrossRef]
  4. Boonraksa, T.; Boonraksa, P.; Marungsri, B. Optimal capacitor location and sizing for reducing the power loss on the power distribution systems due to the dynamic load of the electric buses charging system using the artificial bee colony algorithm. J. Electr. Eng. Technol. 2021, 16, 1821–1831. [Google Scholar] [CrossRef]
  5. Hatziargyriou, N. Microgrids: Architectures and Control; John Wiley & Sons: Hoboken, NJ, USA, 2014. [Google Scholar]
  6. Chen, H.; Gao, L.; Zhang, Z.; Li, H. Optimal Energy Management Strategy for an Islanded Microgrid with Hybrid Energy Storage. J. Electr. Eng. Technol. 2021, 16, 1313–1325. [Google Scholar] [CrossRef]
  7. Chalise, S.; Sternhagen, J.; Hansen, T.M.; Tonkoski, R. Energy management of remote microgrids considering battery lifetime. Electr. J. 2016, 29, 1–10. [Google Scholar]
  8. Pinthurat, W.; Hredzak, B.; Konstantinou, G.; Fletcher, J. Techniques for compensation of unbalanced conditions in LV distribution networks with integrated renewable generation: An overview. Electr. Power Syst. Res. 2023, 214, 108932. [Google Scholar] [CrossRef]
  9. Gao, D.W. Energy Storage for Sustainable Microgrid; Academic Press: Cambridge, MA, USA, 2015. [Google Scholar]
  10. Elbaz, A.; Guneser, M.T. Multi-objective optimization method for proper configuration of grid-connected PV-wind hybrid system in terms of ecological effects, outlay, and reliability. J. Electr. Eng. Technol. 2021, 16, 771–782. [Google Scholar] [CrossRef]
  11. Sansa, I.; Belaaj, N.M.; Villafafilla, R. Optimal sizing design of an isolated micro grid keeping the longest battery lifetime. In Proceedings of the 2017 International Conference on Green Energy Conversion Systems (GECS), Hammamet, Tunisia, 23–25 March 2017; pp. 1–8. [Google Scholar]
  12. Falai, A.; Giuliacci, T.A.; Misul, D.A.; Anselma, P.G. Reducing the Computational Cost for Artificial Intelligence-Based Battery State-of-Health Estimation in Charging Events. Batteries 2022, 8, 209. [Google Scholar] [CrossRef]
  13. Iurilli, P.; Brivio, C.; Carrillo, R.E.; Wood, V. Physics-Based SoH Estimation for Li-Ion Cells. Batteries 2022, 8, 204. [Google Scholar] [CrossRef]
  14. Narayan, N.; Papakosta, T.; Vega-Garita, V.; Qin, Z.; Popovic-Gerber, J.; Bauer, P.; Zeman, M. Estimating battery lifetimes in Solar Home System design using a practical modelling methodology. Appl. Energy 2018, 228, 1629–1639. [Google Scholar] [CrossRef]
  15. Jenkins, D.P.; Fletcher, J.; Kane, D. Lifetime prediction and sizing of lead–acid batteries for microgeneration storage applications. IET Renew. Power Gener. 2008, 2, 191–200. [Google Scholar]
  16. Zhao, B.; Zhang, X.; Chen, J.; Wang, C.; Guo, L. Operation Optimization of Standalone Microgrids Considering Lifetime Characteristics of Battery Energy Storage System. IEEE Trans. Sustain. Energy 2013, 4, 934–943. [Google Scholar] [CrossRef]
  17. Anuphappharadorn, S.; Sukchai, S.; Sirisamphanwong, C.; Ketjoy, N. Comparison the economic analysis of the battery between lithium-ion and lead-acid in PV stand-alone application. Energy Procedia 2014, 56, 352–358. [Google Scholar] [CrossRef] [Green Version]
  18. Kebede, A.A.; Coosemans, T.; Messagie, M.; Jemal, T.; Behabtu, H.A.; Van Mierlo, J.; Berecibar, M. Techno-economic analysis of lithium-ion and lead-acid batteries in stationary energy storage application. J. Energy Storage 2021, 40, 102748. [Google Scholar] [CrossRef]
  19. Chen, S.X.; Gooi, H.B.; Wang, M.Q. Sizing of Energy Storage for Microgrids. IEEE Trans. Smart Grid 2012, 3, 142–151. [Google Scholar] [CrossRef]
  20. Bai, L.; Li, F.; Hu, Q.; Cui, H.; Fang, X. Application of battery-supercapacitor energy storage system for smoothing wind power output: An optimal coordinated control strategy. In Proceedings of the 2016 IEEE Power and Energy Society General Meeting (PESGM), Boston, MA, USA, 17–21 July 2016; pp. 1–5. [Google Scholar]
  21. Den Heeten, T.; Narayan, N.; Diehl, J.C.; Verschelling, J.; Silvester, S.; Popovic-Gerber, J.; Bauer, P.; Zeman, M. Understanding the present and the future electricity needs: Consequences for design of future Solar Home Systems for off-grid rural electrification. In Proceedings of the 2017 International Conference on the Domestic Use of Energy (DUE), Cape Town, South Africa, 4–5 April 2017; pp. 8–15. [Google Scholar]
  22. Narayan, N.; Popovic, J.; Diehl, J.C.; Silvester, S.; Bauer, P.; Zeman, M. Developing for developing nations: Exploring an affordable solar home system design. In Proceedings of the 2016 IEEE Global Humanitarian Technology Conference (GHTC), Seattle, WA, USA, 13–16 October 2016; pp. 474–480. [Google Scholar]
  23. Liu, C.; Wang, X.; Wu, X.; Guo, J. Economic scheduling model of microgrid considering the lifetime of batteries. IET Gener. Transm. Distrib. 2017, 11, 759–767. [Google Scholar] [CrossRef]
  24. Xtender, S. Sun Xtender® PVX-2580L Solar Battery Specifications. Available online: http://www.sunxtender.com/solarbattery.php?id=11 (accessed on 3 November 2022).
  25. Kennedy, J.; Eberhart, R. Particle swarm optimization. In Proceedings of the ICNN’95-International Conference on Neural Networks, Perth, Australia, 27 November–1 December 1995; Volume 4, pp. 1942–1948. [Google Scholar]
  26. Liu, H.; Ji, Y.; Zhuang, H.; Wu, H. Multi-objective dynamic economic dispatch of microgrid systems including vehicle-to-grid. Energies 2015, 8, 4476–4495. [Google Scholar] [CrossRef] [Green Version]
  27. Sreenivasan, G. Solution of dynamic economic load dispatch (DELD) problem with valve point loading effects and ramp rate limits using PSO. Int. J. Electr. Comput. Eng. 2011, 1, 59. [Google Scholar] [CrossRef] [Green Version]
  28. Dhayalini, K.; Sathiyamoorthy, S.; Rajan, C.A. Particle Swarm Optimization Technique for the Coordination of Optimal Wind and Thermal Generation Dispatch. Int. Rev. Electr. Eng. 2013, 8, 1843–1849. [Google Scholar]
  29. Shi, Y.; Eberhart, R.C. Parameter selection in particle swarm optimization. In Proceedings of the International Conference on Evolutionary Programming, San Diego, CA, USA, 25–27 March 1998; Springer: Berlin/Heidelberg, Germany, 1998; pp. 591–600. [Google Scholar]
  30. Birge, B. PSOt-a particle swarm optimization toolbox for use with Matlab. In Proceedings of the 2003 IEEE Swarm Intelligence Symposium. SIS’03 (Cat. No. 03EX706), Indianapolis, IN, USA, 26 April 2003; pp. 182–186. [Google Scholar]
  31. Homer Software. Available online: https://www.homerenergy.com/products/pro/index.html?fbclid=IwAR3mBIFJhn92zNB5JFMnuw96avfNlx1WH_3Px5DsVxO6RlxSclq_kuxcwU8 (accessed on 20 January 2018).
  32. Nguyen, C.L.; Lee, H.H. An optimal hybrid supercapacitor and battery energy storage system in wind power application. In Proceedings of the IECON 2015-41st Annual Conference of the IEEE Industrial Electronics Society, Yokohama, Japan, 9–12 November 2015; pp. 003010–003015. [Google Scholar]
  33. Microgrid, S. The Advantage, and Cost of Lead-Acid and Li-Ion. Available online: https://medium.com/solar-microgrid/battery-showdown-lead-acid-vs-lithium-ion-1d37a1998287 (accessed on 3 November 2022).
  34. Sivaram, A. Cost Comparison of Battery Technology. Available online: https://saurorja.org/2011/08/30/lead-acid-is-the-cheapest-battery-conditions-apply/ (accessed on 3 November 2022).
  35. Ren, M.; Mitchell, C.R.; Mo, W. Managing residential solar photovoltaic-battery systems for grid and life cycle economic and environmental co-benefits under time-of-use rate design. Conserv. Recycl. 2021, 169, 105527. [Google Scholar] [CrossRef]
Figure 1. The relationship between CTF and DOD.
Figure 1. The relationship between CTF and DOD.
Batteries 09 00076 g001
Figure 2. The relationship between energy throughput and DOD.
Figure 2. The relationship between energy throughput and DOD.
Batteries 09 00076 g002
Figure 3. The relationship between SOC and the SOC weighting factor W S O C .
Figure 3. The relationship between SOC and the SOC weighting factor W S O C .
Batteries 09 00076 g003
Figure 4. Flowchart for determination of the optimal BESS capacity. The q B value is changed from 1 to 5 in 0.001 increments.
Figure 4. Flowchart for determination of the optimal BESS capacity. The q B value is changed from 1 to 5 in 0.001 increments.
Batteries 09 00076 g004
Figure 5. Schematic diagram of the standalone microgrid.
Figure 5. Schematic diagram of the standalone microgrid.
Batteries 09 00076 g005
Figure 6. Industrial load, PV source and WT powers taken from [31].
Figure 6. Industrial load, PV source and WT powers taken from [31].
Batteries 09 00076 g006
Figure 7. BESS−power for the case shown in Figure 6.
Figure 7. BESS−power for the case shown in Figure 6.
Batteries 09 00076 g007
Figure 8. SOCs of BESS for q B values of 1, 1.761 and 5.
Figure 8. SOCs of BESS for q B values of 1, 1.761 and 5.
Batteries 09 00076 g008
Figure 9. Relationships between N R B and L B for different values of q B .
Figure 9. Relationships between N R B and L B for different values of q B .
Batteries 09 00076 g009
Figure 10. The costs for different values of q B .
Figure 10. The costs for different values of q B .
Batteries 09 00076 g010
Figure 11. The BESS lifetime and the number of replacements for different values of q B .
Figure 11. The BESS lifetime and the number of replacements for different values of q B .
Batteries 09 00076 g011
Figure 12. Net present value of the system total cost for different values of q B .
Figure 12. Net present value of the system total cost for different values of q B .
Batteries 09 00076 g012
Table 1. Case study parameters.
Table 1. Case study parameters.
ParameterVariableUnitValue
Project life L p years20
BESS Calendar life L B years10
BESS SOC limits S O C m i n - S O C m a x %20–80
Charge/discharge efficiency η B , c h / η B , d %90/90
Initial cost per energy [32,33,34] k B E , i n i t i a l $/kWh183.86
Initial cost per power k B P , i n i t i a l $/kW183.86
Operation & maintenance cost k B , o m $/kWh/year9.19
Discount rate [35]d%5
Table 2. Simulation results for some of the q B values.
Table 2. Simulation results for some of the q B values.
q B L B NR B C initial ($) NPVC om ($) NPVC re ($) NPVC total ($)
1.01.201618,41923.48185,992204,436
1.51.721125,95835.23178,562204,557
1.7612.00929,89641.36170,468200,653
2.02.25833,49746.97168.626202,172
2.22.46836,51351.67176,524213,088
3.03.31648,57370.46171,838220,482
4.04.38463,64993.94153,539217,284
5.05.44378,724117.43142,152220,994
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

Wongdet, P.; Boonraksa, T.; Boonraksa, P.; Pinthurat, W.; Marungsri, B.; Hredzak, B. Optimal Capacity and Cost Analysis of Battery Energy Storage System in Standalone Microgrid Considering Battery Lifetime. Batteries 2023, 9, 76. https://doi.org/10.3390/batteries9020076

AMA Style

Wongdet P, Boonraksa T, Boonraksa P, Pinthurat W, Marungsri B, Hredzak B. Optimal Capacity and Cost Analysis of Battery Energy Storage System in Standalone Microgrid Considering Battery Lifetime. Batteries. 2023; 9(2):76. https://doi.org/10.3390/batteries9020076

Chicago/Turabian Style

Wongdet, Pinit, Terapong Boonraksa, Promphak Boonraksa, Watcharakorn Pinthurat, Boonruang Marungsri, and Branislav Hredzak. 2023. "Optimal Capacity and Cost Analysis of Battery Energy Storage System in Standalone Microgrid Considering Battery Lifetime" Batteries 9, no. 2: 76. https://doi.org/10.3390/batteries9020076

APA Style

Wongdet, P., Boonraksa, T., Boonraksa, P., Pinthurat, W., Marungsri, B., & Hredzak, B. (2023). Optimal Capacity and Cost Analysis of Battery Energy Storage System in Standalone Microgrid Considering Battery Lifetime. Batteries, 9(2), 76. https://doi.org/10.3390/batteries9020076

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