# Improvement of the Voltage Profile and Loss Reduction in Distribution Network Using Moth Flame Algorithm: Wolaita Sodo, Ethiopia

Authors: S. Balakumar; Akililu Getahun; Samuel Kefale; K. Ramash Kumar
Journal: Journal of Electrical and Computer Engineering
 (2021)

---

## Abstract

Voltage stability and line losses are inevitable issues even in modern power systems. There are several techniques that emerged to solve problems in the power system to provide quality and uninterrupted supply to customers. The algorithms used in this paper to determine the appropriate location and size of the Static Var Compensator (SVC) in the Distribution Network (DN) are Moth Flame Optimization (MFO) and Particle Swarm Optimization (PSO). The objective function is defined to minimize voltage deviation and power loss. The burning problem of voltage stability improvement current scenario is because of a rise in electricity demands in all sectors. Paramount duties of power engineers are to keep the system stable and maintain voltage magnitude constant even during peak hours. The results were checked with the aid of MATLAB onWolaita Sodo radial distribution of 34 bus data networks. The potential use of SVC is key to solve distribution system power quality issues and estimating the advantage of the installation. The results obtained from the test system were compared with PSO results. This comparison was done to know the computational time of proposed techniques. The performance of the MFA based SVC was superior in distribution system and highlighted the importance of device.


---

## Body

## 1. Introduction

In modern times, distribution network structure is going through massive transformation. Due to the introduction of new technologies, distribution network structure is becoming more complex, and there is an ever-increasing demand for load. During voltage breakdown, voltage instability occurs, resulting in loss of voltage for a significant part of the system [1]. In order to increase voltage stability margin, FACTS controllers should be installed at the appropriate location to incorporate reactive power. FACTS devices comprise series, hybrid, and shunt devices. When dealing with a specific system voltage instability problem, use the correct component to restore system voltage stability [2]. Network reconfiguration (NR) is known as the process of modifying the network topology by reconfiguring the distribution network (DN).Increasing the power flow tends to reduce the percentage power loss. Because of voltage profile improvements, smooth demand patterns, and improved network reliability, network reconfiguration offers major benefits [3]. The Particle Swarm Optimization (PSO) technique was used to achieve minimization active power losses and total function cost. This involves including the generation costs and the UPFC investment costs, which are deployed by installing the UPFC in the power system network to maximize benefits [4]. One novel optimization, inspired by nature, is presented in this study. The navigational strategy is employed by moths in the universe. The algorithm is known as the moth flame optimization (MFO)—a creative problem that presents engineering design and constrained optimization [5]. Many of the researchers have proved that MFO is one of the suitable algorithms for engineering problems.The authors like Yu et al. had solved power system issues and some of the engineering problems with MFO effectively. The MFO technique was successfully handled in the optimal placements of the DG [6].There are keen researchers using MFA for the engineering problems like power engineering, computer engineering, and mechanical engineering. In this paper, the MFA algorithm was given quick results with less computational time than PSO. The moth flame optimization approach has also been used in medical and related sectors. The following are some of the studies of some academics. For the automatic detection of mitosis in breast cancer tissue pictures, Sayed and Hassanien [7] utilized neutrophils to optimize the moth flame population. Mingjing et al. [8] presented a new learning method for a nuclear limit learning machine using a chaotic moth flame optimization algorithm for medical diagnostics. In this work, a two-stage optimization strategy for optimal distributed generation (DG) planning with energy storage integration is provided.The first stage uses the well-known Loss Sensitivity Factor (LSF) approach to determine the installation locations and initial capacity of DGs, and the second stage determines the optimal installation capacities of DGs to maximize investment benefits and system voltage stability and minimize line losses. The proposed two-stage optimization method outperforms current state-of-the-art DG planning approaches like NSGA-II, MOPSO, and MOHS, implying that our method achieves better outcomes across the range, including higher investment benefits, improved voltage stability, and lower line losses [9]. The Kenyan Power Network has seen significant increases in load demand in recent years. Although growing load demand has been partially met by additional transmission and generating capacity, due to economic, environmental, and geographic limits, the rate of growth has not kept pace with the rate of increase in load demand [10]. The goal of this research is to combine a microgrid system with a Static Synchronous Compensator (STATCOM) controller in order to ensure a higher power flow with a better voltage profile and less power loss. The microgrid and STATCOM controller connected to the grid are modelled using MATLAB/PSAT. The STATCOM is also controlled by Proportional Integral (PI) control and Fuzzy Logic Controller (FLC). By connecting to the main distribution system via a typical IEEE 30-bus connection, the usefulness of STATCOM with microgrid integration is assessed. Finally, it was discovered that STATCOM improves the voltage profile and reduces power loss by increasing the capacity of the distribution line [11].The proposal’s innovation is that it uses IPSO to optimize the model’s parameters with synchro phasors, which enhances the prediction performance of the ELM-based TSP model. Finally, the correctness and validity of the provided approach are confirmed based on test results on both an IEEE 39-bus system and a large-scale real power system. Machine learning has shown to be effective in resolving online TSA issues. The existing PRTSA methodologies, on the other hand, are unable to fulfill the demands of big data administration and usage [12].In this paper, various FACTS devices were tested as fault current limiters. Several FACTS devices are being explored for their performance, including the SSSC, STATCOM, and UPFC. SSSC has little impact on fault current and voltage regulation and instead focuses solely on reactive power flow. UPFC and STATCOM, on the other hand, can reduce fault current while also correcting voltage and regulating current. STATCOM and UPFC absorb reactive power from the system in such a way that fault currents are considerably reduced. UPFC outperforms STATCOM in terms of stability, transients, and voltage control. As a result, the crucial clearing time for the system will be extended [13].


