 Original research
 Open Access
 Published:
Single and multiarea multifuel economic dispatch using a fuzzified squirrel search algorithm
Protection and Control of Modern Power Systems volume 6, Article number: 11 (2021)
Abstract
MultiArea MultiFuel Economic Dispatch (MAMFED) aims to allocate the best generation schedule in each area and to offer the best power transfers between different areas by minimizing the objective functions among the available fuel alternatives for each unit while satisfying various constraints in power systems. In this paper, a Fuzzified Squirrel Search Algorithm (FSSA) algorithm is proposed to solve the singlearea multifuel economic dispatch (SAMFED) and MAMFED problems. Squirrel Search Algorithm (SSA) mimics the foraging behavior of squirrels based on the dynamic jumping and gliding strategies. In the SSA approach, predator presence behavior and a seasonal monitoring condition are employed to increase the search ability of the algorithm, and to balance the exploitation and exploration. The suggested approach considers the line losses, valve point loading impacts, multifuel alternatives, and tieline limits of the power system. Because of the contradicting nature of fuel cost and pollutant emission objectives, weighted sum approach and price penalty factor are used to transfer the biobjective function into a single objective function. Furthermore, a fuzzy decision strategy is introduced to find one of the Pareto optimal fronts as the best compromised solution. The feasibility of the FSSA is tested on a threearea test system for both the SAMFED and MAMFED problems. The results of FSSA approach are compared with other heuristic approaches in the literature. Multiobjective performance indicators such as generational distance, spacing metric and ratio of nondominated individuals are evaluated to validate the effectiveness of FSSA. The results divulge that the FSSA is a promising approach to solve the SAMFED and MAMFED problems while providing a better compromise solution in comparison with other heuristic approaches.
Introduction
The goal of MultiArea MultiFuel Economic Dispatch (MAMFED) is to decide the power delivered by every generator in all the areas and the power flow between the areas in order to lessen the total production cost and pollutants outflows of the interconnected power system considering multifuel alternatives of each generating unit. The total load power is shared among the various areas to decrease the fuel cost as well as the pollutant emission, while the power scheduling between the zones must consider the power balance, generator and tieline limits. Thus, the MAMFEED problem which has received considerable attention in recent years, is addressed as a largescale highly nonlinear multiobjective optimization problem.
Literature review
Over the years, various classical multiarea power generation scheduling methods have been proposed. An efficient approach is applied in [1] for unit commitment and Economic Dispatch (ED) problems with area import/export constraints. The proposed approach is tested on the Texas utilities and Texas municipal power pool systems. In [2], it suggests the use of the DantzigWolfe decomposition principle regarding the revised simplex method and a fastdecoupled power flow algorithm for constrained MultiArea ED (MAED) of power systems. A model of largescale multiarea power generation system is developed in [3] and a rulebased heuristic strategy is used to improve the generation schedule for every zone. Reference [4] proposes a decomposition approach which upgrades the scheduling process and accelerates the execution of a largescale multiarea generation system in a realtime application. The MAED problem is expressed as a capacitated nonlinear network flow problem in [5] and is solved through an Incremental Network Flow Programming (INFP) approach. The results show that the INFP approach is fast, robust, and extendable to largescale systems. An Improved Hopfield Neural Network (IHNN) is presented to solve MAED problems with transmission capacity constraints which shows that IHNN achieves efficient and accurate solutions for the twoarea power system [6]. In [7], a Hopfield Lagrange network is proposed to solve economic emission dispatch (EED) problem with multiple fuel alternatives (MFAs), while the best compromise from the set of obtained solutions is found and compared with that from the LamdaIteration method.
Despite the fact that these traditional strategies can give a strong ideal solution, they experience the disadvantages of having nonsmooth and nonconvex objective functions. Thus, these strategies do not offer reasonable possibilities for dealing with the MAED problem when prohibited operating zones (POZs), valve point loading (VPL), and MFAs are considered [8]. To adapt to these challenges, numerous metaheuristic methodologies have been used to take care of enhancement issues with complicated objectives [9]. The results obtained by evolutionary approaches for complex optimizing problems are good, which means these approaches offer a reasonable possibility for solving MAED problems.
In recent years, swarm intelligence algorithms have been broadly used to overcome the computational unpredictability issues in the MAED problem. Reference [10] proposes an efficient technique for MAED problems using an Evolutionary Programming (EP) approach, while the performance of the various evolutionary algorithms, including the RealCoded Genetic Algorithm (RCGA), Particle Swarm Optimization (PSO), Differential Evolution (DE), and the Covariance Matrix Adapted Evolution Strategy (CMAES), on MAED problems with Karush–Kuhn–Tucker optimality conditions are examined [11]. The simulation results revealed that the CMAES algorithm offers the best results and the Nelder–Mead simplex method offers an optimal solution in the shortest time among the different algorithms considered. The performance of different DE strategies enhanced with timevarying mutation is investigated and analyzed to solve the reserveconstrained MAED problem [12]. The time varying DE approach is found to be capable of determining better solutions than PSO with timevarying acceleration coefficients. In addition, the DE variant approaches are able to find the global best solutions for largescale systems, while basic DE suffers from premature convergence.
A computationally efficient fuzzified PSO algorithm is introduced in [13] to solve the securityconstrained MAED problem of an interconnected power system. The inertia weight of the standard PSO is made adaptive by using fuzzy logic strategy to improve the convergence speed and avoid premature convergence. Artificial bee colony (ABC) optimization is proposed to solve a MAED problem with tieline constraints, transmission losses, multiple fuels, and valve point effects [14]. Evolutionary approaches such as DE, EP, and RCGA are applied to analyze the efficiency of the ABC approach. The simulation results show that the ABC approach converges to a better solution than other comparable approaches. Teaching learningbased optimization (TLBO) has been applied to solve the MAED issue [15]. The approach is tested on three different systems and the results reveal that the presented approach has the ability to provide better solutions and exhibits greater robustness than DE, EP, and RCGA. The different MAED models are solved by employing Fast Convergence Evolutionary Programming (FCEP) in [16] which uses Gaussian and Cauchy mutations to improve the convergence speed and solution quality. The Hybrid Cuckoo Search Algorithm (HCSA) which consolidates the cuckoo search algorithm and TLBO to improve the performance of optimization process, is developed to solve the MAED problem [17]. In [18], a hybrid DEPSO technique is used to address the MAED, reserveconstrained MAELD, and reserveconstrained multiarea ecological/economic dispatch problems. Sum Differential Evolution with Particle Swarm Optimizer (SDEPSO) is used along with the local (Pbest) optimal value in a DE crossover operator. The simulation results show that the hybrid algorithm achieves a good balance between global search ability and provision of better convergence and solution quality.
An Improved Grasshopper Optimization Algorithm (IGOA) is introduced to deal with the MAED problem [19], in which a chaos mechanism is adopted to handle the premature convergence of the grasshopper optimization algorithm. The simulation results demonstrate that IGOA converges better than the other approaches. Backtracking Search Algorithm (BSA) is proposed to solve ED problems considering VPL effects, POZs, and MFAs [20], while the MAED problem is solved by combining Jaya and TLBO approaches to simultaneously solve both nonsmooth and nonconvex characteristics [21]. The Jaya algorithm consistently attempts to find the best solution and avoids the worst solution, and thus is endeavoured to be successful by arriving at the best solution. A Parallel MAED (PMAED) is developed to maintain the independency and transfer the required power in each area [22], which is solved with convex and nonconvex fuel cost functions of different complexities. A Newton Method based Distributed Algorithm (NMDA) and Evolutionary PSO (EPSO) algorithm are employed to solve SingleArea Economic Dispatch (SAED) and MAED problems [23, 24], where the virtual SAED formulation, backtracking line search algorithm, and average consensus theory are embedded in the NMDA approach. In the EPSO algorithm, the GA operators such as mutation, crossover, and selection are employed to enable the search process. A new heuristic approach, the Coulomb’s and Franklin’s lawsbased optimization (CFLBO) algorithm is developed to solve the nonconvex economic and emission dispatch problem [25]. These improved and hybrid forms of heuristic approaches involve complicated computation owing to the use of many control parameters. Thus, an efficient algorithm to solve MAED problems to ascertain the optimum dispatch solutions needs be developed.
Contributions
Recently, a new metaheuristic algorithm, named as Squirrel Search Algorithm (SSA) was proposed in [26]. SSA models the foraging activities of individual squirrels, while each squirrel modifies its position using four processes, namely:

(1)
distributing the population,

(2)
dynamic foraging behavior,

(3)
seasonal adapting intelligence, and.

(4)
random repositioning of individuals at the end of winter season.
The unique features of SSA are as follows:

The gliding constant is used in the location update of squirrels which provides suitable steadiness between exploration and exploitation.

The predator presence behavior is employed to abruptly change the squirrel locations which enhances the exploration ability of the algorithm.

A seasonal monitoring condition is used to prevent the suggested algorithm from being trapped in local optimal solutions.

Levy distribution is used to find new solutions far away from the current best solution which improves the global exploration ability of the algorithm.
These features make SSA capable of overcoming the normal drawbacks of other algorithms such as premature convergence, inadequate ability to find nearby extreme points and absence of efficient constraints handling mechanism. The advantages of the SSA approach are less execution time, ability to solve different complex optimization problems and high capacity in obtaining global optimum solutions. Thus, the SSA approach is able to provide solutions of better quality than the existing heuristic approaches.
The main contributions of this paper are briefly summarized as follows:

For the first time, the application of the SSA approach for solving the SingleArea MultiFuel ED (SAMFED) and MAMFED problems in power system is reported.

A FSSA approach to solve the SAMFED and MAMFED problems is introduced. The proposed algorithm involves the SSA approach, weighted sum approach and fuzzy decision strategy. The weighted sum approach is used to transfer the biobjective functions such as fuel cost and pollutant emission into a single objective function, while the fuzzy decision strategy is used to provide the best compromised solution from the nondominated solution set. The FSSA approach successfully solves the SAMFED problem such that the fuel cost and pollutant emission are simultaneously minimized while fulfilling the power balance and generation limits.

