Optimal Scheduling of Cascaded Hydrothermal Systems Using a New Improved Particle Swarm Optimization Technique

Optimum scheduling of hydrothermal plants generation is of great importance to electric utilities. Many evolutionary techniques such as particle swarm optimization, differential evolution have been applied to solve these problems and found to perform in a better way in comparison with conventional optimization methods. But often these methods converge to a sub-optimal solution prematurely. This paper presents a new improved particle swarm optimization technique called self-organizing hierarchical particle swarm optimization technique with time-varying acceleration coefficients (SOHPSO_TVAC) for solving short-term economic generation scheduling of hydrothermal systems to avoid premature convergence. A multi-reservoir cascaded hydrothermal system with nonlinear relationship between water discharge rate, power generation and net head is considered here. The performance of the proposed method is demonstrated on two test systems comprising of hydro and thermal units. The results obtained by the proposed methods are compared with other methods. The results show that the proposed technique is capable of producing better results.


Introduction
The optimum scheduling of hydrothermal plants is one of the important planning task in power system operation.The generation scheduling problem consists of determining the optimal operation strategy for the next scheduling period, subjected to a variety of constraints.The main objective is to allocate generations of hydroelectric and thermal plants in such a way so as to minimize the total operation cost of the systems subjected to variety of constraints.With the insignificant operating cost of hydroelectric plants, the problem of minimizing the operational cost of a hydrothermal system essentially reduces to that of minimizing the fuel cost for thermal plants under the various constraints on the hydraulic and power system network.It is basically a nonlinear programming problem involving non-linear objective function and a mixture of linear and non-linear constraints.
The main constraints include the cascaded nature of the hydraulic network, the time coupling effect of the hydro sub problem where the water inflow of an earlier time interval affects the discharge capability at a later period of time, the varying hourly reservoir inflows, the physical limitations on the reservoir storage and turbine flow rate, the varying system load demand and the loading limits of both thermal and hydro plants.
The hydrothermal scheduling problem has been the subject of investigation for several decades and many methods have been applied to solve this problem.Some of these solution methods include decomposition techniques [1], dynamic programming [2], semi-definite programming [3] and concept of non-linear network flow [4].In recent times, optimal hydrothermal scheduling problems have been solved by different heuristic techniques such as genetic algorithm [5] simulated annealing technique [6], evolutionary programming technique [7] etc. Yu, et al. used particle swarm optimization technique to solve short-term hydrothermal scheduling problem [8] with an equivalent thermal unit having smooth cost functions.A modified hybrid differentia evolution technique was applied by Lakshminarasimman, et al. [9] to solve short-term hydrothermal generation scheduling problems with promising results.A comparison of particle swarm optimization and dynamic programming for large scale   hydro unit load dispatch was made by Cheng, et al. [10].Recently, Catalao, et al [11] applied mixed-integer quadratic programming method to determine scheduling of head dependent cascaded hydro systems.
Particle swarm optimization (PSO) happens to be a comparatively new combinatorial metaheuristic technique which is based on the social metaphor of bird flocking or fish schooling [12].This algorithm has come to existence in mid 90's and has gained prominence from late 90's.The PSO technique has been applied to various fields of power system optimization.Gaing used PSO to solve economic dispatch problem considering generator constraints [13].Abido proposed a revised PSO technique for optimal design of voltage stabilizer [14].Park, et al. presented a method for solving economic dispatch with non-smooth cost functions [15].A hybrid method for optimal scheduling of short-term electric power generation of cascaded hydroelectric plants based on particle swarm optimization and chance-constrained programming was presented by Jiekang, et al. [16].
A novel parameter automation strategy called self-organizing hierarchical particle swarm optimization technique with time-varying acceleration coefficients (SOH-PSO_TVAC) is applied in this paper for the hydrothermal scheduling to address the problem of premature convergence.In this case, the particle velocities are reinitialized whenever the population stagnates at local optima during the search.A relatively high value of the cognitive component results in excessive wandering of particles while a higher value of the social component causes premature convergence of particles.Hence, time-varying acceleration coefficients (TVAC) [17] are employed to strike a proper balance between the cognitive and social component during the search.The proposed approach was first tested on a simple test system comprising of one equivalent thermal unit and four cascaded hydro unit and then the effectiveness of the SOHPSO_TVAC was demonstrated on a more practical system comprising of six thermal units and four cascaded hydro units.The results have been compared with other evolutionary methods and found to be superior.

