Next Article in Journal
Design Optimization and Coupled Dynamics Analysis of an Offshore Wind Turbine with a Single Swivel Connected Tether
Next Article in Special Issue
Modeling and Parameter Optimization of Grid-Connected Photovoltaic Systems Considering the Low Voltage Ride-through Control
Previous Article in Journal
Performance Assessment of Front-Mounted Beet Topper Machine for Biomass Harvesting
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Probabilistic Power Flow Algorithm Based on Principal Component Analysis and High-Dimensional Model Representation Techniques

State Key Laboratory of Advanced Electromagnetic Engineering and Technology, Huazhong University of Science and Technology, 1037 Luoyu Road, Wuhan 430074, China
*
Author to whom correspondence should be addressed.
Energies 2020, 13(14), 3520; https://doi.org/10.3390/en13143520
Submission received: 6 June 2020 / Revised: 30 June 2020 / Accepted: 7 July 2020 / Published: 8 July 2020
(This article belongs to the Special Issue Electric Power Systems Research 2020)

Abstract

:
Because the penetration level of renewable energy sources has increased rapidly in recent years, uncertainty in power system operation is gradually increasing. As an efficient tool for power system analysis under uncertainty, probabilistic power flow (PPF) is becoming increasingly important. The point-estimate method (PEM) is a well-known PPF algorithm. However, two significant defects limit the practical use of this method. One is that the PEM struggles to estimate high-order moments accurately; this defect makes it difficult for the PEM to describe the distribution of non-Gaussian output random variables (ORVs). The other is that the calculation burden is strongly related to the scale of input random variables (IRVs), which makes the PEM difficult to use in large-scale power systems. A novel approach based on principal component analysis (PCA) and high-dimensional model representation (HDMR) is proposed here to overcome the defects of the traditional PEM. PCA is applied to decrease the dimension scale of IRVs and eliminate correlations. HDMR is applied to estimate the moments of ORVs. Because HDMR considers the cooperative effects of IRVs, it has a significantly smaller estimation error for high-order moments in particular. Case studies show that the proposed method can achieve a better performance in terms of accuracy and efficiency than traditional PEM.

1. Introduction

In recent decades, the penetration level of renewable energy sources (RESs), such as wind and solar power, has increased rapidly. RESs have undispatchable and limited prediction characteristics, which brings challenges for power system operation. For a power system with large-scale RES integration, power generation becomes uncertain, which causes system state variables, such as line flows and bus voltage, to vary. Clearly, in this situation, the state variables should not be seen as having deterministic values, and the traditional deterministic power flow (DPF) is incapable of analyzing system states.
Instead of DPF, the concept of probabilistic power flow (PPF) emerges [1]. In PPF, some power injections, including RES generation and loads, are seen as input random variables (IRVs), while state variables are seen as output random variables (ORVs). Given the probability information of IRVs, the moments or probability distribution of state variables can be obtained by a PPF algorithm. PPF is an effective tool that can provide valuable information for applications, such as operational security assessment [2], planning and economic dispatch [3]. Therefore, the study of the PPF algorithm has seen increased attention.
Since the concept of PPF emerged, a large number of algorithms have been proposed. These algorithms can be roughly classified into three main categories [4]: sampling methods, analytical methods, and approximate methods. The goal of a PPF algorithm is to describe the uncertainty characteristics of state variables precisely and efficiently. However, there are several challenges for an algorithm to achieve this goal, including 1. scales of IRVs—for a large power system, the number of IRVs is high, which may cause the curse of dimensionality problem and lead to an unacceptable calculation burden [5]; 2. the non-Gaussian distribution of IRVs, which means that estimating only mean and variance is not sufficient to describe the shape of the ORVs’ probability distribution, and a PPF algorithm should have the ability to estimate ORVs’ high-order moments accurately; 3. the correlation between IRVs, if not considered, decreases the accuracy of the results.
Monte Carlo simulation (MCS) [6] and its improved variants, such as Latin hypercube sampling [7] and quasi-MCS [8], are representatives of sampling methods. This kind of method requires a large number of samples to be selected from the probability distribution function (PDF) of an IRV and obtains the corresponding number of DPFs. The more samples selected, the more accurate the result is. The mechanism of sampling methods ensures accuracy and is widely regarded as the most accurate method. However, the whole DPF process is quite time-consuming. Therefore, sampling methods are usually a reference method to judge other methods’ accuracy levels.
The conventional method [9] and cumulant method [10,11,12,13,14] are typical representatives of analytical methods. The conventional method for large power systems requires a large number of convolution operations, which is time-consuming. The cumulant method only needs a single DPF, and the computational time is significantly decreased; therefore, this method has received more attention than the conventional method. When using the cumulant method, several assumptions are made, such as the linearization of power flow equations and the independence of IRVs. However, although these assumptions make the algorithm efficient, in cases in which the fluctuation range of IRVs is high, the linearization assumption is violated, which will cause unacceptable calculation error [14].
The Point estimate method [15,16,17,18,19,20,21,22] is a typical representative of approximate methods, which balance accuracy and efficiency well. To the best of the authors’ knowledge, most PEMs are based on Hong’s work [17], which means that 2m+1 DPF calculations are performed when m IRVs are considered. The Nataf transformation is combined with Hong’s PEM to handle correlation issues [20]. Zhao’s PEM scheme [16,21] also called univariate dimension reduction [22], is proposed, which shows better accuracy than Hong’s scheme. Although much work has been performed to improve the performance of PEMs, the existing methods still have two major drawbacks. The first is that the calculation burden is directly related to the number of IRVs, which means that for a system consisting of a large number of IRVs, the computational efficiency will be dragged down, which makes PEMs difficult to apply in practical cases. The second is that it is difficult for the existing methods to estimate high-order moments accurately, such as skewness and kurtosis, which means that existing methods are inappropriate to apply in a power system that has non-Gaussian IRVs.
Evidently, if the drawbacks of the PEM can be eliminated effectively, as an algorithm that can balance the calculation efficiency and accuracy, then it has practical application value. The motivation of this paper is to overcome the drawbacks of existing PEMs; therefore, we propose a novel probabilistic power flow algorithm based on principal component analysis (PCA) and high-dimensional model representation (HDMR). The main contributions of this paper are as follows:
(1) HDMR is applied to obtain the ORVs’ moments in this paper. Unlike existing PEM schemes, HDMR considers the cooperative effects of IRVs on ORVs. Case studies show that the estimation error of high-order moments is significantly reduced compared to that of other PEM schemes;
(2) PCA is applied to handle the correlation of IRVs and improve the efficiency of HDMR. Through PCA, the principal components (PCs) of original IRVs are selected, which contain the salient information of the original IRVs but are much less complex and uncorrelated. By using the PCs instead of the original IRVs, curse of dimensionality problem is avoided, which greatly improves computational efficiency, and correlation is well handled.
The remainder of this paper is organized as follows: Section 2 presents the methodology of the proposed method, followed by the procedure in Section 3. In Section 4, case studies are examined. The proposed method is compared with MCS and Zhao’s PEM. The influence of PCA and the adaptability of the proposed method are also analyzed in Section 4. In Section 5, a discussion of the proposed method is presented. Finally, conclusions are drawn in Section 6.

2. Methodology

2.1. Problem Formation

The load flow equation shows the relationship between IRVs and ORVs, and the Alternating Current (AC) load flow equations are as follows
{ P i = V i j i V j ( G i j cos θ i j + B i j sin θ i j ) Q i = V i j i V j ( G i j sin θ i j B i j cos θ i j )
where P i and Q i are the active and reactive power injections at bus i, respectively, V i and V j are the voltage magnitudes of bus i and j, respectively, θ i j is the voltage angle difference between bus i and j, and G i j and B i j are parameters related to the line that connects bus i and j.
Clearly, the above AC load flow equations can be expressed with a qualitative equation
Y = G ( X )
where X is the vector of IRVs, including active and reactive power injections. Y is the vector of ORVs, including bus voltage magnitude V and angle θ and line active flows P and reactive flows Q, and G ( · ) represents the nonlinear AC power flow equations.
Given the probability information of the IRVs, we want to know the moments of the ORVs. The lth central moment of the ORVs is as follows
Ε ( Y l ) = E [ G l ( X ) ] = G l ( X ) f X ( X ) d x
where E ( · ) is the expectation operator, l denotes the lth central moment, and fX(X) represents the joint probability density function (PDF) of   X .
Here are some challenges for directly solving (3). The first is that the integration expressed in (3) is multi-dimensional; secondly, if X follows non-Gaussian probability distributions, the expression of fX(X) is complicated; moreover, the difficulty of integration is related to the non-linearity of G l ( X ) . Proper transformation and dimensionality reduction are needed to perform this integration. In this paper, PCA and HDMR are combined to perform the integration.

2.2. Principal Component Analysis

PCA [23,24] is an efficient tool for data analysis, especially for high-dimensional data with potential correlations. The mechanism of PCA is orthogonal transformation. Through orthogonal transformation, a high-dimensional dataset is transformed to a low-dimensional dataset while maintaining its main characteristics. In this section, we introduce the procedure of PCA, as it can eliminates the correlations among IRVs, and improve the overall efficiency of proposed algorithm.
In this paper, the power injections of each bus are considered IRVs, specifically, the power injections of RES and load. The distribution model of IRVs can significantly influence the results of PPF, especially the distribution of RESs. However, distribution modelling itself is not the focus of this paper; moreover, the proposed method has good adaptability to any distribution model, thus, samples of IRVs can be derived from commonly used distributions or from historical data.
Suppose m power injection samples are available for each IRV and the dimension of the IRV vector is n; then, an n × m IRV sample matrix X can be constructed:
X = [ x 1 1 x 1 2 x 1 m x 2 1 x 2 2 x 2 m x n 1 x n 2 x n m ]
The first step is centralization. By subtracting the mean value of each IRV μ X , the centred sample matrix XC can be obtained
X C = X μ X
Suppose that the IRVs have linear correlations, and the covariance matrix Σ corresponding to XC is used to describe the correlation
Σ = [ σ 1 2 σ 1 , 2 2 σ 1 , n 2 σ 2 , 1 2 σ 2 2 σ 2 , n 2 σ n , 1 2 σ n , 2 2 σ n , n 2 ]
where σ i 2 is the variance of random variable i and σ i , j 2 is the covariance between random variables i and j. σ i , j 2 is calculated by
σ i , j 2 = 1 m k = 1 m ( x i k μ i ) ( x j k μ j )
where x i k and x j k are the kth samples of IRVs i and j, respectively, and μ i and μ j are the means of IRVs i and j, respectively.
Clearly, the covariance matrix is an n × n symmetric positive semi-definite matrix, and its eigenvalues can be calculated by
det ( Σ λ i I ) = 0     i = 1 , 2 , , n
where det(.) is the determinant, λ i is an eigenvalue, and I is the identity matrix.
Sorting the eigenvalues in descending order, we can calculate the corresponding eigenvectors e i by
Σ e i = λ i e i   i = 1 , 2 , , n
Thus, we can construct the projection matrix E
E = [ e 1   e 2   e n ]
Transforming the dataset XC via E, the PC matrix Z is obtained
Z = E T X c
Z is an n × m matrix that has the same dimension as X or XC. The ith row of the matrix Z is the ith PC, and we use z i to express it. z i is a linear combination of IRVs that has no actual physical meaning.
Here are some characteristics of PCs:
(1) Under the assumption that only linear correlation exists, the covariance between any two PCs is 0; therefore, correlation is eliminated;
(2) The variance of z i is λ i . Clearly, the variance in elements in Z has a descending order. According to PCA, the larger the variance is, the more information is contained in a PC. Therefore, in Z, the importance of z i follows a descending order.
Clearly, matrix Z and matrix X have the same dimension: if all of the components in Z are retained, the scale of the PPF problem is not reduced. Thus, we must determine how many PCs should be retained, which can be achieved by calculating the PCs’ cumulative contributions.
The contribution α i of ith PC is the proportion of its corresponding eigenvalue compared to all eigenvalues. The greater the proportion, the stronger the ability of the PC to integrate the information of the original IRVs.   α i is calculated as follows
α i = λ i λ 1 + λ 2 + λ n
The cumulative contribution of the first k components is calculated as
M k = i = 1 k α i
The greater the value of M k , the more fully the first k PCs represent the information of the original IRVs. When choosing the retained PCs, M k should be compared to a threshold, such as M k     95 % . Suppose the first k components of the projection matrix are chosen as E k = [ e 1 , e 2 , e k ] , then the retained PC matrix is Z k = E k T X c . Z k is a k × m matrix, and it retains the most information of the original IRVs.
For a PEM-like algorithm, the calculation burden is closely related to the dimensions of the IRVs. By performing PCA on the IRV vector X, we can obtain an uncorrelated k × m matrix Z k , where k<<n. Clearly, the calculation burden is not directly related to the dimension of IRVs but relates to the dimension of the retained PC matrix. Using PCA before the moment estimation procedure, the total efficiency of the algorithm is significantly improved; moreover, since there is no correlation between PCs, the impact of the correlations between IRVs is handled well.
A PC is a linear combination of IRVs, and its distribution is arbitrary. Sometimes, the abscissas cannot be directly generated from an arbitrary probability distribution, as it may be outside the region in which the random variable is defined [16]. Here, we perform a Nataf transformation [25] on PCs. This transformation will connect the principal component space to a standard Gaussian space, which will facilitate the subsequent implementation of HDMR and ensure the validity of chosen abscissas.
Assuming that the marginal cumulative distribution function (CDF) Fi(zi), i = 1…k of each PC zi is available, we can transform zi into a standard normal variable ui
{ Φ ( u i ) = F i ( z i ) u i = Φ 1 ( F i ( z i ) ) i = 1 , 2 , , k
where Φ ( ) and Φ - 1 ( ) are the CDF and inverse CDF of the standard normal variable u i .
Therefore, through PCA and marginal transformation, Y can be expressed as a function of the independent standard Gaussian vector U:
Y = G ( X ) = G ( E T Z + μ X ) = G ( E T F 1 ( Φ ( U ) ) + μ X ) = H ( U )

2.3. High-Dimensional Model Representation

HDMR [26,27,28] is an analysis tool for capturing the high-dimensional relationships between IRVs and ORVs. The core idea of HDMR is using several low-dimensional functions to approximate the original high-dimensional functions. Applying HDMR to (14), the function H ( U ) is split into a summation of several sub-function terms
Y = H ( U ) = H 0 + i H i ( u i ) + i 1 < i 2 H i 1 i 2 ( u i 1 , u i 2 ) + + i 1 < i 2 < i l H i 1 i 2 i l ( u i 1 , u i 2 , u i l ) + + H i 1 i 2 i n ( u i 1 , u i 2 , u i n )
H0 is a constant term that represents the mean response on Y. H i ( u i ) is a first-order term that represents the effects of a single IRV u i on Y;   H i 1 , i 2 ( u i 1 , u i 2 ) is a second-order term that represents the cooperative effects of the IRVs u i 1 , u i 2 on Y. Clearly, higher-order terms represent the cooperative effects of multiple IRVs on Y, and H i 1 i 2 i n ( u i 1 , u i 2 , u i n ) represents the cooperative effects of all IRVs u i 1 , u i 2 , u i n on Y.
The detailed expressions of the components in each term are as follows
H 0 = H ( U c ) H i ( u i ) = H ( u i , U i , c ) H 0 H i 1 i 2 ( u i 1 , u i 2 ) = H ( u i 1 , u i 2 , U i 1 i 2 , c ) H i 1 ( u i 1 ) H i 2 ( u i 2 ) H 0 H i 1 i l ( u i 1 , , u i l ) = H i 1 i 2 i l ( u i 1 , u i 2 , u i l , U i 1 i 2 i l , c ) { j 1 j 2 j l 1 } { i 1 i 2 i l } H j 1 j 2 j l 1 ( u j 1 , u j 2 , u j l 1 ) { j 1 j 2 j l 2 } { i 1 i 2 i l } H j 1 j 2 j l 2 ( u j 1 , u j 2 , u j l 2 ) j { i 1 i 2 i l } H j ( u j ) H 0
where U c is the reference vector; as we have transformed the original IRV space to the standard Gaussian space by PCA and marginal transformation, U c is [ 0 , , 0 ] here. U i , c is formulated after eliminating the sub-vectors of U c corresponding to u i ; therefore, H ( u i , U i , c ) = H ( 0 , , 0 , u i , 0 , , 0 ) . The expression of the higher-order component can be obtained recursively.
The expression and calculation process become complex after including the terms in (15). A practical method for doing this is to use a truncated form of (16); given a highest order for term s, only cooperative effects with no more than s variables are involved.
Considering terms up to first order, and combining (16) and (17), H ( U ) is expressed as
H ( U ) = k = 1 n H ( u k , U k , c ) ( n 1 ) H 0 + 2
Similarly, considering terms up to second order, H ( U ) is expressed as
H ( U ) = i 1 < i 2 H ( u i 1 , u i 2 , U i j , c ) ( n 2 ) k = 1 n H ( u k , U k , c ) + ( n 1 ) ( n 2 ) 2 H 0 + 3
where 2 and 3 are the residue errors. Ignoring 2 and 3 , we obtain the first-order and second-order approximations of H ( U ) . Clearly, the first-order approximation only considers a single variable’s individual effect on Y, and all cooperative effects of IRVs are ignored, while the second-order approximation considers the cooperative effects of any two IRVs. It is worth mentioning that the first-order approximation is exactly the same as that used in Zhao’s PEM scheme [15]; therefore, Zhao’s PEM scheme could be considered as a special case of HDMR.
Since load flow equations are non-linear equations, ignoring the cooperative effects of IRVs will inevitably cause estimation errors; therefore, it is necessary to use a form of approximation that considers cooperative effects. According to our experience, ignoring the cooperative effects of IRVs will cause an unacceptable level of estimation error for high-order moments, which will seriously affect the accurate description of the ORVs’ distribution, especially when the IRVs follow a non-Gaussian distribution: considering the cooperative effects of two IRVs can significantly eliminate this error and is sufficient in practical use.
Using the second-order approximation, the lth central moment of Y can be calculated by
E ( Y l ) = i 1 < i 2 E [ H l ( u i 1 , u i 2 , U i 1 i 2 , c ) ] ( n 2 ) k = 1 n E [ H l ( u k , U k , c ) ] + ( n 1 ) ( n 2 ) 2 H 0 l
In (19), the expectation of the univariate and bivariate functions can be calculated by Gauss–Hermite quadrature
E [ H l ( u k , U k , c ) ] = m = 1 r w G H , m π H l ( 2 α G H , m , U k , c ) E [ H l ( u i 1 , u i 2 , U i 1 i 2 , c ) ] = m 1 = 1 r m 2 = 1 r w G H , m 1 w G H , m 2 π H l ( 2 α G H , m 1 , 2 α G H , m 2 , U i 1 i 2 , c )
where r is the number of abscissas, w GH , m is the weight and α GH , m is the abscissa. Table 1 shows the typical abscissas and their corresponding weights.

3. Solution Procedure

In this section, the solution procedure of the proposed method is summarized step by step as follows:
1. Prepare the necessary data, including samples of uncertainty sources (RES and load), which could be generated from distribution models or collected from historical data and the parameters of the power system;
2. Obtain the centralized sample matrix Xc and its corresponding covariance matrix Σ ;
3. Obtain the projection matrix according to (8), (9) and (10);
4. Obtain the retained PC matrix Z k according to (11), (12), and (13) and the predetermined threshold value;
5. According to the dimension of Z k , choose abscissas in standard Gaussian space, then use (11) and (14) to project these abscissas into the original correlated IRV space;
6. Use the projected abscissas as IRVs and calculate the corresponding DPFs;
7. Calculate the moments of the ORVs from (20) and (21).

4. Case Study

The proposed method is tested on two modified test systems. The DLF program is based on MATPOWER 6.0 [29] and is executed on a PC with a 2.8-GHz CPU and 16 GB RAM.
The third-order moment and fourth-order moment, namely skewness and kurtosis, are important parameters to describe the shape of a PDF, we calculate the ORV’s moments up to the fourth-order moment. The performance of the proposed method is compared with that of HDMR alone, Zhao’s PEM proposed in [21], and 50000 iterations of MCS. Additionally, three contribution thresholds of PCA ( M k     99 % , M k     95 % and M k     90 % respectively) are set subjectively to measure the influence of the contribution threshold on the accuracy and efficiency of the proposed method. The number of abscissas of HDMR is chosen as 5, and, for consistency, that of Zhao’s PEM is adjusted to be the same.
The average relative error index (AREI) [11] is adopted to indicate the accuracy of the moment of the ORVs estimated using each method. The average relative error index is defined as follows
ε ¯ i = 1 N j N | γ i , j E V - γ i , j M C γ i , j M C | × 100 %
where i corresponds to the ith central moment; N is the number of types of ORV; j is the corresponding series number; and γ i , j M C and γ i , j E V refer to the results calculated by MCS and the evaluated method, respectively.

4.1. IEEE-30 Test System

A modified IEEE-30 [30] test system is used here to test the accuracy performance of the proposed method. This system contains six generators, 41 lines, and 20 loads.
Four wind farms are integrated in the test system. The probability distribution of wind farm output is assumed to follow the beta distribution, as this distribution is a proper choice to describe the short-term uncertainties of wind power [31], and the power factor of wind farms is assumed to be 0.95. The penetration level of wind power output is approximately 30%. The located bus numbers and uncertainty descriptions of wind farms are shown in Table 2.
The uncertainty of the load at each bus is assumed to follow a Gaussian distribution, with the means of active and reactive power injections equal to those of the base case data and standard deviations equal to 5% of the means. Accordingly, the total number of IRVs is 48 (including the active and reactive injections of four wind farms and 20 loads). All IRVs are assumed to be independent of each other in this case.
The performance of the evaluated methods is shown in Table 3. ε ¯ 3 and ε ¯ 4 , which represent the estimation error of high-order moments of evaluated methods, are shown in Figure 1.
The data presented in Table 3 and Figure 1 show the following results:
(1) HDMR has the minimal estimation error of the first four moments compared to other candidates, but at the cost of a much longer calculation time (more than 50 s);
(2) The estimation accuracy of Zhao’s PEM scheme for the first two orders of moments is satisfactory, but the estimation error is extremely high when estimating third- and fourth-order moments;
(3) Combining PCA with HDMR can be very effective in improving computational efficiency. When the contribution threshold is set at 90%, the calculation time of PCA+HDMR is even shorter than that of the PEM;
(4) Reducing the contribution threshold results in a decrease in computational accuracy. This is because some of the information useful for moment estimation is discarded by PCA. Therefore, in practice, the contribution threshold of PCA should be set carefully. When setting the contribution threshold, we need to balance the conflicting goals of precision and efficiency. The result shows that 99% is a good choice in this case.
Note that here, the estimation error of the third-order moment is relatively high for all evaluated methods. This is because the absolute value of the third-order moment is usually very close to 0; therefore, a small variation will result in a large relative error; however, this kind of large relative error has a minor influence on describing the PDF of a random variable.

4.2. Sensitivity to Correlations

In this section, the outputs of wind farms are assumed to be correlated. We will test the adaptability of the proposed method to different correlation scenarios. Three correlation scenarios are assumed in this section: low correlation, medium correlation and high correlation. The correlation matrix is as follows:
ρ l o w = [ 1 0.25 0.20 0.15 0.25 1 0.23 0.18 0.20 0.23 1 0.30 0.15 0.18 0.30 1 ]
ρ m i d = [ 1 0.43 0.55 0.39 0.43 1 0.60 0.46 0.55 0.60 1 0.52 0.39 0.46 0.52 1 ]
ρ h i g h = [ 1 0.72 0.83 0.69 0.72 1 0.78 0.81 0.83 0.78 1 0.74 0.69 0.81 0.74 1 ]
As seen in Table 3, we found that a contribution threshold setting of 99% is good in terms of both accuracy and efficiency; therefore, the contribution threshold of these three scenarios is set at 99%.
The AREIs of the proposed method in the three correlation scenarios are shown in Table 4. ε ¯ 3 and ε ¯ 4 of these three scenarios are shown in Figure 2.
As the data show in Table 4 and Figure 2, in the three correlation scenarios, the estimated error of the proposed algorithm does not fluctuate greatly, so it can be concluded that the proposed method has good adaptability to different correlation scenarios.

4.3. IEEE-118 Test System

In this section, the adaptability of the proposed method to large power systems is tested. The IEEE 118 system contains 54 generators, 186 lines, and 100 loads. Eight wind farms are integrated in this test system. Accordingly, the total number of IRVs is 216. The located bus numbers and uncertainty descriptions of the wind farms are shown in Table 5. The contribution threshold of PCA is set at 99%, and other assumptions and settings are the same as those of the IEEE-30 test system. The performance of the evaluated methods is shown in Table 6.
The data presented in Table 6 show the following results:
(1) Similar to the results of the IEEE-30 test system, HDMR has the best accuracy performance, followed by PCA+HDMR, and the PEM has the worst performance. Although HDMR has the highest calculation accuracy, its computational burden becomes unacceptable in large systems (up to 406 s in the IEEE-118 test system);
(2) Zhao’s PEM still cannot effectively estimate third- and fourth-order moments. Clearly, ignoring the cooperative effects of the IRVs significantly increases the estimation error of high-order moments;
(3) PCA+HDMR performs high-order moment estimation with a small loss of accuracy. The calculation time in this test system is only 6.52 s. With a comprehensive consideration of calculation accuracy and calculation efficiency, PCA+HDMR is seen to be more suitable for practical use than the other methods.

5. Discussions

5.1. Non-Linear Dependency

In this paper, dependency among IRVs are assumed to be linear. For some cases in which IRVs have non-linear dependencies, using PCA may cause additional errors, since non-linear dependencies cannot be fully eliminated by PCA. In this situation, kernel PCA (kPCA) [32] may be used instead of PCA to increase accuracy. The basic idea of kPCA is to map the data in the input space to a feature space via a non-linear map and then apply PCA in the feature space.

5.2. Probability Distribution Approximation

For practical engineering applications, knowing only the moments of the ORVs may not be sufficient: their PDF or CDF may also be of interest. Gram–Chalier series [10] and Cornish–Fisher series [11] are commonly used in the literature to generate PDF or CDF curves according to moment information. However, there are some defects in such series-based methods. One is that the terms of the series are hard to determine in advance, and the other is that these two series may generate illegitimate negative probability distribution functions [13]. Versatile distributions [33], the Johnson system [8] and Gaussian mixture models [3,34] are good candidates that can provide better performance than series-based methods.

6. Conclusions

To overcome the defects of the traditional PEM, a new PEM approach based on PCA and HDMR is proposed in this study. The proposed method uses PCA to decrease the dimension scale of IRVs and eliminate correlations among IRVs, which makes HDMR adaptable to power systems that have a large number of IRVs. HDMR is used to estimate ORV moments; in this manner, a defect of the traditional PEM, that it has difficulty estimating high-order moments, is eliminated. Case studies show that the proposed method can achieve better performance in terms of accuracy and efficiency than the PEM and MCS when considering a large number of correlated, non-Gaussian IRVs.
In the future, other methods, such as kPCA and parallel computing, can be studied and combined with the proposed method to further enhance its performance. The proposed method may also prove useful for other PPF-based applications.

Author Contributions

H.L. conceived and designed the study; this work was performed under the guidance of Z.Z. and X.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors gratefully acknowledge the support of the National High Technology Research and Development Program of China (863 Program): No. 2015AA050201.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Borkowska, B. Probabilistic load flow. IEEE Trans. Power Appar Syst. 1974, 93, 752–759. [Google Scholar] [CrossRef]
  2. Henneaux, P.; Daniel, S. Probabilistic security analysis of optimal transmission switching. IEEE Trans. Power Syst. 2015, 31, 508–517. [Google Scholar] [CrossRef]
  3. Wang, Z.; Shen, C.; Liu, F.; W, X.; Liu, C.; Gao, F. Chance-constrained economic dispatch with non-Gaussian correlated wind power uncertainty. IEEE Trans. Power Syst. 2017, 32, 4880–4893. [Google Scholar] [CrossRef]
  4. Prusty, B.R.; Jena, D. A critical review on probabilistic load flow studies in uncertainty constrained power systems with photovoltaic generation and a new approach. Renew. Sustain. Energy Rev. 2017, 69, 1286–1302. [Google Scholar] [CrossRef]
  5. Wu, W.; Wang, K.; Li, G.; Jiang, X.; Feng, L.; Crow, M. A multivariate dimension-reduction method for probabilistic power flow calculation. Electr. Power Syst. Res. 2016, 141, 272–280. [Google Scholar] [CrossRef]
  6. Cai, D.; Shi, D.; Chen, J. Probabilistic load flow with correlated input random variables using uniform design sampling. Int. J. Electr. Power Energy Syst. 2014, 63, 105–112. [Google Scholar] [CrossRef]
  7. Chen, Y.; Wen, J.; Cheng, S. Probabilistic load flow method based on Nataf transformation and Latin hypercube sampling. IEEE Trans. Sustain. Energy 2012, 4, 294–301. [Google Scholar] [CrossRef]
  8. Zhang, L.; Cheng, H.; Zhang, S.; Zeng, P.; Yao, L. Probabilistic power flow calculation using the Johnson system and Sobol’s quasi-random numbers. IET Gener. Transm. Distrib. 2016, 10, 3050–3059. [Google Scholar] [CrossRef]
  9. Allan, R.N.; Al-Shakarchi, M. Probabilistic a.c. load flow. Proc. IEE 1976, 123, 531–536. [Google Scholar] [CrossRef]
  10. Zhang, P.; Stephen, 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] [CrossRef]
  11. Usaola, J. Probabilistic load flow with correlated wind power injections. Electr. Power Syst. Res. 2010, 80, 528–536. [Google Scholar] [CrossRef]
  12. Williams, T.; Curran, C. Probabilistic load flow modeling comparing maximum entropy and Gram-Charlier probability density function reconstructions. IEEE Trans. Power Syst. 2012, 28, 272–280. [Google Scholar] [CrossRef]
  13. Fan, M.; Vittal, V.; Heydt, G.; Ayyanar, R. Probabilistic power flow studies for transmission systems with photovoltaic generation using cumulants. IEEE Trans. Power Syst. 2012, 27, 2251–2261. [Google Scholar] [CrossRef]
  14. Wu, W.; Wang, K.; Li, G.; Jiang, X.; Wang, Z. Probabilistic load flow calculation using cumulants and multiple integrals. IET Gener. Transm. Distrib. 2016, 10, 1703–1709. [Google Scholar] [CrossRef]
  15. Rosenblueth, E. Point estimates for probability moments. Proc. Natl. Acad. Sci. 2013, 28, 3114–3125. [Google Scholar] [CrossRef] [Green Version]
  16. Zhao, Y.; Ono, T. New Point Estimates for Probability Moments. J. Eng. Mech. 2000, 126, 433–436. [Google Scholar] [CrossRef]
  17. Hong, H. An efficient point estimate method for probabilistic analysis. Reliab. Eng. Syst. Saf. 1998, 59, 261–267. [Google Scholar] [CrossRef]
  18. Su, C. Probabilistic load-flow computation using point estimate method. IEEE Trans. Power Syst. 2005, 20, 1843–1851. [Google Scholar] [CrossRef]
  19. Morales, J.; Perez-Ruiz, J. Point estimate schemes to solve the probabilistic power flow. IEEE Trans. Power Syst. 2007, 22, 1594–1601. [Google Scholar] [CrossRef]
  20. Liu, Y.; Gao, S.; Cui, H.; Yu, L. Probabilistic load flow considering correlations of input variables following arbitrary distributions. Electr. Power Syst. Res. 2010, 140, 354–362. [Google Scholar] [CrossRef]
  21. Chen, C.; Wu, W.; Zhang, B.; Sun, H. Correlated probabilistic load flow using a point estimate method with Nataf transformation. Int. J. Electr. Power Energy Syst. 2015, 65, 325–333. [Google Scholar] [CrossRef]
  22. Xiao, Q.; He, Y.; Yang, Y.; Lu, Y. Point estimate method based on univariate dimension reduction model for probabilistic power flow computation. IET Gener. Transm. Distrib. 2017, 11, 3522–3531. [Google Scholar] [CrossRef]
  23. Le, D.; Gross, G.; Berizzi, A. Probabilistic modeling of multisite wind farm production for scenario-based applications. IEEE Trans. Sustain. Energy 2016, 6, 748–758. [Google Scholar] [CrossRef] [Green Version]
  24. Saunders, C. Point estimate method addressing correlated wind power for probabilistic optimal power flow. IEEE Trans. Power Syst. 2013, 29, 1045–1054. [Google Scholar] [CrossRef]
  25. Li, H.; Lu, Z.; Yuan, X. Nataf transformation based point estimate method. Chin. Sci. Bull. 2008, 53, 2586–2592. [Google Scholar] [CrossRef]
  26. Xu, H.; Rahman, S. A generalized dimension-reduction method for multidimensional integration in stochastic mechanics. Int. J. Numer. Methods Eng. 2004, 61, 1992–2019. [Google Scholar] [CrossRef] [Green Version]
  27. Chowdhury, R.; Rao, B.; Prasad, A. High-dimensional model representation for structural reliability analysis. Commun. Number. Meth. Eng. 2009, 25, 301–337. [Google Scholar] [CrossRef]
  28. Fan, W.; Wei, J.; Ang, A.; Li, Z. Adaptive estimation of statistical moments of the responses of random systems. Probab. Eng. Eng. Mech. 2016, 43, 50–67. [Google Scholar] [CrossRef]
  29. Zimmerman, R.; Murillo-Sánchez, C.; Thomas, R. MATPOWER: Steady-state operations, planning, and analysis tools for power systems research and education. IEEE Trans. Power Syst. 2010, 26, 12–19. [Google Scholar] [CrossRef] [Green Version]
  30. Power System Test Case Archive. Available online: https://labs.ece.uw.edu/pstca/ (accessed on 1 March 2020).
  31. Bludszuweit, H.; Dominguez-Navarro, J.A.; Llombart, A. Statistical analysis of wind power forecast error. IEEE Trans. Power Syst. 2008, 23, 983–991. [Google Scholar] [CrossRef]
  32. Kwok, J.; Tsang, I. The Pre-Image Problem in Kernel Methods. IEEE Trans. Neural Netw. 2004, 15, 1517–1525. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Zhang, Z.; Sun, Y.; Gao, D.; Lin, J.; Cheng, L. A Versatile Probability Distribution Model for Wind Power Forecast Errors and Its Application in Economic Dispatch. IEEE Trans. Power Syst. 2013, 28, 3114–3125. [Google Scholar] [CrossRef]
  34. Xie, X.; Schenkendorf, R. Robust optimization of a pharmaceutical freeze-drying process under non-Gaussian parameter uncertainties. Chem. Eng. Sci. 2019, 207, 805–819. [Google Scholar] [CrossRef]
Figure 1. (a) ε ¯ 3 of evaluated methods and (b) ε ¯ 4 of evaluated methods.
Figure 1. (a) ε ¯ 3 of evaluated methods and (b) ε ¯ 4 of evaluated methods.
Energies 13 03520 g001
Figure 2. (a) ε ¯ 3 of three correlation scenarios and (b) ε ¯ 4 of three correlation scenarios.
Figure 2. (a) ε ¯ 3 of three correlation scenarios and (b) ε ¯ 4 of three correlation scenarios.
Energies 13 03520 g002
Table 1. Weights and abscissas for the Gauss–Hermite quadrature.
Table 1. Weights and abscissas for the Gauss–Hermite quadrature.
NumbersAbscissasWeights
301.1816
±1.22472.9541 × 10−1
509.4530 × 10−1
±0.95853.9362 × 10−1
±2.02011.9953 × 10−2
708.1027 × 10−1
±0.81624.2561 × 10−1
±1.67355.4512 × 10−2
±2.65199.7178 × 10−4
Table 2. Uncertainty descriptions of wind farms (IEEE-30).
Table 2. Uncertainty descriptions of wind farms (IEEE-30).
WFBusDistribution ParametersFluctuation Range (MW)
115Beta (5.32, 7.34)0–25
216
326Beta (4.18, 1.80)0–18
430
Table 3. Performance of evaluated methods (IEEE-30).
Table 3. Performance of evaluated methods (IEEE-30).
Method AREI Calculation Time (s)
Type of ORV ε ¯ 1 ε ¯ 2 ε ¯ 3 ε ¯ 4
Zhao’s PEMP0.09441.300280.224347.43761.27
Q0.12732.415543.368048.1434
V4.1262 × 10−41.6993189.070169.5160
θ 0.18770.6081125.363245.0952
HDMRP0.05640.241910.48861.699852.13
Q0.07270.09226.83380.9072
V3.6439 × 10−40.112411.45660.6689
θ 0.14170.085918.45501.6844
PCA+HDMR
( M k     99 % )
P0.11321.344712.72401.92754.98
Q0.19860.46414.31471.0406
V5.3721 × 10−40.410711.47070.7315
θ 0.09670.455021.35741.5991
PCA+HDMR
( M k     95 % )
P0.05643.866116.18101.56831.79
Q0.116611.958731.11810.7725
V4.1600 × 10−49.970861.75360.7817
θ 0.19022.624530.04431.4042
PCA+HDMR
( M k     90 % )
P0.07049.201422.74681.39030.89
Q0.054420.547723.70190.9561
V3.4794 × 10−419.526389.38381.0878
θ 0.21066.132222.49201.2454
MCS 164.86
Table 4. Sensitivity to correlation of proposed method.
Table 4. Sensitivity to correlation of proposed method.
MethodAREI
Type of ORV ε ¯ 1 ε ¯ 2 ε ¯ 3 ε ¯ 4
Low correlationP0.10111.37969.52182.6647
Q0.09562.71767.33811.7268
V3.9626 × 10−41.461612.77011.9133
θ 0.27160.176419.18710.6482
Medium correlationP0.03101.915311.00521.5038
Q0.06843.01628.62511.8370
V1.8130 × 10−42.010610.86090.9980
θ 0.08831.256420.14251.6366
High correlationP0.07821.807811.98011.1834
Q0.06583.00017.53881.6983
V4.7408 × 10−41.822210.93580.9288
θ 0.21066.132222.49201.2454
Table 5. Uncertainty descriptions of wind farms (IEEE-118).
Table 5. Uncertainty descriptions of wind farms (IEEE-118).
WFBusDistribution ParametersFluctuation Range (MW)
1–42, 3, 44, 50Beta (3.76, 5.82)0–300
5–882, 88, 98, 115Beta (6.82, 2.44)0–160
Table 6. Performance of evaluated methods (IEEE-118).
Table 6. Performance of evaluated methods (IEEE-118).
MethodAREICalculation Time (s)
Type of ORV ε ¯ 1 ε ¯ 2 ε ¯ 3 ε ¯ 4
Zhao’s PEMP0.03793.3787119.740551.64103.29
Q0.03348.2555163.154052.8882
V2.9447 × 10−44.4120155.389637.7793
θ 0.00881.201278.024066.3490
HDMRP0.02860.17079.36641.8081406.33
Q0.02190.521212.86643.1527
V1.8866 × 10−40.432017.45103.1466
θ 0.00580.144215.23320.5527
PCA+HDMR
( M k     99 % )
P0.03792.412511.47932.08096.52
Q0.03147.923318.90193.6851
V3.3919 × 10−43.422319.65084.8801
θ 0.00851.029817.09020.7055
MCS 391.62 164.86

Share and Cite

MDPI and ACS Style

Li, H.; Zhang, Z.; Yin, X. A Novel Probabilistic Power Flow Algorithm Based on Principal Component Analysis and High-Dimensional Model Representation Techniques. Energies 2020, 13, 3520. https://doi.org/10.3390/en13143520

AMA Style

Li H, Zhang Z, Yin X. A Novel Probabilistic Power Flow Algorithm Based on Principal Component Analysis and High-Dimensional Model Representation Techniques. Energies. 2020; 13(14):3520. https://doi.org/10.3390/en13143520

Chicago/Turabian Style

Li, Hang, Zhe Zhang, and Xianggen Yin. 2020. "A Novel Probabilistic Power Flow Algorithm Based on Principal Component Analysis and High-Dimensional Model Representation Techniques" Energies 13, no. 14: 3520. https://doi.org/10.3390/en13143520

APA Style

Li, H., Zhang, Z., & Yin, X. (2020). A Novel Probabilistic Power Flow Algorithm Based on Principal Component Analysis and High-Dimensional Model Representation Techniques. Energies, 13(14), 3520. https://doi.org/10.3390/en13143520

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