A Refined Teaching-Learning Based Optimization Algorithm for Dynamic Economic Dispatch of Integrated Multiple Fuel and Wind Power Plants

Dynamic economic dispatch problem (DEDP) for a multiple fuel power plant is a nonlinear and nonsmooth optimization problem when valve-point effects, multifuel effects, and ramp-rate limits are considered. Additionally wind energy is also integrated with the DEDP to supply the load for effective utilization of the renewable energy. Since the wind power may not be predicted, a radial basis function network (RBFN) is presented to forecast a one-hour-ahead wind power to plan and ensure a reliable power supply. In this paper, a refined teaching-learning based optimization (TLBO) is applied to minimize the overall cost of operation of windthermal power system. The TLBO is refined by integrating the sequential quadratic programming (SQP) method to fine-tune the better solutions whenever discovered by the former method. To demonstrate the effectiveness of the proposed hybrid TLBO-SQP method, a standard DEDP and one practical DEDP with wind power forecasted are tested based on the practical information of wind speed. Simulation results validate the proposed methodology which is reasonable by ensuring quality solution throughout the scheduling horizon for secure operation of the system.


Introduction
Dynamic economic dispatch problem (DEDP) is indispensable for real-time control of power system operation in power generation system.It comprises allocating the total electricity generation required among the available thermal and other power generating units, assuming that a unit commitment has been previously determined [1].The prime objective is to minimize the cost of generated power subject to physical and operational constraints in the system.The overwhelming majority of published literature, however, deals with static economic dispatch; that is, the dispatch horizon is divided into periods and the dispatch is optimized period by period.On the contrary, a dynamic dispatch process can handle with this dynamic connection both by handling the ramp rate limits of generating units and by modifying the steady state costs to include the extra fuel consumption pertinent to the act of changing the delivered power [2].Thus, the DEDP has been acknowledged as not only a more precise formulation of economic dispatch problem (EDP) but also a complicated dynamic optimization problem [3][4][5][6][7][8][9].
As a potential renewable energy source, the wind power has been drawing enormous consideration, as various environmental and economic apprehensions have been dramatically increased in thermal power generation system [10].Integrating wind energy into existing power system introduces numerous challenges to the operation and planning strategies for the utilities.The fact that wind power is neither easily predictable nor dispatchable is the prime motive behind security and reliability concerns associated with wind-integrated power systems [11].When the wind speed changes randomly, it is more suitable for power systems containing a wind power farm to adopt dynamic models.Using DEDP, it is necessary to know the output data from the wind farm at every moment in the optimization process.Presently, it is rather difficult to forecast the output of a wind farm.The error can be large which further increases the complexity for the DEDP to optimally share the load among the available generators [7].
Subsequently, in a power generating station, generators are possibly operated with different fuels called multiple fuel plants and when such cases exist and are supplied with multiple fuel sources then that leads to an optimization problem of determining the economic fuel to burn in a particular time horizon [3,4].In the case of these generators, unlike the conventional cost function, the cost function of each unit should be presented with a few piecewise functions reflecting the effects of fuel type changes and each segment of the hybrid cost function implies some information about the type of fuel being burned or the operational characteristics of the unit.Thus, incorporating such multiple fuel plants and wind farms for DEDP makes the problem formulation practical and greatly challenging for the solution procedure to achieve a reasonable dispatch solution.Since valve-point effects are considered the fuel cost function of such generators, they will be modeled by superimposing this effect as a rectified sinusoid component into the generating unit (quadratic) fuel cost function [5].Traditional methods like Lagrangian relaxation, gradient projection method, and so forth [3], when used to solve DEDP, suffer from myopia for nonlinear, discontinuous search spaces, leading them to less than desirable performance, often using approximations to limit complexity.When search space is particularly irregular (due to inclusion of valve-point effects), algorithms need to be highly robust to escape from premature convergence.
Over the last few years, evolutionary algorithms, such as the genetic algorithms (GA) [3], evolutionary programming (EP) [1], the simulated annealing (SA), tabu search (TS) [2], differential evolution (DE) [8], and particle swarm optimization (PSO) [6], have been used to solve ED problems.Evolutionary algorithms are appropriate choices for solving DED problems as they prove to be very effective in solving without any restrictions in the shape of the cost curves and also because of their global search ability as well as their robust and effective constraint handling capacity.Perhaps, they do not always guarantee discovering the globally optimal solution in finite time.Therefore, hybrid methods combining two or more optimization methods were introduced [13]: to name a few, hybrid fuzzy PSO and Nelder-Mead (FAPSO-NM), hybrid DE and sequential quadratic programming (DE-SQP), hybrid GA (HGA), evolutionary strategy optimization (ESO), self-organizing hierarchical PSO (SOHPSO), variable scaling hybrid differential evolution (VSHDE), hybrid PSO, and sequential quadratic programming (PSO-SQP) [2] and there is a large still-growing body of literature.To be concise, few are included in [14].More precisely, hybrid methods combining probabilistic methods and deterministic methods are found to be very effective in solving complex optimization problems [2].
Recently, a new optimization technique known as teaching learning based optimization (TLBO) has been developed by Rao et al. [15,16].It is one of the recent evolutionary algorithms and is based on the natural phenomenon of teaching and learning process.It has already proved its superiority over other existing optimization techniques such as GA, ABC, PSO, harmony search (HS), DE, and hybrid-PSO.This research also proposes a hybrid method combining the teaching-learning based optimization (TLBO) and SQP [17].Here, TLBO will be performed as a base level search procedure, which makes a decision to direct the search towards the optimal region.Later the exact method (SQP) will be used to fine-tune that region to get the final solution [4,5].
A practical thermal power system integrated with the wind power generation system is used to demonstrate the application of the proposed hybrid methodology.The wind power system is designed and commissioned by the Suzlon Energy, India, with an installed capacity of a 150 MW (1.5 × 100 mills) wind farm.This wind farm is integrated with the Load Dispatch Center at Erode (LDCE) where a thermal system comprising 7 thermal units with multiple fuel options and valve-point loading effects is controlled.Several experiments are performed to validate the effectiveness of the proposed approach.

