
==== Front
Sensors (Basel)
Sensors (Basel)
sensors
Sensors (Basel, Switzerland)
1424-8220
MDPI

10.3390/s21155246
sensors-21-05246
Article
A New Hybrid Algorithm for Multi-Objective Reactive Power Planning via FACTS Devices and Renewable Wind Resources
https://orcid.org/0000-0003-2232-7189
Syah Rahmad 1
Khorshidian Mianaei Peyman 2
https://orcid.org/0000-0003-2014-6173
Elveny Marischa 3*
Ahmadian Naeim 4
Ramdan Dadan 1
Habibifar Reza 5
https://orcid.org/0000-0002-3697-1808
Davarpanah Afshin 6*
Cano-Ortega Antonio Academic Editor
Sánchez-Sutil Francisco Academic Editor
Laudani Antonino Academic Editor
1 Data Science & Computational Intelligence Research Group, Universitas Medan Area, Medan 20223, Indonesia; rahmadsyah@staff.uma.ac.id (R.S.); dadan@uma.ac.id (D.R.)
2 Department of Electrical and Computer Engineering, Babol Noshirvani University of Technology, Babol 47148-71167, Iran; peyman.khorshidianmianaei@gmail.com
3 Data Science & Computational Intelligence Research Group, Universitas Sumatera Utara, Medan 20155, Indonesia
4 W. P. Carey School of Business, Arizona State University, Tempe, AZ 85281, USA; Naeim.Dian@asu.edu
5 Department of Electrical Engineering, Sharif University of Technology, Tehran 113651-1155, Iran; rhabibifar@alum.sharif.edu
6 Department of Mathematics, Aberystwyth University, Aberystwyth SY23 3BZ, UK
* Correspondence: marischaelveny@usu.ac.id (M.E.); afd6@aber.ac.uk (A.D.)
03 8 2021
8 2021
21 15 524628 5 2021
30 7 2021
© 2021 by the authors.
2021
https://creativecommons.org/licenses/by/4.0/ Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
The power system planning problem considering system loss function, voltage profile function, the cost function of FACTS (flexible alternating current transmission system) devices, and stability function are investigated in this paper. With the growth of electronic technologies, FACTS devices have improved stability and more reliable planning in reactive power (RP) planning. In addition, in modern power systems, renewable resources have an inevitable effect on power system planning. Therefore, wind resources make a complicated problem of planning due to conflicting functions and non-linear constraints. This confliction is the stochastic nature of the cost, loss, and voltage functions that cannot be summarized in function. A multi-objective hybrid algorithm is proposed to solve this problem by considering the linear and non-linear constraints that combine particle swarm optimization (PSO) and the virus colony search (VCS). VCS is a new optimization method based on viruses’ search function to destroy host cells and cause the penetration of the best virus into a cell for reproduction. In the proposed model, the PSO is used to enhance local and global search. In addition, the non-dominated sort of the Pareto criterion is used to sort the data. The optimization results on different scenarios reveal that the combined method of the proposed hybrid algorithm can improve the parameters such as convergence time, index of voltage stability, and absolute magnitude of voltage deviation, and this method can reduce the total transmission line losses. In addition, the presence of wind resources has a positive effect on the mentioned issue.

multi-objective optimization
reactive power (RP) planning
hybrid algorithm
virus colony search
particle swarm optimization
==== Body
1. Introduction

Nowadays, the importance of reactive power (RP) supply is more evident. Even though the RP load flow affects the line losses and the bus voltage and the active power load flow and its costs, unfortunately, the costs of producing and transmitting RP have received little attention [1,2,3,4,5,6,7,8,9,10,11,12]. One reason for this lack of attention is the inherent difficulty in understanding this issue, particularly by economists [12,13,14,15,16,17,18,19,20]. Another reason for the lack of serious attention to RP is the low production price compared to active power [21,22,23,24,25,26,27,28,29,30]. However, according to economics and market calculations, RP is not less valuable than active power [31,32,33,34,35,36,37,38,39,40]. The management and control of RP both in traditional and competitive systems have been one of the main concerns of exploiters [41,42,43,44,45,46,47,48,49,50]. RP in systems was optimized via adjusting the voltage of generators, tap transformers capable of changing tap under load, and the size of parallel capacitors. Due to the expansion of power networks, the exploitation of RP resources for increasing voltage profile and reduction has attracted considerable attention. Although the production of RP does not incur any costs itself, it influences the total cost due to its impact on system losses. In general, RP planning involves the two issues of locating and operating [51,52,53,54,55,56,57,58,59,60]. The issue of locating includes determining the type of RP resources, while the issue of operating emphasizes the optimal arrangement of available RP resources. Parameters that need to ordered RP flow problems are the voltage of controlled buses, tap transformers with the capability to change tap under load, and the size of parallel capacitors. While solving RP flow, it is assumed that economic load flow has been completed, and the amount of real power produced by each generator is determined [61,62,63,64,65,66,67,68]. Therefore, RP planning can be expressed in ways that have functions and constraints [69,70,71,72].

For solving this problem, various objective functions can be considered. RP planning is considered minimizing capital expenditure (CAPEX) and operating expense (OPEX) because their constraints are fulfilled in most papers. CAPEX includes the costs of purchasing, setting up, and conserving novel RP resources in the network. OPEX is considered, and the aim is to minimize these costs [73]. An objective function was considered the primary function and other functions as the constraints, and optimization was conducted to solve this problem. It is challenging to implement such a method to solve multi-objective non-differentiable constraint optimization problems where the solution space is non-convex. In [74], the non-convex and discontinuous constraint in the simplex algorithm is eliminated by implementing the Big M method, which models this disjunctive constraint through a mixed-integer formulation.

On the other hand, due to greenhouse gases and ozone layer depletion, renewable energies such as wind power are currently developing rapidly in most countries. The notion of expansion is synonymous with protection concerning the natural environment and is also discussed in GDP and the use of natural and environmental resources [75,76,77]. The imminent end of fossils has urged researchers to discover renewable energies for substituting the standard energy system [78]. RP control has become more complicated under the influence of renewable resources.

The reasons mentioned above led power industry engineers to revise the existing power systems’ design and present the design of flexible AC transmission systems (FACTS). This design, which involves the use of power electronics devices to increase the transmission capacity and controllability of power systems, has established a new approach to current dynamic problems and control of power systems and has introduced new solutions in this field along with the advancements in power electronics technology and the introduction of the concept of flexible AC transmission systems called FACTS [79]. Planners have found it interesting that concerning RP, FACTS technology provides new chances to control and improve the operational capacity of systems and the new and upgraded lines.

The possibility of controlling the flow within a transmission line at a reasonable cost makes it possible to increase the capacity of the existing lines in the form of lines with more prominent conductors and to increase power flow through such lines under normal and unpredictable conditions by using one of the FACTS controllers [78]. In RP planning, various methods have been proposed that can be divided into mathematics and optimization intelligence.

The first category is based on mathematical models such as non-linear, linear, and quadratic programming [80,81,82]. Unfortunately, the methods described above lose their effectiveness despite the success achieved when the system has many variables in the non-linear structure.

The second category is based on optimization intelligence, which is generally faster and more reliable than the mathematical methods. The volume of computations in this method is less than that of the classical methods. In this category, the following references can be referred to. Shaheen et al. [81] used a modified multi-agent evolutionary model of the RP planning problem. The weakest buses are selected for the installation of RP. Then they are solved using the evolutionary algorithm and target functions based on loss costs and losses. The proposed method in this reference is applied in various standard systems, and the results obtained are compared using genetic and harmonic search algorithms. Bin et al. [82] proposed a multi-objective group search method to accommodate capacitors with different target functions optimally. The voltage drop in lines and the cost of capacitor installation are selected as the target functions. In this paper, the capacitor is used ideally and continuously, and, therefore, the result obtained cannot be applied directly and without changes.

Bhattacharyy et al. [83] addressed the teacher–student learning method for solving RP planning problems by considering linear and non-linear constraints. The proposed method in this reference is applied in a variety of standard systems. Shaheen et al. [84] provide a comprehensive overview of all available methods in the first and the second categories. In this reference, the weaknesses and strengths of the available methods are described meticulously. In [85,86], the gravitational search method is employed to RP problems with FACTS devices presence. Xian et al. [87] modeled the effect of energy on economic distribution for power and pollution. In this model, probabilistic methods are employed for wind energy modeling. El-sobky et al.’s [88] method based on a reliable search area was used for solving the dynamic distribution problem of pollution. Arul et al. [89] used the self-adaptive harmony search algorithm for solving the problem. Control coefficients are focused in a self-regulated manner with program repetitions. The chaos theory based on logistics has also been utilized to improve local search. Despite the methods mentioned earlier, due to the non-linear, discrete, and continuous nature of the constraints in the RP programming problem, a hybrid model that can cover this issue more accurately is required. For overcoming this weakness, the virus colony search optimization algorithm is presented [90]. Because voltage limits and variations are introduced as a constraint, the bus voltage can decrease to a maximum, which increases the probability of voltage instability in those buses. This issue is fundamental in power systems under stress that has been brought on in recent years due to economic issues. Therefore, in this paper, two other target functions, namely, increasing the system security and minimizing the voltage bias of the bus, are considered, and the problem is taken into account as multi-objective programming. The proposed method, coded as a multi-objective method, is one of the highlights that has never been presented in studies thus far.

The model efficiency was investigated on sample systems with FACTS resources formed on proper criteria, and the obtained data were compared to another model. According to the explanations mentioned earlier, the following points can be made about the innovation and highlights of this paper:(A). Improving local and general search via algorithm form on virus search, particle swarm optimization. This hybrid algorithm tries to employ its powerful searching operators in the optimization problem. In addition, it develops the standard particle swarm optimization with the best guiding during the search. Furthermore, it presents a fuzzy mechanism to select the best solution from several solutions.

(B). Modeling the practical constraints in RP programming with continuous and discrete variables in an optimization problem and considering the wind units and FACTS devices to make an accurate model of power system

(C). Consider several analyses and scenarios to evaluate the proposed model and optimization algorithm. In addition, present some comparisons with another optimization algorithm.

(D). The second section deals with modeling RP distribution by considering renewable sources and FACTS devices. The third section is devoted to modeling the proposed multi-objective hybrid algorithm. The fourth section describes how to implement the algorithm on the RP programming, and the fifth section presents simulation data. The sixth section was dedicated to a conclusion.

2. Materials and Methods

