Next Article in Journal
Dynamics Analysis of Misalignment and Stator Short-Circuit Coupling Fault in Electric Vehicle Range Extender
Next Article in Special Issue
A Hybrid LSTM-Based Genetic Programming Approach for Short-Term Prediction of Global Solar Radiation Using Weather Data
Previous Article in Journal
Modeling of Spiral Wound Membranes for Gas Separations—Part II: Data Reconciliation for Online Monitoring
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Non-Convex Combined Heat and Power Economic Dispatch via Improved Artificial Bee Colony Algorithm

1
Department of Electrical Engineering, Faculty of Engineering, University of Zanjan, Zanjan 4537138791, Iran
2
Faculty of Electrical and Computer Engineering, University of Tabriz, Tabriz 5157944533, Iran
3
Institute of Research and Development, Duy Tan University, Da Nang 550000, Vietnam
4
Department of Electrical Engineering, University of Bonab, Bonab 5551761167, Iran
*
Authors to whom correspondence should be addressed.
Processes 2020, 8(9), 1036; https://doi.org/10.3390/pr8091036
Submission received: 14 July 2020 / Revised: 8 August 2020 / Accepted: 17 August 2020 / Published: 25 August 2020
(This article belongs to the Special Issue Bioinspired Computation for Sustainable Energy Systems)

Abstract

:
It is well accepted that combined heat and power (CHP) generation can increase the efficiency of power and heat generation at the same time. With the increasing penetration of CHPs, determination of economic dispatch of power and heat becomes more complex and challenging. The CHP economic dispatch (CHPED) problem is a challenging optimization problem due to non-linearity and non-convexity in both objective function and constraints. Hence, in this paper a novel meta-heuristic algorithm, namely improved artificial bee colony (IABC) algorithm is proposed to solve the CHPED problem. The valve-point effects, power losses as well as the feasible operation region of CHP units are taken into account in the proposed CHPED problem model and the optimal dispatch of power/heat outputs of CHP units is determined via the proposed IABC algorithm. The proposed algorithm is applied on three test systems, in which two of them are large-scale CHPED benchmarks. The obtained results and comprehensive comparison with available methods, demonstrate the superiority of the proposed algorithm for dealing with non-convex and constrained CHPED problem.

1. Introduction

1.1. Motivation and Problem Statement

Heat is considered to be a byproduct of power generation in conventional power generation systems and when it is not fully used that results in lower efficiency. Co-generation systems or combined heat and power (CHP) generation systems use the heat from a power plant and send it around to interested consumers. Thus, co-generation plants can produce both heat and electricity with better energy efficiency and fuel usage [1]. In recent years, CHP systems have attracted more attention due to their higher efficiency (up to 85%), network loss reduction, and rapid return of investment [2,3] compared to conventional systems. The complexity of the economic dispatch problem will be increased by including the CHP systems. Hence, it is necessary to propose appropriate solution procedure to obtain optimal schedules for both heat and power.

1.2. Review of Related Works

Several methods were presented in the literature to solve the CHP economic dispatch (CHPED) problem, which can fall into two main categories: classical optimization methods and heuristic search methods. In the first category, classical Lagrangian-based approaches are used to solve the CHPED problem. In the approximated models, the CHPED can be formulated as a linear programming problem [4]; however, in the case of realistic models the problem is a non-convex and nonlinear model due to the objective function and feasible operation region of the CHP units [5]. Therefore, in the most of cases the CHPED is modeled as a non-convex problem. In these cases, the classical Lagrangian-based approaches may not be able to find the global optimal solution and hence the obtained heat-power schedules may be sub-optimal. Also, the need for gradients and in some cases Hessian matrix of the problem constraints increases the computational burden of these techniques. Some examples are dual partial-separable programming method [6], quadratic program method [6], Lagrange relaxation technique [7], and Lagrangian relaxation with surrogate sub-gradient multiplier updating technique [8]. A more recent approach is semi-definite programming method which was proposed in [9]. In the case of convex problem, this approach gives the optimal solution, and semi-definite programming relaxation of non-convex problem provides a strong calculable bound to the optimal value. A numerical procedure, which uses a direct analytical approach to solve the CHPED problem, was proposed in [10]. Another technique with branch and bound algorithm was recorded in [11] which uses the generalized reduced gradient technique. Most of these numerical methods have approximately the same solutions, but the respective computational loads and CPU times are different.
Most of the recently proposed approaches to solve the CHPED problem are heuristic search methods (second category). In [12] the genetic algorithm was applied to solve the CHPED problem. Genetic algorithm has been used in [13] for daily operation scheduling of CHP units. A solution using a selective particle swarm optimization approach was presented in [14]. In [15] an approach based on time varying acceleration coefficients particle swarm optimization (TVAC-PSO) has been proposed to deal with CHPED problem. Stochastic particle swarm optimization algorithm that takes into account random variations in power and heat demands was used in [16]. Specific evolutionary approaches used to solve the CHPED include harmony search algorithms [3,17], improved ant colony search algorithm [18], enhanced firefly algorithm [19], direct search method [20], artificial immune system [21], bee colony optimization [22], differential evolution [23], an augmented Lagrange combined with Hopfield neural network [24]. Also, [25] optimizes heat and power from CHP units and expected wind power, by stochastic particle swarm optimization approach. These heuristic search methods have the ability of well handling non-convex problems; however, due to the fact that they are population-based and they have stochastic nature, their convergence to the optimal solution is not ensured, and they may be trapped in a local optimal or even a non-optimal solution. The comprehensive review of the application of different stochastic search algorithms for solution of CHPED is provided in [26].
As was previously mentioned, the implemented classical and mathematical-based optimization methods are not efficient for solving nonlinear and non-convex optimization problems. On the other hand, the meta-heuristic algorithms can find better results in comparison with classical optimization methods in non-convex optimization problems. By investigating the literature in CHPED problem solution, it can be observed that different heuristic algorithms yield different solutions. A better solution for CHPED problem has a great economic saving in system operation cost. Hence, it is required to improve the capabilities of heuristic algorithms such that more optimal solutions (i.e., solutions with lower costs) attained for non-convex CHPED problems. It is worth mentioning that some exact gradient-based mathematical programming algorithms, such as [27,28,29,30,31], have not implemented for CHPED in the literature. Therefore, it is not possible to judge their performance in comparison with the meta-heuristic optimization algorithms, and it can be considered in future works.

1.3. Contributions

Many literature works are listed in the previous section that concentrated on the solution of the CHPED problem. Given that the CHPED problem is a non-convex and nonlinear problem of optimization, there is no mathematical or metaheuristic algorithm that can guarantee the optimal global solution to these problems. Because of high economic saving potential of better algorithms, this paper focuses on solution methodology of CHPED problem. In this study, a related problem in the literature [12,15,20,21,22,23,24] is used to compare the results obtained with the methods previously applied to CHPED problem. In this paper a method based on improved artificial bee colony (IABC) algorithm is proposed to solve the CHPED problem. ABC is a heuristic optimization technique, which is based on the intelligent search behavior of honey bee swarm. It provides a population-based search procedure in which individuals (which called foods positions) are modified by the artificial bees, which their aim is to discover the places of food sources with the highest nectar amount [32]. The main contributions of this work can be summarized as follows:
1.
Proposing an improved version of artificial bee colony algorithm for dealing with non-convex optimization problems.
2.
Studying the effectiveness and performance of the proposed algorithm using normal and large-scale test systems and benchmark functions.
3.
Implementation of the proposed algorithm on CHPED problem with different sizes and characteristics.
4.
Compared with available methods in the literature, achieving feasible and better results for large-scale CHPED test systems.

1.4. Paper Organization

The rest of this paper is organized as follows. Section 2 provides the mathematical formulation of the CHPED problem considering valve-point effects, transmission losses and regional heat dispatch. Section 3 describes the proposed IABC algorithm. Section 5 gives the step by step procedure of the proposed IABC algorithm for solving the CHPED problem. Several case studies are presented in Section 6. Finally, conclusions are given in Section 7.

2. Chp Economic Dispatch Problem Formulation

The considered co-generation system in this study consists of power-only units, heat-only units and CHP units. The objective of the CHPED problem is to minimize total cost of serving the heat and power demands. The total cost can be stated as sum of the costs of generating heat and power as follows [20].
min T C = e = 1 N p o C e ( P e p o ) + c = 1 N c h p C c ( P c c h p , H c c h p ) + h = 1 N h o C h ( H h h o ) $ / h
The cost functions of the different units can be expressed using the following quadratic functions.
C e ( P e p o ) = a e p o ( P e p o ) 2 + b e p o P e p o + c e p o $ / h
C c ( P c c h p , H c c h p ) = a c c h p ( P c c h p ) 2 + b c c h p P c c h p + c c c h p + d c c h p ( H c c h p ) 2 + e c c h p H c c h p + f c c h p H c c h p P c c h p $ / h
C h ( H h h o ) = a h h o ( H h h o ) 2 + b h o h H h h o + c h $ / h
The quadratic function approximation (2) is widely used in the literature for modelling the cost function of power-only units [17]. Usually, an absolute sinusoidal term is added to the quadratic cost function for modeling the valve-point effects [33] as follows.
C e ( P e p o ) = a e p o ( P e p o ) 2 + b e p o P e p o + c e p o + | d e p o sin ( f e p o ( P e p o , m i n P e p o ) ) |
Therefore, the CHPED problem becomes non-convex. The objective function (1), should be minimized subject to the following technical constraints [19]:
e = 1 N p o P e p o + c = 1 N c h p P c c h p = P D + P L
P L = e = 1 N p o m = 1 N p o P e p o B e m P m p o + e = 1 N p o c = 1 N c h p P e p o B e c P c c h p + c = 1 N c h p n = 1 N c h p P c c h p B c n P n c h p
c = 1 N c h p H c c h p + k h = 1 N h o H h h o = H D
P e p o , m i n P e p o P e p o , m a x e = 1 , , N p o
H h h o , m i n H h h o H h h o , m a x h = 1 , , N h o
P c c h p , m i n ( H c c h p ) P c c h p P c c h p , m a x ( H c c h p ) c = 1 , , N c h p
H c c h p , m i n ( P c c h p ) H c c h p H c c h p , m a x ( P c c h p ) c = 1 , , N c
where (6) models the power production and consumption balance. The power transmission system loss is calculated using B matrix coefficients using (7). The heat production and demand balance is modeled using (8). The capacity limits of the power-only units and heat-only units are bounded using (9) and (10), respectively. The production limits of heat and power generation of CHP units are modeled using (11) and (12). It is observed from these equations that the upper and lower limits of power generation of CHP units are functions of produced heat (or vice versa). This heat-power dual dependency is presented using feasible operation region (FOR) for a specific CHP. The FOR of CHP units represents either a convex region or non-convex region as described in [34]. In the case of non-convex region, the CHPED problem becomes more complicated, due to non-convexity in both the objective function and the constraints. Typical FORs of CHP units is presented in Figure 1 and Figure 2. As it can be observed, Figure 1 represents a convex region while Figure 2 represents a non-convex region [35].

3. Improved Artificial Bee Colony Algorithm

3.1. Original Abc Algorithm

This algorithm is based on particle swarm intelligence that is inspired by the behavior of honey bees finding food. Bee colony algorithm first was proposed by Karaboga [32]. In this algorithm there are three categories of bees, i.e., employed, onlooker and scout bees. The population of employed bees and onlooker bees are equal (i.e., half of the colony).
Employed bees are responsible for exploiting the nectar sources and providing the waiting bees (onlooker bees) in the hive with information about the nature of the locations of the food source which they exploit. Onlooker bees wait in the hive and decide to exploit a food source based on knowledge exchanged by the bees they are working. Scouts either search the area randomly to find a new food source according to their internal motivation or based on potential external clues [36]. The process of finding food source in honey bee colony can be divided into three parts [32]:
1.
Employed bees discover food sources and determine the quality of nectar and share its location with others bees.
2.
Onlooker bees decide based on the quality of the food sources found by employed bees and follow the location of food sources of employed bees.
3.
If the food source of an employed bee is abandoned, it becomes scout bee and discover new food source randomly.
In ABC algorithm, each multi-dimensional particle (or food source) is shown as follows [32].
X = x 1 , . . . , x j , . . . , x D
Thus, the i-th particle is shown as follows.
X i = x i , 1 , . . . , x i , j , . . . , x i , D
Here, S N is number of artificial bees and D is the number of optimization variables (or problem dimension), i 1 , . . . , S N and j 1 , . . . , D .
That bee employed is associated with a single place of food source. Therefore, the number of places of food supply is equal to the number of bees employed. In each iteration of ABC algorithm, employed bees discover food sources as follows [32].
x i , j n e w = x i , j o l d + r a n d ( 1 , 1 ) × ( x i , j o l d x k , j o l d )
In the above equation k is a random integer that it selected from the set 1 , . . . , S N . After production of new solutions in each iteration, the fitness (nectar) function is calculated from the following expression [32].
F i t ( X i ) = 1 1 + f ( X i ) f ( X i ) 0 1 + a b s ( f ( X i ) ) f ( X i ) < 0
where f ( X i ) is the objective function value for X i to be minimized. After calculation of objective function fitness, if F i t ( X i n e w ) > F i t ( X i o l d ) then X i o l d is replaced with X i n e w .
The onlooker bees select the employed bees location based on the fitness value of their corresponding food sources. For this purpose, the possibility of choosing the food source location is calculated as follows.
P i = F i t ( X i ) k = 1 S N F i t ( X k )
As the nectar quantity of food sources (fitness of solutions) increases, so does the number of onlookers visiting them, which facilitates convergence to the optimal solution [36].
For each iteration a random real number is generated for each source within the range [0,1]. If the probability ( P i in (17) associated with this source is greater than this random number, the onlooker bee modifies the location of this source of food by using (15). After the food source is evaluated from (16), if the fitness value is improved, then the onlooker bee replaces the old food source location by the new one, otherwise it keeps the old location.
If after a certain number of iterations, employed bee’s food source location does not improved, the food source location is abandoned and this location is replaced with a random new location by the scout bee from:
x i , j S c o u t = x j min + r a n d ( 0 , 1 ) × ( x j max x j min )
where x j max and x j min are upper and lower bounds for j-th decision variable x j , respectively.

3.2. Improved Abc (Iabc) Algorithm

The ABC algorithm has been implemented successfully in various optimization problems such as in hydroelectric generation estimation [37] and parameter estimation of solar cells [38]. However, it still attracts the attention of many researchers to improve its performance. Most of these methods modify Equation (15). For example, Kraboga proposed a new search equation for employed bees as follows [36].
x i , j n e w = x i , j o l d + r a n d ( 1 , 1 ) × x i , j o l d x k , j o l d , R i , j M R x i , j o l d , O t h e r w i s e
In the above equation M R is modification rate which is equal to 0.8, R i , j is a uniformly distributed random number in the interval 0 , 1 .
Also, in [39] Gao and Liu proposed a new search equation as follows.
x i , j n e w = x j b e s t + r a n d ( 1 , 1 ) × ( x r 1 , j o l d x r 2 , j o l d )
where r 1 and r 2 are mutually different random integers selected from the set 1 , . . . , S N . x j b e s t is the individual x j corresponds to the particle with the best fitness in the current population.
The results reported using the above modifications indicate that both of the above search rules are very effective approaches in the optimization problems solved by ABC algorithm [39].
In this paper, a hybrid search technique is proposed which combines the above search formulas as follows.
x i , j n e w = x j b e s t + r a n d ( 1 , 1 ) × ( x r 1 , j o l d x r 2 , j o l d ) , R i , j M R x i , j o l d , O t h e r w i s e
The main distinguishing features of the proposed IABC algorithm are as follows:
1.
As it is observed from Equation (21), as the difference between x r 1 , j and x r 2 , j decreases, the disturbance of position x i , j decreases. Therefore, the length of step is adaptively reduced by approaching to an optimal solution, and hence the algorithm converges to the optimal solution.
2.
It is observed from Equation (18) that the algorithm automatically jumps form local optimal or even non-optimal points, since scout bees are generated when no progress made in the search for a specified food source (or solution).
3.
Onlooker bees capability included in this algorithm enables comparison of the behavior of all food sources (or solutions) simultaneously. In other words, it is observed from Equation (17) that if a specified solution (or food source) i has a small P i , then it is a good solution, and hence it is not updated by onlooker bees. Otherwise, it is replaced with new position by onlooker bees.

4. Investigation of Iabc Algorithm on the Benchmark Functions

To investigate the performance of proposed IABC algorithm, two studies are conducted here as follows.

4.1. Study-I: Investigations on Six Benchmark Functions

In this study, six well-known benchmark functions which have different characteristics are examined. These functions which have been employed in [39,40] for evaluation of ABC algorithm and its variants, are given in Table 1.
The proposed IABC algorithm is applied to the above benchmark functions. Similar to the settings of GABC algorithm given in [40], the following settings are considered for the proposed IABC algorithm evaluation: Population size (or SN) is 80, maximum iterations number ( I t e r m a x ) is 5000 and the number of trial runs is 30.
For the above settings, the mean value and standard deviation of the results are presented in Table 2. The obtained results are compared with MABC [39], ABC [40] and GABC [40]. It can be observed from Table 2 that the proposed IABC algorithm converges to better results in comparison with ABC, MABC and GABC algorithms, in terms of the mean and standard deviation of the results.

4.2. Study-Ii: Investigations on Large-Scale Benchmark Functions

To examine the capability of proposed IABC algorithm for solution of large-scale optimization problems, it is implemented on the 300 variables version of the benchmark functions f 4 f 6 [41]. Table 3 compares the mean value obtained by the proposed IABC algorithm with GSO [41], GA [41], PSO [41], EP [41], ES [41], ABC [36], and MABC [39]. It is evidently observed from this table that IABC algorithm outperforms the above existing approaches, since the obtained solution by IABC is very close to the global optimal solution in all considered benchmark functions.
Besides, in Table 4 the performance of IABC algorithm is compared with the basic ABC [36] and MABC [39] in terms of the best and worst obtained solutions for the above three benchmarks. It is observed from this table that the proposed IABC approach gives smaller values for both the best and worst solutions. especially for f 4 and f 5 these values are both zero, which means the algorithm always converges to the global optimal point.

5. Implementation of Iabc on the Chped Problem

To solve the CHPED problem by the IABC algorithm, the following steps are performed.
1.
Step-1: The first stage in IABC algorithm is initialization of the employed bees. Every food source location is a candidate solution of CHPED problem. The position of each food source ( X i ) is a vector of all real power and heat outputs of the units as presented in the following.
X i = [ X i ( 1 ) , X i ( 2 ) , X i ( 3 ) ]
X i ( 1 ) = [ P i , 1 p o , , P i , N p o p o ]
X i ( 2 ) = [ P i , 1 c h p , , P i , N c h p c h p ]
X i ( 3 ) = [ H i , 1 c h p , , H i , N c h p c h p , H i , 1 h o , , H i , N h o h o ]
The initial population of employed bees for power-only and heat-only units are determined from (26) and (27), respectively. Also, the population of employed bees for CHP units is determined from (28) and (29), respectively.
P i , j p o = P j p o , m i n + r a n d ( 0 , 1 ) × ( P j p o , m a x P j p o , m i n )
H i , j h o = H j h o , m i n + r a n d ( 0 , 1 ) × ( H j h o , m a x H j h o , m i n )
P i , j c h p = P j c h p , m i n + r a n d ( 0 , 1 ) × ( P j c h p , m a x P j c h p , m i n )
H i , j c h p = H j c h p , m i n + r a n d ( 0 , 1 ) × ( H j c h p , m a x H j c h p , m i n )
2.
Step-2: By setting I t e r = 1 (where I t e r is the iteration number of the algorithm), discover new food source locations by employed bees using Equation (21).
3.
Step-3: In this step the objective function value for the population of bees are calculated at the current iteration. Since the CHPED is a constrained optimization problem, it is converted to an unconstrained problem using penalty coefficient ( λ ). λ is assumed to be 10000 for all test systems studied in the following section. Hence, the objective function will be as follows.
f ( X i ) = e = 1 N p o C e ( P e p o ) + c = 1 N c h p C c ( P c c h p , H c c h p ) + h = 1 N h o C h ( H h h o ) + λ × ( e = 1 N p o P e p o + c = 1 N c h p P c c h p P D P L ) 2 + λ × ( c = 1 N c h p H c c h p + k h = 1 N h o H h h o H D ) 2
4.
Step-4: Fitness of i-th food source is calculated from Equation (16). If the new food source fitness is better than the old, then the old food source is replaced with the new location (obtained in Step-2, and B S c o u t , i = 0 ( B S c o u t , i is a counter that determines limit value for converting i-th employed bee to scout bee), otherwise old location is preserved and B S c o u t , i = B S c o u t , i + 1 .
5.
Step-5: At this step onlooker bees select food source of employed bees by using the roulette wheel criterion given in (17). Based on the value of P i for each food source, the onlooker bees modify the selected locations of employed bees by using (21), as follows: If P i > r a n d ( 0 , 1 ) , then the fitness of new food source is calculated from (16). If it is better than old fitness value, the old food source is replaced with new location and B S c o u t , i = 0 , otherwise the old food location is kept and B S c o u t , i = B S c o u t , i + 1 .
6.
Step-6: After the completion of the food source update process for employed and onlooker bees, if B S c o u t , i > L i m i t , then that food source is abandoned, and the employed bee is converted to a scout bee. The scout bee selects its new food source randomly in the space via (18).
7.
Step-7: Check the stopping criterion. If the algorithm converged, then go to Step-8, else I t e r = I t e r + 1 and go to Step-2 and repeat the above procedure. In this paper, the stopping criterion is reaching to the maximum number of iterations in each run. In other words, if I t e r = I t e r m a x , then the algorithm stopped.
8.
Step-8: Stop. To clarify the optimization process for energy engineers, the implemented method is presented in Figure 3.

6. Case Studies

In this section, the effectiveness and validity of the proposed method is evaluated by implementing it on three different test systems. The numerical study is performed using MATLAB 7.5 software on a PC with an Intel Core i7, 2.93 GHz CPU and 8 GB of RAM. The obtained results using the proposed IABC are compared with the reported results in the literature. The parameters used in the algorithm for different case studies are presented in Table 5.

6.1. Test System I (7-Unit System)

The first test system consists of four power-only units (units 1–4), two CHP units (units 5–6) and a heat-only unit (unit 7). The valve-point effects and transmission losses are considered in this test system. The cost function parameters of this case along with the feasible region coordinates of CHP units are available in [23]. Data of test system I is provided in Table 6. By investigating the published papers, it was found that there were three different loss matrix data for this system. Hence, we have solved the problem for this system using the available data in three cases.

6.1.1. Case I

In this case, the coefficients of the network loss matrix are provided as follows [15].
B = 10 7 × 49 14 15 15 20 25 14 45 16 20 18 19 15 16 39 10 12 15 15 20 10 40 14 11 20 18 12 14 35 17 25 19 15 11 17 39
The optimal dispatches of the units in this case are provided in Table 7. The obtained results using the proposed IABC algorithm are compared with the results of line-up competition algorithm (LCA) [42], teaching learning-based optimization (TLBO) [42], oppositional teaching learning-based optimization (TLBO) [42], conventional PSO (CPSO) [15] and time varying acceleration coefficients PSO (TVAC-PSO) [15] in Table 7. The distribution of the total cost for 50 independent runs is presented in Figure 4. It is inferred from this figure that in 33 runs the obtained value for cost is less than the mean value, which means the algorithm is capable of attaining a solution better that than the mean value in 66% of trails. The minimum, average and maximum values of the obtained costs for these runs are also provided in Table 7. The convergence characteristics of the proposed method for this case is depicted in Figure 5. It is observed that the proposed IABC algorithm converges quickly in early iterations and hence, the number of maximum runs can be decreased to save the solution time.

6.1.2. Case II

In this case, the coefficients of the network loss matrix are assumed to be as follows [44].
B = 10 6 × 49 14 15 15 20 25 14 45 16 20 18 19 15 16 39 10 12 15 15 20 10 40 14 11 20 18 12 14 35 17 25 19 15 11 17 39 B 0 = 10 3 × 0.3908 0.1297 0.7047 0.0591 0.2161 0.6635 B 00 = 0.056
The optimal solution results using the proposed algorithm are presented in Table 8 and is compared with the real coded genetic algorithm (RCGA) [44] and Grey wolf optimization algorithm (GWO) [45]. It can be observed that the average of the obtained costs using the proposed algorithm for 50 independent runs is lower than the minimum cost obtained using RCGA.

6.1.3. Case III

In this case, the coefficients of the network loss matrix are provided as follows [22].
B = 10 6 × 49 14 15 15 20 25 14 45 16 20 18 19 15 16 39 10 12 15 15 20 10 40 14 11 20 18 12 14 35 17 25 19 15 11 17 39
Using the above data, the problem is solved using the proposed algorithm and the results are presented in Table 9. The obtained results are compared with particle swarm optimization (PSO) [22], real-coded genetic algorithm (RCGA) [22], evolutionary programming (EP) [22], artificial immune system (AIS) [21], bee colony optimization (BCO) [22] and differential evolution (DE) [23] in Table 9. It can be found that the proposed algorithm outperforms all of the previously proposed algorithms in the literature in less time.

6.2. Test System Ii (24-Unit System)

This system is one of the large benchmarks for the CHPED problem, which is proposed in [15]. Data of cost functions of the second test system is presented in Table 10. The valve-point effects are considered in this test system. This test system consists of 13 power only units, 6 CHP units, and 5 heat-only units. The total electricity demand is 2350 MW and the total heat demand is 1250 MWth. The power only units’ data are based on the 13-unit challenging standard economic dispatch test system [46]. Data of the test system can be accessed from [15].
Table 11 shows the obtained results using the proposed algorithm for this system. The obtained results in this table are compared with the recent algorithms such as TLBO [43], OTLBO [43], CPSO [15], TVAC-PSO [15], arithmetic crossover harmony search (ACHS) [47], group search optimization method (GSO) [48], improved GSO (IGSO) [48] and Grey wolf optimization algorithm (GWO) [45]. As it can be observed from this table, the proposed IABC algorithm obtains a better solution compared to other reported algorithms in the literature. It should be mentioned that the minimum, average, and maximum cost of 50 independent runs are also presented in Table 11. Distribution of total costs in these trail runs for this test system is depicted in Figure 6. It is inferred from this figure that in 27 runs the obtained cost by the proposed IABC algorithm is less than the average value of 50 trials. This means that the algorithm is able to find a solution better that than the mean value in 54% of trails. The convergence characteristics of the proposed IABC algorithm for 24-unit test system is provided in Figure 7. As it can be observed from this figure the proposed algorithm is converged to the optimal solution in earlier iterations.

6.3. Test System Iii (48-Unit System)

This system consists of 48 units which includes 26 power-only, 12 CHP, and 10 heat-only units. The data of this system is given in [15]. The total electricity demand is 4700 MW and the total heat demand is 2500 MWth. Table 12 summarizes the optimal heat and power dispatches using the proposed method. The obtained results using the proposed algorithm are compared with the recent algorithms such as TLBO [43], OTLBO [43], CPSO [15], TVAC-PSO [15], group search optimization method (GSO) [48] and improved GSO (IGSO) [48].
It should be mentioned that the total costs presented in this table are directly quoted from the corresponding references. Some of the algorithms have obtained a lower cost, but with the expense of violating some of the technical constraints. Table 13 presents the feasibility analysis of the different algorithms. As it can be observed from this table, the solutions obtained by TLBO, OTLBO, GSO, and IGSO are not feasible. The CHP units which violate their corresponding feasible operation regions are indicate in Table 13 by × symbol. For example, for the TLBO method the CHP unit 38 violates its feasible region and the total power and heat mismatch is 69.9975 . The minimum, mean and maximum obtained total costs are presented for 50 independent runs. Distribution of total costs for 50 independent runs for this test system is presented in Figure 8. It is observed from this figure that in 28 runs the obtained solution is better than the average solution, and if we ignore solution 34, the remaining solutions are also very close to the average solution. This shows good diversity of the obtained solutions by this algorithm. The convergence characteristics of the proposed algorithm for this system is depicted in Figure 9. As it can be observed from this figure the proposed algorithm is converged to the optimal solution in earlier iterations, similar to previous test cases, which confirms the capability of the algorithm in dealing with large scale CHPED problems. To calculate the annual cost-saving, the studied hourly load is considered to be the average load during a year. Therefore, in the 48-unit test system, as an example, the saving will be 6,082,861 ($/year) if the result of the proposed algorithm is compared with the result of the best available feasible algorithm (TVAC-PSO). This comparison shows that the proposed optimization algorithm reduces the system operation cost considerably.

7. Conclusions

A new approach based on the improved artificial bee colony (IABC) algorithm is proposed in this paper for an efficient solution of CHPED problem. Different characteristics and constraints such as valve-point effect, power losses, feasible operation region of CHP units, and capacity limits of units are taken into account in the formulation. The effectiveness of the proposed IABC algorithm is verified using standard benchmark functions and statistical analysis. It is found that the proposed algorithm can find better solutions in terms of the objective function value, convergence speed and the number of solutions with lower objective function than the mean value, compared with other versions of ABC algorithm and other heuristic algorithms. Three test cases with different size and characteristics are used to evaluate the efficiency of the propose algorithm. The obtained results using the proposed IABC algorithm are compared with the most recent proposed algorithms and it is observed that the IABC converges to a feasible solution with the lower total cost in a reasonable time in comparison with the previously reported algorithms. The numerical results substantiate that:
  • The obtained results by the proposed IABC algorithm has small diversity and in most cases the algorithm converges to optimal or near optimal solutions. In other words, the variance of the obtained solutions is small.
  • The algorithm converges in relatively small number of iterations. This means that the algorithm has a good converge speed which enables it to be used in large systems.
  • In test system I, the obtained value for the objective function is less than the average value in 66% of the trial runs. This is 54% for test system II and 56% for test system III, which means that the proposed algorithm is able to attain solutions lower than the mean value, in more than half of the trials.
  • The obtained results are also feasible which indicates that the algorithm has the capability of attaining solutions which are both optimal and feasible.
The better solution results, especially in large test systems, confirms the applicability of the proposed algorithm for dealing with the real world systems. From the application perspective, the proposed method results in an hourly saving of $205.14 per hour which means $1,797,033 saving in each year for small scale 7-unit test system. The hourly saving of $694.4 is obtained for a 48-unit case, which equals more than $6 million saving annually. As a future work, the proposed method could be extended to solve unit commitment problem considering CHP units.

Author Contributions

Conceptualization, A.R. and M.J.; methodology, A.R. and B.M.-I.; software, A.R. and K.J.; validation, A.R., B.M.-I. and A.A.; formal analysis, A.R.; investigation, A.R.; resources, A.R.; data curation, M.J.; writing–original draft preparation, A.R. and M.J.; writing–review and editing, B.M.-I and A.A.; visualization, M.J.; supervision, A.R.; project administration, A.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Parameters

P D
The electric power load of the system [ M W ].
H D
The heat load of the system [ M W t h ].
P L
The electric power loss of the transmission system [ M W ].
P e p o , m a x / m i n
Maximum/minimum generated electric power of eth power-only unit [ M W ].
H h h o , m a x / m i n
Maximum/minimum generated heat of hth boiler [ M W t h ].
P c c h p , m a x ( H c c h p )
Maximum power output of cth CHP unit in M W when generating ( H c c h p ) MWth heat.
P c c h p , m i n ( H c c h p )
Minimum power output of cth CHP unit in M W when generating ( H c c h p ) MWth heat.
H c c h p , m a x ( P c c h p )
Maximum heat output of cth CHP unit in M W t h when generating ( P c c h p ) MW power.
H c c h p , m i n ( P c c h p )
Minimum heat output of cth CHP unit in M W t h when generating ( P c c h p ) MW power.
a e p o
Quadratic cost coefficient of power-only unit e [ $ / ( M W ) 2 h ].
b e p o
Linear cost coefficient of power-only unit e [$/MWh].
c e p o
No-load cost coefficient of power-only unit e [$/h].
d e p o
Magnitude of sinusoidal term in cost function of power-only unit e [$/h].
f e p o
Frequency of sinusoidal term in cost function of power-only unit e [rad/MWh].
a h h o
Quadratic cost coefficient of heat-only unit h $ / ( ( M W t h ) h ) 2 h .
b h h o
Linear cost coefficient of heat-only unit h [$/(MWth)h].
c h h o
No-load cost coefficient of heat-only unit h [$/h].
a c c h p
Quadratic cost coefficient of CHP unit c [ $ / ( M W ) 2 h ].
b c c h p
Linear cost coefficient of CHP unit c [$/MWh].
c c c h p
No-load cost coefficient of CHP unit c [ $ / h ].
d c c h p
Quadratic cost coefficient of CHP unit c [ $ / ( M W t h ) 2 h ].
e c c h p
Linear cost coefficient of CHP unit c [$/(MWth)h].
f c c h p
Quadratic cost coefficient of CHP unit c [ $ / ( M W ) ( M W t h ) h ].

Continuous Variables

P e p o
Power output of ith power-only unit at time t.
P c c h p
Power output of jth CHP unit at time t.
H c c h p
Heat output of jth CHP unit at time t.
H h h o
Heat output of kth heat-only unit at time t.

Functions

C e ( P e p o )
Cost function of e-th power-only unit.
C c ( P c c h p , H c c h p )
Cost function of c-th CHP unit.
C h ( H h h o )
Cost function of h-th heat-only unit.

References

  1. Kazda, K.; Li, X. A Critical Review of the Modeling and Optimization of Combined Heat and Power Dispatch. Processes 2020, 8, 441. [Google Scholar] [CrossRef] [Green Version]
  2. Alipour, M.; Mohammadi-Ivatloo, B.; Zare, K. Stochastic risk-constrained short-term scheduling of industrial cogeneration systems in the presence of demand response programs. Appl. Energy 2014, 136, 393–404. [Google Scholar] [CrossRef]
  3. Nazari-Heris, M.; Mohammadi-Ivatloo, B.; Asadi, S.; Geem, Z.W. Large-scale combined heat and power economic dispatch using a novel multi-player harmony search method. Appl. Therm. Eng. 2019, 154, 493–504. [Google Scholar] [CrossRef]
  4. Rong, A.; Hakonen, H.; Lahdelma, R. An efficient linear model and optimisation algorithm for multi-site combined heat and power production. Eur. J. Oper. Res. 2006, 168, 612–632. [Google Scholar] [CrossRef]
  5. Ghorbani, N. Combined heat and power economic dispatch using exchange market algorithm. Int. J. Electr. Power Energy Syst. 2016, 82, 58–66. [Google Scholar] [CrossRef]
  6. Rooijers, F.; van Amerongen, R. Static economic dispatch for co-generation systems. Power Syst. IEEE Trans. 1994, 9, 1392–1398. [Google Scholar] [CrossRef] [Green Version]
  7. Guo, T.; Henwood, M.I.; van Ooijen, M. An algorithm for heat and power dispatch. IEEE Trans. Power Syst. 1996, 11, 1778–1784. [Google Scholar] [CrossRef]
  8. Sashirekha, A.; Pasupuleti, J.; Moin, N.; Tan, C. Combined heat and power (CHP) economic dispatch solved using Lagrangian relaxation with surrogate subgradient multiplier updates. Int. J. Electr. Power Energy Syst. 2013, 44, 421–430. [Google Scholar] [CrossRef]
  9. Jubril, A.; Adediji, A.; Olaniyan, O. Solving the Combined Heat and Power Dispatch Problem: A Semi-definite Programming Approach. Electr. Power Components Syst. 2012, 40, 1362–1376. [Google Scholar] [CrossRef]
  10. Rao, P. Combined heat and power economic dispatch: A direct solution. Electr. Power Compon. Syst. 2006, 34, 1043–1056. [Google Scholar] [CrossRef]
  11. Geem, Z.W.; Cho, Y.H. Handling non-convex heat-power feasible region in combined heat and power economic dispatch. Int. J. Electr. Power Energy Syst. 2012, 34, 171–173. [Google Scholar] [CrossRef]
  12. Subbaraj, P.; Rengaraj, R.; Salivahanan, S. Enhancement of combined heat and power economic dispatch using self adaptive real-coded genetic algorithm. Appl. Energy 2009, 86, 915–921. [Google Scholar] [CrossRef]
  13. Gopalakrishnan, H.; Kosanovic, D. Operational planning of combined heat and power plants through genetic algorithms for mixed 0–1 nonlinear programming. Comput. Oper. Res. 2015, 56, 51–67. [Google Scholar] [CrossRef] [Green Version]
  14. Ramesh, V.; Jayabaratchi, T.; Shrivastava, N.; Baska, A. A Novel Selective Particle Swarm Optimization Approach for Combined Heat and Power Economic Dispatch. Electr. Power Compon. Syst. 2009, 37, 1231–1240. [Google Scholar] [CrossRef]
  15. Mohammadi-Ivatloo, B.; Moradi-Dalvand, M.; Rabiee, A. Combined heat and power economic dispatch problem solution using particle swarm optimization with time varying acceleration coefficients. Electr. Power Syst. Res. 2013, 95, 9–18. [Google Scholar] [CrossRef]
  16. Wang, L.; Singh, C. Stochastic combined heat and power dispatch based on multi-objective particle swarm optimization. Int. J. Electr. Power Energy Syst. 2008, 30, 226–234. [Google Scholar] [CrossRef]
  17. Khorram, E.; Jaberipour, M. Harmony search algorithm for solving combined heat and power economic dispatch problems. Energy Convers. Manag. 2011, 52, 1550–1554. [Google Scholar] [CrossRef]
  18. Song, Y.; Chou, C.; Stonham, T. Combined heat and power economic dispatch by improved ant colony search algorithm. Electr. Power Syst. Res. 1999, 52, 115–121. [Google Scholar] [CrossRef]
  19. Niknam, T.; Azizipanah-Abarghooee, R.; Roosta, A.; Amiri, B. A new multi-objective reserve constrained combined heat and power dynamic economic emission dispatch. Energy 2012, 42, 530–545. [Google Scholar] [CrossRef]
  20. Chen, C.; Lee, T.; Jan, R.; Lu, C. A novel direct search approach for combined heat and power dispatch. Int. J. Electr. Power Energy Syst. 2012, 43, 766–773. [Google Scholar] [CrossRef]
  21. Basu, M. Artificial immune system for combined heat and power economic dispatch. Int. J. Electr. Power Energy Syst. 2012, 43, 1–5. [Google Scholar] [CrossRef]
  22. Basu, M. Bee colony optimization for combined heat and power economic dispatch. Expert Syst. Appl. 2011, 38, 13527–13531. [Google Scholar] [CrossRef]
  23. Basu, M. Combined Heat and Power Economic Dispatch by Using Differential Evolution. Electr. Power Compon. Syst. 2010, 38, 996–1004. [Google Scholar] [CrossRef]
  24. Dieu, V.; Ongsakul, W. Augmented LagrangeHopfield Network for Economic Load Dispatch with Combined Heat and Power. Electr. Power Compon. Syst. 2009, 37, 1289–1304. [Google Scholar] [CrossRef]
  25. Piperagkas, G.; Anastasiadis, A.; Hatziargyriou, N. Stochastic PSO-based heat and power dispatch under environmental constraints incorporating CHP and wind power units. Electr. Power Syst. Res. 2011, 81, 209–218. [Google Scholar] [CrossRef]
  26. Nazari-Heris, M.; Mohammadi-Ivatloo, B.; Gharehpetian, G. A comprehensive review of heuristic optimization algorithms for optimal combined heat and power dispatch from economic and environmental perspectives. Renew. Sustain. Energy Rev. 2018, 81, 2128–2143. [Google Scholar] [CrossRef]
  27. Neumaier, A.; Shcherbina, O.; Huyer, W.; Vinkó, T. A comparison of complete global optimization solvers. Math. Program. 2005, 103, 335–356. [Google Scholar] [CrossRef]
  28. Bonami, P.; Biegler, L.T.; Conn, A.R.; Cornuéjols, G.; Grossmann, I.E.; Laird, C.D.; Lee, J.; Lodi, A.; Margot, F.; Sawaya, N.; et al. An algorithmic framework for convex mixed integer nonlinear programs. Discret. Optim. 2008, 5, 186–204. [Google Scholar] [CrossRef]
  29. Klanšek, U.; Pšunder, M. Solving the nonlinear transportation problem by global optimization. Transport 2010, 25, 314–324. [Google Scholar] [CrossRef] [Green Version]
  30. Lastusilta, T.; Bussieck, M.R.; Westerlund, T. An experimental study of the GAMS/AlphaECP MINLP solver. Ind. Eng. Chem. Res. 2009, 48, 7337–7345. [Google Scholar] [CrossRef]
  31. Klanšek, U. A comparison between MILP and MINLP approaches to optimal solution of Nonlinear Discrete Transportation Problem. Transport 2015, 30, 135–144. [Google Scholar] [CrossRef] [Green Version]
  32. Karaboga, D.; Basturk, B. On the performance of artificial bee colony (ABC) algorithm. Appl. Soft Comput. 2008, 8, 687–697. [Google Scholar] [CrossRef]
  33. Mohammadi-Ivatloo, B.; Rabiee, A.; Soroudi, A. Nonconvex Dynamic Economic Power Dispatch Problems Solution Using Hybrid Immune-Genetic Algorithm. IEEE Syst. J. 2013, 7, 777–785. [Google Scholar] [CrossRef] [Green Version]
  34. Alipour, M.; Zare, K.; Mohammadi-Ivatloo, B. Short-term scheduling of combined heat and power generation units in the presence of demand response programs. Energy 2014, 71, 289–301. [Google Scholar] [CrossRef]
  35. Alipour, M.; Mohammadi-Ivatloo, B.; Zare, K. Stochastic Scheduling of Renewable and CHP-Based Microgrids. Ind. Inform. IEEE Trans. 2015, 11, 1049–1058. [Google Scholar] [CrossRef]
  36. Akay, B.; Karaboga, D. A modified artificial bee colony algorithm for real-parameter optimization. Inf. Sci. 2012, 192, 120–142. [Google Scholar] [CrossRef]
  37. Uzlu, E.; Akpınar, A.; Özturk, H.T.; Nacar, S.; Kankal, M. Estimates of hydroelectric generation using neural networks with the artificial bee colony algorithm for Turkey. Energy 2014, 69, 638–647. [Google Scholar] [CrossRef]
  38. Oliva, D.; Cuevas, E.; Pajares, G. Parameter identification of solar cells using artificial bee colony optimization. Energy 2014, 72, 93–102. [Google Scholar] [CrossRef]
  39. Gao, W.F.; Liu, S.Y. A modified artificial bee colony algorithm. Comput. Oper. Res. 2012, 39, 687–697. [Google Scholar] [CrossRef]
  40. Zhu, G.; Kwong, S. Gbest-guided artificial bee colony algorithm for numerical function optimization. Appl. Math. Comput. 2010, 217, 3166–3173. [Google Scholar] [CrossRef]
  41. He, S.; Wu, Q.; Saunders, J. Group search optimizer: An optimization algorithm inspired by animal searching behavior. Evol. Comput. IEEE Trans. 2009, 13, 973–990. [Google Scholar] [CrossRef]
  42. Shi, B.; Yan, L.X.; Wu, W. Multi-objective optimization for combined heat and power economic dispatch with power transmission loss and emission reduction. Energy 2013, 56, 135–143. [Google Scholar] [CrossRef]
  43. Roy, P.K.; Paul, C.; Sultana, S. Oppositional teaching learning based optimization approach for combined heat and power dispatch. Int. J. Electr. Power Energy Syst. 2014, 57, 392–403. [Google Scholar] [CrossRef]
  44. Basu, M. Combined heat and power economic emission dispatch using nondominated sorting genetic algorithm-II. Int. J. Electr. Power Energy Syst. 2013, 53, 135–141. [Google Scholar] [CrossRef]
  45. Jayakumar, N.; Subramanian, S.; Ganesan, S.; Elanchezhian, E. Grey wolf optimization for combined heat and power dispatch with cogeneration systems. Int. J. Electr. Power Energy Syst. 2016, 74, 252–264. [Google Scholar] [CrossRef]
  46. Mohammadi-Ivatloo, B.; Rabiee, A.; Soroudi, A.; Ehsan, M. Iteration PSO with time varying acceleration coefficients for solving non-convex economic dispatch problems. Int. J. Electr. Power Energy Syst. 2012, 42, 508–516. [Google Scholar] [CrossRef]
  47. Niu, Q.; Zhang, H.; Wang, X.; Li, K.; Irwin, G.W. A hybrid harmony search with arithmetic crossover operation for economic dispatch. Int. J. Electr. Power Energy Syst. 2014, 62, 237–257. [Google Scholar] [CrossRef]
  48. Hagh, M.T.; Teimourzadeh, S.; Alipour, M.; Aliasghary, P. Improved group search optimization method for solving CHPED in large scale power systems. Energy Convers. Manag. 2014, 80, 446–456. [Google Scholar] [CrossRef]
Figure 1. Convex feasible operation region (FOR) of CHP units.
Figure 1. Convex feasible operation region (FOR) of CHP units.
Processes 08 01036 g001
Figure 2. Non-convex feasible operation region (FOR) of CHP units.
Figure 2. Non-convex feasible operation region (FOR) of CHP units.
Processes 08 01036 g002
Figure 3. Flowchart of implementation process of proposed method on CHPED problem.
Figure 3. Flowchart of implementation process of proposed method on CHPED problem.
Processes 08 01036 g003
Figure 4. Distribution of total costs for 50 independent runs for 7-unit test system (Case I).
Figure 4. Distribution of total costs for 50 independent runs for 7-unit test system (Case I).
Processes 08 01036 g004
Figure 5. Convergence characteristics of the proposed IABC algorithm for 7-unit test system (Case I).
Figure 5. Convergence characteristics of the proposed IABC algorithm for 7-unit test system (Case I).
Processes 08 01036 g005
Figure 6. Distribution of total costs for 50 independent runs for 24-unit test system.
Figure 6. Distribution of total costs for 50 independent runs for 24-unit test system.
Processes 08 01036 g006
Figure 7. Convergence characteristics of the proposed IABC algorithm for 24-unit test system.
Figure 7. Convergence characteristics of the proposed IABC algorithm for 24-unit test system.
Processes 08 01036 g007
Figure 8. Distribution of total costs for 50 independent runs for 48-unit test system.
Figure 8. Distribution of total costs for 50 independent runs for 48-unit test system.
Processes 08 01036 g008
Figure 9. Convergence characteristics of the proposed IABC algorithm for 48-unit test system.
Figure 9. Convergence characteristics of the proposed IABC algorithm for 48-unit test system.
Processes 08 01036 g009
Table 1. Studied benchmark functions.
Table 1. Studied benchmark functions.
NameFormulaD (Problem Dimension)Search SpaceGlobal Minimum
Schaffer [39] f 1 ( x ) = 0.5 + s i n 2 i = 1 n x i 2 0.50 1 + 0.001 × i = 1 n x i 2 2 30 [ 100 , 100 ] n 0
Rosenbrock [39] f 2 ( x ) = i = 1 n 1 ( 100 ( x i + 1 x i 2 ) ) 2 + ( x i 1 ) 2 30 [ 50 , 50 ] n 0
Sphere [39] f 3 ( x ) = i = 1 n x i 2 30 [ 100 , 100 ] n 0
Griewank [39] f 4 ( x ) = 1 + 1 4000 i = 1 n ( x i 100 ) 2 i = 1 n c o s ( x i 100 i ) 30 [ 600 , 600 ] n 0
Rastrigin [40] f 5 ( x ) = i = 1 n ( x i 2 10 cos ( 2 π x i ) + 10 ) 2 30 [ 5.12 , 5.12 ] n 0
Ackly [39] f 6 ( x ) = 20 exp ( 0.2 1 n i = 1 n x i 2 ) exp ( 1 n i = 1 n cos ( 2 π x i ) ) + 20 + e 30 [ 32.768 , 32.768 ] n 0
Table 2. Comparison of the obtained results for benchmark functions.
Table 2. Comparison of the obtained results for benchmark functions.
Benchmark FunctionABC [40]GABC [40]MABC [39]Proposed
NumberNameMeanSDMeanSDMeanSDMeanSD
f 1 [39]Schaffer4.47 × 10 1 2.22 × 10 2 2.81 × 10 1 9.12 × 10 2 2.56 × 10 1 4.65 × 10 2 2.12 × 10 1 2.23 × 10 2
f 2  [39]Rosenbrock3.65 × 10 1 5.04 × 10 1 7.93 × 10 1 1.36 × 10 0 1.73 × 10 1 1.61 × 10 1 1.05 × 10 1 1.45 × 10 1
f 3 [39]Sphere6.38 × 10 16 1.20 × 10 16 4.17 × 10 16 7.36 × 10 17 9.43 × 10 32 6.67 × 10 32 3.21 × 10 35 1.30 × 10 34
f 4 [39]Griewank1.27 × 10 15 1.46 × 10 15 2.96 × 10 17 4.99 × 10 17 0.000.000.000.00
f 5 [40]Rastrigin1.35 × 10 13 7.97 × 10 14 1.32 × 10 14 2.44 × 10 14 0.000.000.000.00
f 6 [39]Ackly4.70 × 10 14 5.95 × 10 15 3.21 × 10 14 3.25 × 10 15 2.98 × 10 14 2.26 × 10 15 2.87 × 10 14 3.65 × 10 15
Table 3. Comparison of the obtained results by IABC with other algorithms for large-scale benchmark functions.
Table 3. Comparison of the obtained results by IABC with other algorithms for large-scale benchmark functions.
#GSO [41]GA [41]PSO [41]EP [41]ES [41]ABC [36]MABC [39]Proposed
f 4 1.82 × 10 7 3.70 × 10 1 1.81 × 10 0 2.84 × 10 2 1.58 × 10 1 1.35 × 10 10 8.35 × 10 15 0
f 5 98.9121.3427.1383.3583.26.82 × 10 14 1.04 × 10 5 0
f 6 1.35 × 10 3 6.24 × 10 0 3.95 × 10 6 2.95 × 10 1 9.62 × 10 0 7.52 × 10 4 9.62 × 10 10 8.93 × 10 11
Table 4. Comparison of the obtained results for the variants of ABC algorithm in large-scale benchmark functions.
Table 4. Comparison of the obtained results for the variants of ABC algorithm in large-scale benchmark functions.
ABC [36]MABC [39]Proposed
#BestWorstBestWorstBestWorst
f 4 06.46 × 10 10 7.55 × 10 15 9.44 × 10 15 00
f 5 01.14 × 10 13 4.27 × 10 11 5.20 × 10 5 00
f 6 5.33 × 10 4 1.10 × 10 3 7.80 × 10 10 1.08 × 10 9 6.82 × 10 11 1.10 × 10 10
Table 5. Parameters of IABC for different test systems.
Table 5. Parameters of IABC for different test systems.
Test System #SNLimit Iter max
I10050300
II200502000
III2001005000
Table 6. Cost function parameters of test system I.
Table 6. Cost function parameters of test system I.
Power Only Units
Unit a e p o b e p o c e p o d e p o f e p o P e p o , m i n P e p o , m a x
10.0082251000.0421075
20.0031.8601400.0420125
30.00122.11001600.03830175
40.00121201800.03740250
CHP Units
a c c h p b c c h p c c c h p d c c h p e c c h p f c c h p feasible region coordinates [ P c c h p , H c c h p ]
50.034514.526500.034.20.031[98.8,0], [81,104.8], [215,180], [247,0]
60.04353612500.0270.60.11[44,0], [44,15.9],[40,75],[110.2,135.6], [125.8,32.4],[125.8,0]
Heat Only Units
a h h o b h h o c h h o H h h o , m i n H h h o , m a x
70.0382.010995002695.20
Table 7. Comparison of the obtained results for 7-unit test system (Case I).
Table 7. Comparison of the obtained results for 7-unit test system (Case I).
Control VariableLCA [42]OTLBO [43]TLBO [43]CPSO [15]TVAC-PSO [15]Proposed
P 1 44.281245.886045.266075.000047.338345.8514
P 2 98.544698.539898.5479112.380098.539898.5388
P 3 112.7192112.6741112.678630.0000112.6735112.6734
P 4 211.4443209.8141209.8284250.0000209.8158209.8169
P 5 93.749493.824994.412193.270192.371893.8594
P 6 40.000040.000240.006240.158540.000040.0000
H 5 29.735829.291425.836532.565537.846729.0616
H 6 74.500075.000274.997072.673874.999974.9839
H 7 45.264145.708449.166644.760637.153245.9542
Minimum Cost ($/h)10,104.3810,094.352910,094.838410,325.333910,100.316410,094.2718
Average Cost ($/h)NA10,099.405710,114.1539NANA10,095.4446
Maximum Cost ($/h)NA10,106.831410,133.6130NANA10,100.9445
CPU time (s)NA3.062.863.293.252.21
Table 8. Comparison of the obtained results for 7-unit test system (Case II).
Table 8. Comparison of the obtained results for 7-unit test system (Case II).
Control VariableRCGA [44]GWO [45]Proposed (IABC)
P 1 74.535752.807445.2848
P 2 99.351898.539898.5507
P 3 174.7196112.6735112.6845
P 4 211.017209.8158209.8439
P 5 100.936393.811593.8194
P 6 44.10364040.0000
H 5 24.367829.370429.3226
H 6 72.5277574.9944
H 7 53.105229.370445.6832
Minimum Cost ($/h)10,712.8610,111.2410,092.9593
Average Cost ($/h)NA10,194.4110,152.5012
Maximuum Cost ($/h)NA10,452.1211,547.5437
CPU time (s)20.34385.26182.21
Table 9. Comparison of the obtained results for 7-unit test system (Case III).
Table 9. Comparison of the obtained results for 7-unit test system (Case III).
Control VariableEP [22]BCO [22]AIS [21]PSO [22]DE [23]RCGA [22]Proposed
P 1 61.36143.945750.132518.462644.211874.683452.5848
P 2 95.120598.588895.5552124.260298.538397.957898.5685
P 3 99.9427112.9320110.7515112.7794112.6913167.2308112.7003
P 4 208.7319209.7719208.7688209.8158209.7741124.9079209.8723
P 5 98.800098.800098.800098.814098.821798.800893.8212
P 6 44.000044.000044.000044.010744.000044.000140.0000
H 5 18.071312.097419.424257.923612.537958.096529.3057
H 6 77.554878.023677.077732.760378.348132.411674.9573
H 7 54.373959.87953.498159.316159.113959.491945.7375
Minimum Cost ($/h)10,390.000010,317.000010,355.000010,613.000010,317.000010,667.000010,111.8592
Average Cost ($/h)NANANANANANA10,656.4161
Maximum Cost ($/h)NANANANANANA13,638.7295
CPU time (s)5.275.165.295.385.266.472.21
Table 10. Cost function parameters of test system II.
Table 10. Cost function parameters of test system II.
Power Only Units
Unit a e p o b e p o c e p o d e p o f e p o P e p o , m i n P e p o , m a x
10.000288.15503000.0350680
20.000568.13092000.0420360
30.000568.13092000.0420360
40.003247.742401500.06360180
50.003247.742401500.06360180
60.003247.742401500.06360180
70.003247.742401500.06360180
80.003247.742401500.06360180
90.003247.742401500.06360180
100.002848.61261000.08440120
110.002848.61261000.08440120
120.002848.61261000.08455120
130.002848.61261000.08455120
CHP Units
a c c h p b c c h p c c c h p d c c h p e c c h p f c c h p feasible region coordinates [ P c c h p , H c c h p ]
140.034514.526500.034.20.031[98.8,0], [81,104.8], [215,180], [247,0]
150.04353612500.0270.60.011[44,0], [44,15.9],[40,75],[110.2,135.6], [125.8,32.4],[125.8,0]
160.034514.526500.034.20.031[98.8,0], [81,104.8], [215,180], [247,0]
170.04353612500.0270.60.011[44,0], [44,15.9],[40,75],[110.2,135.6], [125.8,32.4],[125.8,0]
180.103534.526500.0252.2030.051[20,0],[10,40], [45,55],[60,0]
190.0722015650.022.340.04[35,0],[35,20],[ 90,45],[90,25], [105,0]
Heat Only Units
a h h o b h h o c h h o H h h o , m i n H h h o , m a x
200.0382.010995002695.20
210.0382.0109950060
220.0382.0109950060
230.0523.06514800120
240.0523.06514800120
Table 11. Comparison of the obtained results for 24-unit test system.
Table 11. Comparison of the obtained results for 24-unit test system.
Control VariableTLBO [43]TVAC-PSO [15]ACHS [47]OTLBO [43]CPSO  [15]GSO [48]IGSO [48]GWO  [45]Proposed (IABC)
P 1 538.5656538.5656628.3185538.5656680627.7455628.152538.584628.3185
P 2 299.2123299.2123299.1992299.2123076.2285299.4778299.3423299.1993
P 3 299.122299.122299.199299.1220299.5794154.5535299.3423299.1993
P 4 109.992109.992109.8665109.992180159.438660.846109.9653109.8665
P 5 109.9545109.9545109.8665109.954518061.2378103.8538109.9653109.8666
P 6 110.4042110.404260110.404218060110.0552109.965360
P 7 109.8045109.8045109.8665109.8045180157.1503159.0773109.9653109.8666
P 8 109.6862109.6862109.8665109.6862180107.2654109.8258109.9653109.8665
P 9 109.8992109.8992109.8665109.8992180110.1816159.992109.9653109.8665
P 10 77.399277.39924077.399250.5304113.989441.10377.622340
P 11 77.836477.836476.950577.836450.530479.775577.705577.622376.9498
P 12 55.222555.22255555.22255591.166894.97685555
P 13 55.086155.08615555.086155115.651155.71435555
P 14 81.752481.75248181.7524117.485484.313383.953683.46581
P 15 41.761541.76154041.761545.928140404040
P 16 82.27382.2738182.273117.485481.179685.713382.773281
P 17 40.559940.55994040.559945.928140404040
P 18 10.000210.00021010.000210.001310101010
P 19 31.467931.46793531.467942.110935.0973531.456835
H 14 105.2219105.2219104.8105.2219125.2754106.6588106.4569106.0991104.8
H 15 76.520576.52057576.520580.117574.99874.9987575
H 16 105.5142105.5142104.8105.5142125.2754104.9002107.4073105.789104.8
H 17 75.483375.48337575.483380.117474.99874.9987575
H 18 39.999939.99994039.999940.000540404040
H 19 18.394418.39442018.394423.232219.73852018.378220
H 20 468.9043468.9043470.4468.9043415.9815469.3368466.2575469.7337470.3907
H 21 59.999459.99946059.99946060606060
H 22 59.999959.99996059.99996060606060
H 23 119.9854119.9854120119.9854120119.6511120120120
H 24 119.9768119.9768120119.9768120119.7176119.8823120120
Minimum cost58,006.999258,122.74657,825.436857,856.267659,736.263558,225.74558,049.019757,846.8457,825.2594
Maximum cost58,038.527358,359.552NA57,913.773160,076.690358,318.879258,219.141357,910.9857,857.1058
Mean cost58,006.999258,198.3106NA59,853.47859,853.47858,295.924358,156.519257,873.8657,836.9224
CPU time (s)5.6752.25NA5.8253.3635.5435.545.4849.98
Table 12. Comparison of the obtained results for 48-unit test system.
Table 12. Comparison of the obtained results for 48-unit test system.
Control VariableTVAC-PSO [15]CPSO [15]GSO [48]IGSO [48]TLBO [43]OLTBO [43]Proposed
P 1 538.5587359.0392627.5814629.4952538.5693628.3199628.3071
P 2 75.13474.5831302.5046151.9991225.3021225.3313224.5321
P 3 75.13474.5831225.3696299.2996229.9473223.9653224.6053
P 4 140.6146139.3803178.6488159.2254159.1352159.8516159.7442
P 5 140.6146139.3803178.2134173.6004160.0561109.915109.8049
P 6 140.6146139.3803159.884493.4383109.7821159.7795159.7348
P 7 140.6146139.3803161.4173160.773159.6609109.8946109.9910
P 8 140.6146139.3803108.776159.351159.6492109.9321110.0123
P 9 140.6146139.3803109.0234161.4184109.966159.9569159.7589
P 10 112.199874.7998115.1364115.292740.372640.89740.0033
P 11 112.199874.7998114.2308112.899477.582141.311540.0604
P 12 74.799974.7998107.283997.539492.248955.174855.1632
P 13 74.799974.799893.08115555.175592.400392.3016
P 14 269.2794679.88100448.6854448.8359359.0530
P 15 299.1993148.6585223.7257299.268149.4238225.7871224.3763
P 16 299.1993148.6585356.9056225.4102224.717375.4674.8094
P 17 140.3973139.0809109.2667162.4605109.9355160.1192159.6180
P 18 140.3973139.0809160.4169160.9664159.9052110.3532109.7450
P 19 140.3973139.0809109.6482164.0177159.7255159.819159.7410
P 20 140.3973139.0809160.0005168.4149159.782159.7765159.8047
P 21 140.3973139.0809174.5336159.540260.0777159.737159.6745
P 22 140.3973139.0809118.6394110.8099110.0689160.1751159.6378
P 23 74.799874.799840.06340.639977.681840.11440.0053
P 24 74.799874.799841.2253114.370140.270740.304240.0109
P 25 112.1997112.19935592.327592.410892.414992.1754
P 26 112.1997112.199392.04065555.095692.501292.4037
P 27 86.911992.842381.351282.182181.488285.985790.0393
P 28 56.102798.7199404044.547898.500581.0528
P 29 86.911992.842381.038381.08981.05681.719782.4319
P 30 56.102798.71994040.428191.681948.905581
P 31 10.003110.00021010.691310.54810.083210
P 32 3556.715335.273635.069652.71839.31138.8071
P 33 95.4799109.187782.8788182.152282.023698.9499
P 34 54.923565.60064040.101452.060640.110581.0677
P 35 95.4799109.188181.092282.739481.303998.9518
P 36 54.923565.600640.333640.105645.739845.6747.3001
P 37 23.498110.615810.50871010.007513.870910
P 38 54.088260.59943535.683830.033230.388135.3241
H 27 108.1177111.4458104.9965103.5903105.0678107.5951109.8506
H 28 88.9006125.689874.99874.99878.9162125.4997110.4369
H 29 108.1177111.4458104.8209104.2548104.827105.1942105.5403
H 30 88.9006125.689874.99875.3686119.600682.6853110.3925
H 31 40.001340.000140.00140.099940.234540.034639.9999
H 32 2029.870619.263619.294328.050821.956821.7102
H 33 112.926120.6188105.5564104.8032105.4339105.3622114.8715
H 34 87.882797.099774.99875.085885.4086475.0938110.4267
H 35 112.926120.6188104.8032104.8511105.7694104.9667114.8542
H 36 87.882797.0997332.329375.08679.944779.893681.2985
H 37 45.784940.263939.55144040.000141.655439.9999
H 38 28.676531.63612020.311117.740117.901820.1409
H 39 433.9113357.9456486.4858428.0157394.616445.0937399.5313
H 40 6059.99166059.506159.9359.996760
H 41 6059.99166059.920559.957859.997460
H 42 120120118.7549114.8048118.5797119.8834120
H 43 120120113.2371117.9877118.3425119.5231119.9999
H 44 415.9741370.6214212.5981535.65480.6566428.7605400.9482
H 45 6059.999959.53626059.934659.995760
H 46 6059.999959.91386059.98159.963860
H 47 119.9989119.9856113.9272107.7179117.8207119.5025120
H 48 119.9989119.9856119.2305118.6434119.1898119.444119.9999
Minimum cost117,824.8956119,708.8818117,824.896112,320.4159116,739.364116,579.239117,130.505
Maximum costNANANANA116,756.0057116,613.6505117,182.5525
Mean costNANANANA116,825.8223116,649.4473117,145.5397
CPU time (s)89.6393.3270.6570.6510.3810.9389.51
Table 13. Feasibility analysis of the obtained results for 48-unit test system.
Table 13. Feasibility analysis of the obtained results for 48-unit test system.
CHP Unit No.TVAC-PSO [15]CPSO [15]OLTLBO [43]TLBO [43]GSO [48]IGSO [48]Proposed
273333333
283333333
293333×33
303333333
313333333
323333333
333333333
343333333
353333333
3633333×3
373333333
3833××333
Total Mismatch of
Equations (6) and (8)
−2.00 × 10 4 0.0079 −3.00 × 10 4 69.9975 0.064−1.00 × 10 4 6.23 × 10 4

Share and Cite

MDPI and ACS Style

Rabiee, A.; Jamadi, M.; Mohammadi-Ivatloo, B.; Ahmadian, A. Optimal Non-Convex Combined Heat and Power Economic Dispatch via Improved Artificial Bee Colony Algorithm. Processes 2020, 8, 1036. https://doi.org/10.3390/pr8091036

AMA Style

Rabiee A, Jamadi M, Mohammadi-Ivatloo B, Ahmadian A. Optimal Non-Convex Combined Heat and Power Economic Dispatch via Improved Artificial Bee Colony Algorithm. Processes. 2020; 8(9):1036. https://doi.org/10.3390/pr8091036

Chicago/Turabian Style

Rabiee, Abbas, Mohammad Jamadi, Behnam Mohammadi-Ivatloo, and Ali Ahmadian. 2020. "Optimal Non-Convex Combined Heat and Power Economic Dispatch via Improved Artificial Bee Colony Algorithm" Processes 8, no. 9: 1036. https://doi.org/10.3390/pr8091036

APA Style

Rabiee, A., Jamadi, M., Mohammadi-Ivatloo, B., & Ahmadian, A. (2020). Optimal Non-Convex Combined Heat and Power Economic Dispatch via Improved Artificial Bee Colony Algorithm. Processes, 8(9), 1036. https://doi.org/10.3390/pr8091036

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop