Distributed generations (DGs) consider an important issue for the electric utilities as DGs relieve capacity constraints on the generation, transmission and distribution systems and obviate the need to build new facilities . The increasing penetration of DGs added to the distribution power system create new challenges, which can be classified into three main categories as technical, commercial and regulatory . On the other hand, several benefits accrue by integrating DGs with electric utility, which can be classified into technical and economical. These benefits should be analyzed and quantified in order to increase the potential and value of DGs penetration. The major technical benefits are standby capacity, reliability, security and power quality, environmental and health concern. While, the economic benefits are reduced as costs associated with power losses, deferred investments for generation, transmission or distribution upgrades, lower operating costs due to peak shaving and reduced land use for generation   , in recent years, the integration of the DGs into distribution networks has reached high penetration levels. The proper location and size of the DGs enhance their benefits in the distribution networks. DGs can be powered from either conventional or renewable fuel. In Egypt, wind energy (WE) is considered one of the most promising renewable energy (RE) source as it is renewable, plentiful, clean, and uses little land . Egypt has an excellent wind regime with wind speeds of about 10 m/s at Gulf of Suez coast. While, the large desert regions both to the east and the west of the Nile River, as well as parts of Sinai, have average annual wind speeds of 7 - 8 m/s . The disadvantage of WE is its seasonal variations. If wind power is to supply a significant portion of the demand, either backup power or electrical energy storage (EES) system is needed to ensure that loads will be supplied in reliable way  . Recently greater focus has been placed on the use of new EES technologies such as fuel cell (FC), which depend on Hydrogen (H2) as a fuel   . The electric distribution system (EDS) of Mersi-Mtrouh City as a part of the Egyptian distribution network can be supplied completely from wind farm combined with proton exchange membrane fuel cell (PEMFC) system, in an attempt to isolate it from the grid. The surplus power
generated from the wind farm can be used to produce H2 from the electrolysis of water in an electrolyzer, H2 can be stored in H2 tanks and converted to electrical energy through the PEMFC, for power generation at deficit times . In the literature, the hybrid power systems are discussed, for example, a hybrid photovoltaic fuel-cell generation system employing an electrolyzer for H2 generation and battery for storage purpose is designed and simulated in , power estimation model wind-hydrogen-fuel-cell-inverter system is presented in , photovoltaic/hydrogen/battery hybrid system is simulated in , the simulation results of a grid dependent hybrid system integrating a FC, photovoltaic(PV) and wind for distributed energy applications are shown in , in  a hybrid system consists of solar photovoltaic system and fuel cell system is used to supply the required electricity to a single phase load with the help of Pulse Width Modulation Voltage Source Inverter, a stand-alone hybrid microgrid consisting of wind turbines, photovoltaic arrays and storage battery banks is developed for use in an area in China , a method for optimal sizing of an off-grid hybrid micro-grid system in order to achieve a certain load demand is discussed in . However, various algorithms have been used to solve the DGs placement problem such as fuzzy-genetic algorithm , genetic algorithm , firefly algorithm , bat-inspired algorithm, particle swarm optimization technique , ant colony , and simulated annealing   . In this paper, a hybrid wind energy and proton exchange membrane fuel cell system (WE/PEMFC) is designed for supplying a part of the Egyptian distribution network in a reliable way. This system is dependent on WE as a RE source and PEMFC as the EES. The system load demand is predicted up to year 2022 for Mersi-Matrouh City, and the design of the hybrid WE/PEMFC system is applied. The optimal allocation of inserting this hybrid system in the EDS of Meri-Matrouh City is presented by a two-stage procedure. In the first stage, the sensitivity analysis using three sensitivity factors (SFs) is employed to select the candidate buses for locating the hybrid system. In the second stage, the binary Crow search algorithm (BCSA), discrete Jaya algorithm (DJA) and binary particle Swarm optimization (BPSO) techniques are investigated to find the optimal allocation of the hybrid system. Multi system constraints will be considered to achieve the suggested single and multi-objective function (SOF/MOF). This paper consists of five subsections; Section 1 is the introduction, in Section 2 the proposed procedure for optimal allocation of the hybrid system is discussed, the SFs is illustrated, SOF/MOF and system constraints are discussed, the backward/forward sweep (BFS) method used in load flow calculation is illustrated and a brief description of the three optimization techniques used. Section 3 discusses the design and economics of the wind farm and the PEMFC system. Section 4 is the application and Section 5 is conclusion.
2. Proposed Optimal Allocation of Hybrid WE/PEMFC System Using Parallel Sensitivity Factors
This paper proposes an efficient optimization procedure based on three SFs parallel with three different optimization techniques. First, the SFs will be used to select the candidate system buses for the installation of the hybrid WE/PEMFC system to reduce the search space of the optimization techniques. Second, the three different optimization techniques BCSA, DJA and BPSO are applied to find the optimal allocation of the hybrid WE/PEMFC system by achieving SOF/MOF. Third, the obtained results are compared with each other. The objective functions aim to minimize total active and reactive power losses and to improve voltage profile.
2.1. Sensitivity Analysis
Voltage sensitivity factor (VSF), active losses sensitivity factor (ALSF) and reactive active losses sensitivity factor (RLSF) are presented in this paper to rank the system buses according to their severity, to detect the candidate system buses that will need to install the hybrid WE/PEMFC system.
● VSF is represented by the ratio of the magnitude of the base case voltage to the minimum limit of voltage (0.95 pu) as follow :
The active power losses (APLs) of the distribution line connecting between buses p and q can be formulated as follow:
where; P is the total effective active power supplied beyond the bus q, Q is the total effective reactive power supplied beyond the bus q, Rpq is the resistance of the distribution line connecting between buses p and q and Vq is the voltage at bus q.
● ALSF can be written as:
The reactive power losses (RPLs) of the distribution line connecting between buses p and q can be formulated as follows:
where; Xpq is the reactance of the distribution line connecting between buses p and q.
● RLSF can be written as:
The system buses that have the smallest VSF positive values are considered as the candidate buses that need to install the hybrid system, and will be placed at the top of the VSF list. While, the buses that have the largest ALSF and RLSF positive values are considered as the candidate buses that need to install the hybrid system, where these buses will be placed at the top of the ALSF list and RLSF list, respectively. The candidate system buses for the installation of the hybrid system based on the VSF list, ALSF list and RLSF list are inserted as control variables in the optimization algorithms.
2.2. Single and Multi-Objective Functions
The proposed optimization techniques aim to determine the best allocation of the hybrid WE/PEMFC system in the EDS of Mersi-Matrouh City, by minimizing bus voltage index (BVI), active losses index (ALI) and reactive losses index (RLI) as SOF and MOF.
2.2.1. Single Objective Function (SOF)
By taking the ratio of a measure of an attribute with and without the hybrid WE/PEMFC system at the same loading conditions, some indices can be derived for any attribute as follows:
1) Minimizing BVI:
where; VB and Vb,H are the system bus voltage without and with the hybrid WE/ PEMFC system, respectively.
2) Minimizing ALI:
where, ALsk,H is the APLs of section k of the EDS with considering the hybrid WE/PEMFC system. TALsys is the total active power losses (TAPLs) of the EDS without considering the hybrid WE/PEMFC system.
3) Minimizing RLI:
where, RLsk,H is the RPLs of section k of the EDS with considering the hybrid WE/PEMFC system. TRLsys is the total reactive power losses (TRPLs) of the EDS without considering the hybrid WE/PEMFC system.
2.2.2. Multi-Objective Function (MOF)
The MOF is formulated as:
where; W1,W2 and W3 are the weighting factors (WFs) for BVI%, ALI% and RLI%, respectively. And can be obtained by trial and error.
2.2.3. System constraints
The SOF/MOF are subjected to the following system constraints:
1) System bus voltage constraints
The voltage at each bus VB must be within their minimum and maximum limits (Vbmin and Vbmax) as:
2) Total number of hybrid units (DGs) constraint
The number of DGs (NDGs) must be less than or equal to the maximum permissible number of DGs (NDGmax).
3) Active power flow constraint
The active power flow in the kth line (AFK) should be less than or equal to the maximum permissible active power flow (AFKmax).
4) Reactive power flow constraint
The reactive power flow in line k (RFK) should be less than or equal to the maximum permissible reactive power flow (RFKmax).
2.3. Backward-Forward Sweep (BFS) Load Flow Method
The BFS method is one of the most common ways used for load flow calculation of distribution system as it is fast, simple, robust convergence and low memory requirement for processing with efficiencies and solution accuracies computational . The BFS algorithm involves three iterative basic steps. The three steps are named as the nodal current calculation, the backward sweep and the forward sweep. These steps based on Kirchoff’s current law (KCL) and Kirchoff’s voltage law (KVL).
2.4. Optimization Techniques
In this paper, three optimization techniques are used to determine the best allocation of the hybrid WE/PEMFC system in the EDS of Mersi-Matrouh City. These optimization techniques are the BCSA, DJA and BPSO.
2.4.1. Binary Crow Search Algorithm (BCSA)
CSA is a nature-inspired algorithm presented by Askarzadeh in 2016 , CSA mimics the behavior of Crow birds and their social interaction. Crows watch other birds and observe where they hide their food. Once the other birds leave, they steal their food. The Crows use their intelligence to hide their excess food in a hideout spot and restore this food when they need. If a Crow discovers another one is going after it, it will try to deceive that Crow and will go to another position. CSA can be implemented as follow  :
● Initialization: the Crow in the folk represents a solution to the problem.
● If there is a solution space with dimension d has a Crow folk of N Crows, the position X of Crowi at iteration t can be expressed by the vector:
● Each Crow has a memory to store the best position of its stored food source.
● The memoryM of all Crows at iteration t for dimension d are initialized as follows:
● Solution evaluation (fitness): each solution in the population is evaluated using objective function.
● Position update: to update the position of the Crow between Crows J and i, there are two cases:
Case 1: Crow J does not recognize that Crow i is going after it. Hence, Crow i will get close the storing place of Crow J.
Case 2: Crow J recognizes that Crow i is going after it, so, it will move to another position to deceive Crow i and to save its food.
where; AP is the awareness factor and Fli,t is the flight length of Crow i at iteration t. The main reasons of using CSA are for easy implementation, few control parameters to adjust, fast convergence speed and high efficiency. To further enhance the performance of the classical CSA, BCSA is proposed. At each position, the solution considers values “0” or “1”, where, “0” represents the solution that is not selected and “1” represents the solution that is selected.
2.4.2. Discrete Jaya Algorithm (DJA)
Jaya algorithm is one of the newly developed population based optimization method, which was proposed by Rao in 2016 . The Jaya algorithm is simple in putting into use as a solution is updated only in a single phase by employing a single equation. The basic concept of the Jaya is based on the best and worst solution of objective function. Jaya needs only the common control parameters and does not need any algorithm clearly stated control parameters . A new DJA is the optimization algorithm used here for finding the optimal allocation of the hybrid system problem. In this algorithm, the solution contains two phases, meanwhile the candidate locations are rounded about its real values and the number of candidate buses are determined by binary strategy. The number of the best locations, with their hybrid WE/PEMFC system, are updated by the DJA in each iterations. The algorithm is continued until the termination criterion is met.
2.4.3. Binary Particle Swarm Optimization (BPSO) Technique
The mechanics of PSO technique are inspired by the swarming of biological populations. In PSO technique, the particles change their positions with time . The classical PSO is operated on real values, here the BPSO is proposed which is operated on [0, 1] values.
2.4.4. Proposed Allocation Procedure of Hybrid WE/PEMFC System
The optimal allocation procedure is illustrated below:
1) Input data: input the number of system buses, load demand at each bus, transmission lines’ impedance, control parameters of the used optimization technique.
2) Perform the initial load flow analysis using the BFS method to calculate the SFs, the VSF are arranged in ascending order for all buses of the study system. While, the ALSF and RLSF are arranged in descending order for all system lines. Thus, the candidate buses for inserting the hybrid WE/PEMFC system can be determined.
3) Implement one of the three proposed optimization techniques.
4) Perform the load flow analysis using BFS method again to obtain the voltage at each system bus, APLs and RPLs.
5) Calculate the fitness function of SOF/MOF for each proposed solution.
6) Repeat until the maximum number of iterations is reached or the solution is obtained.
7) The previous steps will be repeated with the other two proposed optimization techniques. Then the obtained results are compared with each other.
3. Optimal hybrid WE/PEMFC System Approach
A hybrid WE/PEMFC system is designed, in which the surplus power from the wind farm is used to generate H2 through the water electrolyzer, the H2 will be stored in the H2 tank, then, H2 will be used to generate the required electricity through the PEMFC to supply loads at deficit times, as so, the loads will be supplied in reliable way, as will be illustrated in the following sections.
3.1. Design of the Wind Turbines Energy
Different wind turbines (WTs) which are suitable for supplying the EDS were optimized considering both technical and economical points of view . The generation model and economics of the WT are discussed as follow:
3.1.1. Generation Model of the Wind Energy
The average output of WT is given as follows :
where; Pave is the average power, Pr is the rated power, Vcin is the cut-in wind speed, Vr is the rated wind speed and Vcout is the cut-out wind speed, C and K are the scale and shape parameters of Weibull probability distribution function, respectively. The average output of the wind farm is given as follows:
3.1.2. The Proposed Economical Model of the Wind Farm
The unit energy cost (UE.C) of WT is calculated dependent on the value of both the annual capital cost (AC.C) and the annual operation cost (AO.C). The UE.C is given as:
where; AEWT is the energy output of WT per year. The total unit energy cost of the wind farm (UWF) is given as:
where; PLand is the price of land used for the installation of the wind farm.
3.1.3. Produced Amount of H2 from the Wind Farm
The amount of H2 that needs to be supplied the FC to produce one kWh of electricity from the surplus power of wind farm can be obtained as follow :
where, nrate is the rate of flow of H2 into the cell in (mol/sec).
3.2. Design of Proton Exchange Membrane Fuel Cell (PEMFC)
In this paper, PEMFC is designed at low operating temperature between [60˚C - 100˚C] that has relatively high electrical efficiency between [40% - 50%]. The design of PEMFC is performed by considering the voltage-current electrical characteristics of the PEMFC as shown in Figure 1. The maximum power is corresponding to the operation of PEMFC between 0.6 and 0.8 Volt per cell . By fitting the V-I curve, the following equation can be obtained as :
where, Vcell, I and A are the single cell voltage, the current and the single cell area, respectively.The stack of the PEMFC can be expressed as follow
Figure 1. Voltage-current curve for a typical FC .
where; Vstack, Pstackand N are the stack voltage, the power of the stack and the number of stack series cells, respectively.
3.2.1. Water Electrolyzer
The electrolyzer size is determined by its production flow rate (H2) as illustrated :
where; H2rate is the H2 flow rate, EL is the average electrical load in Wh/day, ne is the number of electrons per second for 1 ampere, gm is the molar mass of H2 in g/mol, Uf is the utilization factor which takes 0.8, Vcell is the cell voltage, nem is the number of electrons per each molecule of H2 and nmm is the number of molecules per H2 mol. The electrolyzer size can be calculated by dividing the total amount of H2 per year required from the PEMFC by the total windy hours per year as follow:
where; hwind is the total windy hours per year.
3.2.2. Hydrogen Tank
The volume of the H2 tank can be determined from the ideal gas low as follow :
where; VT is the volume of the H2 tank, n is the number of moles, R is the gas constant, T is the temperature in Kelvin and P is pressure.
3.2.3. Calculation of Unit Energy Cost of the Designed PEMFC System
The annual cost (AC) of PEMFC system can be calculated as follow:
where, CO&M is the annual operating and maintenance costs of the PEMFC system, Y is the life time of the PEMFC system, CC is the capital cost, CTank is the annual cost of H2 tank, CElectro is the annual cost of the electrolyzer, Sv is the annual salvage value, Crun is the annual running cost and i is the interest rate. The unit energy cost (COE) can be calculated by dividing the AC of the PEMFC system by the total generated energy, during the system cycle as follow:
where, Ee is the electrical generated energy in kWh.
The optimal hybrid WE/PEMFC system approach discussed in Section 3 is applied for supplying the EDS of Mersi-Matrouh City, in attempt to isolate it from the grid. Also, the proposed optimal allocation procedure of the hybrid WE/ PEMFC system in the EDS of Mersi-Matrouh City discussed in Section 2 is applied.
4.1. Design of Hybrid WE/PEMFC System for Mersi-Matrouh City
Load forecasting of the electrical peak load demand of Mersi-Matrouh City is obtained up to year 2022 , using the historical data shown in Table 1. The peak load demand of Mersi-Matrouh City in year 2022 is equal to 8655.3 MW. The proposed models in Section 3 are used to estimate the optimal technical and economical wind energy which is 3000 kW-WT for Mersi-Matrouh City by considering Equations (14)-(17). Figure 2 shows the EDS of Mersi-Matrouh City up to year 2022. Figure 3 gives the average daily power supplied to the load at Mersi-Matrouh City in year 2022. The characteristics of 3000 kW-WT (E101/3000) model is used to evaluate the output wind energy (EWT). The hybrid WE/PEMFC system consists of 45 × 3000 kW-WT, which gives 239,611.5 MWh per year. The wind farm gives surplus power of about 427.999 MWh which is used to produce 7,275,983 grams of H2 through an electrolyzer, where, QLC-180 PEM electrolyzer cell stack is used . PEMFC can produce electric energy of about 107.88 MWh that can be used to supply the deficit power of 67.6 MWh at the time of shortage power. So, this EDS can be isolated from the grid. These values are illustrated in Table 2.
The proposed economical model in section 3.1.2 is applied to calculate the UE.C and UWF up to year 2022, as shown in Table 3. The proposed economical model of wind farm is performed according to the following assumptions :
● The capital cost (C.C) is assumed as 800 $/kW.
● The annual operation cost (AO.C) is assumed as 2.0 ¢/kWh.
Table 1. Historical electrical peak load data of Mersi-Matrouh City from 2005 to 2019.
Table 2. Energy generation of the hybrid WE/PEMFC system.
Table 3. Capital cost (C.C), the annual capital cost (AC.C), the annual operating cost (AO.C), the unit energy cost (UE.C) and the total unit energy cost of wind farm.
Figure 2. EDS of Mersi-Matrouh City up to 2022.
Figure 3. Wind farm average daily power supplied for Mersi-Matrough City up to year 2022 .
● The interest rate is assumed as 12%, while, the life-time is assumed as 20 years.
● The price of land used for wind farm is equal to 40,060$ per year.
PEMFC system is designed as illustrated in Section 3.2 and the design parameters of the PEMFC stacks are illustrated in Table 4. The design of the PEMFC stack is performed under the following considerations; FC running time is 3000 hrs/year, Vcell = 0.6V, Vstack = 72V, electorlyzer efficiency is 85% and PEMFC efficiency is 45% as Equations (19)-(24).
The COE of the designed PEMFC system is also calculated including H2 tank cost and electrolyzer cost as shown in Table 5 and performed under the following assumptions, as Equation (26):
● The CC of the PEMFC system is 4400$/kW.
● CO&M is assumed as 0.035$/kW.
● i is assumed as 5%.
● Crun and Sv are assumed as 9% of the CC.
To reduce the storage volume of the H2 gas, the pressure needs to increase. In this paper, the type four (type IV) tank is used . The specification of type IV hydrogen tank is shown in the Appendix. The limits of voltage magnitude are taken between 0.95 and 1.05 p.u, the maximum permissible number of DG units is taken as 16 unit, the maximum permissible active power flow is 7.5 MW for line code 1 and 4.9 MW. While, the maximum permissible reactive power flow are 3.63 MVAR for line code1 and 2.37 MVAR for line code 2.
4.2. Results and Discussions of Proposed Procedures
The proposed procedures are applied to the EDS of Mersi-Matrouh City in order to solve the optimal allocation of hybrid WE/PEMFC system problem by achieving the SOF/MOF and satisfying system constraints. Optimal allocating problem is solved with and without considering the SFs. Then these results are compared with each other.
4.2.1. Optimal Candidate Buses Using Sensitivity Factors (SFs)
Figure 3 includes four feeders, each feeder consists of four buses and five connecting cables, except the third feeder which includes only three buses and four connecting cables. This EDS is completely supplied from hybrid WE/PEMFC system that depends on WT which has rated of 3000 kW for each turbine combined with PEMFC stacks which has rated of 5 kW for each unit. Figure 4 and Figure 5 show the values of VSF, ALSF and RLSF, respectively. According to system buses ordering in these figures, the first seven buses of the bus-number are considered as the candidate buses for allocating the hybrid WE/PEMFC system.
Table 4. Design parameters of the PEMFC system for Mersi-Matrouh City.
Table 5. Economics of the PEMFC system for Mersi-Matrouh City.
Figure 4. The values of VSF for each system bus.
Figure 5. The values of both ALSF and RLSF for each transmission line section.
4.2.2. Minimizing the SOF without and with Considering Sensitivity Factors (SFs)
Minimizing the SOF is performed using the three proposed optimization techniques without and with considering the SFs. Table 6 and Table 7 show the optimal allocating for the hybrid system, and the corresponding number of the hybrid units (DGs) using BCSA, DJA and BPSO techniques in year 2022 without and with considering SFs, respectively. While, Table 8 shows a comparison among the fitness values using the three proposed optimization techniques with and without the hybrid system for the SOF without and with considering the SFs in year 2022.
4.2.3. Minimizing of F1 as a SOF
Table 8 illustrates the improvement in F1 without and with considering the SFs using the three proposed optimization techniques in year 2022. From Table 8, the best technique without considering the SFs is DJA that gives improvement in bus voltage about 75%, then, BCSA improves the bus voltage by 72.3%, finally BPSO technique improves the bus voltage by 71.5%. Also, the best technique is
Table 6. Sitting, sizing (S) in MW and corresponding number of DGs (N) using the three proposed optimization techniques without considering the SFs in year 2022 using SOFs.
Table 7. Sitting, sizing (S) in MW and corresponding number of DGs (N) using BCSA, DJA and BPO technique with considering the SFs in year 2022 using SOFs.
Table 8. A comparison of the fitness values using three proposed optimization techniques, with and without the hybrid WE/PEMFC system for the SOFs in year 2022 without and with considering the SFs.
DJA when considering the SFs that improves the bus voltage about 97.5%, then, BCSA improves the bus voltage by 97.3%, finally BPSO technique improves the bus voltage by 94.7%. There is a clear improving in the value of F1 when considering the effect of the SFs as F1 is increased from 75% to 97.5% using DJA. Figure 6 and Figure 7 show the improvement in voltage profile using the three proposed optimization techniques without and with considering the SFs for F1 as SOF, respectively.
4.2.4. Minimizing of F2 as a SOF
Table 8 illustrates the saving in ALI without and with considering the SFs. From this Table, the best technique without considering the SFs is BCSA that gives the highest saving in active power losses (APLs) by 77.2%, then, DJA gives saving in APLs by 77%, finally BPSO technique gives saving in APLs of about 75.5%. Also, the best technique when considering the SFs is DJA that gives the highest saving in APLs by 79.27%, then, both BCSA and BPSO technique that give the same saving in APLs about 79.09%. There is a significant improvement in the value of F2 when considering the SFs that the saving in APLs is increased from 77% to 79.27 using DJA.
4.2.5. Minimizing of F3 as a SOF
Table 8 illustrates the saving in F3 without and with considering the SFs using the three proposed optimization techniques, the best technique without considering the SFs is DJA that gives the highest saving in reactive power losses (RPLs) by 77.8% then, BCSA saves the RPLs by 77.12%, finally BPSO technique saves the RPLs by 76%. Also, the best technique, when considering the SFs, is DJA that gives the highest saving in RPLs by 79.55%, then, both BCSA and BPSO give the same saving in RPLs by 79.32%. There is a significant improving
Figure 6. Voltage improvement using the three optimization techniques without considering the SFs in year 2022 using F1 as SOF.
Figure 7. Voltage improvement using the three optimization techniques with considering the SFs in year 2022 using F1 as SOF.
in the value of F3 after taking the SFs into consideration which the saving in RPLs is increased from 77.8% to 79.55% using DJA.
4.2.6. Minimizing of MOF
MOF are dependent on the values of WFs, which takes the value of W1 as 0.5 to consider the voltage improvement is the first priority, while, the values of W2 and W3 take different values between [0 to 0.5 by step 0.05]. However, the values of W2 and W3 have taken as 0.45 and 0.05, respectively, to optimize the MOF using both BCSA and BPSO technique without and with considering the SFs. While, W2 and W3 are taken as 0.05 and 0.45, respectively, to optimize the MOF using DJA without and with considering the SFs. Table 9 shows the sitting and corresponding number of DGs using BCSA, DJA and BPSO technique using MOF without and with considering the SFs. A comparison of the fitness value using three proposed optimization techniques, with and without the hybrid WE/ PEMFC system using the MOF without and with considering SFs is shown in Table 10. From Table 10, DJA is the best technique without considering the SFs that minimizes the MOF from 67.07% (without inserting the hybrid system) to 18.07%, (after inserting the hybrid system), reduces the TAPLs from 1.2348 MW to 0.313 MW, and, reduces the TRPLs from 1.544 MVAR to 0.3812 MVAR. Also, DJA is the best technique with considering the SFs that minimizes the MOF from 67.07% (without inserting the hybrid system) to 17.22% (after inserting the hybrid system) reduces the TAPLs from 1.2348 MW to 0.2601 MW, and reduces the TRPLs from 1.5421 MVAR to 0.3231 MVAR. There is a significant minimization in the value of MOF with considering the effect of SFs, however, MOF is reduced from 18.07% to 17.22% using DJA.
Table 9. Sitting, sizing (S) in MW and corresponding number of DGs (N) using the three proposed optimization techniques in year 2022 using MOF without and with considering SFs.
Table 10. A comparison of the fitness value using three proposed optimization techniques, with and without the hybrid WE/PEMFC system using the MOF in year 2022 without and with considering SFs.
Figure 8. Voltage improvement using the three optimization techniques without considering the SFs in year 2022 using MOF.
Figure 9. Voltage improvement using the three optimization techniques with considering the SFs in year 2022 using MOF.
In this paper, a proposed hybrid WE/PEMFC system includes WTs and PEMFC has been efficiently presented. Proposed economical and generation models of 3000 kW-WT had been presented. Also, the economics of both wind energy and the PEMFC system had been discussed in year 2022. The BFS method has been used for the load flow calculation. A two-stage procedure has been introduced to allocate a hybrid WE/PEMFC system in the EDS for Mersi-Matrough City. In the first stage, the sensitivities analysis using three SFs have been presented to select the candidate buses for locating the hybrid system. While, in the second stage, the BCSA, DJA and BPSO techniques have been investigated to find the optimal allocation of the hybrid system considering the improvement of voltage profile and the minimization of both APLs and RPLs as objective functions, while system constraints are fully satisfied. Best optimization results have been obtained using DJA compared to BCSA and BPSO technique. Also, the proposed procedure using DJA considering MOF gives the best results in compared with the other two techniques. By considering the effect of SFs, the improvement in the value of the F1 is increased; however, the hybrid system is located at only 5 candidate system buses instead of the 16 buses. So, these SFs decrease the search space for optimization techniques and improve the voltage profile. Also, there are more significant saving in both F2 and F3 when the hybrid system are located at the candidate buses. The proposed procedure represents a potential tool to improve bus voltage magnitude, reduce both APLs and RPLs and helps their operators in smart grid environment.
Table A1. Key features of QLC PEM electrolyzer cell stack .
Table A2. Specification of type four hydrogen tank 500 bar-200L .
 Hoff, T.E., Wenger, H.J. and Framer, B.K. (1996) Distributed Generation an Alternative to Electric Utility Investments in System Capacity. Energy Policy, 24, 137-147.
 Allam, S.M., EL-Zeftawy, A.A. and Shehata, N.K. (2012) Operation of Electrical Utility with Application in Egypt. 8th International Conference on Electrical Engineering, 29-31 May 2012.
 Luo, X., Wang, J., Dooner, M. and Clarke, J. (2015) Overview of Current Development in Electrical Energy Storage Technologies and the Application Potential in Power System Operation. Applied Energy Journal, 137, 511-536.
 Azmy, A.M. and Imam, S.M. (2015) Sizing and Economic Analysis of Standalone PEM Fuel Cell Systems for Residential Utilization. Recent Innovations in Mechatronics, Vol. 2.
 Wang, Y., Diaz, D.F.R., Chen, K.S., Wang, Z. and Adroher, X.C. (2020) Materials, Technological Status, and Fundamentals of PEM Fuel Cells: A Review. Materials Today, 32, 178-203.
 Godwin, U.O., Callistus, O.C. and Paul, I.J. (2013) A Power Estimation Model for Wind—Hydrogen-Fuel Cell—Inverter Hybrid Power Generation System. International Journal of Advanced Research in Electrical, Electronics and Instrumentation Engineering, 2, 2474-2485.
 Das, D., Esmaili, R., Xu, L. and Nichols, D. (2005) An Optimal Design of a Grid Connected Hybrid Wind/Photovoltaic/Fuel Cell System for Distributed Energy Production. 31st Annual Conference of IEEE Industrial Electronic Society, Raleigh, 6-10 November 2005.
 Zakaria, E.M., Arafa, S.H., Nashed, M.F. and Ramadan, S.G. (2020) Fuzzy Logic Control Management with Stand Alone Photovoltaic—Fuel Cell System. ASTES Journal, 5, 424-430.
 Shao, Z. and Lo, K.L. (2017) A Case Study to Determine Optimal Capacity of a Stand-Alone Wind/PV/Battery System. Energy and Power Engineering Journal, 9, 300-307.
 Traore, A., Elgothamy, H. and Zohdy, M.A. (2018) Optimal Sizing of Solar/Wind Hybrid Off-Grid Microgrids Using an Enhanced Genetic Algorithm. Journal of Power and Energy Engineering, 6, 64-77.
 Das, P. (2015) Optimal Sizing and Placement of Distributed Generation in A Radial Distribution System Using Loss Sensitivity Factor and Firefly Algorithm. Intentional Journal of Scientific Research and Management, 3, 2611-2618.
 Becerra, J.E.C. and Riano, H.E.H. (2015) Distributed Generation Placement in Radial Distribution Networks Using a Bat-Inspired Algorithm. Dyna, 82, 60-67.
 Wang, L. and Singh, C. (2008) Reliability-Constrained Optimum Placement of Reclosers and Distributed Generators Distribution Networks Using Ant Colony System Algorithm. IEEE Transactions on Systems, Man, and Cybernetics Part C: Application and Reviews, 38, 757-764.
 Sutthibun, T. and Bhasaputra, P. (2010) Multi-Objective Optimal Distributed Generation Placement Using Simulated Annealing. Proceedings ECTI-CON International Conference on Electrical Engineering/Electronics Computer Telecommunications and Information Technology, 810-813.
 Aly, A.I., Hegazy, Y.G. and Alsharkawy, M.A. (2010) A Simulated Annealing Algorithm for Multi-Objective Distributed Generation Planning. Proceedings IEEE Power and Energy Society General Meeting, Providence, 25-29 July 2010, 1-7.
 Ghadimi, N. and Ghadimi, R. (2012) Optimal Allocation of Distributed Generation and Capacitor Banks in order to Loss Reduction in Reconfigured System. Research Journal of Applied Sciences Engineering and Technology, 4, 1099-1104.
 El-Ela, A.A.A., El-Sehiemy, R.A., Kinawy, A.M. and Mouwafi, M.T. (2016) Optimal Capacitor Placement in Distribution Systems for Power Loss Reduction and Voltage Profile Improvement. IET Generation, Transmission & Distribution, 10, 1209-1221.
 Askarzadeh, A. (2016) A Novel Metaheuristic Method for Solving Constrained Engineering Optimization Problems: Crow Search Algorithm. Computers and Structures, 169, 1-12.
 Souza, R.C.T., Macedo, C.A., Coelho, L.S. and Pierezan, J. (2018) A V-Shaped Binary Crow Search Algorithm for Feature Selection. 2018 IEEE Congress on Evolutionary Computation, Rio de Janeiro, 8-13 July 2018.
 Rao, R.V. (2016) Jaya: Simple and New Optimization Algorithm for Solving Constrained and Unconstrained Optimization Problems. International Journal of Industrial Engineering Computations, 7, 19-34.
 El-Ela, A.A.A., Allam, S.M. and Shehata, N.K. (2019) Optimal Allocation of a Hybrid Renewable Energy-Storage System for Supplying Egyptian Distribution Networks Using Discrete Jaya Algorithm. MEPCON’21, Tanta, 19-21 December 2019.
 El-Ela, A.A.A., Allam, S.M. and Shehata, N.K. (2017) Assessment of Hybrid Renewable Energy with Energy Storage System for Supplying Distribution Networks in Egypt. MEPCON’19, Shibin el Kom, 19-21 December 2017.
 Aly, G.E.M., El-Gharabawy, A.E. and El-Zeftawy, A.A. (2001) Optimizing the Operation of Wind Generators at the Installation Site. Proceedings of the 8th International Middle East Power Conference, MEPCON’2001, Cairo, 29-31 December 2001, 951-955.
 Azmy, A.M. and Imam, S.M. (2015) Sizing and Economic Analysis of Standalone PEM Fuel Residential Utilization. Recent Innovations in Mechatronics, Debrecen, Volume 2.
 El-Ela, A.A.A., Allam, S.M. and Said Ahmed, M.M. (2016) Reliable Operating Assessment for Hybrid PV/Wind Energy Generation System to Supply New City in Egypt. 18th International Middle East Power Systems Conference, MEPCON’18, Helwan, 27-29 December 2016.