Research Article  Open Access
Peng Wei, Yang Liu, "The Integration of WindSolarHydropower Generation in Enabling Economic Robust Dispatch", Mathematical Problems in Engineering, vol. 2019, Article ID 4634131, 12 pages, 2019. https://doi.org/10.1155/2019/4634131
The Integration of WindSolarHydropower Generation in Enabling Economic Robust Dispatch
Abstract
Currently the renewable energies including wind power and photovoltaic power have been increasingly deployed in power system to achieve contamination free and environmentalfriendly power production. However, due to the natural characteristics of wind and solar, both wind power and photovoltaic power contain uncertainty and randomness which may significantly impact the stability, security, and economic efficiency of the conventional power system mainly consisted by hydropower and thermal power. To deal with the issue, this paper presents a twostage robust model which is able to achieve the optimal dayahead dispatch strategy in the worstcase scenario of wind and photovoltaic outputs. Because of the strong interactions between the two stages, the original optimization has been decomposed into the dayahead dispatch master problem and the additional adjustment subproblem considering the uncertainty and randomness of the wind and the photovoltaic outputs. Also, the piecewise linearization technique is employed to convert the original problem into a MILP problem. Afterward, the dualization of the additional adjustment subproblem can be obtained by using linear programming strong duality theory. Additionally, the BigM method enables the linearization of the dual model. The interactediterations between the master problem and the subproblem are successfully implemented which can ultimately figure out the optimal dayahead dispatch strategy of the power system with conventional and renewable energies.
1. Introduction
To achieve the contamination free and environmentalfriendly power generation, clean energies such as wind power generation and photovoltaic power generation have been increasingly deployed in modern power systems [1]. It has been widely admitted that such renewable energies can effectively deal with the increase of the power consumption with less environmental impacts. However, it is also realized that the uncertainty and randomness of the wind and photovoltaic power generations are extremely difficult to be predicted. Therefore, the security and the economicalefficiency including voltage, frequency, and power loss of the conventional power system may be severely affected with high penetration of the wind and photovoltaic power generations. As one of the conventional power generations in power system, hydropower is regarded as a kind of energy source with highly flexible adjustment ability, which has remarkable potential to effectively handle the uncertainty and randomness of the wind and photovoltaic power generations. By coordinating the hydropower, accommodation capability of serving the wind power and photovoltaic power can be significantly improved [2]. Therefore this paper combines the wind power, photovoltaic power, hydropower, and thermal power for coscheduling. Additionally, it should be pointed out that the integration of multiple energy sources may also greatly improve the operating efficiency of power system [3, 4]. Research [3] presents a new mathematical model for the simulation of solid oxide fuel cell based on the combined heat and power system, in which the cogeneration units have high operating efficiency. Research [4] presents a novel polygeneration system, including cogenerative proton exchange membrane fuel cell, concentrated photovoltaicthermal collectors, and so on. The experimental result shows that the total electrical and thermal efficiencies remain high. It can be seen that the integration of multiple energy sources is economically viable. Quite a few researches [5–26] also focus on the integration of the multiple energies in the power system, and also they point out that the optimal dispatch strategy for the coordination of the energies becomes highly important.
To study the integration of the clean energies, researches [5–7] mainly aim at improving the system reliability and flexibility using the optimal coordination of wind power, photovoltaic power, and hydropower. However, they have not considered the uncertainty and randomness of the wind power and the photovoltaic power which may lead to errors in their researches. The deterministic programming and stochastic programming are the common approaches to handle the uncertainties of the renewable energy sources [8]. Researches [9, 10] add the spinning reserve into the deterministic dispatch model. Although their approach improves the system security, it simultaneously limits the economy and the flexibility of the system. Researches [11, 12] establish the uncertainty dispatch models based on the probability distributions, for example, Weibull and Beta distributions. However, it should be realized that these two hypothetical distributions cannot accurately describe the behaviors of wind and solar so that the deviations certainly occur between their dispatch model and the practical requirements. Researches [13, 14] mainly employ the multiscenario approach to convert the stochastic optimization problem into the deterministic problem. The multiscenario approach trades off the scenario reduction in terms of computing efficiency. In the scenario reduction operations, certain extreme scenarios may be lost resulting in the inaccuracy of the optimization [15]. Recently, robust optimization has been widely employed in the economical dispatch as it needs neither precisely modeled probability distribution nor largescale sampling of uncertain variables [16–18]. Research [16] presents a static robust optimization model to deal with the uncertainties. However, the static robust optimization tends to be overconservative. Thus, researches [17, 18] further present the twostage robust optimization model with the “minmaxmin” structure. Additionally, an adjustable robust parameter is introduced to control the conservatism of the scheme. Based on their studies, the original problem can be decomposed into the master problem and the subproblem based on the type of decision variables.
Due to the uncertainty in the renewable energy sources, the adjustment between the base generation and the actual generation should be considered. Researches [19–21] place emphasis on the acquisition of adjustment strategy for unbalanced power caused by the prediction error. Research [19] introduces the multiagent technology in the microgrid to formulate the energy scheduling strategy. Research [20] presents a method that considers variability and uncertainty cost of renewable energy sources between two consecutive scheduling intervals in the realtime. Research [21] further links up the dayahead scheduling and realtime scheduling by a more meticulous generation dispatch mode to formulate a comprehensive strategy. Researches [22–26] introduce the market mechanism to ensure the economic efficiency of power system considering the risk due to uncertainties. The overestimation and underestimation costs of the renewable energy sources are considered in the objective, which is similar to the upregulation and downregulation in the adjustment stage.
Motivated by the previous works, this paper presents a twostage robust model considering the generation of wind power, photovoltaic power and hydropower in enabling the dayahead economical dispatch strategy for power system. In our model the fluctuations of wind output and photovoltaic output caused by the randomness and uncertainty are evened out using hydropower and thermal power. Especially due to the strong adjusting ability of hydropower, the model deals with the fluctuations prioritizing the hydropower and then the thermal power because of its feeblish adjusting ability. Additionally, to improve the wind and photovoltaic power accommodation capability of the power system, this paper regards the wind power and photovoltaic power curtailment costs as penalty function in the presented model. This paper also linearizes the original twostage robust problem and then employs the strong duality theory to handle the maxmin subproblems. At last the ultimate optimal solution can be achieved using the model decomposition and iterations. The effectiveness of the presented approach has been evaluated using the IEEE30 bus system.
The rest of the paper is organized as: Section 2 presents the dayahead dispatch robust model considering the wind power, photovoltaic power, and hydropower; Section 3 presents the linearization of the model; Section 3 solves the preprocessed twostage model; Section 4 evaluates the approach and discusses the simulation results; Section 5 concludes the paper.
2. DayAhead Dispatch Robust Model Considering Wind Power, Photovoltaic Power, and Hydropower
2.1. The Details of the TwoStage Robust Model
Modern power systems contain not only thermal power and hydropower but also wind power and photovoltaic power. In this section, the optimal dispatch model is presented for a power system with thermal power and renewable energies including wind, photovoltaic, and hydropower. At the very beginning of the modeling, without considering the uncertainty and randomness of the wind and photovoltaic power, the dayahead dispatch strategy is an optimization problem based on the identification of the minimum operating cost, which is represented by where (·) is the conventional dayahead operating cost of the power system, which is mainly the cost of the thermal power; x is the dayahead dispatch scheme including the status and the outputs of the thermal units and the hydropower units; (·) and (·) represent the dayahead equality constraint and the inequality constraint representing the power balance constraint, branch power flow constraints, thermal power units constraints, hydropower unit constraint, wind and photovoltaic outputs constraints, and so on.
In (1), the outputs of the wind and photovoltaic generation are based on the prediction. However, precise prediction of the outputs is difficult so that the power imbalance in the system certainly occurs [27]. The imbalance can be compensated by using additional adjustments which can be regarded as a second stage problem named as additional adjustment stage. As a result, in our model, the outer layer is a max problem which is adopted to identify the worstcase scenario (the scenario under which the adjustment cost is the highest) with the x based uncertain parameters. The inner layer is a min problem which is adopted to figure out the minimum additional adjustment costs under the worstcase scenario. The overview of the model is represented by where (·) is the additional adjustment costs based on the dayahead dispatch scheme x; (·) is mainly consisted of units adjustment cost, wind power and photovoltaic power curtailment costs; u is the uncertain parameters; U is the feasibility region of u; y is the adjustment scheme which is able to handle the uncertainty; Ω(·) is the feasibility region based on certain known values (e.g., means that feasibility region of y is constrained by x and u that have been optimized); (·) and (·) represent the power balance constraint, branch power flow constraints, thermal power units constraints, hydropower unit constraint, wind and photovoltaic outputs constraints, and so on in the adjustment stage.
Combining the dayahead operation and the additional adjustment a model is represented by (3). The model can achieve the minimum cost dispatch scheme x under the worstcase scenario considering the uncertainties of wind and photovoltaic outputs.
In (3), x is a type of “hereandnow” variable which represents the optimal dayahead dispatch scheme under the worstcase scenario. In addition, it can hedge against any possible condition within the uncertain sets; y is a type of “waitandsee" variable which represents the optimal adjustment scheme under the worstcase scenario. And it can vary with different operation conditions.
To clarify the model represented by (3), Section 2.2 presents the uncertain parameter u; Section 2.3 presents the objective function considering (·) and (·); Section 2.4 presents the dayahead constraints (·), (·) and the adjustment constraints (·), (·).
2.2. Uncertain Parameter
The uncertain variable u is highly related to the outputs of the wind power and the photovoltaic power. Based on the uncertain set [28, 29], the outputs of the wind power and the photovoltaic power can be represented by the intervals shown by where is the predicted value of the wind power output at time t; is the predicted value of the photovoltaic power output at time t; and are the maximum fluctuation ranges of the wind power output and the photovoltaic power output, respectively.
Additionally, this paper also employs the adjustable robust parameters shown by where T is the simulation period; and are the adjustable robust parameter of the wind power output and the photovoltaic power output, respectively.
2.3. Objective Function
To achieve the optimal economical dayahead dispatch scheme x, costs in different stages are considered. In the first stage the cost is mainly the generation cost of thermal power units . The costs in the second stage mainly include the adjustment cost of thermal units , the adjustment cost of hydrounits , the wind power curtailment cost , and the photovoltaic power curtailment cost . The objective function is represented bywhere a, b, and c are the coefficient of the thermal power units; is the outputs of the thermal power units at time t; and are the penalty prices for the adjustment power of thermal units and hydrounits; and are the penalty prices of wind power curtailment and photovoltaic power curtailment; and are the output adjustment of thermal power and hydropower at time t; and are the wind power and photovoltaic power which are injected into the power grid at time t.
2.4. Constraints
2.4.1. DayAhead Dispatch Constraints
Dayahead dispatch constraints in the model represented by (3) contain the following:
(a) Thermal Power Constraintswhere is the 01 variable which denotes the unit commitment status of the thermal power at time t; and are the upper limit and the lower limit of the outputs of the thermal units; and are the upward and downward ramping capability of the thermal power units.
(b) Hydropower Constraintswhere is the water consumption of the hydrogenerators at time t; is the output of the hydrogenerators at time t; and are the converting coefficients of hydrogenerators; and are the upper and lower limit of water consumption in every subperiod; is the water curtailment of reservoir; U_{t} is the reservoir storage at time t; is the amount of water that flows into the reservoir at time t; and are the upper limit and the lower limit of the reservoir storage. and are the initial value and final value of the reservoir storage.
(c) Power Balance Constraintwhere is the power demand of the load in the system at time t.
(d) Branch Flow Constraintwhere l is the identifier number of the transmission lines; b is the identifier number of the nodes; N_{b} is the total number of the nodes; is the sensitivity factor of the line to the node; is the net active power which is injected into the node at time t; is the maximum transmission power of the line.
2.4.2. Additional Adjustment Constraints
Additional adjustment constraints in the model represented by (3) contain the following:
(a) Thermal Power Units Regulation Constraintswhere is the 01 variable which represents the adjustment status of the thermal power units at time t; and are the upper and lower adjustment limits of the thermal power units. Equation (19) guarantees that the adjusted outputs should still be constrained in the permissible range.
(b) Hydropower Units Regulation Constraintswhere is the binary variable which represents the adjustment status of hydrogenerators at time t; and are the upper and lower limit of the hydrogenerator adjustment power; is the water consumption after the adjustment at time t; is the water curtailment of reservoir after the adjustment at time t; is the reservoir storage after the adjustment at time t. The constraints (21) guarantee that the hydropower should still be constrained in the permissible range after the adjustment.
(c) Power Balance Constraint
Equation (22) indicates that the outputs of multiple powers after adjustment should equal to the power consumption of the load in the power system at time t.
(d) Branch Flow Constraintwhere is the net active power injected into the node after adjustment at time t.
(e) Wind Power Constraintwhere is the actual output of wind turbine generators at time t, which is an uncertainty variable.
(f) Photovoltaic Power Constraintwhere is the actual output of photovoltaic generators at time t, which is an uncertainty variable.
3. Linearization of Thermal Power Generation Cost
The dayahead economical dispatch is a mixed integer quadratic programming (MIQP) problem. Using robust optimization directly may be time consuming to achieve the optimal solution. Therefore, this paper employs piecewise linearization approximation [30] to convert the MIQP problem into a mixed integer linear programming (MILP) problem.
In our model, the nonlinear factor (the first equation in (6)) exists only in thermal power. Figure 1 shows the relationship of the power system operating cost and thermal power output according to the first equation in (6). To linearize the relationship, piecewiselinear approximation is employed. Firstly, we segment the curve into a number of N segments. And then in each segment, by employing the continuous variable P_{j,t} and the state variable S_{j,t}, linear approximation can been done. As a result, the linearized cost of the thermal units is presented by where is the linearized cost function of the thermal units; m_{j} is the slope of segmentation; n_{j} is the intercept of the segmentation; P_{j,t} is the output of the segmentation; S_{j,t} is the status of the segmentation at time t; P_{j} is the segmentation point.
Based on the linearization, the original twostage model can be represented using the matrix shown by where , , and are the coefficient matrices in the objective function; , , , and are the coefficient matrices in the constraints; and are the column vectors; is the vector of all the firststage decisions; is the vector of all the second stage decisions; is the vector of all the uncertainty variables.
4. Solution for the Model
4.1. Model Decomposition
In our model, the worst scenario and its corresponding additional adjustment scheme in the second stage are optimized based on the dayahead dispatch scheme in the firststage optimization. Moreover, the optimization in the second stage also impacts the optimized results of the first stage. The interactions between the two stages result in the employment of C&CG algorithm [18] to decompose the original model into the MP (Master Problem) and SP1 (Subproblem). Additionally, in MP an auxiliary variable is adopted to substitute for the objective of SP1.
4.2. The Equivalent Maximization Problem of SP1
This paper employs strong duality theory to convert the SP1 into an equivalent maximization problem. BigM approach is also adopted to handle the nonlinear terms in the dual problem. The converted SP2 is represented inwhere , , are the dual variables; ,u^{} are the upper and lower limit of uncertainty variable; is the predictive value of uncertainty variable; ^{+}, ^{} are the positive and negative values of ; , are the binary variables introduced by the BigM approach; the adjustable robust parameter Γ includes and .
4.3. Iterations between MP and SP2
Our model is finally converted into the MP and the SP2. The optimizations of OP (28) can be ultimately achieved based on the iterations using the following steps.
Step 1 (initialization). Initialize the worstcase scenario . Set the upper bound of the problem U = +∞ and the lower bound of the problem L=∞. Set the number of iterations k=1 and the convergence gap > 0.
Step 2 (master problem computation). According to the worst scenario (l=1, 2, …, k), solve the MP and achieve the optimal solution (x_{k}, ,y_{1},y_{2}, …,y_{k}), and then set .
Step 3 (subproblem computation). After given the firststage solutions x = x_{k}, the solution (, ) of the SP2 can be obtained. Set the worstcase scenario , .
Step 4 (convergence check). If UL≤λ, return the solution x=x_{k}. Otherwise add the new variable and the following constraints which are associated with the new worstcase scenario, then update k = k+1 and go to Step 2 iteratively.
5. Experimental Result
To evaluate the effectiveness of the presented twostage robust optimal model in enabling the economical dayahead dispatch considering the uncertainties, this paper implements the model using MATLAB. The solution of the model is by employing CPLEX. A modified IEEE30 bus system is employed as the power system model. The load fluctuations are shown in Figure 9. The parameters of the thermal power units are shown in Table 5. This paper only considers one case that the wind power, photovoltaic power and hydropower are connected at the same node (No. 30) in the IEEE30 bus system. The capacity of the hydropower is 18MW while the other parameters are listed in Table 6. The capacity of the wind power is 10MW and the capacity of the photovoltaic power is 10MW. Additionally, according to the research [29], the fluctuations and the uncertain intervals of the wind power and photovoltaic power are shown in Figures 10 and 11. The penalty price of the wind power and photovoltaic power curtailed is 10$/MWh.
5.1. Experimental Result of the TwoStage Robust Optimization
The values of and are set to 10 and 6. During the solution of the twostage robust model, the C&CC approach converges after 6 iterations, which is shown in Figure 2. The final optimized thermal power output and the hydropower output are shown in Figure 3. The worst scenarios of the wind power and the photovoltaic power are shown in Figure 4.
Figure 3 indicates that, compared to the traditional deterministic optimization, the output of the hydropower optimized by the robust optimization is lower in the period of 8h24h. Contrarily, the output of the thermal power becomes higher. The reason is that the robust optimization considers the worst scenario of both wind power and photovoltaic power. In the worst scenario, the output of both wind and photovoltaic power may be lower than the predicted values. As a result, thermal power and hydropower need to output more to compensate the imbalance. Additionally, more adjustment space of hydropower has been reserved at 8h24h to handle the uncertainties of the renewable energies. As the adjustment ability of the hydropower is flexible and also the cost of the hydropower is lower, this point significantly reduces the adjustment demand of the thermal power leading to the economical dispatch scheme.
Under the robust optimized dispatch scheme, the dayahead operation cost is 2359.9 USD while the adjustment cost under the worstcase scenario is 134.8 USD. To evaluate the robustness based on the adjustable robust parameter constraints, this paper also employs Monte Carlo approach to generate a number of 500 random scenarios. The adjustment cost of each scenario is calculated based on the robust optimized dayahead dispatch scheme. Figure 5 shows that the presented approach can identify the worstcase scenario and it can calculate the adjustment cost under the worstcase.
5.2. Analysis of the Economy of the Dispatch Scheme
The economy of the optimized dispatch scheme is evaluated by the number of 500 scenarios generated by Monte Carlo approach. The adjustable robust parameter in the static robust planning is as the same as the twostage robust planning. The comparisons are listed in Table 1.

Three points can be clearly seen in Table 1. (a) Although the twostage robust optimized dayahead operating cost is higher than that of the deterministic optimization, the total cost has been greatly reduced compare to that of the deterministic optimization (no charge for wind power and PV power curtailment). (b) Although the adjustable robust parameter has been introduced, the dayahead operation cost of the static robust planning remains high. In addition, the wind power curtailment and PV power curtailment still occur. Moreover, the total cost for the static robust model is higher than the twostage robust model. (c) Based on the twostage robust optimization, neither wind power curtailment nor PV power curtailment occurs, which greatly improves the accommodation capability of the renewable energies.
As mentioned in the introduction section, the hydropower has strong adjusting ability. And it should be also pointed out that robust model considers the worstcase scenario, in which the output of both wind and photovoltaic power may be lower than the predicted values. Thus, in the robust scheme, hydrounits retain more capacity of upregulation than that in the deterministic scheme to ensure the economic benefits. As a result, though the adjustment cost of hydropower increases, the total adjustment cost is significantly decreased.
In addition, the twostage robust model formulates a better base generation for thermal units and hydrounits than that of the static robust model based on the iterations between the main problem and the subproblem. Consequently, wind power curtailment cost and PV power curtailment cost are reduced to 0. The renewable energy resources accommodation capability of the power system is greatly enhanced.
5.3. The Importance Analysis of the Hydropower
To reveal the importance of the hydropower in the dispatch scheme, extra dispatch scheme optimizations with different hydropower penetrations have been carried out and evaluated. Four different cases are studied as listed below.
Case 1. Not considering the hydropower.
Case 2. The capacity of the hydropower is 9MW.
Case 3. The capacity of the hydropower is 18MW.
Case 4. The capacity of the hydropower is 30MW.
The results are shown in Figures 6 and 7.
It can be observed that, along with the hydropower penetration increases, the output of the thermal power decreases. As a result, the dayahead operation cost can be reduced. In addition, the cost of the dispatch scheme also has been significantly impacted as shown in Table 2.

Table 2 indicates that in Case 1, thermal power has to deal with the uncertainties of the renewable energies. This point leads to higher adjustment cost. Even in some scenarios, adjustment demand exceeds the adjustment capacity which causes wind power curtailment and PV power curtailment. Therefore, the total cost further increases. In Case 2, hydropower participates in managing power imbalance. It shares the adjustment pressure with thermal power. Therefore, the total adjustment cost is greatly reduced. However, PV power curtailment still occurs in this case. In Case 3 and Case 4, neither wind power curtailment nor PV power curtailment occurs. Compared to Case 3, the decline of the adjustment cost in Case 4 is less. It indicates that when the hydropower capacity reaches to 18MW, the power system has a strong ability to cope with the risks due to the uncertainties of the renewable energy sources. Additionally, the experiment also suggests that, to handle the uncertainties of the wind power and PV power, hydropower significantly outperforms the thermal power.
5.4. The Impact of Adjustable Robust Parameter
The previous experiments are based on the values of the adjustable robust parameter =10 and =6. In this section another two groups of adjustable robust parameter =0, =0 and =24, =24 are also employed to study the impact of the optimal dispatch scheme under different values of Γ. Table 3 shows the comparisons.

Table 3 indicates that when the values of the uncertain parameters are 0, the result of the robust optimization is close to that of the deterministic optimization. However, when the values of the uncertain parameters are 24, the dayahead operation cost and total cost become larger than those of =10, =6. This point shows if the values of the adjustable robust parameter are conservative, the performance of the optimization also becomes conservative.
5.5. The Impact of Prediction Error
In order to reveal the influence of the different prediction errors of the renewable energy sources in the twostage robust model, the deviation indicator of the fluctuation range is introduced. To simplify the problem, it assumes that the wind power and PV power have the same deviation indicator. The formula of the deviation is shown as (33):
where and are the fluctuation range of the wind power and PV power illustrated in the Figures 10 and 11.
In this section, five deviation indicators (0, 0.25, 0.5, 0.75, and 1) are employed to the presented model. The dayahead operation cost under the different prediction error is shown in Figure 8. In addition, the adjustment cost is also greatly influenced as shown in Table 4.