The SSA approach is successfully applied to solve the MFMAED such that the total fuel cost is minimized while fulfilling the system and area constraints.

The proposed algorithm is applied to a threearea 10unit system to demonstrate its effectiveness in solving the SAMFED and MAMFED problems. The obtained results are compared with different stateoftheart heuristic approaches.
The structure of this paper is as follows. Section 2 presents the details of the SAMFED and MAMFED models, while the recommended FSSA approach is presented in Section 3. The implementation of the FSSA approach for the ED problems is discussed in Section 4 and in Section 5, the numerical results are discussed and the adequacy of the proposed FSSA approach is demonstrated. Finally, Section 6 concludes the paper.
Problem formulation of multifuel economic dispatch
SAMFED
Fuel cost objective function
The goal of the ED problem is to limit the total expense of thermal units as:
In the ED problem, the fuel cost of each generating unit is communicated as a quadratic function of its capacity yield. As the generating units utilize different fuel choices to create power and consider the valve point impacts, the fuel cost function is expressed in (6).
The multifuel cost curves with and without VPL impacts are shown in Fig. 1.
Emission objective function
The objective of the emission dispatch (EmD) problem is to lessen the hazardous emissions as:
The emission objective function is similar to the fuel cost function while it deals with all discharge types by generation units. The model for emanation function with MFA is introduced in (7).
Economic and emission dispatch (EED) function
The EED problem can be considered as bitarget work with fuel cost and discharge as equalling objectives. This bitarget function can be merged into a solitary target function as:
where w is a primary function of rand [0, 1] which bargains the fuel cost and emanation objectives. So the above condition becomes ED target function when w = 1 and becomes EmD target when w = 0.
Power balance constraint
The total power generated from a set of committed units must fulfil the total load demand, i.e.:
Generator capacity limits
The real output power of thermal units need to be in their minimum and maximum limits as:
MAMFED
Fuel cost objective function
The aim of the MAMFED problem is to find the amount of power that can be efficiently generated in one area and transferred to another area, and to determine the economic fuel choice for each unit. Since generators are provided with multifuel sources, every generator needs to be defined with a few piecewise quadratic capacities superimposed by sine terms mirroring the impact of changes in the type of fuel while the generator must determine the most conservative fuel to consume. The fuel cost function with VPL and MFA impacts [14] is defined in (8).
Power balance constraint
The total power generated by a set of accessible units must satisfy the total load demand, tieline power flow, and transmission losses [9] as:
where i = 1, 2,. . n_{g}, j = 1, 2,. .. M_{i}.
Generator capacity limits
The real output power of the thermal units should be in the minimum and maximum limits as [9]:
Tieline limit
Because of security basis, power transmitted through various lines must not surpass their cutoff points [11]. The power transfer requirement between two regions is characterized by:
Brief overview of the FSSA approach
SSA
The hunt procedure starts when flying squirrels begin scavenging [26]. During fall, squirrels look for nourishment assets by skimming from one tree to the next. At the same time, they change their areas and investigate various regions of woods. As the climatic conditions are sufficiently hot, they can meet their every day vitality needs more rapidly on the eating routine of oak seeds accessible in bounty and thus, they devour oak seeds quickly after discovering them. Subsequent to satisfying their day by day vitality prerequisite, squirrels scan for ideal nourishment hotspot for winter (hickory nuts). Capacity of hickory nuts will help them in keeping up their vitality prerequisites in harsh climate, decrease the expensive searching excursions and increase the likelihood of endurance.
During winter, lost leaves spread in deciduous woodlands result an expanded danger of predation and thus, squirrels become less dynamic but still remain active. Toward the finish of winter season, squirrels again become dynamic. This is monotonous procedure and structures the establishment of SSA. The SSA approach refreshes the places of squirrels as indicated by the ebb and flow season, the sort of squirrels and if chasers showing up.
Instate the population
Assuming the number of squirrels is N, and the upper and lower limits of the pursuit space are X_{U} and X_{L}, the N squirrels are arbitrarily created as:
where X_{i} indicates the i^{th} squirrel, (i = 1: N), rand () is an random number in the range of 0 and 1, and D is the measurement of the issue.
Group the population
SSA requires that there is only a single squirrel at each tree, so for N squirrels, there are N trees in the woods. Among the N trees, there is one hickory tree and N_{a} oak seed trees, while the rests are typical trees having no nourishment. The hickory tree is the best nourishment asset for the squirrels while the oak seed trees come the second. Positioning the fitness estimations of the populace in rising request, the squirrels are separated into three kinds:

Squirrels situated at hickory tree (W_{h});

Squirrels situated at oak seed trees (W_{a});

Squirrels situated at ordinary trees (W_{n}).
Refresh the location of squirrels
The squirrels refresh their situations by skimming to the hickory tree or oak seed trees as follows:
P_{dp} is esteemed at 0.1 and indicates the chaser likelihood. In the event that r > P_{dp}, no chaser shows up, and the squirrels coast in the backwoods to discover the nourishment and are protected. If r < P_{dp}, the chasers show up, and the squirrels are compelled to limit the extent of exercises and are imperiled, and their locations are migrated arbitrarily. d_{g} is the skimming separation that can be determined by:
where h_{g} is the constant estimated 8, tan (ɸ) indicates the coasting point that can be determined by:
The drag power and lift power can be estimated as:
Occasional changeover verdict and arbitrary refreshing
Toward the start of every generation, SSA necessitates that the entire populace is in winter, which implies that the locations of all squirrels are updated by (11) and (12). At the point when the squirrels are refreshed, regardless the season, change is decided by the following formulae:
If S_{tc} < S_{min}, winter is finished and the season goes to summer, otherwise the season is unaltered. At the point when the season goes to summer, the squirrels who float to W_{h} remain at the refreshed area, while the squirrels skimming to W_{a} and not meeting with chasers move their situations as follows:
Le’vy is the arbitrary walk model whose progression complies with the Le’vy appropriation and can be determined by:
where α is determined as:
Fuzzy decision strategy
The objective functions, viz. economic and emission dispatch shall be used in parallel in the multiobjective EELD problem. However, this makes the comparison of the two solutions difficult. Neither solution vector X_{1} nor X_{2} can be superior to each other if they are Pareto optimal, because if superior result is obtained from X_{1} for an objective, X_{2} would offer improved performance for another objective. Obtaining the best solution from multiple nondominated solutions is challenging in multiobjective EELD problem, though it is always possible to collate these outcomes and obtain the best compromised solution. For achieving this, one has to use the proper mechanism to combine both objectives and ensure that it conforms to the target and preference of the decision maker.
Researchers commonly use fuzzy set theory to arrive at the best solution amongst many uncontrolled solutions. It is implausible to achieve both least fuel cost along with least emission as they are contrary to each other. But it is feasible to build a dispatch option that can optimize both. Fuzzy membership functions assign Degree of agreement (DA) to each objective, and merit of the objective is reflected by DA in a linear scale of 0 – 1 (worst to best). F_{j} is a solution in the Paretooptimal set in the j^{th} objective function and is defined by a membership function as:
For each nondominated solution, the normalized membership function \( {\mu}_D^k \) can be calculated as:
The solution that contains the maximum of \( {\upmu}_{\mathrm{D}}^{\mathrm{k}} \) \( {\mu}_D^k \) based on cardinal priority ranking is the best compromised solution, i.e.:
Implementation of FSSA for SAMFED and MAMFED problems
The steps of the solution for MFED problem using FSSA are as follows.

Step 1: Randomly generate the generation values between lower and upper power outputs of each generating unit in all solutions.

Step 2: Evaluate the objective values of all the squirrels using (3).

Step 3: Sort the objective values of each squirrel’s position in ascending order.

Step 4: Declare the flying squirrel with low fitness value as on the hickory nut tree (optimal food source), the next three best flying squirrels on the acorn tree (normal food source) and the rest of the squirrels on the normal trees (no food source).

