Next Article in Journal
Stability Analysis of Stagnation-Point Flow in a Nanofluid over a Stretching/Shrinking Sheet with Second-Order Slip, Soret and Dufour Effects: A Revised Model
Next Article in Special Issue
A Solar–Thermal-Assisted Adiabatic Compressed Air Energy Storage System and Its Efficiency Analysis
Previous Article in Journal
Fusion of Spectra and Texture Data of Hyperspectral Imaging for the Prediction of the Water-Holding Capacity of Fresh Chicken Breast Filets
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Power System Voltage Correction Scheme Based on Adaptive Horizon Model Predictive Control

1
State Grid Shandong Electric Power Research Institute, Jinan 250003, China
2
Key Laboratory of Power System Intelligent Dispatch and Control of Ministry of Education, Shandong University, Jinan 250061, China
3
State Grid Shandong Electric Power Company, Jinan 250001, China
*
Author to whom correspondence should be addressed.
Submission received: 27 February 2018 / Revised: 23 March 2018 / Accepted: 16 April 2018 / Published: 20 April 2018
(This article belongs to the Special Issue Smart Grid and Smart Cities Activities)

Abstract

:
Model predictive control (MPC) is commonly used to compensate for modeling inaccuracies and measurement noise in voltage control problems. The length of the prediction horizon and control horizon of a MPC-based method has significant impact on the control performances. In existing relevant works, those horizon parameters are determined off-line based on experience or enumeration, and keeps constant during the entire receding-horizon optimization process. This paper presents a system voltage correction scheme based on adaptive horizon model predictive control (AH-MPC). The reactive power compensation and voltage regulation devices are coordinated to maintain the system voltages within a desired range. An evaluation index is proposed to determine the horizon parameters, which reflects the maximum voltage regulation ability with the current parameter configuration. Within each sampling interval, the horizon parameters are updated according to the evaluation index and real-time measurements periodically, which comprehensively considers the system uncertainties and voltage recovery speed, and the computational effort is remarkably reduced. The validation and effectiveness of the proposed method is verified by the simulation analysis on the test system.

1. Introduction

The power system voltages should be maintained within predefined limits according to the operation requirements. Nowadays, the increasing penetration of distributed renewable energy has had significant impacts on system behavior due to its variability and uncertainty. Besides, the power systems are operated more closely to their physical limits. The risk of temporary voltage problems or even voltage collapse has increased [1]. When abnormal or unstable voltage profiles are detected, appropriate countermeasures should be taken to enhance voltage quality, power system security and operating economy.
A reliable real-time optimal corrective voltage control scheme requires valid power system prediction models and feedback signals. Due to the inevitable modeling inaccuracies and measurement noise, it is impossible to build an accurate mathematical model to describe all the dynamic features of power systems, thus non-optimal or even inappropriate control decisions may be made. Model predictive control (MPC) is an appealing approach for dealing with the above issues [2], which is especially suitable for complicated systems with changeable dynamics and uncertain factors. The prediction model of a MPC-based method traces the operating state of practical system, and the control sequences are updated and corrected by real-time measurements.
In order to reduce prediction complexity for real-time application, linear prediction models are commonly applied in a MPC-based voltage control scheme. In [3,4], the coordinated voltage control is realized using tree search optimization method based on MPC. The authors of [5,6] apply trajectory sensitivity and MPC to solve the optimal voltage control problem. In [7], the prediction model of a MPC-based real time voltage control scheme is simplified using wide area measurement information in receding horizons. In [8,9], the search solution space of the MPC-based voltage controller is significantly reduced with the prepared control knowledge. A distributed model predictive control scheme for interconnected power systems is presented in [10]. In [11], MPC is applied in a large power system with most effective controls identified.
In recent years, the increasing penetration of renewable power can introduce high variability into transmission and distribution systems, which can be well handled by MPC-based methods. A hierarchical controller based on MPC is proposed to regulate voltages in system with high penetration of wind power in [12]. Model predictive voltage control methods of wind farms are introduced in [13,14,15]. Reference [16] proposes a centralized controller based on MPC to mitigate voltage violations in active distribution networks (ADN) by coordinating power outputs of distributed generation units and the voltage set-point of the transformer load tap changer. The impact of the model predictive control of the ADN voltages to the transmission grid is discussed in [17]. A two-level real-time voltage correction scheme is proposed in [18]. The local level provides the fast response to a disturbance, and the central level coordinates the various distributed generation units relying on the concept of MPC. In [19], the performance of a MPC-based voltage controller is verified by the field test in a real distribution system.
The length of the control horizon and prediction horizon of the existing MPC-based voltage control method has significant impact on the voltage control performances. A larger value of control horizon leads to smoother control but increased computation time. A smaller value of control horizon leads to less control steps and more aggressive control decisions, but it weakens the ability to deal with the modeling inaccuracies and measurement noise problem. The length of the prediction horizon should be equal to the length of the control horizon unless the controller is requested to consider changes happening beyond the control horizon [16]. References [3,4,6,7,8,9,10,11,12,13,14,15,18,19] apply constant horizon parameters, but how those parameters are determined is not discussed. The horizon parameters in [5] is chosen to be the time in which the predicted voltage drops by a predefined percentage at the initial control instant, and the horizon parameters are kept constant during the entire receding-horizon optimization process. In [16,17], the controller uses constant horizon parameters unless the prediction horizon is increased to anticipate load tap changer actions. The control methods discussed so far determines the control horizon parameters off-line based on experience or enumeration without quantitative analysis.
In this paper, a coordinated voltage correction scheme based on adaptive horizon model predictive control (AH-MPC) is proposed to maintain a satisfactory voltage profile while minimizing the control cost. Generator voltage references, shunt capacitors and on-load tap changers are considered as candidate control actions. The proposed scheme is based on the principles of MPC. However, it differs from the standard MPC-based method by way the horizon parameters are chosen and embedded in the optimization procedure. The main contribution of this paper is developing a principle of determine the horizon parameters, which is a key feature missing in previous control schemes. An index is formulated to evaluate and ensure the applicability of the horizon parameter configurations, and barely consume extra computation time. Within each sampling interval, the length of the horizon parameters are updated according to the proposed evaluation index and real-time measurements, aiming at better control performances compared with the traditional MPC-based methods.
This paper is organized as follows: Section 2 introduces the standard MPC-based voltage control model. Section 3 presents the determination of the adaptive horizon parameters. The voltage correction scheme based on AH-MPC is presented in Section 4. The simulation results and conclusions are presented in Section 5 and Section 6 respectively.