Mathematical Formulation of the DEDP
The objective function of DEDP is to minimize the total production cost of a power system over a given dispatch period, while satisfying various constraints [2].Mathematically the objective function is given as Generally, the generator cost function is usually expressed as a quadratic polynomial as Generators with multivalve steam turbines produce ripples like effect on their input-output curves.This effect, known as valve point effect, makes the generator cost function discontinuous and nonconvex.For accurate modeling of the cost function, the valve point effect is considered by superimposing it with the basic cost function Many generating units are supplied with multiple fuel sources and the cost functions of these units are represented with a few or several piecewise quadratic functions.Such a cost function is called a hybrid cost function and each segment of the hybrid cost function gives some information about the fuel burned.The hybrid cost function is given as For more accurate dispatch results, the valve point effect and the multiple fuel options are integrated into the basic cost function.Thus, the basic quadratic cost function given in (2) with  generating units and   fuel options for each unit is given as ( The objective function as given in ( 1) is subject to following equality and inequality constraints.
The power output from all the generating units must satisfy the total demand and the transmission losses of the system.The equality constraint is given as The transmission loss is expressed in a quadratic form as The real power output of each generating unit is limited by the maximum and minimum power limit of the units.It is given as The operating range of the generating units is restricted by their ramp rate limits.This is given as Thus, ( 8) is modified as 2.1.Accounting the Wind Power for the DEDP.Similarly, the nonlinear and fluctuation nature of wind proves to be great challenge for reliability and accuracy of power system that incorporates wind power.With accurate wind speed data, power system operator can predict the required power output.This helps in system planning, scheduling, and storage capacity optimization.To obtain proper and efficient wind power utilization wind speed prediction plays an important factor in forecasting.In the literature, several models have been used for wind speed prediction.The models are mainly physical and statistical models [18,19].Artificial neural network derives its computing power through its massively parallel distributed structure and its ability to learn and generalization for inputs.Neural network models have better performance than other models [12].There are many approaches used based on neural networks such as multilayer perceptron (MLP), ADALINE, back propagation network (BPN), radial basis function network (RBFN), and recurrent neural network (RNN).This paper proposes the radial basis function network for the forecasting of wind speed subsequently with the wind power generated.The wind power is estimated from the forecasted wind speed using the following expressions; as well known the wind power can be harvested only at a particular wind speed; thus, the wind power "" is given by where   is the wind turbine rated power,  is the actual wind speed,   is the wind turbine rated wind speed,  in is the wind turbine cutin speed, and  out is the wind turbine cutout speed.Thus, ( 6) is rewritten as where  ℎ is the wind power generated at time ℎ, and the total power generated from the entire wind farm containing "" wind mills is the summation of the wind power generated by the individual wind turbine.Thus, the DEDP will then be solved for economically dispatching the remaining demanded power using the multiple fuel power plant.

Evaluation Function.
We must define the evaluation function for evaluating the fitness of each candidate in the solution space.It is the sum of the generation cost function  obj and power balance constraint  bnc as in (1) and (3).The evaluation function is as follows: where max is the total fuel cost obtained using  ℎ =   max and   min is the total fuel cost obtained using  ℎ =   min : In order to limit the evaluation value of each candidate of the population within a feasible range, before estimating the fitness value of a candidate, the generated power output must satisfy the constraint given in (8).This evaluation function will be used to find the optimum design variable value for the DEDP.

Wind Power Forecasting Using the Radial Basis Function Network
For the reason that the speed of wind is unpredictable, accurate wind speed forecasting still remains a challenging task.Accurate wind speed prediction aids grid operator to dispatch economically the wind power generated to satisfy the power demand.Wind speed predictions can mainly be done using two models called physical and statistical models.The physical model considers the physical reasoning to get the best results.The statistical model considers online measurements of data.On the other hand, statistical models are more efficient than physical models where real time data is used [20,21].Artificial neural network (ANN) is one such statistical model, which derives its computing power through its massively parallel distributed structure and its ability to learn and generalization for inputs [22][23][24].The choice of ANN is also imperative in accurate wind speed/power forecasting.[12] network uses Gaussian function and has a long history in the applications of recognition and approximating function.RBF networks comprise three layers: the input layer, the hidden layer, and the output layer, where the hidden layer functions as layer of RBF units.The output layer is generally a linear function.The interconnection between input and hidden layer form hypothetical connection and between the hidden and output layers form weighted connections.Each hidden layer unit represents a single radial basis function, with associated center position and width.Each neuron on the hidden layer employs a radial basis function as a nonlinear transfer function to operate on the input data.The most often used RBF is a Gaussian function that is characterized by a center and width.RBF functions by measuring the Euclidean distance between input vector and the radial basis function center.

Radial Basis Function Network. Radial basis function (RBF)
The Gaussian RBF may be tuned by adjusting spread.It is less susceptible to problem with non-stationary input because of the behavior of RBF hidden units.The Gaussian function curve has a peak at zero distance and it decreases as the distance from the centre increases.
The Gaussian function is generally defined as where  in is the net input.
The advantages of RBFN are being more compact and having less training time, while eliminating local minima phenomena.The selection of the centers for the Gaussian function is important for nonlinear approximation.The weights between the hidden and output layer is then updated using the gradient descent rule.
The RBFN design includes structural and parameter design.The structural design involves finding number of neuron.The parameter design involves spread and weight of output mode.The RBFN is much more effective tool which gains the advantage of both generalizing and refining local features over BPN: where  is input vector, "" is number of neurons,   is th centre node in the hidden layer, and   is the weight between hidden and output layer.

The RBF-ANN Model for Forecasting the Wind Power.
This paper proposes a new RBF network architecture for forecasting the wind power by modifying the one designed in [12].Since the method is established as a reliable predictor of wind power, the training data used in this research are different and claimed as original contribution to ensure a more accurate prediction of wind power to assist the DEDP.
The wind speed depends on temperature, pressure difference in the two different parts of the ground, wind direction, and so on.In [12], the authors used three inputs to estimate the wind speed and they are the temperature, wind direction, and the past history of wind speed.Based on the information gathered from the Suzlon Energy Ltd., India, we understand that for a particular wind speed, the temperature and wind vane directions need not be the same at all instances.Perhaps the wind speed is influenced by various other factors.Thus, in this research, unlike [12], three different data of temperature and wind vane directions for the range of wind speed are used to train the RBFN.This improves the accuracy in predicting the wind speed and thereby the wind power with a chance of increasing the reliability of the installed wind power generation.
3.2.1.Architecture.For implementing the architecture, selection of inputs, hidden neuron, and outputs are required.
These selections depend on nature of the problem.The data required for inputs are wind speed, wind direction, and temperature.In spite of these three different data of temperature and wind vane directions, the same range of wind speed is used to train the RBFN.That is, for a particular wind speed there are three different data taken.This is because various other factors also influence the wind speed.Hence, training the RBFN with additional data for the same wind speed increases the prediction of wind speed.This improves the accuracy in predicting the wind speed and thereby the wind power with a chance of increasing the reliability of the installed wind power generation.Neural network with one hidden layer with a sufficient number of hidden neuron is capable of approximating any continuous function.The basic architecture is shown in Figure 1.
The sample of the collected data as shown in Figures 2  and 3 are used as input to the RBFN model.A large number of input parameters can be used for the analysis of the NN models.The input parameters are three different temperatures,  The selection of hidden neurons is very important.Based on trial and error method, 7 hidden neurons are selected with single hidden layer.The constructed architecture includes single hidden layer with 7 hidden neurons, 3 inputs, and 1 output for implementing the model and initializing the weights and epochs.The performance can change by increasing the epochs.The weights are calculated by using gradient descent rule.The error is the difference between actual and target value.The errors are minimized by adjusting the weights.I n s t a l l e dw i n dp o w e r M W 1 0 0× 1.5 MW (150 MW)

Methodology
Step 1 (data collection).The real time data is collected from Suzlon Energy Ltd., India, wind farm.The inputs are wind speed, wind direction, and temperature.The predicted wind speed is an output of the proposed model.The number of samples taken to develop a proposed model is 500, 1000, 2000, 5000, and 10000.The input parameters are as shown in Table 1.The wind energy depends on wind speed, wind direction, and temperature.
Step 2 (data normalization).The data normalization is carried out to improve accuracy of subsequent numeric computation and to obtain better output of model.The min max technique is used for normalization of input data.The advantage is preserving exactly all relationships in the data and it does not introduce bias.The normalization of data is obtained by the following.Normalized input is as follows: Step 3 (designing the neural network).Setup parameter includes learning rate, epoch, and dimensions.The training can be learned from the past data after normalization.The dimensions like number of input, hidden, and output neurons are to be initialized.The three input parameters are temperature, wind direction, and wind speed.The number of hidden layers is one.The number of hidden neurons is selected based on trial and error method.The input signal is processed and net input of model is computed.The net input is the weighted sum of inputs.The activation function is applied over the net input to calculate the output of neural network model.The sample inputs that are applied into the proposed model are as shown in Table 2.The sample inputs used for design of neural network are shown in Table 3.
Step 4 (training of network).For the purpose of developing models, the training, testing, and developing model at end stage for the past years in wind farms are needed.The data required for input are wind speed, wind direction, and temperature.Here, 70% of the total data is used for training the network.
Step 5 (testing of network).Evaluate the performance of network by testing models.Here, all available data is used for testing the network.Finally the predicted wind speed is the output of the NN architectures.
Thus, the RBF network is used to predict the wind speed from the past data and thereby the wind power.Once the wind speed is predicted, the wind power is estimated using the expression given in (11).Once the available wind power is estimated for the given horizon, the DEDP is solved to economically dispatch remaining power demand with the multiple fuel plant.This solves the DEDP using the teachinglearning based optimization refined using the sequential quadratic programming method which is discussed in the following section.

Teaching-Learning Based Optimization
Teaching-learning based optimization (TLBO) is an optimization technique developed by Rao et al. [15,16] based on teaching-learning process in a class among the teacher and the students.Like other nature-inspired algorithms, TLBO is also a population based technique with a predefined population size that uses the population of solutions to arrive at the optimal solution.In this method, population is also called students in a class and design variables are the subjects taken up by the students.Each candidate solution comprises design variables responsible for the knowledge scale of a student and the objective function value symbolizes the knowledge of a particular student.The solution having best fitness in the population (among all students) is considered the teacher.More specifically, an individual student (  ) within the population represents a single possible solution to a particular optimization problem.  is a real-valued vector with  elements, where  is the dimension of the problem and is used to represent the number of subjects that an individual, either student or teacher, enrolls to learn/teach in the TLBO context.The algorithm then tries to improve certain individuals by changing these individuals during the teacher and learner phases, where an individual is only replaced if his/her new solution is better than his/her previous one.The algorithm will repeat itself until it reaches the maximum number of generations.
During the teacher phases, the teaching role is assigned to the best individual ( teacher ).The algorithm attempts to improve other individuals (  ) by moving their position towards the position of the  teacher by referring the current mean value of the individuals ( mean ).This is constructed using the mean values for each parameter within the problem space (dimension) and represents the qualities of all students from the current generation.Equation (1) simulates how student improvement may be influenced by the difference between the teacher's knowledge and the qualities of all students.For stochastic purposes, randomly generated parameters are applied within the equation:  ranges between 0 and 1 and   is a teaching factor which can be either 1 or 2, thus emphasizing the importance of student quality: During the learner phase, student (  ) tries to improve his/her knowledge by peer learning from an arbitrary student   , where  is unequal to .In the case that   is better than   ,   moves towards   (20).Otherwise, it is moved away from   (21).If student  new performs better by following (20) or (21), he/she will be accepted into the population.The algorithm will continue its iterations until reaching the maximum number of generations: Additionally infeasible individuals must be appropriately handled to determine whether one individual is better than the other, when applied to constrained optimization problems.For comparing two individuals, the TLBO algorithm,  according to [15], utilizes Deb's constrained handling method [25]: (i) if both individuals are feasible, the fitter individual (with the better value of fitness function) is preferred; (ii) if one individual is feasible and the other one infeasible, the feasible individual is preferred; (iii) if both individuals are infeasible, the individual having the smaller number of violation (this value is obtained by summing all the normalized constraint violations) is preferred.

Sequential Quadratic Programming [17]
Sequential quadratic programming (SQP) [17] method belongs to the most powerful nonlinear programming algorithms for constrained optimization.Its excellent numerical performance was tested and compared with other nonlinear programming methods and since many years it belongs to the most frequently used algorithms to solve practical constrained optimization problems [17].The method resembles closely Newton's method for constrained optimization just as is done for unconstrained optimization.At each iteration an approximation is made of the Hessian of the Lagrangian function using a Broyden-Fletcher-Goldfarb-Shanno quasi-Newton updating method.This is then used to generate a quadratic programming subproblem whose solution is used to form a search direction for a line search procedure.The SQP subroutine for the EDP is adopted as it is in [17].

The Pseudo Code of the Proposed Refined TLBO Algorithm
The following steps enumerate the step-by-step procedure of the teaching-learning based optimization algorithm refined using the sequential quadratic programming method.(1) Initialize the number of students (population), range of design variables, iteration count, and termination criterion.
(2) Randomly generate the students using the design variables.
(3) Evaluate the fitness function using the generated (new) students.
//teacher phase// (4) Calculate the mean of each design variable in the problem.
(5) Identify the best solution as teacher amongst the students based on their fitness value.Use SQP method to fine-tune the teacher.
(6) Modify all other students with reference to the mean of the teacher identified in step 4.
//learner phase// (7) Evaluate the fitness function using the modified students in step 6.
(8) Randomly select any two students and compare their fitness.Modify the student whose fitness value is better than the other and use again the SQP method to fine-tune the modified student.Reject the unfit student.
(9) Replace the student fitness and its corresponding design variable.
(10) Repeat (test equal to the number of students) step 8, until all the students participate in the test, ensuring that no two students (pair) repeat the test.(11) Ensure that the final modified students' strength equals the original strength, ensuring there is no duplication of the candidates.
(12) Check for termination criterion and repeat from step 4.
The above procedure is used to solve the DEDP once the wind power is forecasted by the RBF network.Here the SQP method will be used to fine-tune the improving (better fitness) solution.This will ensure that the better solution region will not be overrun and will also aid in converging faster towards the possible best solution of the DEDP.

Numerical Experiments
This paper establishes the effectiveness of the proposed solution methodology for DEDP with integrated wind power by solving the following test systems as three different cases: (i) predicting the wind speed using RBFN using the data (wind speed, wind direction, and temperature, as shown in Figures 2 and 3) for a particular day (07 March, 2013; Thursday) provided by Suzlon Energy Ltd., for validating the prediction of wind speed; (ii) solving a standard 10-unit DEDP using the proposed TLBO-SQP to validate the superiority of the method over other existing methods in the literature (this is because the literature review shows that there is not a single paper archived (as per web of knowledge) for DEDP with multiple fuel option considered); (iii) solving a practical 7-unit multifuel DEDP with predicted wind power to share the total load demand on the DEDP (this is to validate the combined performance of the RBFN and the hybrid TLBO-SQP methodologies).
The proposed technique has been implemented in Matlab on a dual-core PC.The performance of the algorithm has been evaluated through simulation.Simulation studies have been carried out on three example test cases, over a scheduling time horizon of 24 hours.

Case 1: Wind Power Predictions.
The data for the wind power prediction is provided by the Suzlon Energy Ltd., India, as they play key role in predicting the wind power for the Load Dispatch Centre, Erode, (LDCE) which is controlling the Neyveli Thermal Power Station (NTPS).The wind farm consists of 100 wind turbines of 1.5 MW each.The LDCE forecasts wind power using Naive method as a generalized method for further scheduling the thermal units to fulfill the balance power demand requirements.Tables 1 and 2 present the training parameters for the back propagation network, RBF network, and improved RBF network.Additionally Table 3 shows the typical sample inputs and wind speed for a 24-hour horizon which is used by the BPN and RBFN [12], for predicting the wind speed.Similarly Figures 2 and 3 show the samples of the training data for the improved RBF network proposed in this research, which is claimed as the new proposal to predict the wind speed as far as this research is concerned.Figure 4 shows the wind speed predicted for 100 trials by the proposed RBFN method for a 24-hour horizon.This shows that the proposed RBFN predicts almost the reasonable and acceptable wind speed (thereby wind power) throughout all the trials.
In order to find the most suited NN architecture, mean square error (MSE) criteria are used.MSE as one of the error indices is used for the evaluation of performance of the networks chosen in this research and is given by where   is predicted output,    is actual output, and  is number of samples.The MSE is used as the criteria for measuring the forecasting performance.
The best suit NN architecture is preferred based on the accuracy of the prediction of wind power.The model accuracy is evaluated by comparing simulation results with actual/measured wind speed at the wind farm by the Suzlon.The experiments were carried out to test the performance of the improved RBF in comparison with that of BPN and RBFN [12].Based on the MSE summarized in Table 4, the improved RBFN proposed in this research proves predicting reasonably closer wind speed prediction than other methods.It is demonstrated that proposed improved RBFN is more accurate than other networks and found suitable for wind speed prediction in wind farms.
Figure 5 shows the comparison plot of the wind speed predicted by various methods for a 24-hour scheduling horizon, where the proposed improved RBFN is very much closely following the actual wind speed.Finally, Table 5 shows  the wind power estimated using (11), from the wind speed predicted using various methods.
Based on the above experiments and results, the proposed IRBFN is proved to be a most suitable alternative for predicting the wind speed and thereby wind power.Thus, this network will be adopted for predicting the wind power for the DEDP for a practical power system.Next the DEDP solution methodology proposed in this research will be validated with a standard 10-unit test system.

Case 2: Validating the Hybrid TLBO-SQP Method for
Solving DEDP.To validate the hybrid TLBO-SQP method for solving the DEDP, this research adopts a 10-unit test system proven as a complex test bench for several solution methods.Also the results are compared with the very recent solutions reported in the literature [14].Here the TLBO uses a population size of 100 and total generation of 10000.A total number of 30 trials were performed on the test system using the proposed hybrid TLBO-SQP method to arrive at the final conclusions about the superiority of the proposed method for solving the DEDP with multiple fuel options.The results are compared with the chaotic self-adaptive differential harmony search algorithm (CSADHS) proposed in [14].This research did not implement the CSADHS method but only compared the final results.
The system data and load demand considered for 24 h time interval is obtained from [2].In this test case, the DEDP only considers valve-point effect in the fuel cost function of generators.The optimal generation cost obtained by the proposed hybrid TLBO-SQP method is $1018679.2135as against the $1018681.872produced by the CSADHS algorithm [14] for a 24 h scheduling horizon and the corresponding production cost for each hour in the horizon is given in Table 6.Table 7 shows the production cost distribution as best, worst, and average production cost obtained for 30 different trial runs.The convergence characteristic of the hybrid TLBO-SQP method is shown in Figure 6.From Figure 6, it is demonstrated that the convergence characteristic of the proposed hybrid TLBO-SQP method is reliable and better than TLBO alone.The standard deviation of generation cost obtained by the proposed TLBO-SQP method is at an average of 1.5 for 30 different trial runs, which is less compared to CSADHS algorithm, thus proving that the proposed method is consistent in reaching the best optimal solution and best suitable for solving the DEDP with valve-point effects.

Case 3: A practical 7-Unit Multifuel DEDP with Integrated
Wind.This is a practical system where 7 thermal units are supplying the power demand and controlled by LDCE.During 2010, wind power generation is penetrated into the system with an installed capacity of 150 MW.As wind power is unpredictable, an improved RBFN is proposed to forecast the wind power and it is very well demonstrated through Case 1 that it is well suited for this system.The data for this practical system is given in Table 8.Similarly power demand for the total scheduling horizon is given in Table 9.The wind power for this entire scheduling horizon is already forecasted and will be used from Table 5 rounded to integer as done in the LDCE.The remaining power demand (difference in Tables 9  and 5) will be scheduled using the proposed hybrid TLBO-SQP method.
The optimal generation cost obtained by the proposed hybrid TLBO-SQP algorithm for 24 h time duration is $9538.1851as against the cost of 9736.1471obtained using the TLBO method alone.The corresponding generation schedule by the hybrid TLBO-SQP method is shown in Table 10.As seen in Table 11, the TLBO-SQP method can obtain the better total fuel cost and best total fuel cost compared to TLBO method, thus resulting in the higher quality solution.Moreover, in all the 30 different trial runs, the TLBO-SQP method produced almost less total fuel costs, thus conforming a better quality solution and convergence characteristic.The reliability of the proposed method in producing quality solutions (total fuel cost less than the other methods and convergence characteristics) is above 90% (as per 30 different trial runs).The average computation time of the proposed TLBO-SQP algorithm for 24 h time duration is 2.53 min as against the 4.21 min taken by the TLBO method.The convergence characteristic is shown in Figure 7. From Figure 7, it is revealed that the convergence characteristic of the proposed hybrid TLBO-SQP algorithm is steady and  fast.The standard deviation of generation cost obtained by proposed hybrid TLBO-SQP algorithm and TLBO are 10.26 and 20.14, respectively.This proves that the proposed algorithm is consistent in reaching the best optimal solution for this practical 7-unit multiple fuel generating station.

Discussion
Due to fluctuation and nonlinearity of wind speed, accurate wind speed prediction plays a major role in operational control of wind farms.Thus, the radial basis function network, a class in artificial neural network, is used to develop architecture for predicting wind speed in wind farms.The evaluation of neural network architecture is carried out for various real time data from wind farm.The performance analysis of modeling frame work was useful for accurate wind speed prediction in wind farm operators and thereby improves the reliability of wind power.Thus, on a particular day the wind speed is predicted by training the network with different samples (the temperature, wind direction, and the past history of wind speed) but data taken for similar wind speeds at various instances.Thus, the predicted wind speed is almost similar to the actual wind speed.
After the wind power is estimated (from the predicted wind speed), the difference in power is dispatched using practical 7-unit multiple fuel DEDP with valve-point effects.Based on the simulation results in Case 2, the hybrid TLBO-SQP method is found suitable producing quality dispatch schedule.Hence, the hybrid TLBO-SQP is used to further dispatch 7-unit multiple fuel DEDP meeting the balance power demand.
The above 7-unit multiple fuel DEDP is considered only with valve-point effects in problem formulation.Several constraints like prohibited operating zones, line flows, voltage limits, and transmission losses will also be included in the problem formulation for better accuracy of the dispatch solution and this will be done once practical data is received from the power station.

Conclusion
Integration of wind power into thermal generation system and its impact on the DEDP are explored in this paper.The DEDP formulation considering valve-point effects and multiple fuel options has been solved using a new hybrid TLBO-SQP algorithm.The DEDP formulation also incorporated wind power to demonstrate and assess the economic benefits of integrating wind power into power system grids.The wind speed prediction method has been analyzed using a radial basis function network and trained using real time wind data with different samples taken for similar wind speeds.Simulation studies revealed that the RBFN is found to be a more reasonable predictor than the existing methods used in load dispatch center.Based on the predicted wind power, the hybrid TLBO-SQP algorithm is used to solve the DEDP to source the remaining power demand with the available thermal generation.The performance of the TLBO algorithm is studied on a standard bench mark DEDP with 10 units, in terms of the convergence rate and accuracy of the solution produced.SQP is used on the run to fine-tune the better solutions produced by the TLBO algorithm.The effectiveness of the TLBO-SQP is finally illustrated on a practical 7-unit multiple fuel DEDP with valve-point effects with practical data.The comprehensive numerical results reveal that the proposed wind speed predictor using RBFN and the hybrid TLBO-SQP method to solve DEDP are most suitable for a more economical power dispatch for practical systems.

Nomenclature
The following nomenclature will be used throughout this paper: Wind turbine rated power : Actual wind speed   : Wind turbine rated wind speed  in : Wind turbine cutin speed  out : Wind turbine cutout speed.

Figure 2 :
Figure 2: Samples of wind vane direction and wind speeds.

Figure 3 :
Figure 3: Samples of temperature and wind speeds.

Figure 4 :
Figure 4: Wind speed predicted for 100 trials by the proposed IRBFN method for a 24-hour horizon (legend WP (hr) stands for wind power at that particular hour).

Figure 5 :
Figure 5: Comparison plot for wind speed predicted by various methods for a 24-hour horizon.

Figure 6 :
Figure 6: Convergence plot for both the methods for Case 2.

Figure 7 :
Figure 7: Convergence plot for both the methods for Case 3.

Table 1 :
Input parameters of the proposed models.

Table 2 :
Parameter selection for designing the various networks.  ,  min ,  max are the actual input data and minimum and maximum input data.  min ,   max are the minimum and maximum target values.

Table 3 :
Typical sample inputs and wind speed for a 24-hour horizon (data provided by Suzlon Energy Ltd).

Table 4 :
Comparison of MSE for various approaches.

Table 5 :
Wind power (for a particular trial) estimated from the wind farm for the predicted wind speed.

Table 6 :
Production cost obtained using the proposed hybrid TLBO-SQP method for each hour for Case 2.

Table 7 :
Summary of best, worst, and mean production cost produced by the various methods for Case 2.

Table 8 :
Data for the 7-unit practical system with multiple fuel cost coefficient.

Table 9 :
Power demand for Case 3.

Table 11 :
Summary of best, worst, and mean production cost produced by the various methods for Case 2.
ℎ ( ℎ ): Incremental fuel cost function ($/hr)  ℎ : R e a lp o w e ro u t p u to ft h eth unit at the ℎth interval (MW) : Number of generating units   : fuel options for each unit : Number of intervals in the given time horizon  ℎ : Power demand at the ℎth interval (MW)  Loss,ℎ : P o w e rL o s sa tt h eℎth interval (MW)   : Transmission loss coefficients  ℎ min / ℎ max : Minimum/maximum limit of the real power of the th unit at the ℎth interval (MW)   /  : Up/down ramp rate limits of the th unit (MW)  0 ℎ : P r e v i o u sp o w e rg e n e r a t e db yt h eth unit at the ℎth interval (MW)   : : Total production cost ($)