Step 5: Update the positions of squirrels located on the oak seed and normal trees using (12), (13) and (14).

Step 6: Randomly relocate the positions of some squirrels when seasonal monitoring condition is satisfied.

Step 7: Repeat steps 2 to 6 until stopping criterion is met.

Step 8: Increment the weight factor in step of 0.1 and repeat steps 2 to 6, until the weight factor reaches unity.
The membership value for each nondominated solution sets which are acquired for different weight factors using (24) is determined, and the position of squirrel on hickory nut tree that procures maximum membership value is chosen as the best compromising solution for the MFED problem.
The flowchart of the FSSA approach to solve EELD problem is shown in Fig. 2.
Numerical results and discussion
To demonstrate the effectiveness of the proposed FSSA approach, it is applied to a threearea, 10unit power system for the SAMFED and MAMFED problems. The coefficients of multifuel cost and emission, generation limits of the generating units, and tieline capacity limits are taken from [7].
The proposed FSSA is executed in Matlab 7.1 on an Intel core i3 processor with 4GB RAM PC, for 50 autonomous trials to assess the solution quality and convergence behavior. The parameter settings of the proposed FBSSA are presented in Table 1.
Three scenarios are researched and the results of the proposed FSSA are compared with Krill Herd Algorithm (KHA), Exchange Market Algorithm (EMA) and other stateofart heuristic approaches.

Scenario 1. Minimization of the total fuel cost and total pollutant emission individually in the SAMFED problem.

Scenario 2. Minimization of both the total fuel cost and the total emission level simultaneously in the SAMFED problem.

