Butterfly optimizer-assisted optimal integration of REDG units in hybrid AC/DC distribution micro-grids based on minimum operational area

This paper presents the impact of optimal location and sizing of renewable and non-renewable-based distributed generators in the AC/DC micro-grid system using the latest optimizer called butterfly optimization algorithm with an aim to minimize power loss. Generally, hybrid AC/DC micro-grids systems are modeled by separating AC and DC feeders with the help of high-power converters (HPC).AC grids sustained by substation and DC grids are maintained by their individual DG units. While planning of DGs in the hybrid AC/DC systems, the power loss incurred by HPCs is not considered avoiding complexity by many authors. In this paper, the sizing of DGs is determined by the operational area required by the type of DG technology as one variable and all possible candidate buses in the respective zones of AC/DC micro-grid system are another variable with due consideration of HPC losses in AC/DC micro-grid system. A hybrid AC/DC MG system is developed by classifying the existing benchmark 33-bus and 69-bus radial distribution systems into various AC/DC zones. To evaluate the proposed approach, it is implemented on aforementioned micro-grid systems and the obtained results are verified with other existing approaches in the literature. The results proved that the proposed approach is better than the other approaches in technical aspects.

grid through HPC. Fuel cells can potentially be integrated with solar-PV technology to provide zero-emissions alternatives to fossil fuels. The energy generation cost of solar-PV systems continues to decrease year by year. By extension, hydrogen technologies that run with solar power also become cheaper to operate. Micro-grids with energy storage devices act as energy hubs that can store both electrical and thermal energy to serve the load demand [26]. A two-stage methodology is used to determine optimal locations based on LRSF (loss reduction sensitive factor) and then optimal sizing of REDGs in prefixed locations by using HNMCS algorithm with an objective of minimization of power loss by considering area required by DG type to calculate DG size as variable in AC/DC micro-grids [27] with neglected HPC losses and operating efficiency. In AC/DC networks, HPCs will have a significant role that can control bidirectional power flows. A typical HPC will operate with an average efficiency of 90% remaining 10% is considered as conversion losses [28].
From the literature, it is evident that the optimal allocation of DGs in a distribution system will definitely improve the system performance. Hence, the selection of suitable DG type followed by identifying the optimal size of DG and its placement is a nonlinear optimization problem that can be solved by proper formulation of the objective function and efficient optimizer. This paper presents a methodology to determine optimal allocation of REDGs in AC/DC micro-grid system with an objective of power loss minimization which includes HPC conversion losses by considering not only area required by REDGs to compute the size of REDG but also locations for the integrating REDGs simultaneously as variables using efficient and novel optimizer called butterfly optimizer (BO) to satisfy the AC/DC loads under various system constraints. The proposed methodology is tested on small-and medium-scale hypothetical AC/DC micro-grids under different cases, and obtained results are compared with the existing results from the literature.

AC load flow model
AC load flow algorithm is performed to calculate the voltage magnitudes of the buses and branch currents of the system. Traditional backward/forward sweep-based load flow has taken as a load flow algorithm.
AC micro-grid with N AC number of buses, ith bus current I i,AC is given by where P Li,AC , Q Li,AC are the active and reactive power demand at ith bus and V i,AC is the voltage magnitude of an ith bus. During backward sweep of the backward/forward based load flow, the branch currents are obtained with the help of bus-injections to branch currents matrix (BIBC) as follows During forward sweep of the backward/forward based load flow, the bus voltages are calculated with the help of branch current to bus voltage matrix (BCBV) as follows: (1) (2) J k,AC = BIBC * I i k = 1, 2, . . . , N AC − 1 Repeat Eqs. (2) and (3) until the difference between the voltages of two adjacent iterations is less than the tolerance value (ε).

DC load flow model
DC micro-grid with N DC number of buses, ith bus current I i,DC is given by where P Li,DC , P REDG , P Li,DCeff are the active power demand, REDG active power and effective active power demand at ith bus, respectively, and V i,DC is the voltage magnitude of an ith bus.
During backward sweep of the backward/forward based load flow, the branch currents are obtained with the help of bus-injections to branch currents matrix (BIBC) as follows: During forward sweep of the backward/forward based load flow, the bus voltages are calculated with the help of branch current to bus voltage matrix (BCBV) as follows: Repeat Eqs. (7) and (8) until the difference between the voltages of two adjacent iterations is less than the tolerance value (ε).

High-power converter (VSC) model
The active ( P AC ) and reactive ( Q AC ) power absorbed by HPC from the AC grid with ignored converter losses can be expressed as follows V AC is the amplitude of AC grid voltage, V C is converter output voltage, X HPC is equivalent reactance of converter, and δ converter modulation angle. However, the converter loss is ignored so that the active power is equal on the AC and DC side. Then the active power can be expressed as V DC is voltage on the DC side and M is a modulation index of the converter. Since the power supply on the DC side of the network is poor inactivity, the losses of the converter can be expressed by the current and resistance of the converter as follows where P C,Loss is power lost in the converter and R is the resistance offered by HPC.

Problem formation
Optimal allocation of REDG units in a hybrid AC/DC micro-grid system is to find the best location as well as the size of REDG units that gives minimum power loss as an objective function with the area required by REDG units as variables while satisfying various operating constraints. The objective function minimization of power loss is described as follows:

Constraints
Power balance constraint Inequality constraints N AC k=1k=1 J 2 K ,AC * R K is power loss in AC micro-grid system, N DC k=1 J 2 K ,DC * R K is power loss in DC micro-grid system, P REDG is power injected by REDG units in DC microgrid system, NB is number of branches in the system, N DG is number of REDG units connected, N is number of buses in the system, i is bus number, k is branch number, A PV , A FC are area required by PV and FC units, n PV , n FC are the number of PV and FC units connected, I PV , η PV , A PV are solar insolation (W/m 2 ), the efficiency of solar PV, the area required by PV unit (m 2 ) and V FC , A FC , J are the output voltage of FC, area required by FC unit (m 2 ), current density (A/m 2 ), P Slack is slack (substation) bus power, and P D is the real power load connected at ith bus. (12)

Butterfly optimization algorithm
Butterfly optimization is based on the ability of the butterflies to locate the source of fragrance accurately. They can also differentiate various fragrances and sense their intensities. In BO algorithm, butterflies are the searching agents. Fitness is correlated with the intensity of fragrance that can be generated by the butterfly. The movement of butterflies in search space will change its fitness. The sharing of information between butterflies is established through the propagation of fragrance. The searching ability of a butterfly depends on the sensing capability of the fragrance. This property will decide the movement of the butterfly towards a global search or local search (random). In BOA, the fragrance is formulated as a function of the physical intensity of stimulus as follows: where f is the perceived magnitude of the fragrance, i.e., fragrance receiving property by other butterflies, c is the sensory modality, I is the stimulus intensity, and a is the power exponent dependent on modality, which accounts the varying degree of absorption.
Most of the cases a&c ∈ [0, 1] . If a = 1, it means there is no absorption of fragrance, i.e., the amount of fragrance emitted by a particular butterfly is sensed in the same capacity by the other butterflies (fragrance propagation in an idealized environment). Thus, a butterfly emitting fragrance can be sensed from anywhere in the domain which in turn helps to reach the global optimum easily. On the other hand, if a = 0, it means that the fragrance emitted by any butterfly cannot be sensed by the other butterflies at all. Another important parameter c ∈ [0, ∞] determines the convergence speed. The values of a and c crucially affect the convergence speed of the algorithm. For the maximization problem, the intensity can be proportional to the objective function [29].
In BO algorithm, the characteristics of butterflies are idealized as follows: 1. Every butterfly is supposed to emit some fragrance which enables the butterflies to attract each other (propagation of information). 2. Every butterfly will move randomly or toward the best butterfly emitting more fragrance. 3. The stimulus intensity of a butterfly is affected or determined by the topography of the objective function.
The detailed steps for implementation of BO algorithm are as follows.
Step 1: Initialize algorithm parameters such as the number of agents N, the dimension of the problem d, the maximum number of iterations Iter max , probability switch P, power exponent PE and sensor modality SM.
Step 2: Generate initial random solution x j i where N is the number of agents and d is the number of decision variables,x j i represents the position of the ith agent in jth dimension generated randomly between the limits as x max,d and x min,d and rand() is a random number between 0 and 1.
where k = 1: N DCMG , N DCMG is number of DC micro-grids, L, A PV and A FC represents location of REDG, area of PV and area of FC, respectively.
Step 3: Evaluate the fitness (objective functions) of agents using Eq. 14. Record the gbest solution so far and set iteration count t as zero.
Step 4: Calculate the fragrance f N for each agent or butterfly using Eq. 20.
Step 5: Perform a global search using Eq. 25 if rand < probability P or local search using Eq. 26 if d > P.
where x d j (t) and x d k (t) are jth and kth butterflies from the solution space which belongs to the same swarm and r is a random number in [0, 1] and then Eq. 27 becomes a local random walk.
Step 6: Evaluate the fitness of each agent in the new population using Eq. 14.
Compare each new solution with the previous solution, if the new solution is better than the previous solution, record the gbest; otherwise, discard the new solution and preserve the previous solution as it is.
Step 8: Stopping criterion. If the maximum number of iterations has reached (iter max ), computation is terminated. Otherwise, Step 4 to Step 7 are repeated.
The implementation flowchart for the BO algorithm is illustrated in Fig. 1.

Results and discussions
Standard radial distribution systems are modified as hybrid AC/DC micro-grid systems for the purpose of the study. The buses in the AC zone will have real and reactive power loads and lines will have resistance and reactance. But in DC zone buses will have only real power loads and lines will have only resistance. Based on the system topology, it In order to demonstrate the effectiveness of the proposed approach for optimal allocation of REDG units in a hybrid AC/DC micro-grid system, it is applied to small-and medium-scale hypothetical 33-and 69-bus AC/DC hybrid distribution micro-grid systems. Tuned algorithm parameters of BO for the implementation are given in Table 1. The details of the REDG units used in this study are furnished in Table 2. All simulations are developed in MATLAB R2017a platform on Intel Core i5 2.7 GHz processor, 8 GB RAM. The optimal allocation of REDG units in a hybrid AC/DC micro-grid system is analyzed through two different cases.
Case 1: Optimal REDG units sizing at locations fixed by LRSF (Eq. 15 in [27]) using the BO algorithm without considering HPC losses.
Case 2: Simultaneous optimal allocation of REDG units using BO algorithm with and without considering HPC losses.

Optimal REDG unit allocation for 33-bus hybrid AC/DC micro-grid system
A hypothetical 33-bus hybrid AC/DC micro-grid system consists of three zones: Zone 1 is ACMG supported by substation and Zone 2, and 3 are DCMG supported by HPC and REDG units. The detailed 33-bus hybrid AC/DC micro-grid system is shown in Fig. 2. The maximum limit for area required by FC and PV REDG is taken as 3000 m 2 and 500 m 2 , respectively. The results obtained for the 33-bus hybrid AC/DC micro-grid system for two cases are given in Table 3. From Table 3, it is observed that the result obtained by BO is better than other results. In Case 1, i.e., without considering HPC losses, the lowest total area required by REDG units is 4537.2 m 2 with minimum system power loss   A detailed comparison of the area required by REDG units versus total power losses of the system by various optimization algorithms is presented in Fig. 3. From Fig. 4, it is clear that in both cases the BO algorithm performed better than other algorithms. Another significant benefit of the REDG unit's optimal allocation in hybrid AC/DC micro-grid is voltage profile improvement which is shown in Fig. 4. From Fig. 6, it is observed that the bus voltage profile has been improved in all three zones of the system for both cases with and without considering HPC losses. Convergence characteristics of the BO algorithm for the minimization of the desired objective function for different cases are given in Fig. 5. From Fig. 6, it is understood that the BO algorithm reached the final solution with good convergence i.e., 13th iteration for Case 1 without HPC losses, at 57th iteration for Case 2 without HPC losses and at 36th iteration for Case 2 with HPC losses. It is identified that the energy supplied through substation is 92514.38 kWh, the amount of coal consumed is 64.76 ton, and CO 2 emission into the atmosphere is 55.14 ton without the integration of REDG units into the system. Energy flow results of the system under REDG units standalone mode of operation are presented in Table 4. From Table 4, it is clear that the proposed approach can improve the power loss reduction percentage, i.e., with HNMCS its value is 67% and with BO is 80%. Energy-saving from the substation is improved from 65 to 67%, coal consumption is reduced from 22 to 19 tons that intern reduced the CO 2 emission from 18.9 tons to 16 tons. And it is also observed that integrated REDG units under the standalone mode of operation in Zone 2 and 3 (DC grid) are producing sufficient   energy to cater to the needs of AC grid without any shortfall and excess energy is useful to support the AC grid.

Optimal REDG unit allocation for 69-bus hybrid AC/DC micro-grid system
69-bus hybrid AC/DC micro-grid system consists of three zones: Zone 1 is ACMG supported by substation and Zone 2, Zone3 and Zone 4 are DCMG supported by HPC and REDG units. The detailed 69-bus hybrid AC/DC micro-grid system is shown in Fig. 6. The maximum limit for area required by FC and PV REDG is taken as 500 m 2 and 800 m 2 , respectively. Numerical outcomes obtained by case-wise based on BO algorithm for 69-b*us hybrid AC/DC micro-grid system are in Table 5. And a comparison of outcomes based on various algorithms is furnished in Table 6. From Table 6, it is observed in Case 1, i.e., without considering HPC losses the lowest total area required by REDG units is 1361.84 m 2 with minimum system power loss of 11.50 kW and the total power loss reduction is 92%. The optimal power required by PV is 15 kW and FC is 1017 kW at 8th bus with 8 PV panels and 40 FC arrays. Similarly, the optimal power required by PV is 8 kW and FC is 277 kW at 28th bus with 4 PV panels and 11 FC arrays and power required by PV is 3 kW and FC is 1656 kW at 61st bus with 2 PV panels and 66 FC arrays. However, the results obtained by BO and HNMCS are almost the same. But the results obtained by BO in Case 2 are surprising. The total area required by REDG units is 1344.2 m 2 with minimum system power loss of 4.65 kW and the total power loss reduction is 98.2%. Optimal locations obtained for the integration of REDG units are 18, 61 and 54 buses. The optimal power produced by PV and FC at these locations is 3.8 kW, 8.03 kW, 24.21 kW, and 470.8 kW, 163.8 kW, and 427.6 kW, respectively. It is observed that the reduction in system power losses is significant in Case 2. The reason is in Case 1, the algorithm does not have the flexibility to choose the optimal locations because they are fixed based on LRSF. At those fixed locations algorithm has to find the best suitable size of REDG units. But, in Case 2 the algorithm has the freedom to choose optimal locations as well as optimal sizes simultaneously. It is worth noting point that the proper size of REDG units at proper locations always guarantees the better performance of the system. It is also observed that the base case power loss is 227.4 kW (including HPC losses). The total area required and power produced by REDG units is 1301 m 2 and 2486.2 kW, respectively. The optimal locations are 19, 61 and 54 buses. In this case, the system power loss is increased from 144 to 227 kW. The reason is power loss due to HPC's will alter the power flow in the respective feeder, and as a result, the change in branch currents may increase the feeder losses. A detailed comparison of the area required by REDG units versus total power losses of the system by various optimization algorithms is presented in Fig. 7. From Fig. 8, it is clear that in both cases the BO algorithm performed better than other algorithms. Another significant benefit of REDG unit's optimal allocation in hybrid AC/DC microgrid is voltage profile improvement which is shown in Fig. 8. From Fig. 9, it is observed that the bus voltage profile has been improved in all three zones of the system for both cases with and without considering HPC losses. Convergence characteristics of the BO algorithm for the minimization of the desired objective function for different cases are given in Fig. 9. From Fig. 9, it is understood that the BO algorithm reached the final solution with smooth convergence. It is identified that the energy supplied through substation is 25848 kWh, the amount of coal consumed is 12.07 ton, and CO 2 emission into the atmosphere is 10.2 ton without the integration of REDG units into the system. Energy flow results of the system with REDG unit's standalone mode of operation are presented in Table 7. From Table 7, it is clear that the proposed approach can improve the power loss reduction percentage i.e., with HNMCS its value is 90.1% and with BO is 97.5%. Energy-saving from the substation is improved from 76 to 80%, coal consumption is reduced from 15.5 to 12 ton that intern reduced the CO 2 emission from 13.2 to 10.2 ton. And it is also observed that integrated REDG units under the standalone mode of operation in Zone 4 (DC grid) are producing excess, i.e., 8601 kWh, energy which is not sufficient to cater the needs of Zone 1 (AC grid), Zones 2 and 3 (DC Grid) and hence it requires load shedding or additional REDG support from other DC grids or through the substation to satisfy the demand. Since BOA is a heuristic search method, its outcome may have certain randomness. In this work, it has been handled carefully by proper selection of decision variable limits. So, BO algorithm is tested by running 100 times for each case. Further, the computational complexity of the proposed BOA is analyzed by popular statistical methods such as Friedman and Quade test. The obtained ranks by BOA are compared with existing results available in the literature and also furnished in Table 8. From Table 8, it is observed that solution-quality and robustness-wise BOA is almost closer to HNMCS. However, HNMCS proves to be better than other algorithms in the race due to its hybrid nature.

Conclusions
In this paper, an efficient approach was proposed for optimal integration of REDG units in hybrid AC/DC distribution micro-grids to maximize the technical benefits of the system with minimum operational area required by REDG units. Optimal locations and sizes for REDG units are determined simultaneously using a butterfly optimizer which is proved as a better approach by obtained results. It is identified that additional REDG units are required in DC zones of the 69-bus system to produce the energy to satisfy the present demand. The proposed approach is an efficient technical tool for achieving better results in an optimal allocation of REDG units in hybrid AC/DC distribution microgrids. The future scope of this work is to include economic and technical benefits of