



Study with the several resources on Docsity
Earn points by helping other students or get them with a premium plan
Prepare for your exams
Study with the several resources on Docsity
Earn points to download
Earn points by helping other students or get them with a premium plan
research on wind turbines and the flow of power
Typology: Study Guides, Projects, Research
1 / 6
This page cannot be seen from the preview
Don't miss anything!




Abstract — Due to uncertainty of wind velocity, wind power generators don’t have deterministic output power. Utilizing wind power generation and thermal power plants together create new concerns for operation engineers of power systems. In this paper, a model is presented to implement the uncertainty of load and generated wind power which can be utilized in power system operation planning. Stochastic behavior of parameters is simulated by generating scenarios that can be solved by deterministic method. A mixed-integer linear programming method is used for solving deterministic generation scheduling problem. The proposed approach is applied to a 12-unit test system including 10 thermal units and 2 wind farms. The results show affectivity of piecewise linear model in unit commitment problems. Also using linear programming causes a considerable reduction in calculation times and guarantees convergence to the global optimum. Neglecting the uncertainty of wind velocity causes higher cost assessment of generation scheduling.
Keywords — Load uncertainty, linear programming, scenario generation, unit commitment, wind farm.
I. INTRODUCTION UE to increasing cost of fuels, there is a profound trend in using wind power. As the amount of wind power generation increases, concerns due to its effects on power system operation and expenses increases. As the wind is an intermittent power supply, developing forecasting tools is so useful. Though the forecasting tools become more accurate day by day, in the absence of a perfect estimation, decision about unit commitment should be done with uncertainties. Researchers have considered uncertainty of imprecise parameters using one of these two methods: fuzzy systems approach and probabilistic (stochastic) approach. Stochastic programming considers stochastic variables and uncertainty in conventional linear and nonlinear programming. Discussed randomness and uncertainty are generally represented using a probability density function (PDF). Unit commitment problem is an important optimization process to determine on/off scheduling of units in power system operation. The purpose is to minimize the total cost of system operation in order to satisfy system load demand and other constraints. Unit commitment problem is generally a nonlinear, large scale, and mixed integer combinatorial problem. This problem has been a dynamic research topic in the recent decades because of potential savings in operation
M. Esmaeeli Shahrakht is with Department of Electrical Engineering, Iran University of science and technology, Tehran, Iran. (Corresponding author to provide phone: +98-915-5623937; e-mail: [email protected]). A. Kazemi is with Department of Electrical Engineering, Iran University of science and technology, Tehran, Iran. (e-mail: [email protected]).
costs. As a result numerous methods including heuristics methods [1], [2], dynamic programming [3], genetic algorithm [4], Neural Networks, Fuzzy Systems [5], mixed integer linear programming [6], and Lagrangian relaxation [7]-[10] have been suggested. Among these methods Lagrangian relaxation is the most widely used approach because of its capability in solving large scaled problems. The most important defect of this approach is that it is required heuristics methods to achieve feasible solutions which might be suboptimal, due to nonconvexities of unit commitment problem. In contrast, mixed integer linear programming guarantees convergence to an optimal solution in a few steps. Moreover, it provides precise and flexible modeling framework [11]. One of the first efforts for solving unit commitment problem using stochastic optimization approach is presented in [12]. In this work the possible future result for the demand are shown with a scenario tree. In this method, it is assumed that only a limited number of scenarios are possible because using all the possible scenarios is computationally impossible. Authors have tested their method using 100 scenarios and a 48-hour time horizon. A stochastic model for long-term security constrained unit commitment is used in [13], which outage of generating units, transmission lines and load prediction errors are considered. In this paper, stochastic unit commitment of thermal power plants and wind farms is implemented with constraints for wind power generation. Wind speed uncertainties during each period of operation time is implemented through several scenario based on the wind speed PDF. To employ the mixed integer linear programming method and achieving global optimum, fuel cost function of thermal units is approximated by a set of piecewise block such that it is not distinguishable from the nonlinear model. The main contributions of this work are as follows:
International Journal of Electrical, Computer, Energetic, Electronic and Communication Engineering Vol:7, No:2, 2013
International Science Index, Electrical and Computer Engineering Vol:7, No:2, 2013 waset.org/Publication/
Load uncertainty and inaccuracy of forecasting for wind speed at all time intervals, faces wind-thermal unit commitment problem to the uncertainty. To model this uncertainty a set of scenarios are generated using Monte Carlo simulation method based on load and wind speed probability density functions. Then cost of operation for these scenarios is evaluated using linear model of mid-term unit commitment. Expected values of parameters are determined according to their values in each scenario and probability of the scenario. A. Load Uncertainty Modeling Due to imprecise load forecasting, there is an uncertainty in expected future load. Typically in literature normal distribution is used for modeling load forecasting error. In this paper, Pd(t) is the forecasted weekly peak load in time interval t and its error assumes that it has a normal distribution of N ( μ σ, ) where μ ( ) t and σ ( ) t are its mean and standard deviation respectively. These quantities in each time interval of t can be determined according to statistic studies. B. Wind Power Uncertainty Modeling Wind speed error can be represented by a suitable PDF. The PDF related to the difference between predicted and measured wind velocity, can be detected with a Gaussian distribution. In [14] the authors show that PDF of wind speed prediction error, in the speeds range suitable for wind energy applications, is a Gaussian distribution. We assume that the uncertainty of wind speed forecasting in all wind farms is a normal distribution with two known parameters, μw (^) i (t) and σ (^) wi (t), respectively for mean and standard deviation of i th wind farm. The generated wind power varies with wind speed at the site of wind farms. The output power of wind turbine could be determined from its power curve which is a plot of output power against wind speed [15]. Fig. 1 shows conventional curve of wind turbine. It is designed to start power generation at cut-in wind speed ( V (^) ci) and be shut down at the cut-out wind speed ( V (^) co) for safety. When the wind speed is higher than rated wind speed ( V (^) r) and lower than cut-out wind speed ( V (^) co), the rated power is generated. As shown in Fig. 1 there is a nonlinear relationship between generated output power and wind speed when wind speed is located between V (^) ci and V (^) r. Therefore the generated wind power by such turbine will be found by:
2
i ci i r i i ci i r r r i co i co
The calculation of constants A, B and C is discussed in [16]. In this paper typical model of wind turbine with cut-in, rated and cut-out wind speed respectively equal to 3, 11.5 and
25 m/s and rated output power equal to 2 MW is used. It is assumed that there are 40 turbines in each wind farm.
Fig. 1 Typical power curve of a wind turbine C. Scenario Generation Inverse transform method is used for scenario generation. It is assumed that PDFs of loads and wind speeds in each time interval is available. To determine a random value for these parameters, cumulative distribution functions (CDFs) is plotted based on their PDFs. Then, a random decimal number between zero and one is selected and corresponding value of the parameter is determined from the cumulative curve. There are 3 stochastic parameter (load and wind speed of 2 wind farms) and 12 time intervals. So each scenario has 36 components (3*12). Each component is a random variable generated based on above approach.
III. UNCERTAINTY MODELING IN UNIT COMMITMENT PROBLEM FORMULATION In this section a model for mid-term unit commitment scheduling is presented in which the uncertainty for main elements is modeled. The time horizon of scheduling for this problem is one season with a weekly interval. Due to time scheduling of this problem in its formulation, ramping and minimum up/down constraints are neglected. A. Problem Formulation The objective function of mid-term unit commitment problem is minimizing the cost of demand supply, while the constraints of required reserve, wind power generation and units capacity are satisfied. The objective function of scenario s is defined as total cost including production cost of thermal units and the operation and maintenance cost of thermal units and wind farms. It can be formulated as:
International Journal of Electrical, Computer, Energetic, Electronic and Communication Engineering Vol:7, No:2, 2013
International Science Index, Electrical and Computer Engineering Vol:7, No:2, 2013 waset.org/Publication/
scenarios.
IV. SIMULATIONS AND NUMERICAL The Proposed optimization algorithm is applied to a 12 test system including 10 thermal units and 2 wind farms. The fuel cost function of thermal units is fitted using a piecewise linear function with four segments. The time horizon of scheduling for this problem is one season with a weekly interval. The input data for thermal units are given in [17] and information of two wind farms are added to them. maximum output power of each wind farm is mean values and standard deviations of weekly peak load and estimated wind speed for each farm in a given in [15]. In this paper, first we compare the results of piecewise linear and nonlinear models in a deterministic unit commitment problem. It will also be shown the effectiveness of the proposed piecewise linear model for unit co problem which considerably reduces the calculation time and guarantees the problem convergence. Then the linear model of stochastic unit commitment will be solved and the results will be compared with the results of problem deterministic solution. The effect of installed wind power is another subject that is analyzed in this paper. A. Deterministic Solution of Unit Commitment In this section a deterministic unit commitment will be considered based on mean values of parameters in each period neglecting their uncertainties. The deterministic mid commitment problem for the test system based on nonlinear model is solved. The results of committing states of units are shown in Table I. Then the problem will be solved based on piecewise linear model and the results of committing states of units are shown in Table II. As seen in these tables, the results of linear programming obtained using piecewise linear approximation of fuel cost function are pretty close to those of nonlinear programming obtained using quadratic fuel cost function. Because suitable precision is considered in the linear approximation and the number of pieces are selected accurately. The cost of operation and calculation time for solving the problem in both cases is presented in Table confirm the effectiveness and usefulness of piecewise linear model. The noticeable reduction of calculation time in linear model is the very important advantage of this model while the results are similar in both models. Also, the effect of the amount of installed wind power on the scheduling cost of UC is studied. For this purpose, once installed wind power is increased to 50% (i.e. installed wind power is 2120 MW) and once it is decreased to 50% (i.e. installed wind power is 240 MW). The cost of UC in the absence of wind farms is also calculated. The results of these studies are given Fig. 3. Although increasing in installed wind power causes reduction in cost of UC, but it does not reduce
Determining desired values of output variables based on their values in different scenarios and the occurrence
The Proposed optimization algorithm is applied to a 12-unit test system including 10 thermal units and 2 wind farms. The fuel cost function of thermal units is fitted using a piecewise linear function with four segments. The time horizon of scheduling for this problem is one season with a weekly input data for thermal units are given in [17] and information of two wind farms are added to them. The power of each wind farm is 80MW. The mean values and standard deviations of weekly peak load and estimated wind speed for each farm in all time intervals are
In this paper, first we compare the results of piecewise linear and nonlinear models in a deterministic unit commitment problem. It will also be shown the effectiveness of the proposed piecewise linear model for unit commitment problem which considerably reduces the calculation time and guarantees the problem convergence. Then the linear model of stochastic unit commitment will be solved and the results will be compared with the results of problem deterministic The effect of installed wind power is another subject
Deterministic Solution of Unit Commitment In this section a deterministic unit commitment will be considered based on mean values of parameters in each period their uncertainties. The deterministic mid-term unit commitment problem for the test system based on nonlinear model is solved. The results of committing states of units are I. Then the problem will be solved based on and the results of committing states of these tables, the results of linear programming obtained using piecewise linear approximation of fuel cost function are pretty close to those of using quadratic fuel cost function. Because suitable precision is considered in the linear approximation and the number of pieces are selected
cost of operation and calculation time for solving the Table III. The results confirm the effectiveness and usefulness of piecewise linear model. The noticeable reduction of calculation time in linear model is the very important advantage of this model while the
ect of the amount of installed wind power on the scheduling cost of UC is studied. For this purpose, once installed wind power is increased to 50% (i.e. installed wind power is 2*120 MW) and once it is decreased to 50% (i.e. ). The cost of UC in the absence of wind farms is also calculated. The results of these
according to the constant rate.
TABLE COMMITTING STATES OF UNITS IN NONLINEAR (1=ON Units 1 2 3 4 5 Th1 1 1 1 1 1 Th2 1 1 1 1 1 Th3 1 1 0 1 1 Th4 1 1 1 1 1 Th5 1 1 1 1 1 Th6 1 1 0 0 1 Th7 0 0 0 0 0 Th8 0 0 0 0 0 Th9 0 0 0 0 0 Th10 0 0 0 0 0
TABLE COMMITTING STATES OF UNITS IN DETERMINISTIC Units 1 2 3 4 5 Th1 1 1 1 1 1 Th2 1 1 1 1 1 Th3 1 1 0 1 1 Th4 1 1 1 1 1 Th5 1 1 1 1 1 Th6 1 1 0 0 1 Th7 0 0 0 0 0 Th8 0 0 0 0 0 Th9 0 0 0 0 0 Th10 0 0 0 0 0
TABLE DETERMINISTIC Modeling method Nonlinear model Piecewise linear model
Fig. 3 The effect of installe B. Stochastic Solution of Unit Commitment The main idea of stochastic optimization is extension and simulation of possible choices. Scenario generation is a conventional way to model these probable conditions. Load and wind speed are two stochastic parameters in unit commitment problem. To simulate these uncertainties scenario generation based on inverse transform is used. Then
56
57
58
0 40 80
Cost of Operation (M$)
Installed Wind Power (MW)
according to the constant rate.
TABLE I ONLINEAR MODEL OF DETERMINISTIC UC ON , 0=OFF) Periods 6 7 8 9 10 11 12 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 1 1 1 0 0 0 0 0 1 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0
TABLE II NITS IN PIECEWISE LINEAR MODEL OF ETERMINISTIC UC (1=ON , 0=OFF) Periods 6 7 8 9 10 11 12 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 1 1 1 1 1 1 0 0 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 1 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0
TABLE III ETERMINISTIC UC RESULTS Cost Calculation time M$ 56.67 6.5 s M$ 56.66 0.5 s
Fig. 3 The effect of installed wind power on operation cost Stochastic Solution of Unit Commitment The main idea of stochastic optimization is extension and simulation of possible choices. Scenario generation is a conventional way to model these probable conditions. Load d are two stochastic parameters in unit commitment problem. To simulate these uncertainties scenario generation based on inverse transform is used. Then
120 160 200 240 280 Installed Wind Power (MW)
International Journal of Electrical, Computer, Energetic, Electronic and Communication Engineering Vol:7, No:2, 2013
International Science Index, Electrical and Computer Engineering Vol:7, No:2, 2013 waset.org/Publication/
deterministic UC problem is solved for each scenario. A combination of results, usually the weighted average of results, is selected to present stochastic solution. In this study, 180 scenarios are generated to simulate the uncertainty of the problem. All scenarios have the same probabilities. The piecewise linear model of fuel cost function is used in this work and linear programming is applied to solve UC problem. For precise analysis of stochastic unit commitment problem, 6 scenario groups each containing 30 scenarios have been considered. The indexes related to total cost of stochastic unit commitment are presented in Table IV. These indexes are including mean value, standard deviation and variation coefficient of total cost.
TABLE IV STOCHASTIC UC RESULTS WITH DIFFERENT SCENARIOS IN LINEAR MODEL Variation coefficient (%)
Number of Total cost (M$) scenarios (^) Standard Deviation
Mean 130 57.414 1.762 0. 230 56.962 1.699 0. 330 56.952 1.562 0. 430 56.827 1.612 0. 530 56.731 1.619 0. 630 56.688 1.638 0.
As seen, the total cost by increasing the number of scenarios has been decreased because the model has became more precise and it more appropriately covers the uncertainties of the problem. In each stage by increasing the number of scenarios, the mean of total operation cost and variation coefficient are decreased. The load supplying contribution of units in stochastic UC is shown in Fig. 4. As seen, the most of units are committed at last week when the system demand is at the highest level. Also, all of generated wind powers have been used in each period. To compare stochastic and deterministic approaches the load supplying contribution of units in deterministic UC is presented in Fig. 5.
Fig. 4 Mid-term stochastic UC results for load contribution
Fig. 5 Mid-term deterministic UC results for load contribution
The calculation time for solving linear model of stochastic unit commitment is 19 minutes. In order to scenarios reduction and hence calculation time reduction, a composed algorithm called backward reduction and forward selection has been used. Three different cases for scenarios reduction have been investigated by this method. The results of these cases are presented in Table V. As seen in Table V, more reduction in the number of scenarios leads to more inaccuracy. But the calculation time is decreased in accordance with percentage of reduction. Therefore, a trade-off between the accuracy and the computation time must be considered.
TABLE V RESULTS OF STOCHASTIC UC IN DIFFERENT PERCENTAGE OF REDUCTION Inaccuracy (%)
Operation Cost
Number of Scenarios
Cases
Base Case 180 56.688 - 85% reduction in base case 27 56.676 0. 90% reduction in base case 18 56.663 0. 95% reduction in base case 9 56.653 0.
In this paper a linear optimization method considering uncertainty in power systems was proposed to solve stochastic UC problem. Utilizing linear optimization method leads to achieving a global optimum without divergence problems seen in nonlinear problems. UC uncertainty is caused by load prediction error and stochastic nature of wind speed in wind farms. In this paper first, mid-term stochastic UC problem as a deterministic optimization problem is solved using linear and nonlinear programming. The results show affectivity of piecewise linear model in unit commitment problems. Also using linear programming causes a considerable reduction in calculation times. Then the uncertainties are simulated by utilizing a set of scenarios. Results show reduction in the number of scenarios leads to inaccuracy. But the calculation time is decreased in accordance with percentage of scenario reduction. In addition, using linear methods causes the average of the scenarios solution, unlike nonlinear methods, does not lead to unfeasible or suboptimal solutions and the global optimum could be accessed.
International Journal of Electrical, Computer, Energetic, Electronic and Communication Engineering Vol:7, No:2, 2013
International Science Index, Electrical and Computer Engineering Vol:7, No:2, 2013 waset.org/Publication/