Evolving ant direction differential evolution for OPF with non-smooth cost functions

https://doi.org/10.1016/j.engappai.2010.10.019Get rights and content

Abstract

In this paper, an effective and reliable algorithm, termed as evolving ant direction differential evolution (EADDE) algorithm, for solving the optimal power flow problem with non-smooth and non-convex generator fuel cost characteristics is presented. In this method, suitable mutation operator for differential evolution (DE) is found by ant colony search. The genetic algorithm evolves the ant colony parameters and the Newton–Raphson method solves the power flow problem. The proposed algorithm has been examined on the standard IEEE 30-bus and IEEE 57-bus systems with three different objective functions. Different cases were considered to investigate the robustness of the proposed method in finding the global solution. The EADDE provides better results compared to classical DE and other methods recently reported in the literature as demonstrated by simulation results.

Introduction

The optimal power flow problem (OPF) problem is an optimization tool through which the electric utilities strive to determine secure and operating conditions for a power system. The OPF solution aims to optimize a selected objective function via optimal adjustment of the power system control variables, while satisfying various equality and inequality constraints. The OPF problem, in general, is a large-scale highly constrained nonlinear non-convex optimization problem. OPF problem has been solved by conventional and evolutionary based algorithms. Many mathematical programming techniques such as linear programming (LP), nonlinear programming (NLP), quadratic programming (QP), Newton method and interior point methods (IPM) have been applied to solve the OPF problem successfully (Momoh et al., 1999, Habiabollahzadeh and Luo, 1989, Burchet et al., 1984, Mota-Palomino and Quintana, 1986, Sun et al., 1984, Yan and Quintana, 1999). Usually, these methods rely on the assumption that the fuel cost characteristic of a generating unit is a smooth, convex function. However, there are situations where it is not possible, or even appropriate, to represent the unit’s fuel cost characteristics as a convex function. This situation arises when valve-points’ units prohibited operating zones and piecewise quadratic cost characteristics are present (Sayah and Zehar, 2008).

In recent years, many heuristic algorithms, such as genetic algorithms (GA) (Lai and Ma, 1997), evolutionary programming (Yuryevich and Wong, 1999), simulated annealing (Roa-Sepulveda and Pavez-Lazo, 2003), tabu search (Abido, 2002) and particle swarm optimization (Abido, 2002) have been proposed for solving the OPF problem, without any restrictions on the shape of the cost curves. Moreover, many hybrid algorithms have been introduced to enhance the search efficiency. For instance, a hybrid tabu search and simulated annealing (TS/TA) (Ongsakul and Bhasaputra, 2002) was applied to solve the OPF with flexible alternating current transmission systems (FACTS) device problem; hybrid evolutionary programming and tabu search or improved tabu search (ITS) (Lin et al., 2002) was used to solve the economic dispatch problem with non-smooth cost functions. Meanwhile, improved evolutionary programming (IEP) (Wang et al., 2002) was successfully used to solve combinatorial optimization problems.

In the recent past, Storn and Price introduced a powerful evolutionary algorithm called differential evolution (DE) to solve numerical global optimization problems, like the OPF problem (Storn and Price, 1997). DE is a numerical optimization approach that is simple, easy to implement, faster than many evolutionary algorithms and robust. DE combines simple arithmetic operators with the classical operators of crossover, mutation and selection to evolve from a randomly generated starting population to a final solution. Each parent competes one-on-one against its corresponding offspring and the fitter of the two becomes a member of the next generation.

The DE algorithm has been successfully applied to various power system optimization problems such as generation expansion planning (Kannan et al., 2005) and hydrothermal scheduling (Lakshminarasimman and Subramanian, 2006). Figueroa and Cederio (2004) applied DE for power system state estimation. Coelho and Mariani (2006) used this algorithm for economic dispatch with valve-point effect. Basu (2008) applied DE for solving the OPF problem incorporating FACTS devices. The hybrid differential evolution (HDE) has been employed for the solution of a large capacitor placement problem (Chiou et al., 2004). The mixed integer hybrid differential evolution (MIHDE) has been employed for hydrothermal coordination (Lakshminarasimman and Subramanian, 2007a), hydrothermal optimal power flow (Lakshminarasimman and Subramanian, 2007b) and the network reconfiguration problem (Su and Lee, 2003). The variable scaling hybrid differential evolution algorithm has been used for solving network reconfiguration of distribution systems (Chiou et al., 2005).

Colorni and Dorigo Maniezzo (1992) proposed the concept of ant system (AS) and applied it to the traveling salesman problem (TSP) (Dorigo and Gambardella, 1997). The ant algorithm has been inspired by the behavior of real ant colonies, in particular, by their foraging behavior. Recently, the ant algorithm has been applied to various optimization problems, such as the short-term generation scheduling problem (Yu et al., 1998), unit commitment (Sisworahardio and El-Keib, 2002) and hydro-electric generation scheduling (Huang, 2001).