## 2. Statement of the Problem

In this paper, Wolaita Sodo substation and its surroundings were taken as case study area. It consisted of several feeders, transformers, lines, and buses. The design of distribution system looks exactly radial which has 34 buses. During the data collection, researchers have noticed several problems associated with substation and around 300 km radius. The researchers aimed to registered current status of the substation by analysing the collected data. The values of voltage variations, power factor, harmonics, etc. had been measured by using MICROVIP3 PLUS. All the collected data were interpreted in the graphical structure. The overall data of the substation were collected and encapsulated in this work for future reference. The following graphs depicted various parameters values like voltage variation and harmonics. Let us view them one by one.


### 2.1. Line Voltage Variation and Power Factor in Distribution Network

This paper’s study area was Wolaita Sodo distribution system which is located in the south part of Ethiopia. The geographical location is southern side of Addis Ababa, the capital city of Ethiopia; it has a latitude and longitude of 6°54′N 37°45′E with an elevation between 1̦600 and 2̦100 meters above sea level. For this study, voltage measurement is taken from Wolaita Sodo substation’s all outgoing feeders control panel. Since the substation supplys 15 kV to the distribution line, the voltage variation is occurring in the line. In Table1, the distribution line voltage (kV) variation for peak time 8:00–21:00 is shown. This observed value leads researchers to focus on voltage profile of the distribution network.Table 1
Distribution line voltage (kV) variation for peak time 8:00–21:00.
Time (s)Gesuba line 1 (kV)Sodo line 2 (kV)Bodit line 4 (kV)Bilate line 5 (kV)33 kV Abala line 2 (kV)33 kV Bedesa line 3 (kV)8:0014.814.614.814.832.732.89:0014.714.514.814.832.832.710:001514.714.714.732.832.711:0014.714.21515333312:001514.68151533.13313:001514.21515333314:001514.651515.132.832.715:0014.814.614.81532.732.716:0015.114.715.114.832.832.817:001514.671515.13333.118:001514.761515333319:0014.814.71515333320:001514.514.914.932.933In Figure1, the voltage variation of Sodo line 2 from the substation control panel is shown. The data was measured in the intervals of one hour. There are drastic changes in the load around Wolaita substation due to residential and industrial enlargement. If observed the obtained graph Figure 1 it was clearly reflecting the variation in the voltage from morning 8 a.m to 8 p.m. the lowest voltage point registered in the graph on particular day was around 14.2 kV. The power factor is a vital issue in the distribution level. The low power factor may lead to line losses and cause deteriorate the performance of the system. This power factor issues can be sort out by external devices like FACTS family because it is local phenomena. In this analysis, researchers were tried to adapt the optimal device to improve voltage profile in distribution system. The power factor for the distribution line of substation observed from the control panel is shown in Table 2.Figure 1
Voltage variation for Sodo line 2 from the substation control panel.Table 2
Power factor for the distribution lines of substation.
Time (s)Gesuba line 1Sodo line 2Bodit line 4Bilate line 533 kV Abala line 233 kV Bedesa line 38:000.90.80.90.90.80.99:000.80.870.80.80.90.810:000.90.850.90.90.80.911:000.90.880.80.80.80.812:000.890.860.90.90.90.913:000.90.870.90.90.90.914:000.90.840.80.80.80.915:000.90.80.840.840.840.816:000.940.870.870.870.870.8417:000.90.890.960.960.960.8718:000.890.890.950.940.950.9119:000.90.80.90.80.940.9620:000.90.810.80.910.90.91In the distribution line, the voltage drop occurs due to long-distance coverage to achieve the utility’s coverage. There is universal truth that power flow always flows from higherδ to lower δ; likewise, reactive power flows from higher voltage to lower voltage—in general, sending end to receiving end—but this will be failed sometimes in case of reactive power due to uneven voltage profile at distribution system. Such common issues could be resolved with aid of FACTS devices. In this work, SVC plays a superior role in enhancing voltage throughout Wolaita Sodo distribution network. There are 34 buses dedicated to distributing power to the consumers. All the buses voltage must be near to knee point (0.95 to 1.05P.u). The moth flame algorithm was used for placing and sizing of the SVC in the network.The main emphasis of the work was voltage profile improvement. The researchers were focused on and investigated reactive power support for the system and worked on maintaining constant voltage with the help of optimization-based SVC. Figure2 shows the photograph while taking measurements by using MICROVIP3 PLUS.Figure 2
Photo during measurements.


## 2.1. Line Voltage Variation and Power Factor in Distribution Network