Problem Formulation
Economic generation scheduling of hydrothermal systems involves the optimization of a problem with nonlinear objective function subject to a mixture of linear, non-linear constraints.As the fuel cost of hydroelectric plants is insignificant in comparison with that of thermal power plants, the objective is to minimize the fuel cost of thermal power plants, while making use of the availability of hydro-resources as much as possible.The objective function and associated constraints are described as follows: where,

 
sit F P is the total fuel cost, is the number of time interval for scheduling horizon, T s N is the number of thermal plants and sit is the power generation by the i-th thermal plants at time t.P Conventionally, the fuel cost curve for any unit can be represented by segments of quadratic functions of the active power output of the generator and can be expressed as where, si , a si b , si c : fuel cost coefficients of the i-th thermal unit.
For more practical and accurate modeling of fuel cost function, the above expression needs to be modified suitably.Modern thermal power plants comprise of generating units having multi-valve steam turbines in order to incorporate flexible operational facilities.The generating units with multi-valve turbines have very different cost curve compared with that defined by (2).The effect of valve-point effect loading may be considered by adding a sinusoidal function [9] to the quadratic cost function described above.Hence, the function described by ( 2 where f P is the fuel cost function of thermal units including the valve point loading effect and , si s i e f are fuel cost coefficients of the i-th thermal generating unit reflecting the valve-point effect. The above objective function described by (3) is to be minimized subject to a variety of constraints as follows: 1) Demand constraints The total power generated must balance the power demand plus losses, at each time interval over the entire scheduling period where hjt is the power generation of jth hydro generating unit at time t, P Lt is power demand at time t and P Lt P is total transmission loss at the corresponding time.
The hydropower generation is a function of water discharge rate and reservoir storage volume, which can be described by (5) as follow: hjt j hjt j hjt j hjt hjt j hjt j hjt j where where is natural inflow of j-th hydro reservoir at time t, hjt is spillage discharge rate of j-th hydro generating unit at time t, mj hjt I S  is the water transport delay from reservoir m to j and uj is the number of upstream hydro generating plants immediately above the j-th reservoir.hjt hj where , are the minimum and maximum storage volume of j th reservoir.hjt hj where, and are the minimum and maximum water discharge rate of the j-th reservoir respectively.

Overview of Some PSO Strategies
There are several variants of PSO.Some of the commonly used PSo are discussed in the following sections.

Classical PSO
The Particle Swarm Optimization (PSO) is one of the recent developments in the category of heuristic optimization technique.Kennedy and Eberhart [12] originally developed the PSO concept based on the behavior of individuals (i.e.particles or agents) of a swarm or group.PSO, as an optimization tool, provides a populationbased search procedure in which individuals called agents or particles change their position with time.In a PSO algorithm, the particles fly around the multidimensional search space in order to find the optimum solution.Each particle adjusts its position according to its own Let in a physical d-dimensional search space, the posi experience and the experience of neighboring particle.
tion and velocity of the i-th particle (i.e.i-th individual in the population of particles) be represented as the vectors is the number of members or elements in a particle, k id is the velocity of individual i at iteration k, w is the eight parameter or swarm inertia, 1 C and 2 are the acceleration constants and The updated velocity can be used to ch is for ange the position of each particle in the swarm as depicted in (12) as: t w provides a Suitable selection of inertia weigh balance between global and local explora ons, thus requiring less iteration on average to find a sufficiently optimal solution.In general, the inertia weight w is set according to the following equation: where is the maximum number of iterations and T he weighting of th max iter the cu iter is rrent number of iterations. he constants 1 C and 2 C represent t e stochastic acc ration rms that pull each particle toward the pbest and ele te gbest positions.

It is
the search toward the opti-

Coefficients (TVAC)
observed from (11) that mum solution is heavily dependent on the two stochastic acceleration components (i.e. the cognitive component and the social component).Thus, it is very important to control these two components properly in order to get optimum solution efficiently and accurately.It has been reported [18] that in PSO, problem-based tuning of parameters is a key factor to find the optimum solution accurately and efficiently.Kennedy and Eberhart [12] re-ed optimization meth ported that a relatively higher value of the cognitive component, compared with the social component, results in excessive roaming of individuals through a larger search space.On the other hand, a relatively high value of the social component may lead particles to rush toward a local optimum prematurely.
In general, for any population-bas od like PSO, it is always desired to encourage the individuals to wander through the entire search space, during the initial part of the search, without clustering around local optima.In contrast, during the latter stages, it is desirable to enhance convergence towards the global optima so that optimum solution can be achieved efficiently.For this, the concept of parameter automation strategy called Time Varying acceleration Coefficients (TVAC) had been introduced [17].The main purpose of this concept is to enhance the global search capability during the early part of the optimization process and to promote the particles to converge toward the global optimum at the end of the search.For this, the cognitive component should be reduced while the social component should be increased during search procedure.In TVAC, this can be achieved by changing the acceleration coefficients with time.With a large cognitive component and small social component at the beginning, the particles are encouraged to move around the search space, instead of moving towards the population best prematurely.On the other hand, during the latter stage of optimization, a small cognitive component and a large social component encourage the particles to converge towards the global optimum.The concept of time varying acceleration coefficients (TVAC) can be introduced mathematically as follows [17].
  where f are constants representing

