Journal of Applied Mathematics

Volume 2014 (2014), Article ID 623216, 11 pages

http://dx.doi.org/10.1155/2014/623216

Research Article

An Environmental-Economic Dispatch Method for Smart Microgrids Using VSS_QGA

Wenxia Liu,¹ Yuying Zhang,¹ Bo Zeng,¹ Shuya Niu,¹ Jianhua Zhang,¹ and Yong Xiao²

¹Transmission and Distribution System Research Institute, North China Electric Power University, Number 2 Beinong Road, Changping District, Beijing 102206, China

²Guizhou Electric Power Test Research Institute, China Southern Power Grid, Number 32 Jiefang Road, Nanming District, Guiyang 550000, China

Received 24 January 2014; Revised 16 March 2014; Accepted 23 March 2014; Published 11 May 2014

Academic Editor: Hongjie Jia

Copyright © 2014 Wenxia Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The increasing penetration of distributed generation resources demands better economic performance of microgrids under the smart-grid era. In this paper, a comprehensive environmental-economic dispatch method for smart microgrids is proposed, with the objective for minimizing the summation of generation and emission costs in the system. As the proposed model belongs to a large-scale nonlinear and nonconvex programming problem, a hybrid heuristic algorithm, named variable step-size chaotic fuzzy quantum genetic algorithm (VSS_QGA), is developed. The algorithm utilizes complementarity among multiple techniques including the variable step size optimization, the rotation mutational angle fuzzy control, and the quantum genetic algorithm and combines them so as to solve problems with superior accuracy and efficiency. The effectiveness of the proposed model is demonstrated through a case study on an actual microgrid system and the advantages in the performance of VSS_QGA is also verified through the comparison with genetic algorithm (GA), the evolutionary programming approach (EP), the quantum genetic algorithm (QGA), and the chaotic quantum genetic algorithm (CQGA).

1. Introduction

A microgrid (MG) is defined as an electrical system that includes multiple loads and distributed energy resources, which can be operated in grid-connected or islanding mode. Due to its superior economic efficiency and operational flexibility [1], MG has been widely regarded as a novel energy supply form under the smart-grid paradigm. If applied properly, MG is beneficial to network investment deferral, pollution reduction, and improvement of energy utilization efficiency and supply reliability [2, 3].

During the 13th five-year plan period, distributed generation (DG) will undergo development in China. With increasing DG penetration into MG, the issue of “how to obtain the optimal dispatch balancing the controversy between economic and environmental benefits while meeting the load demand and network security” has become an urgent problem for MG operators.

Due to the differences in the operation characteristics in bulk power systems and MGs, there are significant distinctions in their objectives and dispatching strategies. First, the dispatch of traditional grid usually focuses on economy for optimization, through dispatching different types of thermal power unit to achieve real-time balance of the network power. Contrast to the main grid, there are lots of uncontrollable renewable distributed generations (RDG) in MG. Only the conventional units and power purchasing from the main grid are control variables. Second, when multiple unfirm generators are connected, the network capacity is shared based on a particular principle of access (PoA) in order to balance economic and environmental benefits and maintain the safe and stable operation of the system. Third, in traditional power system, the single-period based optimal power flow (OPF) is widely used as a suitable tool for power dispatch problem, due to the lack of RDG units [4]. The MG integrates a large number of RDGs and requires dynamic optimal power flow (DOPF) to treat intermittent and temporal correlation of DG. Due to the difference between the dispatching of traditional grid and MG, it is of great significance to investigate the suitable dispatch method for smart MG.

Recently, a number of optimization models and algorithms for economic dispatch (ED) of MG with different types of DG units have been developed [5–7]. In [8], a penalty function-hybrid direct search method (PF-HDSM) was developed to allocate the active power output of multiarea wind-thermal generating units to minimize the total operating cost while fulfilling the equality and inequality constraints. For the same purpose, a simple and efficient harmony search (HS) method with a new pitch adjustment rule (NPAHS) was proposed for dynamic economic dispatch (DED) of electrical power systems in [6]. And a hybrid solution methodology involving modified shuffled frog leaping algorithm (MSFLA) with genetic algorithm (GA) was carried out by [7] in which the valve-point effect was considered. With the increasing public awareness of environmental protection, it has become essential to optimize not only the fuel cost but also the emission levels of harmful pollutants. In this regard, ED problems will become economic emission dispatch (EED) problems when in addition to it there are environmental concerns caused by emissions from fossil-fuel electric power plants [9]. The EED problem is a nonlinear biobjective optimization subject designed to deal with minimization of two conflicting objectives, namely, fuel cost and air pollution from thermal power plants simultaneously [8]. The consideration of environmental emissions in MG operations could render optimization being a nonlinear and nonconvex problem with plenty of local optima. It is difficult to reach the global optimal solution with conventional numerical analysis method, such as Lambda iteration and gradient methods [3]. Recently, some heuristic algorithms have been applied to EED problems in MG, such as virus optimization algorithm (VOA) [10], hybrid particle swarm optimization and gravitational search algorithm (HPSO-GSA) [11], incremental artificial bee colony algorithm with local search (IABC-LS) [4], or ε-multi-objective GA variable (εv-MOGA) [12]. However, the above algorithms lack adaptability for large scale systems, and the calculation speed and convergence accuracy still need to be improved. Therefore, developing more sophisticated algorithms for dispatch optimization of MG is necessary.

GA is a classical algorithm that is used to solve the combinational optimization problem or nonlinear optimization problem with nondifferentiable objective function or complicated constraint condition. Due to its parallel computing character and the ability of adaptive search, GA has been proposed for the solution of EED problems [13]. However, GA also shows more iterative times and slow convergence speed and easily falls into local optimum. In 1996, Ajit Narayanan and Mark Moore proposed quantum genetic algorithm (QGA) which introduced the quantum theory into genetic algorithm. Due to its stronger search capability and more efficient and simple characteristics, QGA has been proven to be a favorable technique to tackle ED/EED problems of power systems [14]. In spite of its improvements, QGA has its shortages like premature convergence and, similarly, being caught in local optimum. To improve this problem, an advanced version with superior accuracy and efficiency is developed.

An EED model for smart microgrids (SMG) is proposed in this paper, which aimed to minimize the generation and emission cost. As the proposed model belongs to a large-scale nonlinear and nonconvex programming problem, a hybrid heuristic algorithm, named variable step-size chaotic fuzzy QGA (VSS_QGA), is developed. The algorithm combines complementarity among several techniques such as the variable step size optimization, the rotation mutation angle fuzzy control, and the quantum genetic algorithm to solve problems with superior accuracy and efficiency. For that very reason, it can be immune to local optima and has an improved searching efficiency compared to cases where the techniques are each used individually. The effectiveness of the proposed model is demonstrated through a case study on an actual MG system and the advantages in the performance of VSS_QGA are also confirmed through comparison with genetic algorithm (GA), evolutionary programming approach (EP), quantum genetic algorithm (QGA), and chaotic quantum genetic algorithm (CQGA).

2. Model Formulation for EED of SMG

In general, it is assumed that the SMG under study is composed of three different types of DG, namely, wind turbines (WT), solar photovoltaic (PV), and diesel internal combustion engine (DICE). As the power output of WT and PV units is mainly dependent on the primary renewable energy supply, the goal of EED is to allocate the active power of controllable DG units (i.e., DICE) that was installed and to achieve minimizing the generation and emission cost of the entire system while meeting customer load demand and network Figure 4.

2.1. Objective Function

The objective function of the proposed EED model for SMG is given as follows:

subject to

