Developing a Hybrid Optimization Algorithm for Optimal Allocation of Renewable DGs in Distribution Network

Distributed generation (DG) is becoming a prominent key spot for research in recent years because it can be utilized in emergency/reserve plans for power systems and power quality improvement issues, besides its drastic impact on the environment as a greenhouse gas (GHG) reducer. For maximizing the benefits from such technology, it is crucial to identify the best size and location for DG that achieves the required goal of installing it. This paper presents an investigation of the optimized allocation of DG in different modes using a proposed hybrid technique, the tunicate swarm algorithm/sine-cosine algorithm (TSA/SCA). This investigation is performed on an IEEE-69 Radial Distribution System (RDS), where the impact of such allocation on the system is evaluated by NEPLAN software.


Introduction
Recently, the integration of distributed generation (DG) in distribution networks is becoming very popular to meet increases in system load [1][2][3][4]. Also, the world is interested in installing several types of DG, especially renewable sources in power grids, such as hydropower, biomass, photovoltaic (PV), and wind turbine [5][6][7] technologies. DG is a newly coined term that describes a technology that has a deep impact on power systems nowadays. A few years ago, distributed generation (DG)-or dispersed generation-started to arise in the world of power systems, aiming to exploit small-scale energy resources to utilize them in electric power generation instead of depending only on centralized largescale power generation stations [1][2][3][4]. DG stations are characterized by being small-scale (usually less than 50 MW) and installed directly to the distribution power system, instead of the traditional transmission power system, enabling the facility that owns such a station to consume a part of the generated power, then export the surplus power [1,2].
Such a technology brings renewable energy resources into action, as most of the renewable-energy-dependent generation stations are small-scale stations, and of course, this will have a drastic effect on the environment. In other words, the more penetration of renewable-energy-dependent generation stations is achieved, the more reduction of greenhouse gasses (GHG) takes place [3][4][5][6][7][8].
Besides the environmental impact of distributed generation, it also has a remarkable impact on power-quality issues in distribution systems. It has a relieving effect on congested transmission and distribution systems due to its unique location: just beside the consumer! Such an advantage results in economic and environmental benefits by reducing the power losses of the system, as there is no need for additional transmission lines, while saving such losses results in a reduction of the GHG effect by about 1% [4].
On the other hand, distributed generation can be utilized for improving power quality in distribution systems, and even in transmission systems, as it has the ability to be con-nected at different voltage levels, i.e., low, medium, and high voltage [1]. It can improve system reliability, maintain voltage stability, and also provide the system with reserve generated power for emergencies [1,5]. Moreover, such technology offers a remarkable contribution to economic investment in energy, inspiring researchers to introduce several techniques for the initiation and penetration of Distributed Energy Resources (DER) technologies and enabling a vision of such penetration regarding economic investment and cost reduction, focusing on biogas and hydrogen cells [9][10][11][12].
Nevertheless, identifying the perfect allocation and size of a distributed generator is a critical issue, as it increases the profit of installing it, which is the goal of installation. Obtaining the best location and size has a deep influence on the impact of installing DG. In [6], the optimal allocation of DG was identified with regards to the operations and investment costs of DG, where the utilized optimization approach tends to decrease the energy loss when the loads are time-varying by determining their generation capacity at different instants. Another optimization method is proposed in [7], regarding the minimization of power losses, where it is simplified so that it does not need excessive computational processes. In [8], a new optimization method, typically the virus colony search (VSC) algorithm was utilized, considering the reliability assessment of the distribution network, where the results are compared with several optimization methods. A hybrid method is proposed in [9] regarding loss minimization and voltage improvement, where the location of DG is identified using an empirical discrete metaheuristic (EDM), while the size is identified by the steepest descent method (SD).
This paper introduces a new hybrid optimization approach wherein the tunicate swarm algorithm (TSA) is merged with the sine-cosine algorithm (SCA), resulting in a novel TSA/SCA hybrid approach. This new approach is used to identify the best size and location of DG in distribution, considering the minimization of system loss. The optimization process is performed in three scenarios: (1) the DG is producing real power only (P-type case), (2) the DG is producing reactive power only (Q-type case), and (3) the DG is generating both reactive and real power (PQ-type case). The performance of the RDS is evaluated after the optimization process using NEPLAN software.
The contributions of this paper are the (1) introduction of a new hybrid approach that consists of the TSA and SCA algorithms, (2) use of an efficient hybrid approach to determine the optimal planning of DG in RDS, and (3) study of the effect of integrating different types of DG in RDS.
The remainder of the paper is organized as follows: the presented problem is explained in Section 2; sensitivity is discussed in Section 3; Section 4 explains the presented algorithm and the obtained results. Section 5 discusses the conclusion. Figure 1 shows the representation of two buses in a distribution system. The system power flow is evaluated in a backwards direction by Equations (1) and (2) [24]:

Problem Formulation
Then, the voltage magnitude of bus (r) can be determined in a forward direction as follows: The problem formulation can be presented as a multiobjective function as follows: where, where, VD(m) represents the voltage deviation at bus (m), VSI(m) represents the voltage stability index at bus (m), S is the total system buses, and B is the total system branches. K 1 , K 2 , and K 3 are weighting factors that are equal to 0.5, 0.25, and 0.25, respectively. The inequality and equality constraints are determined as shown next [19,25,26].
Clean Technol. 2021, 3, FOR PEER REVIEW 3 Figure 1 shows the representation of two buses in a distribution system.

Problem Formulation
Then, the voltage magnitude of bus (r) can be determined in a forward direction as follows: The problem formulation can be presented as a multiobjective function as follows:

Equality Constraints
The power flow balance equation that can be represented by (9) and (10), and the power flow equation that can be represented by (11) and (12), are the inequality constraints as shown next.
where, Q S and P S are the output reactive and real power from the grid, respectively, but Q L,m and P L,m are the reactive and real load demand at bus (m), respectively. M and N are the total number of PVs and capacitors in RDS, respectively. P PV (m) and Q capacitor (m) are the output power of PVs and capacitors at bus (m), respectively.

Inequality Constraints
The system operation constraints are the inequality constraints, which can be represented as follows:

System Voltage Constraints
The bus system voltage is operating between the minimum operating voltage (V Down ) and the maximum operating voltage (V Up ).
DER Sizing Limits: P PV, n ≤ P PV ≤ P PV,a The PV output is operating between the minimum (P PV,n ) and maximum power (P PV,a ) of PVs in RDS.
Capacitor Size Limits:

Line Capacity Limits
The branches current of the system is operating under operating constraints: where, I a,m represents the high operating current in the branch (m).

Tunicate Swarm Algorithm (TSA)
Tunicates are cylinder-shaped, with a gelatinous tunic that is closed at one end and open at the other. Tunicates shine with bioluminescence, generating a faint green-blue light, which can be viewed from several meters away. The size of tunicates is a few millimeters. In the ocean, tunicates absorb water to generate jet propulsion from their open ends using atrial siphons. Tunicates move in water by generating jet propulsion. The updating position of tunicates can be formulated as follow: Tunicates move in a swarm in nature, which can be modeled by the following equation:

Sine-Cosine Algorithm (SCA)
The SCA is derived from the cosine and sine function to create an effective optimization algorithm. The effectiveness of the SCA is based on its exploitation and exploration phases. r 1 is used to balance the exploitation and exploration rates for the SCA over the course of iterations to obtain the global optimum solutions. The position of the SCA is updated as follows: Where, T max and T are the maximum and current iteration, respectively, P i T is the targeted global optimal solution, and X i T represents the current iteration. r 1 , r 2 , r 3 , and r 4 are random numbers.

Improved TSA-SCA Algorithm
The improved TSA-SCA is created by applying the updating position of the SCA to the updating position of the TSA to improve the exploration phase of the TSA. The rest of the pseudo code of the TSA remains the same, as shown in the following equations: The steps of the TSA-SCA to determine the optimal sizes and locations of DG in distribution networks are explained in the following steps: 1.
Read the system data, maximum iteration (I), and number of search agents (S).

2.
Produce the initial population of slime mold between the lower-(w) and upper (p)controlled variables by Equation (36).
where, rand represents a random value between the values of 0 and 1. c and q are the number of tunicates and the problem dimension. 3.
The produced population represents the tunicate position that can be formulated as follows: where, S i,j is the position of the tunicate.

4.
Evaluate the fitness for all locations of tunicates, and obtain the superior position of tunicates and the superior objective function.
Return to step 4 until the final iteration is reached. 7.
Obtain the best location of tunicates (sizes and positions of DG).

Testing and Evaluation
As mentioned in Section 1, the proposed optimization method is applied to the IEEE 69node RDS to identify the optimal sizing and allocation of DG in such a system, considering the minimization of power losses. Such a system is tested by the proposed method when installing one, two, and three DGs to it. Also, the DG mode changes, i.e., DG is tested in three modes: while generating active power only (P-type mode), while generating reactive power only (Q-type mode), and while generating both active and reactive power. The presented test system is a IEEE 69-node RDS that includes load demand of 3801.49 kW and 2694.6 kVAR [23]. This system consists of 69 buses and 68 branches with a base of 12.66 kV and 100 MVA. Figure 2 shows the system implementation by NEPLAN software. The analysis of the IEEE 69-node RDS system in its original case shows the active power losses of the system (PLosses = 225 kW), and the least voltage is in node 65, where the voltage of that node is V65 = 0.9092 pu. Figure 3 shows the loading ability of the system by PV curves, where PV curves are applied on buses 10, 27, 46, and 65. The analysis of the IEEE 69-node RDS system in its original case shows the active power losses of the system (P Losses = 225 kW), and the least voltage is in node 65, where the voltage of that node is V 65 = 0.9092 pu. Figure 3 shows the loading ability of the system by PV curves, where PV curves are applied on buses 10, 27, 46, and 65. The analysis of the IEEE 69-node RDS system in its original case shows the active power losses of the system (PLosses = 225 kW), and the least voltage is in node 65, where the voltage of that node is V65 = 0.9092 pu. Figure 3 shows the loading ability of the system by PV curves, where PV curves are applied on buses 10, 27, 46, and 65. The PV curves show that the voltage of node 65-the weakest node in the RDSstarts to fall below 90% of nominal voltage when loads exceed 110%, and the system collapses when loads exceed 320%. Now the hybrid TSA-SCA method can start to be used to optimize the system by identifying the optimal allocation and sizing of DG installed regarding active power losses. After the optimization process, the system is analyzed to show the impact of the optimization process on the performance of the system. The optimization process is achieved using MATLAB R2019b on a personal computer with 2 GB RAM and Intel(R) Pentium(R) CPU B950, 2.1GHZ, while the analyses are achieved using NEPLAN software. The PV curves show that the voltage of node 65-the weakest node in the RDS-starts to fall below 90% of nominal voltage when loads exceed 110%, and the system collapses when loads exceed 320%. Now the hybrid TSA-SCA method can start to be used to optimize the system by identifying the optimal allocation and sizing of DG installed regarding active power losses. After the optimization process, the system is analyzed to show the impact of the optimization process on the performance of the system. The optimization process is achieved using MATLAB R2019b on a personal computer with 2 GB RAM and Intel(R) Pentium(R) CPU B950, 2.1 GHZ, while the analyses are achieved using NEPLAN software.

System Optimization with Active Power-Generating DG (P-Type Mode)
In this case, the DGs installed are P-type and are configured for system optimization regarding the objective function as it is shown in Table 1. According to the settings of DGs shown in Table 1, the load flow process and PV curves were carried out. Starting with the system with one DG installation, results show power loss (Ploss = 83.2224 kW), and the total power loss reduction = 63.012%. The least voltage is found in node 27, where V 27 = 0.9683 pu. Figure 4 shows the system's PV curves with one P-type DG installation.

Power Losses Reduction (%)
63.012 68.145 69.144 According to the settings of DGs shown in Table 1, the load flow process and PV curves were carried out. Starting with the system with one DG installation, results show power loss (Ploss = 83.2224 kW), and the total power loss reduction = 63.012%. The least voltage is found in node 27, where V27 = 0.9683 pu. Figure 4 shows the system's PV curves with one P-type DG installation. Due to the new DG's location, the PV curve of node 65 improved, as shown in Figure  4. Figure 4 shows that the voltage of node 27-the system's weakest node in that casestarts to fall below 90% of nominal voltage when loads exceed 290%, and the system collapses when loads exceed 750%.
When installing two P-type DGs to the system according to the settings in Table 1, the active power loss of the system is PLoss = 71.6745 kW, and the total active power loss reduction is 68.145%. The least voltage is found in node 65, where V65 = 0.9789 pu. Figure  5 shows the system's PV curves with two P-type DG installations. Due to the new DG's location, the PV curve of node 65 improved, as shown in Figure 4. Figure 4 shows that the voltage of node 27-the system's weakest node in that case-starts to fall below 90% of nominal voltage when loads exceed 290%, and the system collapses when loads exceed 750%.
When installing two P-type DGs to the system according to the settings in Table 1, the active power loss of the system is P Loss = 71.6745 kW, and the total active power loss reduction is 68.145%. The least voltage is found in node 65, where V 65 = 0.9789 pu. Figure 5 shows the system's PV curves with two P-type DG installations. The presence of the two DGs in the estimated locations improves the PV curves even more. Figure 5 shows that the voltage of node 65-the weakest node in that case-starts to fall below 90% of nominal voltage when loads exceed 380%, and the system collapses when loads exceed 750%.
When installing three P-type DGs to the system according to the settings of Table 1, the active power loss of the system is PLoss = 69.4266 kW, and the total active power loss reduction is 69.144%. The least voltage is found in node 65, where V65 = 0.979 pu. Figure 6 presents the system's PV curves with three P-type DG installations. The presence of the two DGs in the estimated locations improves the PV curves even more. Figure 5 shows that the voltage of node 65-the weakest node in that case-starts to fall below 90% of nominal voltage when loads exceed 380%, and the system collapses when loads exceed 750%.
When installing three P-type DGs to the system according to the settings of Table 1, the active power loss of the system is P Loss = 69.4266 kW, and the total active power loss reduction is 69.144%. The least voltage is found in node 65, where V 65 = 0.979 pu. Figure 6 presents the system's PV curves with three P-type DG installations. more. Figure 5 shows that the voltage of node 65-the weakest node in that case-starts to fall below 90% of nominal voltage when loads exceed 380%, and the system collapses when loads exceed 750%.
When installing three P-type DGs to the system according to the settings of Table 1, the active power loss of the system is PLoss = 69.4266 kW, and the total active power loss reduction is 69.144%. The least voltage is found in node 65, where V65 = 0.979 pu. Figure 6 presents the system's PV curves with three P-type DG installations. Figure 6. PV curves of the IEEE 69-node RDS system with 3 P-type DG installations. Figure 6 shows that there is almost no drastic effect on the system's PV curves when installing 3 P-type DGs. The voltage of node 65-the system's weakest node in that casestarts to fall below 90% of nominal voltage when loads exceed 380%, and the system collapses when loads exceed 750%. It is also observed that the presence of the third DG resulted in an improvement on node 10′s PV curve.  Figure 6 shows that there is almost no drastic effect on the system's PV curves when installing 3 P-type DGs. The voltage of node 65-the system's weakest node in that casestarts to fall below 90% of nominal voltage when loads exceed 380%, and the system collapses when loads exceed 750%. It is also observed that the presence of the third DG resulted in an improvement on node 10 s PV curve.

System Optimization with Reactive Power-Generating DGs (Q-Type Mode)
In this case, the DGs installed are Q-type DGs and are configured for system optimization regarding the objective function as it is shown in Table 2. Following the settings of DGs shown in Table 2, the load flow process and PV curves were carried out. When installing one Q-type DG, the power loss of the system is P Loss = 152.041 kW, and the total power loss reduction = 32.426%. The least voltage is found in node 65, where V 27 = 0.9307 pu. Figure 7 presents the system's PV curves with one Q-type DG installation.
In the presence of one Q-type DG, the system's PV curves improved. Figure 7 shows that the voltage of node 65-the system's weakest node in that case-starts to fall below 90% of nominal voltage when loads exceeds 140%, and the system collapses when loads exceeds 400%.
When installing two Q-type DGs to the system according to the settings in Table 2, the active power loss of the system is P Loss = 146.441 kW, and the total active power loss reduction is 34.915%. The least voltage is found in node 65, where V 65 = 0.9311 pu. Figure 8 presents the system's PV curves with two Q-type DG installations.

