MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi 10.1155/2019/4634131 4634131 Research Article The Integration of Wind-Solar-Hydropower Generation in Enabling Economic Robust Dispatch Wei Peng 1 http://orcid.org/0000-0002-7935-7146 Liu Yang 2 Mauro Alessandro 1 Yalong River Hydropower Development Company Ltd. Chengdu 610000 China 2 College of Electrical Engineering and Information Technology Sichuan University Chengdu 610065 China scu.edu.cn 2019 212019 2019 06 08 2018 03 12 2018 17 12 2018 212019 2019 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.

Currently the renewable energies including wind power and photovoltaic power have been increasingly deployed in power system to achieve contamination free and environmental-friendly 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 two-stage robust model which is able to achieve the optimal day-ahead dispatch strategy in the worst-case scenario of wind and photovoltaic outputs. Because of the strong interactions between the two stages, the original optimization has been decomposed into the day-ahead 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 Big-M method enables the linearization of the dual model. The interacted-iterations between the master problem and the subproblem are successfully implemented which can ultimately figure out the optimal day-ahead dispatch strategy of the power system with conventional and renewable energies.

1. Introduction

To achieve the contamination free and environmental-friendly power generation, clean energies such as wind power generation and photovoltaic power generation have been increasingly deployed in modern power systems . 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 economical-efficiency 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 . 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  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  presents a novel poly-generation system, including cogenerative proton exchange membrane fuel cell, concentrated photovoltaic-thermal 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  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  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 . 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 . Recently, robust optimization has been widely employed in the economical dispatch as it needs neither precisely modeled probability distribution nor large-scale sampling of uncertain variables . Research  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 two-stage robust optimization model with the “min-max-min” 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  place emphasis on the acquisition of adjustment strategy for unbalanced power caused by the prediction error. Research  introduces the multiagent technology in the microgrid to formulate the energy scheduling strategy. Research  presents a method that considers variability and uncertainty cost of renewable energy sources between two consecutive scheduling intervals in the real-time. Research  further links up the day-ahead scheduling and real-time scheduling by a more meticulous generation dispatch mode to formulate a comprehensive strategy. Researches  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 two-stage robust model considering the generation of wind power, photovoltaic power and hydropower in enabling the day-ahead 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 two-stage robust problem and then employs the strong duality theory to handle the max-min 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 IEEE-30 bus system.

The rest of the paper is organized as: Section 2 presents the day-ahead dispatch robust model considering the wind power, photovoltaic power, and hydropower; Section 3 presents the linearization of the model; Section 3 solves the preprocessed two-stage model; Section 4 evaluates the approach and discusses the simulation results; Section 5 concludes the paper.

2. Day-Ahead Dispatch Robust Model Considering Wind Power, Photovoltaic Power, and Hydropower 2.1. The Details of the Two-Stage 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 day-ahead dispatch strategy is an optimization problem based on the identification of the minimum operating cost, which is represented by (1)minxCopexs.t.Aopex=0,Bopex0where Cope(·) is the conventional day-ahead operating cost of the power system, which is mainly the cost of the thermal power; x is the day-ahead dispatch scheme including the status and the outputs of the thermal units and the hydropower units; Aope(·) and Bope(·) represent the day-ahead 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.

Combining the day-ahead operation and the additional adjustment a model is represented by (3). The model can achieve the minimum cost dispatch scheme x under the worst-case scenario considering the uncertainties of wind and photovoltaic outputs.(3)minxCopex+maxuUminyΩx,uCope+u,ys.t.Aopex=0,Bopex0Aope+y=0,Bope+y0