In above equation, the objective function is composed of both expected operating cost and emission cost . is a coefficient representing the importance of environmental value and and are penalty factors of and , respectively. Different values of mean different importance between environment and economy, which then lead to different optimization results. Operating personnel can recognize flexible tradeoffs by adjusting these coefficients. is the number of time periods in a scheduling interval. With regard to the regulation of day-ahead generation scheduling, can be considered to be 24 hours. is the number of controllable DICE units in the MG; is the number of distributed renewable energy generation (DREG) in MG. and are the scheduled active power outputs from unit of DICE and unit of DERG in period , respectively; is the active power injected from main grid in period ; and are spot prices of purchasing power from DREG and main grid during period , respectively (/kWh); represents the generation cost function of DICE unit , which is calculated as follows:

In above equation, represents the fuel cost for per kilowatt-hour (/kWh) of electricity produced by DICE unit , where , , and are the coefficients of fuel cost function; is the rate for operation and maintenance of DICE unit (/kWh).

In addition, “” in (3) represents the set of emission pollutants due to generation; and are the emission cost function of DICE and main grid with respect to pollution , respectively, (kg/kWh), which is calculated as follows:

where and are emission intensity from unit of DICE and main grid (t/kWh), respectively, representing the amount of emission induced for producing per kilowatt-hour of pollutant . is the charge rate for emission of pollutant (/t).

In the objective function (1)–(4), despite of the introduction of the coefficient , the tradeoff between economic and environmental motivations can be achieved, which represents the coordinated MG operation towards the low-carbon economy.

2.2. Constraints

The optimization formula of EED in (1) is subject to the following set of equality and inequality constraints, as shown in (7)–(13).

2.2.1. Power Balance Constraint

Consider

where and are predicted value of load demand and network losses in the period ; and represent the expected power output of WG and PV in period , respectively; and are the number of WG and PV units installed in the SMG.

In the SMG, RDG units are integrated into the network based on the “fit-and-forget” strategy, which is lack of controllability and visibility. Optimizing the operation of these resources requires consideration of intertemporal linkages as well as network power flow solutions. But this only optimizes at a single point of time. DOPF is an extension of OPF to cover multiple time periods. DOPF breaks the time-horizon into time-steps and extends all OPF variables and some parameters into time series. The network must obey the typical OPF constraints such as power flow equations during each time-step independently.

2.2.2. Power Output Limits of Controllable DG Units

To guarantee safe and stable operation of a MG, the power outputs of DG units in the MG must fall into their lower and upper limits in the operation capacity, as shown in (6)-(7):

where are the upper and lower boundaries for the active power output of controllable DG , respectively; likewise, denote the upper and lower limits for its reactive power output.

2.2.3. Power Output Limits of RDG Units

As renewable energy sources generally do not provide reactive power to the system, there is only one constraint concerning the active power output for each RDG unit, expressed as follows:

where are the upper and lower limits for the active power output of RDG unit , respectively.

2.2.4. Ramp Rate Limits of Controllable DG

In actual practice, the system operator typically makes energy acquisition plan on an hourly basis. As the load on system changes in subsequent planning intervals, the power output of controllable DG units is changed to match the load. However, this change does not take place instantaneously whenever load changes. Therefore, it is necessary to incorporate the following ramp-rate limits in the EED model:

where represents the duration of time periods; , are the up and down ramp rate limits of the controllable DG , respectively.

2.2.5. Power Interaction Limits of Smart MG with Main Grid

DREG and main grid jointly participate in the optimal operation of power system. When power shortage happens in SMG, it can absorb power from main grid in real time but can be subjected to the maximum capacity constraints of circuits and the system

where are the upper and lower boundaries for the power interacted between SMG and main grid.

Furthermore, the network constraints should be also taken into consideration.

2.2.6. Thermal Limit of Feeders

The magnitude of current passing through the feeders of a system should be within their thermal limit during operation, which is expressed as

where is the permitted maximum current limit in the feeder in the time period .

2.2.7. Node Voltage Constraint

To ensure the quality of electricity supply, the voltage variation for each node in the system must satisfy the following constraint, as is presented below:

where and stand for the minimum and maximum limits with respect to the voltage variation.

3. Variable Step-Size Chaotic Fuzzy Quantum Genetic Algorithm

Chaotic quantum genetic algorithm (CQGA) is an improved version of the well-known quantum evolutionary algorithm (QEA) [15], which incorporates the global searching ability of GA, the local search ability of quantum probability technique, and the adaptability and transverse mobile ability of chaotic algorithm. Hence, CQGA is widely acknowledged as an effective tool to solve EED problem.

To construct VSS_QGA, we take into consideration the effect of initial value selection, the search step size on the convergence speed, and the calculation precision. It combines the variable step size optimal method, the rotation mutational angle fuzzy control method, and the quantum genetic algorithm to obtain the global optimal solution with a fast search speed and higher calculation accuracy. The steps for VSS_QGA in EED are described as follows.

Step 1. Initialize populations with the variable step size chaos algorithm.

( 1) Chaotic System. Chaos is a universal, nonlinear phenomenon in nature. Its behavior is complex and similar to being random. However, due to the existence of the faint positive characteristic, a chaotic system optimizes general chaos variables by using logistic mapping:

where is a chaotic attractor; is the number of iterations. When , the system gets into chaotic state, producing the chaos variables (), between .

After iterations, () are obtained in which is the variable dimension. The second carrier wave can be obtained by using the following formula:

where is the chaos variable in the smaller ergodic interval of the feasible region relative to the optimization objective function TC. is the adjustment coefficient. is an optimal solution in the initial stage of searching.

For the initialization of the chaos system, if , , , optimal solution , and the minimum of optimization objective function , the upper and lower ranges of optimization variables are and , respectively. The chaotic variable is mapped into the optimization variables space. The general method of chaos is used to search by making , , until the local optimal value does not change. At this stage, the step size in search is changed to the following equations:

where is the adjustment coefficient; is the current optimal solution.

The optimization variables are restored to be

Through the new chaotic variable ( is a smaller number), we repeat the above steps of chaotic searching and make until the minimum of optimization objective function does not change. After several iterations, the optimal variable and the corresponding optimal objective function will be obtained by this algorithm.

( 2) Quantum Bit. In CQGA, the smallest information unit is quantum bit (qubit). A qubit state can be or , which can be expressed as

where , are tow complex of probability amplitude of the corresponding state, where ( and represent the occurrence probabilities of qubit in the states of and , respectively).

( 3) Quantum Chromosome. The commonly used coding algorithms in evolutionary algorithm (EA) are binary, decimal, and symbolic coding. In the CQGA, a new coding method based on quantum bit is determined by a pair of complex number. The system with quantum bits can be expressed as follows [16]:

where (). This expression can be used to represent any linear superposition state between and .

( 4) Chaotic Quantum Populations. Using the following units of the logistic mapping produces units of chaotic variables as

where ; is the sequence number of chaotic variables. If , a set of different initial values to a given units of chaotic variables can be obtained. Chaotic variables of units can also be obtained by using (20) (). Using the first qubit of the number chaotic variable initialization populations, by making based on Step 1, we can obtain another solutions. These solutions will consist of the initial populations.

Step 2. The encoding and individual measurements to the power output of generator in populations CQGA and EA are probability algorithms. In the population,   , where is the generation of evolutionary and is the size of population; is the individual of number generation. In , () is the binary code of number unit power generation in number generation. The chromosome is defined as (21)

where represents the length of the quantum chromosome.

In the initial group if ,    and all the in have been initialized, that is, all the possible linear superposition states could emerge with the same probability.

In the process of generating solution sets from , a general solution set, , can be obtained (); for each (), whose length is , there is a binary character string (). Through the range of quantum bits and    (), it can be found that in the binary coding a number between will be randomly generated. In addition, if , 1 will be taken; otherwise, the 0 will be taken.

Step 3. Evaluate individual targets in separately.

