Next Article in Journal
Sustainability Considerations in Digital Fabrication Design Education
Previous Article in Journal
Thermoeconomic Optimization of a Polygeneration System Based on a Solar-Assisted Desiccant Cooling
Previous Article in Special Issue
Linear Permanent Magnet Vernier Generators for Wave Energy Applications: Analysis, Challenges, and Opportunities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Composite Demand-Based Energy Storage Sizing for an Isolated Microgrid System

Department of Electrical Engineering, Faculty of Engineering, Islamic University of Madinah, Madinah 42351, Saudi Arabia
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(2), 1517; https://doi.org/10.3390/su15021517
Submission received: 14 December 2022 / Revised: 7 January 2023 / Accepted: 9 January 2023 / Published: 12 January 2023
(This article belongs to the Special Issue Sustainable Development and Optimisation of Energy Systems)

Abstract

:
This paper presents a comprehensive model for optimal energy storage system (ESS) design for an isolated microgrid. The model presented is a mixed integer linear program (MILP) that considers seasonal varying generation (VG) demand, more specifically seasonal solar cell generator (SCG) demand, SCG maintenance (failure and restoration) rates, and practical operation of conventional generation (CG) while satisfying the required demand and reserve. The model is based on unit commitment (UC) to simulate real operations and physical constraints of CG units, the power balance, and reserve requirements. The objective function aims at minimizing the associated cost of CG, namely, production (fuel), costs of startup and shutdown procedures, and the investment cost of power and energy. The proposed model is assessed on a case study system consisting of multiple SCGs in addition to CG to meet a specific demand. The proposed model showed that the ESS sizing when considering Li-Ion technology and a SCG penetration of 25% was on average approximately 3 MWh and 1.70 MW. Meeting the demand and reserve requirements were the two major constraints when determining the optimal ESS sizing. Moreover, introducing the ESS substantially reduced the operating cost of the system.

1. Introduction