This paper’s study area was Wolaita Sodo distribution system which is located in the south part of Ethiopia. The geographical location is southern side of Addis Ababa, the capital city of Ethiopia; it has a latitude and longitude of 6°54′N 37°45′E with an elevation between 1̦600 and 2̦100 meters above sea level. For this study, voltage measurement is taken from Wolaita Sodo substation’s all outgoing feeders control panel. Since the substation supplys 15 kV to the distribution line, the voltage variation is occurring in the line. In Table1, the distribution line voltage (kV) variation for peak time 8:00–21:00 is shown. This observed value leads researchers to focus on voltage profile of the distribution network.Table 1
Distribution line voltage (kV) variation for peak time 8:00–21:00.
Time (s)Gesuba line 1 (kV)Sodo line 2 (kV)Bodit line 4 (kV)Bilate line 5 (kV)33 kV Abala line 2 (kV)33 kV Bedesa line 3 (kV)8:0014.814.614.814.832.732.89:0014.714.514.814.832.832.710:001514.714.714.732.832.711:0014.714.21515333312:001514.68151533.13313:001514.21515333314:001514.651515.132.832.715:0014.814.614.81532.732.716:0015.114.715.114.832.832.817:001514.671515.13333.118:001514.761515333319:0014.814.71515333320:001514.514.914.932.933In Figure1, the voltage variation of Sodo line 2 from the substation control panel is shown. The data was measured in the intervals of one hour. There are drastic changes in the load around Wolaita substation due to residential and industrial enlargement. If observed the obtained graph Figure 1 it was clearly reflecting the variation in the voltage from morning 8 a.m to 8 p.m. the lowest voltage point registered in the graph on particular day was around 14.2 kV. The power factor is a vital issue in the distribution level. The low power factor may lead to line losses and cause deteriorate the performance of the system. This power factor issues can be sort out by external devices like FACTS family because it is local phenomena. In this analysis, researchers were tried to adapt the optimal device to improve voltage profile in distribution system. The power factor for the distribution line of substation observed from the control panel is shown in Table 2.Figure 1
Voltage variation for Sodo line 2 from the substation control panel.Table 2
Power factor for the distribution lines of substation.
Time (s)Gesuba line 1Sodo line 2Bodit line 4Bilate line 533 kV Abala line 233 kV Bedesa line 38:000.90.80.90.90.80.99:000.80.870.80.80.90.810:000.90.850.90.90.80.911:000.90.880.80.80.80.812:000.890.860.90.90.90.913:000.90.870.90.90.90.914:000.90.840.80.80.80.915:000.90.80.840.840.840.816:000.940.870.870.870.870.8417:000.90.890.960.960.960.8718:000.890.890.950.940.950.9119:000.90.80.90.80.940.9620:000.90.810.80.910.90.91In the distribution line, the voltage drop occurs due to long-distance coverage to achieve the utility’s coverage. There is universal truth that power flow always flows from higherδ to lower δ; likewise, reactive power flows from higher voltage to lower voltage—in general, sending end to receiving end—but this will be failed sometimes in case of reactive power due to uneven voltage profile at distribution system. Such common issues could be resolved with aid of FACTS devices. In this work, SVC plays a superior role in enhancing voltage throughout Wolaita Sodo distribution network. There are 34 buses dedicated to distributing power to the consumers. All the buses voltage must be near to knee point (0.95 to 1.05P.u). The moth flame algorithm was used for placing and sizing of the SVC in the network.The main emphasis of the work was voltage profile improvement. The researchers were focused on and investigated reactive power support for the system and worked on maintaining constant voltage with the help of optimization-based SVC. Figure2 shows the photograph while taking measurements by using MICROVIP3 PLUS.Figure 2
Photo during measurements.


## 3. Research Gap of the Work

Voltage profile improvement and loss minimization were successfully obtained using MFA in the distribution network. Voltage deviation and line losses are indispensable problems even in modern power system. When the undesirable symptoms recurred, the power engineers diagnosed something different to tackle issues. Nowadays, maintaining the distribution system under normal condition is a big challenge due to all countries, particularly for developing countries. All researchers begin the quest for a reliable artificial intelligence method to settle the existing problems. In such a manner, in this work, searching for different and feasible algorithms has intensified. Finally, with deep analysis and literature review, MFA has been selected for optimization techniques. It was tested on Wolaita Sodo 34-bus system successfully and to show greatness of the work, load of the tested system had changed randomly to register changing of the location and size. The proposed algorithm could work on any kind of distribution system; depending on real and reactive power demands, it will optimize the place and size of the FACTS devices. Therefore, research gap of the work is clearly fulfilled and executed in the respective sections of the paper.


## 4. Moth Flame Optimization