Scenario 3. Minimization of the total fuel cost in the MAMFED problem.
Scenario 1
Table 2 presents the optimal generation schedule acquired by the fuel cost minimization of the SAMFED using the proposed FSSA approach. The proposed approach is compared with the known approaches such as EMA, KHA, and BSA [20] in Fig. 3. As can be seen from Fig. 3, the FSSA approach outflanks the aforementioned approaches with the lowest total fuel cost.
The optimal values of generating units obtained by the FSSA for emission minimization are given in Table 3, with the minimum emission value of 6042.5066 kg/h. Figure 4 compares the results from different approaches, which shows that the FSSA approach reduces emission level by 9.31 kg/h and 12.07 kg/h compared to the KHA and EMA approaches, respectively.
Scenario 2
For this scenario, the total fuel cost and emission of the SAMFED are simultaneously minimized, so the objective function of the EED given in (3) is minimized. In order to provide the best tradeoff, the weights w_{1} and w_{2} are varied in the range of 0 to 1 in step of 0.1 so that their sum is one. The nondominated solution sets acquired by the FSSA approach are depicted in Table 4, and the solution corresponding to the maximum membership value estimation is chosen as the BCS. It is seen from Table 4 that the BCS are determined at w_{1} = 0.6 and w_{2} = 0.4 corresponding to μ_{D} = 0.100820. For the FSSA algorithm, the total fuel cost is 641.9474 and the total pollutant emission is 6195.3613 kg/h. Table 5 presents the optimal generation schedule of the SAMFED problem by minimizing both the total fuel cost and emission.
The performance indices of the EED problem such as fuel cost performance index (FCPI) and emission cost performance index (ECPI) are further performed. FCPI and ECPI can be expressed as:
The performance indices of the FSSA approach are compared with the results obtained using the KHA and EMA approaches, as shown in Fig. 5. It is evident that the performance indices acquired using the FSSA approach are superior to the other approaches, with lower divergence. Figure 6 shows the Pareto optimal fronts (POF) curves of different approaches. It is seen that the POF curve of the FSSA algorithm remains closer to the axes, which implies that the FSSA approach offers better results than the KHA and EMA approaches.
Scenario 3
In this scenario, fuel cost minimization of MAMFED is addressed. A 3area, 10unit test system is tested with transmission losses, and VPL and MFA impacts taken into consideration. Areas 1, 2, and 3 comprise four, three and three generating units, respectively as displayed in Fig. 7. The total power demand of the system is 2700 MW, and the power demand shares of areas 1, 2, and 3 are 50%, 25%, and 25% of the total load demand, respectively. The power flow from one area to another is restricted to 100 MW. Table 6 presents the simulation results obtained by the proposed SSA approach. It can be seen that the optimal generation cost obtained by the SSA approach is 654.4665 $/h, which is the lowest among all the approaches. Area 1 imports power from areas 2 and 3, while area 3 also exports power to area 2.
The results of the SSA approach are compared with those of RCGA [14], ABC [14], TLBO [15], EPSO [24], EMA and CFLBO approaches in Fig. 8, and it shows that the proposed strategy outperforms the others with regard to finding the best generation schedule.
Statistical analysis of the FSSA algorithm
Multiobjective performance indicators
The multiobjective performance indicators such as generational distance (GD), spacing metric (SM) and ratio of nondominated individuals (RNI) are used to examine the performance of the proposed approach for Scenario 3.
The smaller estimations of GD and SM demonstrate better union to the POF curve, and better conveyance and assorted variety of the nondominated solutions, respectively. A closer estimation of one for RNI index demonstrates that the greatest number of arrangements in a populace are nondominated. The GD, SM, and RNI measures obtained by the proposed FSSA approach in 50 autonomous trials are shown as box and whisker plots in Fig. 9, which indicate that the FSSA algorithm has estimated the smaller values regarding the GD and SM indices, and closest to one for the RNI indicator.
Fuel cost improvement percentage
Fuel cost improvement percentage (IP) is the ratio of fuel cost difference obtained between two approaches to the higher value of fuel cost obtained, expressed as a percentage as:
The IPs obtained by the FSSA approach and the existing heuristic approaches for scenario 3 are compared in Fig. 10, ranging from 0.113% to 0.5234%. It is noteworthy that the IP of the FSSA approach is high for MAMFED with VPL impacts. Thus, it can conclude that the proposed FSSA approach provides better results than the other approaches.
Computational efficiency
Tables 2, 5, and 6 show that the minimum fuel costs achieved by the FSSA approach are 623.7129, 641.9474, and 654.4665$/h for scenarios 1, 2, and 3, respectively. These costs are lower than the ones presented in recent literature. Figure 11 shows the number of function evaluation adopted by the FSSA and EMA strategies for the various scenarios, indicating lower number of function evaluation with FSSA. Thus, the FSSA technique is more computationally effective than the EMA strategy.
Convergence graph
The convergences of the EMA and FSSA approaches are compared in Fig. 12. As shown, the FSSA approach requires lower number of iterations to converge to the globally optimal solution.
Conclusion
The MAMFED problem, as an imperative issue in modern power system analysis, is addressed in this paper. The optimal generation schedule of all committed units and power transfer between the areas through tielines are determined. The FSSA approach is used to solve the SAMFED and the MAMFED problems, while the transmission line losses, VPL impacts, and tieline limits of the power system are addressed. The biobjective function is transferred into a single objective function by weighted sum approach. Moreover, a fuzzy decision strategy is introduced to find one of the Pareto optimal fronts as the best compromised solution. To demonstrate the effectiveness of the FSSA approach in terms of solution quality, computational efficiency and convergence speed, it is used to solve the MFED problems with three different scenarios and the results are compared with other existing approaches. The results prove that the proposed FSSA approach is an efficient method for solving the MFMAED problem and can provide better compromised solution than the other approaches. This research work not only offers an advanced optimizer for the MFMAED problem but also advances the use of the evolutionary approaches in the energy optimization domain. For future work, it will be intriguing to implement this compelling approach to solve other economic operation problems of hybrid windsolarthermal power systems.
Nomenclature

a_{i}, b_{i}, c_{i} cost coefficients of generator i

a_{ij}, b_{ij}, c_{ij} cost coefficients of generator i in area j

C_{D} drag coefficient

C_{L} lift coefficient

d_{i}, e_{i} cost coefficients of the VPL effect of generator i

d_{g} gliding distance

D drag force

e_{ij}, f_{ij} cost coefficients of the VPL effect of generator i in area j

E_{i} emission of the generator i

F_{i} fuel cost of the generator i

F_{ij}(P_{ij}) fuel cost of the generator i in area j

F_{j}^{max} \( {\mathrm{F}}_{\mathrm{j}}^{\mathrm{max}} \) and F_{j}^{min} \( {\mathrm{F}}_{\mathrm{j}}^{\mathrm{min}} \) maximum and minimum values of jth objective function respectively

F_{bcs} and E_{bcs} fuel cost and emission attained by EED

F_{min} and E_{max} fuel cost and emission attained by ED minimization respectively

F_{max} and E_{min} fuel cost and emission attained by EmD minimization respectively

G_{c} Gliding constant

h price penalty factor in $/h

h_{g} gliding height

IP improvement percentage

k number of fuel alternatives

L lift force

M_{i} number of participated generators in area i

M number of nondominated solutions

n maximum number of objects

ng total number of generating units

P_{D} power demand