A microgrid is defined as “a group of interconnected demands and distributed energy resources with defined electrical boundaries forming a local electric power system at distribution voltage levels, that acts as a single controllable entity and is able to operate in either grid-connected or island mode” [1]. Accordingly, the microgrid can be categorized as either grid-connected or isolated/islanded. Both types can play essential roles in (a) enhancing the stability of power networks, (b) reducing power flow losses in transmission and distribution networks, (c) reducing harmful emissions when integrating renewable energy (RE), (d) providing independent, fully or partially, electric power supply, (e) enhancing the power system reliability and energy quality, and (f) providing back-up supply during outages of the main grid [2]. Microgrid services may be delivered when there are adequate generation resources. These resources may be conventional and/or, more preferably, renewable. RE includes photovoltaic systems (PVs) or solar cell generators (SCGs), wind turbine systems (WTSs), etc. Even though RE resources bring several benefits to power systems, including microgrids, they are characterized by intermittency, which may affect the reliability and economic operation of the power system. Hence, RE is referred to as having variable generation (VG). When planning for an increase in the penetration of VG, it is vital to integrate energy storage systems (ESSs). In recent years, several energy storage technologies have been developed, including electrochemical batteries, compressed air, flywheels, etc. Depending on the technique employed, the ESSs are put forward in different forms and specifications [3,4,5]. Most importantly, ESSs have different capital power and energy costs, which are important factors in determining the required size for use in microgrids. There has been a plethora of research work to determine the optimal sizing of microgrid ESSs depending on whether the microgrid is operating as grid-connected or islanded. Although this work focuses on ESS optimal sizing for islanded microgrids with VG, ESS optimal sizing for large power systems and microgrids (both islanded and grid-connected) involves similar model structures, so those are investigated. The optimal sizing of ESSs is a complex problem as it encompasses several factors and constraints. The authors in [6] proposed a technique for sizing ESSs in the context of WTS generation based on probability theory. A Markov-chain-based stochastic model was used to detect the mismatch between demand and generation. In [7], the authors proposed the optimal sizing of the ESS in a hybrid microgrid with VG using a Monte Carlo simulation (MCS) and considered reliability as a constraint while searching for a pattern optimization model. In [8], the mixed integer linear programming (MILP) approach was used to design a hybrid power system. The method was based on the average time of repair and failure rates of the WTSs/SCGs and relied principally on the MCS to produce their chronological state samples. The non-utilized VG models were highly penalized to avoid possible curtailment, whereas the VG models that either charged the ESSs or directly met the demand were maximized. This model used the deterministic demand forecast. The authors of [9] presented an optimal sizing model for the ESS in the context of WTS and SCG generation and considered the correlation between demand and VG output. All these works proposed general techniques for optimal ESS sizing, rather than specifically for microgrid applications, although they could be used for microgrid applications.
There have been a few research works focused on the optimal sizing of ESSs for islanded microgrids. For example, in [10], the authors presented a long-term methodology for the optimal sizing and planning of the life cycle of the ESS in islanded microgrids, consisting of hybrid SCG-WTS-diesel generators implemented through a multi-scaled decision-making process to meet the demand, while considering capacity fading of ESS modules. Meanwhile, the authors in [11] presented a method to determine the optimal sizing of ESSs in an islanded microgrid based on a two-step cost. The first step involved a unit commitment (UC) to obtain the operation of the microgrid, while the second step used the convex optimization principle, considering different physical and operational constraints, to determine the optimal size of the ESS. Using the incremental cost method based on economic dispatch, the authors of [12] concluded that there was a near-linear relationship between the optimal ESS sizing and ESS efficiency, and the integration of ESS led to a significant reduction in the operation cost. Elsewhere, the authors of [13] proposed an optimization strategy to determine the optimal size of the ESS. A two-layer optimal sizing method combining an iterative method with dynamic programming (DP) was used. Taking a different approach, the authors in [14] used the mixed integer programming (MIP) technique for optimal ESS sizing and considered the reliability of the system as a constraint. The study concluded that a larger size of ESS may have greater costs than benefits for microgrids. Meanwhile, the authors of [15] proposed a size-optimization method based on multi-objective grey wolf (IMOGWO) for a hybrid energy system. The proposed method was aimed at determining the optimal sizes of the different components of systems, including ESSs, while minimizing the annual cost and loss of power supply. In [16], a two-archive many-objective evolutionary algorithm (TA-MaEA) was proposed for the optimal sizing of hybrid microgrids. The aim of this study was to minimize costs, loss of power supply, and emissions, and to maintain the power balance. Elsewhere, the authors in [17] developed a multi-objective optimization model to solve the system-sizing problem for a microgrid while considering the rate of battery degradation, monetary incentives, and PV system azimuth angle. A methodology for optimal sizing of islanded microgrids, including SCGs and ESSs, was then presented in [18] based on a deterministic cost model and incorporating local tax benefits, technical constraints of ESS, and reliability. All aforementioned works proposed different techniques for the optimal sizing of ESSs for microgrid applications. However, the seasonal correlation between demand and VG was only considered by [17], and the failure rate of SCG units was only addressed in [7,8,9,10,12,14]. Furthermore, reserve requirement was used in some of the proposed models, namely [8,9,10,11].
The present research is aims to fill the gaps and address all shortcomings of the previous works related to ESS sizing for microgrid applications. The proposed ESS optimal-sizing model considers: (a) SCGs’ unavailability based on their forced outage rates (FORs), (b) the correlation of seasonal solar radiation with demand, (c) the operational constraints of the CG units, and (d) demand and requirements. Since the demand and SCGs are the main factors for the sizing of the ESS, the seasonal variations SCG and demand will represent the required ESS size accurately in the long term for an isolated microgrid. The approach considers the optimized sizing of the ESS as a MILP model.
For a specific region, the historical solar radiation data are divided into four ranges, one for each yearly season. Then, the expected SCG generation is computed on (a) the generation model, and (b) the availability model, i.e., the collective SCG availability PDFs based on their FORs. Subsequently, an important step is to compute the correlation between each demand and expected output pair for the SGSs. The main aim of optimizing the size of the ESS is to reduce the associated costs with the CG units, namely production (fuel), start-up and shut-down costs, and capital investment costs in the energy and power for the ESS. Figure 1 depicts a flowchart of the proposed model.

2. Research Methodology and Modeling

2.1. Problem Statement

Assuming that the following information is given:
  • Some number (G) of CG units with known specifications.
  • Energy and power capital cost of the ESS, and charging and discharging efficiencies.
  • Historical solar radiation data (Gh) divided into four groups: summer (GSu), fall (GFa), winter (GWi), and spring (GSp).
  • Historical demand data (Dh) for the specific region and season.
  • A solar farm (SF) consists of NSCG SCGs, for which the specifications and FORs are given.
The expected SF power (pSF) values are obtained by interrelating the power with the availability of the SF’s probabilistic model. Usually, at each instant in time, composite demand (CD) is taken as the difference between demand and VG. After computing the CD values for the whole year, these are then used to determine the optimal size of the ESS based on the proposed sizing model, the solution of which gives the desired ESS charging and discharging profile.

2.2. Computation of the Expected Output Power of the SF and Composite Demand

2.2.1. Computation of the Expected Output Power of the SF