Individual targets in are evaluated through a fitness function to ensure that the total generation cost reaches the optimal target. If a satisfactory solution is obtained, the algorithm terminates. Otherwise, it moves on to Step 4.

Step 4. Update by using an appropriate quantum rotating gate .

The QGA applies logic gates into the quantum probability amplitude to maintain diversity of the population. Therefore, updating method using a quantum gate is the key to the QGA. In principle, this updating method applies to find solutions for combinatorial optimization problems through known optimal solution. However, for practical optimization problems, in particular those multivariable continuous function optimization problems, their optimal solutions are not obtained beforehand. Therefore, quantum rotation gates of quantum logic gates are adopted for the QGA:

where is the rotation variation angle of the quantum gate.

Rotation of the control quantum gates on the rotating angle is used to update the quantum bits in this paper, as described in the following:

where is the rotation variation angle, with convergence improved by controlling the amplitude of . If the amplitude of is too high, the problem of premature convergence occurs, and will be constant. To make full use of the trend of fitness function and adjust the scale of MG automatically, needs to change accordingly; therefore, QGA must consider the function value and the rate of change near the individual chromosome in order to adopt the changes of the fitness function into the rotation angle step function, as described in the following:

where is the initial value of iteration. is the control factor. and are the minimum and maximum values of fitness function, respectively; is the gradient at point of the fitness function . and are, respectively, determined by the following equation [17]:

For a discrete case, is the difference at point of the fitness function .

By this method, the quantum rotation gate is applied to all probability amplitudes for individual targets in the population. The following equation can be obtained by using the quantum rotation gate to update the target :

where is the number generation of quantum rotation gates. is the number probability amplitude for a given object. is the number generation probability amplitude of the corresponding object.

Step 5 (perturbation). To solve the problem that CQGA can be prone to be trapped in local optimum, the population requires a perturbation. It is found that through the analysis of CQGA, it is very difficult to prevent the best individual in current generation from becoming the local extremum. Therefore, when the best individual does not change in successive generation, the algorithm is trapped into the local extreme. At this point of time, a perturbation must be applied into the population to help it escape from the local optima and start a new search. The complete process of VSS_QGA is shown in Figure 1.

Figure 1: Flow chart of the VSS_QGA.

4. Case Study

4.1. System Data

In order to verify the effectiveness of the proposed EED and VSS_QGA, a 150-node MG extracted from an actual system in China is used as the test bed, which is shown in Figure 2. The MG system is connected to the main grid via a substation at Bus 0. The scheduling duration of EED under this study is assumed to be 1 day, which is further divided into 24 time intervals. The technical parameters of different DG technologies in MG are presented in Table 1, while the information of spot prices in each day is given in Table 2.

Table 1: Technical data of distributed generation units in SMG.

Type
Operational limits/kW
Climbing Speeds/(kW·min⁻¹)
Number

Lower limit
Upper limit

WT
0
100
—
2

PV
0
20
—
10

#1 DICE
40
100
2
1

#2 DICE
80
200
1
1

Table 2: Information of daily spot prices.

Time/h
Spot price/(¥/(kW·h))

1
0.24

2
0.18

3
0.13

4
0.10

5
0.13

6
0.17

7
0.27

8
0.39

9
0.52

10
0.53

11
0.81

12
1.00

13
0.99

14
1.49

15
0.99

16
0.79

17
0.40

18
0.36

19
0.36

20
0.41

21
0.44

22
0.35

23
0.30

24
0.23

Figure 2: A real 150-node SMG system.