P_{ij} Real power generation of generator j in area i

P_{i},_{min}, P_{i},_{max} minimum and maximum generation of unit i

P_{ij},_{min}, P_{ij},_{max} Minimum and maximum generation j in area i

P_{dp} predator presence probability

r_{a} and r_{b} randomly distributed numbers in [0, 1]

r_{1}, r_{2} and r_{3} random numbers in the range of [0, 1]

S surface area of body

S_{min} minimum value of seasonal constant

t current iteration

t_{max} maximum iteration value

T_{iz} tie line power stream from area i to area z

T_{iz, max} maximum tie line power stream from area i to area z

T_{iz, max} Maximum tie line power stream from area z to area i

V speed

w weight or compromise factor

X_{h} position of squirrel individual which reached the hickory tree

X_{L}, X_{U} lower and upper bounds of squirrel individual

x_{j}^{min} and x_{j}^{max} lower and upper limits of variable j

tan (ɸ) gliding angle

β constant

ρ density of air emission coefficients of generator i
Availability of data and materials
No data has been taken from any sources.
References
 1.
Shoults, R. R., Chang, S. K., Helmick, S., & Grady, W. M. (1980). A practical approach to unit commitment, economic dispatch and savings allocation for multiplearea pool operation with import/export constraints. IEEE Transactions on Power Apparatus and Systems, PAS99(2), 625–635 https://doi.org/10.1109/TPAS.1980.319654.
 2.
Quintana, V. H., Lopez, R., Romano, R., & Valadez, V. (1981). Constrained economic dispatch of multiarea systems using the dantzigwolfe decomposition principle. IEEE Transactions on Power Apparatus and Systems, PAS100(4), 2127–2137 https://doi.org/10.1109/TPAS.1981.316490.
 3.
Ouyang, Z., & Shahidehpour, S. M. (1991). Heuristic multiarea unit commitment with economic dispatch. IEE Proceedings C  Generation. Transmission and Distribution, 242–252 https://doi.org/10.1049/ipc.1991.0030.
 4.
Wang, C., & Shahidehpour, S. M. (1992). A decomposition approach to nonlinear multiarea generation scheduling with tieline constraints using expert systems. IEEE Transactions on Power Systems, 7(4), 1409–1418 https://doi.org/10.1109/59.207362.
 5.
Streiffert, D. (1995). Multiarea economic dispatch with tie line constraints. IEEE Transactions on Power Systems, 10(4), 1946–1951 https://doi.org/10.1109/59.476062.
 6.
Yalcinoz, T., & Short, M. J. (1988). Neural networks approach for solving economic dispatch problem with transmission capacity constraints. IEEE Transactions on Power Systems, 13(2), 307–313 https://doi.org/10.1109/59.667341.
 7.
Thang, N. T. (2013). Economic emission load dispatch with multiple fuel options using hopfield Lagrange network. International Journal of Advanced Science and Technology, 57, 9–24 https://doi.org/10.14257/ijhit.2015.8.1.27.
 8.
Kuo, C. C. (2008). A novel string structure for economic dispatch problems with practical constraints. Energy Conversion and Management, 49, 3571–3577 https://doi.org/10.1016/j.enconman.2008.07.007.
 9.
Li, Y. Z., Jiang, L., Wu, Q. H., et al. (2017). Windthermal power system dispatch using MLSAD model and GSOICLW algorithm. Knowledge Based Systems, 116, 94–101 https://doi.org/10.1016/j.knosys.2016.10.028.
 10.
Jayabarathi, T., Sadasivam, G., & Ramachandran, V. (2000). Evolutionary programming based multiarea economic dispatch with tie line constraints. Electric Machines and Power Systems, 28(12), 1165–1176 https://doi.org/10.1080/073135600449044.
 11.
Manoharan, P. S., Kannan, P. S., Baskar, S., & Iruthayarajan, M. (2009). Evolutionary algorithm solution and KKT based optimality verification to multiarea economic dispatch. International Journal of Electrical Power and Energy Systems, 31(78), 365–373 https://doi.org/10.1016/j.ijepes.2009.03.010.
 12.
Sharma, M., Manjaree, P., & Laxmi, S. (2011). Reserve constrained multiarea economic dispatch employing differential evolution with time varying mutation. International Journal of Electrical Power and Energy Systems, 33(3), 753–766 https://doi.org/10.1016/j.ijepes.2010.12.033.
 13.
Somasundaram, P., & Jothi Swaroopan, N. M. (2011). Fuzzified particle swarm optimization algorithm for multiarea security constrained economic dispatch. Electric Power Components and Systems, 39(10), 979–990 https://www.researchgate.net/deref/http%3A%2F%2Fdx.doi.org%2F10.1080%2F15325008.2011.552094.
 14.