ierarchical PSO with
It is ither based on a con- ognitive and social acceleration factors respectively.

Self-Organizing H TVAC (SOHPSO_TVAC)
seen that the classical PSO is e stant inertia weight factor or a linearly varying inertia weight factor.It is reported that the particles in classical PSO may converge to a local minimum prematurely due to lack of diversity in the population, particularly for complex problems [17].In SOHPSO_TVAC, the previous velocity term in ( 11) is kept at zero.It is observed that in the absence of previous velocity term the particles rapidly rush towards a local optimum solution and then stagnate due to the lack of momentum.In fact in the absence of velocity term, the optimum solution depends highly on the initial population.To overcome this difficulty, the modulus of velocity vector of a particle is reinitialized with a random velocity (called reinitialization velocity) whenever it stagnates in the search space.Stagnation of particles highly influences the performance of PSO in searching global optimum.When a particle is stagnated, it's pbest remains unchanged over a large number of iterat hen more particles are stagnated, the gbest also remains unchanged and the PSO algorithm con to a local minimum prematurely.The necessary momentum is imparted to the particles by reinitialization of velocity vector with a random velocity.The above method can be implemented as follows [17]: If Thus a series of part arm optimizers are at bers ge- , , , rand rand rand rand and rand are num nerated randomly between 0 and 1.A time varying reinitialization strategy is used to overcome the difficulties of selecting appropriate reinitialization velocities.

Development of the Proposed Algori
In this section, an algorithm based on SOHPSO_TVA is described to obtain quality solutions for scheduling problems of hydrothermal systems with cascaded reservoirs.For any population based evolutionary algorithm like PSO, the representation of individuals and their elements is very important.For the present problem, the position of each particle (i.e. each individual in the population of particles) is composed of a set of elements and for the present problem it is the discharge rate of each hydro plant and the power generated by each thermal unit.The algorithm starts with the initialization process.Let tialized randomly satisfying the constraints defined by ( 6) and (10) and can be represented by with   and e elem and , , , hjt Q and (0) sit  must satisfy the water discharge rate and therm enerating capacity constraints as depicted in ( 6) and ( 10) respectively.Assuming the spillage in ( 8) to be zero for simplicity, the water discharge rate of the j-th hydro plant in the dependent interval is then calculated using ( 8) to meet exactly the restrictions on the initial and final reservoir storage.The dependent water discharge rate must satisfy the constraints in (10).At the same time, to meet exactly the power balance constraints, the thermal generation of the dependent thermal generating unit is calculated using (4).Thus, the initial generation is checked against all the constraints.If the constraints are satisfied then movement towards the next step is undertaken.Now, the algorithm can be described as follows: Step 1: Initialize randomly e the al g ach particle according to th value acco fitness value is better than the best fitness va particle ac e of ons reaches the maxim e latest e limit of each unit including individual dimensions, searching points and velocities according to (18).These initial particles must be feasible candidates for solutions that satisfy the practical operating constraints.
Step3: If the lue in history, set current value as the pbest .
Step 4: Modify the member velocity ch of ea cording to (16) and reinitialize it according to (17).
Step 5: Choose the particle with the best fitness valu all the particles as the gbest .
Step 6: If the number rati of ite um, then go to Step 7 else go to Step 4.
Step 7: The individual that generates th gbest is mplemented using in house em-I has been initially applied to a test valuate the performance of the pr performance of PSO algorithm is quite sensitive to the solution of the problem.

Simulation Results
The proposed algorithm was i Matlab code on 3.0 MHz, 2.0 GB RAM PC.It was applied on two illustrative test systems to obtain the simulation results.

