Abstract
This paper proposes an adjustable and distributionally robust chance-constrained (ADRCC) optimal power flow (OPF) model for economic dispatch considering wind power forecasting uncertainty. The proposed ADRCC-OPF model is distributionally robust because the uncertainties of the wind power forecasting are represented only by their first- and second-order moments instead of a specific distribution assumption. The proposed model is adjustable because it is formulated as a second-order cone programming (SOCP) model with an adjustable coefficient. This coefficient can control the robustness of the chance constraints, which may be set up for the Gaussian distribution, symmetrically distributional robustness, or distributionally robust cases considering wind forecasting uncertainty. The conservativeness of the ADRCC-OPF model is analyzed and compared with the actual distribution data of wind forecasting error. The system operators can choose an appropriate adjustable coefficient to tradeoff between the economics and system security.
1 Introduction
Recently, the increased deployment of wind power has brought challenges to electric power systems operation because of its power output uncertainty. To reliably operate the system, stochastic optimization, such as chance-constrained optimal power flow (CC-OPF), has been introduced [1, 2] to minimize the cost of system generation while maintaining reliability and avoiding constraint violations within a defined risk level. Previously, the Gaussian distribution was used to represent wind power forecasting uncertainty [3, 4], but this assumption might lead to an estimation error for the system condition. Distributionally robust CC-OPF (DRCC-OPF) is proposed in [3], which needed only the forecasting mean and variance information, and can lead to a conservative dispatch and high generation cost because the result is robust for any forecasting errors distribution. This violation probability is either higher or lower than the predefined violation level using the predefined distribution assumptions. A higher violation means the system risk level is higher than the preferred level. A lower violation means that the system cost is higher than the optimal value. Therefore, choosing an appropriate coefficient can give a trade-off between the system reliability and cost. In this letter, an adjustable and distributionally robust chance-constrained (ADRCC) optimal power flow (OPF) model is proposed to obtain this optimal coefficient.
When the wind power forecasting uncertainty is taken into consideration, the CC-OPF model is given by:
where \( N \) is the number of generators in the system; \( G_{i} \) is the generation power output; \( D_{i} \) is the load amount; \( P_{i} \) is the forecasted wind power; \( c_{i} \) is the offer price of generation; \( G_{exp,i} \) is the expected generation power output; \( P_{exp,i} \) is the expected wind power; \( \Delta G_{i} \) is the generation response to the wind uncertainty; \( \Delta P_{i} \) is the uncertain output of wind power; \( \beta_{i} \) is the sharing factor of generation for the wind uncertain power output; \( GSF_{l,i} \) is the system generation shift factor; \( LU_{l} \) is the transmission limit; \( G_{{i,{ \hbox{max} }}} \) and \( G_{{i,{ \hbox{min} }}} \) are the generation upper and lower limits, respectively; (3)–(6) are the chance constraints considering the impact of the wind uncertainty on the transmission overload and generation output violation; \( \epsilon \) is the confidence level in the chance constraints (mostly between 1% and 5%). The decision variable in the model is the generation dispatch \( G_{exp,i} \) and the balancing factor \( \beta_{i} \).
Define that the uncertain part of wind power output \( \Delta \varvec{P} \) has the mean as µ and covariance as \( \varvec{\varSigma} \) which is a positive semidefinite matrix. In addition, define \( \varvec{a}_{l} \left(\varvec{\beta}\right) \) as a one-column matrix whose element \( a_{l,i} \left(\varvec{\beta}\right) \) is given by (12), and \( \varvec{b}_{i} \left(\varvec{\beta}\right) \) as a one-column matrix given by (13).
where W is the number of wind power plants.
Then, introduce auxiliary variables \( \eta_{PF,l} \) and \( \eta_{g,i} \) as shown in (14) and (15).
Equations (14) and (15) can be formulated as second-order cone constraints as:
The CC-OPF model is then formulated as (1), (12), (13), (16)–(22).
where \( K_{\epsilon } \) is the coefficient to control the robustness of the chance constraints to wind power forecasting errors. If the forecasting error follows a Gaussian distribution [5], the value of \( K_{\epsilon } \) is decided by (23).
where \( \varPsi \left( x \right) \) is the cumulative distribution function (CDF) of Gaussian distribution.
2 ADRCC-OPF
In (23), \( K_{\epsilon } \) is obtained based on the Gaussian distribution assumption for the wind power forecasting errors, which is not always the case in actual practice, as shown in Fig. 1 [4]. This section discusses the proposed ADRCC-OPF model that does not need the Gaussian distribution assumption for wind power forecasting errors and can control the tradeoff between the reliability and the system cost.
First, the assumption is that the first- and second-order moments of the uncertain variables (mean and covariance) are estimated based on the historical data [2]. A set \( \varvec{p} \) of all the possible distributions \( \varvec{P} \) satisfying the mean and covariance values are represented as:
where \( \varvec{\omega} \) is the uncertain variable (here is the wind power forecasting error); \( P_{0} \left( {{\mathbf{R}}^{\left| v \right|} } \right) \) denotes the set of all of the probabilistic distributions on \( {\mathbf{R}}^{\left| v \right|} \) with mean as µ and covariance matrix as \( \varvec{\varSigma} \); \( \varvec{E}_{\varvec{P}} \left[\varvec{\omega}\right] \) represents the expectation of \( \varvec{\omega} \) [3, 6, 7].
When only the mean and covariance of the distribution are obtained from the historical data, for \( \epsilon \in \left( {0,1} \right) \), the distributionally robust chance constraint is formulated with \( K_{\epsilon } \) value decided by (25) [1].
If the forecasting error distribution is symmetrical with observing the historical forecasting error data, then the DRCC-OPF is formulated [1] with \( K_{\epsilon } \) value decided by (26).
The values of \( K_{\epsilon } \) in the three cases are listed in Table 1 for \( \epsilon \) varying from 2% to 5%.
However, the actual wind power forecasting errors might not be Gaussian or symmetrical. Therefore, using the \( K_{\epsilon } \) listed in Table 1 might lead to over-optimistic or over-conservative results. Using the \( K_{\epsilon } \) values of the distributionally robust case will lead to over-conservative results that has a higher operating cost and is not likely to be adopted by the system operators. Therefore, it is important to choose an appropriate \( K_{\epsilon } \) value leading to a predefined reliability and system cost tradeoff.
For the given \( \epsilon \) value under each distribution assumption (Gaussian, symmetric distributionally robust and distributionally robust), the CC-OPF is performed to obtain the dispatch results and the balancing factors using \( K_{\epsilon } \) values listed in Table 1. With the actual historical wind power forecasting error data, the system actual violations under each distribution assumption are obtained. Comparing the obtained violation value with the predefined value (\( \epsilon \)), the operators can choose the adjustable coefficient as shown in (27) to make the system violation value as close to the predefined value as possible to avoid over-optimistic or over-conservative results shown in Fig. 2.
where \( K_{\epsilon 1} \) and \( K_{\epsilon 2} \) are the coefficients leading to the violation higher and lower than the predefined value, respectively; \( V_{1} \) and \( V_{2} \) are the realized violation values. Here a linear approximation is used for simplicity because the forecasting errors have a linear relationship with \( K_{\epsilon } \) shown in the chance constraints in (19)–(21).
Note that after obtaining the generation dispatch and the balancing factors, the historical forecasting error data are used to calculate the system violation. In this process, the computation procedure can be paralleled to obtain the violation values under different distribution assumption fast and efficiently.
In the system operation, the system operators can choose the appropriate \( K_{\epsilon } \) based on the actual historical data for each operation interval using the actual system realized wind power deviation. Then in the actual operation, they can choose the appropriate coefficients based on the offline simulation. Similar mechanism can be found in Electricity Reliability Council of Texas (ERCOT) dynamic reserve curve [8].
3 Case studies
The test system has been modified from the PJM 5-bus system depicted in Fig. 3. The numbers besides the generators represent the generation offer price ($/MWh) and capacity (MW). The system parameters are from [9]. In this study, the peak load in this system is 1500 MW in the studied interval, and the total load is equally distributed among buses B, C, and D. Two wind power plants with 300 MW as their forecasted mean power output are located at buses B and C. One-year 8760-hour historical wind forecasting error data are utilized from California Independent System Operator (CAISO) to obtain the forecasting error mean and covariance [4], which are also utilized to examine the system violation with DRCC-OPF.
Table 2 lists the system generation cost with different distribution assumptions and confidential levels from 2% to 5%. It is obvious that the cost increases with the confidence level. In Gaussian distribution assumption, the cost increase from $15094.8 to $15170.8 when \( \epsilon \) changes from 5% to 2%. In the distributionally robust case, the cost increases from $15597.7 to $16332.5. In addition, the cost increases with the robustness of the wind power forecasting error distribution assumption. For instance, the cost in the distributionally robust case is the highest, the symmetrically distributionally robust case has the second highest cost, and the Gaussian distribution case has the lowest generation cost. It is clear from Table 1 that higher confidence level (lower \( \epsilon \)) and higher robustness of wind power forecasting error distribution lead to a higher \( K_{\epsilon } \). However, this comes at a cost.
Figure 4 shows the generation dispatch with \( \epsilon \) of 5%. The generation changes with the conservativeness of the wind power forecasting error. In the distributionally robust case, Gen4 and Gen5 reduce their power outputs while Gen3 increases its power to mitigate the wind power uncertainty. This leads to the increment of system generation cost, since Gen3 has a higher cost than Gen5, as shown in Fig. 3.
After the generation dispatch and sharing factors are obtained from the proposed ADRCC-OPF, the historical 8760 wind power forecasting error scenarios are utilized to represent the actual wind power realization and the system violation such as generation violation is calculated [10]. Table 3 shows the output violation of Gen4 under different confidence levels and various distribution assumptions. Figures 5 and 6 are the power flow realization of line 6 and the power output of Gen4 considering the wind power forecasting error.
Table 3 demonstrates that the Gaussian distribution assumption will lead to an optimistic result. For instance, when the confidence level is 1%, Gen4 has a 4.64% probability of violating its generation upper limit. In actual system operation, this will bring risk to system operations. The symmetrically distributionally robust and distributionally robust cases lead to a conservative dispatch because the actual violation possibility of Gen4 is lower than the predefined \( \epsilon \) value in these cases. This over-conservativeness leads to a higher generation cost. Therefore, appropriately choosing the distribution assumption of wind power forecasting error (the value of \( K_{\epsilon } \)) to represent the predefined risk level is of importance for the system operators to tradeoff between the system reliability and the cost. Again, this is also the reason that the proposed method is adjustable in addition to being distributionally robust. In this table, the adjustable coefficients \( K_{\epsilon } \) leading to the violation probability close to the predefined \( \epsilon \) value are also listed. The value of \( K_{\epsilon } \) is obtained using (27). The violation in Gaussian distribution is \( V_{1} \), and the violation in symmetrical distributionally robust case is \( V_{2} \). The values of \( K_{\epsilon } \) in Table 1 for Gaussian and symmetrical distributionally robust are \( K_{\epsilon 1} \) and \( K_{\epsilon 2} \).
In Fig. 5, it is obvious that the transmission power flow will reduce to a narrow range, and that the maximum potential power flow realization reduces, which also leads to lower transmission power flow violations. The black dashed line in Fig. 6 is the upper limit of Gen4. Figure 6 clearly shows that with the conservativeness increases the generation power output will reduce which leads to a lower generation output violation.
The IEEE 118-bus system is tested to verify the proposed model in a large-scale system. The system data is from [9, 11]. Four wind power plants are connected at bus 2, bus 5, bus 53, and bus 86 with the expected power output equal to 300 MW for the studied interval. The confidence level is 5%. The system generation cost and output violation probabilities of Gen2 are listed in Table 4. In this case, the \( K_{\epsilon } \) that leads to the closest value to the defined confidence level should be within the values of the symmetrically distributionally robust case and the distributionally robust case in Table 1.
4 Conclusion
This work proposes a new chance-constrained OPF model called ADRCC-OPF. The influence of the distributional robustness for wind power forecasting on the generation cost and system violations is studied. It shows that the Gaussian distribution assumption for wind power forecasting can lead to overoptimistic results for the system violations and bring risk to system operation. The traditional distributionally robust CC-OPF model will lead to over-conservative results with a high operating cost. It is demonstrated in this work that the proposed ADRCC-OPF can achieve an appropriate tradeoff between the conservativeness and the system cost when an appropriate adjustable coefficient is chosen to control the robustness of the chance constraints.
References
Calafiore GC, El Ghaoui L (2006) On distributionally robust chance-constrained linear programs. J Optim Theory Appl 130(1):1–22
Jiang R, Guan Y (2016) Data-driven chance constrained stochastic program. Math Programm 158(1–2):291–327
Xie W, Ahmed S (2018) Distributionally robust chance constrained optimal power flow with renewables: a conic reformulation. IEEE Trans Power Syst 33(2):1860–1867
Hodge BM, Florita A, Orwig K et al (2012) A comparison of wind power and load forecasting error distributions. In: Proceedings of 2012 world renewable energy forum, Denver, USA, 13–17 May 2012, pp 1–8
Jabr RA (2013) Adjustable robust OPF with renewable energy sources. IEEE Trans Power Syst 28(4):4742–4751
Fang X, Hodge BM, Du E et al (2018) Modelling wind power spatial-temporal correlation in multi-interval optimal power flow: a sparse correlation matrix approach. Appl Energy 230:531–539
Fang X, Hodge BM, Du E et al (2019) Introducing uncertainty components in locational marginal prices for pricing wind power and load uncertainties. IEEE Trans Power Syst. https://doi.org/10.1109/tpwrs.2018.2881131
ERCOT market training website. http://www.ercot.com/content/wcm/training_courses/107/ordc_workshop.pdf. Accessed 1 May 2018
Fang X, Li F, Wei Y et al (2016) Strategic scheduling of energy storage for load serving entities in locational marginal pricing market. IET Gen Transm Distrib 10(5):1258–1267
Hodge BM, Milligan M (2011) Wind power forecasting error distributions over multiple timescales. In: Proceedings of IEEE PES general meeting, Detroit, USA, 24–29 July 2011, pp 1–8
Fang X, Hu Q, Li F et al (2016) Coupon-based demand response considering wind power uncertainty: a strategic bidding model for load serving entities. IEEE Trans Power Syst 31(2):1025–1037
Acknowledgements
This work was co-authored by Alliance for Sustainable Energy, LLC, the manager and operator of the National Renewable Energy Laboratory for the U.S. Department of Energy (DOE) (No. DE-AC36-08GO28308). The funding was provided by U.S. DOE Office of Energy Efficiency and Renewable Energy Wind Energy Technologies Office. The views expressed in the article do not necessarily represent the views of the DOE or the U.S. Government. The U.S. Government retains and the publisher, by accepting the article for publication, acknowledges that the U.S. Government retains a nonexclusive, paid-up, irrevocable, worldwide license to publish or reproduce the published form of this work, or allow others to do so, for U.S. Government purposes.
Author information
Authors and Affiliations
Corresponding author
Additional information
CrossCheck date: 25 January 2019
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
FANG, X., HODGE, BM., LI, F. et al. Adjustable and distributionally robust chance-constrained economic dispatch considering wind power uncertainty. J. Mod. Power Syst. Clean Energy 7, 658–664 (2019). https://doi.org/10.1007/s40565-019-0526-5
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s40565-019-0526-5