Basu, M. (2013). Artificial bee colony optimization for multiarea economic dispatch. International Journal of Electrical Power and Energy Systems, 49(7), 181–187 https://doi.org/10.1016/j.ijepes.2013.01.004.
 15.
Basu, M. (2014). Teaching–learningbased optimization algorithm for multiarea economic dispatch. Energy, 68(4), 21–28 https://doi.org/10.1016/j.energy.2014.02.064.
 16.
Basu, M. (2017). Fast convergence evolutionary programming for multiarea economic dispatch. Electric Power Components and Systems, 45(15), 1629–1637 https://doi.org/10.1080/15325008.2017.1376234.
 17.
Nguyen, K. P., Dinh, N. D., & Fujita, G. (2015). Multiarea economic dispatch using hybrid cuckoo search algorithm. In Proc. 50th International Universities Power Engineering Conference, stoke on Trent, UK, (pp. 1–4) https://doi.org/10.1109/UPEC.2015.7339777.
 18.
Ghasemi, M., Aghaei, J., Akbari, E., et al. (2016). A differential evolution particle swarm optimizer for various types of multiarea economic dispatch problems. Energy, 107(7), 182–195 https://doi.org/10.1016/j.energy.2016.04.002.
 19.
Zhang, P., Ma, W., & Dong, Y. (2019). Multiarea economic dispatching using improved grasshopper optimization algorithm. Evolving Systems. https://doi.org/10.1007/s12530019093206.
 20.
ModiriDelshad, M., Aghay Kaboli, S. H., & TaslimiRenani, E. (2016). Backtracking search algorithm for solving economic dispatch problems with valvepoint effects and multiple fuel options. Energy, 116(1), 637–649 https://doi.org/10.1016/j.energy.2016.09.140.
 21.
Mokarram, M. J., Niknam, T., Aghaei, J., et al. (2019). Hybrid optimization algorithm to solve the nonconvex multiarea economic dispatch problem. IEEE Systems Journal, 13(3), 3400–3409 https://doi.org/10.1109/JSYST.2018.2889988.
 22.
Mokarram, M. J., Gitizadeh, M., Niknam, T., & Niknam, S. (2019). Robust and effective parallel process to coordinate multiarea economic dispatch (MAED) problems in the presence of uncertainty. IET Generation, Transmission and Distribution, 13(18), 4197–4205 https://doi.org/10.1049/ietgtd.2019.0319.
 23.
Qin, J., Wan, Y., Yu, X., & Kang, Y. (2020). A newton methodbased distributed algorithm for multiarea economic dispatch. IEEE Transactions on Power Systems, 35(2), 986–996 https://doi.org/10.1109/TPWRS.2019.2943344.
 24.
Mohammadian, M., Lorestani, A., & Ardehali, M. M. (2018). Optimization of single and multiareas economic dispatch problems based on evolutionary particle swarm optimization algorithm. Energy, 116, 710–724. https://doi.org/10.1016/j.energy.2018.07.167.
 25.
Suman, M., & Sakthivel, V. P. (2020). Coulomb’s and Franklin’s laws based optimization for nonconvex economic and emission dispatch problems. International Energy Journal, 20, 225–238 http://www.rericjournal.ait.ac.th/index.php/reric/article/view/2273.
 26.
Jain, M., Singh, V., & Rani, A. (2019). A novel natureinspired algorithm for optimization: Squirrel search algorithm. Swarm and Evolutionary Computation, 44(2), 148–175 https://doi.org/10.1016/j.swevo.2018.02.013.
Acknowledgements
Authors would like to thank Government College of Engineering Dharmapuri and Annamalai University Chidambaram for providing necessary research facilities.
Funding
No funding received from any agency.
Author information
Affiliations
Contributions
VPS as the corresponding author, contributed significantly to the research, writing and submission of the paper. PDS provided decisive comments and technical supervision for the algorithm. All the authors read and approved the submitted manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare that they have no competing interests.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Sakthivel, V.P., Sathya, P.D. Single and multiarea multifuel economic dispatch using a fuzzified squirrel search algorithm. Prot Control Mod Power Syst 6, 11 (2021). https://doi.org/10.1186/s4160102100188w
Received:
Accepted:
Published:
Keywords
 Fuzzy set theory
 Heuristic optimization
 Multiarea economic dispatch
 Paretooptimal front
 Squirrel search algorithm
 Tieline constraint