2. Voltage Control Model Based on MPC

2.1. Voltage Prediction

MPC finds a sequence of control actions applied within the control horizon in order to achieve the control objectives within the prediction horizon considering equality and inequality constraints. At each sampling point, only the first step of the so-computed sequence is applied. The principle of MPC is shown in Figure 1. The length of the control horizon and prediction horizon are tc and tp respectively. ts is the sampling interval. Nc and Np are the number of sampling intervals within the control horizon and prediction horizon respectively, tc = Ncts, tp = Npts.
Real-time optimization requires reliable and quick decisions. A linear model based on trajectory analysis is applied to predict the trajectory of the system with different control sequences. The discrepancies arising from the model mismatch problems and system uncertainties can be resolved by the closed-loop feature of the centralized MPC-based controller [2]. The trajectory sensitivity analysis can provide valuable insights evaluating the effects of slight control variations on the system evolution. The detailed theory and analytical basis of trajectory sensitivity is discussed in [20].
The system model can be expressed by the differential-algebraic Equations (DAE) (1) and (2) subject to some constraints, the detailed introduction of the DAE model can be found in [1].
x ˙ = f ( x , y , u )
0 = g ( x , y , u )
where x is the vector of state variables, which is related to the generator and its regulation system, the load recovery behavior et al., y is the vector of algebraic variables such as active and reactive powers P and Q, voltage magnitudes V and phase angles θ, y = [P, Q, V, θ]. u is the vector of control variables. In this paper, the set points of automatic voltage regulators (AVR), the switching of shunt capacitors and the position of on-load tap changers (OLTCs) are considered as candidate control actions. The sensitivity Equations (3) and (4) are obtained by differentiating (1) and (2) with respect to u.
x ˙ u ( t ) = f x ( t ) x u ( t ) + f y ( t ) y u ( t ) + f u ( t )
0 = g x ( t ) x u ( t ) + g y ( t ) y u ( t ) + g u ( t )
where yu(t) is composed of the trajectory sensitivity of all the algebraic variables with respect to the candidate control variables at a time instant t, which can be obtained by solving (3) and (4). yu(t) = [Pu(t), Qu(t), Vu(t), θu(t)]. When the time domain simulation of the system (1) and (2) is based on the trapezoidal numerical integration, the Jacobian obtained during the numerical integration process is the same as the coefficient matrix of solving the sensitivity Equations (3) and (4). When the result of time domain simulation is obtained, the calculation of trajectory sensitivity involves little extra computation [20]. When control action Δ u is applied, the corresponding voltage variation at a sampling point k can be approximately expressed as Δ V ^ k = V u ( t n + k t s ) Δ u , V u ( t n + k t s ) is the voltage trajectory sensitivity at sampling point k with respect to control u, and also one of the elements that constitute yu(t).
The voltage prediction based on trajectory sensitivity analysis is graphically illustrated in Figure 2. According to the proportion and superposition characteristics of linear systems, the predicted voltage at a sampling point k is composed of two parts: the voltage without control V ^ k * during the prediction horizon and the voltage variation Δ V ^ k cause by the control actions applied at each sampling points. At the initial control time, V ^ k * is solved based on trapezoidal numerical integration of (1) and (2). Δ V ^ k is solved based on the Jacobian obtained during the numerical integration process. The vector consist of all the predicted voltages at each sampling point can be expressed as:
V ^ k = V ^ k * + Δ V ^ k

2.2. Optimization Model

The proposed scheme aims at minimizing the system voltage deviation and the corresponding control cost. The optimization problem can be formulated as follows:
min ( k = 1 N p ( V r V ^ k ) Q ( V r V ^ k ) T + k = 1 N c Δ u k R Δ u k T )
Subject to:
u k = u 0 + l = 1 k Δ u l
V min V ^ N p V max
u min u k u max
| Δ u k | Δ u m a x
where Vr is the vector consists of the control references. V ^ k is the predicted voltage at sampling point k. Δ u k is the vector consists of the control variable changes at sampling point k. Q and R are both diagonal matrices which define the weight coefficients of the controlled voltages and the cost of the control actions respectively.
The voltage constraints in the prediction horizon is shown in Equation (8), V ^ N p represents the predicted voltage at the sampling point Np. The objective function sets a target of minimizing the voltage deviations at every sampling point under the constraints of control input (9) and (10), and the voltage constraints are only enforced at the end of the prediction horizon. The inequality constraints (9) stipulates the upper and lower bounds of the control inputs.
When undesirable bus voltage magnitudes are detected, appropriate voltage control actions should be applied as soon as possible to reduce the negative influence to the consumers and avoid the deterioration of system voltages. The single-step optimization approach or large control variations in one control step may suffers from deficiencies such as low optimization efficiency or undesirable control results [21]. Thus the constraint of control variation in one control step (10) is included in the model.

3. Adaptive Horizon Parameters