Once the historical solar radiation, Gh, is partitioned according to seasons, it may be regenerated as power using the generation model of the SCG [19], as in (1):
p S C G = { p S C G . r a t e d ( G t 2 G s t d R C   )                 G t [ 0 , R C ) p S C G . r a t e d   ( G t G s t d   )                       G t [ R C , G s t d ]     p S C G . r a t e d                           G t ( G s t d , )  
where Gt is the solar radiation at time t (W/m2), PSCG.rated is the SCG unit rated power (MW), Gstd is the solar radiation in the standard environment (W/m2), and Rc is a certain radiation point set usually at 150 W/m2.
The SCG is represented as a space–time Markov model that is equal to either the SCG rated power, pSCG.rated, or to zero, as depicted in Figure 2. The chance of the ith SCG being unavailable is referred to as the forced outage rate (qSCGi) and is calculated as in Equation (2) [20]:
q S C G i = M T T R S C G i M T T F S C G i + M T T R S C G i
where MTTRSCG/MTTFSCG is the mean time value to maintain the SCG. Both MTTRSCG and MTTFSCG are assumed to follow an exponential distribution. For simplicity, qSCG is taken as identical for all SCGs.
The total law of probability is used and is mathematically represented by Equation (3) and illustrated in Figure 3:
ρ S G   ( p S C G ) = ρ S G   ( p S C G |   S C G   i s   U P ) ( 1 - q S C G ) + ρ S G (   p S C G |   S C G   is D O W N ) ( q S C G )
where ρSG(pSCG) is the PDF of the SCG power output.
Finally, the probability mass function (PMF) of the SF (ρSFA (cSF)) being available can be obtained by the binomial distribution, as given by Equation (4) and shown in Figure 4:
ρ S F A ( c S F ) = r = 0 N S C G ( N S C G r ) ( 1 q S C G ) r ( q S C G ) N S C G r δ ( c S F r   p S C G . r a t e d )
where r is an index for the available SCG units in the SF and NSCG is the number of SCGs in the SF. Thus, after convolving the seasonal partitioned Gh and ρSFA(cSF), the expected SF output (PSF) is as follows:
p S F = ρ S F A ( c S F ) * G S e a s o n Season = Wi ,   Sp ,   Su ,   and   Fa
where * represents the convolution operator.

2.2.2. Composite Demand PDF Computation

The composite demand, CD, is the demand seen by the CG. In other words, CD is the remaining demand after deploying pSF. Hence, seasonal pSF and demand are known, and CD is:
CD = DpSF
Then, the seasonal PDF/CDF (ρCD (cd)/FCD (cd)) of CD is computed and samples will be inputted to the energy sizing model. Note that CD could be negative when there is higher pSF than demand, e.g., when the demand is low during off-peak hours.

2.3. Formulation of the Optimized Sizing of the Energy Storage

The model for optimally sizing the ESS is based on the models proposed in [9] and [21,22,23]. A detailed description is given in the following section. First, it is important to introduce the variables and parameters that appear in the constraints, as listed in Table 1, for both CG and ESS, where T is the simulation time (h), and Δt is the time step (h).
The constraints that relate to the operation of the CG are as follows. The first constraint in (7) determines the CG unit’s on/off status’ variations between time steps. This constraint is important for the determination of the minimum up and down times in addition to the committed CG units. The variables sg,t and zg,t determine the start-up shutdown status and are obtained from determining the binary variables xg,t and xg,t−1. The second constraint in (8) is to determine the limit of the minimum CG. The constraint in (9) sets the maximum limit of the power and reserve provision. The constraint in (10) is the minimum up time, UTg, and (11) is the minimum down time constraint, DTg. The CG unit must be in the on/off state for a time equal to UTg/DTg once it is on/off. The other constraints in (12) and (13) for CG define the limits of the power ramp up/down of a CG unit. These constraints are listed below:
x g , t x g , t 1 = s g , t z g , t g G , t T \ t = 1
p g , t P G m i n   x g , t g G ,   t T  
p g , t p ¯ g , t   P G m a x x g , t + (   S D g P G m a x )   z g , t + 1 g G ,   t T
where p ¯ g , t = p g , t + r g , t .
i = t U T g + 1 t s g , i x g , t t T ,   g G
i = t D T g + 1 t s g , i 1 x g , t D T g t T ,   g G  
p ¯ g , t p g , t 1   S U g   s g , t + R U g   x g , t 1
p g , t 1 p g , t   S D g   z g , t + R D g   x g , t
In this model, the CG units are responsible for providing the reserve, r g , and the ESS is responsible for the reserve (rESS_DN,t\rESS_UP,t). Both the ESS and CG are in the microgrid to satisfy the demand and reserve requirements at each instant of time t. The required reserve for the system, R, is derived as a percentage of the annual or system peak demand, e.g., 10% of the system peak. The ESS reserve provision will be discussed when introducing the ESS constraints. The constraint in (14) ensures that the CG units and the ESS can satisfy both the demand and reserve requirements:
g G p ¯ g , t + p ¯ E S S _ U P , t C D t + p ¯ E S S _ D N , t + R t T
The constraint in (15), related to the power balance, ensures that the available generation (CG and discharging ESS power) meets the demand (CD and charging ESS power) at any given time t.
g G p g , t + p d i s , t = C D t + p c h , t t T  
CDt sampling depends on the season and is calculated by producing random numbers (~unif(0,1)), at any instant of time, t. Then, the inverse transform method (ITM) is applied using the ρCD(cd)/FCD (cd) corresponding to the season in which t falls, as shown in (16):
C D t = F C D 1 ( U 1 , t ) t T  
This process will be repeated for all seasons.
The ESS set of constraints describes the dynamics of the ESS in addition to accounting for reserve provision capability. These operating constraints are the main factors for obtaining the optimal size of the ESS. They determine the ESS charge/discharge schedule and set the power and energy limits. The constraint given in (17) determines the ESS state of charge (SOC) while the constraint in (18) sets the maximum and minimum limits of ESS energy. Then, the constraints in (19) and (20) are the charging and the down reserve and discharging and up reserve power limits, respectively. The constraints in (21) and (22), meanwhile, ensure that no charging and discharging of the ESS take place simultaneously and that there is no simultaneous up and down reserve provision. Finally, the constraints (23) and (24) make sure that the SOC is not exceeded at time t when the ESS provides power and a reserve.
E t + 1 = E t + η c h p c h , t p d i s , t η d i s t T  
0 E t E E S S t T  
0 p c h , t p ¯ E S S _ D N , t P E S S t T  
0 p d i s , t p ¯ E S S _ U P , t P E S S t T  
0 p ¯ E S S _ U P , t α t M t T  
0 p ¯ E S S _ D N , t ( 1 α t ) M t T  
p ¯ E S S _ U P , t E t Δ t t T  
E E S S E t Δ t p ¯ E S S _ D N , t 0 t T  
where p ¯ E S S _ U P , t = p d i s , t + r E S S _ U P , t   and   p ¯ E S S _ D N , t = p c h , t + r E S S _ D N , t .
The linear objective function will minimize the CG production and the capital costs for power and energy for the ESS. It is given by:
m i n   t T g G ( y g , t + S U c o s t g , t + z g , t S D c o s t g ) + d E E S S + e P E S S
where yg,t is the gth CG linearized production cost ($), and SUcost/SDcost are the linearized startup/shutdown costs, as introduced earlier. The production cost of a CG is nonlinear and has the following form:
a g P g , 2 t + b g P g , t + c g
where the parameters ag, bg, and cg are taken as the coefficients cost of the gth CG. For simplicity, a linear model is used in this work; however, when using the quadratic cost function, it can be linearized using piecewise segments, as shown in [20] and explained in [9]. In the same manner, the startup cost will be linearized by an approximate staircase function [9,24,25,26].

3. Results

The hourly solar radiation data used in this study were collected in the City of Madinah in Saudi Arabia [27]. Meanwhile, the load demand data consisting of different building types, e.g., a school and apartment building, were taken from [28]. The demand is fully described below, and was chosen to represent different residential, commercial, health, and education services. A full description of the test system is as follows:
  • The parameter values of the CG units are given in Table 2 [29]. Note that all the parameters of the two CG units are identical except for the cost coefficient. The cost of CG 1 is less than that of CG 2.
  • An SF has NSCG identical SCGs, the specifications for which are shown in Table 3 [30,31].
  • The ESS specifications are listed in Table 4 [32].
  • Demand data are taken from the Open Energy Data Initiative (OEDI) [28], as shown in Table 5. Note that the number of units for each demand is assumed.
Incidentally, the VG penetration level is taken as 25% of the total installed capacity (10 MW). NSCG depends on the penetration level of the VG (for 25% penetration level, NSCG = 50 SCGs). The ESS technology used in the system is a grid-scale Li-Ion battery. The parameters for the simulation in this research are a time of 8760 h, an Δt of 1 h, and an ESS life expectancy of 20 years. e is the energy capital cost, and d is the cost of the power over a period T with a 5% reduction rate.

3.1. Power Outputs of the Expected SF and Composite Demand PDF Results

Samples of the seasonally partitioned historical solar radiation data are depicted in Figure 5. By using the SCG generation models introduced earlier and the historical solar radiation data, the expected seasonal power output of the SF was computed and then convolved with ρSFA(cSF), as described in Equation (5) and shown in Figure 6. Figure 7 shows the seasonal CD calculated using Equation (6). The seasonal CD was then converted to PDF as illustrated in Figure 8. Hourly annual samples were taken from these PDFs, with a total of 8760 samples for each season, which are represented by 2190 samples. The seasonal PDFs highlight the seasonal demand and the expected output variabilities of SFs. Sampling each season’s PDFs detects these variabilities and thus might give an accurate definition of the sizing problem of the ESS. Figure 8 shows each season ρCD (cd). More fluctuations are observed during winter, spring, and fall, whereas there are fewer fluctuations in summer. This may be explained by examining Figure 8 again, in which the summer CD PDF has the highest demand values and resembles a normal distribution, whereas in other seasons, there is a smaller probability of high demand.

3.2. ESS Sizing and CG Operational Cost Results

The ESS sizing model was run multiple times, and Table 6 shows the average ESS sizing and CG operational costs results, with the average ± standard deviation. The average EESS was found to be approximately 3 ± 0.30 MWh, whereas the PESS was 1.70 ± 0.40 MW. These resulted in an ESS capital of $320,514 ± $34,950, which represents about of 40% of the total cost, i.e., $815,4645 ± $22,259. The cost for CG and PV (no ESS) was $718,638 ± $1042. This cost represents the CG cost only. Comparing the CG costs for this case and after the addition of ESS clearly shows that adding the ESS decreased the CG cost (by 31%, to $494,951 ± $20,857). Figure 9 shows the seasonal variation of the ESS SOC.
It is worth noting that the ESS, in addition to CG, could result in a higher ESS sizing and hence a higher investment cost by providing the reserve to meet the reserve requirement. Figure 10 shows the seasonal reserve provided by the ESS and CG. It can be noted that the ESS provided a substantial percentage of the reserve.

4. Conclusions

A mixed integer linear programming (MILP) model for optimal sizing of energy storage for islanded microgrids was presented. The model is a unit commitment model and considers the following: unavailability of solar cell generators, seasonal correlation of solar radiation with demand, operational constraints of CG, and demand and reserve requirements. When running the proposed model, we showed that the ESS sizing, when considering Li-Ion technology and an SCG penetration of 25%, was on average approximately 3 MWh and 1.7 MW. The case study also showed that the CG cost was significantly reduced in the presence of ESS when compared to a scenario without ESS. This reduction could be attributed to a change in the operation of CG when an ESS was present, resulting in lower costs.
While this model considers several aspects, it has some limitations that could be addressed in future work. These include ESS degradation over time, in addition to factors that affect SCG efficiency. In future work, we intend to use the proposed approach for more accurate modeling of long-term ESS operation. The effect of SCG derating factors, such as soiling, conversion, temperature-related factors, and others affecting ESS sizing, are also important and will be considered in future work.

Author Contributions

Conceptualization, A.A. (Abdullah Alamri); methodology, A.A. (Abdullah Alamri); software, A.A. (Abdullah Alamri); validation, A.A. (Abdullah Alamri) and A.A. (Abdulrahman AlKassem); formal analysis, A.A. (Abdullah Alamri); investigation, A.A. (Abdullah Alamri) and A.D.; resources, A.A. (Abdullah Alamri) and A.D.; data curation, A.A. (Abdulrahman AlKassem); writing—original draft preparation, A.A. (Abdullah Alamri); writing—review and editing, A.A. (Abdullah Alamri), A.D. and A.A. (Abdulrahman AlKassem); visualization, A.A. (Abdullah Alamri) and A.A. (Abdulrahman AlKassem); supervision, A.A. (Abdullah Alamri) and A.D.; project administration, A.A. (Abdullah Alamri) and A.A. (Abdulrahman AlKassem); funding acquisition, A.D. All authors have read and agreed to the published version of the manuscript.

Funding

The research was funded by the Islamic University of Madinah (IUM) in Saudi Arabia through the Deanship of Scientific Research (DSR) under grant no. 656 in November 2020 through the “Tamayyuz 2” Financial Program.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. IEEE Std 2030.7-2017; IEEE Standard for the Specification of Microgrid Controllers. IEEE: New York, NY, USA, 2018; pp. 1–43.
  2. Choudhury, S. A comprehensive review on issues, investigations, control and protection trends, technical challenges and future directions for Microgrid technology. Int. Trans. Electr. Energy Syst. 2020, 30, e12446. [Google Scholar] [CrossRef]
  3. The, U.S. Energy Information Administration (EIA). U.S. Battery Storage Market Trends; The U.S. Energy Information Administration (EIA): Washington, DC, USA, 2018. [Google Scholar]
  4. Parhizi, S.; Lotfi, H.; Khodaei, A.; Bahramirad, S. State of the Art in Research on Microgrids: A Review. IEEE Access 2015, 3, 890–925. [Google Scholar] [CrossRef]
  5. Ghiassi-Farrokhfal, Y.; Rosenberg, C.; Keshav, S.; Adjaho, M.-B. Joint Optimal Design and Operation of Hybrid Energy Storage Systems. IEEE J. Sel. Areas Commun. 2016, 34, 639–650. [Google Scholar] [CrossRef] [Green Version]
  6. Dong, J.; Gao, F.; Guan, X.; Zhai, Q.; Wu, J. Storage Sizing with Peak-Shaving Policy for Wind Farm Based on Cyclic Markov Chain Model. IEEE Trans. Sustain. Energy 2016, 8, 978–989. [Google Scholar] [CrossRef]
  7. Arabali, A.; Ghofrani, M.; Etezadi-Amoli, M.; Fadali, M.S. Stochastic Performance Assessment and Sizing for a Hybrid Power System of Solar/Wind/Energy Storage. IEEE Trans. Sustain. Energy 2013, 5, 363–371. [Google Scholar] [CrossRef]
  8. Alamri, A.; AlOwaifeer, M.; Meliopoulos, A.S.; Cokkinides, G.J. Energy Storage Sizing and Reliability Assessment for Power Systems with Variable Generation. In Proceedings of the 2019 IEEE Milan PowerTech, Milano, Italy, 23–27 June 2019; pp. 1–6. [Google Scholar] [CrossRef]
  9. Alamri, A.; Alowaifeer, M.; Meliopoulos, A.P.S. Energy Storage Sizing and Probabilistic Reliability Assessment for Power Systems Based on Composite Demand. IEEE Trans. Power Syst. 2021, 37, 106–117. [Google Scholar] [CrossRef]
  10. Zhang, Y.; Wang, J.; Berizzi, A.; Cao, X. Life cycle planning of battery energy storage system in off-grid wind_solar_diesel microgrid. IET Gener. Transmiss. Distrib. 2018, 12, 4451–4461. [Google Scholar] [CrossRef]
  11. Zolfaghari, M.; Ghaffarzadeh, N.; Ardakani, A.J. Optimal sizing of battery energy storage systems in off-grid micro grids using convex optimization. J. Energy Storage 2019, 23, 44–56. [Google Scholar] [CrossRef]
  12. Hesaroor, K.; Das, D. Optimal sizing of energy storage system in islanded microgrid using incremental cost approach. J. Energy Storage 2019, 24, 100768. [Google Scholar] [CrossRef]
  13. Pham, C.M.; Tran, Q.T.; Bacha, S.; Hably, A.; Nugoc, A.L. Optimal sizing of battery energy storage system for an island microgrid. In Proceedings of the IECON 2018—44th Annual Conference of the IEEE Industrial Electronics Society, Washington, DC, USA, 21–23 October 2018. [Google Scholar]
  14. Bahramirad, S.; Reder, W.; Khodaei, A. Reliability-Constrained Optimal Sizing of Energy Storage System in a Microgrid. IEEE Trans. Smart Grid 2012, 3, 2056–2062. [Google Scholar] [CrossRef]
  15. Zhu, W.; Guo, J.; Zhao, G.; Zeng, B. Optimal Sizing of an Island Hybrid Microgrid Based on Improved Multi-Objective Grey Wolf Optimizer. Processes 2020, 8, 1581. [Google Scholar] [CrossRef]
  16. Cao, B.; Dong, W.; Lv, Z.; Gu, Y.; Singh, S.; Kumar, P. Hybrid Microgrid Many-Objective Sizing Optimization With Fuzzy Decision. IEEE Trans. Fuzzy Syst. 2020, 28, 2702–2710. [Google Scholar] [CrossRef]
  17. Bandyopadhyay, S.; Mouli, G.R.C.; Qin, Z.; Elizondo, L.R.; Bauer, P. Techno-Economical Model Based Optimal Sizing of PV-Battery Systems for Microgrids. IEEE Trans. Sustain. Energy 2019, 11, 1657–1668. [Google Scholar] [CrossRef]
  18. Ropero-Castaño, W.; Muñoz-Galeano, N.; Caicedo-Bravo, E.F.; Maya-Duque, P.; López-Lezama, J.M. Sizing Assessment of Islanded Microgrids Considering Total Investment Cost and Tax Benefits in Colombia. Energies 2022, 15, 5161. [Google Scholar] [CrossRef]
  19. Park, J.; Liang, W.; Choi, J.; El-Keib, A.A.; Shahidehpour, M.; Billinton, R. A probabilistic reliability evaluation of a power system including Solar/Photovoltaic cell generator. In Proceedings of the 2009 IEEE Power & Energy Society General Meeting, Calgary, AB, Canada, 26–30 July 2009; pp. 1–6. [Google Scholar] [CrossRef]
  20. Allan, R.; Billinton, R. Reliability Evaluation of Power Systems, 1st ed.; Plenum: New York, NY, USA, 1996. [Google Scholar]
  21. Ostrowski, J.; Anjos, M.F.; Vannelli, A. Tight Mixed Integer Linear Programming Formulations for the Unit Commitment Problem. IEEE Trans. Power Syst. 2011, 27, 39–46. [Google Scholar] [CrossRef]
  22. Atakan, S.; Lulli, G.; Sen, S. A State Transition MIP Formulation for the Unit Commitment Problem. IEEE Trans. Power Syst. 2017, 33, 736–748. [Google Scholar] [CrossRef] [Green Version]
  23. Baker, K.; Hug, G.; Li, X. Energy Storage Sizing Taking into Account Forecast Uncertainties and Receding Horizon Operation. IEEE Trans. Sustain. Energy 2016, 8, 331–340. [Google Scholar] [CrossRef]
  24. Zimmerman, R.D.; Murillo-Sanchez, C.; Thomas, R.J. MATPOWER’s extensible optimal power flow architecture. In Proceedings of the 2009 IEEE Power & Energy Society General Meeting, Calgary, AB, Canada, 26–30 July 2009; pp. 1–7. [Google Scholar]
  25. Wood, A.J.; Wollenberg, B.F. Power Generation, Operation, and Control, 2nd ed.; Wiley: New York, NY, USA, 1996. [Google Scholar]
  26. Carrion, M.; Arroyo, J. A Computationally Efficient Mixed-Integer Linear Formulation for the Thermal Unit Commitment Problem. IEEE Trans. Power Syst. 2006, 21, 1371–1378. [Google Scholar] [CrossRef]
  27. Renewable Resource Atlas, King Abdullah City for Atomic and Renewable Energy (K.A.CARE). Available online: https://www.energy.gov.sa/en/projects/Pages/atlas.aspx (accessed on 23 August 2022).
  28. Open Energy Data Initiative (OEDI). The U.S. Department of Energy (DOE). Available online: https://data.openei.org (accessed on 20 September 2022).
  29. Bahramirad, S.; Camm, E. Practical modeling of Smart Grid SMS storage management system in a microgrid. In Proceedings of the PES T&D 2012, Orlando, FL, USA, 7–10 May 2012; pp. 1–7. [Google Scholar] [CrossRef]
  30. NREL System Advisor Model (SAM). Available online: https://sam.nrel.gov (accessed on 12 June 2022).
  31. Spertino, F.; Chiodo, E.; Ciocia, A.; Malgaroli, G.; Ratclif, A. Maintenance Activity, Reliability, Availability, and Related Energy Losses in Ten Operating Photovoltaic Systems up to 1.8 MW. IEEE Trans. Ind. Appl. 2020, 57, 83–93. [Google Scholar] [CrossRef]
  32. Schoenung, S. Energy Storage Systems Cost Update: A Study for the DOE Energy Storage Systems Program; Department of Energy: USA, 2011. [Google Scholar] [CrossRef]
Figure 1. Flowchart of the proposed model.
Figure 1. Flowchart of the proposed model.
Sustainability 15 01517 g001
Figure 2. SCG availability model.
Figure 2. SCG availability model.
Sustainability 15 01517 g002
Figure 3. SCG output probability tree.
Figure 3. SCG output probability tree.
Sustainability 15 01517 g003
Figure 4. SF availability PMF (ρSFA(cSF)).
Figure 4. SF availability PMF (ρSFA(cSF)).
Sustainability 15 01517 g004
Figure 5. Solar radiation seasonal samples.
Figure 5. Solar radiation seasonal samples.
Sustainability 15 01517 g005
Figure 6. Expected seasonal power output of SF.
Figure 6. Expected seasonal power output of SF.
Sustainability 15 01517 g006
Figure 7. Seasonal expected composite demand.
Figure 7. Seasonal expected composite demand.
Sustainability 15 01517 g007
Figure 8. Seasonal variations of ρCD (cd).
Figure 8. Seasonal variations of ρCD (cd).
Sustainability 15 01517 g008
Figure 9. Seasonal variation of the ESS SOC.
Figure 9. Seasonal variation of the ESS SOC.
Sustainability 15 01517 g009
Figure 10. Seasonal reserve provision by the ESS and CG.
Figure 10. Seasonal reserve provision by the ESS and CG.
Sustainability 15 01517 g010
Table 1. Variables and parameters used in the model.
Table 1. Variables and parameters used in the model.
Parameter/Variable Description (unit)
Parameters and variables for CG units:
CTgCold start time (h)
DTg/UTgMinimum down/up times of a CG (h)
SUg/SDgStartup/shutdown limit of a CG (MWh−1)
SUcost/SDcostgth CG linearized startup/shutdown costs
sg,tBinary variable of the gth CG startup status (1: turned on, 0: shut down)
rg,tgth CG reserve provision (MW)
RUg/RDgRamp-up/down rates of a CG (MW)
P g m a x / P g m i n Maximum/minimum generation limits of a CG (MW)
pg,t is gth CG produced power (MW)
p ¯ g , t gth CG provided power and reserve (MW)
xg,tBinary variable of the gth CG status (0: off, 1: on)
yg,tgth CG linearized production cost
zg,tBinary variable of the gth CG shutdown status (1: turned off, 0: otherwise)
Parameters and variables for ESS:
αtBinary variable to prevent ESS’s simultaneous charging and discharging
dESS energy capacity capital cost ($/MWh)
eESS power capacity capital cost ($/MW)
EESSEnergy capacity of ESS (MWh)
EtESS energy level or state of charge at time t
ηch/ηdisCharging/discharging efficiency
rESS_DN,tDown reserve provided by the ESS (MW)
rESS_UP,tUp reserve provided by the ESS (MW)
RReserve requirement (MW)
pch,tPower charging of ESS (MW)
pdis,tPower discharging from ESS (MW)
PESSMaximum discharge/charge rate of ESS (MW)
p ¯ E S S _ D N , t ESS charging power and down reserve (MW)
p ¯ E S S _ U P , t ESS discharging power and up reserve (MW)
Table 2. CG units’ characteristics [29].
Table 2. CG units’ characteristics [29].
UnitCost Coeff. ($/MWh)Min. Capacity
(MW)
Max. Capacity
(MW)
Startup Cost ($)
127.71540
239.11540
UnitShutdown Cost ($)Min. Up Time
(h)
Min. Down Time (MW)Ramp Up/Down Rate
(MW/h)
10332.5
20332.5
Table 3. SCGs’ specifications [30,31].
Table 3. SCGs’ specifications [30,31].
SpecificationDescription
pSCG.rated0.05 MW
NSCG50
Gstd1000 W/m2
RC150 W/m2
qSCG0.1667
MTTFSCG1500 h.
MTTRSCG150 h.
Table 4. CG units’ characteristics [32].
Table 4. CG units’ characteristics [32].
SpecificationDescription
ESS technologyLi-Ion
ηchdis85%
Energy capital cost600k $/MWh
Power capital cost400k $/MW
Lifetime20 years
Discount rate5%
e51,814 $/MW
d77,720 $/MWh
Table 5. Demand data [28].
Table 5. Demand data [28].
Demand TypeAverage (KWh/Day)Average (KW)Peak (KW)Demand Factor#Units
Secondary school10,086420.251212.20.351
Primary school2656110.67371.780.31
Midrise apartment building749.8131.2473.640.4220
Medium office2022.584.27254.420.331
Outpatient clinic3956.2164.84360.770.461
Fast food restaurant560.4823.3541.780.565
Large office17,831742.9915310.491
Independent retailer923.938.5104.830.375
Table 6. ESS sizing and CG operational cost results.
Table 6. ESS sizing and CG operational cost results.
CaseEESS (MWh)PESS (MW)Total Cost ($)ESS Cost ($)CG Cost ($)
CG and PV (no ESS)NA NA718,638 ± 1042NA718,638 ± 1042
CG and PV+ESS3.0 ± 0.301.70 ± 0.40815,465 ± 22,259320,514 ± 34,950494,951 ± 20,857
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

Alamri, A.; AlKassem, A.; Draou, A. Composite Demand-Based Energy Storage Sizing for an Isolated Microgrid System. Sustainability 2023, 15, 1517. https://doi.org/10.3390/su15021517

AMA Style

Alamri A, AlKassem A, Draou A. Composite Demand-Based Energy Storage Sizing for an Isolated Microgrid System. Sustainability. 2023; 15(2):1517. https://doi.org/10.3390/su15021517

Chicago/Turabian Style

Alamri, Abdullah, Abdulrahman AlKassem, and Azeddine Draou. 2023. "Composite Demand-Based Energy Storage Sizing for an Isolated Microgrid System" Sustainability 15, no. 2: 1517. https://doi.org/10.3390/su15021517

APA Style

Alamri, A., AlKassem, A., & Draou, A. (2023). Composite Demand-Based Energy Storage Sizing for an Isolated Microgrid System. Sustainability, 15(2), 1517. https://doi.org/10.3390/su15021517

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