In this paper, an efficient evolving ant direction DE based approach is proposed to solve the OPF problem with non-smooth cost functions. Evolving ant direction mutation operator selection is suggested to the original DE algorithm. Though there are five mutation operations stated in this paper, the EADDE uses only one mutation operator in every generation during the solution process. The proposed EADDE method embedded with the ant colony search is able to constantly choose different but most appropriate mutation operators during the solution process to accelerate the search for the global optimum solution. In the proposed algorithm, power flow problem is solved by Newton Raphson method; the control variables for power flow are optimized (overall optimization) by the differential evolution method. The mutation operator (out of five operators) for the differential evolution is selected by ant direction search and the ant parameters are evolved by genetic algorithm method. The proposed approach has been examined and tested on IEEE 30-bus and IEEE 57-bus standard test systems with three different objective functions. Simulation results demonstrate that the EADDE algorithm is superior to the original DE algorithm and provides significantly better results compared to those reported in the literature.

This paper introduces the combined application of the evolving concept and the ant direction hybrid differential evolution for solving the OPF problem, which was not applied to OPF in the literature so far. The remainder of the paper is organized as follows: Section 2 describes the formulation of an optimal power flow problem, while Section 3 explains the standard DE approach. Section 4 then details the procedure of the proposed evolving and direction DE and Section 5 presents the results of the optimization and compares the methods used to solve the case studies of optimal power flow problems with IEEE 30-bus and IEEE 57-bus systems. Lastly, Section 6 outlines the conclusion.

Section snippets

Problem formulation

The main goal of OPF is to optimize a certain objective subject to several equality and inequality constraints. The problem can be mathematically modeled as follows:MinOF(x,u)subject tog(x,u)=0hminh(x,u)hmaxwhere vector x denotes the state variables of a power system network that contains the slack bus real power output (PG1), voltage magnitudes and phase angles of the load buses (Vi,δi) and generator reactive power outputs (QG). Vector u represents control variables that consist of real

Overview of differential evolution

The DE algorithm is a population-based algorithm like genetic algorithms using the similar operators: crossover, mutation and selection. The main difference in constructing better solutions is that DE encodes parameters as floating-point numbers and that its mutation operation does not rely on a predefined probability distribution, but on the differences of randomly sampled pairs of solutions in the population. The algorithm uses the mutation operation as a search mechanism and a selection

Proposed algorithm for OPF by evolving ant direction differential evolution

The main idea of EADDE is to use the ant colony search system to find the proper mutation operator to accelerate the search for the global solution (Wang and Liu, 2006). The optimal values of ant colony parameters are evolved by genetic algorithm (Hozefa and Botee, 1998). The EADDE is discussed in the following:

Step 1: Initialization

The initial population is generated randomly and is given by the following equation:Xi0=Xi,min+rand().(Xi,maxXi,min),i=1,,NPwhere rand() denotes a uniformly

IEEE 30-bus system

The proposed algorithm was implemented in MATLAB computing environment with Pentium-IV, 2.66 GHz computer with 512 MB RAM. The standard IEEE 30-bus test system was used to test effectiveness of EADDE approach. The test system consists of six generating units interconnected with 41 branches of a transmission network with a total load of 283.4 MW and 126.2 MVAR. The bus data and the branch data are taken from Alsac and Stott (1973). The original system has two capacitor banks installed at buses 5

Conclusion

In this paper, a novel evolving ant direction differential evolution (EADDE) approach has been presented for solving the OPF problem with non-smooth and non-convex generator fuel cost curves with different inequality constraints. The EADDE algorithm effectively solves the OPF problem in all the three cases namely quadratic fuel cost, piecewise quadratic cost and quadratic cost with valve point loading. The robustness of EADDE has been demonstrated for different control parameter settings and

References (46)

  • Abou El Ela, A.A., Abido, M.A., Spea, S.R., 2010. Optimal power flow using differential evolution algorithm, Electric...
  • H.M. Botee et al.

    Evolving ant colony optimization

    Advances in Complex Systems

    (1998)
  • R.C. Burchet et al.

    Quadratically convergent optimal power flow

    IEEE Transactions on Power Apparent Systems

    (1984)
  • L.S. Coelho et al.

    Combining of chaotic differential evolution and quadratic programming for economic dispatch optimization with valve-point effect

    IEEE Transactions on Power Systems

    (2006)
  • J.P. Chiou et al.

    Estimation of monod parameters by hybrid differential evolution

    Bioprocess and Biosystems Engineering

    (2001)
  • J.P. Chiou et al.

    Ant direction hybrid differential evolution for solving large capacitor placement problems

    IEEE Transactions on Power Systems

    (2004)
  • Ji-Pyng Chiou et al.

    Variable scaling hybrid differential evolution for solving network reconfiguration of distribution systems

    IEEE Transactions on Power Systems

    (2005)
  • Colorni, A., Dorigo, M., Maniezzo, V., 1992. An investigation of some properties of an ant algorithm. In: Proceedings...
  • M. Dorigo et al.

    Ant colonies for the traveling salesman problem

    BioSystems

    (1997)
  • N.G. Figueroa et al.

    A differential evolution solution approach for power system state estimation

  • H.Y. Fan et al.

    A trigonometric mutation operation to differential evolution

    Journal of Global Optimization

    (2003)
  • H. Habiabollahzadeh et al.

    Hydrothermal optimal power flow based on combined linear and nonlinear programming methodology

    IEEE Transactions on Power Apparent Systems

    (1989)
  • S.J. Huang

    Enhancement of hydroelectric generation scheduling using ant colony system based optimization approaches

    IEEE Transactions on Energy Conversion

    (2001)
  • Cited by (0)

    View full text