The objective of the present paper is to study the optimum installation of Non-dispatchable Distributed Generations (NDG) in the distribution network of given sizes under the given scheme. The uncertainty of various random (uncertain) parameters like load, wind and solar operated DG besides uncertainty of fuel prices has been investigated by the three-point estimate method (3-PEM) and Monte Carlo Simulation (MCS) based methods. Nearly twenty percent of the total number of buses are selected as candidate buses for NDG placement on the basis of system voltage profile to limit the search space. Weibull probability density function (PDF) is considered to address uncertain characteristics of solar radiation and wind speed under different scenarios. Load uncertainty is described by Standard Normal Distribution Function (SNDF). To investigate the solution of optimal probabilistic load flow (OPLF) three-point PEM-based technique was applied. For optimization, Genetic Algorithm (GA), Particle Swarm Optimization (PSO) and GA-PSO hybrid-based Artificial Intelligent (AI) based optimization techniques are employed to achieve the optimum value of the multi-objectives function. The proposed multi-objective function comprises loss and different costs. The proposed methods have been applied to IEEE 33- bus radial distribution network. Simulation results obtained by these techniques are compared based on loss minimization capability, enhancement of system bus voltage profile and reduction of cost and fitness functions. The major findings of the present study are the PEM-based method which provides almost similar results as MCS based method with less computation time and as far as loss minimization capacity, voltage profile improvement etc. is concerned, the hybrid-based optimization methods are compared with GA and PSO based optimization techniques.
Keywords: Non-dispatchable Distributed Generations (NDG); probabilistic load flow (PLF); Monte Carlo Simulation (MCS); point estimate method (PEM); GA-PSO based hybrid approach; optimization
In recent years distributed generations (DGs) are becoming very significant in the area of power engineering because DGs pose enumerable benefits over centralized generating units. Therefore, this research area is becoming more popular for active research. Nowadays, efforts are being made to cut down the line losses and emission of harmful greenhouse gases (GHG). Moreover, DG also improves system voltage profile and also improves system reliability. Additionally, DG units, which are powered by renewable energy sources (RES) such as wind turbine generators (WTG), solar photovoltaic (SPV), biomass, small hydro, and geothermal are extremely useful. Other types of DG technologies that are powered by non-renewable energy sources also exist, such as internal combustion engines, fuel cells etc. Most parameters like electric load demand, power generation of DG, and fuel prices suffer from uncertainty. Therefore, uncertainty is to be investigated thoroughly, which is associated to solar and wind operated DG and load demand. Further, it is established that DGs are universally accepted sustainable power generating units because of their inexhaustible and nonpolluting characteristics. Therefore, optimum placement and sizing of DG in the distribution system is extremely important in order to obtain maximum potential benefits [[
Consequently, in view of these facts, optimum placement of DG (OPDG) is essential. Furthermore, most of the researchers have not addressed the uncertainty of load and power generating units in their research. Hence, OPDG is very much crucial when uncertainties of load, fuel prices, and power generating units like wind turbine generators (WTG) and solar photovoltaic (SPV) are adequately investigated. The MCS-based method is commonly used for addressing the uncertainties of various parameters [[
In the present study, load demand is a random variable. The variation of load demand is considered as a Gaussian distribution function. Further, wind speed and solar insolation are considered as Weibull distribution functions. Further, fuel prices are considered as an uncertain parameter and follow the properties of Geometric Brownian Motion (GBM). To address the uncertainty of these parameters, three PEM-based techniques were applied and the obtained results give the MCS-based technique. The MCS-based method is the benchmark method for the assessment of uncertainty. For optimization GA, PSO and GA-PSO hybridized based optimization techniques are applied. The present research paper is dedicated to the investigation of uncertainty of NDG (SPV and WTG) and fueled-based DG is also considered for installation along with NDG. The main motive of the present study is optimal sizing and siting of non-dispatchable DG in distribution network by considering uncertainty of load demand addressed by Standard Gaussian PDF, and uncertainty of power generation by SPV and WTG by Weibull based PDF. Uncertainty associated in these parameters needs further attention to address. Therefore, the main focus of this study is to minimize the system power loss by placement of NDG in distribution system.
The power generated by solar-based DG and wind-operated DG is an uncertain parameter and has to be addressed carefully. The uncertainty of solar radiation and wind speed is investigated by PEM-based technique. Many researchers have addressed uncertainty of load and power generated by SPV and WTG-based DG. Evangelopoulos and Georgilakis [[
Authors of [[
It is a well-known fact that uncertain inputs give rise to uncertain outputs. Therefore, the methods which are applied to placement of DG units whose output is controllable are not suitable for placement of DG units whose output is intermittent [[
Moreover, uncertainty of non-dispatchable type DG is a crucial feature which needs to be investigated critically as it is evident from available literature [[
Therefore, from the above literature review, none of the research articles address the uncertainty of power generated from wind, solar operated DG, load, and fuel cost. In present manuscript uncertainty of wind and solar operated DG, along with load demand fuel prices uncertainty is investigated by 3-points PEM and MCS based technique. The cost function or objective function is designed by considering different costs that pertain to DG operation under different scenarios. Further, optimum locations and optimum value of cost function is obtained by GA, PSO, and GA-PSO-based hybrid technique.
From the literature it is noticed that very few researchers have addressed the uncertainty of three parameters like load demand and power generated by SPV and WTG based DG by PEM and MCS based method. For optimization of objective function GA, PSO, and GA-PSO-based hybridized methods are not considered. The fueled operated DG is also considered as dispatchable DG for placement.
From the obtained result, it is confirmed that the GA-PSO-based hybrid method is more feasible and provides an optimal solution close to actual results. Moreover, 3-point PEM-based method yields almost similar results as MCS based technique which is the only available technique to investigate the uncertainty of different parameters.
In the present article, three different scenarios and types of DGs (wind, solar, and fueled) have been considered. For optimization GA, PSO, and hybrid (combination of GA and PSO) based heuristic techniques have been applied. The PEM and MCS-based methods have been applied to incorporate the stochastic nature of different parameters. The configuration and properties of computer on which MATLAB programs have been executed is Intel (R) Core (TM) i 5-4570 CPU @ 3.20 GHz 4 GB RAM.
The remaining part of present research article is organized as follows: Section 2 discusses modeling of uncertainties of wind and solar operated DG. Section 3 describes the PEM-based method. Section 4 deals with the problem formulation of present study. Section 5 highlights the scheme and different scenarios of DGs which are to be placed along with solution process. Section 6 discusses the algorithm for selection of initial population of optimization technique. Finally, Section 7 describes the results and discussion.
As solar radiation and wind speed are uncertain parameters, the uncertainty of wind speed and solar radiation followed the Weibull PDF [[
Weibull probability density function (PDF) is generally applied to represent the stochastic nature of wind speed. Many researchers used the same PDF for random behavior of wind speed [[
(
The applied symbols have been explained previously. The output power of WTG is as
(
Here, symbols have their usual meaning.
It is observed that solar radiation also follows the characteristic of Weibull PDF. Therefore, intensity of solar illumination is given by [[
(
(
The relevant symbols of (
The MCS-based method is the only method that is used as uncertainty evaluation. Further, PEM based method is also used to evaluate uncertainty of various parameters. It is an approximate method in which only a few points of the distribution function of a random variable are considered for calculations. These points are called concentrations. By using these points, the information of the output can be obtained [[
In this paper, (2 m + 1) based PEM method is applied because it provides very accurate results for probabilistic load flow (PLF) solutions by running MATLAB program of deterministic load flow only (2m + 1) times [[
Let x
(
Here,
For determination of standard central moments
Here, it is to be noted that m is number of variables to be considered in the PEM method.
Find the standard locations
(
As
Finally, output of any random variable is obtained as follows
(
(
Here the process expected output is E (z), and
In the present article, the problem is formulated for optimum placement of NDG units on the basis of unknown variables. These variables are as:
- The candidate buses are selected on the basis of bus voltage profile i.e., the bus having the lowest value of bus voltage given first priority for DG placement and so on. Thus the sequence of buses is identified for DG placement.
- Nearly 20% of total buses are selected for DG placement on bus voltage priority.
- Size of all DG units (in kW).
- The DG units which are to be placed are considered as WTG and SPV based DG as non-dispatchable types, DG and fueled operated DG (such as biomass and microturbines etc.) considered as dispatchable type DG.
- The placement of DG which is energized by renewable sources depends upon factors such as availability of wind speed and solar radiation.
- These resources are uncertain in nature and this uncertainty is modeled by using Weibull PDF.
- The proposed method is investigated under the scheme presented in Table 2.
- Constant power load is considered and its uncertainty is defined by the standard normal distribution function.
Table 2 describes the ranking of buses for DGs placement and the types of DGs considered for placement. Table 3 describes the technical details of wind turbine generator and solar operated DG in respect of wind speeds i.e., cut-in speed, nominal speed, and cut-out speed and solar radiations. Table 3 represents the data considered for different scenarios. In this table, the different scenarios have been created by considering the different values of shape and scale parameters of Weibull PDF. Moreover, Table 4 represents the three different values of scale and shape parameters of Weibull PDF and respective weight factors.
Optimum placement of DG unit of a given size is codified as a constrained nonlinear optimization problem. Further, the objective is to minimize the system total cost considered as an objective function which includes the investment cost, operating maintenance cost, network loss cost as per [[
(
where w
(
Backward forward sweep method (BFS) based load flow solutions are considered in present study. Following considerations have been assumed:
System under consideration is balanced, radial and fed by a substation.
Voltage of all nodes is close to 1 p.u.
Shunt admittances of lines are negligible.
The base case system power loss are evaluated on the basis of the following formula
(
Here symbols have their usual meanings.
Basically, two types of modeling constraints have been considered in the present study. These are given as below.
In this section, backward forward sweep (BFS) based method is considered for the solution of load flow equations as equality constraints in order to obtain bus voltage.
In the intelligent optimization technique, the inequality constraints should not be violated. Designer of network is intended to form direct relationship with other scientific terms to get feasible solution of power supply and voltage. The inequality constraints are specified as a minimum and maximum value of active and reactive powers produced by DG units. Therefore, inequality is designed in terms of system voltage and output power generated by these DG units.
(
(
(
(
(
where, V
The uncertainties have been addressed by three points PEM and MCS-based approaches. Further, for searching out the different candidate buses for optimal placement, different types of DG units are selected on the basis of increasing order of bus voltage profile to limit the search space. The technical details of wind turbine generator and solar operated DG in respect of wind cut-in speed, nominal speed and cut-out speed and solar radiations and different cost along with scenarios are considered taken from [[
The PEM-based method considered in the present paper is tested on the IEEE 33-bus radial distribution network. The single-line diagram of this system is depicted in Figure 2. The total load on IEEE 33-bus test radial distribution system is (3.715 + j2.300) MVA.
Proper encoding of solution vector is the fundamental tool for efficient application of any meta-heuristic optimization technique. The optimum solution of objective function basically depends upon sizing and siting of DG units. As a result, each solution is represented in the form of chromosomes. The illustration is shown as follows.
The solution vector of different types of DGs are to be allocated as, P = (P
The vector P is zero if there is no DG and its value is 1 or 2 or ......or N
GA is a random population-based search technique in which solutions are obtained in the form of chromosomes. It is an effective technique for both optimization and machine learning techniques. GA is easily parallelized and it is capable to work on discrete and continuous type problem. GA based technique is a time-consuming algorithm, which means its computational time is more compared to other nature-inspired techniques.
PSO is a nature-inspired evolutionary optimization technique. The technique works on the basis of social behavior like flocking of birds and schooling of fish. Basically, PSO starts with randomly generated solutions; these solutions are known as particles. Random velocities are assigned to these particles and are bound to keep track of their positions in search space in which they search the best solution called (pbest). Further, the particle swarm optimizer holds overall best solution, with reference to their position it is known as (gbest). The advantages of PSO are its computational time is less than heuristic techniques. It does not converge at local minima and disadvantages are it is not suitable for the combinatorial type of problems because of its continuous nature. Finally, the following two equations provide the velocity and position of each particle,
GA-PSO based hybridized technique mitigates the disadvantage of GA and PSO and yields better results in the formless computational time and local minima. The process of computation is presented in flowchart as shown in Figure 3.
Generally, the metaheuristic optimization technique starts with a random set of population of chromosomes known as initial solution. The initial population is created whose dimensions are product of (N
- Firstly, an integer is randomly selected between 1 and N
par for each possible solution. - From the vector dimension h with integer elements between 1 and N
par M is randomly selected, i.e., let's assume Npar = 8 and h = 5 is randomly selected, therefore M is filled with integer {1, 2,...,8} and h stand for Ist, IInd, IIIrd, IVth, and Vth gene respectively. - Lastly, a random number M is generated between 1 and N
C for every element of candidate scenarios placed in the gene of M randomly. It can be understood as, say M1 = 1 is first gene of the chromosome and selected randomly between 1 and NC , as M2 = 4 stands for fourth gene of the chromosome and picked up from 1 and NC so on.
Table 5 shows the various is results achieved by PEM based method for scenario-1, and Table 6 represents the simulation results obtained by MCS based method for scenario-1. Finally, Table 7 shows the values of investment cost and fitness function obtained by PEM and MCS based method, respectively, for different scenarios [[
Furthermore, it is observed from Table 5 and Table 6 that the execution time of hybrid optimization technique is more as compared to GA and PSO algorithms. Further, it is observed from Table 5 and Table 6 that the system loss reduction capability hybrid technique is more as compared to GA and PSO. Additionally, the execution time of PSO-based optimization is less than that of GA-based optimization technique.
Figure 4 demonstrates the convergence of fitness function and cost function, respectively, for scenario-1. From these figures, it is observed that the GA-PSO hybrid technique converges very fast while GA based technique takes more time to converge. It is evident from Figure 4 that GA-based technique converges after 46 iterations, the PSO-based technique converges after 15 iterations, and hybrid technique converges after 8 iterations. The same pattern is followed for cost function also as is evident from Figure 5. Furthermore, Figure 6 shows the voltage profile of the 33-bus system before and after DG placement under different techniques. From this figure, it is noticed that GA-PSO hybrid technique is the most powerful technique as far as voltage improvement is concerned among all methods i.e., GA and PSO. Moreover, PSO has better capability to improve the system voltage profile.
From Figure 7 it is observed that the convergence rate of PSO is slightly higher than hybrid-based technique and GA further gives same type results, which means its convergence rate is very slow and requires 49 iterations to finally settle down. From this figure, it is also noticed that PSO technique converges after 21 iterations and hybrid technique converges after 26 iterations. As we see from Figure 8 this figure also follows the same pattern as is shown in Figure 7.
Figure 9 represents the system bus voltage profile for scenario-2 by PEM-based technique. In this figure, it is observed that system voltage improves after DG placement. Further, it is noticed that GA-PSO hybrid technique is most capable of improving the system bus voltage followed by PSO and GA, respectively. Thus from this figure, it is concluded that the hybrid-based method is most suitable for enhancing the system voltage profile.
Figure 10 and Figure 11 are evidence of variation of fitness and cost functions by using GA, PSO, and GA-PSO based hybrid techniques respectively, under scenario-3. From these figures, it is observed that the GA based approach stabilizes and obtains optimized value of functions after 49 iterations and PSO provides optimized value after 45 iterations. Finally, GA-PSO based hybrid method provides an optimized value of objective after 34 iterations. Therefore, it is concluded that the GA-PSO-based optimization technique provides quick and better results among GA and PSO-based methods. Finally, Figure 12 shows the system voltage profile before and after DG placement by GA, PSO, and GA-PSO-based optimization technique in this method GA based technique provides better improvement in voltage profile than PSO and hybrid GA-PSO based technique.
The simulation results shown in Table 7, Table 8, Table 9 and Table 10 for PEM and MCS-based method for different optimization techniques for scenario-2 and scenario-3.
In Figure 6, Figure 9, and Figure 12 represent the system voltage profile (in p.u.) at different buses with and without DG using PEM-based method by incorporating three optimization technique such as GA, PSO, and hybrid (combination of GA and PSO). Further, it is observed from these figures that the hybrid-based optimization technique is more capable of enhancing the system voltage profile as compared to the other two techniques for scenario-1 and 2. Moreover, the same is applicable for MCS based method with and without DG units.
The proposed algorithm was applied on the IEEE 33-bus radial distribution system with the base power of 1 MVA and the base voltage of 12.66 kV. The total load of the distribution system is 3.715 MW and 2.3 MVAr. The GA-PEM, PSO-PEM, and hybrid (combination of GA and PSO) -PEM along with GA-MCS, PSO-MCS, and hybrid (combination of GA and PSO)-MCS-based method were applied to address uncertainty. The probabilistic load solutions were obtained on the basis of MCS and PEM-based techniques. To implement programs for PLF and GA, PSO and hybrid-based optimization techniques were developed in MATLAB environment. The load is considered as an uncertain parameter and varies according to the standard normal distribution function at each bus. Further, the uncertainties of wind speed and solar radiation are considered in the present study as following the properties of Weibull PDF. The candidate buses for DG placement are selected on the basis of voltage profiles of different buses and are put in ascending order. Further DG sizes and respective candidate buses are shown in Table 2. The fuel prices are also considered an uncertain parameter and follow the characteristic of geometric Brownian motion (GBM). The magnitude of voltage profile at each bus does not exceed ±5% of the nominal grid voltage. Further, it is assumed that the maximum limit of power in distribution lines is 4 MVA. The details of simulated results obtained in this method are presented in Table 5, Table 6, Table 7, Table 8, Table 9, Table 10 and Table 11, for PEM and MCS-based methods for different scenarios. From these results, it has been observed that the proposed method is competent enough to reduce the system loss by DG installation. Moreover, the voltage profile of the system with and without DG has been shown in Figure 6, Figure 9, and Figure 12. From these figures, it was observed that system voltage was enhanced significantly by placement of DG. Moreover, it is also observed that the hybrid-based optimization technique is more capable of enhancing the system voltage profile as compared to GA and PSO-based optimization techniques. Additionally, the convergence of fitness function is depicted in Figure 4, Figure 7, and Figure 10, for PEM-based method embedded with GA, PSO, and hybrid-based optimization techniques for three scenarios. The cost functions are shown in Figure 5, Figure 8, and Figure 11 for the PEM-based method embedded with GA, PSO and hybrid-based optimization techniques for scenario-3. Different tables show the simulation results are obtained by the above-mentioned method. It is concluded that the loss reduction capacity of different methods for different scenarios is almost the same for GA, PSO, and hybrid-based optimization techniques.
In this paper (2m + 1) PEM method is presented, which addresses the uncertainty of load, power generated by solar and wind operated DG. The probabilistic load flow (PLF) solution is carried out to investigate the uncertainty of these parameters. The PEM and MCS-based methods were applied to achieve solutions of PLF. Further, for optimization of fitness function and minimization of investment cost with DG placement, GA, PSO, and Hybrid (combination of GA and PSO) based methods were applied with PEM and MCS. The simulation results obtained by the PEM method were compared with the MCS method, which is used as a benchmark method for uncertainty evaluation. Moreover, as it is evident from the results obtained by PEM-GA, PEM-PSO and PEM-Hybrid based methods are almost similar to MCS-GA, MCS-PSO, and PEM-Hybrid method for three different scenarios. On the basis of obtained results, it is noticed that the simulation time of MCS based method is longer, as compared to PEM based method. Further, it is observed that for the same scenario, the hybrid-based optimization technique provides a better improvement in system voltage profile. Finally, it is concluded that PEM based method provides similar results as MCS with a lesser computational burden. For scenario-3, the system voltage profile is enhanced significantly by GA.
Graph: Figure 1 Flow chart for probabilistic load flow considering point estimate method.
DIAGRAM: Figure 2 Single line diagram for IEEE 33-bus radial distribution network.
Graph: Figure 3 Flow chart for a hybrid GA-PSO based approach.
Graph: Figure 4 Variation of fitness function for GA, PSO, and hybrid-based optimization technique using PEM based technique for scenario-1.
Graph: Figure 5 The cost function for GA, PSO, and hybrid-based optimization technique using PEM based method for scenario-1.
Graph: Figure 6 System voltage profile for scenario-1 with and without DG by PEM based method in case of GA, PSO and hybrid GA-PSO based optimization technique.
Graph: Figure 7 Fitness function for GA, PSO and hybrid optimization technique using PEM based method for scenario-2.
Graph: Figure 8 The cost function for GA, PSO, and hybrid optimization technique using MCS based method for scenario-2.
Graph: Figure 9 System voltage profile using PEM based load flow method and GA, PSO, and hybrid based optimization technique for scenario-2 with and without DG unit.
Graph: Figure 10 Fitness value for GA, PSO, and hybrid-based optimization technique using PEM for scenario-3.
Graph: Figure 11 Cost function for GA, PSO, and hybrid-based optimization technique using PEM method for scenario-3.
Graph: Figure 12 System voltage profile using PEM scenario-3 load flow method and GA, PSO and hybrid-based optimization technique with and without DG unit.
Table 1 Brief detail of extensive literature review.
Reference Number and Year of Publication Type of Application Main Contribution and Characteristics Analyzed Method/Remarks Recommended System/Similarity with Current Approach [ Optimal placement of DGs (Solar, wind and fueled) based DG in distribution system Design objective function by combining different costs associated to DGs operation and maintenance and designed objective function is optimized by GA based heuristic technique Application of 3-point PEM based to address uncertainty of different parameters and results are compared with MCS based technique The current approach is applied on standard IEEE 33-bus test distribution system and designed objective function is optimized by application of GA, PSO and GA-PSO based hybrid technique [ Extends the application PEM Main contribution of paper to address the random parameters of DG by PEM based method In proposed paper, the PEM -based method is highlighted to investigate the uncertainty of unsymmetrical system in presence of DG. Only probabilistic load flow concept is similar [ Hong's PEM based method for PLF To address the uncertainty associated with PLF, Hong's PEM-based methods are applied. The authors considered two scenarios as IEEE 14-bus and IEEE 118-bus distribution network. Similarity only in this way that PLF by PEM based method [ New algorithm for PLF Uncertainties of bus injection and line parameters is investigated by PEM based method The proposed algorithm is implemented on IEEE test system and deterministic load flow is used to address the PLF PLF is similar parameter to the current approach [ Applied 2-PEM based method for investigating the uncertain parameters 2-PEM based method for investigation of uncertainties of optimal power flow (OPF) 2-PEM and MCS based methods are applied to address the optimal power flow (OPF) and proposed method is implemented on 129-bus test system No similarity except PEM and MCS based technique [ Evolutionary or swarm-based search algorithms Three nature inspired techniques such as IGA, IPSO and ICSO are applied. The cost function which addressing the benefit of network by optimal sizing and siting of DGs and SCs is optimized by IGA, IPSO and ICSO optimization approaches No similarity: IGA, IPSO and ICSO are applied in different perspective [ Main focus on optimal power flow (OPF) by considering wind uncertainty The objective function comprising economic aspect of system is optimized by best guided artificial bee colony optimization algorithm (GABC) The uncertainty of wind speed is addressed by Weibull density function for optimization GABC technique is applied and the proposed technique is tested on IEEE 30-bus test system Applied methods are different, optimized function is entirely different test system is different but uncertainty of wind speed is addressed by same Weibull density function [ The chance constrained programming (CCP) based mathematical model proposed Authors developed a CCP based mathematical model comprising various costs including reliability cost which is optimized by considering different operating constrains of DGs. An objective function which includes all the attributes of optimal DG sizing and placement are included is optimized by Monte-Carlo simulation embedded GA is applied in order to optimize the considered objective function and the proposed algorithm is tested on IEEE 37-node test feeder Considered objective function is different and is optimized by GA based technique, therefore GA is one of the similar technique which is applied but objective function is entirely different [ Optimum placement and sizing of DG in distribution network by considering a novel method. Authors formulated a multi-objective functions function which comprising of different costs and optimized by an improved NSGA-II based technique. The proposed algorithm is tested on IEEE 69-bus test radial distribution system the objective function is optimized by NSGA-II for optimal placement and sizing of DG No similarity [ The paper proposed a meta-heuristic voltage stability index based method for optimal siting and sizing of DG A multi-objective based objective function comprising functions which leads towards maximize the benefits by multi-objective particle swarm optimization (MOPSO) based method The proposed method is tested on 30-bus test system and 41-bus Indian distribution system in order to improve the voltage profile and minimize cost No similarity [ A methodology has been proposed for optimally allocating wind-based DG The authors formulated a multi-objective function which combines all possible operating conditions of load and wind operated DGs with their probabilities. The objective function is optimized by MINLP based technique for optimal allocation of wind based DG in order to minimize annual energy losses. Further, in present reference uncertainty of load is not considered Not similarity [ Analysis of probabilistic based load flow by 2-PEM based technique Comparison of obtained results in terms of uncertainty of wind based DG by proposed method and MCS based method 2-point PEM based method is analyzed for uncertainty with MCS based and tested on IEEE 9-bus and 57-bus test systems Only probabilistic load flow concept is similar. In proposed method 3- PEM based load flow is used but in [ [ Analysis of probabilistic power flow PEM based method for analysis of uncertainty of wind To address the uncertainty of wind the PEM and MCS based methods are compared Only probabilistic load flow concept is similar [ Probabilistic approach for uncertain parameters The proposed paper addresses the time dependent uncertainty of load by considering PEM based technique PEM based technique is considered to address the uncertainty of load of different nature and an objective function is constituted by active power, reactive power and voltage profile indices The objective function is optimized by metaheuristic optimization technique Invasive Weed Optimization (IWO) based technique for size and location of DGs [ Evolutionary programming (EP) based method Address the uncertainty of wind and solar based DGs by considering probabilistic based method Uncertainty of wind and solar based DGs by probabilistic method and model the uncertainties of wind speed, solar radiation and load the method is implemented on IEEE 33-bus and 69-bus test distribution system No similarity [ Conservation voltage reduction (CVR) and distributed-generation (DG) integration This paper investigates the interactions between CVR and DG placement to minimize load consumption in distribution networks, while keeping the lowest voltage level within the predefined range The optimal placement of DG units is formulated as a stochastic optimization problem considering the uncertainty of DG outputs and load consumptions for this purpose A sample average approximation algorithm-based technique is developed and proposed technique is tested on IEEE 37-bus test system No similarity [ Novel approach for PLF based PEM PEM and MCS based techniques are applied to investigate different uncertain parameters The proposed algorithm is implemented on IEEE 14-bus test system In this paper only PLF is investigated by considering uncertainty of load, while current approach is applied on IEEE 33-bus test system to address uncertainty parameters like Solar, wind, and load and objective function is constituted by considering different weight [ First paper that address the concept of PLF In this paper concept of PLF was introduced to address the uncertainty of node data PLF solution is obtained by application of deterministic load flow PLF is addressed by author to investigate the uncertainty of node data but in current research paper the uncertainty of many parameters is addressed [ The paper deals Energy storage system (ESS) in active distribution network The objective of the problem is to minimize the yearly cost ESSs. The authors investigated the AC power flow based method to address the active and reactive power using PSO based optimization method. No similarity: only metaheuristic based PSO optimization is applied to optimized for different objective function [ Probabilistic optimal power flow (POPF) by considering wind speed To obtain the solution of POPF model for integration of wind operated DG in distribution system two kinds of 2m, 2 The uncertainty of load is investigated by two kinds of 2m, 2 Only (2m+1) based PEM is provide better solution it is confirmed by this research paper. In our research article the uncertainty of DGs power, load and fuel costs are investigated by the same method along with objective function is optimized by hybrid GA-PSO, PSO and GA based optimization methods. [ Rosenbleuth's based PEM technique is analyzed Author investigated the Rosenbleuth's based PEM technique to address the uncertainty of different random parameters Investigate the uncertainty of different random parameters by MCS and Rosenbleuth's based PEM technique and observed that PEM based technique reduces the computational burden 2m point based PEM is employed but in our paper (2 [ Weibull based PDF is applied to investigate the uncertainty To assess the wind energy potential, Weibull based probability density functions is employed to investigate the uncertainty of wind speed In this research article an investigation of uncertainty of wind operated DG Weibull based PDF is employed on real time data of 9 power stations which are located in the United Arab Emirates (UAE). No similarity [ Monte Carlo Simulation (MCS) based probabilistic load flow is investigated Uncertain nature of wind operated DG is investigated by MCS based probabilistic load flow is addressed to optimal placement of DG and capacitor to improve system voltage profile The optimal sizing and siting of wind power DG (WPDG) and capacitor in distribution system is optimized by modified PSO (MPSO) based optimization technique is employed and proposed technique is tested on IEEE 33-bus distribution system No similarity [ AC probabilistic optimal power flow (P-OPF) is applied to address the uncertainty using MCS. Applied PEM based method and results are compared with MCS based method for wind speed and line outage factor uncertainty PEM and MCS based technique to address the wind speed and line outage factor uncertainty and proposed technique is applied on IEEE 30-bus test system No similarity except PEM and MCS based technique [ In proposed technique integrates the diagonal band Copula and sequential Monte Carlo The authors Considered the multivariable -based stochastic method to investigate uncertainty of solar based systems. Objective function is optimized by Big Bang-Big crunch method for optimal placement of DGs and proposed method is applied on IEEE 37-bus test system MCS based method is only similar but applied in different objective function [ Optimal allocation of biomass based DGs. Authors suggested an adaptive equilibrium optimizer technique (EO) to enhance performance of biomass DGs. EO techniques is applied to IEEE-33-bus and practical large-scale 141-bus system of AES-Venuzuala in metropolitan area of Caracas. No similarity [ MOPSO-based optimization for optimum sizing and allocation of shunt capacitor banks (SCBs) and DGs. Multi-objective function comprising the balancing current in different sections along with bus voltage stability and system power loss is optimized by SPEA, NSGA, MODE and combination of ICA/GA. The uncertainty of loads is modeled by using fuzzy data theory, and objective function is optimized by SPEA, NSGA, MODE and combination of ICA/GA, the proposed algorithm is tested on IEEE 33-bus and an actual realistic 94 bus Portuguese test systems No similarity: different techniques are applied for uncertainty and optimization of objective function [ An Improved Sunflower Optimization Algorithm (ISFOA)-based Monte Carlo simulations Authors Investigated performance enhancement of smart distribution system using ISFOA technique. The technique is tested on IEEE-33 and real time 84-bus distribution system No similarity [ Optimal distributed generation units are correlated with fault current limiter sites. Different types of DGs are allocated and correlated in a single stage with fault current limiters (FCLs). A fuzzy-based multiobjective (FBMO) formulation is provided for performance improvement both in normal and faulty operating conditions. The proposed method is studied for 33-bus, 69-bus, and the Egyptian East Delta distribution systems. No similarity [ A modified marine predators optimizer (MMPO) is used for simultaneous distribution network reconfiguration (DNR). Authors have shown the superiority of the proposed MMPO for simultaneous DNR and DG allocation. This technique is tested on IEEE 33-bus and 69-bus distribution systems. No Similarity [ A methodology developed on the basis of probabilistic generation-load model The authors applied a combination of all working conditions of renewable DGs-based methodology to assess the probabilities by MINLP based technique. The paper addresses the optimal sizing and siting of all types of DGs (dispatchable and non-dispatchable) with the objective to minimize total annual energy loss of the system No similarity with proposed method. [ Analytical and MCS based technique Investigation of system reliability when wind based DG is integrated in distribution system MCS and analytical based methods are employed in order to calculate the system reliability in the form of LOEE and LOLE in presence of wind operated DG No similarity
Table 2 Candidate schemes for a different type of DGs placement and their sizing.
Ranking of Buses for DG Placement Installed Capacity of DG (in kW) DG Type * 18 20 40 60 80 100 1, 2 17 40 80 120 160 200 1, 2, 3 16 40 80 120 160 200 1, 2 33 20 40 60 80 100 1, 2 32 20 40 60 80 100 1, 2, 3 15 100 200 300 400 500 1, 2, 3 31 100 200 300 400 500 1, 2, 3 14 40 80 120 160 200 1, 2 13 40 80 120 160 200 1, 2, 3
Table 3 Technical specifications of DGs for different type of DGs.
S.No. DG Type Technical Specifications 1. Wind turbines 2. 3. 4. Power factor = 0.9 lagging 5. Solar photovoltaic 6. Power factor = 1 7. Fueled DGs Stable power 8. Power factor = 0.9 lagging
Table 4 Data for wind and solar-based DG for different scenarios.
S.No. Scenarios Different Parameters ( Different Parameters ( Weight for Different Objective Functions 1. 1 b1 = 0.1, b2 = 0.11, b3 = 0.34, b4 = 0.34, b5 = 0.11 2. 2 b1 = 0.1, b2 = 0.11, b3 = 0.34, b4 = 0.34, b5 = 0.11 3. 3 b1 = 0.34, b2 = 0.11, b3 = 0.34, b4 = 0.11, b5 = 0.10
Table 5 Simulation results for optimal DGs size and their corresponding optimum location by PEM based method for scenario-1 using different metaheuristic optimization techniques.
Particle Swarm Optimization (PSO) Genetic Algorithm (GA) Hybrid Optimization (GA and PSO) Execution time (in seconds) Power loss without DG (in kW) Power loss with DG (in kW) Execution time (in seconds) Power loss without DG (in kW) Power loss with DG (in kW) Execution time (in seconds) Power loss without DG (in kW) Power loss with DG (in kW) 156.824654 202.6771 93.3447 238.258334 202.6771 89.0562 372.409068 202.6771 98.6957 DG TYPE DG TYPE DG TYPE Wind DG size (in kW) and location (bus) Solar DG size (in kW) and location (bus) Fueled DG size (in kW) and location (bus) Wind DG size (in kW) and location (bus) Solar DG size (in kW) and location (bus) Fueled DG size (in kW) and location (bus) Wind DG size (in kW) and location (bus) Solar DG size (in kW) and location (bus) Fueled DG size (in kW) and location (bus) 0 (0) 40 (13) 120 (13) 500 (15) 100 (15) 300 (15) 300 (15) 60 (32) 100 (15)
Table 6 Simulation results for optimal DGs size and their corresponding optimum location by MCS for scenario-1 using different metaheuristic optimization techniques.
Particle Swarm Optimization (PSO) Genetic Algorithm (GA) Hybrid optimization (GA and PSO) Execution time (in seconds) Power loss without DG (in kW) Power loss with DG (in kW) Execution time (in seconds) Power loss without DG (in kW) Power loss with DG (in kW) Execution time (in seconds) Power loss without DG (in kW) Power loss with DG (in kW) 1054.205281 202.6771 89.9687 1615.485610 202.6771 91.4555 2656.724803 202.6771 95.2132 DG TYPE DG TYPE DG TYPE Wind DG size (in kW) and location (bus) Solar DG size (in kW) and location (bus) Fueled DG size (in kW) and location (bus) Wind DG size (in kW) and location (bus) Solar DG size (in kW) and location (bus) Fueled DG size (in kW) and location (bus) Wind DG size (in kW) and location (bus) Solar DG size (in kW) and location (bus) Fueled DG size (in kW) and location (bus) 400 (15) 40 (14) 400 (15) 500 (15) 100 (15) 300 (15) 400 (31) 400 (14) 400 (15)
Table 7 Comparison of different parameters using PEM and MCS based method for uncertainty for different scenarios of various optimization techniques.
Probabilistic Technique Different Intelligent Technique Investment Cost before DG (in USD) Investment Cost after DG (in USD) Fitness Value of Function PEM scenario 1 GA 1.4425 × 106 7.3094 × 105 406.0774 PSO 7.6799 × 105 426.6635 HYBRID 7.7463 × 105 430.3498 PEM scenario 2 GA 7.4988 × 105 416.5972 PSO 7.4988 × 105 429.2726 HYBRID 7.8558 × 105 436.4327 PEM scenario 3 GA 8.6499 × 105 480.5495 PSO 9.3932 × 105 521.8421 HYBRID 1.0225 × 106 568.0348 MCS scenario 1 GA 7.3850 × 105 410.2800 PSO 7.8040 × 105 433.5537 HYBRID 7.7596 × 105 431.0876 MCS scenario 2 GA 7.4615 × 105 414.5283 PSO 7.7363 × 105 429.7956 HYBRID 7.8163 × 105 434.2403 MCS scenario 3 GA 8.4815 × 105 471.1919 PSO 1.0075 × 106 559.7419 HYBRID 9.8414 × 105 546.7418
Table 8 Results for optimal DGs size and their corresponding optimum location using by PEM for scenario-2 using different metaheuristic optimization techniques.
Particle Swarm Optimization (PSO) Genetic Algorithm (GA) Hybrid optimization (GA and PSO) Execution time for algorithm (in seconds) Power loss without DG (in kW) Power loss with DG (in kW) Execution time for algorithm (in seconds) Power loss without DG (in kW) Power loss with DG (in kW) Execution time for algorithm (in seconds) Power loss without DG (in kW) Power loss with DG (in kW) 154. 045316 202.6771 92.7158 235.338242 202.6771 93.6265 396.141804 202.6771 93.9004 DG TYPE DG TYPE DG TYPE Wind DG size (in kW) and location (bus) Solar DG size (in kW) and location (bus) Fueled DG size (in kW) and location (bus) Wind DG size (in kW) and location (bus) Solar DG size (in kW) and location (bus) Fueled DG size (in kW) and location (bus) Wind DG size (in kW) and location (bus) Solar DG size (in kW) and location (bus) Fueled DG size (in kW) and location (bus) 200 (15) 300 (15) 120 (13) 0 (0) 500 (15) 100 (31) 200 (15) 80 (14) 400 (15)
Table 9 Results for optimal DG sizes and their corresponding optimum location using by PEM for scenario-3 using different metaheuristic optimization techniques.
Particle Swarm Optimization (PSO) Genetic Algorithm (GA) Hybrid optimization (GA and PSO) Execution time for algorithm (in seconds) Power loss without DG (in kW) Power loss with DG (in kW) Execution time for algorithm (in seconds) Power loss without DG (in kW) Power loss with DG (in kW) Execution time for algorithm (in seconds) Power loss without DG (in kW) Power loss with DG (in kW) 168.890842 202.6771 94.5629 248.739512 202.6771 93.7272 384.948469 202.6771 93.4615 DG TYPE DG TYPE DG TYPE Wind DG size (in kW) and location (bus) Solar DG size (in kW) and location (bus) Fueled DG size (in kW) and location (bus) Wind DG size (in kW) and location (bus) Solar DG size (in kW) and location (bus) Fueled DG size (in kW) and location (bus) Wind DG size (in kW) and location (bus) Solar DG size (in kW) and location (bus) Fueled DG size (in kW) and location (bus) 160 (16) 300 (15) 500 (15) 0 (0) 500 (15) 300 (15) 300 (15) 120 (14) 200 (15)
Table 10 Results for optimal DGs size and their corresponding optimum location using by MCS for scenario-2 using different metaheuristic optimization techniques.
Particle Swarm Optimization (PSO) Genetic Algorithm (GA) Hybrid optimization (GA and PSO) Execution time (in seconds) Power loss without DG (in kW) Power loss with DG (in kW) Execution time (in seconds) Power loss without DG (in kW) Power loss with DG (in kW) Execution time (in seconds) Power loss without DG (in kW) Power loss with DG (in kW) 1047. 508148 202.6771 92.7158 1647.367694 202.6771 93.7272 2733.378447 202.6771 94.7006 DG TYPE DG TYPE DG TYPE Wind DG size (in kW) and location (bus) Solar DG size (in kW) and location (bus) Fueled DG size (in kW) and location (bus) Wind DG size (in kW) and location (bus) Solar DG size (in kW) and location (bus) Fueled DG size (in kW) and location (bus) Wind DG size (in kW) and location (bus) Solar DG size (in kW) and location (bus) Fueled DG size (in kW) and location (bus) 300 (15) 0 (0) 500 (15) 0 (0) 40 (32) 100 (31) 200 (14) 80 (13) 300 (15)
Table 11 Results for optimal DG sizes and their corresponding optimum location using MCS for scenario-3 using different metaheuristic optimization techniques.
Particle Swarm Optimization (PSO) Genetic Algorithm (GA) Hybrid Optimization (GA and PSO) Execution time for algorithm (in seconds) Power loss without DG (in kW) Power loss with DG (in kW) Execution time for algorithm (in seconds) Power loss without DG (in kW) Power loss with DG (in kW) Execution time for algorithm (in seconds) Power loss without DG (in kW) Power loss with DG (in kW) 1063.508148 202.6771 103.7904 1716.338242 202.6771 93.6265 2773.036281 202.6771 100.2172 DG TYPE DG TYPE DG TYPE Wind DG size (in kW) and location (bus) Solar DG size (in kW) and location (bus) Fueled DG size (in kW) and location (bus) Wind DG size (in kW) and location (bus) Solar DG size (in kW) and location (bus) Fueled DG size (in kW) and location (bus) Wind DG size (in kW) and location (bus) Solar DG size (in kW) and location (bus) Fueled DG size (in kW) and location (bus) 200 (14) 120 (16) 400 (15) 0 (0) 100 (31) 100 (15) 80 (16) 160 (14) 200 (15)
Conceptualization, P.P., D.C.M., H.M., M.A.A.; methodology, P.P., D.C.M., H.M.; software, M.A.A., I.A.K.; validation, P.P., D.C.M., H.M., M.A.A., I.A.K.; formal analysis, P.P., D.C.M., H.M., M.A.A., I.A.K.; investigation, P.P., D.C.M., H.M., M.A.A., I.A.K.; resources, D.C.M., H.M.; data curation, M.A.A., I.A.K.; writing—original draft preparation, P.P., D.C.M.; writing—review and editing, M.A.A., I.A.K.; visualization, P.P., D.C.M., H.M., M.A.A., I.A.K.; supervision, D.C.M., I.A.K.; project administration, H.M., M.A.A.; funding acquisition, H.M., M.A.A. All authors have read and agreed to the published version of the manuscript.
The authors extend their appreciation to the Researchers Supporting Project at King Saud University, Riyadh, Saudi Arabia, for funding this research work through the project number RSP-2021/278.
Not applicable.
Not applicable.
Not applicable.
The authors declare no conflict of interest.
The authors would like to acknowledge the technical and non-technical support from Intelligent Prognostic Private Limited Delhi, India Researcher's Supporting Project, Intelligent Prognostic Private Limited Delhi, India.
By Prem Prakash; Duli Chand Meena; Hasmat Malik; Majed A. Alotaibi and Irfan Ahmad Khan
Reported by Author; Author; Author; Author; Author