
Citation: Rachel R. Sleeter, William Acevedo, Christopher E. Soulard, Benjamin M. Sleeter. Methods used to parameterize the spatially-explicit components of a state-and-transition simulation model[J]. AIMS Environmental Science, 2015, 2(3): 668-693. doi: 10.3934/environsci.2015.3.668
[1] | Mikaeel Ahmadi, Mir Sayed Shah Danish, Mohammed Elsayed Lotfy, Atsushi Yona, Ying-Yi Hong, Tomonobu Senjyu . Multi-objective time-variant optimum automatic and fixed type of capacitor bank allocation considering minimization of switching steps. AIMS Energy, 2019, 7(6): 792-818. doi: 10.3934/energy.2019.6.792 |
[2] | Santosh Kumar Sharma, D. K. Palwalia, Vivek Shrivastava . Distributed generation integration optimization using fuzzy logic controller. AIMS Energy, 2019, 7(3): 337-348. doi: 10.3934/energy.2019.3.337 |
[3] | Abdollah Kavousi-Fard, Amin Khodaei . Multi-objective optimal operation of smart reconfigurable distribution grids. AIMS Energy, 2016, 4(2): 206-221. doi: 10.3934/energy.2016.2.206 |
[4] | Minh Y Nguyen . Optimal voltage controls of distribution systems with OLTC and shunt capacitors by modified particle swarm optimization: A case study. AIMS Energy, 2019, 7(6): 883-900. doi: 10.3934/energy.2019.6.883 |
[5] | HVV Priyadarshana, MA Kalhan Sandaru, KTMU Hemapala, WDAS Wijayapala . A review on Multi-Agent system based energy management systems for micro grids. AIMS Energy, 2019, 7(6): 924-943. doi: 10.3934/energy.2019.6.924 |
[6] | V. V. Thang, Thanhtung Ha . Optimal siting and sizing of renewable sources in distribution system planning based on life cycle cost and considering uncertainties. AIMS Energy, 2019, 7(2): 211-226. doi: 10.3934/energy.2019.2.211 |
[7] | Md Mashud Hyder, Kaushik Mahata . Reconfiguration of distribution system using a binary programming model. AIMS Energy, 2016, 4(3): 461-480. doi: 10.3934/energy.2016.3.461 |
[8] | Kushan Tharuka Lulbadda, K. T. M. U. Hemapala . The additional functions of smart inverters. AIMS Energy, 2019, 7(6): 971-988. doi: 10.3934/energy.2019.6.971 |
[9] | Baseem Khan, Hassan Haes Alhelou, Fsaha Mebrahtu . A holistic analysis of distribution system reliability assessment methods with conventional and renewable energy sources. AIMS Energy, 2019, 7(4): 413-429. doi: 10.3934/energy.2019.4.413 |
[10] | Zijian Hu, Hong Zhu, Chen Deng . Distribution network reconfiguration optimization method based on undirected-graph isolation group detection and the whale optimization algorithm. AIMS Energy, 2024, 12(2): 484-504. doi: 10.3934/energy.2024023 |
Distribution systems usually consist of feeders with a radial configuration in which approximately 13% of the generated power is wasted as power losses [1]. The increasing demand and load led to the development of distribution systems which causes a further voltage drop, increases losses, as a result reduction of the bus voltage stability and load imbalance. Therefore, an appropriate installation of DGs and SCBs can reduce network losses, improve network performance, and postpone investment.
Many meta-heuristics swarm techniques, which mimic the social behavior of swarms, herds, flocks, or schools of creatures in nature, have been proposed to resolve the issue of placement and sizing of DGs and SCBs. A Chaotic Artificial Bee Colony (CABC) algorithm was offered to find the best candidate locations and optimum sizes for DG units [2]. A Cuckoo Search Algorithm (CSA) was proposed to find optimal DG allocation and sizing in a radial distribution system [3]. A Whale Optimization Algorithm (WOA) was adopted for optimal allocation and sizing of DGs in RDSs [4]. In [5], the researchers used Particle Swarm Optimization (PSO), the Firefly algorithm (FA), and the Novel Bat Algorithm (NBA) methods to solve a multiobjective problem for various DG cases. A hybrid method of an improved PSO and Bat-Inspired Algorithm (BIA) for the optimal placement and sizing of DGs in RDSs was proposed in [6].
An Artificial Bee Colony (ABC) method was proposed for solving the SCBs problem [7,8,9]. A hybrid technique of the fuzzy and ABC approaches was introduced to obtain the optimal locations and sizing of SCBs [10]. Two stage approach of the fuzzy and Bat Algorithm (BA) methods was proposed to optimize the locations and sizes of SCBs in a radial distribution network [11]. The loss sensitivity and Cuckoo Search Algorithm (CSA) methods were suggested for solving the SCBs problem [12]. A Firefly Algorithm (FFA) method was offered [13,14]. The loss sensitivity and Fruit Fly optimization algorithm (FOA) methods were proposed in [15]. An Oppositional Krill Herd (OKH) approach was suggested to solve the SCBs problem in [16]. A PSO was proposed by other researchers in [17,18].
A fuzzy-dragonfly and fuzzy multi-verse optimizer methods were proposed to solve the SCBs problem in RDSs, and also a hybrid of the fuzzy and genetic algorithms was proposed to study the impact of temperature on SCBs placement in distorted RDSs [19,20,21,22]. A Multiobjective Particle Swarm Optimization (MOPSO) was suggested to determine optimal candidate locations and size of DGs and SCBs simultaneously in RDSs [23]. FFA and PSO methods were proposed to solve the DGs and SCBs placements problem in distribution systems [24,25].
Besides these swarm techniques, many remarkable researches have been carried out to solve either DGs/ SCBs placement problem or both stimulatingly. For example, A harmony search algorithm [26] was introduced to optimize the placement and sizing of DGs problem, in distribution network, aiming to improve the voltage profile, to reduce the power loss and to minimize the Total harmonic Distortion (THD). A hybrid of fuzzy decision-making tool and Multi-Objective Non-Sorting Genetic Algorithm (MONSGA) [27] was proposed to find the optimal placement and sizing of DGs in distribution networks through considering three objective functions i.e. power loss reduction, minimization of the comprehensive costs, and minimization of stability and security indices. In [28], two scenarios for DG placement in a distributions system were presented. In the first scenario, only real power loss reduction was considered. Both the optimal size and location of DGs were obtained as outputs from the system loss formula. In the next scenario, the voltage stability index was considered.
A water cycle algorithm [29] was presented to find the optimal placement and sizing of DGs and SCBs in RDSs aiming to minimize power losses, voltage deviation, total electrical energy cost, total emissions produced by generation sources and to improve the voltage stability. An elephant herding optimization algorithm [30] was introduced to find the optimal placement and sizing of DGs in distribution systems to improve the voltage profile and stability of system buses, to reduce the system loses, and to maximize the cost-savings.
PSO [31] was proposed to solve a multi-objective DGs sizing and placement problem for optimizing benefit-cost analysis of their installation and the total power losses during failures in distribution systems. Also, adaptive PSO [32] was introduced to efficiently tackle the problem of simultaneous placement of DG and CB to revamp voltage magnitude and reduce power loses in RDS. A Salp Swarm Algorithm (SSA) [33] was presented to find the optimal allocation of DGs to reduce active power loss and SCBs to improve bus voltages in RDS. Also, SSA [34] was proposed for locating and optimal sizing of DGs and SCBs on RDS for assuring the power quality (PQ) through enhancing the voltage level, minimizing the power losses, and the whole operating cost of the grid.
A hybrid of weight improved PSO and GSA algorithm [35] was proposed for optimizing the DGs and SCBs problem to achieve the reduction of power loss and enhancement of load ability in RDSs. A hybrid of fuzzy logic controller and ant-lion optimization algorithm with PSO based combination [36] was proposed to allocate the optimal placement of DGs in radial distribution systems. An enhanced genetic algorithm [37] was introduced to identify the optimal location and capacity of the simultaneous DGs/SCBs placement in the radial systems. A spring search algorithm [38] was proposed to find the optimal sizing and placement of SCBs and DGs in RDS with aiming to attain economic, technical, and environmental advantages.
In this study, a Multiobjective Grey Wolf Optimizer (MOGWO) is proposed to find the optimal DGs and SCBs placement locations and sizing simultaneously in distorted RDSs. The contributions proposed by this study are as follows:
1. To solve the DGs and SCBs placement problem by considering a a multi-objective method with Pareto fronts solutions.
2. To investigate the impact of ambient temperature on solving the DGs and SCBs placement problem.
3. To apply fuzzy set theory to select the optimal solution of non-dominated solutions.
4. To perform a MOGWO method on standard and real case studies.
5. To compare the obtained results with the MONSGA and MOPSO methods to show the effectiveness and capability of the proposed method in solving multi-objective problems when considering the temperature.
The rest of the paper is organized as follows. Section 2 describes a problem formulation which consist of objective functions, constrains, cost and temperature modelling, and a fuzzy decision-making mechanism. Section 3 introduces the proposed solution technique used in this paper. Section 4 describe results and discussion of standard and real case studies. Finally, Section 5 presents the conclusion, acknowledgments, and references.
The optimizing of DGs and SCBs placement problem involves identifying the number and locations of DGs and SCBs to be placed simultaneously on distorted RDSs. The goals and objectives of this study are to reduce the total power loss ($ P_{loss} $), enhance the Voltage Stability Index (VSI), balance the loading of the Section Currents Index (SCI), reduce the Total Harmonic Distortion (THD) and maximize the cost-savings. A simplified system model that used for obtaining the objective functions i.e., the total $ P_{loss} $, VSI, and SCI values is
The optimization problem ($ f(X) $) that need to minimized is
$ minf(X)=min[f1(X),f2(X),f3(X)] $
|
(1) |
$ X = (x_1, x_2, ..., x_n)\in R^n $ |
$ s.t. {gi(X)≥0,i=0,...,nhi(X)=0,j=0,...,n $
|
Here, $ f_1(X) = Total P_{loss} $, $ f_2(X) = SCI $, $ f_3(X) = VSI $, and $ X $ is a solution vector subjected to $ g_i(X) $ and $ h_i(X) $ which are inequality and equality constraints, and $ R^n $ is the decision variable space.
The objective functions i.e., the total $ P_{loss} $, VSI, and SCI were calculated using the load follow studies. A load flow tool called a backward-forward sweep algorithm, was used to compute the fundamental system power loss [39].
$ Pfloss=Nb∑i=2(Pgqi−Pdqi−VpiVqiYqicos(δpi−δqi+θqi)) $
|
(2) |
Here, $ P_{loss}^f $ is the total fundamental active system loss, $ Pg_{qi} $ is the active power output of the DG at bus qi, $ Pd_{qi} $ is the active power demand at bus qi, $ V_{pi} $ and $ V_{qi} $ are the voltages at buses pi and qi, respectively. $ N_{b} $ is the number of system buses, $ Y_{qi} $ is the admittance between bus pi and bus qi, $ \delta_{pi} $ is the phase angle of voltage at bus pi, $ \delta_{qi} $ is the phase angle of voltage at bus qi and $ \theta_{qi} $ is the phase angle of $ Y_i $ ($ Y_{qi}e^{j\theta_{qi}} $)
The harmonic power flow method [40] was utilized to determine the harmonic power loss for the distorted RDSs:
$ Phloss=H∑h=2Nb−1∑qi=0Rqi,qi+1[Vhqi+1,1−Vhqi][yhqi,qi+1]2, $
|
(3) |
$ TotalPloss=Phloss+Pfloss, $
|
(4) |
The system loading is calculated as follows:
$ SCI=Nb−1∑k=1|Vpi−VqiZk|, $
|
(5) |
The stability index (SI)can be evaluated as follows [41]:
$ SI(qi)=|Vpi|4−4[Pqi(qi)Xk−Qqi(qi)Rk]2−4[Pqi(qi)Rk+Qqi(qi)Xk]|Vpi|2, $
|
(6) |
$ SI(qi) $ should be maximized. Therefore, the voltage stability index (VSI) will be:
$ VSI=Nb∑i=2|1SI(qi)|, $
|
(7) |
Equations (4–7), which are used as the objective functions for solving the DGs and SCBs placement problem are subject to the following constraints:
● Bus voltage:$ V_{i}^{min}\leq V_{i}\leq V_{i}^{max} $
● Power flow:$ |PF_i| \leq PF_{i}^{max} $
● System power factor: $ |p_{f_{system}} \ \ \ | \geq p_{f_{all}}^{min} $
● Installed reactive power: $ (Q_{C}^{Total}) \leq Q_{L}^{Total} $
● The DG source utilized must be permissible in the range of the size and power factor: $ pf_{min}^{DG}\leq (pf_{qi}^{DG}\leq pf_{max}^{DG} $
● The THD must be within the allowable limit [42].
The IEEE 69-bus and 38-bus SEC RDSs (real data) were considered for this study. This was evaluated based on weather information using a least square method in which the power consumption increased by 300 MW per $ C^\circ $ based on the following equation derived from Figure 2 [43]:
$ y=300x+10400, $
|
(8) |
For an IEEE 69-bus RDS with $ P_d $ = 3.80 MW, the slope of increasing power consumption can be found by ($ 300\times $ 3.80 $ MW $)/(10.4 $ GW $) = 110 $ kW/C^\circ $. Then, the relationship between power consumed by this system and the ambient temperature is
$ y=0.110x+3.8, $
|
(9) |
For a 38-bus SEC RDS with $ P_d $ = 8.34 MW,
$ y=0.241x+8.34, $
|
(10) |
$ P_d $ and y are in MW, and $ x $ is in the interval [$ 25^\circ $ – $ 50^\circ $].
The cost of installation of SCBs and DG can be calculated as follows [44]:
$ CSCBs=NSCBs∑n=1KcQc, $
|
(11) |
$ CDGs=NDGs∑n=1CDGi+KIDG, $
|
(12) |
The cost of power from substation before $ C_{bef.plac} $ and after $ C_{aft.plac} $ installation of SCBs and DG can be calculated as follows:
$ Cbef.plac=Nyr∑n=1PWnKEMP×(real(VSI∗S))bef.×T, $
|
(13) |
$ Caft.plac=Nyr∑n=1PWnKEMP×(real(VSI∗S))aft.×T, $
|
(14) |
The presenting worth (PW) factor is formulated as follows:
$ PW=1+InfR1+IntR, $
|
(15) |
The following equations for calculating the cost saving (S) and Benefit-to-Cost Ratio (BCR) are used:
$ S=Cbef.plac−(Caft.plac+CSCBs+CDGs), $
|
(16) |
$ BCR=SCSCBs +CDGs, $
|
(17) |
$ C_{SCBs} $ is the installation cost of SCBs, $ K_C $ is the cost of commercial SCBs [45], and $ Q_C $ is the injected reactive power in $ kVAr $. $ C_{DG} $ is the cost of DG installation, $ {C_{DG}}_i $ is the installed capacity of a DG and $ K_{IDG} $ is the cost of DG investment ($318,000/MW). $ K_{EMP} $ is the energy cost ($49/MWh), $ V_s $ is the slack bus voltage, and $ I_s $ is the current fed from the slack bus to the entire system. $ PW $ is the worth factor, $ Inf_R $ is the inflation rate (2.2%), $ Int_R $ is interest rate (3%), T is 2190 h per year, and $ N_{yr} $ is 10 years.
A fuzzy decision-maker is utilized to select the optimal solution from the non-dominated solutions of MOGWO output. In the membership function $ M_{i}^{k} $, its largest value is the best optimal non-dominated solution, which computed as follows:
$ Mki=∑mi=1 μki∑Nndk=1∑mi=1 μki $
|
(18) |
$ \mu_{i}^{k} $ is a membership function for each Pareto front solution $ k $, $ N_{nd} $ is the number of non-dominated solutions, and $ m $ is the number of targeted objective functions.
The MOGWO method is suggested to identify the optimal location placement and sizing of DGs and SCBs simultaneously in a distorted RDS. This method mimics the leadership hierarchy and hunting tactic of grey wolves in nature and is guided by paired leaders ($ \alpha $), subordinate wolves ($ \beta $), subordinate wolves ($ \delta $), and scapegoat wolves ($ \omega $)[46]. The mathematical model of the suggested approach is as follow:
$ →D=|→C→Xp(t)−→X(t)|, $
|
(19) |
$ →X(t+1)=|→Xp(t)−→A⋅→D|, $
|
(20) |
$ t $ is the current iteration, $ \vec{A} $ and $ \vec{D} $ are coefficient vectors that help the GWO algorithm to explore and exploit the search space, $ \vec{X_p}(t) $ is the position vector of the prey, and $ \vec{X}(t) $ is the position vector of the grey wolf. The vectors $ \vec{A} $ and $ \vec{C} $ are calculated as follows:
$ →A=2α⋅→r1−α, $
|
(21) |
$ →C=2→r2, $
|
(22) |
The $ \alpha $ value linearly decreases from 2 to 0 during the course of the iterations and $ r_1 $ and $ r_2 $ are random vectors in [0, 1]. The range of $ \vec{C} $ is $ 2 \geq \vec{C} \geq 0 $ which improves the exploration when $ \vec{C} > 1 $ and the exploitation when $ \vec{C} < 1 $. The $ \vec{C} $ vector is generated randomly whose aim is to ensure that the exploration/exploitation at any stage can obtain the global optima while avoiding local optima.
The grey wolf optimizer begins by generating a set of random solutions as the first population. During the optimization, the three best-obtained solutions are saved and treated as $ \alpha $, $ \beta $, and $ \delta $ solutions. The next equations (21–23) update $ \alpha $, $ \beta $, and $ \delta $ wolf positions from the prey in order to simulate the hunting and discovery of promising areas of the search space:
$ →Dα=|→C1→Xα−→X|, $
|
(23) |
$ →Dβ=|→C2→Xβ−→X|, $
|
(24) |
$ →Dδ=|→C3→Xδ−→X|, $
|
(25) |
The positions of wolves in the promising areas of the search space are modelled as follows:
$ →X1=→Xα−→A1⋅→Dα, $
|
(26) |
$ →X2=→Xβ−→A2⋅→Dβ, $
|
(27) |
$ →X3=→Xδ−→A3⋅→Dδ, $
|
(28) |
$ →X(t+1)=→X1+→X2+→X33 $
|
(29) |
$ \vec{X_{\alpha}} $, $ \vec{X_{\beta}} $, and $ \vec{X_{\delta}} $ indicate wolf positions in the search space. $ \vec{X}(t+1) $ represents the random position of a wolf. The MOGWO is integrated with a repository and leader selection in order to perform the proposed multi-objective optimization. The repository is used to store the best solution (non-dominated solutions) obtained from the optimization algorithm. The leader strategy selects the better $ \alpha $, $ \beta $, and $ \delta $ wolves from the repository. Figure 3 shows a flowchart of the MOGWO technique for optimum multi-objective placement and sizing of DGs and SCBs with a fuzzy decision-maker.
The injection of harmonic currents is in the order of odd numbers (3–49) at a nominal frequency of 60 Hz, where 15% of each load at each load bus is assumed to be non-linear. The proposed method applied to two cases study with two scenarios as explained in Table 1.
No. | Scenario |
1 | No consideration of ambient temperature |
2 | Consideration of ambient temperature |
A typical 69-bus and 68-branch RDS is an IEEE standard RDS as shown in Figure 4. The system voltage and total load are 12.66 kV and 4.66 MVA, respectively [47].
In this scenario, the impact of ambient temperature was not considered. The MONSGA, MOPSO and MOGWO algorithms discussed above were implemented for this system while considering all the objective functions and constraints. Various optimal results such as the optimal location, size, VSI, SCI and THD values, real power losses, reduction cost of power and energy losses and benefit-cost ratio (BCR) before and after installing of DGs and SCBs in the distorted 69-bus RDS were calculated. It can be observed that for the base case (before the installation of DGs and SCBs), the real power loss was 232.92 kW. By applying the MOGWO method, this loss was reduced to 30.13 kW, which was less than that in other methods as shown in Table 3. The VSI, SCI, and THD had permissible values which were 1.058 p.u., 0.331 p.u., and 2.20, respectively. The total annual reduction cost of the power losses was $ 1,265,500 for a cost-saving of $ 1,000,000. It can be seen that the calculated BCR was 1.22.
Method | Base Case | MONSGA | MOPSO | Proposed MOGWO | ||||
$ T_{No} $ | $ T_{With} $ | $ T_{No} $ | $ T_{With} $ | $ T_{No} $ | $ T_{With} $ | $ T_{No} $ | $ T_{With} $ | |
Location bus for DG | 0 | 0 | 62 | 61 | 8 | 17 | 54 | 61 |
0 | 0 | 61 | 62 | 62 | 63 | 62 | 62 | |
0 | 0 | 53 | — | 11 | — | 19 | — | |
Size of DG (MW) | 0 | 0 | 1.00 | 1.22 | 1.01 | 1.24 | 0.601 | 1.15 |
0 | 0 | 0.80 | 1.10 | 1.24 | 0.762 | 1.22 | 1.22 | |
0 | 0 | 0.985 | — | 0.672 | — | 0.54 | — | |
Total MW | 0 | 0 | 2.79 | 2.32 | 2.93 | 2.00 | 2.58 | 2.37 |
Location bus for SCB | 0 | 0 | 37 | 32 | 8 | 62 | 46 | 61 |
0 | 0 | 61 | 34 | 40 | 63 | 37 | 65 | |
0 | 0 | 28 | — | 17 | — | 62 | — | |
Size of SCB (MVAr) | 0 | 0 | 1.10 | 0.184 | 2.20 | 0.615 | 0.470 | 2.07 |
0 | 0 | 1.17 | 1.122 | 0.150 | 0.710 | 0.841 | 0.370 | |
0 | 0 | 0.448 | — | 0.550 | — | 1.25 | — | |
Total MVAr | 0 | 0 | 2.71 | 1.406 | 2.90 | 1.33 | 2.56 | 2.443 |
Max. voltage (p.u) | 1.032 | 1.026 | 1.01 | 1.02 | 1.004 | 1.01 | 1.01 | 1.01 |
Min. voltage (p.u) | 0.9438 | 0.9010 | 0.980 | 0.9537 | 0.9758 | 0.9386 | 0.9863 | 0.9689 |
VSI (p.u) | 1.283 | 1.581 | 1.092 | 1.16 | 1.102 | 1.30 | 1.058 | 1.15 |
SCI (p.u) | 0.728 | 0.976 | 0.310 | 0.611 | 0.327 | 0.615 | 0.331 | 0.449 |
THD | 9.92 | 9.23 | 1.96 | 6.10 | 1.82 | 4.28 | 2.20 | 2.95 |
Power losses (kW) | 232.9 | 438.48 | 31.43 | 155.73 | 69.19 | 152.87 | 30.13 | 105 |
Reduction (%) | 0 | 0 | 87.54 | 64.50 | 70.30 | 65.14 | 87.1 | 76.10 |
DG cost ($K/y) | 0 | 0 | 887.790 | 738.480 | 930.970 | 635.810 | 820.460 | 752.520 |
SCB cost ($K/y) | 0 | 0 | 3.522 | 2.384 | 3.518 | 2.414 | 3.444 | 2.483 |
Power cost ($M/y) | 3.100 | 3.980 | 1.159 | 2.715 | 1.062 | 2.795 | 1.266 | 2.196 |
Benefit ($M/y) | 0 | 0 | 1.045 | 0.522 | 1.090 | 0.544 | 1.0 | 1.027 |
BCR | 0 | 0 | 1.17 | 0.70 | 1.18 | 0.85 | 1.22 | 1.36 |
The impact of ambient temperature on this system was studied using Eq 9. One observes that in the base case, the real power loss was 438.48 kW and reduced to 105 kW after applying the MOGWO method which was less than in the other methods as shown in Table 3. The performance of the MOGWO method was compared to other methods through its elapsed time needed to converge as explained in Table 2 for both scenarios. The indices of VSI, SCI, and THD had permissible values comparing to other methods; these values were 1.15 p.u., 0.449 p.u., and 2.95, respectively. Figure 5 shows the voltage profile of this system. The figure compares the magnitudes of the voltage profile before and after installing DGs and SCBs for both scenarios. In the base case, the voltage profile was poor, which violates the lower limit of the voltage constraint. However, this voltage profile was improved after the installation of DGs and SCBs. Figure 6 shows VSI values on buses of the system, and Figure 7 shows the section current for also both scenarios. Almost all currents in all sections were decreased after installing of DGs and SCBs. The total annual reduction cost of the power losses was $ 2,195,500 for a cost-saving of $ 1,026,800. The calculated BCR was 1.36.
Scenario | MONSGA | MOPSO | Proposed MOGWO |
No Temperature ($ T_{No} $) | 621.38 | 584.76 | 565.18 |
With Temperature ($ T_{With} $) | 673 | 648.92 | 589.78 |
A typical 38-bus and 37-branch RDS as shown in Figure 8, is a SEC distribution system. The system voltage and total load are 13.8 kV and 8.77 MVA, respectively.
The MONSGA, MOPSO and MOGWO algorithms were also implemented for this system while considering all the objective functions and constraints. Various optimal results such as the optimal location, size, VSI, SCI and THD values, real power losses, reduction cost of power and energy losses and BCR before and after installing of DGs and SCBs in the distorted 38-bus RDS were calculated. It can be observed that for the base case (before the installation of DGs and SCBs), the real power loss was 194.11 kW as shown in Table 5. By applying the MOGWO method, this loss was reduced to 128.12 kW, which was less than that in other methods. The VSI, SCI, and THD had permissible values which were 1.10 p.u., 0.220 p.u., and 2.29, respectively. The total annual reduction cost of the power losses was $ 2,298,500 for a cost-saving of $ 473,030. It can be seen that the calculated BCR was 1.46
Method | Base Case | MONSGA | MOPSO | Proposed MOGWO | ||||
$ T_{No} $ | $ T_{With} $ | $ T_{No} $ | $ T_{With} $ | $ T_{No} $ | $ T_{With} $ | $ T_{No} $ | $ T_{With} $ | |
Location bus for DG | 0 | 0 | 27 | 27 | 6 | 6 | 25 | 24 |
0 | 0 | 21 | 26 | 27 | 22 | 26 | 26 | |
Size of DG (MW) | 0 | 0 | 0.509 | 0.415 | 0.505 | 0.509 | 0.509 | 0.509 |
0 | 0 | 0.4710 | 0.368 | 0.509 | 0.509 | 0.505 | 0.509 | |
Total MW | 0 | 0 | 0.976 | 0.783 | 1.13 | 1.017 | 1.13 | 1.017 |
Location bus for SCB | 0 | 0 | 12 | 23 | 4 | 25 | 33 | 33 |
0 | 0 | 23 | 17 | 38 | 34 | 21 | 21 | |
Size of SCB (MVAr) | 0 | 0 | 0.962 | 1.434 | 1.64 | 1.63 | 0.824 | 0.98 |
0 | 0 | 1.700 | 1.310 | 1.10 | 0.960 | 1.54 | 1.70 | |
Total MVAr | 0 | 0 | 2.66 | 2.744 | 2.74 | 2.59 | 2.36 | 2.68 |
Max. voltage (p.u) | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 |
Min. voltage (p.u) | 0.9603 | 0.9466 | 0.9765 | 0.9666 | 0.9764 | 0.9690 | 0.9764 | 0.9687 |
VSI (p.u) | 1.19 | 1.27 | 1.10 | 1.16 | 1.10 | 1.14 | 1.10 | 1.14 |
SCI (p.u) | 0.314 | 0.425 | 0.218 | 0.327 | 0.213 | 0.30 | 0.220 | 0.31 |
THD | 3.17 | 3.97 | 2.41 | 3.31 | 2.47 | 3.43 | 2.29 | 3.05 |
Power losses (kW) | 194.11 | 403.83 | 132 | 318.76 | 128.36 | 304.18 | 128.12 | 300 |
Reduction (%) | 0 | 0 | 32 | 21.1 | 34.24 | 24.68 | 33..95 | 25.6 |
DG cost ($K/y) | 0 | 0 | 310.030 | 248.980 | 322.360 | 323.410 | 320.087 | 323.410 |
SCB cost ($K/y) | 0 | 0 | 2.576 | 2.581 | 2.522 | 2.558 | 2.576 | 2.576 |
Power cost ($M/y) | 3.900 | 5.470 | 2.416 | 4.222 | 2.439 | 4 | 2.197 | 3.943 |
Benefit ($M/y) | 0 | 0 | 0.365 | 1 | 0.331 | 1.155 | 0.473 | 1.204 |
BCR | 0 | 0 | 1.17 | 3.97 | 1.02 | 3.54 | 1.46 | 3.70 |
In the same way, the temperature effects on this system were studied using Eq 10. In the base case, the real power loss was 403.83 kW and reduced to 300 kW after applying the MOGWO method which was less than in the other methods as shown in Table 5. The performance of the MOGWO method was compared to other methods through its elapsed time needed to converge as explained in Table 4 for both scenarios. The indices of VSI, SCI, and THD had permissible values comparing to other methods; these values were 1.11 p.u., 0.310 p.u., and 3.05, respectively. Figure 9 shows the voltage profile of this system. The figure compares the magnitudes of the voltage profile before and after installing DGs and SCBs for both scenarios. In the base case, the voltage profile was poor, which violates the lower limit of the voltage constraint. However, this voltage profile was improved after the installation of DGs and SCBs. Figure 10 shows VSI values on buses of the system, and Figure 11 shows the section current for also both scenarios. Almost all currents in all sections were decreased after installing of DGs and SCBs. The total annual reduction cost of the power losses was $ 3.942,600 for a cost-saving of $ 1,204,100. The calculated BCR was 3.70.
Scenario | MONSGA | MOPSO | Proposed MOGWO |
No Temperature ($ T_{No} $) | 243.65 | 238.26 | 221.34 |
With Temperature ($ T_{With} $) | 296.52 | 289.28 | 269.98 |
The MOGWO technique was proposed to identify the best candidate locations and optimum sizing of DGs and SCBs simultaneously under the impact of ambient temperature. The power loss reduction, voltage stability improvement of the system nodes, and load balancing in the sections were simulated as a multi-objective function. The best optimal solution was chosen via a fuzzy-based mechanism. The proposed method was applied to an IEEE 69-bus and 38-bus RDSs (real data from SEC). The obtained results were compared with the MONSGA and MOPSO methods, and showed better results in terms of system power loss reduction, voltage stability and voltage profile improvement, THD level, and load balancing. Additionally, economic benefits such as cost savings and BCR were evaluated. It is recommended that the power utilities adopt this technique to improve their technical capability and financial figures.
The authors extend their appreciation to the Deanship of Scientific Research, King Saud University for funding this work through research group no. (RG-1439-028).
The authors declare no conflict of interest in this paper.
[1] |
Vitousek PM, Mooney HA, Lubchenco J, et al. (1997) Human domination of earth's ecosystems. Science 277: 494-499. doi: 10.1126/science.277.5325.494
![]() |
[2] | DeFries RS, Foley JA, Asner GP (2004) Land-use choices: Balancing human needs and ecosystem function. Front Ecol Environ 2: 249-257. |
[3] | IPCC (2000) In: Watson, R.T., Noble, I.R., Bolin, B., Ravindranath, N.H., Verardo, D.J., Dokken, D.J. (Eds.), Land Use, Land-Use Change, and Forestry. Special Report of the Intergovernmental Panel on Climate Change. Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, p. 377. |
[4] | Foley JA, Defries R, Asner GP, et al. (2005) Global consequences of land use. Science 309: 570-574. |
[5] |
Pielke RA (2005) Land use and climate change. Science 310: 1625-1626. doi: 10.1126/science.1120529
![]() |
[6] | Alcamo J (2008) Searching for the future of land: scenarios from the local to global scale. In: Alcamo J (Ed.), Environmental Futures: The Practice of Environmental Scenario Analysis. Elsevier, Amsterdam, The Netherlands. |
[7] | USGCRP, The National Global Change Research Plan 2012-2021: A Strategic Plan for the U.S. Global Change Research Program. 132 pp. The U.S. Global Change Research Program Washington, D.C., 2012, Available from: http://downloads.globalchange.gov/strategic-plan/2012/usgcrp-strategic-plan-2012.pdf |
[8] |
Strengers B, Leemans R, Eickhout B, et al. (2004) The land-use projections and resulting emissions in the IPCC SRES scenarios as simulated by the IMAGE 2.2 model. GeoJournal 61: 381-393. doi: 10.1007/s10708-004-5054-8
![]() |
[9] | National Research Council (2014) Advancing Land Change Modeling: Opportunities and Research Requirements. Washington, DC: The National Academies Press. |
[10] | Soares-Filho BS, Nepstad DC, Curran LM, et al. (2006) Modeling conservation in the Amazon basin. Nature 440: 520-523. |
[11] | Eastman JR (2007) The Land Change Modeler, a software extension for ArcGIS. Worcester, Mass.: Clark University. |
[12] | Matthews RB, Gilbert NG, Roach A, et al. (2007) Agent-based land-use models: A review of applications. Landscape Ecol 22: 1447-1459. |
[13] |
Acosta-Michlik L, Espaldon V (2008) Assessing vulnerability of selected farming communities in the Philippines based on behavioural model of agent's adaptation to global environmental change. Glob Environ Chang 18: 554-563. doi: 10.1016/j.gloenvcha.2008.08.006
![]() |
[14] | Brown DG, Page SE, Riolo R, et al. (2004) Agent based and analytical modeling to evaluate the effectiveness of greenbelts. Environ Modell Softw 19: 1097-1109. |
[15] |
Verburg PH (2002) Land use change modelling at the regional scale: The CLUE-S model. Environ Manage 30: 391-405. doi: 10.1007/s00267-002-2630-x
![]() |
[16] | Burnham BO (1973) Markov intertemporal land use simulation model. South J Agr Econ 5: 253-258. |
[17] | Baker WL (1989) A review of models of landscape change. Landscape Ecol 2: 111-133. |
[18] | Turner M (1987) Spatial simulation of landscape changes in Georgia: A comparison of 3 transition models. Landscape Ecol 1: 29-36. |
[19] | Daniels CJ, Frid L (2011) Predicting landscape vegetation dynamics using state-and-transition simulation models, In Proceedings of the First Landscape State-and-Transition Simulation Modeling Conference, Portland, OR, USA, 14-16 June 2011; Kerns, B.K., Shlisky, A.J., Daniel, C.J., Eds.; U.S. Department of Agriculture, Forest Service, Pacific Northwest Research Station: Portland, OR, USA, 2012; pp.5-22. |
[20] | Sleeter BM, Sohl T, Bouchard MA, et al. (2012) Scenarios of land use and land cover change in the conterminous U.S.: Utilizing the special report on emission scenarios at ecoregional scales. Glob Environ Chang 22: 896-914. |
[21] |
Omernik JM (1987) Ecoregions of the conterminous U.S. Ann Assoc Am Geogr 77: 118-125. doi: 10.1111/j.1467-8306.1987.tb00149.x
![]() |
[22] | Sleeter BM, Liu J, Daniel C, et al. (2015) An integrated approach to modeling changes in land use, land cover, and disturbance and their impact on ecosystem carbon dynamics: a case study in the Sierra Nevada Mountains of California. AIMS Environ Sci 2: 577-606. |
[23] | Nakicenovic N, Swart R (Eds.) (2000) IPCC Special Report on Emission Scenarios; Cambridge University Press; Cambridge, UK; p. 570. |
[24] |
Sleeter BM, Sohl T, Loveland T, et al. (2013) Land-cover change in the conterminous United States from 1973 to 2000. Glob Environ Chang 23: 733-748. doi: 10.1016/j.gloenvcha.2013.03.006
![]() |
[25] | EPA (U.S. Environmental Protection Agency), Primary distinguishing characteristics of Level III ecoregions of the continental United States. Environmental Protection Agency, 1999. Available from: http://www.epa.gov/wed/pages/ecoregions/level_iii.htm |
[26] | Gallant AL, Loveland T, Sohl T, et al. (2004) Using an ecoregion framework to analyze land-cover and land-use dynamics. Environ Manage 34: S89-S110. |
[27] | Soulard CE and Acevedo W, Multi-temporal harmonization of independent land-use/land-cover datasets for the conterminous U.S. American Geophysical Union, 2013. Available from: http://adsabs.harvard.edu/abs/2013AGUFM.B41E0448S |
[28] | McConnell WJ, Moran EF (Eds.) (2001) Meeting in the middle: the challenge of meso-level integration. An international workshop on the harmonization of land-use and land-cover classification. LUCC Report Series No. 5. Anthropological Center for Training and Research on Global Environmental Change - Indiana University and LUCC International Project Office, Louvain-la-Neuve. |
[29] |
Jansen LJM, Groom GB, Carrai G (2008) Land-cover harmonisation and semantic similarity: some methodological issues. J Land Use Sci 3:131-160. doi: 10.1080/17474230802332076
![]() |
[30] | Homer C, Dewitz J, Fry J, et al. (2007) Completion of the 2001 national land cover database for the conterminous U.S. Photogramm Eng Remote Sens 73: 337-341. |
[31] | Vogelmann JE, Howard SM, Yang L, et al. (2001) Completion of the 1990's national land cover data set for the conterminous U.S. Photogramm Eng Remote Sens 67: 650-662. |
[32] | Fry J, Xian G, Jin S, et al. (2011) Completion of the 2006 national land cover database for the conterminous U.S. Photogramm Eng Remote Sens 77:858-864. |
[33] |
Jin S, Yang L, Danielson P (2013) A comprehensive change detection method for updating the National Land Cover Database to circa 2011. Remote Sens Environ 132: 159-175. doi: 10.1016/j.rse.2013.01.012
![]() |
[34] | LANDFIRE: LANDFIRE Existing Vegetation Type layer. U.S. Department of Interior, Geological Survey, 2013 June last update. Available from: http://landfire.cr.usgs.gov/viewer/ (accessed 1 October 2012) |
[35] | U.S Department of Agriculture, National Agricultural Statistics Service (2011) Cropland Data Layer. Available from http://nassgeodata.gmu.edu/CropScape/ (accessed on 1 October 2012). |
[36] | U.S. Geological Survey, Gap Analysis Program (GAP). May 2011. National Land Cover, Version 2 |
[37] |
Hansen MC, Potapov PV, Moore R, et al. (2013) High-Resolution Global Maps of 21st-Century Forest Cover Change. Science 342: 850-853. doi: 10.1126/science.1244693
![]() |
[38] | LANDFIRE: LANDFIRE Disturbance (1999-2010). U.S. Department of Interior, Geological Survey, 2013 June last update. Available from: http://landfire.cr.usgs.gov/viewer/ (accessed on 1 October 2012). |
[39] |
Roy DP, Ju J, Kline K, et al. (2010) Web-enabled Landsat data (WELD): Landsat ETM+ composited mosaics of the conterminous U.S. Remote Sens Environ 114: 35-49. doi: 10.1016/j.rse.2009.08.011
![]() |
[40] | Eidenshink J, Schwind B, Brewer K, et al. (2007) A project for monitoring trends in burn severity. Fire Ecol Spec Issue 3: 3-21. |
[41] | Anderson JR, Hardy E, Roach JT, et al. (1976) A land use and land cover classification scheme for use with remote sensor data. U.S. Geological Survey, Reston, VA. USA, Professional Paper 964. |
[42] |
Pan Y, Chen JM, Birdsey R, et al. (2011) Age structure and disturbance legacy of North American forests. Biogeosciences 8: 715-732. doi: 10.5194/bg-8-715-2011
![]() |
[43] |
Masek JG, Huang R, Wolfe R, et al. (2008) North American forest disturbance mapped from a decadal Landsat record. Remote Sens Environ 112: 2914- 2926. doi: 10.1016/j.rse.2008.02.010
![]() |
[44] |
Huang C, Goward SN, Masek JG, et al. (2010) An automated approach for reconstructing recent forest disturbance history using dense Landsat time series stacks. Remote Sens Environ 114: 183-198. doi: 10.1016/j.rse.2009.08.017
![]() |
[45] |
Tobler W (1970) A computer movie simulating urban growth in the Detroit region. Econ Geogr 46: 234-240. doi: 10.2307/143141
![]() |
[46] | U.S. Geological Survey GAP. Protected Areas Database of the United States (PAD-US), version 1.3 Combined Feature Class. 2012. Available from: http://gapanalysis.usgs.gov/padus/ (accessed on 13 November 2013). |
[47] | Natural Resources Conservation Service. Soil Survey Geographic (SSURGO) Database. U.S. Department of Agriculture, 2011. Available from: http://soildatamart.nrcs.usda.gov (accessed on 22 November 2013). |
[48] | Sleeter R, Gould M (2007) Geographic Information System Software to Remodel Population Data Using Dasymetric Mapping Methods; U.S. Geological Survey Techniques and Methods 11-C2; U.S. Geological Survey, Reston, VA. USA, p. 15. |
[49] | Wilson BT, Woodall CW, DM Griffith (2013) Imputing forest carbon stock estimates from inventory plots to a nationally continuous coverage. Carbon Balance Manag 8: 1-15. |
[50] | IMAGE team (2001) The IMAGE 2.2 implementation of the SRES scenarios: climate change scenarios resulting from runs with several GCMs. RIVM CD-ROM Publication 481508019, National Institute of Public Health and the Environment, Bilthoven. |
[51] |
Van Vuuren D, Edmonds JA, Kainuma M, et al. (2011) The Representative Concentration Pathways: An Overview. Climatic Change 109: 5-31. doi: 10.1007/s10584-011-0148-z
![]() |
[52] |
Wilson TS, Sleeter BM, Sleeter RR, et al. (2014) Land use threats and protected areas: a scenario-based, landscape level approach. Land 3: 362-389. doi: 10.3390/land3020362
![]() |
[53] | Sleeter BM (2008) Late 20th century land change in the Central California Valley Ecoregion. California Geographer 48:27-60 |
1. | Waleed Fadel, A Review of Distributed Generation Optimization with Shunt Capacitors in Reconfigured Distribution Networks, 2024, 136, 0929-6212, 1, 10.1007/s11277-024-11062-x |
No. | Scenario |
1 | No consideration of ambient temperature |
2 | Consideration of ambient temperature |
Method | Base Case | MONSGA | MOPSO | Proposed MOGWO | ||||
$ T_{No} $ | $ T_{With} $ | $ T_{No} $ | $ T_{With} $ | $ T_{No} $ | $ T_{With} $ | $ T_{No} $ | $ T_{With} $ | |
Location bus for DG | 0 | 0 | 62 | 61 | 8 | 17 | 54 | 61 |
0 | 0 | 61 | 62 | 62 | 63 | 62 | 62 | |
0 | 0 | 53 | — | 11 | — | 19 | — | |
Size of DG (MW) | 0 | 0 | 1.00 | 1.22 | 1.01 | 1.24 | 0.601 | 1.15 |
0 | 0 | 0.80 | 1.10 | 1.24 | 0.762 | 1.22 | 1.22 | |
0 | 0 | 0.985 | — | 0.672 | — | 0.54 | — | |
Total MW | 0 | 0 | 2.79 | 2.32 | 2.93 | 2.00 | 2.58 | 2.37 |
Location bus for SCB | 0 | 0 | 37 | 32 | 8 | 62 | 46 | 61 |
0 | 0 | 61 | 34 | 40 | 63 | 37 | 65 | |
0 | 0 | 28 | — | 17 | — | 62 | — | |
Size of SCB (MVAr) | 0 | 0 | 1.10 | 0.184 | 2.20 | 0.615 | 0.470 | 2.07 |
0 | 0 | 1.17 | 1.122 | 0.150 | 0.710 | 0.841 | 0.370 | |
0 | 0 | 0.448 | — | 0.550 | — | 1.25 | — | |
Total MVAr | 0 | 0 | 2.71 | 1.406 | 2.90 | 1.33 | 2.56 | 2.443 |
Max. voltage (p.u) | 1.032 | 1.026 | 1.01 | 1.02 | 1.004 | 1.01 | 1.01 | 1.01 |
Min. voltage (p.u) | 0.9438 | 0.9010 | 0.980 | 0.9537 | 0.9758 | 0.9386 | 0.9863 | 0.9689 |
VSI (p.u) | 1.283 | 1.581 | 1.092 | 1.16 | 1.102 | 1.30 | 1.058 | 1.15 |
SCI (p.u) | 0.728 | 0.976 | 0.310 | 0.611 | 0.327 | 0.615 | 0.331 | 0.449 |
THD | 9.92 | 9.23 | 1.96 | 6.10 | 1.82 | 4.28 | 2.20 | 2.95 |
Power losses (kW) | 232.9 | 438.48 | 31.43 | 155.73 | 69.19 | 152.87 | 30.13 | 105 |
Reduction (%) | 0 | 0 | 87.54 | 64.50 | 70.30 | 65.14 | 87.1 | 76.10 |
DG cost ($K/y) | 0 | 0 | 887.790 | 738.480 | 930.970 | 635.810 | 820.460 | 752.520 |
SCB cost ($K/y) | 0 | 0 | 3.522 | 2.384 | 3.518 | 2.414 | 3.444 | 2.483 |
Power cost ($M/y) | 3.100 | 3.980 | 1.159 | 2.715 | 1.062 | 2.795 | 1.266 | 2.196 |
Benefit ($M/y) | 0 | 0 | 1.045 | 0.522 | 1.090 | 0.544 | 1.0 | 1.027 |
BCR | 0 | 0 | 1.17 | 0.70 | 1.18 | 0.85 | 1.22 | 1.36 |
Scenario | MONSGA | MOPSO | Proposed MOGWO |
No Temperature ($ T_{No} $) | 621.38 | 584.76 | 565.18 |
With Temperature ($ T_{With} $) | 673 | 648.92 | 589.78 |
Method | Base Case | MONSGA | MOPSO | Proposed MOGWO | ||||
$ T_{No} $ | $ T_{With} $ | $ T_{No} $ | $ T_{With} $ | $ T_{No} $ | $ T_{With} $ | $ T_{No} $ | $ T_{With} $ | |
Location bus for DG | 0 | 0 | 27 | 27 | 6 | 6 | 25 | 24 |
0 | 0 | 21 | 26 | 27 | 22 | 26 | 26 | |
Size of DG (MW) | 0 | 0 | 0.509 | 0.415 | 0.505 | 0.509 | 0.509 | 0.509 |
0 | 0 | 0.4710 | 0.368 | 0.509 | 0.509 | 0.505 | 0.509 | |
Total MW | 0 | 0 | 0.976 | 0.783 | 1.13 | 1.017 | 1.13 | 1.017 |
Location bus for SCB | 0 | 0 | 12 | 23 | 4 | 25 | 33 | 33 |
0 | 0 | 23 | 17 | 38 | 34 | 21 | 21 | |
Size of SCB (MVAr) | 0 | 0 | 0.962 | 1.434 | 1.64 | 1.63 | 0.824 | 0.98 |
0 | 0 | 1.700 | 1.310 | 1.10 | 0.960 | 1.54 | 1.70 | |
Total MVAr | 0 | 0 | 2.66 | 2.744 | 2.74 | 2.59 | 2.36 | 2.68 |
Max. voltage (p.u) | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 |
Min. voltage (p.u) | 0.9603 | 0.9466 | 0.9765 | 0.9666 | 0.9764 | 0.9690 | 0.9764 | 0.9687 |
VSI (p.u) | 1.19 | 1.27 | 1.10 | 1.16 | 1.10 | 1.14 | 1.10 | 1.14 |
SCI (p.u) | 0.314 | 0.425 | 0.218 | 0.327 | 0.213 | 0.30 | 0.220 | 0.31 |
THD | 3.17 | 3.97 | 2.41 | 3.31 | 2.47 | 3.43 | 2.29 | 3.05 |
Power losses (kW) | 194.11 | 403.83 | 132 | 318.76 | 128.36 | 304.18 | 128.12 | 300 |
Reduction (%) | 0 | 0 | 32 | 21.1 | 34.24 | 24.68 | 33..95 | 25.6 |
DG cost ($K/y) | 0 | 0 | 310.030 | 248.980 | 322.360 | 323.410 | 320.087 | 323.410 |
SCB cost ($K/y) | 0 | 0 | 2.576 | 2.581 | 2.522 | 2.558 | 2.576 | 2.576 |
Power cost ($M/y) | 3.900 | 5.470 | 2.416 | 4.222 | 2.439 | 4 | 2.197 | 3.943 |
Benefit ($M/y) | 0 | 0 | 0.365 | 1 | 0.331 | 1.155 | 0.473 | 1.204 |
BCR | 0 | 0 | 1.17 | 3.97 | 1.02 | 3.54 | 1.46 | 3.70 |
Scenario | MONSGA | MOPSO | Proposed MOGWO |
No Temperature ($ T_{No} $) | 243.65 | 238.26 | 221.34 |
With Temperature ($ T_{With} $) | 296.52 | 289.28 | 269.98 |
No. | Scenario |
1 | No consideration of ambient temperature |
2 | Consideration of ambient temperature |
Method | Base Case | MONSGA | MOPSO | Proposed MOGWO | ||||
$ T_{No} $ | $ T_{With} $ | $ T_{No} $ | $ T_{With} $ | $ T_{No} $ | $ T_{With} $ | $ T_{No} $ | $ T_{With} $ | |
Location bus for DG | 0 | 0 | 62 | 61 | 8 | 17 | 54 | 61 |
0 | 0 | 61 | 62 | 62 | 63 | 62 | 62 | |
0 | 0 | 53 | — | 11 | — | 19 | — | |
Size of DG (MW) | 0 | 0 | 1.00 | 1.22 | 1.01 | 1.24 | 0.601 | 1.15 |
0 | 0 | 0.80 | 1.10 | 1.24 | 0.762 | 1.22 | 1.22 | |
0 | 0 | 0.985 | — | 0.672 | — | 0.54 | — | |
Total MW | 0 | 0 | 2.79 | 2.32 | 2.93 | 2.00 | 2.58 | 2.37 |
Location bus for SCB | 0 | 0 | 37 | 32 | 8 | 62 | 46 | 61 |
0 | 0 | 61 | 34 | 40 | 63 | 37 | 65 | |
0 | 0 | 28 | — | 17 | — | 62 | — | |
Size of SCB (MVAr) | 0 | 0 | 1.10 | 0.184 | 2.20 | 0.615 | 0.470 | 2.07 |
0 | 0 | 1.17 | 1.122 | 0.150 | 0.710 | 0.841 | 0.370 | |
0 | 0 | 0.448 | — | 0.550 | — | 1.25 | — | |
Total MVAr | 0 | 0 | 2.71 | 1.406 | 2.90 | 1.33 | 2.56 | 2.443 |
Max. voltage (p.u) | 1.032 | 1.026 | 1.01 | 1.02 | 1.004 | 1.01 | 1.01 | 1.01 |
Min. voltage (p.u) | 0.9438 | 0.9010 | 0.980 | 0.9537 | 0.9758 | 0.9386 | 0.9863 | 0.9689 |
VSI (p.u) | 1.283 | 1.581 | 1.092 | 1.16 | 1.102 | 1.30 | 1.058 | 1.15 |
SCI (p.u) | 0.728 | 0.976 | 0.310 | 0.611 | 0.327 | 0.615 | 0.331 | 0.449 |
THD | 9.92 | 9.23 | 1.96 | 6.10 | 1.82 | 4.28 | 2.20 | 2.95 |
Power losses (kW) | 232.9 | 438.48 | 31.43 | 155.73 | 69.19 | 152.87 | 30.13 | 105 |
Reduction (%) | 0 | 0 | 87.54 | 64.50 | 70.30 | 65.14 | 87.1 | 76.10 |
DG cost ($K/y) | 0 | 0 | 887.790 | 738.480 | 930.970 | 635.810 | 820.460 | 752.520 |
SCB cost ($K/y) | 0 | 0 | 3.522 | 2.384 | 3.518 | 2.414 | 3.444 | 2.483 |
Power cost ($M/y) | 3.100 | 3.980 | 1.159 | 2.715 | 1.062 | 2.795 | 1.266 | 2.196 |
Benefit ($M/y) | 0 | 0 | 1.045 | 0.522 | 1.090 | 0.544 | 1.0 | 1.027 |
BCR | 0 | 0 | 1.17 | 0.70 | 1.18 | 0.85 | 1.22 | 1.36 |
Scenario | MONSGA | MOPSO | Proposed MOGWO |
No Temperature ($ T_{No} $) | 621.38 | 584.76 | 565.18 |
With Temperature ($ T_{With} $) | 673 | 648.92 | 589.78 |
Method | Base Case | MONSGA | MOPSO | Proposed MOGWO | ||||
$ T_{No} $ | $ T_{With} $ | $ T_{No} $ | $ T_{With} $ | $ T_{No} $ | $ T_{With} $ | $ T_{No} $ | $ T_{With} $ | |
Location bus for DG | 0 | 0 | 27 | 27 | 6 | 6 | 25 | 24 |
0 | 0 | 21 | 26 | 27 | 22 | 26 | 26 | |
Size of DG (MW) | 0 | 0 | 0.509 | 0.415 | 0.505 | 0.509 | 0.509 | 0.509 |
0 | 0 | 0.4710 | 0.368 | 0.509 | 0.509 | 0.505 | 0.509 | |
Total MW | 0 | 0 | 0.976 | 0.783 | 1.13 | 1.017 | 1.13 | 1.017 |
Location bus for SCB | 0 | 0 | 12 | 23 | 4 | 25 | 33 | 33 |
0 | 0 | 23 | 17 | 38 | 34 | 21 | 21 | |
Size of SCB (MVAr) | 0 | 0 | 0.962 | 1.434 | 1.64 | 1.63 | 0.824 | 0.98 |
0 | 0 | 1.700 | 1.310 | 1.10 | 0.960 | 1.54 | 1.70 | |
Total MVAr | 0 | 0 | 2.66 | 2.744 | 2.74 | 2.59 | 2.36 | 2.68 |
Max. voltage (p.u) | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 |
Min. voltage (p.u) | 0.9603 | 0.9466 | 0.9765 | 0.9666 | 0.9764 | 0.9690 | 0.9764 | 0.9687 |
VSI (p.u) | 1.19 | 1.27 | 1.10 | 1.16 | 1.10 | 1.14 | 1.10 | 1.14 |
SCI (p.u) | 0.314 | 0.425 | 0.218 | 0.327 | 0.213 | 0.30 | 0.220 | 0.31 |
THD | 3.17 | 3.97 | 2.41 | 3.31 | 2.47 | 3.43 | 2.29 | 3.05 |
Power losses (kW) | 194.11 | 403.83 | 132 | 318.76 | 128.36 | 304.18 | 128.12 | 300 |
Reduction (%) | 0 | 0 | 32 | 21.1 | 34.24 | 24.68 | 33..95 | 25.6 |
DG cost ($K/y) | 0 | 0 | 310.030 | 248.980 | 322.360 | 323.410 | 320.087 | 323.410 |
SCB cost ($K/y) | 0 | 0 | 2.576 | 2.581 | 2.522 | 2.558 | 2.576 | 2.576 |
Power cost ($M/y) | 3.900 | 5.470 | 2.416 | 4.222 | 2.439 | 4 | 2.197 | 3.943 |
Benefit ($M/y) | 0 | 0 | 0.365 | 1 | 0.331 | 1.155 | 0.473 | 1.204 |
BCR | 0 | 0 | 1.17 | 3.97 | 1.02 | 3.54 | 1.46 | 3.70 |
Scenario | MONSGA | MOPSO | Proposed MOGWO |
No Temperature ($ T_{No} $) | 243.65 | 238.26 | 221.34 |
With Temperature ($ T_{With} $) | 296.52 | 289.28 | 269.98 |