As can be seen from Figure 8, the dayahead operation cost increases with the increasing deviation indicator value. This is because that the worse scenarios are considered in the twostage robust model. Additionally, the increasing adjustment cost also suggests that more generation adjustments are needed.
6. Conclusion
This paper presents a twostage robust model in enabling economical dayahead dispatch considering wind power, photovoltaic power, hydropower, and thermal power. Based on the model decomposition and linearization, the original model can be converted into MP and SP, which are finally handled by C&CG approach. Based on the experimental results, the contributions of the presented model are as follows:(1)The twostage robust model is able to achieve the optimal dayahead dispatch scheme under the worst wind and photovoltaic output scenarios. The scheme is proved to be economy and has great potential to improve the capability of the power system to accommodate wind and photovoltaic power.(2)Hydropower can effectively handle the uncertainties caused by wind power and photovoltaic power. In our model by considering the adjusting ability of the hydropower, the adjustment demand of the thermal unit can be significantly reduced. Additionally, the wind power curtailment and PV power curtailment can also be potentially avoided.(3)The adjustable robust parameter enables the performance adjustment of the optimized dispatch scheme by tuning their values. As a result, the compromise between economy and robust could be achieved under different dispatching requirements.
Taking more factors (demand side management, battery storage system, electric vehicle [31], and so on) into account in windsolarhydropower system is our followup work in the future.
Appendix
See Figures 9, 10, and 11 and Tables 5 and 6.
Data Availability
The IEEE bus model used to support the findings of this study has been deposited in the Department of Electrical Engineering at the University of Washington. The details of the modifications of the models and the customized parameters are stated in the paper. IEEE bus model site is available at http://www.ee.washington.edu/research/pstca/.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this article.
References
 Y. Bai, Y. Wang, Q. Xia, X. Sun, M. Yang, and J. Zhang, “A fullscenario SCED with coordinative optimization of hydrothermalwind power,” Proceedings of the Chinese Society of Electrical Engineering, vol. 33, no. 13, pp. 1–9, 2013. View at: Google Scholar
 L. Ye, X. Qu, Y. Yao et al., “Analysis on intraday operation characteristics of hybrid windsolarhydro power generation system,” Automation of Electric Power System, vol. 42, no. 04, pp. 158–164, 2018. View at: Google Scholar
 F. Arpino, M. Dell'Isola, D. Maugeri, N. Massarotti, and A. Mauro, “A new model for the analysis of operating conditions of microcogenerative SOFC units,” International Journal of Hydrogen Energy, vol. 38, no. 1, pp. 336–344, 2013. View at: Publisher Site  Google Scholar
 F. Calise, R. D. Figaj, N. Massarotti, A. Mauro, and L. Vanoli, “Polygeneration system based on PEMFC, CPVT and electrolyzer: Dynamic simulation and energetic and economic analysis,” Applied Energy, vol. 192, pp. 530–542, 2017. View at: Publisher Site  Google Scholar
 L. Chen, J. Chen, Z. Li et al., “Scheduling strategy of hybrid windphotovoltaichydro power generation system,” Electric Power Construction, vol. 34, no. 12, pp. 1–6, 2013. View at: Publisher Site  Google Scholar
 X. Yang, Q. Chen, M. Wang, and L. Zhang, “Cooperating control for wind farm and hydro power plant based on the fruit fly optimization,” Proceedings of the Chinese Society of Electrical Engineering, vol. 37, no. 18, pp. 5286–5293, 2017. View at: Google Scholar
 X. Yang, W. Wang, B. Xue et al., “On shortterm united optimal operation of wind power, thermal power and water power,” Journal of Hydroelectric Engineering, vol. 37, no. 18, pp. 5286–5293, 2017. View at: Google Scholar
 J. A. Momoh and S. S. Reddy, “Review of optimization techniques for Renewable Energy Resources,” in Proceedings of the 2014 IEEE Symposium on Power Electronics and Machines for Wind and Water Applications, PEMWA 2014, Milwaukee, WI, USA, July 2014. View at: Google Scholar
 W. Gan, X. Ai, J. Fang, Y. Ta, and J. Wen, “Coordinated optimal operation of the wind, coal, hydro, gas units with energy storage,” Transactions of China Electrotechnical Society, vol. 32, pp. 11–20, 2017. View at: Google Scholar
 H. Liu, P. Daoxin, H. Xin et al., “Joint scheduling optimization model of multitype generator with uncertainty,” Water Resources and Power, vol. 34, no. 02, pp. 153–158, 2016. View at: Google Scholar
 L. Li, W. Wu, and D. Shen, “Modeling and research of windpvwater hybrid power generation system,” East China Electric Power, vol. 40, no. 07, pp. 1157–1160, 2012. View at: Google Scholar
 X. Liu and W. Xu, “Economic load dispatch constrained by wind power availability: A hereandnow approach,” IEEE Transactions on Sustainable Energy, vol. 1, no. 1, pp. 2–9, 2010. View at: Publisher Site  Google Scholar
 W. Dong, Q. Wang, and L. Yang, “A coordinated dispatching model for a distribution utility and virtual power plants with wind/photovoltaic/hydro generators,” Automation of Electric Power Systems, vol. 39, no. 9, pp. 75–207, 2015. View at: Google Scholar
 Y. Zou and L. Yang, “Synergetic dispatch models of a wind/PV/hydro virtual power plant based on representative scenario set,” Power System Technology, vol. 39, no. 7, pp. 1855–1859, 2015. View at: Publisher Site  Google Scholar
 Y. Liu, C. Jiang, J. Shen, and J. Hu, “Coordination of Hydro Units with Wind Power Generation Using Interval Optimization,” IEEE Transactions on Sustainable Energy, vol. 6, no. 2, pp. 443–453, 2015. View at: Publisher Site  Google Scholar
 C. Peng, P. Xie, J. Zhan, and H. Sun, “Robust economic dispatch of microgrid using improved bacterial foraging algorithm,” Power System Technology, vol. 38, no. 9, pp. 2392–2398, 2014. View at: Google Scholar
 Y. Yang, W. Pei, W. Deng, H. Xiao, and Z. Qi, “Optimal operation of gridconnected microgrid based on Benders decomposition,” Electric Power Automation Equipment, vol. 34, no. 10, pp. 21–27, 2014. View at: Publisher Site  Google Scholar
 Y. Liu, L. Guo, and C. Wang, “Economic dispatch of microgrid based on two stage robust optimization,” Proceedings of the CSEE, vol. 38, no. 14, pp. 4013–4022, 2018. View at: Google Scholar
 S. Wang, Y. Niu, L. Fang et al., “Dual stage scheduling strategy for microgrid community considering uncertainty of renewable energy,” Power System Protection and Control, vol. 46, no. 17, pp. 89–98, 2018. View at: Google Scholar
 S. Surender Reddy, P. R. Bijwe, and A. R. Abhyankar, “Realtime economic dispatch considering renewable power generation variability and uncertainty over scheduling period,” IEEE Systems Journal, vol. 99, pp. 1–12, 2015. View at: Publisher Site  Google Scholar
 S. Surender Reddy and P. R. Bijwe, “Dayahead and real time optimal power flow considering renewable energy resources,” International Journal of Electrical Power & Energy Systems, vol. 82, pp. 400–408, 2016. View at: Publisher Site  Google Scholar
 C. Ma, C. Dong, Z. Lü et al., “Shortterm trading and optimal operation strategy for commercial virtual power plant considering uncertainties,” Power System Technology, vol. 40, no. 5, pp. 1543–1549, 2016. View at: Google Scholar
 S. S. Reddy, A. R. Abhyankar, and P. R. Bijwe, “Market clearing for a windthermal power system incorporating wind generation and load forecast uncertainties,” in Proceedings of the 2012 IEEE Power and Energy Society General Meeting, PES 2012, USA, July 2012. View at: Google Scholar
 Y. Wang and J. Zhao, “Optimal spinning reserve of power system with wind power penetrated under power market environment,” Power System and Clean Energy, vol. 33, no. 7, pp. 123–128, 2017. View at: Google Scholar
 S. S. Reddy, “Optimal scheduling of windthermal power system using clustered adaptive teaching learning based optimization,” Electrical Engineering, vol. 99, no. 2, pp. 535–550, 2017. View at: Publisher Site  Google Scholar
 S. S. Reddy, P. R. Bijwe, and A. R. Abhyankar, “Optimum dayahead clearing of energy and reserve markets with wind power generation using anticipated realtime adjustment costs,” International Journal of Electrical Power & Energy Systems, vol. 71, pp. 242–253, 2015. View at: Publisher Site  Google Scholar
 C. Zhang, T. Liu, and T. Bu, “Dynamic optimal dispatch of windhydrothermal power system,” Modern Electric Power, vol. 34, no. 04, pp. 8–14, 2017. View at: Google Scholar
 S. Dehghan, N. Amjady, and A. Kazemi, “Twostage robust generation expansion planning: A mixed integer linear programming model,” IEEE Transactions on Power Systems, vol. 29, no. 2, pp. 584–597, 2014. View at: Publisher Site  Google Scholar
 H. Gao, J. Liu, Z. Wei, and Y. Su, “A bilevel robust planning model of active distribution network and its solution method,” Proceedings of the Chinese Society of Electrical Engineering, vol. 37, no. 5, pp. 1389–1400, 2017. View at: Google Scholar
 X. Wu, X. Wang, J. Wang et al., “Economic generation scheduling of a microgrid using mixed integer programming,” Proceedings of the CSEE, vol. 33, no. 28, pp. 1–9, 2013. View at: Google Scholar
 S. Surender Reddy, J. Y. Park, and C. M. Jung, “Optimal operation of microgrid using hybrid differential evolution and harmony search algorithm,” Frontiers in Energy, vol. 10, no. 3, pp. 355–362, 2016. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 Peng Wei and Yang Liu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.