Test Syst
The proposed method system consisting of four cascaded hydro units and an equivalent thermal plant.This has been done to e oposed method in comparison to other population based methods [8,9], on the same test system.The scheduling period has been kept to 24 hours with one hour time interval.The water transport delay between connected reservoirs is also considered.For a direct comparison, effect of valve point loading is not considered in this case.The hydraulic network is shown in Figure 1.The load demand, hydro unit power generation coefficients, reservoir inflows, reservoir limits, the generation limits, cost coefficients of thermal unit and other data are taken from [8,9] and hence they are not shown in this paper.
The the various parameter settings.Tuning of parameters is essential in all PSO based methods.Based on empirical studies on a number of mathematical benchmark functions [17] it has been reported the best range of variation as 2.5 -0.5 for 1 C and 0.5 -2.5 for 2 C .The idea is to use a high initial value of the cognitive coefficient to make use of full range of the search space and to avoid premature convergence with a low social coefficient.We experimented with the same range and several combinations of the values of 1 C and 2 C were tested.The best results were obtained for 2.5 -1.2 for 1 C and 0.8 -2.5 for 2 C out of 50 trial runs.The optimiza n is done with a randomly initialized population of 30 swarms.The maximum iteration was set at 500.shows the optimal water discharge the proposed method.The optimal hydrothermal generation schedule along with demand for 24 hours is shown in Table 2. From the Table 2, it is clearly seen that total demand is met by the total power (both hydro and thermal power) generated in every scheduling interval.The proposed method converges to the optimal solution in 6.65 seconds.The optimal cost is found to be $ 922018.24.The results of the proposed method are compared with the results obtained by various versions of particle swarm optimization (PSO) techniques and genetic algorithm [8], modified hybrid differential evolution [9], improved PSO technique [19] and are shown in Table 3.It is clearly seen that the proposed method based SOHPSO_TVAC yields comparable results.

Test System-II
To evaluate the perform based on SOHPSO_TVAC, it was further applied to a test system that consists of a multi-chain cascade of four hydro units and six thermal units.The effect of valve point loading has been taken into account in this case to illustrate the robustness of the proposed algorithm.The scheduling period has been kept to 24 hours with one hour time interval.The water transport delay between connected reservoirs is also considered.The hydro subsystem configuration, hydro unit power generation coefficients, reservoir inflows, reservoir limits and other data related to hydro sub-system are same as that of test system-I [8,9].The load demand, generation limits and cost coefficients of six thermal units are given in Tables 4  and 5 respectively.In this case also, the parameters were selected separately.Out of 50 trial runs, best results were obtained for 2.5 -1.0 for 1 and 1.08 -2.5 for 2 .The optimization is done with a randomly initialized population of 50 swarms.The maximum iteration for this case was set at 1000.The iteration number was increased at a step of fifty and beyond 1000 no improvement in results was obtained.Hence, maximum iteration was set at this value.

C C
Optimal hourly water discharge rate obtained by the proposed algorithm for this system is shown in Table 6.Table 7 presents the optimal hydro-generation schedule including total hydropower generation while Table 8 presents the optimal thermal-generation along with total thermal-power generation.From the Tables 7 and 8, it is clearly seen that total demand is exactly met by the total power (both hydro and thermal power) generated in   [9] 921893.948 HDE [9] 922872.68 9 MDE [9] 922555.44 45 DE [9] 923574.3150 PSO [19] 922553.498.46 922018.246.65   4. Computation time for optimal solution for this case is found to be 76.25 seconds and optimal fuel cost is found to be $104232.48.The same problem was also solved by the so called classical PSO described in section 3.1.For this case optimal cost was found to be $106322.23. Figure 2 compares the convergence characteristic of fuel cost classical PSO and pro-sub-optimal solution prematurely.It also seen that the proposed SOHPSO_TVAC algorithm successfully addresses the problem of premature convergence and produces better results.

Conclusions
Optimum scheduling of hydrothermal plants generation is of great importanc parti-ration coefficients (SOHPSO_TVAC) for solvi t-term economic generation scheduling of hydro-sh thermal systems to avoid premature convergence has been proposed and successfully applied to solve daily hydrothermal scheduling problem.To evaluate the performance of the proposed algorithm, it has been applied on two test systems comprising of a multi-chain cascade of hydro units and several thermal units and results are presented.The effect of valve point loading is also considered.The results obtained by the proposed method have been compared with other evolutionary algorithms like improved PSO, GA and modified hybrid differential evolution (MHDE).It is found that proposed method SOHPSO_TVAC can address the problem of premature convergence and produce better results.

R 4 )
Reservoir storage volume constraints min max hj d gbest .The modified velocity and position of each par n be calculated using the current velocity and the distance from automatically inside the main PSO according to the behavior of the particles in the search space until the convergence criteria is satisfied.The variables be the initial popup N er o number of particles.For a system with lation of h N numb f hydro units and s N number of thermal s, position of k-th individual a population is ini-unit of Copyright © 2011 SciRes.
the i-th thermal unit at time t.The range of the elements(0)
of j-th hydro unit at time w an o -th reservoir at time r ate of j-th reservoir minimum and maximum storage volume of t min hj ma hj  : lo er d upper generation limits for j th hydro unit  , x hjt Q : water discharge rate f j t hjt V : storage volume of j-th reservoir at time t of thermal generating units h N : number of hydro generating units Copyright © 2011 SciRes.SGRE

Figure 2. Convergence Characteristics for fuel cost.
novel algorithm called self-organizing hierarchical cle swarm optimization technique with time-varying ac