Next Article in Journal
Assessment of Credible Capacity for Intermittent Distributed Energy Resources in Active Distribution Network
Next Article in Special Issue
A Data Analysis Technique to Estimate the Thermal Characteristics of a House
Previous Article in Journal
A New Methodology for Assessing the Energy Consumption of Building Stocks
Previous Article in Special Issue
The Comparison Study of Short-Term Prediction Methods to Enhance the Model Predictive Controller Applied to Microgrid Energy Management
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A PMU-Based Method for Smart Transmission Grid Voltage Security Visualization and Monitoring

Department of Electrical Engineering, Feng Chia University (FCU), No. 100, Wenhwa Road, Seatwen, Taichung 40724, Taiwan
*
Author to whom correspondence should be addressed.
Energies 2017, 10(8), 1103; https://doi.org/10.3390/en10081103
Submission received: 26 June 2017 / Revised: 18 July 2017 / Accepted: 25 July 2017 / Published: 27 July 2017
(This article belongs to the Collection Smart Grid)

Abstract

:
With the rapid growth of usage of phasor measurement units (PMUs) for modern power grids, the application of synchronized phasors (synchrophasors) to real-time voltage security monitoring has become an active research area. This paper presents a novel approach for fast determination of loading margin using PMU data from a wide-area monitoring system (WAMS) to construct the voltage stability boundary (VSB) of a transmission grid. Specifically, a new approach for online loading margin estimation that considers system load trends is proposed based on the Thevenin equivalent (TE) technique and the Mobius transformation (MT) technique. A VSB is then computed by means of real-time PMU measurements and is presented in a complex load power space. VSB can be utilized as a visualization tool that is able to provide real-time visualization of the current voltage stability situation. The proposed method is fast and adequate for online voltage security assessment. Furthermore, it enables us to significantly increase a system operator’s situational awareness for operational decision making. Simulation studies were carried out using different sized power grid models under various operating conditions. The simulation results are shown to validate the capability of the proposed method.

1. Introduction

Recently, voltage stability has become a challenging issue in power industries since several major blackouts worldwide have been mainly attributed to voltage collapse [1,2,3]. Based on this concern, much attention has been devoted to understanding the voltage collapse phenomenon. Voltage collapse is described as the phenomenon in which the events in sequential operations are accompanied by voltage instability. This may lead to blackouts or abnormally low voltage distribution in a vital part of the power grid [4,5]. Therefore, it is essential that system operators easily identify system voltage instability to prevent possible voltage collapse.
In view of this, a variety of methods have been proposed for static voltage security analysis, such as sensitivity methods [6,7], continuation power flow (CPF) methods [8,9], and impedance match methods [10,11,12,13,14,15]. Indeed, the Power-Voltage (P-V) curve is widely utilized for voltage stability analysis. Generally, such a curve is produced by calculating power flow solutions for successively increasing load levels. To cope with the divergence issue of power flow calculation near the point of voltage collapse, the CPF methods are presented [8,9]. Since the CPF methods belong to model-based tools, a large amount of computation time for a bulk power grid model is required. In addition, an accurate system model is needed for this kind of model-based approach. According to the description above, CPF methods are not suitable for online voltage stability analysis in its current setting. In contrast, impedance match methods are measurement-based tools for online voltage stability monitoring [10,11,12,13,14,15]. Such methods are based on Thevenin equivalent (TE) impedance matching criterion using measured voltage and current phasors. As voltage instability take places at a load bus, the load impedance is equal to the TE impedance in magnitude. As a consequence, the difference between the load impedance and TE impedance can be utilized as an index for voltage stability monitoring. Even though the impedance margin can be employed to determine whether the present operating condition is voltage secure, loading margin information, which is more useful for system planners and operators, becomes obscure. Hence, an alternative approach is required.
The wide area monitoring system (WAMS) is expected to be a core technology in future smart grids [16,17]. With the gradual installation of phasor measurement units (PMUs) in many power grids, the PMU-based WAMS is becoming indispensable to sustaining reliability and security for modern power grids. PMUs, which are an advanced and fully developed technology, have been produced by many international manufacturers. Furthermore, most of the PMUs available on the market follow the IEEE standard C37.118, which defines frequency, rate of change of frequency, and synchronized phasors (synchrophasors) under all power system operating conditions. A more detailed description about this standard can be found in [18]. Figure 1 illustrates a typical structure of a PMU-based WAMS in a smart grid. Indeed, the wide-area measurements from PMUs, which is coherent and time-synchronized, can help increase real-time situational awareness for system operators to carry out time-critical tasks under critical conditions. This indicates that the PMU-based WAMS serves great opportunity to reconfigure the power system before reaching the voltage collapse point and avoid power system blackouts.
This research is mainly devoted to smart transmission grid voltage security visualization and monitoring. The contribution of this paper is the development of a novel method which aims to advance wide-area situational awareness enhanced with voltage stability monitoring. To this end, a new algorithm based on the Thevenin equivalent (TE) technique and the Mobius transformation (MT) technique for online system loading margin estimation is proposed. The approximation of a voltage stability boundary (VSB) can then be obtained by a quadratic curve using PMU data. The computed VSB can provide real-time visualization of the current voltage stability situation. This enables a significant increase in the system operator’s situational awareness for operational decision making. Furthermore, the proposed method is simple and fast, which makes it adequate for online applications.
The rest of this work is as follows: The impedance match method to static voltage stability analysis is briefly reviewed in Section 2. Section 3 describes the proposed method, and some computational details are discussed. The proposed method is tested on several Electrical and Electronics Engineers (IEEE) power grid models under various operating conditions in Section 4. Section 5 gives some concluding remarks.

2. Review of the Impedance Match Method to Static Voltage Stability Analysis

Impedance match methods are measurement-based tools for static voltage stability analysis [10,11,12,13,14,15]. Such methods use the concept of tracking the Thevenin equivalent (TE) parameters for a power system at a node using local voltage and current phasor measurements.
A complex power grid seen at a load bus i can be simplified to a TE network, shown in Figure 2, in which E i t h is the TE voltage, Z i t h is the TE impedance, and Z i L is the load impedance at load bus i. Note that E i t h , Z i t h , and Z i L are all in phasor representation,
E i t h = | E i t h | β i Z i t h = | Z i t h | α i Z i L = | Z i L | θ i
In this research, the phasors ( E i t h , Z i t h , Z i L ) are in per units, the magnitudes ( | E i t h | , | Z i t h | , | Z i L | ) are in per units, and the angles ( β i , α i , θ i ) are in degrees.
A key idea of the impedance match method is summarized as follows: when voltage collapse occurs, | Z i L | = | Z i t h | . Consequently, the ratio of | Z i t h | / | Z i L | can be utilized as an index for voltage stability monitoring, as shown in Figure 3. If the ratio is smaller than 1, the system will be stable. On the other hand, if the ratio is higher than 1, the system is unstable.
Due to the elegance and simplicity of this impedance-matching concept, the measurement-based real-time voltage stability analysis becomes possible [10,11,12,13,14,15]. The impedance-based index methods can be adopted to identify if a given operating condition is within a secure voltage margin or not. The information of the loading margin, however, becomes obscure. To address these difficulties, a new method, which is able to provide loading margin, is proposed in this research.

3. Proposed Method

This section presents a new method to transmission grid voltage stability visualization and monitoring utilizing PMU measurements. The main steps in the proposed method are shown in Figure 4 and are described in the following subsections.

3.1. Determination of Load Increase Direction

Loading margin estimations depend on load levels and load increase directions (LIDs). To illustrate this, Figure 5 shows the voltage collapse (critical point) surface in a load power parameter space with active and reactive load powers as coordinates. In Figure 5, let S 0 be the initial point, and Δ S 1 and Δ S 2 be two distinct LIDs, respectively. If S 0 increases along the direction of Δ S 1 , the maximum loading point would be S 1 max . Instead, if S 0 follows Δ S 2 , the critical point would be S 2 max . Clearly, distinct LIDs result in distinct voltage collapse points. Thus, it is essential to estimate the voltage collapse point with consideration of the LID.
The complex load power change is typically expressed as,
Δ S i t = S i t S i t 1
where
S i t = V i t ( I i t )
S i t 1 = V i t 1 ( I i t 1 )
In the above equations, V i and I i denote voltage and current measurements at load bus i, and t denotes the tth sampling point. In this research, the LID is accomplished by using two sets of consecutive PMU measurements.

3.2. Estimation of Thevenin Equivalent Network Parameters

In Figure 2, application of Kirchoff’s voltage law to TE network yields,
E i t h = V i + Z i t h I i
where V i and I i denote the load voltage and current phasors. Suppose that the data points of V i and I i are available from the installed PMU. The unknown variables in Equation (5) are E i t h and Z i t h . Since the Thevenin equivalent parameters is approximately constant under various loading conditions [10,11,12,13,14,15], E i t h and Z i t h can be determined by using two data pairs ( V i t 1 , I i t 1 ) and ( V i t , I i t ) as,
E i t h = V i t 1 I i t V i t I i t 1 I i t I i t 1
Z i t h = V i t 1 V i t I i t I i t 1
From the preceding equations, TE network for a power system at a bus can be acquired. Moreover, the load impedance can be obtained by,
Z i L = | V i | 2 S i
where S i denotes the conjugate of the complex load power S i = V i I i * .
According to Equations (6)–(8), all parameters in TE network can be obtained by synchrophasor measurements.

3.3. Estimation of Loading Margin

In Thevenin equivalent circuit depicted in Figure 2, the real load power at bus i is given by,
P i = | Z i L | | E i t h | 2 cos θ i | Z i t h | 2 + 2 | Z i L | | Z i t h | cos ( α i β i ) + | Z i L | 2
As P i increases, the trajectory of | Z i L | and | Z i t h | is illustrated in Figure 6, in which Z i M denotes the impedance margin, and it is defined to be,
Z i M = | Z i L | | Z i t h |
Once | Z i L | and | Z i t h | are obtained, Z i M can be determined. From Figure 6, it is clear that Z i M = 0 when | Z i L | = | Z i t h | , which means that the maximum power transfer occurs at bus i. Furthermore, if the maximum loading point P i max is determined, the loading margin can be obtained by,
P i M = P i max P i 0
where P i 0 denotes the current operating point. With P i M , the voltage stability margin (VSM) is defined as,
V S M i = P i M P i 0
Since Z i M = | Z i L | | Z i t h | can be obtained by utilizing PMU measurements, the aim is to transform the real-time monitored Z i M into direct information about loading margin P i M . In this research, the Mobius transformation (MT) technique [19] is employed to find a transformation T such that P i M = T ( Z i M ) . Based on the MT technique, the transformation T ( Z i M ) of Z i M is a map of the form,
T ( Z i M ) = a × Z i M + b c × Z i M + d , a , b , c , d
where the coefficients a, b, c, and d have to be determined. Indeed, the transformation depends on the 2 × 2 complex matrix,
g = [ a b c d ]
and let T g be the transformation. This means that g is the matrix associated to T g . It is noteworthy that the composition of two transformations T g 1 and T g 2 is associative [19], i.e., the matrix associated to T g 1 T g 2 is given by,
g 1 2 = g 1 g 2 T g 1 T g 2 = T g 1 g 2
where the arrow “ ” means “maps to”. This technique is used to construct a transformation that maps Z i M to P i M .
To start with, let ϒ be the ratio of Z i M and V S M i , i.e.,
ϒ = Z i M V S M i
Substituting from Equation (12) for V S M i , P i M can be expressed as,
P i M = P i 0 Z i M ϒ
In addition, according to Equation (9) and the maximum power transfer theorem, Equation (16) becomes,
ϒ = 2 | Z i L | | Z i t h | [ 1 + cos ( α i θ i ) ] Z i M
From Equations (17) and (18), the transformation T g 2 such that P i M = T g 2 ( ϒ ) and the transformation T g 1 such that ϒ = T g 1 ( Z i M ) can be obtained. Thus, the matrices g 1 and g 2 associated to T g 1 and T g 2 respectively can be written as,
g 1 = [ 0 2 | Z i L | | Z i t h | [ 1 + cos ( α i θ i ) ] 1 0 ]
g 2 = [ 0 P i 0 Z i M 1 0 ]
Since g 1 g 2 T g 2 T g 1 ( Z i M ) , the composite function T g 1 g 2 = T g 2 ( T g 1 ( Z i M ) ) is given by,
T g 1 g 2 ( Z i M ) = a × Z i M + b c × Z i M + d
where
[ a b c d ] = [ 0 P i 0 Z i M 1 0 ] [ 0 2 | Z i L | | Z i t h | [ 1 + cos ( α i θ i ) ] 1 0 ] = [ P i 0 Z i M 0 0 2 | Z i L | | Z i t h | [ 1 + cos ( α i θ i ) ] ]
Figure 7 illustrates the construction of T g 1 g 2 such that P i M = T g 1 g 2 ( Z i M ) . Once P i M is determined, S i max can be calculated by S i max = P i max + j Q i max , where P i max = P i M + P i 0 and Q i max = P i max tan θ i . Thus, the maximum loading point of the entire system S L max is obtained by S L max = i S i max , where i { PQ } .

3.4. Construction of Voltage Stability Boundary

Based on the approximate quadratic property of the VSB in a complex load power space [20,21], the VSB can be constructed by means of the quadratic curve fitting technique. A quadratic function for the active and reactive load power is,
Q L ( P L ) = i = 0 2 c i ( P L ) i
where c i for i = 0 , 1 , 2 are parameters to be determined. As stressed earlier, distinct LIDs lead to distinct maximum loading points S L max . Given three different look-ahead LIDs { L I D k } k = 1 3 , the corresponding maximum loading point { S L k max = P L k max + j Q L k max } k = 1 3 can be quickly estimated by using the method proposed in Section 3.3. Thus, the quadratic curve parameters in Equation (23) can be found by,
[ Q L 1 max Q L 2 max Q L 3 max ] = [ 1 P L 1 max ( P L 1 max ) 2 1 P L 2 max ( P L 2 max ) 2 1 P L 3 max ( P L 3 max ) 2 ] [ c 0 c 1 c 2 ]
Using the values of c 0 , c 1 , and c 2 , a VSB in a complex load power space, as shown in Figure 8, can be obtained. Such a VSB is able to give system operators a global view of power system voltage stability. The computation process to construct a VSB curve in P-Q plane is given in Algorithm 1.
Algorithm 1. Construction of a VSB curve in P-Q plane.
1: Input: { ( V i t , I i t ) } t = 1 2 and { L I D k } k = 1 3 ;
2: for k = 1 , 3 do
3: for i { P Q } do
4: P i 0 Re ( V i I i ) ;
5: θ i Arg ( V i / I i ) ;
6: Z i M ImpedanceMargin( V i , I i );
7: P i M PowerMargin( Z i M , P i 0 );
8: P i max P i M + P i 0 ;
9: Q i max P i max tan θ i ;
10: S i max P i max + j Q i max ;
11: end for
12: S L k max i S i max
13: end for
14: Determine { c j } j = 0 2 by Equation (24);
15: Construct a VSB curve by Equation (23);
16: return VSB curve
17: function ImpedanceMargin( V i , I i )
18: Compute Z i t h by Equation (7);
19: Compute Z i L by Equation (8);
20: Z i M | Z i L | | Z i t h | ;
21: return Z i M
22: end function
23: function PowerMargin( Z i M , P i 0 )
24: Compute [ a b c d ] by Equation (22);
25: Compute P i M by Equation (21);
26: return P i M
27: end function

4. Case Studies

The capability of the proposed algorithm is further demonstrated by utilizing different-size test systems. The system data and the simulation scenarios are listed in Table 1. Two metrics, accuracy and execution time, are utilized to evaluate the performance of the proposed method. Furthermore, its performance is also compared to the CPF method [8] that is widely used in voltage stability assessment. All computations are performed on an Intel® Core™ i5, 1.7 GHz computer.

4.1. Load Change Cases

In this research, extensive simulation studies for load change cases have been performed, including various load levels and various load patterns. Table 2 lists some selected cases, where λ 0 represents the initial load level, ( Δ P L E v e n , Δ Q L E v e n ) denotes load increases at even load buses, ( Δ P L O d d , Δ Q L O d d ) denotes load increases at odd load buses, and ( Δ P L A l l , Δ Q L A l l ) denotes load increases at all load buses.

4.1.1. IEEE 14-Bus Model

In order to test the algorithm, the IEEE 14-bus model is first utilized. The diagram of the test model is depicted in Figure 9. This model involves 14 buses, 5 generation units, 9 loads, and 20 transmission lines [22].
The results of system VSM estimated by the proposed method under the selected test cases is illustrated in Figure 10, where the values of VSM corresponding to the test cases are depicted on the Y-axis of the figure. From Figure 10, one can see that the estimated VSMs are almost the same as the actual ones.

4.1.2. IEEE 30-Bus Model

The proposed method is applied to the IEEE 30-bus model, depicted in Figure 11. This test model contains 30 buses, 6 generation units, 24 loads, and 41 transmission lines [22].
Figure 12 depicts the simulation results under the test cases listed in Table 2. In Figure 12, one can see that the proposed algorithm is able to provide acceptable accuracy for VSM identification.

4.1.3. Statistical Evaluation

To test the capability of the proposed algorithm for online voltage stability assessment, many simulations were conducted. These include different-size test systems, various load levels, and various load patterns. Among the simulation studies, Table 3 compares the selected results of VSM, along with the execution time for the tested methods. Table 3 shows that the VSMs obtained by the proposed approach get closer to the ones obtained by the CPF approach. However, the proposed approach has less execution time in all cases compared to the CPF approach.
In order to compare the overall efficiency of all the tested methods, an index called the efficiency coefficient (EC) is defined as,
EC = | log 10 ( V S M ) | Time
This index combines execution time and accuracy, represented by the VSM based on the given LID. A larger value of EC indicates high efficiency, in terms of less execution time and higher accuracy. Table 3 shows the values of EC for the compared methods. Since a larger system requires more execution time, the value of EC decreases when the system size increases. In addition, Table 3 clearly shows that the proposed approach is significantly more efficient than the conventional CPF approach.

4.2. Topology Change Cases

The aim of this case study is to demonstrate the capability of the proposed method to deal with topology changes in the operating condition. The considered topology change cases include transmission line and generator outages.
Table 4 compares the results in terms of EC for the studied systems. In addition, the table shows several out-of-service cases for each test system. Table 4 shows that the efficiency of the proposed method is significantly higher than the CPF method. Indeed, the proposed method yields results with acceptable accuracy, despite its ease of implementation and low computational cost. Furthermore, we observe, based on extensive simulations carried out, that the proposed method is easily applicable to any power system, regardless of system size and configuration.

4.3. Visualization of VSB in P-Q Plane

To examine the effects of different look-ahead LIDs on voltage stability assessment, we performed many load patterns on the studied systems. Figure 13 depicts an illustrative case of VSB determination in PQ plane for the IEEE 14-bus model. The maximum transferable load powers (P, Q, and |S|) and the corresponding operating point can be easily identified in Figure 13 (i.e., a VSB curve can give more meaningful information to power engineers). This real-time visualization of voltage stability is achieved via the identification of a VSB curve using PMU measurements.
The criterion used for performance evaluation of VSB approximation is expressed as,
EC avg = 1 n c p i = 1 n c p EC i
where EC i is defined previously in Equation (25) and n c p denotes the total number of compared points on the VSB. In this study, n c p = 100 for each test system. Also note that the exact VSB is computed by point-by-point voltage collapse point computation. Table 5 gives the values of EC avg for the compared methods. Through quantitative comparisons with the CPF method, it can be seen that the proposed algorithm is capable of providing acceptable accuracy, but with much less execution time and higher efficiency. This means that a possible VSB, which can improve real-time situational awareness to enhance voltage stability monitoring, can be given quickly by the proposed method.

4.4. Effects of Measurement Inaccuracies

To verify the effects of measurement inaccuracies on VSB determination, we carried out a lot of test cases by simply adding random noise to the original measurements. The total number of compared points on the VSB is set to be 100 for each test system. The results in terms of average error are summarized in Table 6, showing that measurement errors indeed degrade the performance of the VSB determination. However, the issue of data pre-processing can be easily solved by different filtering methods available in studies [23,24].

5. Conclusions

A new PMU-based algorithm for online voltage security monitoring is developed to increase the system operator’s situational awareness for operational decision making. Based on the TE technique and the MT technique for determining loading margin of a transmission grid, the approximation of VSB can be obtained by a quadratic curve using PMU data. Visualization of VSB in a complex load power space gives system operators a global view of power system voltage stability. The proposed algorithm can cope with different kinds of load increase cases. In addition, the algorithm is simple and fast, which makes it adequate for online applications. Numerical test results, using different-size power grid models, illustrate the effectiveness, flexibility, and capability of the proposed algorithm.

Acknowledgments

This work was supported by the Ministry of Science and Technology of Taiwan, under Contract MOST 104-2221-E-035-041 and Contract MOST 106-3113-E-002-012.

Author Contributions

Heng-Yi Su proposed the approach and wrote the manuscript with input from Tzu-Yi Liu.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Andersson, G.; Donalek, P.; Farmer, R.; Hatziargyriou, N.; Kamwa, I.; Kundur, P.; Martins, N.; Paserba, J.; Pourbeik, P.; Sanchez-Gasca, J.; et al. Causes of the 2003 major grid blackouts in north America and Europe, and recommended means to improve system dynamic performance. IEEE Trans. Power Syst. 2005, 20, 1922–1928. [Google Scholar] [CrossRef]
  2. Simpson-Porco, J.W.; Dörfler, F.; Bullo, F. Voltage collapse in complex power grids. Nat. Commun. 2016. [Google Scholar] [CrossRef] [PubMed]
  3. Ajjarapu, V. Computational Techniques for Voltage Stability Assessment and Control; Springer-Verlag: New York, NY, USA, 2006. [Google Scholar]
  4. Su, H.-Y.; Hsu, Y.-L.; Chen, Y.-C. PSO-Based Voltage Control Strategy for Loadability Enhancement in Smart Power Grids. Appl. Sci. 2016, 6, 449. [Google Scholar] [CrossRef]
  5. Deng, W.; Zhang, B.; Ding, H.; Li, H. Risk-Based Probabilistic Voltage Stability Assessment in Uncertain Power System. Energies 2017, 10, 180. [Google Scholar] [CrossRef]
  6. Flatabo, N.; Ognedal, R.; Carlsen, T. Voltage stability condition in a power transmission system calculated by sensitivity methods. IEEE Trans. Power Syst. 1990, 5, 1286–1293. [Google Scholar] [CrossRef]
  7. Flatabo, N.; Fosso, O.; Ognedal, R.; Carlsen, T. A method for calculation of margins to voltage instability applied on the Norwegian system for maintaining required security level. IEEE Trans. Power Syst. 1993, 8, 920–928. [Google Scholar] [CrossRef]
  8. Ajjarapu, V.; Christy, C. The continuation power flow: A tool for steady state voltage stability analysis. IEEE Trans. Power Syst. 1992, 7, 416–423. [Google Scholar] [CrossRef]
  9. Chiang, H.D.; Flueck, A.J.; Shah, K.S.; Balu, N. CPFLOW: A practical tool for tracing power system steady-state stationary behavior due to load and generation variations. IEEE Trans. Power Syst. 1995, 10, 623–634. [Google Scholar] [CrossRef]
  10. Vu, K.; Begovic, M.M.; Novosel, D.; Saha, M.M. Use of local measurements to estimate voltage-stability margin. IEEE Trans. Power Syst. 1999, 14, 1029–1035. [Google Scholar] [CrossRef]
  11. Smon, I.; Verbic, G.; Gubina, F. Local voltage-stability index using Tllegen’s theorem. IEEE Trans. Power Syst. 2006, 21, 1267–1275. [Google Scholar] [CrossRef]
  12. Corsi, S.; Taranto, G.N. A real-time voltage instability identification algorithm based on local phasor measurements. IEEE Trans. Power Syst. 2008, 23, 1271–1279. [Google Scholar] [CrossRef]
  13. Verbic, G.; Gubina, F. A new concept of voltage-collapse protection based on local phasors. IEEE Trans. Power Del. 2004, 19, 576–581. [Google Scholar] [CrossRef]
  14. Parniani, M.; Vanouni, M. A fast local index for online estimation of closeness to loadability limit. IEEE Trans. Power Syst. 2010, 25, 584–585. [Google Scholar] [CrossRef]
  15. Sun, T.; Li, Z.; Rong, S.; Lu, J.; Li, W. Effect of Load Change on the Thevenin Equivalent Impedance of Power System. Energies 2017, 10, 330. [Google Scholar] [CrossRef]
  16. Phadke, A.G.; Thorp, J.S. Synchronized Phasor Measurements and Their Applications; Springer-Verlag: New York, NY, USA, 2008. [Google Scholar]
  17. Aminifar, F.; Fotuhi-Firuzabad, M.; Safdarian, A.; Davoudi, A.; Shahidehpour, M. Synchrophasor measurement technology in power systems: Panorama and state-of-the-art. IEEE Access 2014, 2, 1607–1628. [Google Scholar] [CrossRef]
  18. Martin, K.E. Synchrophasor standards development-IEEE C37.118 & IEC 61850. In Proceedings of the 44th Hawaii International Conference on System Sciences, Kauai, HI, USA, 4–7 January 2011; pp. 1–8. [Google Scholar]
  19. Dubrovin, B.A.; Fomenko, A.T.; Novikov, S.P. Modern Geometry-Methods and Applications; Springer-Verlag: New York, NY, USA, 1984. [Google Scholar]
  20. Taylor, C.W. Power System Voltage Stability; McGraw-Hill: New York, NY, USA, 1994. [Google Scholar]
  21. Van Cutsem, T.; Vournas, C. Voltage Stability of Electric Power Systems; Kluwer: Norwell, MA, USA, 1998. [Google Scholar]
  22. Power Systems Test Case Archive, University of Washington College of Engineering. Available online: http://www.ee.washington.edu/re-serach/pstcal/ (accessed on 14 April 2017).
  23. Tate, J.E.; Van, Cutsem T. Extracting steady state values from phasor measurement unit data using FIR and median filters. In Proceedings of the 2009 IEEE/PES Power Systems Conference and Exposition, Seattle, WA, USA, 15–18 March 2009; pp. 1–6. [Google Scholar]
  24. Glavic, M.; Van Cutsem, T. Wide-area detection of voltage instability from synchronized phasor measurements, part I: Principle. IEEE Trans. Power Syst. 2009, 24, 1408–1416. [Google Scholar] [CrossRef]
Figure 1. Illustration of phasor measurement unit (PMU)–based wide-area monitoring system (WAMS) in a smart grid.
Figure 1. Illustration of phasor measurement unit (PMU)–based wide-area monitoring system (WAMS) in a smart grid.
Energies 10 01103 g001
Figure 2. Thevenin equivalent (TE) network for a power system seen at a load bus.
Figure 2. Thevenin equivalent (TE) network for a power system seen at a load bus.
Energies 10 01103 g002
Figure 3. Illustration of the impedance match method for voltage stability monitoring.
Figure 3. Illustration of the impedance match method for voltage stability monitoring.
Energies 10 01103 g003
Figure 4. Main steps in the proposed algorithm for voltage stability visualization and monitoring.
Figure 4. Main steps in the proposed algorithm for voltage stability visualization and monitoring.
Energies 10 01103 g004
Figure 5. Load power parameter space.
Figure 5. Load power parameter space.
Energies 10 01103 g005
Figure 6. Trajectory of | Z i L | and | Z i t h | as P i increases.
Figure 6. Trajectory of | Z i L | and | Z i t h | as P i increases.
Energies 10 01103 g006
Figure 7. Illustration of the transformation between | Z i M | and P i M .
Figure 7. Illustration of the transformation between | Z i M | and P i M .
Energies 10 01103 g007
Figure 8. Voltage stability boundary (VSB) in a complex load power space.
Figure 8. Voltage stability boundary (VSB) in a complex load power space.
Energies 10 01103 g008
Figure 9. IEEE 14-bus test system.
Figure 9. IEEE 14-bus test system.
Energies 10 01103 g009
Figure 10. Estimation of system VSM for IEEE 14-bus model.
Figure 10. Estimation of system VSM for IEEE 14-bus model.
Energies 10 01103 g010
Figure 11. IEEE 30-bus test system.
Figure 11. IEEE 30-bus test system.
Energies 10 01103 g011
Figure 12. Estimation of system VSM for IEEE 30-bus model.
Figure 12. Estimation of system VSM for IEEE 30-bus model.
Energies 10 01103 g012
Figure 13. VSB in P-Q plane for IEEE 14-bus model. The point at (1.5, 1.35) is the initial operating point.
Figure 13. VSB in P-Q plane for IEEE 14-bus model. The point at (1.5, 1.35) is the initial operating point.
Energies 10 01103 g013
Table 1. System data and the simulation scenarios.
Table 1. System data and the simulation scenarios.
Number ofIEEE 14-BusIEEE 30-BusIEEE 57-BusIEEE 118-BusIEEE 300-Bus
Generators5675469
Loads9245064231
Lines204180186411
ScenarioDescription
Load changeVarious load levels together with various load patterns
Topology changeTransmission line and generator outages
Table 2. Selected load increase cases.
Table 2. Selected load increase cases.
CaseLoad LevelLoad Pattern
1 λ 0 = 1.2 ( Δ P L E v e n , Δ Q L E v e n ) = { ( Δ P i , Δ Q i ) | i { even   load   buses } }
2 λ 0 = 1.4
3 λ 0 = 1.6
4 λ 0 = 1.2 ( Δ P L O d d , Δ Q L O d d ) = { ( Δ P i , Δ Q i ) | i { odd   load   buses } }
5 λ 0 = 1.4
6 λ 0 = 1.6
7 λ 0 = 1.2 ( Δ P L A l l , Δ Q L A l l ) = { ( Δ P i , Δ Q i ) | i { all   load   buses } }
8 λ 0 = 1.4
9 λ 0 = 1.6
Table 3. Comparison of the results for the considered test systems in load change cases.
Table 3. Comparison of the results for the considered test systems in load change cases.
Test SystemVSM (%)Time (s)EC
ProposedCPFProposedCPFProposedCPF
IEEE 14-bus22.1722.480.251.365.380.99
IEEE 30-bus34.0134.630.412.593.740.59
IEEE 57-bus17.5717.961.066.041.170.21
IEEE 118-bus9.069.411.6210.340.590.09
IEEE 300-bus12.9413.2712.4473.420.090.02
Table 4. Comparison of the results for the considered test systems in topology change cases.
Table 4. Comparison of the results for the considered test systems in topology change cases.
Test SystemOut of ServiceEC
ProposedCPF
IEEE 14-busLine 9–105.251.36
Line 12–134.940.98
G25.031.01
IEEE 30-busLine 21–223.680.63
Line 10–173.310.47
G133.970.72
IEEE 57-busLine 23–241.780.59
Line 9–551.040.15
G81.360.34
IEEE 118-busLine 18–190.630.19
Line 63–640.580.11
G240.600.14
IEEE 300-busLine 11–130.140.05
Line 15–370.080.02
G100.110.04
Table 5. Comparison of the results in terms of ECavg for the considered test systems.
Table 5. Comparison of the results in terms of ECavg for the considered test systems.
Test System EC avg
ProposedCPF
IEEE 14-bus5.071.12
IEEE 30-bus3.650.64
IEEE 57-bus1.390.37
IEEE 118-bus0.610.15
IEEE 300-bus0.120.04
Table 6. Effects of measurement errors to VSB determination.
Table 6. Effects of measurement errors to VSB determination.
Test SystemError (%)
No Measurement ErrorsWith Measurement Errors
IEEE 14-bus−1.56−4.66
IEEE 30-bus−2.05−5.24
IEEE 57-bus−1.79−4.98
IEEE 118-bus−2.46−5.47
IEEE 300-bus−2.34−5.31

Share and Cite

MDPI and ACS Style

Su, H.-Y.; Liu, T.-Y. A PMU-Based Method for Smart Transmission Grid Voltage Security Visualization and Monitoring. Energies 2017, 10, 1103. https://doi.org/10.3390/en10081103

AMA Style

Su H-Y, Liu T-Y. A PMU-Based Method for Smart Transmission Grid Voltage Security Visualization and Monitoring. Energies. 2017; 10(8):1103. https://doi.org/10.3390/en10081103

Chicago/Turabian Style

Su, Heng-Yi, and Tzu-Yi Liu. 2017. "A PMU-Based Method for Smart Transmission Grid Voltage Security Visualization and Monitoring" Energies 10, no. 8: 1103. https://doi.org/10.3390/en10081103

APA Style

Su, H. -Y., & Liu, T. -Y. (2017). A PMU-Based Method for Smart Transmission Grid Voltage Security Visualization and Monitoring. Energies, 10(8), 1103. https://doi.org/10.3390/en10081103

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