The prediction of power output for a single WG and PV is shown in Figure 2. Here, it is assumed that the cutin/rated/cutout speeds of WG units are  m/s,  m/s, and  m/s, respectively. The variation of wind speed is assumed to follow a Weibull distribution, with the shape parameter and the scale parameter , where is the mean value of wind speed in the site. The illumination intensity is simulated based on the Beta distribution, with the parameters and . The daily duration of system demand is shown in Figure 3. We take 0.38/kW·h, which is consistent with the average feed-in price implemented in China; the emission cost for SO₂, NOx, and CO₂ is assumed to be 6.57, 8.57, and 0.02 /kg, respectively. The upper and lower boundaries and for the power interacted between DREG and main grid are set to 0 kW and 200 kW, respectively.

Figure 3: Power forecast of renewable energy generation.

Figure 4: Power forecast of load demand.

Initial parameter setting has a significant impact on behavior of the algorithm. Therefore, appropriate parameters of VSS_QGA have to be found correctly. According to [14], a large population size could slightly improve the mean cost value, but it will increase the CPU time. Therefore, population size of VSS_QGA is set to 20. Since large angles may cause premature convergence and small angles generally produce better solution, the magnitude of rotation angle is set to 0.02 π. Besides, maximum iteration is , qubit length is 20, and mutation rate is 0.05.

4.2. Optimization Results

In order to balance the economy and environmental protection in smart MG, is set to , which means the contribution rates of both operation and emission cost to the object function are 50%. VSS_QGA is applied to solve the model of economic and environment dispatch.

Figure 5 gives the optimization results of EED in above system. As can be seen from the power output curves of WT, PV, DICE, and main grid, the wind generation and main grid are the main sources of power supply from 1 a.m. to 7 a.m. After 8 a.m, the energy contribution of PV rises significantly and reaches its peak from 12 a.m. to 3 p.m. At the same time, due to high feed-in price in these periods, DICE units replace the main grid and become the main source to MG. As evening (5 p.m.–7 p.m.) approaches, the energy mix changes again. Wind generation and main grid play dominate roles as wind is strong and the solar output is weakened due to sunset. Last, there is an obvious increase of power generation from DICE between 7 p.m. to 9 p.m., implying that these units are put into operation due to the high spot prices in these periods.

Figure 5: Environmental and economic dispatch schemes.

When is 0, the model in this paper evolves into traditional economic dispatch model with an objective to minimize the expectation operation cost. While takes a greater value, the model evolves into environmental dispatch section with an objective to minimize the emission Table 3. In order to verify that the algorithm in this paper is applicable for all kinds of scheduling problems, though defining and , economic dispatch and environmental dispatch model are built and optimized, respectively. Optimization results under different MG scheduling paradigms are as shown in Figure 6.

Table 3: Emission characteristics of both traditional power plant and distributed generation (g/kW·h).

Generation Type
SO₂

CO₂

Main grid
0.330
2.880
643.890

DG
WT
0
0
0

PV
0
0
0

DICE
0.464
4.331
232.037

(a)

(b)

(c)

(d)

Figure 6: Power sour output under two different paradigms with and .

The output of power sources under two different operation modes is shown in Figure 6. As can be seen from Figures 5 and 6, with the growing of , DERG and main grid output are on the rise, while there is significant fall of power injection from DICE. This is mainly because when , that is, only the operating cost is considered, importing energy from the main grid is preferred to the local production because the price of the former is lower in valley periods; thus, DG units offer the extra power that is required. And in the flat and peak periods, the energy contribution of DICE and DREG increases gradually due to their superiority in prices. When , that is, only emission cost is considered, the emission-free DREG units are preferred, which is consistent with our expectation. The optimization results under different MG scheduling paradigms are shown in Table 4.

Table 4: Optimization results under different MG scheduling paradigms.

Dispatch paradigm
Energy purchase cost
DG generation cost
Main grid

emission cost
DG emission cost
Total cost

Economic Dispatch
782.46
1942.77
198.72
168.59
3092.54

Environmental Dispatch
2259.33
1368.22
146.43
44.11
3818.09

Economic-Environmental Dispatch
910.46
2140.91
96.23
112.06
3259.66