In (3), x is a type of “here-and-now” variable which represents the optimal day-ahead dispatch scheme under the worst-case scenario. In addition, it can hedge against any possible condition within the uncertain sets; y is a type of “wait-and-see" variable which represents the optimal adjustment scheme under the worst-case 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 Cope(·) and Cope+(·); Section 2.4 presents the day-ahead constraints Aope(·), Bope(·) and the adjustment constraints Aope+(·), Bope+(·).

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 (4)PWTG,tPWTG,tpre-PWTG,tflu,PWTG,tpre+PWTG,tfluPPVG,tPPVG,tpre-PPVG,tflu,PPVG,tpre+PPVG,tfluwhere PWTG,tpre is the predicted value of the wind power output at time t; PPVG,tpre is the predicted value of the photovoltaic power output at time t; PWTG,tflu and PPVG,tflu 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 (5)t=1TPWTG,t-PWTG,tprePWTG,tpreΓWTGt=1TPPVG,t-PPVG,tprePPVG,tpreΓPVGwhere T is the simulation period; ΓWTG and ΓPVG 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 day-ahead dispatch scheme x, costs in different stages are considered. In the first stage the cost is mainly the generation cost of thermal power units CTMG. The costs in the second stage mainly include the adjustment cost of thermal units CTMGdelta, the adjustment cost of hydrounits CHDGdelta, the wind power curtailment cost CWTGloss, and the photovoltaic power curtailment cost CPVGloss. The objective function is represented by(6)Cope=CTMG=t=1Ta·PTMG,t2+b·PTMG,t+cCope+=CTMGdelta+CHDGdelta+CWTGloss+CPVGloss=t=1TλTMG·ΔPTMG,t+λHDG·ΔPHDG,t+λWTGloss·PWTG,t-PWTG,tget+λPVGloss·PPVG,t-PPVG,tgetwhere a, b, and c are the coefficient of the thermal power units; PTMG,t is the outputs of the thermal power units at time t; λTMG and λHDG are the penalty prices for the adjustment power of thermal units and hydrounits; λWTGloss and λPVGloss are the penalty prices of wind power curtailment and photovoltaic power curtailment; ΔPTMG,t and ΔPTMG,t are the output adjustment of thermal power and hydropower at time t; PWTG,tget and PPVG,tget are the wind power and photovoltaic power which are injected into the power grid at time t.

2.4. Constraints 2.4.1. Day-Ahead Dispatch Constraints

Day-ahead dispatch constraints in the model represented by (3) contain the following:

(a) Thermal Power Constraints (7) S T M G , t · P T M G m i n P T M G , t S T M G , t · P T M G m a x (8) - R T M G d o w n P T M G , t + 1 - P T M G , t R T M G u p where STMG,t is the 0-1 variable which denotes the unit commitment status of the thermal power at time t; PTMGmax and PTMGmin are the upper limit and the lower limit of the outputs of the thermal units; RTMGup and RTMGdown are the upward and downward ramping capability of the thermal power units.

(b) Hydropower Constraints(9)PHDGminPHDG,tPHDGmax(10)UHDG,t=λHDGa·PHDG,t+λHDGb(11)vflowminUHDG,tvflowmax(12)vdpmminUHDG,t+UHDG,tlossvdpmmax(13)Ut+1=Ut+vin,t-UHDG,t-UHDG,tloss(14)UminUtUmax(15)U1=Uini(16)UT=Utermwhere UHDG,t is the water consumption of the hydrogenerators at time t; PHDG,t is the output of the hydrogenerators at time t; λHDGa and λHDGb are the converting coefficients of hydrogenerators; vflowmax and vflowmin are the upper and lower limit of water consumption in every subperiod; UHDG,tloss is the water curtailment of reservoir; Ut is the reservoir storage at time t; vin,t is the amount of water that flows into the reservoir at time t; Umax and Umin are the upper limit and the lower limit of the reservoir storage. Uini and Uterm are the initial value and final value of the reservoir storage.

(c) Power Balance Constraint (17) P T M G , t + P H D G , t + P W T G , t p r e + P P V G , t p r e = P L O A D , t where PLOAD,t is the power demand of the load in the system at time t.