Seyedali Mirjalili introduces the Moth-Flame Optimizer in 2015. MFO is an algorithm based upon population and represents the moths in the set of matrix:(1)M=m11m11…m1d⋮⋮⋮⋮m21m22…m2d⋮⋮⋮⋮mn1mn2…mnd.The numbersn and d are the number of moths (dimension) variables. In this paper, assume that the corresponding fitness values are stored in the following way for all moths [5]: where n is the moth count,(2)OM=om1⋮om2omn.Let us remember that the fitness value is the fitness return value (target) for each moth function. The fitness function and the output of the fitness function are passed to each moth position vector (first row in matrixM, for example) assigned as a fitness function to the corresponding moth (for example, OM Matrix OM). It is the moth set in a matrix. Moth flame optimization algorithm is conceived as transverse orientation by moth motion in nature. Moth flames could go right along the long distances of the night by keeping a fixed angle to their moon, but because of the artificial light they may move in a deadly spiral path. The moths are candidates for solutions in the proposed method, whereas the dependent variables or control variables are their position in space. The biggest distinction between the moths and the flames is that the moths are considered to be search agents migrating around the search area and that the flames portrayed by moth are the ultimate solution acquired. The flames are therefore the flags for moths; every moth moves around the flames and updates its position when a stronger decision is obtained, where n is the moth number and d the variable numbers (dimensions). Across all moths, researchers also assume that there is an array for the storing corresponding fitness values saved in array (M). We understand that for every moth the value of fitness is the fitness return value (objective) function. Each fitness function is transferred to the position vector (e.g., first row in matrix M) of each moth and the fitness function output is delegated as the fitness function of the respective moth (for example, OMI in matrix and OM). Flames are also the main components of the proposed algorithm. We consider a moth matrix-like matrix.Researchers understand exactly that theM and F arrays dimensions are the same. For flames, authors also believe that the fitness values are saved in an array:(3)OF=Of1Of2⋮Ofn,where n is the moth number. Moths and flames are both solutions. The difference is how each iteration is handled and updated. The moths are original search agents travelling around the search area, while flames are the best position to get moths up to now. Each moth is therefore looking for and updating a flame if a better solution is found. A moth has not ever lost its optimal answer with this technique.


## 5. Moth Flame Algorithm for Placing and Sizing SVC

### 5.1. Step-by-Step Algorithm of MFO

Step1: Read the information on the line and bus.Step2: Compute power loss using a backward/forward sweeping load flow method, without an SVC connection. Deal with the power loss as the standard model.Step3: Initialize moth population variables upper and lower limits, highest iterations and browse variable dimensions.Step4: The random population (array) of the moth places continues to generate an initial array by equating appropriate with an array of the SVC sizes within defined by default top and bottom borders. Set the number of iteration (ki = Zero).Step5: Inspect whether the limitations are fulfilled or not in the moth positions.Step6: If limitations have not been violated, compute the value of objective function (Complete loss of power) with load flow.Step7: Save as best flame the minimum value of the target function and the related to moth position as optimal flame position.Step8: Using equation (A)(i) = (Li.ebt.cos (2πt) + Bj) Review the moth position to enable moths to converge to an optimal solution.Step9: If exceed the maximum threshold of iteration, go to step5.Step10: If not, update theK =K + 1 iteration count and return to step5.Step11: Break the best solution for the problem of optimization. The finest possible solution is to optimize SVC’s location and size. The respective objective function value which embodies the lowest possible real power loss.


### 5.2. Flowchart of MFA