MPC is a receding-horizon multi-step optimization method. The computation time and control performance are closely related to the horizon parameter selection, which determines when and how the system variables will reach their target states. This section presents the determination of adaptive horizon parameters, which are periodically updated according to the real-time measurements and system trajectory prediction during the entire optimization process of MPC.

3.1. Principles for the Horizon Parameter Settings

At each sampling point, the MPC-based method obtains an optimal control sequence [u1, u2, …, uNc] distributed within the control horizon. The control variables of the control sequence keeps constant between sampling intervals, and only changes at sampling points within the control horizon. The basic principles for the horizon parameters are listed as follows:
  • Sampling interval ts. The length of sampling interval ts should fulfill the requirements of measurement collection, data processing and optimization, which change little during the entire receding horizon optimization process. A constant value ts = 10 s is applied in the proposed scheme, which also matches the value used in the coordinated secondary voltage control in operation in practical systems [22].
  • Control horizon tc. With a constant sampling interval, the optimal length of the control horizon is related to the voltage evolution and response feature, which varies over the receding horizon optimization process. In this paper, the system voltage is predicted using a linear model as is introduced in Section 2. The sensitivity feature changes little if no disturbances occur during the optimization process [21], and the solved control variation is proportional to the detected voltage deviation. When severe voltage deviation is detected, the length of the control horizons should be short enough to meet the requirement of the voltage recovery speed and online computation time, and should be long enough to deal with modeling inaccuracies and measurement noise. As the receding horizon optimization progressed, the controlled voltages are gradually closer to its desired range, and smaller control variations are needed to accomplish the control objective. In this situation, shorter control horizon parameters could reduce settling times of the optimization process and computational effort.
  • Prediction horizon tp. The system response caused by the control actions are predicted within the prediction horizon, thus the length of the prediction horizon should be greater or equal to the length of the control horizon. When the controlled voltages are gradually stabilized to the steady state values, the necessity to consider changes happening beyond the control horizon is reduced from a computational viewpoint. Therefore, at the beginning of the MPC implementation, the length of the prediction horizon is set longer than the length of control horizon. As the receding horizon optimization progressed, the length of the prediction horizon gradually approaches the length of the control horizon, which could also save unnecessary computation time

3.2. Voltage Regulation Ability Evaluation