Due to the increase in electric power consumption over the past decades, electric power supply systems have been expanded to such an extent that today the optimal distribution of RP for optimal programming and utilization of power systems between energy production units with the lowest cost has become widespread and complex in power system issues. The formulation of the problem of RP programming by considering linear and non-linear constraints is defined as follows.

2.1. RP with FACTS without Wind Unit

In this paper, to achieve more comprehensive RP programming, different target functions have been applied according to the following descriptions.

2.1.1. Loss Function

It is an essential function in RP programming and aims to reduce network losses, increasing network efficiency. To calculate network losses, the Newton–Raphson method has been employed [48]:(1) PL=∑k=1NLgk[Vi2+Vj2−2ViVjcos(θi−θj)

Vi and Vj represent voltage for the beginning and last of lines i and j, and the θi and θj show the angles at the beginning and end of the line, respectively.

2.1.2. Voltage Equalization

Common issues in RP programming are network voltage stability. As the nature of the network voltage is different from that of the loss function, a new function called network voltage equalization, which indicates difference for desired per bus and obtained after each run of the load distribution program, is as follows:(2) VIndex=1N∑i=1N|Vact−Vdes|

In the above relation, N is the candidate bus number for installation in RP sources, Vact is the value obtained from the load distribution problem with RP sources, and Vdes is the expected value according to the constraints on the power system network.

2.1.3. The Cost of FACTS Devices

The cost of equipment performance consists of two parts: the cost of the active power and the cost including shunt equipment. The target is obtained via:(3) CTotal=CE+CSh

CE is the cost of losses ($), and CSh is the cost of shunt capacitors ($). CE generally increases by considering system losses. In this paper, installing shunt equipment cost, installing equipment cost, and the cost of energy losses have been considered to be equal to 3 $/kVar, 1000 $, and 0.06 $/KWh, respectively. Now, the target function for RP programming is expressed by considering the FACTS equipment as follows:(4) CFACTs=CTCSC+CSVC

where CTCSC and CSVC show the cost of equipment performance of TCSC and SVC in RP programming, respectively, as follows [17]:(5) CTCSC=0.0015(TCSCvalue)2−0.7130(TCSCvalue)+153.75(US$/kvar)

(6) CSVC=0.0003(SVCvalue)2−0.3051(SVCvalue)+127.38(US$/kvar)

2.1.4. The Objective Function

As seen above, the final target function is presented separately for each section. The target should be considered as consisting of all of the above points. The final target function proposed in this paper is based on the value of losses, voltage fluctuations, and the assignment of FACTS devices and the related constraints as follows:(7) Minimizeα∑k=1,…,NLPloss(k)+β∑i=1,…,NQciyi+γ∑i=1Nd|Vi−Vi0|subject.to.∑i=1Nyi=NcPGi(k)−PLi(k)−Vi(k)∑jwiVj(k)×(Gijcosθij(k)+Bijsinθij(k))=0QGi(k)+QCi(k)(yi)−QLi(k)−Vi(k)∑jwiVj(k)×(Gijsinθij(k)−Bijcosθij(k))=0Vimin(k)≤Vi(k)≤Vimax(k)PGimin(k)≤PGi(k)≤PGimax(k),QCimin(k)≤QCi(k)≤QCimax(k),i∈NG−Sijmax(k)≤Sij(k)≤−Sijmax(k),Ptieline(k)≤TTC(QC(k))

where NG, N, NL, and Nc represent generators, buses, load levels, and Var resource numbers. k shows load levels, and k, T(k), α, β, and yi represent kth load level, cost of energy, coefficient of the maintenance cost, cost of allocating the equipment, and the binary variable, respectively. If yi = 1, the Var equipment in that busbar is installed; otherwise, yi = 0. PGi, QGi, PLi, and QLi represent the output of the active power of generators, the output of the RP, the active load, and load in request in ith, respectively. Sijk is the current flowing from the i-j line. When Sijk is not equal to zero, then Qci will be a specific capacity of Var for the bus i, and Qcik is the value of Var required at k. Ptieline is power flowing with communication lines between them. TTC(Qc(k)) is a one-way interpolation following Formula (4) for voltage stability. Qc is the total set of Qci candidates. (8) p(x,y)=∑i=0m∑j=0nf(xi,yi)∏l≠i(x−xlxi−xl)∏k≠j(y−ykyj−yk)

The y and x were Var capacities, p(x,y), TTC, and Var function capacity. Equation (7) addresses the minimization of cost, voltage stability, and system losses. As stated previously, one of the goals of the present paper is to ensure the system voltage stability during optimal load distribution. The most crucial issue is voltage stability for both standard and critical (unwanted) conditions. For an extensive system with a probability of various errors, these two modes are a challenge. In order to address this problem, the following formulation are able to guarantee system stability under different working conditions:(9) mink{max∑i∈Source Areaj∈Sink Area(Pij(k)−Pij0)}subject.to.∑i=1Nyi=NcPGi(k)−PLi(k)−Ui(k)∑jwiUj(k)×(Gijcosθij(k)+Bijsinθij(k))=0QGi(k)+QCi(k)(yi)−QLi(k)−Ui(k)∑jwiUj(k)×(Gijsinθij(k)−Bijcosθij(k))=0Uimin(k)≤Ui(k)≤Uimax(k)PGimin(k)≤PGi(k)≤PGimax(k),QCimin(k)≤QCi(k)≤QCimax(k),i∈Source Area−Sijmax(k)≤Sij(k)≤−Sijmax(k),PLi(k)PLi0=QLi(k)QLi0,i∈Sink Area

where Q0Li and P0Li are the initial demand for active RP in ith busbar, respectively. Pij0, Pij(k) were the initial power and the current power flowing in the line between the i-j busbars. K = 1, …, Ncntg shows the different working conditions, including normal conditions and those containing the error.

2.1.5. Fuzzy Classification

One of the essential points in RP was investigated optimal location for FACTS devices. In other words, the most crucial point in RP was determining the proper location of the Var equipment in order to reduce calculation time and ensure voltage stability. In this paper, interactive fuzzy classification is proposed to obtain FACTS location devices. Fuzzy classification combines standard fuzzy classification and regression [49]. In the proposed fuzzy method, the distance (dij) between each data (xi) and the center of the class (v) is expressed based on the standard fuzzy distance and regression:(10) d2((xj,yj),(vi,βi))=||xi−vi||⏟fuzzy distance+(yj−βiTxj)2⏟fuzzy−regression distance

where xj∈X and yj∈Y are input and output data, respectively. nd in the index j = 1, …, nd is the data number in a vector for training.

c in the index i = 1, …, c is the number of centers of the classes. y^j=βiTx is the polynomial estimation, where β is a polynomial coefficient. Based on the above definitions, the membership function uj(xi) can be expressed as follows:(11) uik=(∑j=1cdik2−(mindik2−η)djk2−(mindjk2−η))−1,η>0

where η>0 is considered as a constant value. In order to optimize fuzzy members, the error of the regression is considered in the following objective function:(12) JmMFC=∑i=1c∑k=1ndμikmdik2⏟FCM+∑i=1c∑j=1ndμikm(yk−fi(τik))2⏟error of fuzzy regression

where μik is the undefined function from the kth input vector in the ith fuzzy category. m is the degree of fuzzification that expresses the amount of overlap between the categories. As the abovementioned relationship shows, JmMFC consider two parts. The first part, which is similar to the standard fuzzy classification mechanism (FCM), controls the categorization accuracy. The second part is the sum of the error squares (f(τi)) derived from the ith category. This section checks the error rate due to the fuzzy classification based on the fuzzy function, the matrix (τ), and the “or” and “and” constraints. Based on the new fuzzy distance defined in the above relation, we can introduce the new fuzzy membership function as follows:(13) μik1<i<c1<k<nd=(∑j=1c[dik2+(yk−fi(τik))2dik2+(yk−fi(τjk))2]1m−1)−1

In the new membership function, the centrality of each category can be expressed as follows:(14) ∀1<i<cvi=∑k=1nd(μik)m(xkyk)∑k=1n(μik)m

The above formula states that the membership function for xi depends on the distance between that variable and all vi categories. The categories created in this paper can be defined based on the loss factor, voltage stability, and the cost of installing FACTS equipment.

2.2. RP with FACTS and Wind Unit

Rapid progress in wind power generation has led to a significant increase in the installation of wind turbines in the world. Research centers worldwide are extensively studying and researching various aspects of the optimal utilization of this energy resource. Simulation of wind turbines by software is a research tool for investigating the behavior of wind turbines. In addition, increasing the simulation parameters and considering the system details increase the accuracy of the simulations and the simulation time. Therefore, various dynamic models of generators for power generation from wind have been studied, one of which is electricity generation from the double-fed induction generator (DFIG) [6]. DFIG can be a good option for systems with variable speed (approximately 30% synchronous speed). Figure 1 shows the schema of a DFIG [91].

In a DFIG, the injection of RP for the network depends on the control system and the size of the converter. Figure 2 shows the P-Q characteristic for wind turbines (WTs) used in this paper. Specifications of this turbine (Gamesa WT G80-2.0MW) can be obtained from the reference [92]. This turbine provides a 0.98 power factor in the pre-phase mode and 0.96 in the post-phase mode. Therefore, its reactivity depends on the generation of active power, as shown in Figure 3.

The green line in Figure 3 shows the P-Q of the wind turbine (WT) at the power factor 1. For power less than 10 megawatts, the cables have a more significant effect than the transformer, while the transformer is more influential in the high-power WT. In addition, if the wind farm (WF) receives the RP of a capacitor lower than the range of the active power generation, then WTs are set to the selfie mode. Based on the above explanations, since WFs can inject or receive RP, we need to adjust them optimally. In order to optimize RP, STATCOM or capacitor is used. The ultimate target function, by taking into account the loss function, is as follows:(15) Minimize J(Varx,Vary)=Min Plosses

where Vary shows the location of the tapped transformer. The optimal configuration of STATCOM and capacitor is performed independently, which includes the following constraints:(16) QWTimin≤QWTi≤QWTimax,i=1,2,…,NG

(17) Timin≤Ti≤Timax

(18) QStatcommin≤QStatcom≤QStatcommax

There is also a non-linear constraint. RP demand for PCC to adjust voltage is modeled as follows:(19) QPCC∗=QPCCmeas

In this paper, the condition of obtaining the possible response is given in Equation (20) to guarantee that the set of solutions covers the limits imposed on the system:(20) Sik+1={Sik+vik+1,Simin≤Sik+vik+1≤SimaxSimin         ,Simin>Sik+vik+1Simax         ,Sik+vik+1>Simax

Since Constraint (19) is an equality constraint, to apply it non-linearly, we can use the following equation based on error (ε):(21) QPCC*−QPCCmeas|<ε

Based on the explanations mentioned earlier, Figure 4 shows the strategy used in this paper.

3. The Hybrid Algorithm of Virus Colony Search and Modified Particle Swarm Optimization

3.1. Virus Colony Search Algorithm

Definitions and Basic Concepts

VCS algorithm form on viruses population, and it was founded for two virus behaviors. Viruses require host cells to ensure their survival, production, and plurality. In what follows, the strategy of viruses has been described, including host cell contamination [22].

The function of the immune system is also considered as the opposite relationship between these two stages. These stages are explained in more detail as follows:(1) Virus replication: Viruses randomly scan the cells to find the raw materials needed for survival. In this process, a random walk method can be the best mathematical model for expressing the function of viruses to detect host cells.

(2) The infection or influence of host cells: When a virus selects a cell, it tries to replicate itself in the best manner. In other words, the virus replicates itself based on the essential materials in the host cell until the host cell dies and acts like a virus. To model this behavior, the CMA-ES mathematical model is a covariance matrix-based method.

(3) Function of the immune system: The immune system has the essential task of protecting the cell against the replication and spread of the virus. In addition, more powerful viruses save themselves from reproducing in the proper position.

Matching mathematical models

Virus algorithm has three models: (1) Gaussian model to virus propagation that utilized for recognition improving; (2) CMA-ES model to cell damage modeling that utilized for social behavior improving; and (3) evolution strategy to resist a virus that utilized for improving two other strategies [93]. Therefore, the virus colony algorithms have five steps:Vpop virus and Hpop cells were utilized in the VCS model.

Virus for transferring has unique behavior.

Viruses can desire to infect the cell.

Virus reproduction was formed on cell destruction to find survival.

Best viruses can remain for replication.

Virus transmission

Gaussian model is formulated via this equation:(22) Vpopi’=Gaussian(Gbestg,τ)+(r1.Gbestg−r2.Vpopi)

The ɩ standard deviation was calculated via log(g)/g.(Vpopi−Gbestg)the log(g)/g is used to increase virus search. The (Gbestg,τ) is produced the best answer for virus searching via vector (r1.Gbestg−r2.Vpopi).

The cell influence

The cell was infected, and it can be modeled via CMA-ES in these steps:

Step 1: Update the Hpop via:(23) Hpopig=Xmeang+σig×Ni(0,Cg)

where σg > 0. Xmeang was obtained via this equation:(24) Xmean0=1N∑i=1NVpopi

Step 2: Select vector γ and obtain parents vector via this equation:(25) Xmeang+1=1γ∑i=1γωi.Vpopiγbest|ωi=ln(γ+1)/(∑j=1γ(ln(γ+1)−ln(j))

In the above relation, γ=[N/2], wi is the coefficient of the combination, and two mentioned paths are obtained via:(26) Pσg+1=(1−Cσ)Pσg+Ca(2−Cσ)γω1σg(Cg)−12(Xmeang+1−Xmeang)

(27) PCg+1=(1−Cc)PCg+hσCC(2−CC)γω1σg(Xmeang+1−Xmeang)

The parameters were set via Cσ=(γω+2)/(N+γω+3), s.

Step 3: Update σg+1 and Cg+1 via:(28) σg+1=σg×exp(Cσdσ(||Pσg+1||E||N(0,1)||−1)) 

(29) Cg+1=(1−C1−Cγ)Cg+C1PCg+1(PCg+1)T+Cγ∑i=1γwi(Vpopiγbest−Xmeang)σg.(Vpopiγbest−Xmeang)Tσg

where Cγ is acted via this equation:(30) C1=1γw((1−1γw)min{1,2γw−1(N+2)2+γw}+1γw2(N+2)2),Cγ=(γw−1)C1

Immune system

The evolutionary movement is obtained via these steps:

Step 1: Pr to Vpop population:(31) Prrank(i)=N−i+1N

Step 2: Each individual grows individually to Vpop via this equation:(32) {Vpopi,j"=Vpopk,j−rand.(Vpoph,j−Vpopi,j)                                                       if r>Prrank(i)Vpopi,j"=Vpopi,j                                                       otherwise

The response is tried to save the best stage value, and other unusual response can be obtained via this equation:(33) xij=lowij+rand×(upij−lowij)

3.2. Modified Particle Swarm Optimization

Eberhart has introduced a swarm optimization model by the group search for food by birds or fish [94]. In the past decade, a variety of topologies has been proposed for the PSO algorithm. Topologies used to exchange information between particles in the PSO algorithm include star, ring, and square topologies [95]. Figure 5 shows the relationship between particles in different topologies.

In the star topology for D-dimensional space, the best personal status for particle i is indicated by p→i=(pi1,pi2,…,piD) and the group’s best position g→=(g1,g2,…,gD). The relations between the velocity and movement of particle i in the next moment or repetition are obtained by the following relations [29]. (34) vid(t+1)=ωvid(t)+c1rand1(pid(t)−xid(t))+c2rand2(gd(t)−xid(t))

(35) x→(t+1)=x→(t)+v→(t+1)

c1 and c2 are called the cognitive and the social parameters, respectively, which are usually set to 2. To randomize the nature of velocity, the coefficients c1 and c2 are multiplied by the random numbers rand1 and rand2. Usually, in the implementation of the PSO, the value of ω decreases linearly from one to zero. Generally, the inertia coefficient ω is governed by the following equation. (36) ω=ωmax−ωmax−ωminitermax.iter

The value of vmax is usually chosen in such a manner that vmax = kxmax and 0.1 < k < 1. xmax specifies the length of the search. As is clear from Formula (34), the coefficients c1 and c2 are usually considered constant, which is the dark point in the local and final search for the particle swarm model. These coefficients can be suggested to improve the efficiency of virus search:(37) vid(t+1)=ωvid(t)+c1irand1(pid(t)−xid(t))+c2irand2(gd(t)−xid(t))

The correlation coefficient c1i is updated in each replication. If c1i has a small value, then c1ir1 is small and local search is strengthened. Conversely, it has an immense value if it has a significant value, thus improving the overall search. To select the best value for c1i, two thresholds of T1 < 0 and T2 > 0 and two variables of R1 in the range (0, T1) and R2 in the range (0, T2) are used and defined as c1i=2R1 and c2i=2R2. As a result, the two vector populations with coefficients of c1i and c2i are generated. When T1 < 0, then R1 is negative, and c1i has a small value, thus improving local search.

3.3. The Proposed Hybrid Algorithm

This section introduces the structure of the hybrid modified particle swarm optimization and HMPSO-VCS model. Figure 6 shows a schema of the HMPSO-VCS. It also shows that the hybrid VCS and particle swarm optimization algorithm is started for an initial population. When N-dimensional problem is obtained, then the hybrid algorithm has 4N members generated entirely randomly. Next, 4N members are sorted according to their eligibility, and the upper 2N members enter the VCS as a virus. VCS formulations create the new 2N member population. The mechanism of optimizing the particle swarm is applied to the 2N lower members as particles. In applying the mechanism of particle swarm optimization, the population is obtained via VCS. The Pgbest is the corresponding member is used in Equation (37). The population generated by applying the mechanism of the particle swarm optimization to the population obtained via virus colony search algorithm was merged, and the 4N new member population is sorted according to the eligibility, and the previous process is repeated until convergence is achieved.

The steps for the HMPSO-VCS algorithm are as follows:(1) Random generation of the initial population with 4N members as initial responses

(2) Evaluating and sorting the population according to its eligibility

(3) Applying the VCS algorithm to the 2N upper members of the population based on the mutation and crossover of the generations Selection: For the target population, the best 2N members are selected based on their eligibility.

Crossover: For a well-chosen population, we use the crossover of two parents to produce a new generation.

Mutations: 20% of the new population is mutated.

(4) The particle swarm optimization algorithm is applied to another 2N population based on the production relations of the new population, and the new population is produced. Next, 2N population is combined with the 2N population generated by the virus colony search algorithm.

(5) Repeat the previous steps from Step (2) until the convergence or termination requirements are met.

3.4. The Proposed Multi-Objective Algorithm

The proposed algorithm has been used based on the entropy accumulation to label each objective function. Entropy optimization is based on the Pareto criterion. In this method, for optimization in each step, a set is labeled as a better response and enters the next step. The criterion for the multi-objective problem can be obtained via:

The X1 was no worse than X2, and X1 were better than X2, and mathematical terms are stated via:(38) X1≺X2⇔(∀i∈{1,2,…,n}:fi(X1)≤fi(X2))^(∃i∈{1,2,…,n}:fi(X1)≤fi(X2))

In addition, the X∈Xf is stated non-dominated of relation to A⊆Xf if (39) ∃a∈A:X≺a

Therefore, we can assume that the vector X is optimal [96]. The non-dominated decision vectors sets is obtained via:(40) P(A)={a∈A|aisNon−dominatedA}

Set A was equal to Xf, then P(A) is produced Pareto optimal front. An ideal point and anti-ideal point, objective functions are minimal and maximal, respectively.

3.5. Combining Fuzzy Logic with the Proposed Algorithm

Figure 7 indicates membership μc [97] for fuzzy that states cost concept.

The fi(Pgi) is stated via μi(Pgi) in this equation:(41) μi=fimax−fifimax−fimin

The lower limits of ith objective function:(42) FDMi={0μi1μi≤00<μi<1μi≥1

The FDMk function is obtained via:(43) FDMk=∑i=12FDMik∑j=1M∑i=12FDMij

The decision maker is selected the best solution by considering the function in actual conditions.

4. Implementation of Problem of RP Planning

To solve RP with FACTS presence, the following steps are performed:

Step 1: The set initial parameters of the suggested algorithm are studied. The initial population with N units is produced as follows:(44) P=[Var1,Var2,…,Varn  Q1,Q2,Q3,…,Qm,T1,T2,…,To]

In this regard, the number of FACTS devices, tap transformers, and shunt capacitors is encoded as mentioned above. Subnets n, m, and o represent the number of variables used. Variables are also generated based on the upper and lower bounds:(45) Ti=Ti,min+rand×(Ti,max−Ti,min),i=1,2,3,…,oQi=Qi,min+rand×(Qi,max−Qi,min),i=1,2,3,…,mVari=Vari,min+rand×(Vari,max−Vari,min),i=1,2,3,…,n

The initial responses should be placed within the constraints to create a better search for the algorithm at this stage.

Step 2: Apply the limitations on the problem, arrange them on the cost function.

Step 3: Use the proposed fuzzy cluster algorithm to select the best load bus bars for installing the FACTS devices.

Step 4: Choose the best response, and these sets of responses are replaced in initial responses that are considered null.

Step 5: Use the target function introduced for the loss function, voltage, and cost of FACTs installation, calculate the fitness of responses and select the best response.

Step 6: Compare good response and stored response.

Step 7: Sort data based on the Pareto criterion and choose the best response sets as the best current response.

Step 8: Replace the worst responses by random data and use a modified particle swarm to find the best response.

Step 9: Is termination completed? In this case, leave the program; if not, go to Step 3.

Figure 8 shows the flowchart of the proposed algorithm.

5. Results and Simulation Analysis

5.1. Determination of Optimal Location by the Proposed Fuzzy Method

First, we propose a fuzzy method for obtaining FACTS installation. The purpose of this action is to determine the most talented buses for installing these resources. To determine these locations, criteria such as L-index, U/U0, and voltage stability have been used. (A). U/U0 criterion: This criterion is considered based on the voltage range for all buses as follows:

(46) UU0=[U1U01,…,UkU0k]

The U and U0 are the voltage range vector of the normal state for all buses and the voltage range vector of all bass when all loads are zero. (B). L-Index criterion: Generally, increasing load and the optimal performance of the power system are considered more than before. A voltage failure can be created suddenly in the system, resulting from weaker buses. For the jth bus, the rate of voltage drop or failure is expressed as follows:

(47) {Lj=|1−∑i=1NPVFjiViVj−1,j=1,2,…,NPQFji=−[Y1]−1[Y2]

Y1 and Y2 matrices were below the invertible submatrices for YBUS matrix, which are obtained by sorting them according to the following matrix:(48) [IPQIPV]=[Y1Y2Y3Y4][VPQVPV]

(C). Voltage stability criterion is expressed in Formula (2).

5.2. IEEE 30-Bus Standard

The 30-bus system is the most popular in power engineering. This system consists of 6 power plants and 46 transmission lines. Table 1 presents the results of applying the fuzzy method and determines the most talented bus for SVC. As the numerical results show, according to the L-Index, the most talented buses are {3, 7, 9, 14, 15, 17, 23, and 24}. According to the U/U0 criterion, these talented buses are {4, 7, 15, 17, 20, 21, and 25}. Finally, based on the voltage criterion, the best buses for installation are {4, 7, 15, 16, 17, 21, and 24}. Based on the criteria described, the candidate buses are not the same in all three modes. By fuzzy classification, the best buses are classified as follows: {7, 15, 17, 21}, {4, 24}, and {9, 16}. Similarly, fuzzy classification is performed for TCSC, and the best candidates are lines 5, 25, 28, and 41. Now, we should consider the RP planning with the presence FACTS devices in the system under study. Table 2, Table 3 and Table 4 shows these results with the presence of FACTS devices. As mentioned previously, two types of SVC equipment have been used as shunt and TCSC series in line with the present study. Table 3 indicates a comparison of the costs and losses of the 30-bus standard system in RP planning. According to this table, it is clear that when FACTS equipment operates in sync with RP devices such as tap transformers and RP generators, the system operation costs and losses considerably improve compared with when FACTS devices are not present. In addition, the results reveal that the mentioned model has better performance than another models.

As can be seen from the results of Table 2, Table 3 and Table 4, the proposed method has been able to obtain far fewer answers compared to other methods in terms of cost criteria, loss rate, and operating cost. On average, it can be seen that the proposed method worked about 5% better.

5.3. IEEE 57-Bus Standard System

Here, to illustrate the proposed model, a more extensive 57-bus system is selected with the information mentioned in other papers. Similar to the previous section, in this section, the fuzzy classification method is used first. Based on this method, buses 25, 38, and 49 are the best option for SVC installation, and lines 13, 37, 57, and 61 are the best candidates for TCSC installation. Table 5 shows the comparison between system costs and losses of the 57-bus standard system in RP planning. The results indicate mentioned model is more successful than other models. This represents the better local and overall search of this algorithm to find the final result. This algorithm has shown better performance despite the increasing number of input variables.

5.4. RP Planning for Wind Turbines (WTs)

RP planning performance in the Mongolia system is studied. There are approximately 200 wind turbines installed in this system. Capacity is over 300 megawatts. With the help of T2 transmissions, the output of these turbines is increased to 35 kV and is transmitted through the transmission lines to the distribution site. This wind farm contains 403 nodes, and more information is available in reference [99]. Figure 9 shows a section of this network. To compensate for RP, the SVC compensator is installed in the nodes. The reference power in this system is 100 megawatts, and the most significant revenue generated from RP is approximately $500 million. Table 6 provides the comparison between the performance of the mentioned model and other models for RP planning.

As the results presented in the following tables demonstrate, the proposed method has produced significant improvements in responses to the problem. It is also clear from Table 6 that when the wind speed changes from 4 to 12 m/s, the loss rate increases significantly, but compared to other optimization methods, the proposed method has fewer losses.

5.5. RP Planning Based on the Possibility of Controlling Wind Turbines (WTs)

As the last analysis of the system under study, this section investigates RP planning by fully covering Section 2.2. RP planning is applied to the wind farm (WF) for this model, according to Figure 10. The short circuit current for the external network is 400 megawatts. In the wind farm, 12 connected wind turbines are considered.

In addition, in this system, the STATCOM compensator is used at PCC. The objective in the scenario is to reduce the losses and keep the system error in the safe range based on setting the PCC setting point. Usually, such a function is considered for the demand for a large amount of RP. In addition, the tapped transformer and RP compensator are used as other control items. In order to control RP, six different control models are applied to control RP in the wind farm. Figure 11 shows the structure of these control systems. The maximum RP error based on the PCC adjustment system is approximately 5%. The strategies adopted are as follows:

Strategy 1: The variables are RP for wind turbine QWTi, location of the tapped transformer for PCC, and capacitor banks. The capacitor capacity was 0.96 MVAr. Figure 11a shows the structure of this strategy.

Strategy 2: RP of wind turbine QWTi is the only control parameter. Power production maximum, amount of RP was 4 MVAr. Figure 11b shows the structure of this strategy.

Strategy 3: In this case, RP planning was formed on the RP amount of wind turbine QWTi injected into the capacitor system and bank. Power production maximum, amount of RP was 4 MVAr and Q*PCC = 4 MVAr. Figure 11c shows the structure of this strategy.

Strategy 4: In this case, RP planning is based on adjusting the tapped transformer and the generation of RP from the wind turbine. Therefore, two control parameters in this system are QWTi and tap transformer. RP value is equal to Q*PCC = 4 MVAr. Figure 11d shows the structure of this strategy.

Strategy 5: STATCOM compensator was installed for PCC. Therefore, in addition to the two previous variables, the STATCOM compensator is considered a new control parameter. Figure 11e shows the structure of this strategy. The STATCOM system with 30% of the total power and wind turbine capacity is selected. Therefore, the capacity of RP with the STATCOM compensator is approximately 1.22 MVAr.

Strategy 6 is similar to Strategy 5, but the difference is that the constraint Q*PCC is not considered. Virus search is employed, and Figure 11 shows the structure of this strategy.

Table 7 and Table 8 show RP data by the proposed algorithm based on the possible response under different working conditions and the corresponding compensation. It is found that different individual adjustment points are considered for WT because the distance and feature of PCC vary from one wind turbine to another.

Table 9 represents the simulation results for strategies 2 to 6 (Figure 11). It is clear that lower power is obtained via STATCOM, and error is lower with Strategy 5 using the STATCOM and the proposed HMPSO-VCS. In Strategies 5 and 6, the results depict that RP error is more significant than PSO than the proposed HMPSO-VCS method.

In the following, the proposed algorithms with other optimization methods are compared. For this purpose, Table 10 reports the list of different test functions which used in this simulation and Table 11 shows the results of the simulation. The proposed method is compared with many other methods and results shown that the proposed algorithm has worked much better than other methods [101,102,103,104,105,106,107].

6. Conclusions

RP distribution in power systems is optimized by adjusting the voltage of generators, tap transformers with the ability to change under load, and the size of parallel capacitors. With the expansion of power networks and the development of renewable resources, considerable attention is devoted to RP resources in the network. With this end in view, objective functions based on losses, voltage, and cost of RP sources were applied in this paper. This paper presents a new model for optimal reactive power distribution. To bring the issue to a more realistic level, renewable resources and facsimiles have also been considered. On the other hand, by increasing the numbers of these types of equipment, if it is not possible to make an optimal and appropriate choice from them, the system’s stability is not adequately established, but it also leads to more errors. For this reason, solving this problem has become an optimization, and a new hybrid algorithm based on virus search and particle clustering has been used to solve it. Using this hybrid model has been able to improve local and global search dramatically. The fuzzy theory was employed to install RP resources of shunt and the SVC and TCSC series. In this model, various criteria were used to determine the best position. In the method used in this scheme, the discrete variables are considered continuous and treated like any other continuous variables in the algorithm’s space. At each stage, after reaching the new position, the discrete variables that have become continuous are now rounded off to the appropriate number before applying to the load program. For example, tapping tap changers is a discrete quantity. In the process of running the algorithm, when a value between the two tap changers is obtained for this quantity, this value is rounded off to the nearest step and then applied to the load program. Finally, the problem of RP planning was solved using a novel algorithm from on virus search. Appropriate search, choosing the best response at each stage, improving the set of responses, and the crossover of generations to give the best responses are among prominent features for the proposed method. The results obtained in different sections show that the proposed method has an average performance of 5% better.

Author Contributions

Conceptualization, R.S. and R.H.; methodology, P.K.M.; software, M.E.; validation, N.A. and D.R.; formal analysis, N.A.; investigation, A.D.; resources, R.H.; data curation, P.K.M.; writing—original draft preparation, D.R.; writing—review and editing, A.D.; visualization, R.H.; supervision, A.D.; project administration, R.S.; funding acquisition, M.E. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Figure 1 Variable speed wind turbine with a DFIG.

Figure 2 Q of WT in G80 (-2.0 MW).

Figure 3 Q of WF in twelve G80 (-2.0MWWTs).

Figure 4 Solution search process.

Figure 5 Topologies used in the PSO algorithm.

Figure 6 Combining the VCS (⟶) and swarm optimization (-→).

Figure 7 The membership function for the cost and loss function and the voltage synchronization.

Figure 8 Applying the proposed algorithm and FACTS presence.

Figure 9 The structure of the system under study with the presence of wind turbines.

Figure 10 Structure of the tested wind farm.

Figure 11 RP strategies, the proposed HMPSO-VCS, and basic VCS.

sensors-21-05246-t001_Table 1 Table 1 The results of determining the most talented bus in SVC installation.

PQ Bus No.	L-Index	U/U0	Voltage Fluctuation (%)	
3	1.2222	1.227	1.1108	
4	1.1307	1.189	1.0667	
6	0.9732	0.9599	1.1627	
7	0.996	1.1638	1.1512	
9	1.2546	1.1074	1.0168	
10	0.9115	1.2892	1.0727	
12	1.096	1.1596	0.9062	
14	0.9672	1.2201	1.2936	
15	1.2915	1.0815	0.9669	
16	1.1851	1.073	0.9425	
17	1.1002	1.2301	1.049	
18	1.0884	0.9334	0.9792	
20	0.9238	0.9533	1.0959	
21	1.1728	0.9694	1.0358	
23	0.917	1.0564	1.2807	
24	0.9286	1.2326	1.2681	
25	1.1087	1.2213	0.9211	
28	0.9387	0.9242	1.1951	
30	1.2273	1.0597	1.0076	

sensors-21-05246-t002_Table 2 Table 2 The optimum range for RP resources for the 30-bus standard system without the presence of FACTS devices.

Method	Generation of RP with Generators (Perunit)	Adjustment of Tap Transformer (Perunit)	Shunt Capacitor (Perunit)	
The proposed method	0.1601 (2)
0.2611 (5)
0.3034 (8)
0.0105 (11)
0.2243 (13)	0.9046 (11)
0.9014 (12)
0.9032 (15)
0.8732 (36)	0.0243 (7)
0.0167 (15)
0.0065 (17)
0.0134 (21)	
SPSO [98]	0.1705 (2)
0.2655 (5)
0.3028 (8)
−0.0121 (11)
0.2467 (13)	0.9 (11)
0.9 (12)
0.9019 (15)
0.9 (36)	0.0330 (7)
0.0527 (15)
0.0 (17)
0.0 (21)	
APSO [98]	0.1629 (2)
0.2671 (5)
0.2964 (8)
0.0563 (11)
0.1797 (13)	0.9 (11)
0.9133 (12)
0.9 (15)
0.9010 (36)	0.0174 (7)
0.0468 (15)
0.0025 (17)
0.0209 (21)	
EPSO [98]	0.1462 (2)
0.2642 (5)
0.2977 (8)
0.1061 (11)
0.2213 (13)	0.9026 (11)
0.9 (12)
0.9 (15)
0.9 (36)	0.0 (7)
0.0 (15)
0.0295 (17)
0.0 (21)	

sensors-21-05246-t003_Table 3 Table 3 The optimal range for RP sources for the 30-bus system with the presence of FACTS devices.

Method	Generation of RP with Generators (Perunit)	Adjustment of Tap Transformer (Perunit)	SVC (Perunit)	TCSC	
The proposed method	0.532 (2)
0.065 (5)
0.032 (8)
0.344 (11)
0.002 (13)	0.843 (11)
0.703 (12)
0.984 (15)
0.898 (36)	0.0 (7)
0.0 (15)
0.0 (17)
0.04 (21)	0.1298 (25)
0.0837 (41)
0.2312 (28)
0.2109 (5)	
SPSO [98]	0.6 (2)
0 (5)
0 (8)
0.4 (11)
0 (13)	0.9 (11)
0.9 (12)
0.9 (15)
0.9223 (36)	0.0 (7)
0.0 (15)
0.0 (17)
0.0840 (21)	0.1463 (25)
0.0419 (41)
0.1049 (28)
0.1388 (5)	
APSO [98]	0 (2)
0 (5)
0 (8)
0.4 (11)
0 (13)	0.9 (11)
0.9501 (12)
0.9180 (15)
0.9330 (36)	0.0 (7)
0.0 (15)
0.0 (17)
0.0768 (21)	0.1463 (25)
0.0419 (41)
0.1049 (28)
0.1388 (5)	
EPSO [98]	0.6 (2)
0 (5)
0 (8)
0.4 (11)
0 (13)	0.9439 (11)
0.9 (12)
0.9 (15)
0.9326 (36)	0.0 (7)
0.0 (15)
0.0 (17)
0.0 (21)	0.1463 (25)
0.0419 (41)
0.1049 (28)
0.1368 (5)	

sensors-21-05246-t004_Table 4 Table 4 Active power losses and operation costs for the 30-bus system with different optimization methods.

FACTS	Method	Active Power Losses	Operation Costs after RP Planning
(×106)	Reduction in Active Power Losses	Reduction in Operation Costs
(×106)	
without the presence of FACTS devices	The proposed method	0.0756	3.403	0.0024	1.5103	
SPSO [98]	0.0684	3.5951	0.0027	1.41916	
APSO [98]	0.0684	3.5966	0.0027	1.40416	
EPSO [98]	0.0685	3.6000	0.0026	1.37016	
with the presence of FACTS devices	The proposed method	0.0431	2.3543	0.0279	1.38283	
SPSO [98]	0.0435	2.3622	0.0276	1.37481	
APSO
[98]	0.0434	2.3558	0.0277	1.38121	
EPSO [98]	0.0438	2.3671	0.0273	1.36991	

sensors-21-05246-t005_Table 5 Table 5 Active power losses and operation costs for the 57-bus system with various optimization methods.

FACTS	Method	Active Power Losses	Operation Costs after RP Planning
(×107)	Reduction in Active Power Losses	Reduction in Operation Costs
(×106)	
without the presence of FACTS devices	The proposed method	0.222	1.413	0.0821	1. 43	
SPSO [98]	0.2522	1.325	0.0277	1.46	
APSO [98]	0.2495	1.311	0.0304	1.60	
EPSO [98]	0.2526	1.327	0.0273	1.44	
with the presence of FACTS devices	The proposed method	0.2205	1.163	0.0798	3.05	
SPSO [98]	0.2210	1.168	0.0589	3.03	
APSO [98]	0.2231	1.179	0.0568	2.92	
EPSO [98]	0.2275	1.203	0.0524	2.68	

sensors-21-05246-t006_Table 6 Table 6 Numerical comparison of the losses of the system under study with the presence of wind resources.

Project	Investment of RP Compensation (Million Yuan)	The System Loss (kW)	
V = 4 m/s	V = 8 m/s	V = 12 m/s	
TGA [100]	338	1872	2480	3129	
IGA [100]	336	1731	2292	2892	
PSO	332	1729	2290	2687	
VCS	327	1723	2288	2676	
The proposed method	311	1698	2261	2623	

sensors-21-05246-t007_Table 7 Table 7 Data of option I.

Q*PCC	PWF
100%	PWF 80%	PWF
50%	PWF 20%	PWF 10%	
4	3.5	3.5	2	3	1	0.5	
QWT1	0.0316	0.0203	0.0232	0.0077	0.0469	0.0089	0.0802	
QWT2	0.0845	0.0154	0.0368	0.0236	0.0208	0.0236	0.0301	
QWT3	0.0788	0.0205	0.0535	0.0721	0.044	0.0302	0.0629	
QWT4	0.0495	0.0392	0.0236	0.0026	0.0562	0.0612	0.0178	
QWT5	0.056	0.028	0.0543	0.0836	0.0611	0.0123	0.0027	
QWT6	0.0528	0.0831	0.064	0.0657	0.0356	0.0649	0.067	
QWT7	0.0187	0.0387	0.022	0.044	0.0331	0.0096	0.045	
QWT8	0.0271	0.0166	0.0106	0.0521	0.0889	0.0588	0.0432	
QWT9	0.0424	0.0814	0.0267	0.0214	0.0034	0.0445	0.0814	
QWT10	0.0207	0.0882	0.0287	0.0413	0.0797	0.0701	0.0549	
QWT11	0.076	0.0395	0.0382	0.0867	0.0822	0.0644	0.0556	
QWT12	0.0175	0.01	0.0457	0.0492	0.0717	0.0813	0.0773	
Tab	−2	−2	−2	−2	−2	−2	−2	
Comp	ON	ON	ON	ON	ON	OFF	OFF	

sensors-21-05246-t008_Table 8 Table 8 Data of strategy I.

PWF	Proportional Distribution (PD)	
Q*PCC	Plosses
(MVAr)	QPCC*-QPCCmeas
(%)	
100%	4	0.0222	14.287	
100%	3.5	0.1743	10.725	
80%	3.5	0.1239	8.109	
50%	2	0.1478	2.312	
50%	3	0.1623	33.029	
20%	1	0.0820	11.023	
10%	0.5	0.0046	28.635	
PWF	HMPSO-VCS	
Plosses
(MVAr)	QPCC*−QPCCmeas
(%)	Reduction
Plosses %	
100%	0.1102	4.2122	0.0723	
100%	0.1101	4.2355	0.7019	
80%	0.0712	4.2197	1.7081	
50%	0.0266	4.2283	2.0564	
50%	0.0285	4.1091	6.1498	
20%	0.0023	4.2401	10.216	
10%	0.0010	4.2222	6.522	

sensors-21-05246-t009_Table 9 Table 9 Data for an option I.

WT Units	Strategy	
2	3	4	5	6	
QWT1	0.3756	0.4071	0.4586	0.2844	0.3446	
QWT2	0.1275	0.1218	0.1429	0.2347	0.3741	
QWT3	0.253	0.4646	0.3786	0.006	0.2253	
QWT4	0.3495	0.175	0.3769	0.1686	0.0419	
QWT5	0.4455	0.0983	0.1902	0.0811	0.1145	
QWT6	0.4796	0.1255	0.2839	0.3971	0.4567	
QWT7	0.2736	0.308	0.0379	0.1556	0.0762	
QWT8	0.0693	0.2366	0.027	0.2643	0.4129	
QWT9	0.0746	0.1758	0.2654	0.0828	0.2692	
QWT10	0.1288	0.4154	0.3896	0.301	0.4981	
QWT11	0.4204	0.2926	0.467	0.1315	0.0391	
QWT12	0.1271	0.2749	0.065	0.327	0.2213	
Comp	–	1	–	–	–	
Tab	–	–	−2	−2	−2	
QST	–	–	–	1.218	1.219	
Plosses	0.1209	0.1212	0.1021	0.1118	0.1119	
QPCC*−QPCCmeas (%)	4.9602	4.9401	4.9316	4.0345	40.721	

sensors-21-05246-t010_Table 10 Table 10 The mathematical detailed of employed benchmark functions, D: dimension, [L,U]: lower and upper bands, Fun: function name, No: number, Min: minimum value.

Many Local Minima Group	
No	Fun	[L,U]	D	Formulation	Min	
1	Ackley	[−32, 32]	30	f1(x)=−20exp(−0.21D∑i=1Dxi2)−exp(1D∑i=1Dcos(2πxi))+20+exp(1)	0	
2	Bukin N. 6	x1∈[−15,−5]x2∈[−3,3]	2	f2(x)=100|x2−0.01x12|+0.01|x1+10|	0	
3	Cross-in-tray function	[−10, 10]	2	f3(x)=−0.0001(|sin(x1)sin(x2)exp(|100−x12+x22π|)|+1)0.1	−2.06261	
4	Drop-wave	[−5.12, 5.12]	2	f4(x)=1+cos(12x12+x22)0.5(x12+x22)+2	−1	
5	Eggholder	[−512, 512]	2	f5(x)=−(x2+47)sin(|x2+x12+47|)−x1sin(|x1−x2−47|)	−959.6407	
6	Gramacy and Lee (2012)	[0.5, 2.5]	1	f6(x)=sin(10πx)2x+(x−1)4	−0.8690	
7	Griewank	[−600, 600]	30	f7(x)=∑i=1Dxi24000−∏i=1Dcos(xii)+1	0	

sensors-21-05246-t011_Table 11 Table 11 Statistical results obtained by GA, PSO, GSA, ABC, HBMO, GWO, COA, LCA, SSA, and VCS through 30 independent runs on mentioned benchmark functions in Table 1, Best, Worst, Mean, and STD denotes the best solution, the worst solution, the mean solution, and the standard deviation, respectively.

Fun. No		Algorithms	GA	PSO	GSA	ABC	COA	Proposed	
Indices		
f1	Best	2.142	7.053 × 10−3	6.983 × 10−6	5.982 × 10−10	7.873 × 10−6	2.948 × 10−12	
Worst	3.525	1.279 × 10−2	3.982 × 10−2	3.043 × 10−6	5.563 × 10−2	3.128 × 10−5	
Mean	3.240	9.858 × 10−1	0.195 × 10−1	0.989 × 10−9	2.098 × 10−4	1.029 × 10−9	
STD	3.058 × 10−1	1.586 × 10−2	1.329 × 10−1	1.093 × 10−2	9.340 × 10−2	1.928 × 10−2	
f2	Best	7.837 × 10−1	5.38	2.098	4.652	5.127 × 10−2	3.822 × 10−3	
Worst	9.651 × 10+2	9.043 × 10−1	8.675	8.341	3.428	5.938 × 10−1	
Mean	8.564 × 10−1	7.542 × 10−1	3.342	6.783	1.090 × 10−1	1.625 × 10−2	
STD	1.243 × 10−1	1.054 × 10+1	1.764 × 10−1	3.762	7.893	2.736 × 10−2	
f3	Best	−0.323	−0.645	−1.972	−1.234	−1.867	−1.948	
Worst	1.248	1.121	−1.146	0.653	−1.023	−1.226	
Mean	1.032	0.342 × 10−1	−1.345	−8.837 × 10−1	−1.648	−1.873	
STD	8.637 × 10−1	9.748	1.097 × 10−1	2.019 × 10−1	3.265 × 10−2	2.039 × 10−2	
f4	Best	−1.425 × 10−1	−4.837 × 10−1	−1.000	−8.948 × 10−1	−9.763 × 10−1	−1.000	
Worst	3.243	−1.093 × 10−1	−9.849 × 10−1	−2.192 × 10−1	−1.039 × 10−1	−0.827	
Mean	1.907 × 10−1	−1.623 × 10−1	−9.992 × 10−1	−7.765 × 10−1	−6.543 × 10−1	−0.928	
STD	1.783 × 10−1	5.425 × 10−1	−6.938 × 10−5	1.492 × 10-2	2.327 × 10−1	1.413 × 10−4	
f5	Best	−7.653 × 10−2	−9.321 × 10+2	−9.585 × 10+2	−9.596 × 10-2	−9.564 × 10+2	−9.596 × 10+2	
Worst	−6.536 × 10−2	−9.025 × 10+2	−9.413 × 10+2	−9.513 × 10-2	−9.323 × 10+2	−9.151 × 10+2	
Mean	−7.192 × 10+1	−9.294 × 10+2	−9.526 × 10+2	−9.576 × 10-2	−9.488 × 10+2	−9.029 × 10+2	
STD	3.201	3.052 × 10−1	1.0451 × 10-2	3.656 × 10-2	2.114 × 10−4	1.983 × 10−4	
f6	Best	−8.690 × 10−1	−8.690 × 10−1	−8.690 × 10−1	−8.690 × 10−1	−8.690 × 10−1	−8.690 × 10−1	
Worst	−8.689 × 10−1	−8.690 × 10−1	−8.690 × 10−1	−8.690 × 10−1	−8.690 × 10−1	−8.690 × 10−1	
Mean	−8.690 × 10−1	−8.6900 × 10−1	−8.690 × 10−1	−8.690 × 10−1	−8.690 × 10−1	−8.690 × 10−1	
STD	0	0	0	0	0	0	
f7	Best	1.208 × 10−1	0	0	4.837 × 10−8	2.897 × 10−8	0	
Worst	2.685 × 10−1	1.029	2.546 × 10−2	2.039 × 10−1	8.675 × 10−1	0	
Mean	1.952 × 10−1	1.938 × 10−2	1.645 × 10−5	5.038 × 10−7	2.653 × 10−4	0	
STD	3.915 × 10−2	1.523 × 10−3	5.320 × 10−4	9.456 × 10−14	3.332 × 10−5	0	

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Mehdinejad M. Mohammadi-Ivatloo B. Dadashzadeh-Bonab R. Zare K. Solution of optimal reactive power dispatch of power systems using hybrid particle swaarm optimization and imperialist competitive algorithms Electr. Power Energy Syst. 2016 83 104 116 10.1016/j.ijepes.2016.03.039
2. Guo W. Liu T. Dai F. Xu P. An improved whale optimization algorithm for feature selection Comput. Mater. Contin. 2020 62 337 354 10.32604/cmc.2020.06411
3. Zhao Y. Lu J. Yan Q. Lai L. Xu L. Research on cell manufacturing facility layout problem based on improved nsga-ii Comput. Mater. Contin. 2020 62 355 364 10.32604/cmc.2020.06396
4. Zhu J. Zhou B. Optimization design of rc ribbed floor system using eagle strategy with particle swarm optimization Comput. Mater. Contin. 2020 62 365 383 10.32604/cmc.2020.06655
5. Liao Z. Wang J. Zhang S. Cao J. Min G. Minimizing movement for target coverage and network connectivity in mobile sensor networks IEEE Trans. Parallel Distrib. Syst. 2014 26 1971 1983 10.1109/TPDS.2014.2333011
6. Gu K. Wu N. Yin B. Jia W. Secure data query framework for cloud and fog computing IEEE Trans. Netw. Serv. Manag. 2019 17 332 345 10.1109/TNSM.2019.2941869
7. Basu M. Quasi-oppositional differential evolution for optimal reactive power dispatch Electr. Power Energy Syst. 2016 78 29 40 10.1016/j.ijepes.2015.11.067
8. Siddiqui I.F. Lee S.U. Abbas A. A novel knowledge-based battery drain reducer for smart meters Intell. Autom. Soft Comput. 2020 26 107 119 10.31209/2019.100000132
9. Mezhuyev V. Gunchenko Y. Shvorov S. Chyrchenko D. A method for planning the routes of harvesting equipment using unmanned aerial vehicles Intell. Autom. Soft Comput. 2020 26 121 132
10. Chao L. Zhang K. Li Z. Zhu Y. Wang J. Yu Z. Geographically weighted regression based methods for merging satellite and gauge precipitation J. Hydrol. 2018 558 275 289 10.1016/j.jhydrol.2018.01.042
11. Zhang K. Wang Q. Chao L. Ye J. Li Z. Yu Z. Yang T. Ju Q. Ground observation-based analysis of soil moisture spatiotemporal variability across a humid to semi-humid transitional zone in China J. Hydrol. 2019 574 903 914 10.1016/j.jhydrol.2019.04.087
12. Zhao X. Gu B. Gao F. Chen S. Matching Model of Energy Supply and Demand of the Integrated Energy System in Coastal Areas J. Coast. Res. 2020 103 983 989 10.2112/SI103-205.1
13. Liu Y. Zhang B. Feng Y. Lv X. Ji D. Niu Z. Yang Y. Zhao X. Fan Y. Development of 340-GHz Transceiver Front End Based on GaAs Monolithic Integration Technology for THz Active Imaging Array Appl. Sci. 2020 10 7924 10.3390/app10217924
14. Niu Z. Zhang B. Wang J. Liu K. Chen Z. Yang K. Zhou Z. Fan Y. Zhang Y. Ji D. The research on 220GHz multicarrier high-speed communication system China Commun. 2020 17 131 139 10.23919/JCC.2020.03.011
15. Zhang Z. Liu M. Zhou M. Chen J. Dynamic reliability analysis of nonlinear structures using a Duffing-system-based equivalent nonlinear system method Int. J. Approx. Reason. 2020 126 84 97 10.1016/j.ijar.2020.08.006
16. Hu J. Zhang H. Li Z. Zhao C. Xu Z. Pan Q. Object traversing by monocular UAV in outdoor environment Asian J. Control. 2020 10.1002/asjc.2415
17. Hu J. Wang M. Zhao C. Pan Q. Du C. Formation control and collision avoidance for multi-UAV systems based on Voronoi partition Sci. China Ser. E Technol. Sci. 2019 63 65 72 10.1007/s11431-018-9449-9
18. Zuo X. Dong M. Gao F. Tian S. The Modeling of the Electric Heating and Cooling System of the Integrated Energy System in the Coastal Area J. Coast. Res. 2020 103 1022 1029 10.2112/SI103-213.1
19. Wang B. Ma F. Ge L. Ma H. Wang H. Mohamed M.A. Icing-EdgeNet: A Pruning Lightweight Edge Intelligent Method of Discriminative Driving Channel for Ice Thickness of Transmission Lines IEEE Trans. Instrum. Meas. 2020 70 1 12 10.1109/tim.2020.3018831 33776080
20. Research on evaluating vulnerability of integrated electricity-heat-gas systems based on high-dimensional random matrix theory CSEE J. Power Energy Syst. 2019 6 878 889 10.17775/cseejpes.2019.00440
21. Wang B. Zhang L. Ma H. Wang H. Wan S. Parallel LSTM-Based Regional Integrated Energy System Multienergy Source-Load Information Interactive Energy Prediction Complex. 2019 2019 7414318 10.1155/2019/7414318
22. Yin F. Xue X. Zhang C. Zhang K. Han J. Liu B. Wang J. Yao J. Multifidelity Genetic Transfer: An Efficient Framework for Production Optimization SPE J. 2021 1 22 10.2118/205013-PA
23. Xue X. Zhang K. Tan K.C. Feng L. Wang J. Chen G. Zhao X. Zhang L. Yao J. Affine Transformation-Enhanced Multifactorial Optimization for Heterogeneous Problems IEEE Trans. Cybern. 2020 1 15 10.1109/TCYB.2020.3005453 32701457
24. Deng R. Li M. Linghu S. Research on Calculation Method of Steam Absorption in Steam Injection Thermal Recovery Technology Fresenius Environ. Bull. 2021 30 5362 5369
25. Zhang L. Zheng H. Wan T. Shi D. Lyu L. Cai G. An Integrated Control Algorithm of Power Distribution for Islanded Microgrid Based on Improved Virtual Synchronous Generator IET Renew. Power Gener. 2021 10.1049/rpg2.12191
26. Zhang X. Wang Y. Wang C. Su C.-Y. Li Z. Chen X. Adaptive Estimated Inverse Output-Feedback Quantized Control for Piezoelectric Positioning Stage IEEE Trans. Cybern. 2018 49 2106 2118 10.1109/TCYB.2018.2826519 29994043
27. Cai X. Wang J. Zhong S. Shi K. Tang Y. Fuzzy quantized sampled-data control for extended dissipative analysis of T–S fuzzy system and its application to WPGSs J. Frankl. Inst. 2021 358 1350 1375 10.1016/j.jfranklin.2020.12.002
28. Cai X. Zhong S. Wang J. Shi K. Robust H∞ control for uncertain delayed T-S fuzzy systems with stochastic packet dropouts Appl. Math. Comput. 2020 385 125432 10.1016/j.amc.2020.125432
29. Cai X. Shi K. Zhong S. Wang J. Tang Y. Dissipative analysis for high speed train systems via looped-functional and relaxed condition methods Appl. Math. Model. 2021 96 570 583 10.1016/j.apm.2021.03.042
30. Cai X. Shi K. Zhong S. Pang X. Dissipative Sampled-Data Control for High-Speed Train Systems With Quantized Measurements IEEE Trans. Intell. Transp. Syst. 2021 1 12 10.1109/tits.2021.3052940
31. Hua L. Zhu H. Shi K. Zhong S. Tang Y. Liu Y. Novel Finite-Time Reliable Control Design for Memristor-Based Inertial Neural Networks With Mixed Time-Varying Delays IEEE Trans. Circuits Syst. I Regul. Pap. 2021 68 1599 1609 10.1109/TCSI.2021.3052210
32. Dong S. Zhu H. Zhong S. Shi K. Liu Y. New study on fixed-time synchronization control of delayed inertial memristive neural networks Appl. Math. Comput. 2021 399 126035 10.1016/j.amc.2021.126035
33. Liu C. Deng F. Heng Q. Cai X. Zhu R. Liserre M. Crossing Thyristor Branches-Based Hybrid Modular Multilevel Converters for DC Line Faults IEEE Trans. Ind. Electron. 2020 68 9719 9730 10.1109/TIE.2020.3026277
34. Yu B. Urban spatial structure and total-factor energy efficiency in Chinese provinces Ecol. Indic. 2021 126 107662 10.1016/j.ecolind.2021.107662
35. Xiao G. Song K. He Y. Wang W. Zhang Y. Dai W. Prediction and experimental research of abrasive belt grinding residual stress for titanium alloy based on analytical method Int. J. Adv. Manuf. Technol. 2021 1 15 10.1007/s00170-021-07272-3
36. Ni T. Liu D. Xu Q. Huang Z. Liang H. Yan A. Architecture of Cobweb-Based Redundant TSV for Clustered Faults IEEE Trans. Very Large Scale Integr. VLSI Syst. 2020 28 1736 1739 10.1109/TVLSI.2020.2995094
37. Li Y. Jia D. Rui Z. Peng J. Fu C. Zhang J. Evaluation method of rock brittleness based on statistical constitutive relations for rock damage J. Pet. Sci. Eng. 2017 153 123 132 10.1016/j.petrol.2017.03.041
38. Nejad R.M. Berto F. Wheatley G. Tohidi M. Ma W. On fatigue life prediction of Al-alloy 2024 plates in riveted joints Structures 2021 33 1715 1720 10.1016/j.istruc.2021.05.055
39. Huang J. Zhang J. Ren J. Chen H. Anti-rutting performance of the damping asphalt mixtures (DAMs) made with a high content of asphalt rubber (AR) Constr. Build. Mater. 2021 271 121878 10.1016/j.conbuildmat.2020.121878
40. Huang J. Kumar G.S. Sun Y. Evaluation of workability and mechanical properties of asphalt binder and mixture modified with waste toner Constr. Build. Mater. 2021 276 122230 10.1016/j.conbuildmat.2020.122230
41. Huang J. Wang Q.-A. Influence of crumb rubber particle sizes on rutting, low temperature cracking, fracture, and bond strength properties of asphalt binder Mater. Struct. 2021 54 1 15 10.1617/s11527-021-01647-4
42. Huang J. Sun Y. Zhang J. Reduction of computational error by optimizing SVR kernel coefficients to simulate concrete compressive strength through the use of a human learning optimization algorithm Eng. Comput. 2021 1 18 10.1007/s00366-021-01305-x
43. Huang J. Zhang Y. Sun Y. Ren J. Zhao Z. Zhang J. Evaluation of pore size distribution and permeability reduction behavior in pervious concrete Constr. Build. Mater. 2021 290 123228 10.1016/j.conbuildmat.2021.123228
44. Huang J. Kumar G.S. Ren J. Sun Y. Li Y. Wang C. Towards the potential usage of eggshell powder as bio-modifier for asphalt binder and mixture: Workability and mechanical properties Int. J. Pavement Eng. 2021 1 13 10.1080/10298436.2021.1905809
45. Lee B.K. Lee Y.S. Distinction between real faces and photos by analysis of face data Intell. Autom. Soft Comput. 2020 26 133 139 10.31209/2019.100000134
46. Tang Q. Wang K. Song Y. Li F. Park J.H. Waiting time minimized charging and discharging strategy based on mobile edge computing supported by software-defined network IEEE Internet Things J. 2019 7 6088 6101 10.1109/JIOT.2019.2957124
47. Tang Q. Wang K. Yang K. Luo Y.S. Congestion-balanced and welfare-maximized charging strategies for electric vehicles IEEE Trans. Parallel Distrib. Syst. 2020 31 2882 2895 10.1109/TPDS.2020.3003270
48. Huang C.-M. Huang Y.-C. Combined differential evolution algorithm and ant system for optimal reactive power dispatch Energy Procedia 2012 14 1238 1243 10.1016/j.egypro.2011.12.1082
49. He S. Xie K. Xie K. Xu C. Wang J. Interference-aware multisource transmission in multiradio and multichannel wireless network IEEE Syst. J. 2019 13 2507 2518 10.1109/JSYST.2019.2910409
50. Li W. Xu H. Li H. Yang Y. Sharma P.K. Wang J. Singh S. Complexity and algorithms for superposed data uploading problem in networks with smart devices IEEE Internet Things J. 2019 7 5882 5891 10.1109/JIOT.2019.2949352
51. Ping H. Study of assessment method based on coupling factor of casualty in earthquake disasters in guangdong area Comput. Syst. Sci. Eng. 2020 35 191 199 10.32604/csse.2020.35.191
52. Chen J. Luo Y. Du R. The impact of privacy seal on users’ perception in network transactions Comput. Syst. Sci. Eng. 2020 35 199 206 10.32604/csse.2020.35.199
53. Zhang Y. The implementation of an English word learning system feedback system and smartphone app Comput. Syst. Sci. Eng. 2020 35 207 214 10.32604/csse.2020.35.207
54. Davarpanah A. Parametric study of polymer-nanoparticles-assisted injectivity performance for axisymmetric two-phase flow in EOR processes Nanomaterials 2020 10 1818 10.3390/nano10091818 32932654
55. Zhou S. Tan B. Electrocardiogram soft computing using hybrid deep learning CNN-ELM Appl. Soft Comput. 2020 86 105778 10.1016/j.asoc.2019.105778
56. Long M. Chen Y. Peng F. Simple and accurate analysis of BER performance for DCSK chaotic communication IEEE Commun. Lett. 2011 15 1175 1177 10.1109/LCOMM.2011.092011.111088
57. Bahramian F. Akbari A. Nabavi M. Esfandi S. Naeiji E. Issakhov A. Design and tri-objective optimization of an energy plant integrated with near-zero energy building including energy storage: An application of dynamic simulation Sustain. Energy Technol. Assess. 2021 47 101419 10.1016/j.seta.2021.101419
58. Nabavi M. Elveny M. Danshina S.D. Behroyan I. Babanezhad M. Velocity prediction of Cu/water nanofluid convective flow in a circular tube: Learning CFD data by differential evolution algorithm based fuzzy inference system (DEFIS) Int. Commun. Heat Mass Transfer. 2021 126 105373 10.1016/j.icheatmasstransfer.2021.105373
59. Rezaei M. Farahanipad F. Dillhoff A. Elmasri R. Athitsos V. Weakly-supervised hand part segmentation from depth images Proceedings of the 14th PErvasive Technologies Related to Assistive Environments Conference New York, NY, USA 29 June 2021 218 225
60. Farahanipad F. Rezaei M. Dillhoff A. Kamangar F. Athitsos V. A Pipeline for Hand 2-D Keypoint Localization Using Unpaired Image to Image Translation Proceedings of the 14th PErvasive Technologies Related to Assistive Environments Conference New York, NY, USA 29 June 2021 226 233
61. Ehyaei M.A. Ahmadi A. Rosen M.A. Davarpanah A. Thermodynamic Optimization of a Geothermal Power Plant with a Genetic Algorithm in Two Stages Processes 2020 8 1277 10.3390/pr8101277
62. Habibifar R. Karimi M.R. Ranjbar H. Ehsan M. Economically based distributed battery energy storage systems planning in microgrids Proceedings of the Electrical Engineering (ICEE), Iranian Conference Mashhad, Iran 8–10 May 2018 1257 1263
63. Habibifar R. Khoshjahan M. Ghasemi M.A. Optimal Scheduling of Multi-Carrier Energy System Based on Energy Hub Concept Considering Power-to-Gas Storage Proceedings of the 2020 IEEE Power & Energy Society Innovative Smart Grid Technologies Conference (ISGT) Washington, DC, USA 17–20 February 2020 1 5
64. Sujod M. Erlich I. Engelhardt S. Improving the reactive power capability of the DFIG-based wind turbine during operation around the synchronous speed IEEE Trans. Energy Convers. 2013 28 736 745 10.1109/TEC.2013.2272975
65. Jafarishiadeh F. Mohammadi F. Sahraei-Ardakani M. Preventive Dispatch for Transmission De-icing IEEE Trans. Power Syst. 2020 35 4104 4107 10.1109/TPWRS.2020.3004079
66. Jin Y. Davarpanah A. Using photo-fenton and floatation techniques for the sustainable management of flow-back produced water reuse in shale reservoirs exploration Water Air Soil Pollut. 2020 231 1 8 10.1007/s11270-020-04812-7
67. Hu X. Xie J. Cai W. Wang R. Davarpanah A. Thermodynamic effects of cycling carbon dioxide injectivity in shale reservoirs J. Pet. Sci. Eng. 2020 195 107717 10.1016/j.petrol.2020.107717
68. Davarpanah A. Shirmohammadi R. Mirshekari B. Aslani A. Analysis of hydraulic fracturing techniques: Hybrid fuzzy approaches Arab. J. Geosci. 2019 12 402 10.1007/s12517-019-4567-x
69. Razmjoo A.A. Davarpanah A. Zargarian A. The Role of Renewable Energy to Achieve Energy Sustainability in Iran. An Economic and Technical Analysis of the Hybrid Power System Technol. Econ. Smart Grids Sustain. Energy 2019 4 7 10.1007/s40866-019-0063-3
70. Dibazar S.Y. Salehi G. Davarpanah A. Comparison of Exergy and Advanced Exergy Analysis in Three Different Organic Rankine Cycles Process. 2020 8 586 10.3390/pr8050586
71. Razmjoo A. Aliehyaei M. Ahmadi A. Pazhoohesh M. Marzband M. Khosravi M.M. Shahhoseini A. Davarpanah A. Implementation of energy sustainability using hybrid power systems, a case study Energy Sources Part A Recover. Util. Environ. Eff. 2019 1 14 10.1080/15567036.2019.1687623
72. Razmjoo A.A. Sumper A. Davarpanah A. Energy sustainability analysis based on SDGs for developing countries Energy Sources, Part A Recover. Util. Environ. Eff. 2019 42 1041 1056 10.1080/15567036.2019.1602215
73. Azma A. Narreie E. Shojaaddini A. Kianfar N. Kiyanfar R. Seyed Alizadeh S.M. Davarpanah A. Statistical modeling for spatial groundwater potential map based on gis technique Sustainability 2021 13 3788 10.3390/su13073788
74. Nagarajan K. Parvathy A.K. Rajagopalan A. Multi-Objective Optimal Reactive Power Dispatch using Levy Interior Search Algorithm Int. J. Electr. Eng. Inform. 2020 12 547 570 10.15676/ijeei.2020.12.3.8
75. Ghennam T. Aliouane K. Akel F. Francois B. Berkouk E.M. Advanced control system of DFIG based wind generators for reactive power production and integration in a wind farm dispatching Energy Convers. Manag. 2015 105 240 250 10.1016/j.enconman.2015.07.058
76. Wang S. Chen N. Yu D. Foley A. Zhu L. Li K. Yu J. Flexible fault ride through strategy for wind farm clusters in power systems with high wind power penetration Energy Convers. Manag. 2015 93 239 248 10.1016/j.enconman.2015.01.022
77. Sarker J. Goswami S.K. Solution of multiple UPFC placement problems using Gravitational Search Algorithm Electr. Power Energy Syst. 2014 55 531 541 10.1016/j.ijepes.2013.10.008
78. Gitizadeh M. Shidpilehvar M. Mardaneh M.d. A new method for SVC placement considering FSS limit and SVC investment cost Electr. Power Energy Syst. 2013 53 900 908 10.1016/j.ijepes.2013.06.009
79. Lu F.C. Hsu Y.Y. Reactive power/voltage control in a distribution substation using dynamic programming IEE Proc. Gener. Transm. Distrib. 1995 142 639 645 10.1049/ip-gtd:19952210
80. Quintana V.H. Santos-Nieto M. Reactive-power dispatch by successive quadratic programming IEEE Trans. Energy Conv. 1989 4 425 435 10.1109/60.43245
81. Granville S. Optimal reactive power dispatch through interior point methods IEEE Trans. Power Syst. 1994 9 136 146 10.1109/59.317548
82. Shaheen A.M. El-Sehiemy R.A. Farrag S.M. A novel adequate bi-level reactive power planning strategy Electr. Power Energy Syst. 2016 78 897 909 10.1016/j.ijepes.2015.12.004
83. Bin Z. Ka W.C. Tao Y. Hua W. Jie T. Strength Pareto Multigroup Search Optimizer for Multiobjective Optimal Reactive Power Dispatch IEEE Trans. Ind. Inform. 2014 10 1012 1022
84. Bhattacharyy B. Babu R. Teaching Learning Based Optimization algorithm for reactive power planning Electr. Power Energy Syst. 2016 81 248 253 10.1016/j.ijepes.2016.02.042
85. Shaheen A.M. Spea S.R. Farrag S.M. Abido M.A. A review of meta-heuristic algorithms for reactive power planning problem Ain Shams Eng. J. 2015 9 215 231 10.1016/j.asej.2015.12.003
86. Bhattacharyya B. Kumar S. Reactive power planning with FACTS devices using gravitational search algorithm Ain Shams Eng. J. 2015 6 865 871 10.1016/j.asej.2015.02.002
87. Kanna B. Singh S.N. Towards reactive power dispatch within a wind farm using hybrid PSO Electr. Power Energy Syst. 2015 69 232 240 10.1016/j.ijepes.2015.01.021
88. Xian L. Wilsun X. Minimum Emission Dispatch Constrained by Stochastic Wind Power Availability and Cost IEEE Trans. Power Syst. 2010 25 1705 1713 10.1109/TPWRS.2010.2042085
89. El-sobky B. Abo-elnaga Y. Multi-objective economic emission load dispatch problem with trust-region strategy Electr. Power Syst. Res. 2014 108 254 259 10.1016/j.epsr.2013.11.003
90. Arul R. Ravi G. Velusami S. Chaotic self-adaptive differential harmony search algorithm based dynamic economic dispatch Electr. Power Energy Syst. 2013 50 85 96 10.1016/j.ijepes.2013.02.017
91. Li M.D. Zhao H. Weng X.W. Han T. A novel nature-inspired algorithm for optimization: Virus colony search Adv. Eng. Softw. 2016 92 65 88 10.1016/j.advengsoft.2015.11.004
92. Ghasemi A. Valipour K. Tohidi A. Multi-objective optimal reactive power dispatch using a new multi-objective strategy Electr. Power Energy Syst. 2014 57 318 334 10.1016/j.ijepes.2013.11.049
93. Höppner F. Klawonn F. Obtaining interpretable fuzzy models from fuzzy clustering and fuzzy regression Proceedings of the Fourth International Conference on Knowledge-Based Intelligent Engineering Systems and Allied Technologies Salt Lake City, UT, USA 9–10 October 2000 162 165
94. Chowdhury M.A. Shen W.X. Hosseinzadeh N. Pota H.R. A novel aggregated DFIG wind farm model using mechanical torque compensating factor Energy Convers Manag. 2013 67 265 274 10.1016/j.enconman.2012.12.001
95. Hansen N. Müller S.D. Koumoutsakos P. Reducing the time complexity of the derandomized evolution strategy with covariance matrix adaptation (CMAES) Evol. Comput. 2003 11 1 18 10.1162/106365603321828970 12804094
96. Kennedy J. Eberhart R. Particle swarm optimization IEEE Int. Conf. Neural Netw. 1995 4 1942 1948
97. Zou D. Li S. Li Z. Kong X. A new global particle swarm optimization for the economic emission dispatch with or without transmission losses Energy Convers. Manag. 2017 139 45 70 10.1016/j.enconman.2017.02.035
98. Khorram E. Khaledian K. Khaledyan M. A numerical method for constructing the Pareto front of multi-objective optimization problems J. Comput. Appl. Math. 2014 261 158 171 10.1016/j.cam.2013.11.007
99. Srivastava L. Singh H. Hybrid multi-swarm particle swarm optimization based multi-objective reactive power dispatch IET Gener. Transm. Distrib. 2015 9 727 739 10.1049/iet-gtd.2014.0469
100. Bhattacharyya B. Raj S. Swarm intelligence based algorithms for reactive power planning with Flexible AC transmission system devices Electr. Power Energy Syst. 2016 78 158 164 10.1016/j.ijepes.2015.11.086
101. Wiik J. Gjerde J.O. Gjengedal T. Steady-state power system issues when planning large wind farms IEEE Power Eng. Soc. Winter Meet. 2002 1 657 661
102. Zeng X.J. Tao J. Zhang P. Pan H. Wang Y.Y. Reactive Power Optimization of Wind Farm based on Improved Genetic Algorithm Energy Procedia 2012 14 1362 1367 10.1016/j.egypro.2011.12.1102
103. Martinez-Rojas M. Sumper A. Gomis-Bellmunt O. Sudrià-Andreu A. Reactive power dispatch in wind farms using particle swarm optimization technique and feasible solutions search Appl. Energy 2011 88 4678 4686 10.1016/j.apenergy.2011.06.010
104. Cao Y. Doustgani A. Salehi A. Nemati M. Ghasemi A. Koohshekan O. The economic evaluation of establishing a plant for producing biodiesel from edible oil wastes in oil-rich countries: Case study Iran Energy 2020 213 118760 10.1016/j.energy.2020.118760
105. Karim S.H.T. Tofiq T.A. Shariati M. Rad H.N. Ghasemi A. 4E analyses and multi-objective optimization of a solar-based combined cooling, heating, and power system for residential applications Energy Rep. 2021 7 1780 1797 10.1016/j.egyr.2021.03.020
106. Ghasemi A. Moghaddam M. Thermodynamic and environmental comparative investigation and optimization of landfill vs. incineration for municipal solid waste: A case study in varamin, Iran J. Therm. Eng. 2020 6 226 246 10.18186/thermal.820234
107. Ghasemi A. Shayesteh A.A. Doustgani A. Pazoki M. Thermodynamic assessment and optimization of a novel trigeneration energy system based on solar energy and MSW gasification using energy and exergy concept J. Therm. Eng. 2020 7 349 366 10.18186/thermal.850819