Figure3 shows the exact function of the MFA. The coding of the algorithm was developed in the m-script of the MATLAB platform. In this section, crucial execution of the program was explained step by step. The results of the MFA had been compared with PSO. The system data could be read by the function [d, V_base] = data (). The maximum number of the bus and line could be taken by the command nbus = max (max (d (: , 1). The results of load flow data were input to placing and sizing of SVC. In this optimization algorithm, maximum number of search agent is 30. The maximum number of iterations is 100. N is denoted as number of search agent. The SVC rating is in kVAR. The upper limit and lower limit were decided based on total reactive power demand of the system. The reactive power ranges from −1kVAR to +1kVAR. The dimension of the problem was taken as 1. Then initialize the positions of the moths using (4) and search randomly to fix best fitness.(4)Moth_pos:,1=SVCmin+SVCmax−SVCmin∗randN,1.Figure 3
Flowchart of MFA.The number of the flames in the work was decided by the following:(5)Flame_no=roundN−Iter∗N−1Max_iteration.It is mandatory to check whether the moths are out of search space or not. If the moths are out of search space, it should be brought into search boundaries. The fitness of moths was calculated by MATLAB function:(6)Moth_fitness1,i=Obj_funcMoth_posi,:.The next steps after the fitness were sorting the first population of the moths based on best fitness. After the sorting of the flame, update the best position of the flames. The 8th step in the algorithm gives the absolute value of the flame distance. The algorithm had converged after flame scored best value. The converged speed of the algorithm was faster than PSO in this work. In this paper, MFA dominated over the PSO for the placing and sizing the SVC in the test system.


## 5.1. Step-by-Step Algorithm of MFO

Step1: Read the information on the line and bus.Step2: Compute power loss using a backward/forward sweeping load flow method, without an SVC connection. Deal with the power loss as the standard model.Step3: Initialize moth population variables upper and lower limits, highest iterations and browse variable dimensions.Step4: The random population (array) of the moth places continues to generate an initial array by equating appropriate with an array of the SVC sizes within defined by default top and bottom borders. Set the number of iteration (ki = Zero).Step5: Inspect whether the limitations are fulfilled or not in the moth positions.Step6: If limitations have not been violated, compute the value of objective function (Complete loss of power) with load flow.Step7: Save as best flame the minimum value of the target function and the related to moth position as optimal flame position.Step8: Using equation (A)(i) = (Li.ebt.cos (2πt) + Bj) Review the moth position to enable moths to converge to an optimal solution.Step9: If exceed the maximum threshold of iteration, go to step5.Step10: If not, update theK =K + 1 iteration count and return to step5.Step11: Break the best solution for the problem of optimization. The finest possible solution is to optimize SVC’s location and size. The respective objective function value which embodies the lowest possible real power loss.


## 5.2. Flowchart of MFA

Figure3 shows the exact function of the MFA. The coding of the algorithm was developed in the m-script of the MATLAB platform. In this section, crucial execution of the program was explained step by step. The results of the MFA had been compared with PSO. The system data could be read by the function [d, V_base] = data (). The maximum number of the bus and line could be taken by the command nbus = max (max (d (: , 1). The results of load flow data were input to placing and sizing of SVC. In this optimization algorithm, maximum number of search agent is 30. The maximum number of iterations is 100. N is denoted as number of search agent. The SVC rating is in kVAR. The upper limit and lower limit were decided based on total reactive power demand of the system. The reactive power ranges from −1kVAR to +1kVAR. The dimension of the problem was taken as 1. Then initialize the positions of the moths using (4) and search randomly to fix best fitness.(4)Moth_pos:,1=SVCmin+SVCmax−SVCmin∗randN,1.Figure 3
Flowchart of MFA.The number of the flames in the work was decided by the following:(5)Flame_no=roundN−Iter∗N−1Max_iteration.It is mandatory to check whether the moths are out of search space or not. If the moths are out of search space, it should be brought into search boundaries. The fitness of moths was calculated by MATLAB function:(6)Moth_fitness1,i=Obj_funcMoth_posi,:.The next steps after the fitness were sorting the first population of the moths based on best fitness. After the sorting of the flame, update the best position of the flames. The 8th step in the algorithm gives the absolute value of the flame distance. The algorithm had converged after flame scored best value. The converged speed of the algorithm was faster than PSO in this work. In this paper, MFA dominated over the PSO for the placing and sizing the SVC in the test system.


## 6. Results and Discussion

### 6.1. Scenario 1

In this section, two different optimization techniques were used for sizing and placing the SVC in the Wolaita Sodo network. To ensure the suitable optimization techniques for the distributed network, different scenarios were formed and the effectiveness of the techniques was observed. The computational time of the algorithm was also taken as a crucial parameter to fix the superiority in the Wolaita Sodo distribution network. The suggested reasoning was tested on 34 radial bus systems and the results obtained were improved with MFA and PSO optimizations. The test system had a demand of 3.72 MW and 2.3 MVAR for real power and reactive capacity, respectively. The PSO technique is carried out on a radial 34-bus network, as per the standard steps. This MATLAB simulation program was developed and run on Intel Pentium dual core personal computer with 2 GB RAM. The final result obtained from the MATLAB software is shown in Table3. It is evident that the PSO approach shows its advancement and benefits over the existing method in terms of line loss decrement and computational time.Table 3
Comparison results show location and rating of SVC and power loss reduction in 34-bus system.
TechniquesMFAPSOOptimum placing of deviceBus 7Bus 11Optimal size of SVC23001352 kVARActive power loss (kW)Without SVC210.99369.57 kWWith SVC168.57270.01 kW% reduction in real power losses42.42 kW27.6Reactive power loss (kVAR)Without SVC143.13356.07 kVARWith SVC126.26260.12 kVAR% reduction in reactive power losses16.87 kVAR27.4 kVARComputational time in sec0.65 sec1.145873Figure4 clearly reiterates voltage deviation in the system while absence of SVC and voltage magnitude was restored in its optimal range with presence of SVC. In Figure 4, red colour line shows the results with SVC. Almost all 34 buses’ voltage was maintained constant with reactive power support. The convergence line of PSO drawn on 34-bus system is also depicted in Figure 5. The iteration-1 starts objective function of 298.7986 and objective function is converged at iteration-100 value of 269.32 and elapsed time is 1.145873 seconds. Before SVC placement on optimal location, real power loss (kW) is 369.57 and reactive power loss (kVAR) is 356.07. After optimal location and optimal sizing of SVC/, the real power and reactive power loss were reduced into 270.01 kW and 260.12 kVAR, respectively.Figure 4
Voltage profile with and without SVC using PSO.Figure 5
Convergence property of PSO algorithm.


### 6.2. Scenario 2

The first scenario was tested successfully with PSO algorithm. In the second scenario, recent optimization technique, moth flame algorithm, was introduced for sizing and placing the SVC. This technique coding was developed in MATLAB and run on Intel Pentium dual core like PSO code. In Figure6, results of voltage profile are constant in the entire distribution network. The range of the distribution system voltage started from 0.91 to 1.05 pu. The range in Figure 4 was accepted and matched with IEEE standard. This is one of the evidences that line losses would be minimized because of the constant voltage maintained in the system. The bus 7 was selected as weakest bus in the next for moment. The optimal size of the SVC was designed with MFA and placed at bus 7 for reactive power support. During the operation, bus 7 suffered more than the rest of the buses to maintain constant voltage at knee point.Figure 6
Convergence property of MFO.The optimal size of the SVC had been decided value of 2300 kVAR. The active power loss in kW without SVC was 210.99 kW and with SVC 168.57 kW. The percentage of the losses reduction in the system was 42.43 kW. The reactive power losses without and with SVC were 143.13 and 126.26 kVAR, respectively. The percentage of reactive power loss observed was exactly 16.87 kVAR. The overall results of the work are showcased in Table3. The comparisons were made between MFA and PSO in the all level performance. Figure 7 shows voltage profile of the 34 buses and all the bus voltage was maintained constant around 1pu after inclusion of SVC with MFA.Figure 7
Voltage profile of 34-bus system with MFA.The MFO technique had taken slight edge in terms of computational time and was given proper placing and sizing to fix the problem in the system. The power system has been more dynamical operation and uncertainty may occur any time. To prevent the faults, stability issues, compensation, and protection of the devices, the computerized monitoring system must act accordingly as fast as possible to take decision over the uncertainties. In that sense, MFA was shown its superiority in the case of improving voltage profile of the system. The big contribution of the work was helping to choose the rating of SVC during needed conditions. The operation of these proposed algorithms, MFA and PSO, would deliver exact sizing and placing of the SVC in the system. The detailed information of the collected and calculated data of Wolaita Sodo substation is recorded in Tables4 to 7. The data used in the work was measured from the case study substation. To implement the algorithm-based SVC in the system proper data from concern substation is mandatory. The MICROVIP3 PLUS device was used to measure harmonics and voltage variations in different places of the substations.Table 4
Bus data for radial distribution system.
Branch no.Nominal load at receiving busActive load power (kW)Reactive power (kVAR)S/S 100224.414.5376.837.3474.946.2540.943.2635.913.5755.91088.953.39970.142.81078.520.81143.225.71243.6123.61370.320.91478.29.561528.610.51614.87.231724.511.21825.38.621956.124.72036.114.32148.222.52220.515.62333.617.6243015.52570.220.92655.9102738.237.12843.6123.62941.820.73014085.33137.922.93225283390.836.2342110.2Table 5
Transformer data of the distribution system.
Branch no.Transformer (kVA)No-load losses (W)Full load losses (W)Z (%)S/S 10002315110036305.231002802000645035012304531582041006620089027005.2720090030004.581002802050593159503650410100520180041163045625204.5122001002200413630130065004.5141005601820415315730365041620033028605.51710044017604.5183151050420041910050018506201004401700421200120260042220060027004.5231006001550424251103504253151300430042650220100042720060027004.52810032017004.529315110037805.3301004001550431100272150543220033037004.133350220100434252106954.5Table 6
Line impedance and admittance per unit on a distribution line.
Branch no.Sending endReceiving endImpedance (Z pu)R (pu)X (pu)|y|<Y1120.012196 + 0.05756i0.01220.0575616.9959782230.02591 + 0.03344i0.025910.0334423.6389−52.23340.11651 + 0.07165i0.116510.071657.3111−31.64450.03048 + 0.05504i0.030480.0550415.8942−61.025560.16952 + 0.17598i0.169520.175984.09252−46.076670.15212 + 0.07592i0.152120.075925.88191−26.527780.00204 + 0.07014i0.002040.0701414.2512−88.338890.26552 + 0.17672i0.265520.176723.13526−33.6499100.01855 + 0.00397i0.018550.0039752.7146−121010110.09104 + 0.077876i0.091040.077888.34698−40.51111120.16519 + 0.0776i0.165190.07765.47919−25.161212130.02263 + 0.072987i0.022630.0729913.0865−72.771313140.002089 + 0.01968i0.002090.0196850.5291−83.941414150.091004 + 0.088996i0.0910.0897.85625−44.361515160.082178 + 0.02799i0.082180.0279911.5189−18.81616170.169369 + 0.070098i0.169370.07015.45548−22.481717180.15014 + 0.085191i0.150140.085195.79289−29.57182190.04984 + 0.090347i0.049840.090359.69158−61.111919200.096827 + 0.08694i0.096830.086947.68457−41.922020210.09178 + 0.135156i0.091780.135166.12097−55.82214220.09178 + 0.136356i0.091780.136366.08394−56.052222230.0251 + 0.07067i0.02510.0706713.3342−70.452323240.08867 + 0.10267i0.088670.102677.3714−49.182424250.02511 + 0.07067i0.025110.0706713.3336−70.44257260.07015 + 0.07015i0.070150.0701510.0799−452626270.173004 + 0.175156i0.1730.175164.06189−45.352727280.08999 + 0.010676i0.089990.0106811.035−6.772828290.071867 + 0.062764i0.071870.0627610.480441.132911300.082142 + 0.09494i0.082140.094947.96543−49.133230310.0929067 + 0.08735i0.092910.087357.84183−43.233114320.08735 + 0.08589i0.087350.085898.16303−44.513232330.088996 + 0.081876i0.0890.081888.2692742.613333340.152498 + 0.090676i0.15250.090685.63635−30.73Table 7
Each line real and reactive power in the system run in backward/forward MATLAB.
LinePower at busFromTokWkVAr12423.5981999.437231342.361732.404341725.7871061.46545843.74971523.622561188.4481233.752671328.319662.91777828.9178996.463389962.906640.9058910853.0899182.71741011728.5451623.1971112559.4409262.65181213170.2558549.0849131445.2561426.36981415111.1573108.7036151691.66231.2216161769.125728.6095171822.945513.019321968.2395123.7001192070.263.0337202125.536537.6061422109.423162.5714222335.884100.9834232460.74670.3374242514.39340.5043726147.9189193.52032627112.8311114.234272891.46810.8508282934.057529.7441130130.5708150.9173303131.660529.76731432107.5726105.7773323384.373677.6234333420.031611.9109The collected 34-bus data is listed in Table4. In Wolaita Sodo region, there are 34 distribution transformers under working conditions with different rating. All the transformer data are placed in Table 5. The line data of the system is catalogued in Table 6 and real and reactive power data of each bus during load flow execution is showcased in Table 7. The overall data given in Tables 4–7 could be useful for future enhancement of the realistic work in the same substation.


## 6.1. Scenario 1

In this section, two different optimization techniques were used for sizing and placing the SVC in the Wolaita Sodo network. To ensure the suitable optimization techniques for the distributed network, different scenarios were formed and the effectiveness of the techniques was observed. The computational time of the algorithm was also taken as a crucial parameter to fix the superiority in the Wolaita Sodo distribution network. The suggested reasoning was tested on 34 radial bus systems and the results obtained were improved with MFA and PSO optimizations. The test system had a demand of 3.72 MW and 2.3 MVAR for real power and reactive capacity, respectively. The PSO technique is carried out on a radial 34-bus network, as per the standard steps. This MATLAB simulation program was developed and run on Intel Pentium dual core personal computer with 2 GB RAM. The final result obtained from the MATLAB software is shown in Table3. It is evident that the PSO approach shows its advancement and benefits over the existing method in terms of line loss decrement and computational time.Table 3
Comparison results show location and rating of SVC and power loss reduction in 34-bus system.
TechniquesMFAPSOOptimum placing of deviceBus 7Bus 11Optimal size of SVC23001352 kVARActive power loss (kW)Without SVC210.99369.57 kWWith SVC168.57270.01 kW% reduction in real power losses42.42 kW27.6Reactive power loss (kVAR)Without SVC143.13356.07 kVARWith SVC126.26260.12 kVAR% reduction in reactive power losses16.87 kVAR27.4 kVARComputational time in sec0.65 sec1.145873Figure4 clearly reiterates voltage deviation in the system while absence of SVC and voltage magnitude was restored in its optimal range with presence of SVC. In Figure 4, red colour line shows the results with SVC. Almost all 34 buses’ voltage was maintained constant with reactive power support. The convergence line of PSO drawn on 34-bus system is also depicted in Figure 5. The iteration-1 starts objective function of 298.7986 and objective function is converged at iteration-100 value of 269.32 and elapsed time is 1.145873 seconds. Before SVC placement on optimal location, real power loss (kW) is 369.57 and reactive power loss (kVAR) is 356.07. After optimal location and optimal sizing of SVC/, the real power and reactive power loss were reduced into 270.01 kW and 260.12 kVAR, respectively.Figure 4
Voltage profile with and without SVC using PSO.Figure 5
Convergence property of PSO algorithm.


## 6.2. Scenario 2

The first scenario was tested successfully with PSO algorithm. In the second scenario, recent optimization technique, moth flame algorithm, was introduced for sizing and placing the SVC. This technique coding was developed in MATLAB and run on Intel Pentium dual core like PSO code. In Figure6, results of voltage profile are constant in the entire distribution network. The range of the distribution system voltage started from 0.91 to 1.05 pu. The range in Figure 4 was accepted and matched with IEEE standard. This is one of the evidences that line losses would be minimized because of the constant voltage maintained in the system. The bus 7 was selected as weakest bus in the next for moment. The optimal size of the SVC was designed with MFA and placed at bus 7 for reactive power support. During the operation, bus 7 suffered more than the rest of the buses to maintain constant voltage at knee point.Figure 6
Convergence property of MFO.The optimal size of the SVC had been decided value of 2300 kVAR. The active power loss in kW without SVC was 210.99 kW and with SVC 168.57 kW. The percentage of the losses reduction in the system was 42.43 kW. The reactive power losses without and with SVC were 143.13 and 126.26 kVAR, respectively. The percentage of reactive power loss observed was exactly 16.87 kVAR. The overall results of the work are showcased in Table3. The comparisons were made between MFA and PSO in the all level performance. Figure 7 shows voltage profile of the 34 buses and all the bus voltage was maintained constant around 1pu after inclusion of SVC with MFA.Figure 7
Voltage profile of 34-bus system with MFA.The MFO technique had taken slight edge in terms of computational time and was given proper placing and sizing to fix the problem in the system. The power system has been more dynamical operation and uncertainty may occur any time. To prevent the faults, stability issues, compensation, and protection of the devices, the computerized monitoring system must act accordingly as fast as possible to take decision over the uncertainties. In that sense, MFA was shown its superiority in the case of improving voltage profile of the system. The big contribution of the work was helping to choose the rating of SVC during needed conditions. The operation of these proposed algorithms, MFA and PSO, would deliver exact sizing and placing of the SVC in the system. The detailed information of the collected and calculated data of Wolaita Sodo substation is recorded in Tables4 to 7. The data used in the work was measured from the case study substation. To implement the algorithm-based SVC in the system proper data from concern substation is mandatory. The MICROVIP3 PLUS device was used to measure harmonics and voltage variations in different places of the substations.Table 4
Bus data for radial distribution system.
Branch no.Nominal load at receiving busActive load power (kW)Reactive power (kVAR)S/S 100224.414.5376.837.3474.946.2540.943.2635.913.5755.91088.953.39970.142.81078.520.81143.225.71243.6123.61370.320.91478.29.561528.610.51614.87.231724.511.21825.38.621956.124.72036.114.32148.222.52220.515.62333.617.6243015.52570.220.92655.9102738.237.12843.6123.62941.820.73014085.33137.922.93225283390.836.2342110.2Table 5
Transformer data of the distribution system.
Branch no.Transformer (kVA)No-load losses (W)Full load losses (W)Z (%)S/S 10002315110036305.231002802000645035012304531582041006620089027005.2720090030004.581002802050593159503650410100520180041163045625204.5122001002200413630130065004.5141005601820415315730365041620033028605.51710044017604.5183151050420041910050018506201004401700421200120260042220060027004.5231006001550424251103504253151300430042650220100042720060027004.52810032017004.529315110037805.3301004001550431100272150543220033037004.133350220100434252106954.5Table 6
Line impedance and admittance per unit on a distribution line.
Branch no.Sending endReceiving endImpedance (Z pu)R (pu)X (pu)|y|<Y1120.012196 + 0.05756i0.01220.0575616.9959782230.02591 + 0.03344i0.025910.0334423.6389−52.23340.11651 + 0.07165i0.116510.071657.3111−31.64450.03048 + 0.05504i0.030480.0550415.8942−61.025560.16952 + 0.17598i0.169520.175984.09252−46.076670.15212 + 0.07592i0.152120.075925.88191−26.527780.00204 + 0.07014i0.002040.0701414.2512−88.338890.26552 + 0.17672i0.265520.176723.13526−33.6499100.01855 + 0.00397i0.018550.0039752.7146−121010110.09104 + 0.077876i0.091040.077888.34698−40.51111120.16519 + 0.0776i0.165190.07765.47919−25.161212130.02263 + 0.072987i0.022630.0729913.0865−72.771313140.002089 + 0.01968i0.002090.0196850.5291−83.941414150.091004 + 0.088996i0.0910.0897.85625−44.361515160.082178 + 0.02799i0.082180.0279911.5189−18.81616170.169369 + 0.070098i0.169370.07015.45548−22.481717180.15014 + 0.085191i0.150140.085195.79289−29.57182190.04984 + 0.090347i0.049840.090359.69158−61.111919200.096827 + 0.08694i0.096830.086947.68457−41.922020210.09178 + 0.135156i0.091780.135166.12097−55.82214220.09178 + 0.136356i0.091780.136366.08394−56.052222230.0251 + 0.07067i0.02510.0706713.3342−70.452323240.08867 + 0.10267i0.088670.102677.3714−49.182424250.02511 + 0.07067i0.025110.0706713.3336−70.44257260.07015 + 0.07015i0.070150.0701510.0799−452626270.173004 + 0.175156i0.1730.175164.06189−45.352727280.08999 + 0.010676i0.089990.0106811.035−6.772828290.071867 + 0.062764i0.071870.0627610.480441.132911300.082142 + 0.09494i0.082140.094947.96543−49.133230310.0929067 + 0.08735i0.092910.087357.84183−43.233114320.08735 + 0.08589i0.087350.085898.16303−44.513232330.088996 + 0.081876i0.0890.081888.2692742.613333340.152498 + 0.090676i0.15250.090685.63635−30.73Table 7
Each line real and reactive power in the system run in backward/forward MATLAB.
LinePower at busFromTokWkVAr12423.5981999.437231342.361732.404341725.7871061.46545843.74971523.622561188.4481233.752671328.319662.91777828.9178996.463389962.906640.9058910853.0899182.71741011728.5451623.1971112559.4409262.65181213170.2558549.0849131445.2561426.36981415111.1573108.7036151691.66231.2216161769.125728.6095171822.945513.019321968.2395123.7001192070.263.0337202125.536537.6061422109.423162.5714222335.884100.9834232460.74670.3374242514.39340.5043726147.9189193.52032627112.8311114.234272891.46810.8508282934.057529.7441130130.5708150.9173303131.660529.76731432107.5726105.7773323384.373677.6234333420.031611.9109The collected 34-bus data is listed in Table4. In Wolaita Sodo region, there are 34 distribution transformers under working conditions with different rating. All the transformer data are placed in Table 5. The line data of the system is catalogued in Table 6 and real and reactive power data of each bus during load flow execution is showcased in Table 7. The overall data given in Tables 4–7 could be useful for future enhancement of the realistic work in the same substation.


## 7. Conclusion

MFO algorithm is used in this article to find the optimum location, size, and SVC units to improve voltage profile and minimize the power loss in a Wolaita Sodo distribution of the 34-bus system. The findings indicated the efficiency of the algorithm, since the enforced MFO algorithm implies a lower power loss value than the literary results reported and PSO algorithm. The system's voltage profile is also improved immensely with SVCs being placed at relevant locations where more reactive power is required for voltage stable. The validity of the proposed system was noteworthy and results in different places in the articles were tabled. In this paper, effectiveness of the algorithm was compared with one of the popular algorithms PSO. Both algorithms hav effectively helped with maintaining constant voltage profile over long distances. The weakest zone of the network was bus 11 and bus 7 under different situations. The prolonged changes in the reactive power cause massive changes in the voltage which sometimes leads to brownout or blockout in the distribution ring. The moth flame optimization had been sized and placed the SVC in the optimal place for correcting the voltage deviations. The consequence of the constant voltage at all 34 buses in Wolaita Sodo leads to minimization in the line losses. The obtained results from the proposed system have proved that, under any uncertainty, conditions in Wolaita Sodo substation would be restored as soon as possible. The computational times of the algorithms were registered as 0.65 sec for MFO and 1.14 sec for PSO.



---

*Source: 9987304-2021-12-01.xml*