As is shown in Table 4, the optimal scheme derived from the economic dispatch strategy has the minimum economic cost, which is only 2725.23 . However, due to lack of consideration on environmental burden in the optimization, the emission cost in such system is dramatically high, which poses a negative effect on the overall attractiveness of the case. When the environmental objective is solely considered, the environmental benefits of MG can be significantly improved. This leads to a drop in the emission cost from DG and the grid by 52.29 and 124.48 , which corresponds to a decrease of 26.5% and 73.84%, respectively, compared to the economic dispatch scheme.

Under the proposed EED paradigm, the environmental and economic objectives are collectively optimized. The overall cost is calculated to be 3259.66 , which is the lowest among all scenarios. Although the scheme is not the best in either economic or environmental criteria, the coordination embodied in EED scheme enables superior performance of the plan in terms of comprehensive benefits.

4.3. Discussion

In this section, the convergence and robustness of the proposed optimization method VSS_QGA is evaluated.

To achieve this end, 50 trials are performed with different initial values generated randomly. The objective function value of final schemes in these trails is recorded. Table 5 shows the best, worst, average, and standard deviation (SD), rational standard deviation (%RSD), and corrected rational standard deviation (%CRSD) of the results.

Table 5: Summary of different trials in VSS_QGA.

Parameter
Value

Best Cost (¥)
3204.42

Worst Cost (¥)
3283.97

Average Cost (¥)
3259.66

SD
533.21

%RSD
0.975

%CRSD
0.987

As can be seen from Table 5, there is little difference in the fitness values in the 50 trials and all of them are close to the average. Furthermore, it is also found that %RSD and %CRSD are 0.975% and 0.987%, respectively, which are less than 1%. The above facts demonstrate that the proposed method is not sensitive to initial solution and is able to provide valid dispatch decisions for smart MG close to the global optima in real applications. Through above comparison, the robustness of VSS_QGA is validated.

To evaluate the convergence characteristic of VSS_QGA, different optimization methods are employed to solve EED problem. The comparisons of GA, EP, QGA, CQGA, and VSS_QGA with 500 iterations are given in Figure 7.

Figure 7: Calculation results of GA, EP, QGA, CQGA, and VSS_QGA.

As observed, among all above algorithms, VSS_QGA exhibits the best performance. It is obvious that CQGA and VSS_GA have faster convergence speed compared to the rest of methods. This is because unlike the searching strategies of GA, EP, and QGA (which is typically based on the probability to exclude poor solutions by escaping local optimum traps), VSS_QGA and CQGA obtain the optimal solutions by improving chaos motion ergodicity and local search strategy, which is more suitable and efficient for solving the scheduling optimization problem. Different from CQGA, VSS_QGA takes advantage of the variable step size searching and the rotated angle of the fuzzy control strategies, which in turn leads to a superior searching efficiency and convergence performance to reach the global optimal solution.

5. Conclusions

In this study, an application of a novel VSS_QGA for EED problem in SMG is developed. The proposed EED model for SMG considers the minimization of generation and emission costs as the objective while taking different constraints such as energy balance, operation limits of power sources, and network characteristics into account. To get faster and better optimization performance, without changing the searching mechanism of the QGA, VSS_QGA introduces the variable step size optimization and the rotation mutational angle fuzzy control method into the algorithm, which improves the evolution speed and enables the optimization immune to the local optima. The effectiveness of the proposed methods is confirmed by comparing the results with the most recently reported literatures, including QGA, EP, and GA. The corresponding results show that the optimal scheduling plan obtained by the proposed EED model will produce greater economic benefit and social benefits, particularly for the large-scale applications. Furthermore, VSS_QGA has superior convergence, robustness, and less computational complexities as compared to other methods.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This work was supported by the National Key Technology R&D Program (2013BAA02B02) of China.

References

C. Chen, S. Duan, J. Yin, and T. Cai, “Energy management system of distributed generation based on power forecasting,” Transactions of China Electrotechnical Society, vol. 25, no. 3, pp. 150–156, 2010.

X. Lin, Y. Lu, and L. Wang, “New fault region location scheme in distribution system with DGs,” Transactions of China Electrotechnical Society, vol. 23, no. 11, pp. 139–165, 2008.

Y. Xu, L. Zhang, Z. Wang, and X. Li, “Algorithm of service restoration for large area blackout in distribution network with distributed generators,” Transactions of China Electrotechnical Society, vol. 25, no. 4, pp. 135–141, 2010.

D. Aydin, S. Özyön, C. Yaşar, and T. Liao, “Artificial bee colony algorithm with dynamic population size to combined economic and emission dispatch problem,” International Journal of Electrical Power & Energy Systems, vol. 54, pp. 144–153, 2014.

C. L. Chen, Z. Y. Chen, and T. Y. Lee, “Multi-area economic generation and reserve dispatch considering large-scale integration of wind power,” International Journal of Electrical Power & Energy Systems, vol. 55, pp. 171–178, 2014.

Q. Niu, H. Zhang, K. Li, and G. W. Irwin, “An efficient harmony search with new pitch adjustment for dynamic economic dispatch,” Energy, vol. 65, no. 25, p. 43, 2013.

P. Roy, P. Roy, and A. Chakrabarti, “Modified shuffled frog leaping algorithm with genetic algorithm crossover for solving economic load dispatch problem with valve-point effect,” Applied Soft Computing, vol. 13, no. 11, pp. 4244–4252, 2013.

H. T. Jadhav and R. Roy, “Gbest guided artificial bee colony algorithm for environmental/economic dispatch considering wind power,” Expert Systems with Applications, vol. 40, no. 16, pp. 6385–6399, 2013.

S. Jiang, Z. Ji, and Y. Shen, “A novel hybrid particle swarm optimization and gravitational search algorithm for solving economic emission load dispatch problems with various practical constraints,” International Journal of Electrical Power & Energy Systems, vol. 55, pp. 629–644, 2014.

Y. C. Liang and J. R. Cuevas Juarez, “A normalization method for solving the combined economic and emission dispatch problem with meta-heuristic algorithms,” International Journal of Electrical Power & Energy Systems, vol. 54, pp. 163–186, 2014.

S. Jiang, Z. Ji, and Y. Shen, “A novel hybrid particle swarm optimization and gravitational search algorithm for solving economic emission load dispatch problems with various practical constraints,” International Journal of Electrical Power & Energy Systems, vol. 55, pp. 628–644, 2014.

E. Afzalan and M. Joorabian, “Emission, reserve and economic load dispatch problem with non-smooth and non-convex cost functions using epsilon-multi-objective genetic algorithm variable,” International Journal of Electrical Power & Energy Systems, vol. 52, pp. 55–67, 2013.

C. Yaşar and S. Özyön, “Solution to scalarized environmental economic power dispatch problem by using genetic algorithm,” International Journal of Electrical Power & Energy Systems, vol. 38, no. 1, pp. 54–62, 2012.

G.-C. Liao, “Solve environmental economic dispatch of Smart MicroGrid containing distributed generation system-Using chaotic quantum genetic algorithm,” International Journal of Electrical Power & Energy Systems, vol. 43, pp. 779–787, 2012.

A. D. S. Nicolau, R. Schirru, and A. Alvarenga De Moura Meneses, “Quantum evolutionary algorithm applied to transient identification of a nuclear power plant,” Progress in Nuclear Energy, vol. 53, no. 1, pp. 86–91, 2011.

G. S. Sailesh Babu, D. Bhagwan Das, and C. Patvardhan, “Real-parameter quantum evolutionary algorithm for economic load dispatch,” IET Generation, Transmission and Distribution, vol. 2, no. 1, pp. 22–31, 2008.

A. M. Martínez, “Recognizing imprecisely localized, partially occluded, and expression variant faces from a single sample per class,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 24, no. 6, pp. 748–763, 2002.