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

: 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 reconﬁguration, 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, ﬁrst, 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.

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.

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.

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.

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].

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.

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). Line2 Line3 Line4 Then the bus voltage matrix is formed by using Equation (3) By expanding Equation (3), the node voltages can be obtained by using Equations (4) and (5).
where i is 1,2,. . . ,n bus . Z ai , Z bi and I i are the element of impedance matrix that signify the a th row and i th column, b th row and i th column, and current injection at bus-i, respectively.
Current Injection, where P in,i and jQ in,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).
For an electrical system with n l number of branch/lines, the line loss is given by (8).
where B LOSS l is line loss, S in,i is apparent power injected at bus-i, l is 1 to n l . [ where load factor (LF) is given by . LF li is load factor of the lth line due to the ith bus injection (it is non-zero if lth line is in the path of ith 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).   (13) is reduced according to the Figure 4 where power injection at all the buses except the source bus are available.

r a c k er-s of t w a r e . c o m
The calculation of effective power injections is done as given in (16).
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 th bus from the source bus can be observed in the LF li as guided in Equation (10). If the i th bus is not near, the number of lines between i th bus and source bus is being plenty and the corresponding Z ai , Z bi , and Y l parameters will account in the electrical loss component. Furthermore, if the node voltage is high, the value of LF li 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 th 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.
S in e f f i is the effective injection of complex power, which is the sum of injected powers from other buses connected to ith bus as shown in Figure 4.
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.

Power Balance
The AP and RP balance expressions are shown in Equations (19) and (20).
where P net i = 0 and Q net i = 0 are the net AP and RP at i-bus, respectively. P dg i and Q dg i represent DG AP and RP at i-bus, respectively. Active and reactive load demands are mentioned by P dem i and Q dem 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.

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).

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).
where P real i , P real j , Q real i , Q real j are the AP and RP at i and j-buses, respectively. N bus = number of buses or nodes, C 1 ij and C 2 ij are defined as follows.
where V i , δ i and V j , δ j are the voltages and corresponding angles at ith and jth buses, respectively, R ij = resistance of a branch between i and j-buses.

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).
where Q gen i and Q dem i are the RP generation and demand at the ith bus (including the slack bus), respectively. Q dem i = RP demand at the ith bus.

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].

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].

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.
where AP gen i = AP generated by the generation units at ith bus, AP dem i = AP demand at ith bus.
where RP gen i = RP generated by the generation units at ith bus, RP dem i = RP demand at ith bus.

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].
where AP max a ij = the highest permissible apparent powers (AP a ) for lines i to j, AP a ij = the actual AP 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].
where p = 1, 2,. . . . . . , n DG , AP min DG p (set to zero) and AP max DG p (from Equation (28)) are the lower and upper AP outputs of DG unit p, respectively. RP min DG p (set to zero) and RP max DG p (from Equation (29)) are the lower and upper RP outputs of DG unit p, respectively. n DG = 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].
where |V i minimum | and |V i maximum | = the lower and upper boundaries of the bus voltage |V i | which are set to 95% and 105%, respectively.

D. Branch current
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].

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).
where V p+1 n = n th particle velocity at (p + 1) th iteration, W' = particle inertial weight, V p n = n th particle velocity at p th iteration, C 1 , C 2 = constants (0, 2.5), R 1 , R 2 = numbers generated randomly (0, 1), Personal BEST i = the n th particle's best position considering its own property, Global BEST = the n th particle's best position considering the whole population, X p+1 n , X p n = n th particle position at (p + 1) th and p th 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.

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).

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 num k × F k failure per year (36) Outage duration (annual); U p = ∑ k=n F k D pk hour per year (37) Outage duration (average); D p = U p λ p hour (38) where F k = failure rate (average) of the kth element, n = number of elements in the EDS, num k = number of kth elements in the EDS, D pk = duration of failure at p th load point due to k th failed element. The calculation equations of λ p and U p are given in Appendix A.2.

System-Based Indices
These indices are further categorized in load-oriented indices and customer oriented indices as given in Equations (39) where P p = demand/load (average) of the pth load point, EENS p = expected ENS at Pth load or customer point. where n p = total load points, N p = number of customers at Pth load point.

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.
where ECOST p = expected interrupted cost at Pth load point, IEAR p = interrupted energy assessment rate at Pth load point, LOEE p = loss of expected energy, f (D pk ) = system composite customer damage function ($ per kilowatt) as provided in Also, ASAI can be derived as follows.
where ASUI = Average Service Unavailability Index.

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 Sections 4.1-4.3.

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 mech ) is a function of generator rotor speed and wind speed as formulated in Equations (47)-(51) [70].
where ρ = air density, A s = area swept by the turbine rotor blades, Vwind = speed of the wind, C p is the non-linear function of the tip speed ratio (λ) and pitch angle (θ), ω rotor = generator rotor speed, GR = gear ratio, R rotor = rotor radius at the turbine blades, C 1 to C 6 and x are constants and computed in [70].
where P WTG = output WTG power, P WTG,rated = rated output WTG power, V rated = rated wind speed.

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].
where r s,STC = normalized series resistance under STC, and all other parameters in Equations (61)- (65) are evaluated at STC.

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].
where P Battery(DC) (t) = DC charging/discharging power of the battery in ∆t interval (W), E Battery (t) = energy of battery (Wh), P Battery(AC) (t) = AC power discharged/charged state of battery, η Battery = efficiency of the battery, PF Inverter = inverter power factor, SOC(t) = state of charge of the battery, η Charging = charging efficiency of the battery, C Battery = battery capacity (Ah), V = nominal voltage of the battery (V).

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 Tables 4-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.

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 lth line depends on whether the lth is in the path of the ith bus to the source node or not. The multiplication of Load Factor and injected apparent power provides the APL and RPL of lth line due to the ith bus AP a injection. Thus, the maximum value of Index-1 at ith 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.

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.

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 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.

Effect on Load-Oriented Indices
The EENS and AENS are obtained and tabulated in Tables 10 and 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.

Effect on System-Oriented Indices
The SAIDI, and SAIFI are obtained and tabulated in Tables 12 and 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.  Figure 11. System-oriented indices for different DG reliability data (a) ASAI at different λ p , (b) ASAI at different RT.

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.

•
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.

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.