(d) Branch Flow Constraint (18) b = 1 N b κ l b · P b , t P l m a x where l is the identifier number of the transmission lines; b is the identifier number of the nodes; Nb is the total number of the nodes; κlb is the sensitivity factor of the lth line to the bth node; Pb,t is the net active power which is injected into the bth node at time t; Plmax is the maximum transmission power of the lth line.

Additional adjustment constraints in the model represented by (3) contain the following:

(a) Thermal Power Units Regulation Constraints (19) Δ S T M G , t · Δ P T M G m i n Δ P T M G , t Δ S T M G , t · Δ P T M G m a x S T M G , t · P T M G m i n P T M G , t + Δ P T M G , t S T M G , t · P T M G m a x - R T M G d o w n P T M G , t + 1 + Δ P T M G , t + 1 - P T M G , t - Δ P T M G , t R T M G u p where ΔSTMG,t is the 0-1 variable which represents the adjustment status of the thermal power units at time t; ΔPTMGmax and ΔPTMGmin 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 Constraints (20) Δ S H D G , t · Δ P H D G m i n Δ P H D G , t Δ S H D G , t · Δ P H D G m a x (21) P H D G m i n P H D G , t + Δ P H D G , t P H D G m a x U H D G , r t , t = λ H D G a · P H D G , t + Δ P H D G , t + λ H D G b v f l o w m i n U H D G , r t , t v f l o w m a x v d p m m i n U H D G , r t , t + U H D G , r t , t l o s s v d p m m a x U r t , t + 1 = U r t , t + v i n , t - U H D G , r t , t - U H D G , r t , t l o s s U m i n U r t , t U m a x where ΔSHDG,t is the binary variable which represents the adjustment status of hydrogenerators at time t; ΔPHDGmax and ΔPHDGmin are the upper and lower limit of the hydrogenerator adjustment power; UHDG,rt,t is the water consumption after the adjustment at time t; UHDG,rt,tloss is the water curtailment of reservoir after the adjustment at time t; Urt,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 (22) P T M G , t + Δ P T M G , t + P H D G , t + Δ P H D G , t + P W T G , t g e t + P P V G , t g e t = P L O A D , t

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 Constraint (23) b = 1 N b κ l b · P r t , b , t P l m a x where Prt,b,t is the net active power injected into the bth node after adjustment at time t.

(e) Wind Power Constraint (24) 0 P W T G , t g e t P W T G , t where PWTG,t is the actual output of wind turbine generators at time t, which is an uncertainty variable.

(f) Photovoltaic Power Constraint (25) 0 P P V G , t g e t P P V G , t where PPVG,t is the actual output of photovoltaic generators at time t, which is an uncertainty variable.

3. Linearization of Thermal Power Generation Cost

The day-ahead 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  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, piecewise-linear approximation is employed. Firstly, we segment the curve into a number of N segments. And then in each segment, by employing the continuous variable Pj,t and the state variable Sj,t, linear approximation can been done. As a result, the linearized cost of the thermal units is presented by (26)CTMG,tap=j=1NmjPj,t+njSj,t(27)mj=CTMG,t+1ap-CTMG,tapPj,t+1-Pj,tnj=CTMG,tap-Pj,t·mjPj·Sj,tPj,tPj+1Sj,tSTMG,t=j=1Nnj,t1where CTMG,tap is the linearized cost function of the thermal units; mj is the slope of jth segmentation; nj is the intercept of the jth segmentation; Pj,t is the output of the jth segmentation; Sj,t is the status of the jth segmentation at time t; Pj is the jth segmentation point.

Piecewise linearizing of generation cost curve of thermal units.

Based on the linearization, the original two-stage model can be represented using the matrix shown by (28)OP:minxdTx+maxuUminyeTy+fTus.t.Dxg,Fx+Gyh,Iuyuwhere d, e, and f are the coefficient matrices in the objective function; D, F, G, and Iu are the coefficient matrices in the constraints; g and h are the column vectors; x is the vector of all the first-stage decisions; y is the vector of all the second stage decisions; u 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 day-ahead dispatch scheme in the first-stage 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  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. (29)MP:minxdTx+ηs.t.ηeTy+fTuDxg,Fx+Gyh,Iuyu(30)SP1:maxuUminyeTy+fTus.t.Fx+Gyh,Iuyu

4.2. The Equivalent Maximization Problem of SP1

This paper employs strong duality theory to convert the SP1 into an equivalent maximization problem. Big-M approach is also adopted to handle the nonlinear terms in the dual problem. The converted SP2 is represented in(31)SP2:maxu,α,βhTα-xTFTα+u+θ++u-θ-+u01-θ+-θ-s.t.GTα+IuTβ=e,α0,β0,θ=f+β-M1-μt++θtθt+M1-μt++θt-M1-μt-+θtθt-M1-μt-+θt-Mμt+θt+Mμt+,-Mμt-θt-Mμt-μt++μt-1,tμt++μt-Γwhere α, β, γ are the dual variables; u+,u- are the upper and lower limit of uncertainty variable; u0 is the predictive value of uncertainty variable; θ+, θ- are the positive and negative values of θ; μt+, μt- are the binary variables introduced by the Big-M approach; the adjustable robust parameter Γ includes ΓWTG and ΓPVG.

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 worst-case scenario u1. 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 ul (l=1, 2, …, k), solve the MP and achieve the optimal solution (xk, η,y1,y2, …,yk), and then set L=η.

Step 3 (subproblem computation).

After given the first-stage solutions x = xk, the solution (uko, yko) of the SP2 can be obtained. Set the worst-case scenario uk+1=uko, U=eTyko+fTuko.

Step 4 (convergence check).

If U-Lλ, return the solution x=xk. Otherwise add the new variable yk+1 and the following constraints (32) which are associated with the new worst-case scenario, then update k = k+1 and go to Step 2 iteratively.(32)Add:ηeTyk+1+fTuk+1Fx+Gyk+1h,Iuyk+1uk+1

5. Experimental Result

To evaluate the effectiveness of the presented two-stage robust optimal model in enabling the economical day-ahead dispatch considering the uncertainties, this paper implements the model using MATLAB. The solution of the model is by employing CPLEX. A modified IEEE-30 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 IEEE-30 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 , 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 Two-Stage Robust Optimization

The values of ΓWTG and ΓPVG are set to 10 and 6. During the solution of the two-stage 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.

Convergence of C&CG approach.

The output power of thermal units and hydrounits.

The worst-case scenario of WTG and PVG.

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 8h-24h. 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 8h-24h 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 day-ahead operation cost is 2359.9 USD while the adjustment cost under the worst-case 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 day-ahead dispatch scheme. Figure 5 shows that the presented approach can identify the worst-case scenario and it can calculate the adjustment cost under the worst-case.

Robustness analysis.

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 two-stage robust planning. The comparisons are listed in Table 1.

Comparison for operation costs.

Category Two-stage robust planning Static robust planning Deterministic planning
Avg. Max Avg. Max Avg. Max
Day-ahead operation cost /\$ 2359.9 2374.2 2264.1

Adjustment cost /\$ Thermal units 2.7 30.8 85.0 251.5 375.9 1153.7
Hydro units 123.2 169.4 153.9 221.8 51.1 145.3
Wind power curtailment 0 0 73.0 500.9 314.9 1053.1
PV power curtailment 0 0 53.1 473.2 225.2 813.9

Total cost /\$ 2485.8 2541.4 2739.2 3348.3 3231.2 4098.5

Three points can be clearly seen in Table 1. (a) Although the two-stage robust optimized day-ahead 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 day-ahead 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 two-stage robust model. (c) Based on the two-stage 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 worst-case 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 two-stage 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.

The output power of thermal units in different cases.

Day-ahead operation cost in different cases.

It can be observed that, along with the hydropower penetration increases, the output of the thermal power decreases. As a result, the day-ahead operation cost can be reduced. In addition, the cost of the dispatch scheme also has been significantly impacted as shown in Table 2.

The effect of hydropower.

Category Case 1 Case 2 Case 3 Case 4
Avg. Max Avg. Max Avg. Max Avg. Max
Adjustment cost of thermal units /\$ 404.1 1109.5 21.2 107.9 2.7 30.8 0.5 7.2

Adjustment cost of hydro units /\$ 110.3 143.35 123.2 169.4 124.3 174.0

Wind power curtailment cost /\$ 154.2 1089.0 0 0 0 0 0 0

PV power curtailment cost /\$ 193.3 763.0 4.9 37.3 0 0 0 0

Total cost /\$ 751.6 1852.0 136.34 245.4 125.9 181.5 124.8 176.7

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 ΓWTG=10 and ΓPVG=6. In this section another two groups of adjustable robust parameter ΓWTG=0, ΓPVG=0 and ΓWTG=24, ΓPVG=24 are also employed to study the impact of the optimal dispatch scheme under different values of Γ. Table 3 shows the comparisons.

Optimization results and economic efficiency with different adjustable robust parameter Γ.

Category Γ W T G =0, ΓPVG=0 Γ W T G =10, ΓPVG=6 Γ W T G =24, ΓPVG=24
Avg. Max Avg. Max Avg. Max
Day-ahead operation cost /\$ 2264.1 2359.9 2430.0

Adjustment cost /\$ 969.6 1831.8 125.9 181.5 123.4 169.4

Total cost /\$ 3233.7 4095.9 2485.8 2541.4 2553.4 2599.4

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 day-ahead operation cost and total cost become larger than those of ΓWTG=10, ΓPVG=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 two-stage 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):(33)P°WTG,tflu=α·P¯WTG,tfluP°PVG,tflu=α·P¯PVG,tflu

where P¯WTG,tflu and P¯PVG,tflu 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 day-ahead 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.

Category α=0 α=0.25 α=0.5 α=0.75 α=1
Avg. Max Avg. Max Avg. Max Avg. Max Avg. Max
Adjustment cost of thermal units /\$ 0 0 0 0 0.6 10.6 2.6 28.1 2.7 30.8

Adjustment cost of hydro units /\$ 0 0 31.3 42.3 61.2 83.6 91.2 122.4 123.2 169.4

Wind power curtailment cost /\$ 0 0 0 0 0 0 0 0 0 0

PV power curtailment cost /\$ 0 0 0 0 0 0 0 0 0 0

total adjustment cost /\$ 0 0 31.3 42.3 61.9 83.6 93.8 148.0 125.9 181.5

The parameters of thermal generator.

Category Node P T M G m a x (MWh) P T M G m i n (MWh) R T M G u p (MW/h) R T M G d o w n (MW/h) a b c λ T M G (\$/Mwh)
1 1 40 10 15 10 0.02 2 0 8.152
2 2 40 10 20 15 0.0175 1.75 0 12.170
3 22 25 5 10 10 0.0625 1 0 7.771
4 27 22.5 5 10 8 0.00834 2.25 0 6.899
5 23 15 5 8 5 0.025 2 0 5.024
6 13 20 5 8 5 0.025 2 0 6.020

The parameters of hydrogenerator.

λ H D G a (m3/(Mwh)) λ H D G b (m3/(Mwh)) v f l o w m a x (104m3/h) v f l o w m i n (104m3/h) v d p m m a x (104m3/h) v d p m m i n (104m3/h)
2500 0 4.86 0 7.00 3.25

v i n , t (104m3/h) U m a x (108m3) U m i n (108m3) U i n i (108m3) U t e r m (108m3) λ H D G (\$/Mwh)

4.10 1.3586 1.1623 1.3366 1.3366 3.5

Comparison for the planning results under different prediction error.

Wind power fluctuation interval.

PV fluctuation interval.

As can be seen from Figure 8, the day-ahead operation cost increases with the increasing deviation indicator value. This is because that the worse scenarios are considered in the two-stage robust model. Additionally, the increasing adjustment cost also suggests that more generation adjustments are needed.

6. Conclusion

This paper presents a two-stage robust model in enabling economical day-ahead 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:

The two-stage robust model is able to achieve the optimal day-ahead 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.

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.

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 , and so on) into account in wind-solar-hydropower system is our follow-up 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.

Bai Y. Wang Y. Xia Q. Sun X. Yang M. Zhang J. A full-scenario SCED with coordinative optimization of hydro-thermal-wind power Proceedings of the Chinese Society of Electrical Engineering 2013 33 13 1 9 2-s2.0-84877910440 Ye L. Qu X. Yao Y. Analysis on intraday operation characteristics of hybrid wind-solar-hydro power generation system Automation of Electric Power System 2018 42 04 158 164 Arpino F. Dell'Isola M. Maugeri D. Massarotti N. Mauro A. A new model for the analysis of operating conditions of micro-cogenerative SOFC units International Journal of Hydrogen Energy 2013 38 1 336 344 2-s2.0-84871985173 10.1016/j.ijhydene.2012.10.003 Calise F. Figaj R. D. Massarotti N. Mauro A. Vanoli L. Polygeneration system based on PEMFC, CPVT and electrolyzer: Dynamic simulation and energetic and economic analysis Applied Energy 2017 192 530 542 2-s2.0-84995704429 10.1016/j.apenergy.2016.08.018 Chen L. Chen J. Li Z. Scheduling strategy of hybrid wind-photovoltaic-hydro power generation system Electric Power Construction 2013 34 12 1 6 10.1049/cp.2012.1783 Yang X. Chen Q. Wang M. Zhang L. Cooperating control for wind farm and hydro power plant based on the fruit fly optimization Proceedings of the Chinese Society of Electrical Engineering 2017 37 18 5286 5293 2-s2.0-85037811040 Yang X. Wang W. Xue B. On short-term united optimal operation of wind power, thermal power and water power Journal of Hydroelectric Engineering 2017 37 18 5286 5293 Momoh J. A. Reddy S. S. Review of optimization techniques for Renewable Energy Resources Proceedings of the 2014 IEEE Symposium on Power Electronics and Machines for Wind and Water Applications, PEMWA 2014 July 2014 Milwaukee, WI, USA 2-s2.0-84928612103 Gan W. Ai X. Fang J. Ta Y. Wen J. Coordinated optimal operation of the wind, coal, hydro, gas units with energy storage Transactions of China Electrotechnical Society 2017 32 11 20 2-s2.0-85029680874 Liu H. Daoxin P. Xin H. Joint scheduling optimization model of multi-type generator with uncertainty Water Resources and Power 2016 34 02 153 158 Li L. Wu W. Shen D. Modeling and research of wind-pv-water hybrid power generation system East China Electric Power 2012 40 07 1157 1160 Liu X. Xu W. Economic load dispatch constrained by wind power availability: A here-and-now approach IEEE Transactions on Sustainable Energy 2010 1 1 2 9 2-s2.0-77955689529 10.1109/TSTE.2010.2044817 Dong W. Wang Q. Yang L. A coordinated dispatching model for a distribution utility and virtual power plants with wind/photovoltaic/hydro generators Automation of Electric Power Systems 2015 39 9 75 207 2-s2.0-84929757265 Zou Y. Yang L. Synergetic dispatch models of a wind/PV/hydro virtual power plant based on representative scenario set Power System Technology 2015 39 7 1855 1859 10.13335/j.1000-3673.pst.2015.07.014 2-s2.0-84937402573 Liu Y. Jiang C. Shen J. Hu J. Coordination of Hydro Units with Wind Power Generation Using Interval Optimization IEEE Transactions on Sustainable Energy 2015 6 2 443 453 2-s2.0-85027945493 10.1109/TSTE.2014.2382122 Peng C. Xie P. Zhan J. Sun H. Robust economic dispatch of microgrid using improved bacterial foraging algorithm Power System Technology 2014 38 9 2392 2398 2-s2.0-84907255087 Yang Y. Pei W. Deng W. Xiao H. Qi Z. Optimal operation of grid-connected microgrid based on Benders decomposition Electric Power Automation Equipment 2014 34 10 21 27 10.3969/j.issn.1006-6047.2014.10.004 2-s2.0-84908276480 Liu Y. Guo L. Wang C. Economic dispatch of microgrid based on two stage robust optimization Proceedings of the CSEE 2018 38 14 4013 4022 Wang S. Niu Y. Fang L. Dual stage scheduling strategy for microgrid community considering uncertainty of renewable energy Power System Protection and Control 2018 46 17 89 98 Surender Reddy S. Bijwe P. R. Abhyankar A. R. Real-time economic dispatch considering renewable power generation variability and uncertainty over scheduling period IEEE Systems Journal 2015 99 1 12 10.1109/JSYST.2014.2325967 2-s2.0-84960373595 Surender Reddy S. Bijwe P. R. Day-ahead and real time optimal power flow considering renewable energy resources International Journal of Electrical Power & Energy Systems 2016 82 400 408 10.1016/j.ijepes.2016.03.033 2-s2.0-84964049728 Ma C. Dong C. Z. Duan Q. Wang K. Ai X. Li Y. Short-term trading and optimal operation strategy for commercial virtual power plant considering uncertainties Power System Technology 2016 40 5 1543 1549 2-s2.0-84971287164 Reddy S. S. Abhyankar A. R. Bijwe P. R. Market clearing for a wind-thermal power system incorporating wind generation and load forecast uncertainties Proceedings of the 2012 IEEE Power and Energy Society General Meeting, PES 2012 July 2012 USA 2-s2.0-84870596591 Wang Y. Zhao J. Optimal spinning reserve of power system with wind power penetrated under power market environment Power System and Clean Energy 2017 33 7 123 128 Reddy S. S. Optimal scheduling of wind-thermal power system using clustered adaptive teaching learning based optimization Electrical Engineering 2017 99 2 535 550 2-s2.0-84975503656 10.1007/s00202-016-0382-5 Reddy S. S. Bijwe P. R. Abhyankar A. R. Optimum day-ahead clearing of energy and reserve markets with wind power generation using anticipated real-time adjustment costs International Journal of Electrical Power & Energy Systems 2015 71 242 253 2-s2.0-84926471648 10.1016/j.ijepes.2015.03.002 Zhang C. Liu T. Bu T. Dynamic optimal dispatch of wind-hydro-thermal power system Modern Electric Power 2017 34 04 8 14 Dehghan S. Amjady N. Kazemi A. Two-stage robust generation expansion planning: A mixed integer linear programming model IEEE Transactions on Power Systems 2014 29 2 584 597 2-s2.0-84897661673 10.1109/TPWRS.2013.2287457 Gao H. Liu J. Wei Z. Su Y. A bi-level robust planning model of active distribution network and its solution method Proceedings of the Chinese Society of Electrical Engineering 2017 37 5 1389 1400 2-s2.0-85020217000 Wu X. Wang X. Wang J. Economic generation scheduling of a microgrid using mixed integer programming Proceedings of the CSEE 2013 33 28 1 9 Surender Reddy S. Park J. Y. Jung C. M. Optimal operation of microgrid using hybrid differential evolution and harmony search algorithm Frontiers in Energy 2016 10 3 355 362 2-s2.0-84975260802 10.1007/s11708-016-0414-x