Next Article in Journal
Erratum: Algan, E.K.; Ummanel, A. Toward Sustainable Schools: A Mixed Methods Approach to Investigating Distributed Leadership, Organizational Happiness, and Quality of Work Life in Preschools. Sustainability 2019, 11, 5489
Previous Article in Journal
Motivations and Attitudes: An Empirical Study on DIY (Do-It-Yourself) Consumers in Hungary
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Probabilistic Optimal Power Flow for Day-Ahead Dispatching of Power Systems with High-Proportion Renewable Power Sources

1
School of Electrical Engineering and Automation, Harbin Institute of Technology, Harbin 150001, China
2
Electric Power Research Institute of HITZ, Harbin Institute of Technology at Zhangjiakou, Zhangjiakou 075000, China
3
Shenyang Power Supply Company, State Grid Liaoning Electric Power Supply Co., Ltd., Shenyang 110000, China
*
Author to whom correspondence should be addressed.
Sustainability 2020, 12(2), 518; https://doi.org/10.3390/su12020518
Submission received: 1 December 2019 / Revised: 30 December 2019 / Accepted: 8 January 2020 / Published: 9 January 2020
(This article belongs to the Section Energy Sustainability)

Abstract

:
With the increasing proportion of uncertain power sources in the power grid; such as wind and solar power sources; the probabilistic optimal power flow (POPF) is more suitable for the steady state analysis (SSA) of power systems with high proportions of renewable power sources (PSHPRPSs). Moreover; PSHPRPSs have large uncertain power generation prediction error in day-ahead dispatching; which is accommodated by real-time dispatching and automatic generation control (AGC). In summary; this paper proposes a once-iterative probabilistic optimal power flow (OIPOPF) method for the SSA of day-ahead dispatching in PSHPRPSs. To verify the feasibility of the OIPOPF model and its solution algorithm; the OIPOPF was applied to a modified Institute of Electrical and Electronic Engineers (IEEE) 39-bus test system and modified IEEE 300-bus test system. Based on a comparison between the simulation results of the OIPOPF and AC power flow models; the OIPOPF model was found to ensure the accuracy of the power flow results and simplify the power flow model. The OIPOPF was solved using the point estimate method based on Gram–Charlier expansion; and the numerical characteristics of the line power were obtained. Compared with the simulation results of the Monte Carlo method; the point estimation method based on Gram–Charlier expansion can accurately solve the proposed OIPOPF model

1. Introduction

The steady state analysis (SSA) of power systems is the most fundamental method for the investigation of power system planning and operations. In particular, a calculation method for the steady state operation of power systems is based on power flow. The injection power and voltage of traditional power flow calculation are both determinate quantities, and the bus voltage and line power obtained are also determinate quantities [1,2]. With an increase in proportion of uncertain power sources such as wind and solar power generation in the power grid, there is a corresponding increase in the number of indeterminate variables in the power grid. A grid with a high proportion of uncertain power sources is referred to as a power system with high proportions of renewable power sources. Due to the indeterminate variables, in accordance with the law of probability and statistics, the probabilistic optimal power flow (POPF) model is suitable for the SSA of power systems with high proportions of renewable power sources (PSHPRPSs) [3,4,5,6,7].
After the proposal of the POPF, significant research attention was directed toward the development of the model and calculation method [8,9,10,11,12,13,14]. The POPF model mainly includes the direct-current (DC) POPF model [3], the linear alternating-current (AC) POPF model [8], the piecewise linear AC POPF model [9], and the nonlinear AC POPF model [13]. The calculation methods of POPF models include cumulant, point estimation method, and Monte Carlo simulation method. Moreover, a method of combined cumulants and the Gram–Charlier expansion was proposed. For the DC power flow model, the relationship between the target and control variables can be directly obtained, and the operation speed can be improved [15]. The point estimation method was proposed to solve the AC probabilistic power flow model [16]. Literature [17] applies the two-point estimation method to account for the uncertainty in the optimal power flow problem in the context of the competitive power market. Although the Monte Carlo method demonstrates a high accuracy, it is typically employed as a verification method for other probabilistic power flow algorithms, due to the significant computational load required [18]. Literature [19] proposes a new possibilistic probabilistic power flow that exploits artificial neural networks to model the mathematical relationship between input and output variables and solve the non-linear equation set of power flow problems. Literature [20] proposes a novel method that utilizes radial basis functions (RBF) neural networks in order to find the input/output relationship of power-flow nonlinear equation sets for the microgrids, which does not depend on the partial derivatives of equations and Jacobian matrix (JM) inversions. Literature [21] proposes a novel method suitable for calculating the ill-conditioned network power flow, which utilizes neural networks to find the relationship between inputs and outputs of load-flow equations. Furthermore, the POPF is highly applicable to the optimization of grid dispatching control systems and the reduction of the generation cost or line loss [22,23,24]. In several studies, the use of the cumulants method for the solution of the POPF was proposed [25,26,27,28]. Uncertain renewable energy sources such as solar power and wind power increase the uncertain variables in the grid. For grid research with wind or solar power, probabilistic power flow models are commonly used [29,30,31]. Literature [32] reviews the latest research progress in design, planning, and control problems in the field of renewable and sustainable energy, and concludes that heuristic approaches, Pareto-based multi-objective optimization, and parallel processing are promising research areas. Literature [33] presents the model for finding an optimal joint bid of hydroelectric systems and wind parks on a day-ahead electricity market, which can be combined with hydropower system to make day-ahead plan for the power grid with wind power. In the PSHPRPSs, the power prediction error of uncertain power sources is significantly high. The POPF model of PSHPRPSs should consider the actual situation of significantly high errors in the power generation prediction of uncertain power sources.
The maintenance of power balance is the main objective of grid dispatching control. Because the PSHPRPSs contains a large number of uncertain power sources, only relying on automatic generation control (AGC) cannot balance the power forecast error of day-ahead dispatching. For the safe and steady operation of the power grid, the addition of real-time dispatching links between the day-ahead dispatching and AGC is required, to accommodate uncertain power sources and limit the errors within the range of the AGC [34,35,36,37,38]. Due to the large day-ahead power prediction error of PSHPRPSs, the day-ahead generation plan is inaccurate, and the day-ahead active power prediction error is accommodated by the real-time dispatching and AGC. The reason for the power forecast error of PSHPRPSs is the inherent uncertainty of uncertain renewable energy sources, not the power flow model. Therefore, the use of AC power flow model for day-ahead dispatching of PSHPRPSs cannot improve the calculation accuracy but also complicates the power flow models. However, the probabilistic DC power flow model will lead to large calculation errors. On this basis, in this paper, the once-iterative probabilistic optimal power flow (OIPOPF) model is proposed as the POPF model for the day-ahead dispatching SSA of PSHPRPSs. The power flow model proposed in this paper not only retains the same calculation accuracy as AC power flow, but also has the same operation efficiency as the DC power flow model. The model was solved in this study using the point estimate method based on the Gram–Charlier expansion, and the numerical characteristics of the target variables were obtained. The OIPOPF mathematical model and its solution algorithm were validated on the modified Institute of Electrical and Electronic Engineers (IEEE) 39-bus test system and modified IEEE 300-bus test system. Moreover, the simulation results verify the accuracy of the model and the effectiveness of the method.
The remainder of this paper is structured as follows. Section 2 presents the mathematical model of OIPOPF and details the solution of the proposed OIPOPF using the point estimation method based on the Gram–Charlier expansion. Section 3 presents the simulation results for the OIPOPF, as applied to the IEEE 39-bus test system, to verify the accuracy of the OIPOPF mathematical model. Thereafter, the conclusions of this study are presented in Section 4.

2. Materials and Methods

2.1. POPF Mathematical Model

Renewable energy sources such as wind power and solar power have uncertainties based on their randomness and intermittent nature. For the PSHPRPSs whose bus variables obey probability and statistics law, it is suitable to use the POPF for grid SSA. In this paper, an OIPOPF mathematical model is used to formulate a day-head dispatching plan for the PSHPRPSs.
For a PSHPRPSs with n buses, the bus active power difference, Δ P , is defined as expressed by Equation (1):
Δ P = P S P e 0 ,
where Δ P is a 1 × n 1 dimensional vector; P S is the bus active power; and P e 0 is the initial calculated value of the bus active power. P S = ( P S , 1 , P S , 2 , , P S , i , , P S , n ) , P S , i is the i th variable of P S vector. If P S , i is an uncertain renewable energy source, its mean value is μ S , i , and its standard deviation is σ S , i . Moreover, P e 0 is calculated using AC power flow equations, which can reduce the error in the calculation of Δ P .
The voltage phase angle error calculated by the modified active power error is more accurate. The grid bus active power difference, Δ P , and the bus voltage phase angle difference, Δ θ , satisfy the functional relationship expressed by Equation (2):
Δ P = B Δ θ ,
where Δ θ is a 1 × n 1 dimensional vector. The bus voltage phase angle difference Δ θ is defined as expressed by Equation (3):
Δ θ = θ S θ e ,
where θ S is the bus voltage phase angle, and θ e is the calculated value of bus voltage phase angle.
In Equation (2), the coefficient matrix, B , is an n 1 × n 1 dimensional matrix with elements as expressed by Equation (4):
{ B i j = 1 x i j B i i = j i j i 1 x i j .
Given that B is an invertible matrix, Equation (2) is transformed into Equation (5):
Δ θ I = E · Δ P I .
In Equation (5), E = B 1 , and E is a n 1 × n 1 dimension matrix.
For any branch, l , the from-bus is the i bus, and the to-bus is the j bus. The corresponding voltage phase angle of bus i is θ i , and that of bus j is θ j . Moreover, θ i and θ j are the i th and j th elements in the θ e vector, respectively. The active power of the line between bus i and bus j is defined as expressed by Equation (6):
P i j = V i 2 y i 0 + V i 2 G i j V i V j ( G i j cos θ i j + B i j sin θ i j ) ,
where B i j represents the reactance between bus i and bus j ; G i j represents the conductance between bus i and bus j ; and y i 0 represents the admittance between bus i and the zero-potential bus.
For a grid with L lines, the line power is P L = ( p l 1 , p l 2 , , p i j , , p L ) T .
The day-ahead dispatching of PSHPRPSs considers the distribution of active power in the power grid. Therefore, according to Equation (6), the line power of the OIPOPF is defined as expressed by Equation (7):
P L = f ( θ 1 , θ 2 , , θ i , θ n ) .
The day-ahead dispatching of PSHPRPSs not only needs to set generation plan and reserves accommodation reserve power for the grid, but also needs to reserve line active power margin for the grid to prevent line power congestion during real-time dispatching. Hence, the line power in the day-ahead dispatching plan of PSHPRPSs should be minimized, so as to reserve line power margin. The OIPOPF model should be defined as expressed by Equation (8), to constrain the minimization of line power.
min h ( P i ) = min l = 1 L P l ,
where P i represents the power of power sources, and P l represents the line power.
The OIPOPF mathematical model for the day-ahead dispatching of PSHPRPSs is constrained by power balance constraints, tie-line power constraints, and reserve capacity constraints.

2.1.1. Power Balance Constraints

P d = P u + P c
Equation (9) is the power balance expression; where P d represents the total load power of the power grid, P u represents the total power of the uncertain power generation, and P c represents the total power generated by conventional power sources such as thermal power and hydropower sources.

2.1.2. Tie-Line Power Constraints

The maximum line power is P l max , and the power of each line in Equation (7) is P l . Moreover, P l represents the l th line power. Equation (7) satisfies the line power constraint expression as shown in Equation (10).
P l P l max ( i = 1 , 2 , , L )

2.1.3. Reserve Capacity Constraints

The number of power sources in the grid is assumed to be m . Moreover, P i represents the i th power source. The power source described in P i satisfies the reserve capacity constraints in Equation (11).
P i min P i P i max ( i = 1 , 2 , m )

2.2. Point Estimation Method Based on Gram–Charlier Expansion

2.2.1. Moments, Central Moments, and Their Relationship

For a continuous random variable X of any positive integer k , the x k function can be integrated on ( , + ) with respect to f ( x ) . The k order moment of X can be expressed by Equation (12) [39]:
α k = + x k f ( x ) d x ,
where α k represents the k order moment.
Moreover, the mean value of the continuous random variable X is μ , and the k order central moment of X can then be defined as expressed by Equation (13):
β k = + ( x μ ) k f ( x ) d x ,
where β k represents the k order central moment.
The k order moment, α k , and the k order central moment, β k , of the continuous random variable X satisfy the relationship expressed by Equation (14) [40]:
β k = j = 0 k C k j α k j ( α 1 ) j , ( k = 2 , 3 , , n ) ,
where the zero-order origin moment α 0 = 1 . Moreover, the coefficient C k j satisfies Equation (15):
C k j = k ! j ! ( k j ) ! .
According to Equation (14), the k order central moment, β k , can be obtained by the k order moment, α k . The relationship between the central moment and the moment satisfies the relationship expressed by Equation (16):
{ β 1 = 0 β 2 = α 2 α 1 2 β 3 = α 3 3 α 2 α 1 + 2 α 1 3 β 4 = α 4 4 α 3 α 1 + 6 α 2 α 1 2 3 α 1 3 β 5 = α 5 5 α 4 α 1 + 10 α 3 α 2 10 α 2 α 1 + 5 α 1 5 .

2.2.2. Point Estimation Method

For a grid with n buses, P i represents the i th bus active power; with a mean of μ i standard deviation of σ i . In the point estimation method, for n 1 input variables ( P 1 , P 2 , , P i , , P n 1 ) , each variable contains j points. From the first variable to the n 1 th variable, the POPF is calculated when only one point is considered for one variable at a time and the other variables are mean. For example, for the i th input variable, the POPF calculated at the j th estimation point is expressed by Equation (17) [41]:
P L i , j = F ( μ 1 , μ 2 , , P i , j , , μ n 1 ) .
The relationship between P i and μ i , σ i satisfies the relationship expressed by Equation (18).
P i , j = μ i + η i , j σ i
In Equation (18), η i , j represents the standard central moments, which are defined as expressed by Equation (19):
η i , j = { λ i , 3 2 + ( 1 ) 3 j λ i , 4 3 4 λ i , 3 2   ,   j = 1 , 2         0        ,   j = 3 .
In Equation (19), λ is calculated by the ratio of the central moment, M ( P i ) , to the standard deviation, σ i , which is defined as expressed by Equation (20).
λ i , t = M t ( P i ) ( σ i ) t   ( i = 1 , 2 , 3 , , n )   ,
where M ( P i ) represents the central moment, which is defined as expressed by Equation (13).
After calculating each point in turn, the grid line power, P L , and the line power, P L i , j , calculated by each point estimation satisfy the relationship expressed by Equation (21).
P L t = i = 1 n 1 j = 1 J w i , j ( P L i , j ) t ,
where w i , j represents the weight coefficient, which is defined as expressed by Equation (22).
w i , j = { ( 1 ) 3 j η i , j ( η i , 1 η i , 2 )   ,   j = 1 , 2 1 n 1 η i , 4 η i , 3 2   ,   j = 3

2.2.3. Gram–Charlier Expansion

Based on the Gram–Charlier expansion, the cumulative distribution function (CDF) of a continuous random variable, X , is expressed by Equation (23), and the probability distribution function (PDF) is expressed by Equation (24) [42].
H ( x ) = Φ ( x ) + ξ 1 1 ! Φ ( x ) + ξ 2 2 ! Φ ( x ) + ξ 3 3 ! Φ ( x ) + ,
h ( x ) = φ ( x ) + ξ 1 1 ! φ ( x ) + ξ 2 2 ! φ ( x ) + ξ 3 3 ! φ ( x ) + ,
where Φ ( x ) and φ ( x ) are the CDF and the PDF of the standard normal distribution, respectively; and ξ i is a constant coefficient as shown in Equation (25).
{ ξ 1 = ξ 2 = 0 ξ 3 = β 3 σ 3 ξ 4 = β 4 σ 4 3

2.2.4. Algorithm Calculation Steps

The steps for calculating the OIPOPF using the point estimation method based on the Gram–Charlier expansion are as follows:
(1)
The standard central moment, λ , of the probability variable is calculated according to Equation (20).
(2)
The standard location, η , of the probability variable is calculated according to Equation (19).
(3)
The weight, w , of the probability variable is calculated according to Equation (22).
(4)
The estimated point, P i , j , is calculated according to Equation (18).
(5)
The OIPOPF, P L i , j , is calculated according to Equation (17) and its constraints.
(6)
After calculating the OIPOPF for all the estimated points, the moment, α , of each line power is calculated according to Equation (21).
(7)
The central moment, β , of the line power is calculated according to Equation (14) based on the moment, α , of each line power.
(8)
The Gram–Charlier series coefficient, ξ , is calculated according to Equation (25) based on the central moment, β , of each line power.
(9)
The CDF and PDF of the power of each line are calculated according to Equations (23) and (24).
A flow chart for the calculation of the OIPOPF based on the Gram–Charlier expansion point estimation method is presented in Figure 1.

3. Results

In this study, the OIPOPF mathematical model of PSHPRPSs was applied to the modified IEEE 39-bus test system and the modified IEEE 300-bus test system for verification. In this section, a comparison of the OIPOPF calculation results with those of the DA power flow model and AC power flow model is presented, to verify the accuracy of the OIPOPF mathematical model. The OIPOPF model proposed in this paper is solved by the point estimation method based on the Gram–Charlier expansion. The results were compared with those obtained using the Monte Carlo method. The simulations were conducted using MATLAB 8.3.0.532 (R2014a, TheMathWorks, Inc, Natick, MA, USA) on a personal computer with an i3-4170 CPU, 3.70 GHz processor, and 4.0 GB of random-access memory (RAM).

3.1. Modified Power Test System Raw Data

3.1.1. Modified IEEE 39-Bus Test System

The modified IEEE 39-bus test system diagram is presented in Figure 2. As shown in Figure 2, the modified IEEE 39-bus test system adds uncertain renewable energy sources such as wind and photovoltaic power.
The IEEE 39 bus test power system contains 10 power sources. In the modified test system, four of the power sources were replaced by uncertain renewable energy sources. The numerical characteristics such as the mean and standard deviation of each uncertain renewable energy sources prediction error are shown in Table 1.

3.1.2. Modified IEEE 300-Bus Test System

The IEEE 300 bus test power system contains 69 power sources. In the modified test system, 28 power sources replaced were by uncertain renewable energy sources. The numerical characteristics such as the mean and standard deviation of each uncertain renewable energy source prediction error are shown in Table 2.

3.2. Comparison of Simulation Results with DC Power Flow Model Results

The day-ahead dispatching of PSHPRPSs considers the distribution of active power in the power grid. In this study, the once-iteration power flow model, DC power flow model, and AC power flow model were used to determine the active power of power lines in the power grid, and the results were compared to verify the accuracy of the OIPOPF. In order to verify the accuracy of the calculation results of the power flow model proposed in this paper, the calculation results of the OIPOPF model and the DC power flow model are compared with those of the AC power flow model. The relative error of the simulation results of the two models relative to the AC power flow model is defined by Equation (26).
e r r o r = | P L A C P L 1 | P L A C × 100 % ,
where e r r o r represents the percentage error of the simulation results; P L A C represents the simulation results of the AC power flow model; P L 1 represents the simulation results of the OIPOPF model or DC power flow model.
Table A1 (Appendix A) presents the percentage errors of the OIPOPF mathematical model and AC power flow simulation results for the modified IEEE 39-bus test system. Table A2 (Appendix A) presents the percentage errors of the DC power flow model and AC power flow model simulation results for the modified IEEE 39-bus test system.
Based on the percentage error, the average error was calculated. The average error e r r o r ¯ is defined as expressed by Equation (27), and N L represents the total number of percentage errors.
e r r o r ¯ = e r r o r N L
According to Table 3, the average relative error of the AC power flow and DC power flow calculation result for the modified IEEE 39-bus test system is 4.45%; the average relative error of the AC power flow and the OIPOPF model calculation result for the modified IEEE 39-bus test system is 0.64%. The average relative error of the AC power flow and DC power flow calculation result for the modified IEEE 300-bus test system is 11.91%; the average relative error of the AC power flow and the OIPOPF model calculation result for the modified IEEE 300-bus test system is 0.90%. Regardless of whether the modified IEEE 39-bus test system or the modified IEEE 300-bus test system was used, the simulation results show that the results of the optimal power flow model proposed in this paper are more accurate than those of the DC power flow model, and similar to those of the AC power flow model. Furthermore, in the PSHPRPSs, the uncertain power source accounts for a relatively high proportion, and its power prediction error is generally 25–40% [43]. The relative average error of the OIPOPF and AC power flow models for the modified IEEE 39-bus test system and the modified IEEE 300-bus test system are 0.64% and 0.90%, respectively, which is significantly less than the power prediction error of PSHPRPSs. Therefore, the simulation results prove the accuracy of the OIPOPF model proposed in this paper.
In order to verify that the probabilistic power flow model proposed in this paper can improve computing efficiency, the running time of the OIPOPF mathematical model, the DC power flow model and the AC power flow model are compared. The comparison results are shown in Table 4.
As shown in Table 4, regardless of whether the modified IEEE 39-bus test system or the modified IEEE 300-bus test system was used, the operation time of the OIPOPF model is similar to that of the DC power flow, and shorter than that of the AC power flow. Therefore, the simulation results show that the proposed model can maintain high operation efficiency.
To sum up, using the OIPOPF mathematical model to set the day-ahead dispatching plan of PSHPRPSs can not only maintain high calculation accuracy but also have high operation efficiency.

3.3. Simulation Results of OIPOPF Model

In the modified IEEE 39-bus test system and the modified IEEE 300-bus test system, point estimation methods based on Gram–Charlier expansion and Monte Carlo simulation method were employed for the calculation of the OIPOPF model proposed in this paper. The Monte Carlo simulation method was found to be more accurate, and it is generally used to verify the accuracy of the simulation results of other models or algorithms. In this study, the simulation results of the two abovementioned methods were compared to verify that the point estimation method based on Gram–Charlier expansion can be used to solve the proposed OIPOPF model. The mean and standard deviation of the line power calculated using these two simulation methods are shown in Table A3 (Appendix B).
The simulation results of the point estimation method based on Gram–Charlier expansion are represented by A , and the simulation results of the Monte Carlo simulation method are represented by A m c s . The relative error, ε , of the simulation result is expressed by Equation (28):
ε = | A A m c s | | A m c s | .
The relative error, ε , can represent the mean relative error, ε μ , or the standard deviation relative error, ε σ .
To verify the accuracy of the point estimation method based on Gram–Charlier expansion, the cumulative distribution average root mean square error (CDARMSE) was defined as expressed by Equation (29).
C D A R M S E = i = 1 N ( A i A m c s , i ) 2 N ,
where N represents the number of point estimates; A i represents the cumulative distribution value obtained using the point estimation method based on Gram–Charlier expansion for the i th point estimation; and A m c s represents the cumulative distribution value obtained using the Monte Carlo simulation method for the i th point estimation.
Table A4 (Appendix B) presents the relative errors and CDARMSE values of the line power for the modified IEEE 39-bus test system calculated using the Gram–Charlier point estimation and Monte Carlo simulation methods. As seen in Table A4, the CDARMSE values of the two methods were generally less than 0.001. Hence, the point estimation method based on Gram–Charlier expansion can be used to solve the proposed OIPOPF model. Table A4 indicates that the CDARMSE of the line power between busses 31 and 6 was the highest. The PDF and CDF of the line power between busses 31 and 6, as obtained using the two abovementioned methods, are shown in Figure 3 and Figure 4. As shown in Figure 3 and Figure 4, P E C C S represents the simulation results of the point estimation method based on Gram–Charlier expansion; M C represents the simulation results of the Monte Carlo method. Although the CDARMSE of the line power between busses 31 and 6 was the highest, the simulation results of the two methods were very similar; thus, verifying the accuracy of the point estimation method based on Gram–Charlier expansion.
This paper uses 2 m + 1 point estimation method based on Gram–Charlier expansion to calculate the OIPOPF mathematical model proposed in this paper. In order to obtain the calculation result, 2 m + 1 simulation calculations are required, where m is the number of indeterminate variables. The modified IEEE 39-bus test system has four indeterminate variables, and the modified IEEE 300-bus test system has 28 indeterminate variables.
According to Table 5, the number of simulation calculations of different sizes of test systems is related to the number of uncertain variables.
This paper compares the simulation results of the point estimation method based on Gram–Charlier expansion and Monte Carlo simulation method in the modified IEEE 300-bus test system to verify that the point estimation method based on Gram–Charlier expansion can calculate the OIPOPF model in large test power systems. In the simulation results of the modified IEEE 300-bus test system, the CDARMSE of the line power between busses 123 and 124 was the highest. The PDF and CDF of the line power between busses 123 and 124, as obtained using the two abovementioned methods, are shown in Figure 5 and Figure 6. As shown in Figure 5 and Figure 6, the simulation results of the two methods were very similar. Therefore, it is verified that the point estimation method based on Gram–Charlier expansion can calculate the OIPOPF model in large test power system.

3.4. Discussion

Because the PSHPRPSs contain a high proportion of uncertain renewable energy sources, the prediction error of the day-ahead power is large, and the day-ahead dispatching plan cannot be accurate. Therefore, using the probabilistic AC power flow model not only cannot effectively improve the calculation accuracy, but will also obviously reduce the calculation efficiency. The use of probabilistic DC power flow model can improve the calculation efficiency but reduce the calculation accuracy. In order to verify that the probabilistic power flow model proposed in this paper can not only maintain high calculation accuracy but also maintain high calculation efficiency, this section compares the OIPOPF calculation results with those of the DC power flow model and AC power flow model.
According to the simulation results in Section 3.2, regardless of whether the modified IEEE 39-bus test system or the modified IEEE 300-bus test system was used, the calculation accuracy of OIPOPF model is similar to that of AC power flow model, and higher than that of DC power flow model; the calculation efficiency of OIPOPF model is similar to that of DC power flow model, and higher than that of AC power flow model. To sum up, the OIPOPF model proposed in this paper not only maintains high calculation accuracy, but also maintains high calculation efficiency when setting the day-ahead dispatching plan for the PSHPRPSs.
In this paper, the point estimation method based on Gram–Charlier expansion is used to solve the OIPOPF mathematical model. In order to verify that the point estimation method based on Gram–Charlier expansion can solve the mathematical model of probabilistic optimal power flow proposed in this paper, the simulation results of point estimation method based on Gram–Charlier expansion and Monte Carlo simulation method are compared. According to the simulation results in Section 3.3, regardless of whether the modified IEEE 39-bus small test system or the modified IEEE 300-bus large test system was used, the CDARMSE of the two calculation methods are generally less than 0.001. The figures were drawn by selecting the line power corresponding to the maximum error of the simulation results in the modified IEEE 39-bus small test system and the modified IEEE 300-bus large test system. As shown in the figures in Section 3.3, the simulation results of point estimation method based on Gram–Charlier expansion are similar to those of Monte Carlo simulation method. In summary, by comparing the simulation results of the two algorithms, it is verified that the point estimation method based on Gram–Charlier expansion can accurately solve the OIPOPF mathematical model.
Based on the comparison of the simulation results, the accuracy and efficiency of the proposed OIPOPF model of PSHPRPSs was verified. Moreover, based on the comparison of the calculation results obtained using the two methods, it was included that the point estimation method based on Gram–Charlier expansion can accurately calculate the proposed OIPOPF model.

4. Conclusions

Given that PSHPRPSs contain a large number of uncertain power sources, the day-ahead power prediction error is significantly high. During the dispatch cycle, it is difficult for the day-ahead generation plan to meet the actual active power demand of the power grid. To maintain the power balance of the power grid, the day-ahead power prediction error is accommodated by real-time dispatching and AGC. Due to the significant uncertainty of the power generated by the uncertain power sources, the conventional power flow calculation is not suitable for the SSA of PSHPRPSs. Based on the above reasons, this paper proposes the OIPOPF mathematical model as the power flow model for the SSA of PSHPRPSs. In this paper, the OIPOPF and AC power flow models were applied to a modified IEEE 39-bus test system and modified IEEE 300-bus test system. The simulation results verify that the OIPOPF model can not only ensure the accuracy of power flow results, but also simplify the power flow model. The OIPOPF mathematical model was solved using the point estimation method based on the Gram–Charlier expansion. Compared with the Monte Carlo simulation method, the results indicate that the point estimation method based on Gram–Charlier expansion can be employed to accurately calculate the proposed OIPOPF model. Furthermore, the simulation results verify the accuracy of the OIPOPF model and algorithm.

Author Contributions

Conceptualization, Y.C. and Z.G.; methodology, Y.C.; software, Y.C.; validation, Y.C., and Z.G.; formal analysis, Y.Y.; investigation, H.L.; resources, G.W.; data curation, Y.H.; writing—original draft preparation, Y.C.; writing—review and editing, A.T.T.; visualization, Y.H.; supervision, H.L.; project administration, G.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded in part by Key R&D Projects in the Hebei Province of China (grant number 19212103D), the Innovation capability improvement project in the Hebei Province of China (grant number 19962113D), and the National Natural Science Foundation of China (grant number 51877047).

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

POPFProbabilistic optimal power flow
SSASteady state analysis
PSHPRPSPower systems with high proportions of renewable power sources
AGCAutomatic generation control
OIPOPFOnce-iterative probabilistic optimal power flow
DCDirect current
ACAlternating current
CDFCumulative distribution function
PDFProbability distribution function
CDARMSECumulative distribution average root mean square error

Appendix A

Appendix A presents the simulation results of OIPOPF mathematical model compared with DC power flow model for the modified IEEE 39-bus test system (see Section 3.2 for details). Table A1 presents the percentage errors of the OIPOPF mathematical model and AC power flow simulation results for the modified IEEE 39-bus test system. Table A2 presents the percentage errors of the DC power flow model and AC power flow model simulation results for the modified IEEE 39-bus test system.
Table A1. Percentage errors of the OIPOPF and AC power flow simulation results for the modified IEEE 39-bus test system.
Table A1. Percentage errors of the OIPOPF and AC power flow simulation results for the modified IEEE 39-bus test system.
From BusTo Bus e r r o r From BusTo Bus e r r o r
120.3714150.48
1390.8315160.09
230.8116170.01
2253.0416190.03
2300.3016210.02
343.5916240.14
3182.8817180.59
450.5317274.91
4140.1019200.14
560.2319330.09
580.0620340.01
670.0921220.14
6110.1622230.50
6310.2922350.11
780.1823240.09
891.7523360.09
9392.2525260.06
10110.1425371.10
10130.1426270.48
10320.1726280.02
12111.1126290.03
12131.0128290.08
13140.1629380.03
Table A2. Percentage errors of the DC power flow model and AC power flow model simulation results for the modified IEEE 39-bus test system.
Table A2. Percentage errors of the DC power flow model and AC power flow model simulation results for the modified IEEE 39-bus test system.
From BusTo Bus e r r o r From BusTo Bus e r r o r
122.68141510.30
1396.1115165.63
234.2216170.87
2257.0216191.93
2301.0216211.57
3414.9216245.73
3184.7217180.83
4510.0117272.61
4141.0519201.56
564.1319330.46
580.6220340.50
671.1721220.72
6114.8222233.67
6316.5322353.49
781.8623240.03
8914.5423360.25
93916.93252617.11
10113.9625370.31
10134.0326270.61
10320.1726283.23
121113.5926292.60
121310.4628291.08
13144.3929380.63

Appendix B

Appendix B presents the simulation results of the point estimation method for the modified test power system (see Section 3.3 for details). Table A3 presents the numerical characteristics of the line power simulation results for the modified IEEE 39-bus test system. Table A4 presents the relative error and CDARMSE of the line power for the modified IEEE 39-bus test system.
Table A3. Numerical characteristics of line power simulation results for the modified IEEE 39-bus test system.
Table A3. Numerical characteristics of line power simulation results for the modified IEEE 39-bus test system.
From BusTo BusMonte Carlo MethodPoint Estimation Method
MeanStandard DeviationMeanStandard Deviation
21−1.8310.693.1215.32
139−1.8310.693.1215.32
23−8.7851.42−9.9053.49
252−10.6161.25−6.7864.98
3020000
34−4.2924.44−2.7625.84
1834.4927.137.1428.23
540.202.51−2.245.40
1444.0923.605.0024.57
650.975.780.835.94
580.775.293.086.89
670.815.292.656.38
1161.8410.853.4911.53
316−0.061.4600
780.815.292.656.38
891.5810.575.7313.24
9391.5810.575.7313.24
10111.689.883.1710.49
1013−1.649.91−3.1710.49
32100.042.1800
1112−0.160.97−0.311.04
13120.160.970.311.04
1314−1.8010.87−3.4911.53
1415−5.8934.36−8.4935.88
16155.8934.368.4935.88
161725.92147.0926.29153.49
191610.7793.3015.2496.52
211612.6298.0511.9499.97
24168.4162.917.6063.71
17184.4927.137.1428.23
172721.43120.1619.15125.65
1920−10.7593.30−15.2496.52
33190.010.2700
342010.7593.3015.2496.52
222112.6298.0511.9499.97
22230.3067.021.8066.54
352212.92134.6613.74137.40
23248.4162.917.6063.71
36238.1188.605.8087.00
252611.3261.177.7864.02
37250.719.781.0010.00
2627−21.43120.16−19.15125.65
2826−16.3790.60−13.4694.80
2926−16.3790.60−13.4694.80
2928−16.3790.60−13.4694.80
3829−32.75181.19−26.93189.60
Table A4. Relative error and cumulative distribution average root mean square error (CDARMSE) of line power for the modified IEEE 39-bus test system.
Table A4. Relative error and cumulative distribution average root mean square error (CDARMSE) of line power for the modified IEEE 39-bus test system.
From BusTo Bus ε μ ε σ C D A R M S E
212.700.430.0027
1392.700.430.0027
230.130.040.0002
2520.360.060.0004
3021.710.180.2274
340.360.060.0004
1830.590.040.0005
5411.971.150.0053
1440.230.040.0003
650.140.030.0002
583.000.300.0025
672.280.210.0020
1160.890.060.0008
3161.001.000.4965
782.280.210.0020
892.630.250.0022
9392.630.250.0022
10110.890.060.0008
10130.940.060.0009
32101.001.000.4957
11120.910.070.0009
13120.910.070.0009
13140.940.060.0009
14150.440.040.0004
16150.440.040.0004
16170.010.040.0002
19160.420.030.0003
21160.050.020.0001
24160.100.010.0001
17180.590.040.0005
17270.110.050.0002
19200.420.030.0003
33191.001.000.4939
34200.420.030.0003
22210.051.000.0001
22234.940.030.0001
35220.060.020.0001
23240.100.010.0001
36230.280.020.0002
25260.310.050.0004
37250.420.020.0002
26270.110.050.0002
28260.180.050.0002
29260.180.050.0002
29280.180.050.0002
38290.180.050.0002

References

  1. Semlyen, A.; Leon, F.D. Quasi-Newton power flow using partial Jacobian updates. IEEE Trans. Power Syst. 2001, 16, 332–339. [Google Scholar] [CrossRef]
  2. Vanamerongen, R.A.M. A general-purpose version of the fast decoupled loadflow. IEEE Trans. Power Syst. 1989, 4, 760–770. [Google Scholar] [CrossRef]
  3. Borkowska, B. Probabilistic Load Flow. IEEE Trans. Power Appar. Syst. 1974, 3, 752–759. [Google Scholar] [CrossRef]
  4. Constante-Flores, G.E.; Illindala, M.S. Data-Driven Probabilistic Power Flow Analysis for a Distribution System with Renewable Energy Sources Using Monte Carlo Simulation. IEEE Trans. Ind. Appl. 2019, 55, 174–181. [Google Scholar] [CrossRef]
  5. Peng, S.; Chen, H.; Lin, Y.; Shu, T.; Lin, X.; Tang, J.; Li, W.; Wu, W. Probabilistic Power Flow for Hybrid AC/DC Grids with Ninth-Order Polynomial Normal Transformation and Inherited Latin Hypercube Sampling. Energies 2019, 12, 3088. [Google Scholar] [CrossRef] [Green Version]
  6. Xie, Z.Q.; Ji, T.Y.; Li, M.S.; Wu, Q.H. Quasi-Monte Carlo Based Probabilistic Optimal Power Flow Considering the Correlation of Wind Speeds Using Copula Function. IEEE Trans. Power Syst. 2018, 33, 2239–2247. [Google Scholar] [CrossRef]
  7. Zhu, J.; Zhang, Y.; Chen, H. Probabilistic Load Flow with Wind Farms Using a Frequency and Duration Method. IEEE Access 2018, 6, 74812–74827. [Google Scholar] [CrossRef]
  8. Allan, R.N.; Al-Shakarchi, M.R.G. Probabilistic a.c. load flow. Proc. Inst. Electr. Eng. 1976, 123, 531–536. [Google Scholar] [CrossRef]
  9. Allan, R.N.; Silva, A.M.L.D. Probabilistic load flow using multilinearisations. IEE Proc. C Gener. Transm. Distrib. 1981, 128, 280–287. [Google Scholar] [CrossRef]
  10. Deng, X.; He, J.; Zhang, P. A Novel Probabilistic Optimal Power Flow Method to Handle Large Fluctuations of Stochastic Variables. Energies 2017, 10, 1623. [Google Scholar] [CrossRef] [Green Version]
  11. Huynh, V.K.; Ngo, V.D.; Le, D.D.; Nguyen, N.T.A. Probabilistic Power Flow Methodology for Large-Scale Power Systems Incorporating Renewable Energy Sources. Energies 2018, 11, 2624. [Google Scholar] [CrossRef] [Green Version]
  12. Lyu, J.-K.; Heo, J.-H.; Park, J.-K.; Kang, Y.-C. Probabilistic Approach to Optimizing Active and Reactive Power Flow in Wind Farms Considering Wake Effects. Energies 2013, 6, 5717–5737. [Google Scholar] [CrossRef] [Green Version]
  13. Sobierajski, M. Method of stochastic load flow calculation. Arch. Elektrotech. 1978, 60, 37–40. [Google Scholar] [CrossRef]
  14. Zhang, X.; Guo, Z.; Chen, W. Probabilistic Power Flow Method Considering Continuous and Discrete Variables. Energies 2017, 10, 590. [Google Scholar] [CrossRef] [Green Version]
  15. Shi, D.; Cai, D.; Chen, J.; Duan, X.; Li, H.; Yao, M. Probabilistic load flow calculation based on cumulant method considering correlation between input variables. Proc. Chin. Soc. Electr. Eng. 2012, 32, 104–113. [Google Scholar]
  16. Saunders, C.S. Point Estimate Method Addressing Correlated Wind Power for Probabilistic Optimal Power Flow. IEEE Trans. Power Syst. 2014, 29, 1045–1054. [Google Scholar] [CrossRef]
  17. Verbic, G.; Canizares, C.A. Probabilistic Optimal Power Flow in Electricity Markets Based on a Two-Point Estimate Method. IEEE Trans. Power Syst. 2006, 21, 1883–1893. [Google Scholar] [CrossRef]
  18. Li, J.; Zhang, B.; Liu, Y. Data Mining in Nonlinear Probabilistic Load Flow Based on Monte Carlo Simulation. In Proceedings of the 1st International Conference on Information Science and Engineering, ICISE2009, Nanjing, China, 26–28 December 2009; IEEE Computer Society: Nanjing, China, 2009; pp. 833–836. [Google Scholar]
  19. Baghaee, H.R.; Mirsalim, M.; Gharehpetian, G.B.; Talebi, H.A. Fuzzy unscented transform for uncertainty quantification of correlated wind/PV microgrids: Possibilistic–probabilistic power flow based on RBFNNs. IET Renew. Power Gener. 2017, 11, 867–877. [Google Scholar] [CrossRef]
  20. Baghaee, H.R.; Mirsalim, M.; Gharehpetian, G.B.; Talebi, H.A. Three Phase AC/DC Power-flow for Balanced/Unbalanced Microgrids Including Wind/Solar, Droop-Controlled and Electronically-Coupled Distributed Energy Resources using RBF Neural Networks. IET Power Electron. 2016, 10, 313–328. [Google Scholar] [CrossRef]
  21. Baghaee, H.R.; Mirsalim, M.; Gharehpetian, G.B.; Talebi, H.A. Generalized three phase robust load-flow for radial and meshed power systems with and without uncertainty in energy resources using dynamic radial basis functions neural networks. J. Clean. Prod. 2018, 174, 96–113. [Google Scholar] [CrossRef]
  22. Pan, W.; Liu, W.-Y.; Yang, Y.-H. Point estimate method for probabilistically optimal power flow computation. Proc. Chin. Soc. Electr. Eng. 2008, 28, 28–33. [Google Scholar]
  23. Wang, R.; Gu, W.; Sun, R.; Li, Q. Analysis on wind power penetration limit based on probabilistically optimal power flow. Power Syst. Technol. 2011, 35, 214–220. [Google Scholar]
  24. Xue, M.-J.; Yang, X.-P.; Ma, X.-Y. The united optimal operation of wind power and cascaded hydropower based on the optimal power flow theory. J. Hydraul. Eng. 2014, 45, 368–375. [Google Scholar]
  25. Li, X.; Li, Y.; Zhang, S. Analysis of Probabilistic Optimal Power Flow Taking Account of the Variation of Load Power. IEEE Trans. Power Syst. 2008, 23, 992–999. [Google Scholar]
  26. Schellenberg, A.; Rosehart, W.; Aguado, J. Cumulant-based probabilistic optimal power flow (P-OPF) with Gaussian and gamma distributions. IEEE Trans. Power Syst. 2005, 20, 773–781. [Google Scholar] [CrossRef] [Green Version]
  27. Tamtum, A.; Schellenberg, A.; Rosehart, W.D. Enhancements to the Cumulant Method for Probabilistic Optimal Power Flow Studies. IEEE Trans. Power Syst. 2009, 24, 1739–1746. [Google Scholar] [CrossRef]
  28. Zou, B.; Xiao, Q. Solving Probabilistic Optimal Power Flow Problem Using Quasi Monte Carlo Method and Ninth-Order Polynomial Normal Transformation. IEEE Trans. Power Syst. 2014, 29, 300–306. [Google Scholar] [CrossRef]
  29. Ke, D.; Chung, C.Y.; Sun, Y. A Novel Probabilistic Optimal Power Flow Model with Uncertain Wind Power Generation Described by Customized Gaussian Mixture Model. IEEE Trans. Sustain. Energy 2016, 7, 200–212. [Google Scholar] [CrossRef]
  30. Li, Y.Z.; Wu, Q.H. Downside Risk Constrained Probabilistic Optimal Power Flow with Wind Power Integrated. IEEE Trans. Power Syst. 2016, 31, 1649–1650. [Google Scholar] [CrossRef]
  31. Sun, W.; Zamani, M.; Zhang, H.; Li, Y. Probabilistic Optimal Power Flow with Correlated Wind Power Uncertainty via Markov Chain Quasi-Monte-Carlo Sampling. IEEE Trans. Ind. Inform. 2019, 15, 6058–6069. [Google Scholar] [CrossRef]
  32. Banos, R.; Manzano-Agugliaro, F.; Montoya, F.G.; Gil, C.; Alcayde, A.; Gomez, J. Optimization methods applied to renewable and sustainable energy: A review. Renew. Sustain. Energy Rev. 2011, 15, 1753–1766. [Google Scholar] [CrossRef]
  33. Kneevi, G.; Topi, D.; Juri, M.; Nikolovski, S. Joint market bid of a hydroelectric system and wind parks. Comput. Electr. Eng. 2019, 74, 138–148. [Google Scholar] [CrossRef]
  34. Liu, C.; Chao, Q.; Wei, L. Wind-storage coupling based on actual data and fuzzy control in multiple time scales for real-time rolling smoothing of fluctuation. Electr. Power Autom. Equip. 2015, 35, 35–41. [Google Scholar]
  35. Wang, K.; Zhang, B.; Yan, D.; Li, Y.; Luo, T. A multi-time scale rolling coordination scheduling method for power grid integrated with large scale wind farm. Power Syst. Technol. 2014, 38, 2434–2440. [Google Scholar]
  36. Wang, Z.; Guo, Z. On critical timescale of real-time power balancing in power systems with intermittent power sources. Electr. Power Syst. Res. 2018, 155, 246–253. [Google Scholar] [CrossRef]
  37. Zhai, J.; Ren, J.; Zhou, M.; Li, Z. Multi-time scale fuzzy chance constrained dynamic economic dispatch model for power system with wind power. Power Syst. Technol. 2016, 40, 1094–1099. [Google Scholar]
  38. Zhang, B.; Wu, W.; Zheng, T.; Sun, H. Design of a multi-time scale coordinated active power dispatching system for accommodating large scale wind power penetration. Autom. Electr. Power Syst. 2011, 35, 1–6. [Google Scholar]
  39. Pei, Z.; Lee, S.T. Probabilistic load flow computation using the method of combined cumulants and Gram-Charlier expansion. IEEE Trans. Power Syst. 2004, 19, 676–682. [Google Scholar]
  40. Liu, D. About the Researches on Numerical Characteristics of Geometric Distribution. Master’s Thesis, Northeastern University, Shenyang, China, December 2007. [Google Scholar]
  41. Morales, J.M.; Perez-Ruiz, J. Point Estimate Schemes to Solve the Probabilistic Power Flow. IEEE Trans. Power Syst. 2007, 22, 1594–1601. [Google Scholar]
  42. Pender, J. Gram charlier expansion for time varying multiserver queues with abandonment. SIAM J. Appl. Math. 2014, 74, 1238–1265. [Google Scholar]
  43. Thomas, A. Wind Power in Power System; John Wiley & Sons Ltd.: London, UK, 2012; pp. 753–769. [Google Scholar]
Figure 1. Flow chart for the calculation of the once-iterative probabilistic optimal power flow (OIPOPF) based on the Gram–Charlier expansion point estimation method. CDF—cumulative distribution function, PDF—probability distribution function.
Figure 1. Flow chart for the calculation of the once-iterative probabilistic optimal power flow (OIPOPF) based on the Gram–Charlier expansion point estimation method. CDF—cumulative distribution function, PDF—probability distribution function.
Sustainability 12 00518 g001
Figure 2. Modified Institute of Electrical and Electronic Engineers (IEEE) 39-bus test system diagram.
Figure 2. Modified Institute of Electrical and Electronic Engineers (IEEE) 39-bus test system diagram.
Sustainability 12 00518 g002
Figure 3. PDF of line with the largest cumulative distribution average root mean square error (CDARMSE) for the modified IEEE 39-bus test system.
Figure 3. PDF of line with the largest cumulative distribution average root mean square error (CDARMSE) for the modified IEEE 39-bus test system.
Sustainability 12 00518 g003
Figure 4. CDF of line with the largest CDARMSE for the modified IEEE 39-bus test system.
Figure 4. CDF of line with the largest CDARMSE for the modified IEEE 39-bus test system.
Sustainability 12 00518 g004
Figure 5. PDF of line with the largest CDARMSE for the modified IEEE 300-bus test system.
Figure 5. PDF of line with the largest CDARMSE for the modified IEEE 300-bus test system.
Sustainability 12 00518 g005
Figure 6. CDF of line with the largest CDARMSE for the modified IEEE 300-bus test system.
Figure 6. CDF of line with the largest CDARMSE for the modified IEEE 300-bus test system.
Sustainability 12 00518 g006
Table 1. Numerical characteristics of uncertain renewable energy sources prediction errors of the modified IEEE 39-bus test system.
Table 1. Numerical characteristics of uncertain renewable energy sources prediction errors of the modified IEEE 39-bus test system.
BusMeanStandard Deviation
340.15240.9652
350.13741.374
360.0580.87
370.010.1
Table 2. Numerical characteristics of uncertain renewable energy sources prediction errors of the modified IEEE 300-bus test system.
Table 2. Numerical characteristics of uncertain renewable energy sources prediction errors of the modified IEEE 300-bus test system.
BusMeanStandard DeviationBusMeanStandard Deviation
840.0480.4752130.0370.372
910.0260.2552200.020.2
920.0390.392210.0550.55
980.0170.1682220.0350.35
1080.0220.2172270.040.403
1240.0340.3470010.0570.567
1410.0380.38170020.0720.723
1430.080.79670110.0330.334
1460.0180.18470120.0470.472
1470.0320.31770170.0430.43
1490.020.20370230.0290.285
1520.0470.47270240.0510.51
1530.0320.31670710.0220.216
1980.0520.52490540.0150.15
Table 3. Average errors of simulation results for the modified IEEE 39-bus test system and the modified IEEE 300-bus test system.
Table 3. Average errors of simulation results for the modified IEEE 39-bus test system and the modified IEEE 300-bus test system.
Scale of Test Power SystemOnce-Iterative Probabilistic Optimal Power Flow (OIPOPF) ModelDirect Current (DC) Power Flow Model
IEEE 39-bus test system0.64%4.45%
IEEE 300-bus test system0.90%11.91%
Table 4. Simulation operation time of the modified IEEE 39-bus test system and the modified IEEE 300-bus test system.
Table 4. Simulation operation time of the modified IEEE 39-bus test system and the modified IEEE 300-bus test system.
Scale of Test Power SystemOIPOPFDC Power FlowAC Power Flow
Operation time of IEEE 39-bus test system (s)0.0120.0090.15
Operation time of IEEE 300-bus test system (s)0.0370.0310.55
Table 5. Number of once-iterative probabilistic optimal power flow OIPOPF mathematical model simulation calculations.
Table 5. Number of once-iterative probabilistic optimal power flow OIPOPF mathematical model simulation calculations.
Scale of Test Power SystemVariablesPECCS Simulation TimesMC Simulation Times
Operation time of IEEE 39-bus test system491000
Operation time of IEEE 300-bus test system28571000

Share and Cite

MDPI and ACS Style

Chen, Y.; Guo, Z.; Li, H.; Yang, Y.; Tadie, A.T.; Wang, G.; Hou, Y. Probabilistic Optimal Power Flow for Day-Ahead Dispatching of Power Systems with High-Proportion Renewable Power Sources. Sustainability 2020, 12, 518. https://doi.org/10.3390/su12020518

AMA Style

Chen Y, Guo Z, Li H, Yang Y, Tadie AT, Wang G, Hou Y. Probabilistic Optimal Power Flow for Day-Ahead Dispatching of Power Systems with High-Proportion Renewable Power Sources. Sustainability. 2020; 12(2):518. https://doi.org/10.3390/su12020518

Chicago/Turabian Style

Chen, Yue, Zhizhong Guo, Hongbo Li, Yi Yang, Abebe Tilahun Tadie, Guizhong Wang, and Yingwei Hou. 2020. "Probabilistic Optimal Power Flow for Day-Ahead Dispatching of Power Systems with High-Proportion Renewable Power Sources" Sustainability 12, no. 2: 518. https://doi.org/10.3390/su12020518

APA Style

Chen, Y., Guo, Z., Li, H., Yang, Y., Tadie, A. T., Wang, G., & Hou, Y. (2020). Probabilistic Optimal Power Flow for Day-Ahead Dispatching of Power Systems with High-Proportion Renewable Power Sources. Sustainability, 12(2), 518. https://doi.org/10.3390/su12020518

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop