A Synchrophasor Based Optimal Voltage Control Scheme with Successive Voltage Stability Margin Improvement

1 Department of Electrical Engineering, Feng Chia University (FCU), No.100, Wenhwa Road, Seatwen, Taichung 40724, Taiwan 2 Department of Information Engineering and Computer Science, Feng Chia University (FCU), No.100, Wenhwa Road, Seatwen, Taichung 40724, Taiwan; chenyic@fcu.edu.tw 3 Department of Automatic Control Engineering, Feng Chia University (FCU), No.100, Wenhwa Road, Seatwen, Taichung 40724, Taiwan; hsuyl@fcu.edu.tw * Correspondence: hengyisu@fcu.edu.tw; Tel.: +86-4-2451-7250 (ext. 3822); Fax: +86-4-2451-6842


Introduction
Voltage instability has been regarded as one of the primary threats to the security of modern power network operation during the past few decades.Power system disturbances such as a continuous load increase and/or a major change in network topology can result in voltage collapse.The voltage collapse problem, which is characterized by the loss of voltage magnitude at certain locations of the power grid, has caused several severe blackout events worldwide [1,2].
A number of planning and operation technologies have been proposed to mitigate the risk of voltage collapse [3].Among these technologies, phasor measurement unit (PMU) based schemes to secure power systems have become one of the enabling techniques which are under active investigations.Indeed, the use of PMUs in modern power systems becomes popular [4,5].
In order to provide a better voltage support in transmission networks, the coordinated voltage control has been developed.It is organized as a hierarchical structure with three levels: the primary, secondary, and tertiary voltage control.Significant attention has been given to the study of the secondary level, which is an automatic regulation of voltage and reactive power for power systems [6][7][8][9][10][11][12][13][14][15][16][17][18][19].The task of the secondary voltage control (SVC) is to regulate system voltage profile and to protect against potential voltage instability scenario at an early stage.The conventional approaches use only the voltage information at the observed buses (the so-called pilot nodes) as a triggering signal, but just monitoring voltage magnitudes alone may give the wrong indication in static voltage stability study [20].In other words, voltage stability issues cannot be fully prevented by the traditional SVC methods; therefore, an additional accepted measure of static voltage stability is required for a more reliable SVC scheme.The voltage stability margin (VSM), which is defined to be the distance between the given operating point and the voltage-collapse point, can serve as the static voltage stability indicator.A wide variety of approaches have been proposed for static VSM evaluation [21][22][23][24][25][26], wherein the measurement-based methods [25,26] address such problem by using the impedance-matching concept, i.e., at the maximum loading condition, load impedance is equal to Thevenin equivalent impedance in magnitude.Due to the elegance and simplicity of the measurement-based method, it becomes an attractive alternative in static voltage stability monitoring.In this paper, a VSM index computed from synchrophasor data is utilized.
This paper is concerned with designing an optimal voltage control scheme using synchrophasor measurements.In the proposed technique, the voltage magnitude and static voltage stability margin are considered as two key criteria for activating the proposed control strategy; therefore, a more efficient secondary voltage control can be accomplished.The rest of the paper is organized as follows: In Section 2, the proposed scheme based on the PMU technique is presented.Section 3 describes the fundamental theories and mathematical principles of the proposed scheme.Numerical simulations and test results are given and discussed in Section 4. Section 5 concludes the paper.

Overview of the Proposed Scheme
The overall architecture of the proposed scheme is shown in Figure 1.The functions and principles of the scheme are briefly described in this section.
stability study [20].In other words, voltage stability issues cannot be fully prevented by the traditional SVC methods; therefore, an additional accepted measure of static voltage stability is required for a more reliable SVC scheme.
The voltage stability margin (VSM), which is defined to be the distance between the given operating point and the voltage-collapse point, can serve as the static voltage stability indicator.A wide variety of approaches have been proposed for static VSM evaluation [21][22][23][24][25][26], wherein the measurement-based methods [25,26] address such problem by using the impedance-matching concept, i.e., at the maximum loading condition, load impedance is equal to Thevenin equivalent impedance in magnitude.Due to the elegance and simplicity of the measurement-based method, it becomes an attractive alternative in static voltage stability monitoring.In this paper, a VSM index computed from synchrophasor data is utilized.
This paper is concerned with designing an optimal voltage control scheme using synchrophasor measurements.In the proposed technique, the voltage magnitude and static voltage stability margin are considered as two key criteria for activating the proposed control strategy; therefore, a more efficient secondary voltage control can be accomplished.The rest of the paper is organized as follows: In Section 2, the proposed scheme based on the PMU technique is presented.Section 3 describes the fundamental theories and mathematical principles of the proposed scheme.Numerical simulations and test results are given and discussed in Section 4. Section 5 concludes the paper.

Overview of the Proposed Scheme
The overall architecture of the proposed scheme is shown in Figure 1.The functions and principles of the scheme are briefly described in this section.The designed scheme is illustrated with the part enclosed by the dashed line shown in Figure 1.Indeed, it consists of two primary functions: (i) voltage profile and static VSM monitoring function, which monitors not only voltage magnitudes but static VSM as well using synchrophasor data from PMUs; and (ii) adaptive secondary voltage control (ASVC) function which provides appropriate control actions to the controllers.The proposed scheme aims to regulate voltage profiles and to enlarge static voltage stability margin by means of a proper adjustment of VAR sources when the power grid has poor voltage profile and/or shows the risk to voltage collapse.
The flowchart of the proposed scheme is illustrated in Figure 2 and is described in the following: Step 1: In this study, we assume that PMUs are installed at the selected pilot buses.Using the precise timing signal provided by GPS as the common time base for PMUs, both magnitude and phase angle of voltage and current signals at different PMU locations can be measured, at exactly the same time instant from all observable system buses.Once the PMU measurements with time tags are entered into the designed scheme, these coherent and real-time measured quantities are then used to process the tasks of the proposed scheme.
Appl.Sci.2016, 6, 14 3 of 12 The designed scheme is illustrated with the part enclosed by the dashed line shown in Figure 1.Indeed, it consists of two primary functions: (i) voltage profile and static VSM monitoring function, which monitors not only voltage magnitudes but static VSM as well using synchrophasor data from PMUs; and (ii) adaptive secondary voltage control (ASVC) function which provides appropriate control actions to the controllers.The proposed scheme aims to regulate voltage profiles and to enlarge static voltage stability margin by means of a proper adjustment of VAR sources when the power grid has poor voltage profile and/or shows the risk to voltage collapse.
The flowchart of the proposed scheme is illustrated in Figure 2 and is described in the following: Step 1: In this study, we assume that PMUs are installed at the selected pilot buses.Using the precise timing signal provided by GPS as the common time base for PMUs, both magnitude and phase angle of voltage and current signals at different PMU locations can be measured, at exactly the same time instant from all observable system buses.Once the PMU measurements with time tags are entered into the designed scheme, these coherent and real-time measured quantities are then used to process the tasks of the proposed scheme.Step 2: At this stage, the proposed scheme performs voltage profile evaluation and static voltage stability margin estimation simultaneously.
(1) The voltage magnitudes, which are extracted from the PMUs in step 1, are used to evaluate the voltage magnitudes of the monitored buses.As long as the measured voltages are not within the specified limit, voltage violation is confirmed.Then, the scheme will issue a trigger signal to the ASVC function.Otherwise, the scheme will go back to Step 1. (2) As mentioned before, depending on bus voltages for the activation of control actions is not a very reliable strategy to be adopted.Another important index, which measures the Step 2: At this stage, the proposed scheme performs voltage profile evaluation and static voltage stability margin estimation simultaneously.
(1) The voltage magnitudes, which are extracted from the PMUs in step 1, are used to evaluate the voltage magnitudes of the monitored buses.As long as the measured voltages are not within the specified limit, voltage violation is confirmed.Then, the scheme will issue a trigger signal to the ASVC function.Otherwise, the scheme will go back to Step 1.
(2) As mentioned before, depending on bus voltages for the activation of control actions is not a very reliable strategy to be adopted.Another important index, which measures the proximity of an operating point to voltage instability, is the voltage stability margin (VSM).Therefore, the scheme needs to further check to see whether the current power system state is operated within a sufficient safety margin or not.If the computed value VSM is less than the predetermined security margin, inadequate VSM is identified.This will initiate the ASVC action accordingly.
Step 3: When critical voltage level or inadequate VSM is detected, the ASVC method will be activated automatically.The proposed control strategy can provide appropriate control actions to improve voltage profiles as well as voltage stability margin of the power system.That is, the proposed scheme is capable of preventing the system from possible voltage collapse.
The following section shows the detailed principles and methods used in the proposed scheme.

Principles of the Scheme
The proposed scheme is intended to automatically improve power system voltage profiles and static voltage stability margin simultaneously.The basic principles of the scheme will be derived in the following subsections.

Adaptive Secondary Voltage Control (ASVC) Function
First, consider that the approximate model of the small disturbance voltage-var control is represented by where ∆ Q and ∆ |V| are the reactive power and voltage magnitude change vectors; B stands for the system susceptance matrix; and the subscripts L and G denote the load and the voltage-controlled buses, respectively.In the above matrix equation, load voltage changes ∆ |V| can be expressed as where Notice that q and u are considered as reactive power load disturbances and control variables, respectively.
In this study, the feedback control law is used to develop the secondary voltage controller.The control vector is determined by using the worst-case design, the technique which addresses the problem by applying the minimization of the maximum load voltage deviation as an objective function for the optimal control model.Mathematically, the problem is to search an optimal set of control solutions such that the l 8 norm of load voltage changes is minimized.
Suppose that PMUs are placed at the chosen pilot nodes; the voltage deviations at the observed buses can be elaborated by ∆ ˇˇV p ˇˇ" ˇˇV * p ˇˇ´ˇˇVp ˇˇ (5) where ˇˇV * p ˇˇand ˇˇV p ˇˇdenote the pilot-node set-point voltages determined by tertiary level and pilot-node measured voltages obtained from installed PMUs, respectively.In the proposed secondary voltage control, both voltage profile and static voltage stability margin are adopted to serve as trigger Based on different numbers of pilot buses being used, the following conditions are considered.
(1) All Load Buses as Pilot Buses: This means that voltage measurements are available at all load buses.Under this circumstance, the first term 1 J q in Equation ( 2) can be replaced by p Δ V ; therefore, the problem formulation of the voltage control strategy is stated as where min G V and max G V are the lower and upper bounds of generator voltage magnitudes, respectively.
(2) Partial Load Buses as Pilot Buses: This means that voltage measurements are only available at the given pilot buses, and thus only has some elements of the vector under this situation.Indeed, can be related to q by the equation , in which Jp is the matrix with the rows of corresponding to the pilot points selected.
Since there are fewer measurements than variables to be estimated, the equation is underdetermined.In the proposed method, however, any reactive power disturbance variable vector q that satisfies can be utilized to determine feasible control signals that minimize . The least-norm technique [27], which is the most commonly used for solving an underdetermined set of linear equations, is applied to approximate q.In this case, the problem is to find q that satisfies and minimizes 2 q .Thus, the estimate of q is given by ( ) With the optimal q * , the proposed control strategy can be formulated in the following search problem: The constrained optimization problems in Equations ( 6) and ( 8) can be reformulated as linear programming problems [27] and be solved by using a linear programing solver such as linprog function in the MATLAB optimization toolbox [28].The optimal control actions, which are obtained Based on different numbers of pilot buses being used, the following conditions are considered.
(1) All Load Buses as Pilot Buses: This means that voltage measurements are available at all load buses.Under this circumstance, the first term J 1 q in Equation ( 2) can be replaced by ∆ ˇˇV p ˇˇ; therefore, the problem formulation of the voltage control strategy is stated as where ˇˇV min G ˇˇand ˇˇV max G ˇˇare the lower and upper bounds of generator voltage magnitudes, respectively.(2) Partial Load Buses as Pilot Buses: This means that voltage measurements are only available at the given pilot buses, and thus ∆ ˇˇV p ˇˇonly has some elements of the vector J 1 q under this situation.Indeed, ∆ ˇˇV p ˇˇcan be related to q by the equation J p q " ∆ ˇˇV p ˇˇ, in which J p is the matrix with the rows of J 1 corresponding to the pilot points selected.
Since there are fewer measurements than variables to be estimated, the equation J p q " ∆ ˇˇV p ˇˇis underdetermined.In the proposed method, however, any reactive power disturbance variable vector q that satisfies J p q " ∆ ˇˇV p ˇˇcan be utilized to determine feasible control signals that minimize || ∆ |V L | || 8 .The least-norm technique [27], which is the most commonly used for solving an underdetermined set of linear equations, is applied to approximate q.In this case, the problem is to find q that satisfies J p q " ∆ ˇˇV p ˇˇand minimizes ||q|| 2 .Thus, the estimate of q is given by With the optimal q*, the proposed control strategy can be formulated in the following search problem: min The constrained optimization problems in Equations ( 6) and ( 8) can be reformulated as linear programming problems [27] and be solved by using a linear programing solver such as linprog function in the MATLAB optimization toolbox [28].The optimal control actions, which are obtained from Equations ( 6) or (8), ensure that the resulting worst load voltage change can precisely remain within the range of the predefined constraint regardless of any unexpected load disturbances acting on the system.

Voltage Stability Margin Estimation Function
A common static voltage stability index is expressed by the voltage stability margin which shows how far the system is away from a possible instability event.In order to rapidly assess static voltage stability of a power grid, the local measurement-based methods have been presented in the works [25,26].The key idea of these techniques is provided here.
To start with, consider a load at bus i connected to a complex power system, which can be simplified to a single-machine-infinite-bus system by an estimated Thevenin equivalent network as shown in Figure 4, where E th i , Z th i , and Z L i correspond to Thevenin equivalent voltage, Thevenin equivalent impedance, and load impedance in phasor representation at bus i, respectively.When the power transmitted is maximum at bus i, ˇˇZ L i ˇˇis identical to ˇˇZ th i ˇˇ.Based on the impedance match theory, tracking ˇˇZ L i ˇˇand ˇˇZ th i ˇˇplays an important role in real-time voltage instability detection.
Appl.Sci.2016, 6, 14 6 of 12 from Equations ( 6) or (8), ensure that the resulting worst load voltage change can precisely remain within the range of the predefined constraint regardless of any unexpected load disturbances acting on the system.

Voltage Stability Margin Estimation Function
A common static voltage stability index is expressed by the voltage stability margin which shows how far the system is away from a possible instability event.In order to rapidly assess static voltage stability of a power grid, the local measurement-based methods have been presented in the works [25,26].The key idea of these techniques is provided here.
To start with, consider a load at bus i connected to a complex power system, which can be simplified to a single-machine-infinite-bus system by an estimated Thevenin equivalent network as shown in Figure 4, where th i E , th i Z , and L i Z correspond to Thevenin equivalent voltage, Thevenin equivalent impedance, and load impedance in phasor representation at bus i, respectively.When the power transmitted is maximum at bus i,  Using the approach above, the impedance based VSM index of bus i is defined as where Z can be obtained by using the local measured voltage and current phasors which are available from the installed PMU [25,26].The static voltage stability margin for the entire power system is defined to be where p represents the set of numbers of pilot buses.Note that the value of VSM is between 0 and 1.At the voltage-collapse point, the VSM is equal to 0. Although the approach proposed in [25] was employed to estimate VSM, it is noteworthy that any other measurement-based method of VSM analysis can be applied for this purpose, for example, the one proposed in [26].
The index VSM incorporated with the voltage magnitudes at the monitored buses will be used as the triggering signals for activating the proposed ASVC approach.Figure 5 illustrates the proposed control strategy, where 1 ε and 2 ε are the threshold levels for voltage change and security margin, respectively.In Figure 5, if the power system is operated at the critical state with poor voltage profile or insufficient VSM, then the ASVC function will be automatically initiated to steer the power system away from the point which is prone to voltage collapse.Using the approach above, the impedance based VSM index of bus i is defined as where ˇˇZ L i ˇˇand ˇˇZ th i ˇˇcan be obtained by using the local measured voltage and current phasors which are available from the installed PMU [25,26].The static voltage stability margin for the entire power system is defined to be VSM " min tVSM i | i " 1, 2, . . ., pu where p represents the set of numbers of pilot buses.Note that the value of VSM is between 0 and 1.At the voltage-collapse point, the VSM is equal to 0. Although the approach proposed in [25] was employed to estimate VSM, it is noteworthy that any other measurement-based method of VSM analysis can be applied for this purpose, for example, the one proposed in [26].
The index VSM incorporated with the voltage magnitudes at the monitored buses will be used as the triggering signals for activating the proposed ASVC approach.Figure 5 illustrates the proposed control strategy, where ε 1 and ε 2 are the threshold levels for voltage change and security margin, respectively.In Figure 5, if the power system is operated at the critical state with poor voltage profile or insufficient VSM, then the ASVC function will be automatically initiated to steer the power system away from the point which is prone to voltage collapse.

Simulation Results
This section presents numerical examples of the developed control scheme, i.e., simultaneous system voltage profiles improvement as well as static voltage stability margin enlargement, on a sample power network.The simulation program is coded using MATLAB ® and implemented on a personal computer which has a CPU with Intel ® Core™2 Duo 2.66 GHz and a memory with 4 GB.
The IEEE 30-bus system is used as an example to show the effectiveness of the proposed scheme.This sample system consists of 41 transmission lines, 6 VAR sources, and 24 loads.The system data including line parameters and bus data are given in [29].For all examples, the permissible value adopted for voltage magnitude is set at 0.95 1.05 V ≤ ≤ p.u., and the security margin for the test system is set as 0.3 VSM ≥ p.u..These are the threshold levels for triggering the proposed control algorithm.The voltage profile improvement index rms x used in this research is expressed by the root mean square value of voltage changes at all load buses, i.e., where m is the number of load buses in the test system, and j x for each 1, 2, , j m =  denotes the resulting voltage change at each of the load buses.In this test system, m is set to be 24.Note that the smaller the value of rms x , the better the voltage profile will be.In order to verify the performance of the presented methodology to power system secondary voltage control, we have studied a lot of experiments.These simulation cases include different load levels, different load patterns, various pilot-bus selections, and various branch outage contingencies.Among those investigated cases, some typical test results for the scenarios given in Table 1 are briefly summarized in the following.

Simulation Results
This section presents numerical examples of the developed control scheme, i.e., simultaneous system voltage profiles improvement as well as static voltage stability margin enlargement, on a sample power network.The simulation program is coded using MATLAB ® and implemented on a personal computer which has a CPU with Intel ® Core™2 Duo 2.66 GHz and a memory with 4 GB.
The IEEE 30-bus system is used as an example to show the effectiveness of the proposed scheme.This sample system consists of 41 transmission lines, 6 VAR sources, and 24 loads.The system data including line parameters and bus data are given in [29].For all examples, the permissible value adopted for voltage magnitude is set at 0.95 ď |V| ď 1.05 p.u., and the security margin for the test system is set as VSM ě 0.3 p.u..These are the threshold levels for triggering the proposed control algorithm.The voltage profile improvement index x rms used in this research is expressed by the root mean square value of voltage changes at all load buses, i.e., where m is the number of load buses in the test system, and x j for each j " 1, 2, . . ., m denotes the resulting voltage change at each of the load buses.In this test system, m is set to be 24.Note that the smaller the value of x rms , the better the voltage profile will be.
In order to verify the performance of the presented methodology to power system secondary voltage control, we have studied a lot of experiments.These simulation cases include different load levels, different load patterns, various pilot-bus selections, and various branch outage contingencies.Among those investigated cases, some typical test results for the scenarios given in Table 1 are briefly summarized in the following.

Case I
In the first case, bus load change pattern is considered.The load at bus #30 is gradually increased by 15% of its initial load level from the time instant 1 min to 2.5 min of the simulation.
In this test, the voltage magnitude |V 30 | " 0.9388 p.u. of pilot bus #30 drops below the constraint of 0.95 p.u. at t " 2.5 min.This low voltage violation activates the proposed ASVC immediately, and Figure 6a illustrates the voltage trace at pilot bus #30 during the simulation test of case I.When the proposed scheme is used to carry out the secondary voltage control, voltage violation will be eliminated effectively.In addition, the value of VSM is increased from 0.51 to 0.72 p.u., as shown in Figure 6b.This demonstrates that larger VSM can be achieved simultaneously with improvement of overall network voltage.

Case I
In the first case, bus load change pattern is considered.The load at bus #30 is gradually increased by 15% of its initial load level from the time instant 1 min to 2.5 min of the simulation.
In this test, the voltage magnitude p.u. of pilot bus #30 drops below the constraint of 0.95 p.u. at min.This low voltage violation activates the proposed ASVC immediately, and Figure 6a illustrates the voltage trace at pilot bus #30 during the simulation test of case I.When the proposed scheme is used to carry out the secondary voltage control, voltage violation will be eliminated effectively.In addition, the value of VSM is increased from 0.51 to 0.72 p.u., as shown in Figure 6b.This demonstrates that larger VSM can be achieved simultaneously with improvement of overall network voltage.

Case II
In case II, a new operating condition and new load change pattern are investigated.From the time instant 0.5 min to 2 min of the simulation, all the loads in the IEEE 30-bus system are increased by 20% based on the initial load levels.
During the load increase, the voltage magnitudes at the monitored buses (bus #3 and #12 in this case) are within the pre-determined voltage limits, but static voltage stability margin of the power system is less than the threshold level of 0.3 p.u. at 2 t = min.This means that the system has a great potential of voltage collapse without any remedial control.Under such situation, however, no control actions will be initiated by the traditional SVC methods which consider only bus voltages as trigger signal.Instead, the proposed scheme employs not only voltage magnitudes but also static voltage stability margin as criteria; therefore, inadequate VSM ( 0.23 0.3 VSM = < p.u.) is identified.After applying the proposed ASVC, all the voltage magnitudes are maintained within the range of predefined constraint and the static voltage stability margin is enlarged as well.The simulation result is illustrated in Figure 7, showing that the VSM is significantly increased from 0.23 to 0.57 p.u..The enlargement of the static voltage stability margin is a result of the improvement in the system voltage profile.

Case II
In case II, a new operating condition and new load change pattern are investigated.From the time instant 0.5 min to 2 min of the simulation, all the loads in the IEEE 30-bus system are increased by 20% based on the initial load levels.
During the load increase, the voltage magnitudes at the monitored buses (bus #3 and #12 in this case) are within the pre-determined voltage limits, but static voltage stability margin of the power system is less than the threshold level of 0.3 p.u. at t " 2 min.This means that the system has a great potential of voltage collapse without any remedial control.Under such situation, however, no control actions will be initiated by the traditional SVC methods which consider only bus voltages as trigger signal.Instead, the proposed scheme employs not only voltage magnitudes but also static voltage stability margin as criteria; therefore, inadequate VSM (VSM " 0.23 ă 0.3 p.u.) is identified.
After applying the proposed ASVC, all the voltage magnitudes are maintained within the range of predefined constraint and the static voltage stability margin is enlarged as well.The simulation result is illustrated in Figure 7, showing that the VSM is significantly increased from 0.23 to 0.57 p.u..The enlargement of the static voltage stability margin is a result of the improvement in the system voltage profile.

Case III
This example is to demonstrate the performance of the proposed scheme with respect to branch outage contingency following several loads increase condition.In this case, the loads at bus #12, #24, and #29 are increased by 20% from the time instant 0.5 min to 2 min, and the transmission line connected between bus #27 and #28 is tripped in a contingency at 2 min.
As can be expected, the system has poor voltage profile and unsecure margin due to the disturbances of both load increase and line outage.Figure 8 shows the simulation result of case III.One can observe that the proposed scheme works satisfactorily in this case.

Case III
This example is to demonstrate the performance of the proposed scheme with respect to branch outage contingency following several loads increase condition.In this case, the loads at bus #12, #24, and #29 are increased by 20% from the time instant 0.5 min to 2 min, and the transmission line connected between bus #27 and #28 is tripped in a contingency at 2 min.
As can be expected, the system has poor voltage profile and unsecure margin due to the disturbances of both load increase and line outage.Figure 8 shows the simulation result of case III.One can observe that the proposed scheme works satisfactorily in this case.

Case III
This example is to demonstrate the performance of the proposed scheme with respect to branch outage contingency following several loads increase condition.In this case, the loads at bus #12, #24, and #29 are increased by 20% from the time instant 0.5 min to 2 min, and the transmission line connected between bus #27 and #28 is tripped in a contingency at 2 min.
As can be expected, the system has poor voltage profile and unsecure margin due to the disturbances of both load increase and line outage.Figure 8 shows the simulation result of case III.One can observe that the proposed scheme works satisfactorily in this case.

Case IV
In the fourth case, we assume that PMUs are installed at all load buses.In addition, several loads change is considered.The loads at bus #3, #16, and #20 are gradually increased by 15% from the time instant 0.5 min to 1.5 min of the simulation.
Figure 9 shows the traces of the selected pilot buses (bus #23 and #26 in this case) and VSM during the simulation of case IV.Since insufficient VSM ( 0.26 0.3 VSM = < ) is detected at 1.5 t = min, the control strategy will be activated accordingly.Thus, the voltage levels and voltage stability margin are considerably enhanced after ASVC.This can be regarded clearly in Figure 9.

Case IV
In the fourth case, we assume that PMUs are installed at all load buses.In addition, several loads change is considered.The loads at bus #3, #16, and #20 are gradually increased by 15% from the time instant 0.5 min to 1.5 min of the simulation.
Figure 9 shows the traces of the selected pilot buses (bus #23 and #26 in this case) and VSM during the simulation of case IV.Since insufficient VSM (VSM " 0.26 ă 0.3) is detected at t " 1.5 min, the control strategy will be activated accordingly.Thus, the voltage levels and voltage stability margin are considerably enhanced after ASVC.This can be regarded clearly in Figure 9.

Case IV
In the fourth case, we assume that PMUs are installed at all load buses.In addition, several loads change is considered.The loads at bus #3, #16, and #20 are gradually increased by 15% from the time instant 0.5 min to 1.5 min of the simulation.
Figure 9 shows the traces of the selected pilot buses (bus #23 and #26 in this case) and VSM during the simulation of case IV.Since insufficient VSM ( 0.26 0.3 VSM = < ) is detected at 1.5 t = min, the control strategy will be activated accordingly.Thus, the voltage levels and voltage stability margin are considerably enhanced after ASVC.This can be regarded clearly in Figure 9.   Table 2 summarizes the test results for the above three cases.An inspection from the shown table indicates that both voltage profile and static voltage stability margin are greatly improved after applying the proposed control scheme.

Conclusions
A synchrophasor based optimal voltage control scheme, which considers both voltage profile and static voltage stability margin, is developed in order to achieve secure grid operations.A detailed derivation of the principles used in the scheme is presented with illustrated figures.An extensive simulation studies on the IEEE 30-bus test system is carried out to demonstrate the feasibility and effectiveness of the proposed scheme.

Figure 1 .
Figure 1.Schematic view of the proposed scheme.Figure 1.Schematic view of the proposed scheme.

Figure 1 .
Figure 1.Schematic view of the proposed scheme.Figure 1.Schematic view of the proposed scheme.

Figure 2 .
Figure 2. Flowchart of the proposed scheme.

Figure 2 .
Figure 2. Flowchart of the proposed scheme.

Figure 3 .
Figure 3. Block diagram of the proposed ASVC method.Once voltage limit violation or insufficient VSM is identified, the ASVC algorithm will be triggered automatically.

Figure 3 .
Figure 3. Block diagram of the proposed ASVC method.Once voltage limit violation or insufficient VSM is identified, the ASVC algorithm will be triggered automatically.

Figure 4 .
Figure 4. Thevenin equivalent network at bus i with PMU installed.

Figure 4 .
Figure 4. Thevenin equivalent network at bus i with PMU installed.

Figure 5 .
Figure 5. Illustration of the proposed control scheme.

Figure 5 .
Figure 5. Illustration of the proposed control scheme.

Figure 6 .
Figure 6.Simulation result of case I. (a) The trace of 30 V ; (b) The trace of VSM.

Figure 6 .
Figure 6.Simulation result of case I. (a) The trace of |V 30 |; (b) The trace of VSM.

Figure 7 .
Figure 7. Simulation result of case II.(a) The trace of 3 V ; (b) The trace of 12 V ; (c) The trace of VSM.

Figure 7 .
Figure 7. Simulation result of case II.(a) The trace of |V 3 |; (b) The trace of |V 12 |; (c) The trace of VSM.

Figure 7 .
Figure 7. Simulation result of case II.(a) The trace of 3 V ; (b) The trace of 12 V ; (c) The trace of VSM.

Figure 8 .
Figure 8. Simulation result of case III.(a) The trace of 24 V ;(b) The trace of 29 V ; (c) The trace of VSM.

Figure 8 .
Figure 8. Simulation result of case III.(a) The trace of |V 24 |; (b) The trace of |V 29 |; (c) The trace of VSM.

Figure 8 .
Figure 8. Simulation result of case III.(a) The trace of 24 V ;(b) The trace of 29 V ; (c) The trace of VSM.

Figure 9 .
Figure 9. Simulation result of case IV.(a) The trace of |V 23 |; (b) The trace of |V 26 |; (c) The trace of VSM.

Table 2
summarizes the test results for the above three cases.An inspection from the shown table indicates that both voltage profile and static voltage stability margin are greatly improved after applying the proposed control scheme.

Table 2 .
Performance evaluation for the proposed scheme on the IEEE 30-bus system.

Table 2 .
Performance evaluation for the proposed scheme on the IEEE 30-bus system.