Next Article in Journal
Functional Properties and Microbiological Stability of Fatty Acid Methyl Esters (FAME) under Different Storage Conditions
Next Article in Special Issue
High-Performance 3-Phase 5-Level E-Type Multilevel–Multicell Converters for Microgrids
Previous Article in Journal
Glucose-Oxygen Biofuel Cell with Biotic and Abiotic Catalysts: Experimental Research and Mathematical Modeling
Previous Article in Special Issue
Impact of Submodule Faults on the Performance of Modular Multilevel Converters
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Reliability Assessment of Wind-Solar PV Integrated Distribution System Using Electrical Loss Minimization Technique

by
Sachin Kumar
1,
Kumari Sarita
1,
Akanksha Singh S Vardhan
2,
Rajvikram Madurai Elavarasan
3,*,
R. K. Saket
1 and
Narottam Das
4,5,*
1
Department of Electrical Engineering, Indian Institute of Technology (BHU), Varanasi 221005, India
2
Department of Electrical Engineering, Shri G.S. Institute of Technology and Science, Indore 452003, Madhya Pradesh, India
3
Electrical and Automotive Parts Manufacturing Unit, AA Industries, Chennai 600123, Tamilnadu, India
4
School of Engineering and Technology, Central Queensland University, Melbourne, VIC 3000, Australia
5
Centre for Intelligent Systems, School of Engineering and Technology, Central Queensland University, Brisbane, QLD 4000, Australia
*
Authors to whom correspondence should be addressed.
Energies 2020, 13(21), 5631; https://doi.org/10.3390/en13215631
Submission received: 31 August 2020 / Revised: 16 October 2020 / Accepted: 22 October 2020 / Published: 28 October 2020

Abstract

:
This article presents the Reliability Assessment (RA) of renewable energy interfaced Electrical Distribution System (EDS) considering the electrical loss minimization (ELM). ELM aims at minimizing the detrimental effect of real power and reactive power losses in the EDS. Some techniques, including integration of Renewable Energy Source (RES), network reconfiguration, and expansion planning, have been suggested in the literature for achieving ELM. The optimal RES integration (also referred to as Distributed Generation (DG)) is one of the globally accepted techniques to achieve minimization of electrical losses. Therefore, first, the locations to accommodate these DGs are obtained by implementing two indexes, namely Index-1 for single DG and Index-2 for multiple DGs. Second, a Constriction Factor-based Particle Swarm Optimization (CF-PSO) technique is applied to obtain an optimal sizing(s) of the DGs for achieving the ELM. Third, the RA of the EDS is performed using the optimal location(s) and sizing(s) of the RESs (i.e., Solar photovoltaic (SPV) and Wind Turbine Generator (WTG)). Moreover, a Battery Storage System (BSS) is also incorporated optimally with the RESs to further achieve the ELM and to improve the system’s reliability. The result analysis is performed by considering the power output rating of WTG-GE’s V162-5.6MW (IECS), SPV-Sunpower’s SPR-P5-545-UPP, and BSS-Freqcon’s BESS-3000 (i.e., Battery Energy Storage System 3000), which are provided by the corresponding manufacturers. According to the outcomes of the study, the results are found to be coherent with those obtained using other techniques that are available in the literature. These results are considered for the RA of the EDS. RA is further analyzed considering the uncertainties in reliability data of WTG and SPV, including the failure rate and the repair time. The RA of optimally placed DGs is performed by considering the electrical loss minimization. It is inferred that the reliability of the EDS improves by contemplating suitable reliability data of optimally integrated DGs.

Graphical Abstract

1. Introduction

1.1. Literature Survey

The incorporation of distributed power sources in the EDS can support the increasing load demands. To satisfy the increasing load demands, an extremely thorough comparison and assessment has been performed on exploring the potentiality of the distributed energy sources, including the wind and solar energies [1,2]. The EDS is the combination of electrical loads and Distributed Generation (DG) (especially, WTG and SPV) integrated with electrical storage system (namely BSS) [3]. In this regard, several comprehensive studies are available on sustainable energy production and management. The state policies, renewable energy plants, and the development of renewable energy sources including solar energy, wind energy, small hydroelectric energy, biomass energy, tidal energy, and geothermal energy have been analyzed for different states of India [4,5,6]. As the energy crisis deepens day by day, the EDS accompanied by RESs is a better solution and also it acts as a complement to the central power grid system [7]. Therefore, a study on renewable energy potential has been performed for five countries, namely China, Iceland, India, Sweden, and the United States of America. The strengths, weaknesses, opportunities, and threats related to green energy generation and use have been reviewed for these countries using a SWOT analysis [8,9]. Sometimes an arrangement including the battery storage is needed due to high penetration of renewable resources into an electrical network to fulfill the load demands [10,11]. Thus, integration of BSS in the distribution system prevails as a better solution for obtaining a steady power output, especially from WTG and SPV owing to the uncertainties involved in the energy harvesting sources such as wind speed [1], solar irradiation, and ambient temperature [12]. Therefore, it is a requisite to determine the siting(s) and sizing(s) of WTG, SPV, and BSS. The system reliability, stability, and power quality are thus, improved substantially. Therefore, to obtain improvements in system’s attributes, the multi-objective optimization methods for problems including operational cost, siting and sizing of distributed energy resources (DERs), Carbon Dioxide (CO 2 ) emission, total power loss, voltage deviation, demand-side management, charging-discharging of BSS, total harmonic distortion, and system reliability are also established [13,14,15,16,17].
Improvement in network reliability is observed when congestion management algorithms are implemented to identify the transmission line congestion [18]. The DGs and BSSs are scheduled optimally to alleviate this transmission line congestion. A two-step optimization approach is used for solving the congestion problem [19]. In this optimization method, the optimum location and the size of the SPV array are observed, and then the BSS size with location is determined for accomplishing further reduction in electrical losses and voltage deviation. The authors in [20] have proposed an index to obtain the optimal siting of DGs in EDS. This index is implemented to resolve the multiple problems, including total ELM, ENS, and voltage deviation. Another objective, namely loss of yearly energy minimization, is observed through the integration of DERs-based DGs and network reconfiguration [21]. Furthermore, a novel two-stage stochastic programming is proposed, and the uncertainty considerations together with the load variation are also studied especially for wind energy and solar power generation [22]. In this method, the total cost is reduced by incorporating BSS into the EDS and by considering the demand response programs in planning. Simultaneously, the enhancement of the power system reliability is achieved due to the obtained optimal BSS size and location. The reliability improvement and reduction in network losses are also observed using compound co-optimization strategic plan [23]. The reliability indices such as expected ENS, EIR, LOLE, and LOLP, are defined in the co-optimization strategy. Furthermore, a moth-flame optimization [24], Olympic games ranking process [25], firefly algorithm [26], lightning search algorithm [27], crow search optimization [28], and an improved variant PSO [29] techniques have been implemented and discussed in the literature so as to obtain the optimal site, optimal size, optimal parameters of DGs, and ELM.
The RA provides a better evaluation of any power system’s performance [30]. For assessing the power system’s reliability, some indices have been mentioned in the available literature. These indices are categorized into load-based and system-based indices. The indices values decrease if the ageing of the sub-components is considered [31]. Reliability is considered to be a primary requirement in the designing phase of EDS. Thus, the optimal siting and optimal sizing of RESs for further RA have also been considered. To fulfill the RA in EDS [15], has introduced a restoration strategy for ENS calculation. The optimization of reliability indices has been considered in [20,23], and the improvement in the system’s reliability is observed.

1.2. Motivation

The complexity of the electrical power system relies in its ability to tackle the unexpected load variations. These load variations are managed by the electrical distribution system (EDS) for continuous electrical power supply to the loads. This system is categorized into radial and loop structured distribution systems. The use of radial DS (RDS) is preferred in this study as it is simple, cheap, and mostly applicable to sparsely distributed loads. Furthermore, to meet the rapid growth in demand as well as to promote sustainability, integration of renewable energy source (RES) in the distribution system is desideratum. The RESs, including Wind Turbine Generator (WTG) and Solar photovoltaic (SPV), provide minimal electrical losses, improves system bus voltages, possesses less operation costs, and more significantly, it emits less CO 2 emissions. Thus, the optimization of these parameters would ultimately improves the reliability of the EDS. Reliability in an electrical power system is defined as the ability to provide adequate, stable, and reliable power for a particular distribution system. Therefore, the study on Reliability Assessment (RA) and performance analysis of RDS considering the optimal RES siting and sizing is crucial. It is also observed that the integration of renewable sources in EDS has considerable impacts on systems’ operation and planning strategy [32]. Some issues co-related to RES integration are explained as follows.
  • Electrical loss minimization using system reconfiguration [13].
  • Electrical loss minimization [33,34],
  • Reduction in investment during system capacity enhancement [35],
  • Improvement in bus voltage [36],
  • Mitigation of greenhouse gases [37],
  • Improvement in voltage stability [38],
  • Improvement in system’s reliability by considering the reliability indices [39,40,41],
  • Enhancement in system security [42],
  • Facilitate system restoration [43],
  • Reduction in harmonic distortion [44],
  • Optimal load management strategy [45].
The mentioned issues have been analyzed individually that is without considering the effect on another issue(s). The concept of analyzing the effect of the mentioned issues motivates the authors for carrying out this research work. In the view of strong motivation, this work is concentrated on the reliability assessment of IEEE 33 bus distribution system by considering the effect of electrical loss minimization through optimal siting and sizing of distributed generation system.

1.3. Contribution

The available literature focus on optimal siting and sizing of DGs obtaining a better power system performance. Also, the DG integration of SPV, WTG, and BSS in the distribution system has not been discussed as of authors knowledge. The system RA is discussed individually without analyzing the impact of optimal DG integration. Thus, the two indexes have been implemented for obtaining the optimal location(s) of DG(s) in the selected EDS. Then, a Constriction Factor (CF)-based PSO technique for optimal sizing(s) of DG(s) is presented. Simultaneously, all the parameters considered in this study (as given in Table 1) are not discussed in the existing literature and thus, this work provides completeness in contribution. According to the outcomes of the study, the results are found coherent with those obtained using other techniques that are available in the literature. However, the results can be claimed to be even better since the proposed approach considers all the parameters given in Table 1, including the ELM and the reliability aspects. The comparative analysis for a reduction in electrical power loss and voltage deviation is performed by integrating the WTG, SPV, and BSS in IEEE 33 bus EDS. The results obtained are compared with the results in the literature, as described in Table 2. After obtaining the reduced power loss and improved bus voltages, the reliability assessment is accomplished. The performance analysis is done for the selected distribution system and then the system’s RA is carried out with and without the integration of DGs.
Though the research on DG siting and sizing is abundant in the literature, the RA in the EDS still prevails as an emerging area of research. Thus, the work mainly focuses on the reliability assessment of an IEEE 33 bus system is carried out in the view of loss minimization technique of WTG, SPV, and BSS integration. For the accomplishment of the mentioned tasks, the research work is performed according to the framed flowchart as shown in Figure 1. Furthermore, the following contributions are highlighted in this work.
  • The optimal locations for SPV, WTG, and BSS is obtained by considering Index-1 and Index-2.
  • Optimal sizes of SPV, WTG, and BSS are derived by employing CF-PSO technique.
  • Distribution system performance, including minimization of I 2 R loss and minimization of deviation in bus voltages, is analyzed with and without DGs.
  • Comparison of CF-PSO technique with other nature-inspired optimization techniques for achieving a sound reliability assessment.
  • A brief study is done on the inclusion of uncertainties in WTG and SPV reliability data such as failure rate ( λ p ) and time to repair (RT).
  • Reliability assessment is done by evaluating the five indices namely EENS, AENS, SAIDI, SAIFI, and ASAI for Case 1, Case 2, and Case 3; where Case 1 is for integrating WTG only, Case 2 is for WTG+SPV, and Case 3 is for WTG+SPV+BSS (adding BSS optimally).
The work is further continued to accomplish the mentioned contributions as follows. The problem statement, objective function, and methodologies implemented, including Index-1, Index-2, and CF-based PSO techniques are delineated in Section 2. The basics of EDS reliability assessment is formulated and explained in Section 3. In Section 4, an overview of SPV, WTG, and BSS is explained using mathematical expressions. Furthermore, the results of several cases and scenarios are elucidated with supporting graphical depictions and tables in Section 5. Finally, the conclusions with future scopes are drawn in Section 6. The Appendix A of some important data and figure is provided at the end of the work.

1.4. Parameters Considered for the Study

The parameters considered for the study are described as follows.
  • DG siting and DG sizing: The determination of bus voltages and the flow of powers is done by an Optimal Power Flow (OPF) method. The optimal siting of DGs is required for ELM during this power flow results. The performance of the power network is affected by an inappropriate location of DG. The IEEE 1547 standards for integration and operation of DG into EDSs are presented in [49].
  • Power loss: The occurrence of Active Power Loss (APL) is greater than Reactive Power Loss (RPL) in EDS. Hence, distribution companies should reduce these losses and this can be accomplished by means of reconfiguration of feeder, capacitor allocation, high voltage distribution system, grading of the conductor, DG placement, and many other methods.
  • Bus voltage: It is expected to maintain bus voltages nearly 1 pu with an angle of 0 0 . The power loss occurring in EDS during OPF creates a voltage drop at each bus of the system. Therefore, the DG integration technique is implemented for voltage profile (VP) improvement.
  • DG type: The three DGs have been considered, which are categorized as WTG, SPV, and BSS. The classification of the several DG technologies is based on the generation of active power ‘P’ and reactive power ‘Q’, as illustrated in Figure 2.
  • Reliability: The reliability indices considered for the distribution system reliability are as follows.
    -
    Expected Energy Not Supplied (EENS); MWh per year
    -
    Average Energy Not Supplied (AENS); MWh per customer per year
    -
    System Average Interruption Duration Index (SAIDI); hour per customer per year
    -
    System Average Interruption Frequency Index (SAIFI); failure per customer per year
    -
    Average System Availability Index (ASAI); pu
The calculations of these five indices are performed using Equations (A1)–(A7) of the Appendix A.2. The reliability indices considered in this work can be used in obtaining other indices, including IEAR, CAIDI and ASUI, as illustrated in Equations (42), (45) and (46c).

2. Problem Formulation, Objective Function (OF), and Methodology

The bus voltage and system reliability are the most affecting factors for the EDS when losses are considered. It becomes necessary to improve these two factors by implementing distributed generation (DG) into the EDS. DG siting is one of the favored techniques used in the EDS for the improvement of the system’s reliability and bus Voltage Profile (VP), including ELM. The DG location, DG size, DG power factor (pf), DG penetration, and DG type are required for the effective implementation of DGs in the EDS. Simultaneously, it is required to study the mathematical expressions, and modeling of related parameters and DGs integrated into the system. An overview of parameters considered for the ELM and mathematical modeling are elaborated in the upcoming subsections.

2.1. Optimal Location

Obtaining the optimal location of DGs is crucial part of EDS. To identify the optimal locations, two indexes are used. Index-1 is implemented only for placing the single DG and Index-2 is incorporated for placing more than one DGs in the EDS. Power loss is minimized by using Index-1 for placing 1DG (viz. Case 1). However, Index-2 provides minimum power loss for multiple DGs (viz. Case 2 and Case 3). These two indicators are represented by Equations (17) and (18), respectively [50,51]. It can be observed from the equation of Index-1 that the large index value depicts the weakest node of the system because the complex power injected at bus i is large. It implies that the single DG can be placed at this particular bus. On the other hand, Index-2 shows the voltage stability, which concludes that the reduced values of this index give the weakest bus of the EDS. Table 3 shows the values of both the indexes with corresponding five buses to arrange DG optimally in a given EDS. Therefore, DGs can be placed hierarchically at these buses.
In Figure 3, V a and V b are the magnitudes of the voltages at buses a and b, respectively. δ a and δ b are the phase angles of the voltages at buses a and b, respectively. Z l and Y l are the impedance and admittance of l-line, respectively. R l and X l are the resistance and reactance of a l-line. I l is the current in the l-line. The electrical power loss in the l-line is given by (1).
S L O S S l = ( V a V b ) × I l *
I l = ( V a V b ) × Y l
Then the bus voltage matrix is formed by using Equation (3) where [ Z b u s ] is the network impedance matrix, [ I b u s ] is the bus injection matrix, and n b u s is the number of total buses in EDS.
[ V b u s ] n b u s × 1 = [ Z b u s ] n b u s × n b u s [ I b u s ] n b u s × 1
By expanding Equation (3), the node voltages can be obtained by using Equations (4) and (5).
V a = i = 1 n b u s Z a i × I i
V b = i = 1 n b u s Z b i × I i
where i is 1,2,…, n b u s . Z a i , Z b i and I i are the element of impedance matrix that signify the a t h row and i t h column, b t h row and i t h column, and current injection at bus-i, respectively.
Current Injection , I i = ( P i n , i + j Q i n , i ) * V i *
where P i n , i and j Q i n , i are active power and reactive power injected at bus-i, respectively. V i * is the voltage at bus-i. Now, put (2) and (4) to (6) in (1) then the electrical loss of the line-l is derived as (7).
S L O S S , l = [ V a V b ] n b u s ( Z a i Z b i ) Y l V i * * [ P i n , i + j Q i n , i ]
For an electrical system with n l number of branch/lines, the line loss is given by (8).
[ B L O S S l ] = i = 1 n b u s ( V a V b ) ( Z a i Z b i ) * Y l * V i S i n , i
where B L O S S l is line loss, S i n , i is apparent power injected at bus-i, l is 1 to n l .
[ B L O S S l ] = i = 1 n b u s [ L F l i ] [ S i n , i ]
L F l i = i = 1 n b u s ( V a V b ) ( Z a i Z b i ) * Y l * V i
L F l i = non - zero , if l - line is in the path of bus - i 0 , Else
where load factor (LF) is given by ( V a V b ) ( Z a i Z b i ) * Y l * V i . L F l i is load factor of the l t h line due to the i t h bus injection (it is non-zero if l t h line is in the path of i t h bus else zero) as described in (11). For example, a 6-bus distribution is taken for explanation as shown in Figure 4 and a general branch loss formula is derived as (12).
B L O S S 1 B L O S S 2 B L O S S 3 B L O S S 4 B L O S S 5 5 × 1 = L F 11 L F 12 L F 13 L F 14 L F 15 L F 16 L F 21 L F 22 L F 23 L F 24 L F 25 L F 26 L F 31 L F 32 L F 33 L F 34 L F 35 L F 36 L F 41 L F 42 L F 43 L F 44 L F 45 L F 46 L F 51 L F 52 L F 53 L F 54 L F 55 L F 56 5 × 6 × S i n 1 S i n 2 S i n 3 S i n 4 S i n 5 S i n 6 6 × 1
[ B L O S S M ] 5 × 1 = [ L F M ] 5 × 6 × [ S i n M ] 6 × 1
where [ B L O S S M ] is branch/line loss matrix, [ L F M ] is the load factor matrix, [ S i n M ] is complex power injection matrix. Equation (13) is reduced according to the Figure 4 where power injection at all the buses except the source bus are available.
B L O S S 1 B L O S S 2 B L O S S 3 B L O S S 4 B L O S S 5 5 × 1 = 0 L F 12 L F 13 L F 14 L F 15 L F 16 0 0 L F 23 L F 24 L F 25 L F 26 0 0 0 L F 34 L F 35 0 0 0 0 0 L F 45 0 0 0 0 0 0 L F 56 5 × 6 × S i n 1 S i n 2 S i n 3 S i n 4 S i n 5 S i n 6 6 × 1
A i = l = 1 n l L F l i
The calculation of effective power injections is done as given in (16).
c S e f f , 6 = S i n , 6 S e f f , 5 = S i n , 5 S e f f , 4 = S i n , 4 + S e f f , 5 S e f f , 3 = S i n , 3 + S e f f , 4 + S e f f , 6 S e f f , 2 = S i n , 2 + S e f f , 3 S e f f , 1 = S i n , 1 + S e f f , 2
Equation (17) shows the final equation for the calculation of Index-1. The index is implemented for attaining the optimal position of one renewable energy source as DG. In Equation (15), | A i | is fully dependent on LF values of all the branches (or lines) connected between bus and the source bus (main station). The closeness of the i t h bus from the source bus can be observed in the L F l i as guided in Equation (10). If the i t h bus is not near, the number of lines between i t h bus and source bus is being plenty and the corresponding Z a i , Z b i , and Y l parameters will account in the electrical loss component. Furthermore, if the node voltage is high, the value of L F l i will be small and vice-versa as observed in the derived equation. The equation of Index-1 is also accounted for the effective complex power supplied by the i t h bus. The Index-1 will be high only when both the terms are high in (17). Thus, the value of Index-1 represents its main contribution in the total electrical loss and hence, in finding the optimal siting of one RES.
( I n d e x 1 ) i = | A i | × | S i n e f f i |
S i n e f f i is the effective injection of complex power, which is the sum of injected powers from other buses connected to i t h bus as shown in Figure 4.
( I n d e x 2 ) i + 1 = | V i | 4 4 ( P i + 1 X j Q i + 1 R j ) 2 4 ( P i + 1 R j Q i + 1 X j ) | V i | 2
where j is branch number, V i is sending bus voltage, P i + 1 and Q i + 1 are the Active Power (AP) and Reactive Power (RP) at the receiving end bus, respectively. R j and X j are the resistance and reactance between sending and receiving end bus, respectively.

2.2. Power Balance

The AP and RP balance expressions are shown in Equations (19) and (20).
P n e t i = P d g i P d e m i V i j = 1 N b u s V j Y i , j c o s ( δ i δ j θ i + θ j )
Q n e t i = Q d g i Q d e m i V i j = 1 N b u s V j Y i , j s i n ( δ i δ j θ i + θ j )
where P n e t i = 0 and Q n e t i = 0 are the net AP and RP at i-bus, respectively. P d g i and Q d g i represent DG AP and RP at i-bus, respectively. Active and reactive load demands are mentioned by P d e m i and Q d e m i , respectively. V j is the bus voltage at j-bus, Y i , j is the branch admittance between i, j-buses, δ i and δ j represent the phase angles of i-bus and j-bus voltages, respectively. ( θ i θ j ) are the impedance angle of branch connected between i and j-buses.

2.3. Objective Function (OF)

In this research the OF is considered to be APL minimization in the EDS. The reliability indices are then evaluated by fixing the optimal location and size of the DGs. The OF of the problem is given in Equation (21).

2.3.1. Active Power Loss (APL)

The minimization of APL occurred in EDS is the OF considered. The primary aim of the OPF technique is to minimize the system APL as given as Equation (21).
min A P l o s s = i = 1 N b u s j = 1 N b u s C 1 i j ( P r e a l i P r e a l j + Q r e a l i Q r e a l j ) + C 2 i j ( Q r e a l i P r e a l j P r e a l i Q r e a l j )
where P r e a l i , P r e a l j , Q r e a l i , Q r e a l j are the AP and RP at i and j-buses, respectively. N b u s = number of buses or nodes, C 1 i j and C 2 i j are defined as follows.
C 1 i j = R i j V i V j c o s ( δ i δ j )
C 2 i j = R i j V i V j s i n ( δ i δ j )
where V i , δ i and V j , δ j are the voltages and corresponding angles at i t h and j t h buses, respectively, R i j = resistance of a branch between i and j-buses.

2.3.2. Reactive Power Loss (RPL)

The availability of RP ensures the AP transmission from source to load. Voltage stability margin or bus voltages are also dependent on this RP support. The RPL is obtained at different pf of DG using Equation (23).
R P l o s s = i = 1 N b u s Q g e n i i = 1 N b u s Q d e m i
where Q g e n i and Q d e m i are the RP generation and demand at the i t h bus (including the slack bus), respectively. Q d e m i = RP demand at the i t h bus.

2.3.3. Reliability Indices

The indices are assessed for divergent DG reliability data by fixing the site and size of DGs. Furthermore, the reliability improvement of distribution network has been accomplished by integrating one DG and multiple DGs in the EDS. Several reliability indices exist to observe the system’s reliability such as EENS, AENS, SAIDI, SAIFI, and ASAI which are also used in this study to analyze the reliability improvement. The indices of the network reliability are dependent function of failure rate ( λ p ) and repair time (RT) as given in Equation (24) [20,23,52].
Reliability Indices = f ( λ p , R T )

2.4. Constraints

The OF minimization is a primary task to obtain the optimal results. OF minimization is subjected to design the constraints so that the requirements of the EDS must be satisfied with DG operation. Thus, the constraints are discussed in the succeeding subsection [53,54].

2.4.1. Equality Constraints

These constraints follow the Kirchhoff’s current rule as the algebraic sum of powers in and powers out should be equal in an EDS [54,55]. Two of these constraints are described as follows.
i = 1 N b u s A P g e n i = i = 1 N b u s A P d e m i + A P l o s s
where A P g e n i = AP generated by the generation units at i t h bus, A P d e m i = AP demand at i t h bus.
i = 1 N b u s R P g e n i = i = 1 N b u s R P d e m i + R P l o s s
where R P g e n i = RP generated by the generation units at i t h bus, R P d e m i = RP demand at i t h bus.

2.4.2. Inequality Constraints

These constraints are associated with the limits applied to the system parameters for the operation of EDS. Some of these constraints are described as follows.
A. Power flow
To maintain the line capacity within limits, these constraints ensure the apparent power to be within limits at the ends of a line [53,54].
A P a i j A P a i j m a x
where A P a i j m a x = the highest permissible apparent powers ( A P a ) for lines i to j, A P a i j = the actual A P a transmitted from i to j.
B. DG capacity
These limits ensure the non-reversal of power flow. The power from the substation is provided to the EDS must be greater than the DG power. Also, the DG has the minimum and maximum power generation boundaries [56].
i = 1 n D G A P D G i 0.85 × i = 1 n b u s A P d e m i + A P l o s s
i = 1 n D G R P D G i 0.85 × i = 1 n b u s R P d e m i + R P l o s s
A P D G p m i n A P D G p A P D G p m a x
R P D G p m i n R P D G p R P D G p m a x
where p = 1, 2,……, n D G , A P D G p m i n (set to zero) and A P D G p m a x (from Equation (28)) are the lower and upper AP outputs of DG unit p, respectively. R P D G p m i n (set to zero) and R P D G p m a x (from Equation (29)) are the lower and upper RP outputs of DG unit p, respectively. n D G = number of DGs present in the distribution network.
C. Bus voltage
The voltages at buses present in the EDS must be limited within minimum and maximum limits [57,58].
| V i m i n i m u m | | V i | | V i m a x i m u m |
where | V i m i n i m u m | and | V i m a x i m u m | = the lower and upper boundaries of the bus voltage | V i | which are set to 95% and 105%, respectively.
D. Bus voltage
It refers as thermal capacity of the EDS lines. The current in the distribution lines must be within limits and should exceed the maximum current as given in Equation (33) [20].
I i I i m a x

2.5. Constriction Factor-Based PSO (CF-PSO) Technique

PSO is a novel progression computational technique which is in the frame since 1995. The use of this method is seen in RP dispatch [59], generation scheduling [60], renewable source integrated power system [61], and cost analysis [62]. In basic PSO method, the candidate solution is improved iteratively under any given constraint. The PSO algorithm is shown in Figure 5. Due to the reduction in computational time and requirement of less memory, PSO has overtaken many algorithms including the Genetic algorithm (GA) as PSO is mutation free. It searches the optimized value globally with the help of several particles present in a swarm based on specific constraints. As all particles have its local and global best values because of its own and global positions. This method updates the particle position and velocity as described in Equations (34) and (35).
V n p + 1 = W × V n p + C 1 × R 1 × ( P e r s o n a l B E S T i X n p ) + C 2 × R 2 × ( G l o b a l B E S T X n p )
X n p + 1 = X n p + C f × V n p + 1
where V n p + 1 = n t h particle velocity at ( p + 1 ) t h iteration, W’ = particle inertial weight, V n p = n t h particle velocity at p t h iteration, C 1 , C 2 = constants (0, 2.5), R 1 , R 2 = numbers generated randomly (0, 1), P e r s o n a l B E S T i = the n t h particle’s best position considering its own property, G l o b a l B E S T = the n t h particle’s best position considering the whole population, X n p + 1 , X n p = n t h particle position at ( p + 1 ) t h and p t h iterations, respectively. C f = constriction factor (CF) assures efficient convergence [63,64].
Due to faster convergence to the global point, the basic PSO faces the difficulty of premature convergence. The particles have started oscillating around the optimal point without providing any type of restriction to the highest velocity of the particles available in swarm. Therefore, the optimal global solution is rare to obtain. The use of properly defined CF is briefly described for advance convergence of the PSO [65]. This can also be applied for the DG siting and DG size in the EDS. It reduces the computation time and requires little memory. Although this technique suffers from partial optimization, by altering its parameter during problem solving will produce an improved result [66,67,68]. To obtain the improved result, a constriction factor (CF) is used and thus, the method is known as CF-based PSO technique. The parameters set for the CF-PSO are as follows. The values of initial weight, final weight, C 1 , C 2 , R 1 , R 2 , and CF are considered to be 9 × 10 1 , 4 × 10 1 , 201 × 10 2 , 201 × 10 2 , 0 to 1, 0 to 1, and 729 × 10 3 , respectively. A flowchart is provided in Figure 5 to obtain the DG location, DG sizing and system reliability of the DGs in 33 bus EDS.

3. Reliability Assessment of Distribution System

The RA of an EDS is as important as contrasted to other components and parts of the EDS. The IEEE guide for EDS reliability is given by standard number 1366-2012 [69]. According to given standard, reliability of an EDS can be analyzed using some reliability indices. The reliability indices considered for EDS reliability assessment include EENS, AENS, SAIDI, SAIFI, and ASAI. These indices are mainly classified in two categories which are elaborated in Equations (39)–(46c).

3.1. Reliability Parameters at Load Point ‘p’

The reliability indices are the function of reliability parameters mentioned in Equations (36)–(38). The reliability parameters have been calculated at load point ’p’ as follows.
Failure rate ( average ) ; λ p = k n n u m k × F k failure per year
Outage duration ( annual ) ; U p = k = n F k D p k hour per year
Outage duration ( average ) ; D p = U p λ p hour
where F k = failure rate (average) of the k t h element, n = number of elements in the EDS, n u m k = number of k t h elements in the EDS, D p k = duration of failure at p t h load point due to k t h failed element. The calculation equations of λ p and U p are given in Appendix A.2.

3.2. System-Based Indices

These indices are further categorized in load-oriented indices and customer oriented indices as given in Equations (39)–(40) and Equations (41)–(46c), respectively.

3.2.1. Load-Oriented Indices

Load-Oriented Indices have been calculated at load point ’p’, as mentioned in Equations (39)–(40).
E E N S p = P p U p megawatt hour per year
A E N S p = p = 1 n p E E N S p p = 1 n p N p megawatt hour per customer per year
where P p = demand/load (average) of the p t h load point, E E N S p = expected ENS at P t h load or customer point. where n p = total load points, N p = number of customers at P t h load point.

3.2.2. Customer Oriented Indices

These indices have allowed to enhance the EDS’s reliability related to the improvement of customer or load services. The two of the indices namely ECOST and IEAR are related to the cost reliability and thus, termed as reliability worth of the system.
E C O S T p ( = L O E E p ) = P p k = n f ( D p k ) F k k $ per year
I E A R p = E C O S T p E E N S p $ per kilowatt hour
where E C O S T p = expected interrupted cost at P t h load point, I E A R p = interrupted energy assessment rate at P t h load point, L O E E p = loss of expected energy, f ( D p k ) = system composite customer damage function ($ per kilowatt) as provided in Table A3 of Appendix A.
S A I F I = p = 1 n p λ p N p p = 1 n p N p failure per customer per year
S A I D I = p = 1 n p U p N p p = 1 n p N p hour per customer per year
CAIDI ( = S A I D I S A I F I ) = p = 1 n p U p N p p = 1 n p λ p N p hour per customer per interruption
where CAIDI = Customer Average Interruption Duration Index.
A S A I = 8760 p = 1 n p N p p = 1 n p U p N p 8760 p = 1 n p N p per unit
Also, ASAI can be derived as follows.
A S A I = 1 S A I D I 8760
A S U I = 1 A S A I per unit
where ASUI = Average Service Unavailability Index.

4. Modeling of WTG, SPV, and BSS

The reliability assessment of the IEEE 33 bus EDS is accomplished, considering the optimal siting(s) and sizing(s) of SPV, WTG, and BSS. In this regard, a brief modeling and specifications of these RESs are illustrated in Section 4.1, Section 4.2 and Section 4.3.

4.1. Wind Turbine Generator

The V162-5.6MW(IECS based on IEC IIB), a WTG, manufactured by General Electric Company is considered for its output power rating. The specifications of the WTG considered in this study are provided in Table 4. The mechanical power of WTG ( P m e c h ) is a function of generator rotor speed and wind speed as formulated in Equations (47)–(51) [70].
P m e c h ( v w i n d , ω r o t o r ) = 1 2 × ρ × v w i n d 3 × C p ( λ , θ )
λ = ω r o t o r × G R × R r o t o r v w i n d
C p ( v w i n d , ω r o t o r , θ ) = C 1 ( C 2 1 α C 3 θ C 4 θ x C 5 ) × e x p ( C 6 α )
1 α = 1 ( λ + 0.08 θ ) 0.035 1 + θ 3
where ρ = air density, A s = area swept by the turbine rotor blades, V w i n d = speed of the wind, C p is the non-linear function of the tip speed ratio ( λ ) and pitch angle ( θ ), ω r o t o r = generator rotor speed, GR = gear ratio, R r o t o r = rotor radius at the turbine blades, C 1 to C 6 and x are constants and computed in [70].
P W T G = 0 ; 0 V V c i n or V V c o u t P W T G , r a t e d × V V c i n V r a t e d V c i n ; V c i n V V r a t e d P W T G , r a t e d ; V r a t e d V V c o u t
where P W T G = output WTG power, P W T G , r a t e d = rated output WTG power, V r a t e d = rated wind speed.

4.2. Solar Photovoltaic

The SPR-P5-545-UPP, a Solar PV Module, manufactured by Sunpower Company is considered for its output power rating. The specifications of the SPV module considered in this study are provided in Table 5. The SPV module is developed by implementing several cells. The power output for the SPV module can be derived as described by Equations (52)–(65) [71,72,73].
P S P V ( A C ) ( t ) = P S P V ( o u t ) ( t ) × η i n v e r t e r
P S P V ( o u t ) ( t ) = F F A ( t ) × I s h o r t ( t ) × V o p e n ( t )
I s h o r t ( t ) = S R 1000 [ I s h o r t , S T C + C I ( T S M ( t ) 25 ) ]
V o p e n ( t ) = V o p e n , S T C + C V [ T S M ( t ) 25 ]
F F A ( t ) = F F i ( t ) × [ 1 R s ( t ) ]
F F i ( t ) = V o p e n , 0 ( t ) l n [ V o p e n , 0 ( t ) + 0.72 ] V o p e n , 0 ( t ) + 1
V o p e n , 0 ( t ) = V o p e n ( t ) × Q N × K [ T S M ( t ) + 273.15 ]
R s ( t ) = R s × I s h o r t ( t ) V o p e n ( t )
T S M ( t ) = T a m b ( t ) + S R ( T n o m 20 ) 0.8
where t = time instant, P S P V ( A C ) ( t ) = AC output power, P S P V ( o u t ) ( t ) = maximum output power, η i n v e r t e r = inverter efficiency, F F A ( t ) = fill factor actual, I s h o r t ( t ) = short circuit current under operating conditions, V o p e n ( t ) = open circuit voltage under operating conditions, SR = solar radiation intensity (W/m 2 ), STC = standard test conditions, C I = temperature coefficient for current (A/ C), T S M ( t ) = solar module temperature ( C), C V = temperature coefficient for voltage (V/ C), F F i ( t ) = fill factor ideal, R s ( t ) = normalized series resistance of solar module, V o p e n , 0 ( t ) = normalized open circuit voltage, Q = an electron charge, N(≈1) = ideality factor, K = Boltzmann’s constant, R s = series resistance of SPV module, T a m b ( t ) = ambient temperature of SPV module, T n o m = nominal operating cell temperature.
The R s is being evaluated as follows [74].
R s = R s , S T C = r s , S T C × V o p e n , S T C I s h o r t , S T C
R s , S T C = 1 F F A , S T C F F i , S T C
F F A , S T C = V m p p , S T C × I m p p , S T C V o p e n , S T C × I s h o r t , S T C
F F i , S T C = V o p e n , 0 , S T C l n [ V o p e n , 0 , S T C + 0.72 ] V o p e n , 0 , S T C + 1
V o p e n , 0 , S T C = V o p e n , S T C × Q N × K [ T S M , S T C ( t ) + 273.15 ]
where r s , S T C = normalized series resistance under STC, and all other parameters in Equations (61)–(65) are evaluated at STC.

4.3. Battery Storage System

The BESS 3000, a Lithium-Ion Battery System, manufactured by Freqcon Company is considered for its output power rating. The specifications of the BSS considered in this study are provided in Table 6. The flowchart of BSS dispatch modeling is illustrated in Figure 6. Mathematical modeling is further defined. The BSS dispatch strategy starts functioning by monitoring the peak load hours. If the peak load is greater/less than the capacity of WTG and SPV, the BSS discharges/charges to support the distribution system; otherwise the BSS operates as a neutral device. Also, the BSS SOC will decide the charging or neutral operation during the off-peak hours. The battery model is derived with the help of Equations (66)–(68) assuming that no battery is self-discharging [75].
P B a t t e r y ( D C ) ( t ) = E B a t t e r y ( t ) E B a t t e r y ( t Δ t ) Δ t
P B a t t e r y ( A C ) ( t ) = P B a t t e r y ( D C ) ( t ) η B a t t e r y ; P B a t t e r y ( D C ) ( t ) > 0 P B a t t e r y ( D C ) ( t ) × η B a t t e r y × P F I n v e r t e r ; O t h e r w i s e
S O C ( t ) = S O C ( t Δ t ) + η C h a r g i n g × P B a t t e r y ( D C ) ( t ) C B a t t e r y × V × Δ t
where P B a t t e r y ( D C ) ( t ) = DC charging/discharging power of the battery in Δt interval (W), E B a t t e r y ( t ) = energy of battery (Wh), P B a t t e r y ( A C ) ( t ) = AC power discharged/charged state of battery, η B a t t e r y = efficiency of the battery, P F I n v e r t e r = inverter power factor, S O C ( t ) = state of charge of the battery, η C h a r g i n g = charging efficiency of the battery, C B a t t e r y = battery capacity (Ah), V = nominal voltage of the battery (V).

5. Results and Discussion

The DG location, DG size, and EDS reliability are obtained and analyzed. The 33 bus EDS (Figure A1 of Appendix A.1) is considered in this study. The branch and load data for this EDS are adopted from [76]. It contains 33 buses and 32-branches with a total of 3.715 MW and 2.3 MVAr AP and RP loads, respectively. This EDS operates at 12.66 kV, 100 MVA base values. The APL and RPL for without DG case are obtained as 0.211009 MW and 0.143056 MVAr, respectively. The objective function (OF) minimization is done by implementing the CF-based PSO, as explained in Section 2. The following steps are followed to obtain the results.
Step 1:
Optimal siting(s) and sizing(s) of WTG, SPV, and BSS are evaluated considering electrical loss minimization (ELM). The technical ratings of WTG, SPV, and BSS have been illustrated in Table 4, Table 5 and Table 6, respectively. The BSS is assumed to be fully charged and produces its rated output power.
Step 2:
APL, RPL, and bus voltages are obtained by integrating WTG, WTG+SPV, and WTG+SPV+BSS (referred as Case 1, Case 2, and Case 3, respectively) in the EDS to analyze the results obtained in Step 1.
Step 3:
Reliability indices are estimated for EDS considering two different WTG and SPV reliability data, including λ p and RT (for Scenario 1 to Scenario 6).
Step 4:
Furthermore, the reliability improvement is analyzed by adding BSS (considering 100% reliable) to the EDS in the presence of WTG and SPV. All related reliability data used are mentioned in Table A2 of the Appendix A.

5.1. DG Location and DG Rating

The bus number is obtained to allocate WTG, SPV, and BSS. The two indexes are implemented to obtain the locations as described in Section 2.1. It is observed from the analysis that Index-1 provides the location suitable for a single DG. This index examines the effective apparent power injection to the buses. The Load Factor value of l t h line depends on whether the l t h is in the path of the i t h bus to the source node or not. The multiplication of Load Factor and injected apparent power provides the APL and RPL of l t h line due to the i t h bus A P a injection. Thus, the maximum value of Index-1 at i t h bus indicates the candidate bus to place one DG. The first six values of this index for respective buses are provided in Table 3.
The optimal siting for several DGs are found by implementing the Index-2. This index provides the hierarchy of weak buses in the EDS. It shows the sensitivity of the bus towards voltage collapse. The value of Index-2 must be greater than or equals to zero for ensuring the stable operation of the distribution system. The minimum value of this index depicts more sensitivity to the voltage collapse and thus, referred to as the weak bus. The optimal location(s) and size(s) of WTG, SPV, and BSS for three cases are obtained and reproduced in Table 7. The DGs are accommodated according to the locations obtained from the indexes, as mentioned in Section 2.1. The DG size and minimum APL are then evaluated, implementing CF-PSO as described in Section 2.

5.2. APL, RPL, and Bus Voltages

The accommodation of DG at an optimal location with optimal size reflects in VP improvement and minimization of APL and RPL. Most of the research has concentrated on APL minimization because of the dominance of I 2 R losses in the EDS. In contrast, the RPL minimization for overall voltage improvement of all the 32 buses of the EDS has also been observed. The estimation of APL, RPL, and VP is considered before analyzing the system’s reliability. This is done to analyze the system’s reliability with the optimal DG size, DG location, minimum power loss, and better VP.
Several kinds of research are performed to obtain APL minimization, which is cited in Table 2 for single and multiple DGs, respectively. It is observed from Table 2 that the authors of the mentioned literature have not dealt with the cases considering various DGs combination. Therefore, it is vital to observe that the results obtained considering several types of DGs are compared with the conventional DGs [77]. The output results obtained for the APL are tabulated in Table 8. This table shows that the APL value is better for Case 1, and comparable for Case 2 and Case 3. The slight variations in APL for Case 2 and Case 3 are observed because the authors have considered the pf of WTG only. Furthermore, the bus voltage profile with one DG and multiple DGs is drawn in Figure 7. The voltages at all the buses vary according to the real power loss and reactive power loss in the distribution system. Therefore, the system requires real power support for active power loss minimization, which improves the bus voltages by compensating the I 2 R losses. Furthermore, it is concluded from Table 8 that the active power loss minimization is not reduced significantly for Case 3 as compared to Case 2. Thus, there is a marginal improvement in bus voltage profile for Case 3 as compared to Case 2, which is depicted in Figure 7. Also, the improvement in bus voltages is observed when multiple DGs are placed. This VP is further improved at 0.85 and 0.82 pfs. This is because of increment in reactive power support at system buses. It is the point of interest to know about the two voltage peaks when the system is operated with single DG. The two voltage peaks appear at bus number 7 and 26 because these buses are directly connected to bus 6, at which single DG is placed optimally. Also, the size of the single DG is greater than the sum of the size of two DG and slightly lesser than the sum of the size of three DG, as obtained in Table 8. A comparison between present work and the best available method is made for ELM. Simultaneously, from Table 9 it can be inferred that the minimum bus voltage is improved and RPL is minimized with the implementation of multiple DGs at different pfs as illustrated. The graphical representation of APL and RPL for without DG, one DG, and multiple DGs are represented in Figure 8.

5.3. Reliability Assessment

The indices obtained show the improvement in the system’s reliability. The indices are calculated for two different reliability data of DGs. It is observed that the indices are dependent on two reliability data, namely λ p and RT of the system’s elements. The present work has considered different reliability data for DG only. The best reliability improvement is observed for 0.2 of λ p and 12 h of RT. A detailed description of the DG reliability data effect on indices is given as per the following Scenarios.
  • Scenario 1: 0.2 f/yr and 12 h (as provided in Table A2 of Appendix A.2)
  • Scenario 2: 0.4 f/yr and 12 h
  • Scenario 3: 0.6 f/yr and 12 h
  • Scenario 4: 0.2 f/yr and 24 h
  • Scenario 5: 0.2 f/yr and 48 h
  • Scenario 6: No failure
Scenarios 1, 2, and 3 are considered by fixing the RT and varying λ p . The appropriate case from the first three cases is then considered for variable RT to extract the best case from the top Five Scenarios. The values of these reliability data are fed manually in the optimization technique to get the values of indices for the system’s reliability improvement. Furthermore, the following key assumptions are considered to assess the reliability of the EDS.
  • Circuit breakers, distribution lines, and potential transformers are available throughout with 100% reliability.
  • The λ p and RT of DG, Buses, feeders, and substations are given in Table A2.
  • RT for each distribution branch = 10 h.

5.3.1. Effect on Load-Oriented Indices

The EENS and AENS are obtained and tabulated in Table 10 and Table 11 for all cases as illustrated in Figure 9a–d, respectively. It is important to note that the EENS and AENS decrease with the number of DGs, and these are also decreased with decreasing values of λ p and RT. As the increasing number of DGs are integrated into an EDS, the supplied energy is improved in the EDS, and thus, the indices related to the energy not supplied are reduced. This reduction is more while integrating the DGs with lesser λ p and RT values. The reducing EENS and AENS are desirable, and thus, the EDS reliability enhances with the integration of DGs with appropriate reliability data values.

5.3.2. Effect on System-Oriented Indices

The SAIDI, and SAIFI are obtained and tabulated in Table 12 and Table 13 for all cases considering all scenarios as shown in Figure 10a–d, respectively. The important point to be noted here that the SAIDI, and SAIFI, decreases with the increasing number of DGs; SAIDI is also decreased with the decreasing values of λ p and RT. It is worthy to note that the SAIFI is not affected by the RT of the DG. It is because this index is independent of RT as shown in Equation (43). As the increasing number of DGs are incorporated into the EDS, the duration of the interruption and the number of interruptions occurred are reduced in the system. Thus, the SAIDI, and SAIFI are reduced. Customer Average Interruption Duration Index can be determined using the ratio of SAIDI, and SAIFI as shown in Equation (45). The reduction in the value of indices is more while integrating the DGs with lesser λ p and RT values. Moreover, reducing SAIDI and SAIFI are desirable for EDS reliability enhancement.
The ASAI is determined and tabulated in Table 14 for all cases considering six scenarios as illustrated in Figure 11a,b. The electrical power service availability for all loads increases with the integration of multiple DGs. This index is further increased when DGs have a lower λ p and RT values. The increment in ASAI increases leads to the decrement in ASUI as given in Equations (46a) and (46c) which is desirable for the EDS reliability improvement.

6. Conclusions and Scope for Future Work

The optimal location(s) and size(s) of distributed generations derived for three cases ensure an enhancement in the electrical loss minimization and also, improves the bus voltage profile when compared to a system without distributed generations. While analyzing all the cases, it is observed that Case 3 has provided the best result. It was observed that the APL value is reduced at UPF by 0.00003 MW, 0.014472 MW, and 0.021307 MW for Case 1, Case 2, and Case 3, respectively when compared to [77]. The minimum value of bus voltage is improved by 6.42% for Case 3 at UPF when compared to No DG Case. The bus voltage profile was further improved at 0.85 and 0.82 power factors by 7.71% and 7.78%, respectively. After achieving the acceptable results, the reliability analysis is performed for the distributed generation integrated distribution system. It is concluded that the combination of Case 3 with Scenario 6 provides the best system’s reliability results. However, Scenario 6 being an ideal one, Case 3 with Scenario 1 is considered to yield better results for the improvement in the system’s reliability. This research analysis can be extended in the future for a greater number of standard distribution systems considering the following aspects.
  • Reliability assessment of larger EDS.
  • Inclusion of reliability data of subsystems.
  • System reconfiguration.
  • Considering CO 2 emission.
  • Economical aspects related to the system’s reliability, including net present value, Levelized cost of energy, and many other aspects.
  • Security.
The future aspects can be dealt with large radial distribution systems, including IEEE 69 and 118 bus systems. Inclusion of reliability data in terms of electrical, mechanical, and structural subsystems which would be favorable in obtaining accurate reliability of the system. Distribution system reliability improvement can also be achieved by adjusting the number of branches which is termed as a system reconfiguration technique. Furthermore, the dependency of reliability analysis on CO 2 emission, economics, system protection, and system security can also be considered for further research works.

Author Contributions

Conceptualization, Methodology, and Draft writing: S.K.; Software Implementation and Review: K.S.; Mathematical formulations: A.S.S.V.; Investigation, Results validation, and Data curation: R.M.E.; Supervision, Results validation, and Review: R.K.S.; Review and Editing: N.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors wholeheartedly thank the Department of Electrical Engineering, Indian Institute of Technology (Banaras Hindu University), Varanasi for providing the laboratory facilities and technical support to accomplish the research work promptly. The first author expresses his gratitude towards Govind Ballabh Pant Institute of Engineering and Technology, Pauri Garhwal, Uttarakhand, India for allowing him to pursue PhD from IIT(BHU) Varanasi, Uttar Pradesh, India.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

This work has used the following abbreviations:
AENSAverage Energy Not Supplied (MWh per customer per year)
APLActive Power Loss (MW)
ASAIAverage Service Availability Index (pu)
BSSBattery Storage System
CF-PSOConstriction Factor-based Particle Swarm Optimization
DGDistributed Generation
DGenDiesel Generator
EDSElectrical Distribution System
EENSExpected Energy Not Supplied (MWh per year)
EIREnergy index of reliability (pu)
ELMElectrical Loss Minimization
ENSEnergy Not Supplied (MWh)
GAGenetic Algorithm
GEGeneral Electric
IEEEInstitution of Electrical and Electronics Engineers
LOLELoss of Load Expectation (hour)
LOLPLoss of Load Probability (pu)
MOGAMulti-objective Genetic Algorithm
MWMegawatt
OPFOptimal Power Flow
pfPower Factor (pu)
RAReliability Assessment
RDSRadial Distribution System
RESRenewable Energy Source
RPLReactive Power Loss (MW)
SAIDISystem Average Interruption Duration Interruption (hour per customer per year)
SAIFISystem Average Interruption Frequency Interruption (failure per customer per year)
SPRSurface Plasmon Resonance
SPVSolar photovoltaic
WTGWind Turbine Generator
VPVoltage Profile

Appendix A

Appendix A.1

The IEEE 33 bus EDS as shown in Figure A1 has been selected for analysis. As mentioned in Section 5, it is inferred that the 33-bus distribution system is analyzed for obtaining the locations and sizes. The branch and load data have been taken from [76].
Figure A1. The IEEE 33 bus EDS.
Figure A1. The IEEE 33 bus EDS.
Energies 13 05631 g0a1

Appendix A.2

Table A1. Load distribution for 33 bus [52].
Table A1. Load distribution for 33 bus [52].
Type of Load
Bus Number (or Load Point)Number of LoadsMixedSame Type of Loads
2–5148Industrial (I)CIR
6–910Commercial (C)
11, 12132
13–15110Residential (R)
162
17–20118
21–26126
27–31108
32, 3358
Table A2. Reliability data for 33 bus [52].
Table A2. Reliability data for 33 bus [52].
Reliability Data for All Loads, Feeders, etc.
Bus, Feeder, etc. λ p (f/yr)RT (h)
Load@40.32111.04
Load@(5, 7–12, 29,0.30111.44
30, 14, 16, 18–22, 25–28)
Load@(13, 15)0.31411.17
Load@(17, 23, 24)0.2081.75
Load@(31–33)0.32710.96
substation0.15
feeder (2, 3, 6)0.23
DG0.212
From Table A2, average failure rate, and unavailability can be determined as given in Equations (A1) and (A2). The system reliability indices are thus evaluated for all six scenarios using the values of λ p and U p .
average failure rate , λ p = ( no . of loads × failure rates ) + ( no . of substations × failure rate ) + ( no . of feeders × failure rate ) + ( no . of DG × failure rate )
outage duration or unavailability , U p = ( no . of loads × failure rates × R T ) + ( no . of substations × failure rate × R T ) + ( no . of feeders × failure rate × R T ) + ( no . of DG × failure rate × R T )
EENS = [ ( Demand or load at P t h load point ) × ( Annual outage duration at P t h load point ) ]
AENS = ( EENS at P t h load point ) Total number of customers at all load points
SAIDI = [ ( annual outage duration at P t h load point ) × ( Number of customers at P t h load point ) ] Total number of customers at all load points
SAIFI = [ ( average failure rate at P t h load point ) × ( Number of customers at P t h load point ) ] Total number of customers at all load points
ASAI = 1 S A I D I 8760
Table A3. Cost per kilowatt for Reliability worth estimation [78].
Table A3. Cost per kilowatt for Reliability worth estimation [78].
Type of LoadInterruption Duration (minutes)Cost ($/kW)
10.38
202.97
Commercial608.55
24031.32
48083.01
11.63
203.87
Industrial609.09
24025.16
48055.81
10
200.09
Residential600.48
2404.91
48015.69

References

  1. Raju, K.; Madurai Elavarasan, R.; Mihet-Popa, L. An Assessment of Onshore and Offshore Wind Energy Potential in India Using Moth Flame Optimization. Energies 2020, 13, 3063. [Google Scholar]
  2. Elavarasan, R.M. The motivation for renewable energy and its comparison with other energy sources: A review. Eur. J. Sustain. Dev. Res. 2019, 3, em0076. [Google Scholar] [CrossRef]
  3. Kumar, N.M.; Chopra, S.S.; Chand, A.A.; Elavarasan, R.M.; Shafiullah, G. Hybrid renewable energy microgrid for a residential community: A techno-economic and environmental perspective in the context of the SDG7. Sustainability 2020, 12, 3944. [Google Scholar] [CrossRef]
  4. Elavarasan, R.M.; Shafiullah, G.; Padmanaban, S.; Kumar, N.M.; Annam, A.; Vetrichelvan, A.M.; Mihet-Popa, L.; Holm-Nielsen, J.B. A comprehensive review on renewable energy development, Challenges, and policies of leading Indian states with an international perspective. IEEE Access 2020, 8, 74432–74457. [Google Scholar] [CrossRef]
  5. Elavarasan, R.M.; Shafiullah, G.; Manoj Kumar, N.; Padmanaban, S. A State-of-the-Art review on the drive of renewables in Gujarat, State of India: Present situation, barriers and future initiatives. Energies 2020, 13, 40. [Google Scholar] [CrossRef] [Green Version]
  6. Madurai Elavarasan, R.; Selvamanohar, L.; Raju, K.; Vijayaraghavan, R.R.; Subburaj, R.; Nurunnabi, M.; Khan, I.A.; Afridhis, S.; Hariharan, A.; Pugazhendhi, R. A Holistic Review of the Present and Future Drivers of the Renewable Energy Mix in Maharashtra, State of India. Sustainability 2020, 12, 6596. [Google Scholar] [CrossRef]
  7. Tan, X.; Li, Q.; Wang, H. Advances and trends of energy storage technology in microgrid. Int. J. Electr. Power Energy Syst. 2013, 44, 179–191. [Google Scholar] [CrossRef]
  8. Elavarasan, R.M.; Afridhis, S.; Vijayaraghavan, R.R.; Subramaniam, U.; Nurunnabi, M. SWOT analysis: A framework for comprehensive evaluation of drivers and barriers for renewable energy development in significant countries. Energy Rep. 2020, 6, 1838–1864. [Google Scholar]
  9. Elavarasan, R.M. Comprehensive review on India’s growth in renewable energy technologies in comparison with other prominent renewable energy based countries. J. Sol. Energy Eng. 2020, 142, 030801. [Google Scholar] [CrossRef]
  10. Sharma, S.; Bhattacharjee, S.; Bhattacharya, A. Grey wolf optimisation for optimal sizing of battery energy storage device to minimise operation cost of microgrid. IET Gener. Transm. Distrib. 2016, 10, 625–637. [Google Scholar] [CrossRef]
  11. Energy Storage System: Roadmap for India 2019–32. 2019. Available online: http://www.indiasmartgrid.org/reports/ISGF_Report_Energy_Storage_System_RoadmapforIndia_2019to2032_11July2019_Draft.pdf (accessed on 29 September 2020).
  12. Hassan, T.; Abbassi, R.; Jerbi, H.; Mehmood, K.; Tahir, M.F.; Cheema, K.M.; Elavarasan, R.M.; Ali, F.; Khan, I.A. A Novel Algorithm for MPPT of an Isolated PV System Using Push Pull Converter with Fuzzy Logic Controller. Energies 2020, 13, 4007. [Google Scholar] [CrossRef]
  13. Ebrahimi, H.; Marjani, S.R.; Talavat, V. Optimal planning in active distribution networks considering nonlinear loads using the MOPSO algorithm in the TOPSIS framework. Int. Trans. Electr. Energy Syst. 2020, 30, e12244. [Google Scholar] [CrossRef]
  14. Liu, L.; Zhang, Y.; Da, C.; Huang, Z.; Wang, M. Optimal allocation of distributed generation and electric vehicle charging stations based on intelligent algorithm and bi-level programming. Int. Trans. Electr. Energy Syst. 2020, 30, e12366. [Google Scholar] [CrossRef]
  15. da Silva Seta, F.; de Oliveira, L.W.; de Oliveira, E.J. Comprehensive approach for distribution system planning with uncertainties. IET Gener. Transm. Distrib. 2019, 13, 5467–5477. [Google Scholar] [CrossRef]
  16. Ahmadi, M.; Lotfy, M.E.; Howlader, A.M.; Yona, A.; Senjyu, T. Centralised multi-objective integration of wind farm and battery energy storage system in real-distribution network considering environmental, technical and economic perspective. IET Gener. Transm. Distrib. 2019, 13, 5207–5217. [Google Scholar] [CrossRef]
  17. Deb, G.; Chakraborty, K.; Deb, S. Spider monkey optimization technique–based allocation of distributed generation for demand side management. Int. Trans. Electr. Energy Syst. 2019, 29, e12009. [Google Scholar] [CrossRef]
  18. Dehnavi, E.; Aminifar, F.; Afsharnia, S. Congestion management through distributed generations and energy storage systems. Int. Trans. Electr. Energy Syst. 2019, 29, e12018. [Google Scholar] [CrossRef]
  19. Chedid, R.; Sawwas, A. Optimal placement and sizing of photovoltaics and battery storage in distribution networks. Energy Storage 2019, 1, e46. [Google Scholar] [CrossRef] [Green Version]
  20. Babu, K.B.; Maheswarapu, S. A solution to multi-objective optimal accommodation of distributed generation problem of power distribution networks: An analytical approach. Int. Trans. Electr. Energy Syst. 2019, 29. [Google Scholar] [CrossRef]
  21. Hesaroor, K.; Das, D. Annual energy loss reduction of distribution network through reconfiguration and renewable energy sources. Int. Trans. Electr. Energy Syst. 2019, 29, e12099. [Google Scholar] [CrossRef]
  22. BiazarGhadikolaei, M.; Shahabi, M.; Barforoushi, T. Expansion planning of energy storages in microgrid under uncertainties and demand response. Int. Trans. Electr. Energy Syst. 2019, 29, e12110. [Google Scholar] [CrossRef]
  23. Amir, V.; Azimian, M.; Razavizadeh, A.S. Reliability-constrained optimal design of multicarrier microgrid. Int. Trans. Electr. Energy Syst. 2019, 29, e12131. [Google Scholar] [CrossRef]
  24. RaguRaman, L.; Ravindran, M. MFLRS-RDF technique for optimal sizing and performance analysis of HRES. Int. J. Numer. Model. Electr. Netw. Devices Fields 2020, 33, e2675. [Google Scholar] [CrossRef]
  25. Gholami, M.; Zakariazadeh, A. Olympic ranking–based allocation of distributed generation units in distribution networks. Int. Trans. Electr. Energy Syst. 2020, 30, e12220. [Google Scholar] [CrossRef]
  26. Muhammad, M.A.; Mokhlis, H.; Amin, A.; Naidu, K.; Franco, J.F.; Wang, L.; Othman, M. Enhancement of simultaneous network reconfiguration and DG sizing via Hamming dataset approach and firefly algorithm. IET Gener. Transm. Distrib. 2019, 13, 5071–5082. [Google Scholar] [CrossRef]
  27. Manna, D.; Goswami, S.K. Optimum placement of distributed generation considering economics as well as operational issues. Int. Trans. Electr. Energy Syst. 2020, 30, e12246. [Google Scholar] [CrossRef]
  28. Pandey, A.K.; Kirmani, S. Optimal location and sizing of hybrid system by analytical crow search optimization algorithm. Int. Trans. Electr. Energy Syst. 2020, 30, e12327. [Google Scholar] [CrossRef]
  29. Kiani, A.T.; Nadeem, M.F.; Ahmed, A.; Khan, I.; Elavarasan, R.M.; Das, N. Optimal PV Parameter Estimation via Double Exponential Function-Based Dynamic Inertia Weight Particle Swarm Optimization. Energies 2020, 13, 4037. [Google Scholar] [CrossRef]
  30. Kumar, S.; Saket, R.; Dheer, D.K.; Holm-Nielsen, J.; Sanjeevikumar, P. Reliability enhancement of electrical power system including impacts of renewable energy sources: A comprehensive review. IET Gener. Transm. Distrib. 2020, 14, 1799–1815. [Google Scholar] [CrossRef]
  31. Samrout, M.; Yalaoui, F.; Châtelet, E.; Chebbo, N. New methods to minimize the preventive maintenance cost of series–parallel systems using ant colony optimization. Reliab. Eng. Syst. Saf. 2005, 89, 346–354. [Google Scholar] [CrossRef]
  32. Shahzad, M.; Ahmad, I.; Gawlik, W.; Palensky, P. Load concentration factor based analytical method for optimal placement of multiple distribution generators for loss minimization and voltage profile improvement. Energies 2016, 9, 287. [Google Scholar] [CrossRef] [Green Version]
  33. Parihar, S.S.; Malik, N. Optimal allocation of renewable DGs in a radial distribution system based on new voltage stability index. Int. Trans. Electr. Energy Syst. 2020, 30, e12295. [Google Scholar] [CrossRef]
  34. Hassan, A.A.; Fahmy, F.H.; Nafeh, A.E.S.A.; Abu-elmagd, M.A. Hybrid genetic multi objective/fuzzy algorithm for optimal sizing and allocation of renewable DG systems. Int. Trans. Electr. Energy Syst. 2016, 26, 2588–2617. [Google Scholar] [CrossRef]
  35. Atteya, I.I.; Ashour, H.; Fahmi, N.; Strickland, D. Radial distribution network reconfiguration for power losses reduction using a modified particle swarm optimisation. CIRED Open Access Proc. J. 2017, 2017, 2505–2508. [Google Scholar] [CrossRef] [Green Version]
  36. Wazir, A.; Arbab, N. Analysis and optimization of IEEE 33 bus radial distributed system using optimization algorithm. JETAE J. Emerg. Trends Appl. Eng. 2016, 1, 2518–4059. [Google Scholar]
  37. Paliwal, N.K.; Singh, A.K.; Singh, N.K.; Kumar, P. Optimal sizing and operation of battery storage for economic operation of hybrid power system using artificial bee colony algorithm. Int. Trans. Electr. Energy Syst. 2019, 29, e2685. [Google Scholar] [CrossRef] [Green Version]
  38. Saric, M.; Hivziefendic, J.; Konjic, T.; Ktena, A. Distributed generation allocation considering uncertainties. Int. Trans. Electr. Energy Syst. 2018, 28, e2585. [Google Scholar] [CrossRef]
  39. Baran, M.; Wu, F. Network reconfiguration in distribution systems for loss reduction and load balancing. IEEE Trans. Power Deliv. 1989, 4, 1401–1407. [Google Scholar] [CrossRef]
  40. Chong, B.; Zhang, X.; Godfrey, K.; Yao, L.; Bazargan, M. Optimal location of unified power flow controller for congestion management. Eur. Trans. Electr. Power 2010, 20, 600–610. [Google Scholar] [CrossRef]
  41. Aghajani, A.; Kazemzadeh, R.; Ebrahimi, A. Optimal energy storage sizing and offering strategy for the presence of wind power plant with energy storage in the electricity market. Int. Trans. Electr. Energy Syst. 2018, 28, e2621. [Google Scholar] [CrossRef]
  42. Charfi, S.; Atieh, A.; Chaabene, M. Optimal sizing of a hybrid solar energy system using particle swarm optimization algorithm based on cost and pollution criteria. Environ. Prog. Sustain. Energy 2019, 38, e13055. [Google Scholar] [CrossRef]
  43. Wu, C.; Lou, Y.; Lou, P.; Xiao, H. DG location and capacity optimization considering several objectives with cloud theory adapted GA. Int. Trans. Electr. Energy Syst. 2014, 24, 1076–1088. [Google Scholar] [CrossRef]
  44. Duong, M.Q.; Pham, T.D.; Nguyen, T.T.; Doan, A.T.; Tran, H.V. Determination of optimal location and sizing of solar photovoltaic distribution generation units in radial distribution systems. Energies 2019, 12, 174. [Google Scholar] [CrossRef] [Green Version]
  45. Madurai Elavarasan, R.; Ghosh, A.; K Mallick, T.; Krishnamurthy, A.; Saravanan, M. Investigations on performance enhancement measures of the bidirectional converter in PV–wind interconnected microgrid system. Energies 2019, 12, 2672. [Google Scholar] [CrossRef] [Green Version]
  46. Natarajan, M.; Ramadoss, B.; Lakshmanarao, L. Optimal location and sizing of MW and MVAR based DG units to improve voltage stability margin in distribution system using a chaotic artificial bee colony algorithm. Int. Trans. Electr. Energy Syst. 2017, 27, e2287. [Google Scholar] [CrossRef]
  47. Nawaz, S.; Bansal, A.K.; Sharma, M.P. Allocation of DG and capacitor units for power loss reduction in radial distribution system. In Proceedings of the 2016 International Conference on Recent Advances and Innovations in Engineering (ICRAIE), Jaipur, India, 23–25 December 2016; pp. 1–6. [Google Scholar]
  48. Vita, V. Development of a decision-making algorithm for the optimum size and placement of distributed generation units in distribution networks. Energies 2017, 10, 1433. [Google Scholar] [CrossRef] [Green Version]
  49. Basso, T. IEEE 1547 and 2030 Standards for Distributed Energy Resources Interconnection and Interoperability with the Electricity Grid; Technical report; National Renewable Energy Lab. (NREL): Golden, CO, USA, 2014. [Google Scholar]
  50. Huang, W.; Zhang, N.; Yang, J.; Wang, Y.; Kang, C. Optimal configuration planning of multi-energy systems considering distributed renewable energy. IEEE Trans. Smart Grid 2017, 10, 1452–1464. [Google Scholar] [CrossRef]
  51. Bahramirad, S.; Daneshi, H. Optimal Sizing of Smart Grid Storage Management System in a Microgrid. In Proceedings of the 2012 IEEE PES Innovative Smart Grid Technologies, ISGT 2012, Washington DC, USA, 16–20 January 2012; pp. 1–6. [Google Scholar]
  52. Kumar, D.; Samantaray, S.; Kamwa, I.; Sahoo, N. Reliability-constrained based optimal placement and sizing of multiple distributed generators in power distribution network using cat swarm optimization. Electr. Power Compon. Syst. 2014, 42, 149–164. [Google Scholar] [CrossRef]
  53. Soliman, S.A.H.; Mantawy, A.A.H. Modern Optimization Techniques with Applications in Electric Power Systems; Springer Science & Business Media: New York, NY, USA, 2011. [Google Scholar]
  54. Bhumkittipich, K.; Phuangpornpitak, W. Optimal placement and sizing of distributed generation for power loss reduction using particle swarm optimization. Energy Procedia 2013, 34, 307–317. [Google Scholar] [CrossRef] [Green Version]
  55. Aman, M.; Jasmon, G.; Mokhlis, H.; Bakar, A. Optimal placement and sizing of a DG based on a new power stability index and line losses. Int. Trans. Electr. Energy Syst. 2012, 43, 1296–1304. [Google Scholar] [CrossRef]
  56. Shukla, T.; Singh, S.; Srinivasarao, V.; Naik, K. Optimal sizing of distributed generation placed on radial distribution systems. Electr. Power Compon. Syst. 2010, 38, 260–274. [Google Scholar] [CrossRef]
  57. Abdel-Aal, H.; Bassyouni, M.; Abdelkreem, M.; Abdel-Hamid, S.; Zohdy, K. Feasibility Study for a Solar-Energy Stand-Alone System:(SESAS). Smart Grid Renew. Energy 2012, 3, 204. [Google Scholar] [CrossRef] [Green Version]
  58. Hamouda, Y.A. Wind energy in Egypt: Economic feasibility for Cairo. Renew. Sustain. Energy Rev. 2012, 16, 3312–3319. [Google Scholar] [CrossRef] [Green Version]
  59. Zhao, B.; Guo, C.; Cao, Y. A multiagent-based particle swarm optimization approach for optimal reactive power dispatch. IEEE Trans. Power Syst. 2005, 20, 1070–1078. [Google Scholar] [CrossRef]
  60. Kusiak, A.; Zhang, Z.; Xu, G. Minimization of wind farm operational cost based on data-driven models. IEEE Trans. Sustain. Energy 2013, 4, 756–764. [Google Scholar] [CrossRef] [Green Version]
  61. Wang, L.; Singh, C. Multicriteria design of hybrid power generation systems based on a modified particle swarm optimization algorithm. IEEE Trans. Energy Convers. 2009, 24, 163–172. [Google Scholar] [CrossRef]
  62. Ramezani, M.; Haghifam, M.R.; Singh, C.; Seifi, H.; Moghaddam, M.P. Determination of capacity benefit margin in multiarea power systems using particle swarm optimization. IEEE Trans. Power Syst. 2008, 24, 631–641. [Google Scholar] [CrossRef]
  63. Yang, H.; Xie, K.; Tai, H.M.; Chai, Y. Wind farm layout optimization and its application to power system reliability analysis. IEEE Trans. Power Syst. 2015, 31, 2135–2143. [Google Scholar] [CrossRef]
  64. Sunil Joseph, P.; DineshBalaji, C. Transmission loss minimization using optimization technique Based on PSO. 2011 IEEE Symp. Ind. Electron. Appl. 2013, 1, 1–5. [Google Scholar] [CrossRef]
  65. Clerc, M.; Kennedy, J. The particle swarm-explosion, stability, and convergence in a multidimensional complex space. IEEE Trans. Evol. Comput. 2002, 6, 58–73. [Google Scholar] [CrossRef] [Green Version]
  66. El-Zonkoly, A. Optimal placement of multi-distributed generation units including different load models using particle swarm optimisation. IET Gener. Transm. Distrib. 2011, 5, 760–771. [Google Scholar] [CrossRef]
  67. Kansal, S.; Kumar, V.; Tyagi, B. Optimal placement of different type of DG sources in distribution networks. Int. J. Electr. Power Energy Syst. 2013, 53, 752–760. [Google Scholar] [CrossRef]
  68. Kansal, S.; Sai, B.; Tyagi, B.; Kumar, V. Optimal placement of distributed generation in distribution networks. Int. J. Eng. Sci. Technol. 2011, 3, 47–55. [Google Scholar] [CrossRef]
  69. Subcommittee, D. IEEE Guide for Electric Power Distribution Reliability Indices. Distribution 2012, 1997, 1–43. [Google Scholar]
  70. Heier, S. Grid Integration of Wind Energy: Onshore and Offshore Conversion Systems; John Wiley & Sons: Sussex, UK, 2014. [Google Scholar]
  71. Koutroulis, E.; Kolokotsa, D.; Potirakis, A.; Kalaitzakis, K. Methodology for optimal sizing of stand-alone photovoltaic/wind-generator systems using genetic algorithms. Sol. Energy 2006, 80, 1072–1088. [Google Scholar] [CrossRef]
  72. van Sark, W.G. Teaching the relation between solar cell efficiency and annual energy yield. Eur. J. Phys. 2007, 28, 415. [Google Scholar] [CrossRef]
  73. Singh, R.; Sharma, M.; Rawat, R.; Banerjee, C. An assessment of series resistance estimation techniques for different silicon based SPV modules. Renew. Sustain. Energy Rev. 2018, 98, 199–216. [Google Scholar] [CrossRef]
  74. Yang, J.; Sun, Y.; Xu, Y. Modeling impact of environmental factors on photovoltaic array performance. Int. J. Energy Environ. 2013, 4, 955–968. [Google Scholar]
  75. Gunawardana, A. Proper sizing of energy storage for grid connected photovoltaic system. Master’s Thesis, Department of Engineering Science Faculty of Engineering and Science, University of Agder, Kristiansand, Norway, 2014. [Google Scholar]
  76. Sultana, S.; Roy, P.K. Multi-objective quasi-oppositional teaching learning based optimization for optimal location of distributed generator in radial distribution systems. Int. J. Electr. Power Energy Syst. 2014, 63, 534–545. [Google Scholar] [CrossRef]
  77. Mahmoud, K.; Yorino, N.; Ahmed, A. Optimal distributed generation allocation in distribution systems for loss minimization. IEEE Trans. Power Syst. 2015, 31, 960–969. [Google Scholar] [CrossRef]
  78. Billington, R.; Allan, R.N. Reliability Evaluation of Power Systems; Springer: New York, NY, USA, 1996. [Google Scholar]
Figure 1. Workflow of the research work.
Figure 1. Workflow of the research work.
Energies 13 05631 g001
Figure 2. Models of several DG technologies.
Figure 2. Models of several DG technologies.
Energies 13 05631 g002
Figure 3. General 2-bus system to formulate the line loss and load factor.
Figure 3. General 2-bus system to formulate the line loss and load factor.
Energies 13 05631 g003
Figure 4. 6-bus EDS for Index-1 calculation.
Figure 4. 6-bus EDS for Index-1 calculation.
Energies 13 05631 g004
Figure 5. Algorithm implemented for the research work.
Figure 5. Algorithm implemented for the research work.
Energies 13 05631 g005
Figure 6. General BSS dispatch strategy opted for battery operation.
Figure 6. General BSS dispatch strategy opted for battery operation.
Energies 13 05631 g006
Figure 7. Voltage profile for 33 bus system (WTG at 0.9 pf).
Figure 7. Voltage profile for 33 bus system (WTG at 0.9 pf).
Energies 13 05631 g007
Figure 8. Active and Reactive power losses (WTG at UPF).
Figure 8. Active and Reactive power losses (WTG at UPF).
Energies 13 05631 g008
Figure 9. Load-oriented indices for different DG reliability data (a) EENS at different λ p , (b) EENS at different RT, (c) AENS at different λ p , and (d) AENS at different RT.
Figure 9. Load-oriented indices for different DG reliability data (a) EENS at different λ p , (b) EENS at different RT, (c) AENS at different λ p , and (d) AENS at different RT.
Energies 13 05631 g009
Figure 10. System-oriented indices for different DG reliability data (a) SAIDI at different λ p , (b) SAIDI at different RT, (c) SAIFI at different λ p , and (d) SAIFI at different RT.
Figure 10. System-oriented indices for different DG reliability data (a) SAIDI at different λ p , (b) SAIDI at different RT, (c) SAIFI at different λ p , and (d) SAIFI at different RT.
Energies 13 05631 g010
Figure 11. System-oriented indices for different DG reliability data (a) ASAI at different λ p , (b) ASAI at different RT.
Figure 11. System-oriented indices for different DG reliability data (a) ASAI at different λ p , (b) ASAI at different RT.
Energies 13 05631 g011
Table 1. Published work.
Table 1. Published work.
No. of BusParameters Considered
SizeLocationVoltageLossReliabilityPower FactorDGReference
34, 69 PV/WTG[33]
33, 69, 119 PV/WTG[34]
12 DGen[38]
33, 69 PV[44]
13 PV/BSS[19]
33, 118 PV/WTG[21]
38 PV/WTG[27]
33, 69 PV/ESS[28]
38, 69 GT/WTG[46]
69, 118 DGen/Cap[47]
33 PV/WTG[48]
33SPV/WTG/BSSPresent Work
Table 2. Work related to IEEE 33 bus with multiple DGs.
Table 2. Work related to IEEE 33 bus with multiple DGs.
Method# of DGDG PositionTotal DG Size (MW)Loss (MW)Reference
MOGA1(SPV)81.63330.113[34]
2(SPV)14, 300.8337, 0.998510.08435
1(WTG)81.850.08556
2(WTG)14, 301.1, 0.750.04791
GA3(SPV)14, 24, 280.6947, 1.1844, 1.46280.0756[44]
ABC3(SPV)9, 24, 321.1372, 1.0674, 0.80310.0752
PSO3(SPV)9, 24, 301.0625, 1.0447, 0.95180.0744
BBO3(SPV)14, 24, 300.7539, 1.0994, 1.07140.0715
CSO5(BSS)1, 4, 11, 12, 180.15, 0.4117, 0.6705, 0.1, 8.90550.02379[28]
1(SPV)620.0908
DMA1810.1175[48]
1(WTG)331.650.1068
Table 3. Values of indexes with corresponding bus number [17,20].
Table 3. Values of indexes with corresponding bus number [17,20].
Index-1Index-2
Value Bus No. Value Bus No.
1.350 × 10 3 6 41.52 × 10 3 30
0.928 × 10 3 29 16.44 × 10 3 13
0.867 × 10 3 30 16.43 × 10 3 24
0.864 × 10 3 5 7.36 × 10 3 31
0.735 × 10 3 28 6.49 × 10 3 20
Table 4. Wind Turbine (V162-5.6 MW) specifications.
Table 4. Wind Turbine (V162-5.6 MW) specifications.
ParameterRating (Unit)
Rated output power5.6 MW
Cut-in Speed ( V c i n )3 m/s
Cut-out Speed ( V c o u t )25 m/s
Temperature 20 C to 45 C
Diameter162 m
Swept Area20612 m 2
Frequency50/60 Hz
Hub Height119 m, 125 m, 148 m, 149 m, and 166 m
Table 5. Bifacial Solar Panel (SPR-P5-545-UPP) specifications.
Table 5. Bifacial Solar Panel (SPR-P5-545-UPP) specifications.
ParameterRating (Unit)ParameterRating (Unit)
nominal power545 WMaximum Series Fuse25 A
Tolerance of Power ± 3 / 0 % Temperature−40–85  C
Efficiency21.1%Power Temperature Coefficient 0.34 % / C
Rated voltage46.1 VVoltage Temperature Coefficient 0.28 % / C
Rated current11.84 ACurrent Temperature Coefficient 0.06 % / C
Open circuit voltage55.8 VWeight31.5 kg
Short circuit current12.62 ASolar CellsMono-crystalline
Maximum System Voltage (IEC)1500 VL × B × H mm 3 2362 × 1092 × 35
Table 6. Battery Storage System (BESS 3000) specifications.
Table 6. Battery Storage System (BESS 3000) specifications.
ParameterRating
Rated output power3000 kW
Storage Capacity1000 kWh
Rated output current2795 A
Rated output AC voltage620 V
Power factor0.95 Cap …0.95 Ind
Total harmonic distortion<3%
Efficiency>98%
TypeLithium-ion
IGBT Switching Frequency (Converter)2–4 kHz
Table 7. DG location and DG size obtained.
Table 7. DG location and DG size obtained.
Case 1Case 2Case 3
Location (bus no.)630, 1330, 13, 24
size@upf (MW)2.5641.148, 0.8431.048, 0.801, 1.105
Table 8. APL (MW) obtained considering WTG power factor.
Table 8. APL (MW) obtained considering WTG power factor.
pfNo DGCase 1Case 2Case 3
Unity 0.111040.07270.05148
Present work0.850.211010.068310.045390.02795
0.82 0.068310.04440.02702
APL obtained considering power factor of all Conventional DGs
pfNo DGSingle DGTwo DGThree DG
Unity 0.111070.0871720.072787
EA [77]0.850.2110.0681700.031190.01552
0.82 0.0678700.030410.01514
Table 9. Minimum voltage, DG location, and RPL (MVAr) obtained.
Table 9. Minimum voltage, DG location, and RPL (MVAr) obtained.
Minimum Voltage (%)RPL
pf No DG Case 1 Case 2 Case 3 No DG Case 1 Case 2 Case 3
Unity 94.2696.8896.86 0.081680.051210.03848
0.8590.4495.7498.1298.150.143060.055040.032570.02185
0.82 96.098.2098.22 0.055040.031950.02119
Table 10. EENS (MWh per year) evaluated for different Scenarios.
Table 10. EENS (MWh per year) evaluated for different Scenarios.
Scenario 1Scenario 2Scenario 3Scenario 4Scenario 5Scenario 6
No DG82.76382.76382.76382.76382.76382.763
Case 165.53368.46573.39768.46578.32958.601
Case 231.81729.24931.81729.24934.38524.113
Case 330.13527.56730.13527.56732.70322.431
Table 11. AENS (MWh per customer per year) evaluated for different Scenarios.
Table 11. AENS (MWh per customer per year) evaluated for different Scenarios.
Scenario 1Scenario 2Scenario 3Scenario 4Scenario 5Scenario 6
No DG0.02550.02550.02550.02550.02550.0255
Case 10.01960.02030.02260.02110.02410.0181
Case 20.00820.0090.00980.0090.01060.0074
Case 30.00770.00850.00930.00850.01010.0069
Table 12. SAIDI (hour per customer per year) evaluated for different Scenarios.
Table 12. SAIDI (hour per customer per year) evaluated for different Scenarios.
Scenario 1Scenario 2Scenario 3Scenario 4Scenario 5Scenario 6
No DG24.01224.01224.01224.01224.01224.012
Case 118.76420.08521.40620.08522.72817.442
Case 27.3888.0538.7198.0539.3856.722
Case 37.2017.8668.5327.8669.1986.535
Table 13. SAIFI (failure per customer per year) evaluated for different Scenarios.
Table 13. SAIFI (failure per customer per year) evaluated for different Scenarios.
Scenario 1Scenario 2Scenario 3Scenario 4Scenario 5Scenario 6
No DG3.1793.1793.1793.1793.1793.179
Case 12.1092.2192.3292.1092.1091.999
Case 20.9150.9701.0260.9150.9150.859
Case 30.8420.8970.9530.8420.8420.786
Table 14. ASAI (pu) evaluated for different Scenarios.
Table 14. ASAI (pu) evaluated for different Scenarios.
Scenario 1Scenario 2Scenario 3Scenario 4Scenario 5Scenario 6
No DG0.997260.997260.997260.997260.997260.99726
Case 10.997860.997710.997560.997710.997410.99801
Case 20.999160.999080.999000.999080.998930.99923
Case 30.999180.999100.999030.999100.998950.99925
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kumar, S.; Sarita, K.; Vardhan, A.S.S.; Elavarasan, R.M.; Saket, R.K.; Das, N. Reliability Assessment of Wind-Solar PV Integrated Distribution System Using Electrical Loss Minimization Technique. Energies 2020, 13, 5631. https://doi.org/10.3390/en13215631

AMA Style

Kumar S, Sarita K, Vardhan ASS, Elavarasan RM, Saket RK, Das N. Reliability Assessment of Wind-Solar PV Integrated Distribution System Using Electrical Loss Minimization Technique. Energies. 2020; 13(21):5631. https://doi.org/10.3390/en13215631

Chicago/Turabian Style

Kumar, Sachin, Kumari Sarita, Akanksha Singh S Vardhan, Rajvikram Madurai Elavarasan, R. K. Saket, and Narottam Das. 2020. "Reliability Assessment of Wind-Solar PV Integrated Distribution System Using Electrical Loss Minimization Technique" Energies 13, no. 21: 5631. https://doi.org/10.3390/en13215631

APA Style

Kumar, S., Sarita, K., Vardhan, A. S. S., Elavarasan, R. M., Saket, R. K., & Das, N. (2020). Reliability Assessment of Wind-Solar PV Integrated Distribution System Using Electrical Loss Minimization Technique. Energies, 13(21), 5631. https://doi.org/10.3390/en13215631

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