Power Losses Reduction (%)
32.426 34.915 35.498 Following the settings of DGs shown in Table 2, the load flow process and PV curves were carried out. When installing one Q-type DG, the power loss of the system is PLoss = 152.041 kW, and the total power loss reduction = 32.426%. The least voltage is found in node 65, where V27 = 0.9307 pu. Figure 7 presents the system's PV curves with one Q-type DG installation. In the presence of one Q-type DG, the system's PV curves improved. Figure 7 shows that the voltage of node 65-the system's weakest node in that case-starts to fall below 90% of nominal voltage when loads exceeds 140%, and the system collapses when loads exceeds 400%.
When installing two Q-type DGs to the system according to the settings in Table 2, the active power loss of the system is PLoss = 146.441 kW, and the total active power loss reduction is 34.915%. The least voltage is found in node 65, where V65 = 0.9311 pu. Figure  8 presents the system's PV curves with two Q-type DG installations.  Figure 8 shows that there is almost no drastic effect on the system's PV curves when installing 2 P-type DGs. The voltage of node 65-the system's weakest node in that casestarts to fall below 90% of nominal voltage when loads exceed 140%, and the system collapses when loads exceed 400%.
When installing three Q-type DGs to the system according to the settings of Table 2, the active power loss of the system is PLoss = 145.129 kW, and the total active power loss reduction is 35.498%. The least voltage is found in node 65, where V65 = 0.9314 pu. Figure  9 presents the system's PV curves with three Q-type DG installations.  Figure 8 shows that there is almost no drastic effect on the system's PV curves when installing 2 Q-type DGs. The voltage of node 65-the system's weakest node in that case-starts to fall below 90% of nominal voltage when loads exceed 140%, and the system collapses when loads exceed 400%.
When installing three Q-type DGs to the system according to the settings of Table 2, the active power loss of the system is P Loss = 145.129 kW, and the total active power loss reduction is 35.498%. The least voltage is found in node 65, where V 65 = 0.9314 pu. Figure 9 presents the system's PV curves with three Q-type DG installations.
Again, Figure 9 shows that there is almost no drastic effect on the system's PV curves when installing 3 Q-type DGs. The voltage of node 65-the system's weakest node in that case-starts to fall below 90% of nominal voltage when loads exceed 140%, and the system collapses when loads exceed 400%. installing 2 P-type DGs. The voltage of node 65-the system's weakest node in that casestarts to fall below 90% of nominal voltage when loads exceed 140%, and the system collapses when loads exceed 400%.
When installing three Q-type DGs to the system according to the settings of Table 2, the active power loss of the system is PLoss = 145.129 kW, and the total active power loss reduction is 35.498%. The least voltage is found in node 65, where V65 = 0.9314 pu. Figure  9 presents the system's PV curves with three Q-type DG installations. Again, Figure 9 shows that there is almost no drastic effect on the system's PV curves when installing 3 Q-type DGs. The voltage of node 65-the system's weakest node in that case-starts to fall below 90% of nominal voltage when loads exceed 140%, and the system collapses when loads exceed 400%.

System Optimization with Active and Reactive Power-Generating DG (PQ-Type Mode):
In this case, the DGs installed are PQ-type DGs and are configured for system optimization regarding the objective function as it is shown in Table 3.