Under the condition of fulfilling the control variation constraint (10), large control variations are forced to be applied in multi-steps. The limited candidate control actions in one control step have theoretical maximum voltage regulation ability for every bus. The procedure to estimate the ability is summarized as follows:
(1)
At tn, the voltage evolution trajectories are predicted within a predetermined prediction horizon [tn, tn + tpn] based on the trapezoidal numerical simulation. tpn is the length of prediction horizon applied in the n-th optimization. If the predicted voltage deviation of load bus i satisfies the inequality constraint (11) at the initial time instant t1, the voltage magnitude of bus i is determined as the evaluation target.
| V r V ^ ( i ) * ( t 1 + t p 1 ) | D B / 2
where V ^ ( i ) * ( t 1 + t p 1 ) is the voltage magnitude of bus i without control at the prediction time instant (t1 + tp1) obtained by numerical simulation. DB is the control deadband. Vmax and Vmin initially presented in Equation (8) are subject to: Vmax = Vr + DB/2, Vmin = VrDB/2.
(2)
The trajectory sensitivity matrix of the evaluation targets at time instant tn + tpn with respect to the control action of each step can be expressed as:
S k = [ S k ( 1 , 1 ) ( t n + t p n ) S k ( 1 , 2 ) ( t n + t p n ) S k ( 1 , N j ) ( t n + t p n ) S k ( 2 , 1 ) ( t n + t p n ) S k ( 2 , 2 ) ( t n + t p n ) S k ( 2 , N j ) ( t n + t p n ) S k ( N i , 1 ) ( t n + t p n ) S k ( N i , 2 ) ( t n + t p n ) S k ( N i , N j ) ( t n + t p n ) ]
where Sk(i, j) represents the trajectory sensitivity of the i-th evaluation target with respect to the j-th candidate control in the k-th control step. Ni and Nj are the total number of evaluation targets and available control actions respectively.
(3)
Based on the voltage prediction and trajectory sensitivity analysis, the maximum control target variation caused by controls applied at sampling point k is evaluated as follows.
When S k ( i , j ) ( V r V ^ ( i ) * ( t n + t p n ) ) 0 , the optimization tends to increase the j-th control in order to raise (when S k ( i , j ) 0 and V r V ^ ( i ) * ( t n + t p n ) ) or reduce (when S k ( i , j ) 0 and V r V ^ ( i ) * ( t n + t p n ) ) the voltage magnitude of bus i towards its desired value. Under the inequality constraints (9) and (10), the predicted voltage magnitude variation of bus i reaches its maximum value when the variation of control j equals its maximum positive value:
Δ V k ( i , j ) max ( t n + t p n ) = { S k ( i , j ) ( t n + t p n ) Δ u ( j ) max      i f      u k ( j ) + Δ u ( j ) max u ( j ) max S k ( i , j ) ( t n + t p n ) ( u ( j ) max u k ( j ) )      i f      u k ( j ) + Δ u ( j ) max u ( j ) max
where Δ V k ( i , j ) max ( t n + t p n ) is the predicted maximum variation of the evaluation target i at a future time instant tn + tpn, which is caused by the j-th control action applied at the sampling point k. Δ u ( j ) max is the j-th element in vector Δ u m a x , representing the maximum value of the corresponding control variation. u ( j ) max is the j-th element in vector umax.
Similarly, when S k ( i , j ) ( V r V ^ ( i ) * ( t n + t p n ) ) < 0 , the optimization tends to reduce the j-th control in order to reduce (when S k ( i , j ) > 0 and V r < V ^ ( i ) * ( t n + t p n ) ) or raise (when S k ( i , j ) < 0 and V r > V ^ ( i ) * ( t n + t p n ) ) the voltage magnitude of bus i towards its desired value. The maximum voltage variation value of bus i:
Δ V k ( i , j ) max ( t n + t p n ) = { S k ( i , j ) ( t n + t p n ) ( Δ u ( j ) max )      i f      u k ( j ) Δ u ( j ) max u ( j ) min S k ( i , j ) ( t n + t p n ) ( u ( j ) min u k ( j ) )      i f      u k ( j ) Δ u ( j ) max u ( j ) min
where u ( j ) min is the j-th element in vector umin.
The maximum voltage variation of the i-th evaluation target considering all candidate control actions applied at sampling point k can be expressed as:
Δ V k ( i ) max ( t n + t p n ) = j = 1 N j Δ V k ( i , j ) max ( t n + t p n )
(4)
Considering multi control steps applied at the first m sampling points, The maximum voltage regulation ability at time instant tn + tpn caused by the control sequence ( Δ u 1 , Δ u 2 , …, Δ u m ) is evaluated as follows:
V a ( i ) m = V ^ ( i ) * ( t n + t p n ) + k = 1 m Δ V k ( i ) max ( t n + t p n )
where Va(i)m is the evaluated voltage regulation ability of the i-th evaluation target when the length of control horizon is mts.

3.3. Determination of the Horizon Parameters

Based on the voltage regulation ability evaluation, an applicability judgment index for control horizon parameter settings is proposed:
I s ( i ) m = D B 2 + V r V ^ ( i ) * ( t n + t p n ) | V r V ^ ( i ) * ( t n + t p n ) | ( V a ( i ) m V r )
where Is(i)m is the applicability judgment index of the i-th evaluation target when the control horizon is set to mts. It should be noted that the proposed scheme will not be triggered if V r = V ^ ( i ) * ( t n + t p n ) .
When the predicted voltage without control is lower than its reference value ( V r > V ^ ( i ) * ( t n + t p n ) ), Equation (13) can be written as:
I s ( i ) m = V a ( i ) m ( V r D B 2 ) = V a ( i ) m V min
Similarly, when V r < V ^ ( i ) * ( t n + t p n ) , Equation (13) can be written as:
I s ( i ) m = ( V r + D B 2 ) V a ( i ) m = V max V a ( i ) m
If I s ( i ) m > 0 , it indicates that the controller is able to bring the voltage magnitude of bus i back to its desired range with the current horizon parameter settings. Let Ism be a vector that consists the applicability judgment index of all the evaluation targets, Ism = (Is(1)m, Is(2)m, ……, Is(Ni)m). The initial value of m is set to 1. If ∀Is(i)m ≤ 0 in vector Ism, the proposed algorithm will keep increasing m until all the elements in Ism meet Is(i)m > 0, then mts is determined as the control horizon of the n-th optimization and the prediction horizon of the (n + 1)-th optimization, mts = tcn = tp(n+1).
If it is evaluated that the candidate control variables cannot achieve ∀Is(i)m > 0 according to Equations (13)–(16), m will be reset to 1, and control actions with high cost such as load shedding will be added to the candidate control actions. In this paper, it is assumed that the power system is equipped with adequate reactive power compensation and voltage regulation devices for a voltage correction scenario.
As is presented in Figure 3, the length of the control horizon and prediction horizon are adaptively adjusted during the receding-horizon optimization process.

4. Voltage Correction Scheme Based on AH-MPC

The structure of the proposed scheme is shown in Figure 4. The system measurements are collected and processed by a state estimator in the control center. The results are then used to initialize the state variable of the system model and predict the future system trajectory. Compared with a standard MPC controller, the byproduct of the voltage prediction results are used to update the horizon parameter settings. With the numerical simulation results of V ^ * and the trajectory sensitivity information already available, the calculation of the horizon parameters based on Equations (13)–(17) requires solving a set of linear algebraic equations, which involves little extra computation.
With adaptive horizons applied, the objective function (6) and constraint (8) in the standard MPC-based optimization model are rewritten as:
min ( k = 1 N p n ( V r V ^ k ) Q ( V r V ^ k ) T + k = 1 N c n Δ u k R Δ u k T )
V min V ^ N p n V max
N p ( n + 1 ) = N c n = t c n / t s
where Npn and Ncn are the number of prediction steps and control steps applied in the n-th optimization.
Figure 5 shows the flow chart for the voltage correction scheme based on AH-MPC. The procedure of the proposed scheme in the n-th optimization include the following.
  • At a sampling point tn, trapezoidal numerical simulation is used to obtain the predicted system trajectory within the prediction horizon [tn, tn + tpn].
  • The trajectory sensitivity is obtained as a by-product of the first step.
  • Based on the predicted system trajectory and the trajectory sensitivity results, the adaptive horizon parameters are determined as is introduced in Section 3.
  • The optimal control model is formulated according to Equations (7), (9), (10), (20), (21) and (22).
  • By solving the optimal control model, the first step of the so computed control sequence is applied in the system at tn + ts.
  • The above procedure is repeated at the initial time of the (n + 1)-th optimization.
If disturbances are detected by phase measurement unit (PMU) measurements and monitoring systems, all the parameters are initialized at the next sampling point.

5. Simulations

The proposed method is applied to the nine-bus system and the New England test system. Dynamic recovery load models [23,24] are used in the simulation, with αs = βs = 0, αt = βt = 2 and Tp = Tq = 60. The voltage reference Vr = 1, and the voltage control deadband is 0.1 p.u., which indicates that the voltage magnitude should be kept within [0.95 p.u., 1.05 p.u.] with control. The sampling interval ts is 10 s. The initial prediction horizon tp1 is 60 seconds. The diagonal elements of Q and R are set to be 1 because the correction of each bus voltage is equally important and no expensive control actions such as load shedding are considered in the simulation cases. The trajectory prediction and adaptive horizon parameter determination is realized based on the quasi-steady state model [1], and the actual control performance is illustrated by the full-time domain simulation on a detailed system model.
The model prediction was realized with the Matlab/Simulink-based tool. The optimization problem was solved by GAMS [25]. The calculations were completed using a windows machine with Intel Core 2.50 GHz CPU and 8 GB of RAM.

5.1. Case A: Nine-Bus Test System

Simulation on a simple test system is presented in this section to demonstrate the basic idea and implementation procedure of the proposed scheme. The topology of the test system is shown in Figure 6. The available candidate control actions include the shunt capacitors at bus6, the AVR set points of generator at bus2 and bus3, the OLTCs of the transformer between bus9 and bus3. The constraint of control variation in one control step for the shunt capacitors, AVR set points and the OLTCs are 0.1 p.u., 0.04 p.u. and 1.67% respectively.
Line between bus 4 and bus 6 are tripped at t = 10 s. The evolution of voltages at load bus 5, 6 and 8 are illustrated in Figure 7. Due to the loss of the transmission system, the voltage magnitudes drops below the lower bound if no control actions are applied.
The voltage evolution with the proposed control scheme is shown in Figure 8. The step-like voltage variations are caused by the discrete control feature of the reactive power compensation and voltage regulation devices, which is accordant with the voltage response feature in real-life applications. At the initial time instant t1 = 10 s, the proposed controller predicts the system trajectory over a predefined prediction horizon tp1 = 60 s. The voltage magnitude of load bus 6 and bus 8 fall out of the lower limit 0.95 p.u at the end of the prediction horizon t1 + tp1 = 70 s, thus the voltage magnitudes of these two buses are determined as the evaluation targets. The horizon parameters are updated at the initial time of each optimization according to the applicability judgment index Ism. The value of Ism and the applied horizon parameters in this case is listed in Table 1, where n represents the n-th optimization.
When n = 1 and m = 1, the applicability judgment index Is1 is obtain according to Equation (17). As is shown in Table 1, the value of the element Is(1)1 = −0.00195 < 0, thus the proposed controller increases m and calculate Is2. Since all the elements of Is2 are positive, the control horizon of the current optimization tc1 and the prediction horizon of the next optimization tp2 are both determined as 20 s. The control model of the current optimization is then formulated and solved. Considering the time consuming of communication and calculation, the above procedure is started at 10 s, and the first step of the so-computed control sequence is applied a sampling interval later. When n = 2 and m = 1, all the elements of Is1 are positive, thus tc2 = tp3 = 10 s. According to the system trajectory prediction implemented during the third optimization, the voltage magnitudes of all the system buses has restored to the desired range and stops declining, the receding-horizon optimization is then completed. The adaptive horizons applied in the simulation is graphically illustrated in Figure 9. The control actions are listed in Table 2, where Uref,Gi represents the AVR set point of generator i, Ci is the applied shunt capacitor at bus i, nt(ij) is the ratio of the transformer connecting bus i and j.

5.2. Case B: New England Test System

In this case, the validation and effectiveness of the proposed method further verified by the simulation on the New England 39-bus test system. Figure 10 shows the topology of the test system. The candidate control actions include the shunt capacitors at bus 8, the AVR set points of all the 10 generators, the OLTCs of the transformers connecting bus12 and 11, 12 and 13, 19 and 20. The control parameter settings are the same as case A.
At t = 10 s, generator 3 are tripped. As is shown in Figure 11.The voltage magnitude of load bus 4, 7, 8 and 12 drops below 0.95 p.u. after the fault occurs.
According to Equation (11), the voltage magnitudes of load bus 4, 7, 8 and 12 are determined as the evaluation targets. Figure 12 presents the corrected voltage trajectory with the proposed scheme started at t = 10 s. The applied control actions are listed in Table 3. As the receding horizon optimization progresses, the system voltages are gradually brought back inside the limits at t = 40 s. The applicability judgment index Ism is listed in Table 4. When large voltage deviation of the control targets is detected at the initial control instant t = 10 s, the proposed scheme applies three-step optimization according to the value of Is3. As the voltages gradually approaching the desired range, the declining voltage trajectories are stabilized, and smaller control variations are needed to complete the optimization process. In the third optimization, the control variation needed is small enough to be applied in one step without causing negative effect listed in [21], and shorter prediction length beyond the control horizon is needed due to the steadier voltage evolution trend, which reduces the computation time and settling times of the entire optimization process. The applied decreasing prediction horizon and control horizon are graphically illustrated in Figure 13.

5.3. Simulating Modeling Inaccuracies

In this section, the modeling inaccuracies and disturbance are considered in the course of controlling the system. The fault scenario, system parameter and other controller parameter settings are the same as case B. At t = 25 s, A disturbance is imposed to the load at bus 3, bus 4, bus 6 and bus7 with 5% increase of the active and reactive power.
If the disturbance can be measured and sent to the control center, the voltage behavior based on AH-MPC is shown in Figure 14. The system voltage drops dramatically due to the disturbance at t = 25 s. At the next sampling point t = 30 s, the optimization is carried out based on the updated system measurements.
If the disturbance is ignored by the control center, there will be modeling inaccuracy problems at all future control steps after t = 25 s. The voltage behavior based on AH-MPC is shown in Figure 15. Although the prediction model is not completely updated after the disturbance, the controller will get the actual voltage responses according to measured feedback signals at each future sampling point, and the structure of receding-horizon multi-step optimization offers chance to correct the consequences of modeling errors.
From Figure 14 and Figure 15, we can see that the proposed control scheme is still valid. This is mainly because the disturbance and modeling inaccuracy problems are compensated by the closed-loop feature of the centralized controller.

5.4. Control Performance Comparison

The control performance of the proposed method is compared with the ordinary MPC-based method using constant horizon parameters (tc = 30 s, tp = 60 s). The fault scenario, system parameter and other controller parameter settings are the same as case B. The control result is shown in Figure 16. All the voltages meet the constraint requirements after 4 optimizations are implemented. The performance comparison is listed in Table 5. An index of average voltage deviation of all the load buses is introduced as:
Δ V offset = 1 19 i = 1 19 ( 1 Δ t t 1 t 1 + Δ t | V r V ^ i ( t ) | d t )
The integration time Δ t = 300 s. Smaller value of Δ V offset indicates faster voltage recovery speed and better control performance.
The comparison results of Figure 12 and Figure 16 and Table 5 show that both methods are able to stabilize the system voltage within a desire range, and the AH-MPC-based method outperforms MPC-based scheme with constant horizon in computation effort and voltage recovery speed.
The proposed scheme applies decreasing prediction parameters and control horizon. When n = 1, both methods applies the same horizon parameters and cost almost the same computation time. With the increase of n, the time consuming per optimization of the AH-MPC-based scheme is remarkably reduced. This is mainly because the numerical simulation needed in the model prediction is implemented within a shorter prediction horizon, and the control sequence consists of less control steps is computed. Moreover, the proposed scheme takes less settling times, thus there are less optimal control model to be formulated and solved, which also leads to decreased computational effort.

5.5. Discussion

A principle of determine the horizon parameters based on measurements and prediction is developed in this paper, which is integrated with the standard MPC-based method as is graphically illustrated in Figure 4. The determination of horizon parameters costs little computation effort, and the performance of the controller can be improved, which has been verified by the simulations. It should be noted that the basic idea of adaptive horizon parameter settings is not restricted to the method presented in Section 2, but can be combined with any MPC-based method using linear prediction model.
The time-consuming nature of the optimization is affected by many aspects, such as hardware configuration, system scale and prediction method. Simulation results of a load increase scenario in IEEE 118 bus test system [26] is presented to further verify the practicability of the proposed method. The active and reactive powers of loads at buses 82, 83, 84, 85 and 86 are linearly increased at a rate of 10% per second until t = 20 s. The voltage trajectories of bus 83 with and without control are shown in Figure 17. With the actions of the AVR set points of generators at bus 77, 85, 87, 89 and 90, the system voltages are kept within a desired range. The maximum computation time for one optimization (when n = 1) is 5.18 s.
The measurement collection and optimization starts at each sampling point, and the control decisions are applied to the practical system a sampling interval (ts = 10 s) later, which offers sufficient time for the telecommunication and computation. In a real-life application, it is confidently expected that the computation time can be further reduced:
  • The simulations were carried out on a personal computer with the Matlab/Simulink-based tool. Both the hardware and software configuration can be significantly improved in a real-life application [27].
  • The optimization problem can be solved concentrating on areas close to the disturbance with a simple representation used for other areas [21].

6. Conclusions

This paper proposes a voltage correction scheme based on AH-MPC. By coordinating the control actions of the AVR set points, OLTCs and shunt capacitors, the system voltages are kept within a desired range after the disturbance occurs. The inevitable prediction model mismatch problem is compensated by the closed-loop feature of the centralized controller. This paper also presents a real-time determination method for horizon parameters which is commonly ignored in the related research. The horizon parameters are updated at each sampling point according to the real-time feedback signals and system prediction results, which effectively improves the control performance compared with the traditional MPC-based method. The information needed to determine the horizon parameters can be obtained as by-products of the prediction process in a standard MPC-based method, thus little extra computation time is needed.

Author Contributions

Yan Zhang and Meng Liu designed the scheme and wrote the paper; Wen Zhang offered constructive suggestions for the algorithm design; Wenchuan Sun, Xingwang Hu and Gang Kong helped perform the simulations and analyze the results.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Van Cutsem, T.; Vournas, C. Voltage Stability of Electric Power Systems; Springer: Boston, MA, USA, 1998; pp. 193–210. [Google Scholar]
  2. Rawlings, J.B. Tutorial overview of model predictive control. IEEE Control Syst. Mag. 2000, 20, 38–52. [Google Scholar] [CrossRef]
  3. Larsson, M.; Hill, D.J.; Olsson, G. Emergency voltage control using search and predictive control. Int. J. Electr. Power Energy Syst. 2002, 24, 121–130. [Google Scholar] [CrossRef]
  4. Larsson, M.; Karlsson, D. Coordinated system protection scheme against voltage collapse using heuristic search and predictive control. IEEE Trans. Power Syst. 2003, 18, 1001–1006. [Google Scholar] [CrossRef]
  5. Jin, L.; Kumar, R.; Elia, N. Model predictive control-based real-time power system protection schemes. IEEE Trans. Power Syst. 2010, 25, 988–998. [Google Scholar] [CrossRef]
  6. Hiskens, I.A.; Gong, B. Voltage stability enhancement via model predictive control of load. Intell. Autom. Soft Comput. 2006, 12, 117–124. [Google Scholar] [CrossRef]
  7. Zhang, Y.; Zhang, W.; Chu, X.; Liu, Y. Real-time optimal voltage control using measurement-based aggregate load model. Electr. Power Syst. Res. 2014, 116, 293–300. [Google Scholar] [CrossRef]
  8. Ma, H.; Hill, D.J. Adaptive coordinated voltage control-Part I: Basic scheme. IEEE Trans. Power Syst. 2014, 29, 1546–1553. [Google Scholar] [CrossRef]
  9. Ma, H.; Hill, D.J. Adaptive coordinated voltage control—Part II: Use of learning for rapid response. IEEE Trans. Power Syst. 2014, 29, 1554–1561. [Google Scholar] [CrossRef]
  10. Moradzadeh, M.; Boel, R.; Vandevelde, L. Anticipating and coordinating voltage control for interconnected power systems. Energies 2014, 7, 1027–1047. [Google Scholar] [CrossRef] [Green Version]
  11. Martin, J.A.; Hiskens, I.A. Corrective model-predictive control in large electric power systems. IEEE Trans. Power Syst. 2017, 32, 1651–1662. [Google Scholar]
  12. Tang, Z.; Hill, D.J.; Liu, T.; Ma, H. Hierarchical voltage control of weak subtransmission networks with high penetration of wind power. IEEE Trans. Power Syst. 2018, 33, 187–197. [Google Scholar] [CrossRef]
  13. Zhao, H.; Wu, Q.; Wang, J.; Liu, Z.; Shahidehpour, M.; Xue, Y. Combined active and reactive power control of wind farms based on model predictive control. IEEE Trans. Energy Convers. 2017, 32, 1177–1186. [Google Scholar] [CrossRef]
  14. Zhao, H.; Wu, Q.; Guo, Q.; Sun, H.; Huang, S.; Xue, Y. Coordinated Voltage Control of a Wind Farm Based on Model Predictive Control. IEEE Trans. Sustain. Energy. 2016, 7, 1440–1451. [Google Scholar] [CrossRef]
  15. Guo, Y.; Gao, H.; Wu, Q.; Zhao, H.; Østergaard, J. Coordinated voltage control scheme for VSC-HVDC connected wind power plants. IET Renew. Power Gener. 2018, 12, 198–206. [Google Scholar] [CrossRef]
  16. Valverde, G.; Van Cutsem, T. Model predictive control of voltages in active distribution networks. IEEE Trans. Smart Grid. 2013, 4, 2152–2161. [Google Scholar] [CrossRef]
  17. Aristidou, T.; Valverde, G.; Van Cutsem, T. Contribution of distribution network control to voltage stability: A case study. IEEE Trans. Smart Grid. 2017, 8, 106–116. [Google Scholar] [CrossRef]
  18. Soleimani Bidgoli, H.; Van Cutsem, T. Combined Local and Centralized Voltage Control in Active Distribution Networks. IEEE Trans. Power Syst. 2018, 33, 1374–1384. [Google Scholar] [CrossRef]
  19. Balram, P.; Carlson, O.; Tuan, L.A. Demonstration of voltage control in a real distribution system using model predictive control. IET Gener. Transm. Distrib. 2017, 11, 3922–3929. [Google Scholar] [CrossRef]
  20. Hiskens, I.A.; Pai, M. Trajectory sensitivity analysis of hybrid systems. IEEE Trans. Power Syst. 2000, 47, 204–220. [Google Scholar] [CrossRef]
  21. Glavic, M.; Hajian, M.; Rosehart, W.; Van Cutsem, T. Receding-horizon multi-step optimization to correct nonviable or unstable transmission voltages. IEEE Trans. Power Syst. 2011, 26, 1641–1650. [Google Scholar] [CrossRef]
  22. Vu, H.; Pruvot, P.; Launay, C.; Harmand, Y. An improved voltage control on large-scale power system. IEEE Trans. Power Syst. 1996, 11, 1295–1303. [Google Scholar] [CrossRef]
  23. Hill, D.J. Nonlinear dynamic load models with recovery for voltage stability studies. IEEE Trans. Power Syst. 1993, 8, 166–176. [Google Scholar] [CrossRef]
  24. Karlsson, D.; Hill, D.J. Modelling and identification of nonlinear dynamic loads in power systems. IEEE Trans. Power Syst. 1994, 9, 157–166. [Google Scholar] [CrossRef]
  25. Rosenthal, R.E. GAMS—A User’s Guide; GAMS Development Corporation: Washington, DC, USA, 2010; pp. 1–3. [Google Scholar]
  26. Power System Test Case Archive. Available online: http://www.ee.washington.edu/research/pstca/ (accessed on 20 April 2018).
  27. Wen, J.; Wu, Q.; Turner, D.; Cheng, S.; Fitch, J. Optimal coordinated voltage control for power system voltage stability. IEEE Trans. Power Syst. 2004, 19, 1115–1122. [Google Scholar] [CrossRef]
Figure 1. The principle of model predictive control.
Figure 1. The principle of model predictive control.
Applsci 08 00641 g001
Figure 2. Voltage prediction based on trajectory sensitivity analysis.
Figure 2. Voltage prediction based on trajectory sensitivity analysis.
Applsci 08 00641 g002
Figure 3. Determination of adaptive horizon parameters.
Figure 3. Determination of adaptive horizon parameters.
Applsci 08 00641 g003
Figure 4. Structure of the voltage correction scheme based on adaptive horizon model predictive control (AH-MPC).
Figure 4. Structure of the voltage correction scheme based on adaptive horizon model predictive control (AH-MPC).
Applsci 08 00641 g004
Figure 5. Flow chart for the voltage correction scheme based on AH-MPC.
Figure 5. Flow chart for the voltage correction scheme based on AH-MPC.
Applsci 08 00641 g005
Figure 6. Nine-bus test system.
Figure 6. Nine-bus test system.
Applsci 08 00641 g006
Figure 7. Voltage evolution without control in case A.
Figure 7. Voltage evolution without control in case A.
Applsci 08 00641 g007
Figure 8. Voltage evolution using AH-MPC in case A.
Figure 8. Voltage evolution using AH-MPC in case A.
Applsci 08 00641 g008
Figure 9. Adaptive horizons applied in case A.
Figure 9. Adaptive horizons applied in case A.
Applsci 08 00641 g009
Figure 10. New England 39-bus test system.
Figure 10. New England 39-bus test system.
Applsci 08 00641 g010
Figure 11. Voltage evolution without control in case B.
Figure 11. Voltage evolution without control in case B.
Applsci 08 00641 g011
Figure 12. Voltage evolution using AH-MPC in case B.
Figure 12. Voltage evolution using AH-MPC in case B.
Applsci 08 00641 g012
Figure 13. Adaptive horizons applied in case B.
Figure 13. Adaptive horizons applied in case B.
Applsci 08 00641 g013
Figure 14. Voltage evolution using AH-MPC with disturbance occurred during the control process.
Figure 14. Voltage evolution using AH-MPC with disturbance occurred during the control process.
Applsci 08 00641 g014
Figure 15. Voltage evolution using AH-MPC with both disturbance and model mismatch occurred during the control process.
Figure 15. Voltage evolution using AH-MPC with both disturbance and model mismatch occurred during the control process.
Applsci 08 00641 g015
Figure 16. Voltage evolution with control in case B using constant horizon parameters.
Figure 16. Voltage evolution with control in case B using constant horizon parameters.
Applsci 08 00641 g016
Figure 17. Voltage evolution of bus 83 in IEEE 118 test system using AH-MPC.
Figure 17. Voltage evolution of bus 83 in IEEE 118 test system using AH-MPC.
Applsci 08 00641 g017
Table 1. The applicability judgment index Ism in case A.
Table 1. The applicability judgment index Ism in case A.
nIs1 (p.u.)Is2 (p.u.)Control Horizon (s)Prediction Horizon (s)
1(−0.0195, 0.0128)(0.0298, 0.0483)2060
2(0.0285, 0.0429) 1020
Table 2. Control actions in case A.
Table 2. Control actions in case A.
Time (s)Control
20Uref,G2 = 1.042 p.u., Uref,G3 = 1.033 p.u., C6 = 0.1 p.u., nt(12-11) = 1.0167
30Uref,G2 = 1.064, C6 = 0.2 p.u.
Table 3. Control actions in case B.
Table 3. Control actions in case B.
Time (s)Control
20Uref,G1 = 1.062 p.u., Uref,G2 = 1.025 p.u., Uref,G4 = 1.007 p.u., Uref,G5 = 1.036 p.u., Uref,G6 = 1.051 p.u., Uref,G7 = 1.069 p.u.,
Uref,G8 = 1.017 p.u., Uref,G9 = 1.025p.u., Uref,G10 = 1.023 p.u., C8 = 0.1 p.u., nt(12-11) = 1.0227, nt(12-13) = 1.0227, nt(19-20) = 1.0767
30Uref,G1 = 1.102 p.u., Uref,G2 = 1.065 p.u., Uref,G4 = 1.001 p.u., Uref,G5 = 1.017 p.u., Uref,G6 = 1.029 p.u., Uref,G7 = 1.061 p.u.,
Uref,G8 = 1.021 p.u., Uref,G9 = 1.002 p.u., Uref,G10 = 1.055 p.u., C8 = 0.2 p.u., nt(12-11) = 1.0394, nt(12-13) = 1.0227, nt(19-20) = 1.0934
40Uref,G1 = 1.128 p.u., Uref,G2 = 1.084 p.u., Uref,G4 = 0.998 p.u., Uref,G5 = 1.010 p.u., Uref,G6 = 1.004 p.u., Uref,G7 = 1.054 p.u.,
Uref,G8 = 1.031 p.u., Uref,G9 = 0.989 p.u., Uref,G10 = 1.079 p.u., C8 = 0.3 p.u., nt(12-11) = 1.0394, nt(12-11) = 1.0227, nt(12-11) = 1.1101
Table 4. The applicability judgment index Ism in case B.
Table 4. The applicability judgment index Ism in case B.
nIs1 (p.u.)Is2 (p.u.)Is3 (p.u.)Control Horizon (s)Prediction Horizon (s)
1(−0.0211, −0.0301,
−0.0275, −0.0198)
(0.0013, −0.0045,
−0.0019, 0.0201)
(0.0196, 0.0152,
0.0177, 0.0510)
3060
2(0.0019, −0.0028,
−0.0012, 0.0225)
(0.0201, 0.0169,
0.0184, 0.0534)
2030
3(0.0148, 0.0144,
0.0177,0.0496)
1020
Table 5. Control performance comparison.
Table 5. Control performance comparison.
Control MethodnTime Consuming per Optimization (s)tcn (s)tpn (s)Total Time Consuming (s)ΔVoffset (p.u.)
AH-MPC11.91930603.5673.4210 × 10−2
21.0792030
30.5691020
MPC11.91430607.5253.6610 × 10−2
21.876
31.902
41.833

Share and Cite

MDPI and ACS Style

Zhang, Y.; Liu, M.; Zhang, W.; Sun, W.; Hu, X.; Kong, G. Power System Voltage Correction Scheme Based on Adaptive Horizon Model Predictive Control. Appl. Sci. 2018, 8, 641. https://0-doi-org.brum.beds.ac.uk/10.3390/app8040641

AMA Style

Zhang Y, Liu M, Zhang W, Sun W, Hu X, Kong G. Power System Voltage Correction Scheme Based on Adaptive Horizon Model Predictive Control. Applied Sciences. 2018; 8(4):641. https://0-doi-org.brum.beds.ac.uk/10.3390/app8040641

Chicago/Turabian Style

Zhang, Yan, Meng Liu, Wen Zhang, Wenchuan Sun, Xingwang Hu, and Gang Kong. 2018. "Power System Voltage Correction Scheme Based on Adaptive Horizon Model Predictive Control" Applied Sciences 8, no. 4: 641. https://0-doi-org.brum.beds.ac.uk/10.3390/app8040641

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