System Optimization with Active and Reactive Power-Generating DG (PQ-Type Mode)
In this case, the DGs installed are PQ-type DGs and are configured for system optimization regarding the objective function as it is shown in Table 3. Following the settings of DGs shown in Table 3, the load flow process and PV curves were carried out. When installing one PQ-type DG, the power losses, P Loss = 23.169 kW, and the total power loss reduction = 89.702%. The least voltage is found in node 27, where V 27 = 0.9725 pu. Figure 10 presents the system's PV curves with one PQ-type DG installation.
Clean Technol. 2021, 3, FOR PEER REVIEW 12  Following the settings of DGs shown in Table 3, the load flow process and PV curves were carried out. When installing one PQ-type DG, the power losses, PLoss = 23.169 kW, and the total power loss reduction = 89.702%. The least voltage is found in node 27, where V27 = 0.9725 pu. Figure 10 presents the system's PV curves with one PQ-type DG installation. Due to the presence of the powerful PQ-type DG, the system's PV curves acquired superior improvement, as the DG provides both active and reactive power to the RDS. Figure 10 shows that the voltage of node 27-the system's weakest node in that casestarts to fall below 90% of nominal voltage when loads exceed 310%, and the system col- Figure 10. PV curves of the IEEE 69-node RDS system with 1 PQ-type DG installation.
Due to the presence of the powerful PQ-type DG, the system's PV curves acquired superior improvement, as the DG provides both active and reactive power to the RDS. Figure 10 shows that the voltage of node 27-the system's weakest node in that case-starts to fall below 90% of nominal voltage when loads exceed 310%, and the system collapses when loads exceed 850%.
When installing two PQ-type DGs to the system according to the settings in Table 3, the active power loss of the system is P Loss = 7.2013 kW, and the total active power loss reduction is 96.799%. The least voltage is found in node 69, where V 69 = 0.9943 pu. Figure 11 presents the system's PV curves with two PQ-type DG installations.
Due to the presence of the two PQ-type DGs, the system's PV curves were improved even more. Figure 11 shows that the voltage of node 65-the system's weakest node in that case-start to fall below 90% of nominal voltage when loads exceed 450%, and the system collapses when loads exceed 850%.
When installing three PQ-type DGs to the system according to the settings of Table  3, the active power loss of the system is PLoss = 4.2665 kW, and the total active power loss reduction is 98.104%. The least voltage is found in node 65, where V65 = 0.997 pu. Figure  12 presents the system's PV curves with three PQ-type DG installations. Figure 12. PV curves of the IEEE 69-node RDS system with 3 PQ-type DG installations. Figure 12 shows that a slight increase in the system's load capacity happens when installing 3 PQ-type DGs in the system. The voltage of node 65-the system's weakest node in that case-starts to fall below 90% of nominal voltage when loads exceed 460%, Figure 11. PV curves of the IEEE 69-node RDS system with 2 PQ-type DG installations.
Due to the presence of the two PQ-type DGs, the system's PV curves were improved even more. Figure 11 shows that the voltage of node 65-the system's weakest node in that case-start to fall below 90% of nominal voltage when loads exceed 450%, and the system collapses when loads exceed 850%.
When installing three PQ-type DGs to the system according to the settings of Table 3, the active power loss of the system is P Loss = 4.2665 kW, and the total active power loss reduction is 98.104%. The least voltage is found in node 65, where V 65 = 0.997 pu. Figure 12 presents the system's PV curves with three PQ-type DG installations. Figure 12 shows that a slight increase in the system's load capacity happens when installing 3 PQ-type DGs in the system. The voltage of node 65-the system's weakest node in that case-starts to fall below 90% of nominal voltage when loads exceed 460%, and the system collapses when loads exceed 850%. Even more, the PV curve of node 10 is slightly improved more than in the case of two PQ-mode DGs. Table 4 shows the comprehensive results between TSA-SCA techniques and several other techniques. even more. Figure 11 shows that the voltage of node 65-the system's weakest node in that case-start to fall below 90% of nominal voltage when loads exceed 450%, and the system collapses when loads exceed 850%. When installing three PQ-type DGs to the system according to the settings of Table  3, the active power loss of the system is PLoss = 4.2665 kW, and the total active power loss reduction is 98.104%. The least voltage is found in node 65, where V65 = 0.997 pu. Figure  12 presents the system's PV curves with three PQ-type DG installations.  Figure 12 shows that a slight increase in the system's load capacity happens when installing 3 PQ-type DGs in the system. The voltage of node 65-the system's weakest node in that case-starts to fall below 90% of nominal voltage when loads exceed 460%, and the system collapses when loads exceed 850%. Even more, the PV curve of node 10 is slightly improved more than in the case of two PQ-mode DGs. From Table 4, the improved TSA-SCA has efficient characteristics to obtain the best results when compared to other efficient algorithms.

Conclusions
This paper introduces a new hybrid approach, a mixture between the tunicate swarm algorithm (TSA) and the sine-cosine algorithm (SCA) optimization technique. The new TSA/SCA optimization technique was tested on an IEEE 69-node RDS system, where the fitness is decreasing the active power loss through identifying the optimal sizing and allocation of DGs installed on the system. The optimization process took place in three modes: with active power DGs (P-mode DGs), reactive power DGs (Q-mode DGs), and with both active and reactive power DGs (PQ-mode DGs). The performance of the optimized system was evaluated by NEPLAN software to show the impact of the optimization process on the system. The analyses shows that the objective function was successfully achieved by the proposed technique, the active power loss was obviously minimized and the load demand of the system was greatly increased so that it can withstand more loads. The results proved that integration of multiple DGs gives better results than the integration of a single DG in a distribution network. Installing PQ-type optimization gives better results than the integration of P-type or Q-type.