Next Article in Journal
Suborganizations of Institutions in Library and Information Science Journals
Previous Article in Journal
The SP Theory of Intelligence: An Overview
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Passivity-Based Nonlinear Excitation Control of Power Systems with Structure Matrix Reassignment

1
School of Computer Science, Zhengzhou University of Light Industry, Zhengzhou 450002, China
2
School of Electrical Engineering, Zhengzhou University, Zhengzhou 450001, China
3
Department of Automation, Tsinghua University, Beijing 100084, China
4
Department of Electronics & Computer Science, University of Southampton, Southampton SO17 1BJ, UK
*
Author to whom correspondence should be addressed.
Information 2013, 4(3), 342-350; https://doi.org/10.3390/info4030342
Submission received: 10 April 2013 / Revised: 12 August 2013 / Accepted: 12 August 2013 / Published: 20 August 2013

Abstract

:
Passivity-based control is widely used in electronic circuit systems because it can utilize their internal structures to facilitate the controller design. In this paper, we first propose a dissipative Hamiltonian realization of power systems and discuss the disadvantages of the traditional passivity-based excitation controller. Then, a novel excitation controller is put forward to reassign the interconnection and dissipative matrix, and the corresponding Hamiltonian function. Simulation results verify that the proposed controller can effectively improve the transient stability of the power system.

1. Introduction

The port-controlled Hamiltonian (PCH) system has gained increasing interest in the control and energy community [1]. The main characteristic of the PCH system is that the structural information is available in the model description: the total energy of the system is represented by the Hamiltonian function, which can be usually used as a Lyapunov function in stability analysis, while the interconnection property is explicitly described by the interconnection and damping matrix.
By utilizing the energy concept and the internal structural property, several methods have been developed for the controller design of PCH systems [2,3,4,5]. Interconnection and damping assignment passivity-based control (IDA-PBC) is a widely used state-feedback technique which can stabilize the considered nonlinear systems by shaping the energy and re-constructing the structure matrices. So far, IDA-PBC has been successfully applied to electromechanical systems [6,7], power electronics [8], and power systems [9,10,11,12,13,14,15].
For power systems, the energy shaping is generally achieved by modifying the energy transfer pattern between the mechanical and the electrical components, which is obtained by injecting damping into the electrical dynamic. The conventional IDA-PBC excitation controller cannot directly reassign the mechanical damping in the swing equation, though the damping in the swing equation is very important for transient stability improvement. In this paper, we propose a novel IDA-PBC excitation controller to enhance the transient stability of power systems by choosing a convenient interconnection matrix and find a solution for the matched partial derivative equation. The proposed excitation controller not only provides a compensation damping into electrical dynamic equation, but also into the mechanical swing equation. Simulation results verify that, compared to the conventional energy-shaping method [14] and PSS+AVR controller [16], the proposed control scheme can effectively improve the transient stability of the power system.
The rest of the paper is organized as follows. In section 2 we provide some preliminary results on the IDA-PBC method. In Section 3 we propose the dynamic model of the single machine infinite bus power systems and put forward a novel IDA-PBC excitation controller. In Section 4, we simulate the power system, to illustrate the effectiveness of the proposed control scheme. Section 5 summarizes the results and concludes.

2. IDA-PBC Methodology of Nonlinear Systems

Consider the following nonlinear system:
x ˙ = f ( x ) + g ( x ) u
where x R n is the state vector and u R m is the control input. f ( x ) and g ( x ) are smooth vector functions in the appropriate dimensions.
Suppose the system can be represented as the following PCH formulation,
x ˙ = [ J ( x ) R ( x ) ] H ( x ) + g ( x ) u
where J ( x ) = J T ( x ) and R ( x ) = R T ( x ) 0 are the interconnection and damping matrices respectively. H ( x ) is the Hamiltonian function, which represents the total energy stored in the system. H ( x ) is the gradient vector function of H ( x ) with respect to x.
If the Hamiltonian function H ( x ) has a strict minimum at the equilibrium point, it can be chosen as the Lyapunov candidate for the stability analysis and passivity-based controller design. Unfortunately, for a physical system the Hamiltonian function does not always have a minimum at the desired equilibrium or the dynamic performance of the system is not satisfactory; so the IDA-PBC is always needed to shape the energy function, the interconnection, and the damping matrix to get a desired Hamiltonian function or to improve the dynamic performance of the system.
Generally, for a dissipative Hamiltonian realization of the system in Equation (1), we have following result [4]:
Lemma 1
Consider the system in Equation (1). Assume there exist matrices g ( x ) , J d ( x ) = J d T ( x ) , R d ( x ) = R d T ( x ) 0 and differentiable scalar function H d ( x ) such that the following partial differential equation (PDE) holds
g ( x ) f ( x ) = g ( x ) [ J d ( x ) R d ( x ) ] H d ( x )
where g is a left annihilator of g ( x ) , i.e., g ( x ) g ( x ) = 0 and H d ( x ) is such that
x * = arg min H d ( x )
with x * the desired equilibrium point to be stabilized. Then, under the following feedback controller
u ( x ) = g ( x ) g ( x ) 1 g J d ( x ) R d ( x ) H d f ( x )
the closed loop system can be represented as
x ˙ = [ J d R d ( x ) ] H d ( x )
with x * being a local stable equilibrium. Furthermore, if the largest invariant set
x R n | | H d ( x ) | T R d ( x ) H d = 0 = { x * }
the closed loop system is asymptotically stable.
According to Lemma 1, the main procedure of the IDA-PBC methodology can be briefly described as follows:
(i)
Fix the desired structure of the interconnection and damping matrix.
(ii)
Derive a partial differential equation (PDE) parameterized by the chosen matrices whose solutions characterize all the energy functions that can be assigned.
(iii)
Choose one solution that satisfies the minimum requirement and computer the control.
Thus, in order to formulate the considered nonlinear system within a PCH system, the key step is to find a solution to the PDE in Equation (3). In the next section, we construct the pon-parameterized IDA and parameterized IDA methods proposed in [4], in order to reassign the interconnection and dissipative matrices of the single machine infinite bus power system and formulate a feedback controller to improve its transient stability.

3. IDA-PBC Excitation Control of Power Systems

Consider a single machine infinite bus power system as shown in Figure 1.
Figure 1. the configuration of a single machine infinite bus power system.
Figure 1. the configuration of a single machine infinite bus power system.
Information 04 00342 g001
The dynamic of the system can be presented as
δ ˙ = ω ω 0 ω ˙ = ω 0 M P m D M ( ω ω 0 ) ω 0 E q V s sin δ M x d Σ E ˙ q = E q T d + x d x d T d 0 x d Σ V s cos δ + 1 T d 0 V f
where δ is the power angle of the generator, in radians (rad); ω is the rotor speed of the generator, ω 0 = 2 π f 0 , in rad/s; E q is the q-axis internal transient voltage of the generator, in per unit; V f is the voltage of the field circuit of the generator—the control input, in per unit; P m is the mechanical power, assumed to be constant, in per unit; V s is the voltage of the infinite bus, in per unit; x d is the d-axis transient reactance of the generator, in per unit; x d is the d-axis reactance, in per unit; M is the inertia coefficient of the generator, in s; D is the damping constant, in per unit; T d 0 is the d-axis transient open-circuit time constant, in s; T d is the d-axis transient closed-circuit time constant, in s; x d Σ = x d + 1 2 x L + x T , where x L is the reactance of the transmission line, in per unit; x T is the reactance of the transformers, in per unit.
Let x 1 = δ , x 2 = ω ω 0 , x 3 = E q and u = V f . Denote a = ω 0 P m M , b = D M , c = ω 0 V s M x d Σ , d = 1 T d , and e = x d x d T d 0 x d Σ V s . Direct calculation shows that the system can be Hamiltonian realized as follows:
x ˙ = [ J ( x ) R ( x ) ] H + g ( x ) u
where g ( x ) = 0 , 0 , 1 T d 0 T . The interconnection and dissipative matrices are
J ( x ) = 0 1 0 1 0 0 0 0 0 , R ( x ) = 0 0 0 0 b 0 0 0 e c
respectively. The corresponding Hamiltonian function is
H ( x ) = 1 2 x 2 2 + c d 2 e x 3 2 c x 3 cos x 1 a x 1
It can be verified that H ( x ) possesses a strict local minimum at the equilibrium point x 1 * = 1 2 arcsin 2 a d c e , x 2 * = 0 , and x 3 * = e d cos x 1 * .
Noticing that the largest invariant set { H ˙ ( x ) = 0 } contains only the equilibrium point, the power system can be stabilized by a damping injection controller as follows:
u = k L g H ( x )
where k > 0 is the feedback gain.
The closed loop system can be written as
x ˙ = [ J ( x ) R c ( x ) ] H
where the dissipative matrix is
R c ( x ) = 0 0 0 0 b 0 0 0 e c + k
It can be seen from the matrix R c ( x ) that the damping injection excitation controller cannot change the mechanic-electric dynamic of the system and thus cannot effectively attenuate the oscillation.
In the following, we use the IDA-PBC method to design an effective excitation controller to improve the performance of the system. First, assume the desired dissipative Hamiltonian realization of the power system can be formulated as
x ˙ = ( J d ( x ) R d ( x ) ) H d
where the dissipative matrix R d = R c . Choose the interconnection matrix as
J d ( x ) = 0 1 α 1 0 0 α 0 0
Suppose the desired Hamiltonian function can be chosen as
H d ( x ) = H ( x ) + H i n j ( x )
Let the feedback stabilization controller u = ϕ ( x ) . Then the closed loop system satisfies
x ˙ = 0 1 α 1 b 0 α 0 ( e c + k ) ( H ( x ) + H i n j ( x ) ) = 0 1 0 1 b 0 0 0 e c H ( x ) + 0 0 1 T d 0 ϕ ( x )
Thus we have
H i n j x 3 = 0 H x 1 b H i n j x 2 = 0 u = α H x 1 + H i n j x 1 k H x 3
From above equations we can see that H i n j does not depend on x 3 . Taking into consideration that
H x 1 = c x 3 sin x 1 a H x 2 = x 2 H x 3 = c d e x 3 c cos x 1
we have
H i n j = ψ ( x 1 ) + 1 b ( a c x 3 sin x 1 ) x 2 u = α c x 3 sin x 1 a + H i n j x 1 k ( c d d x 3 c cos x 1 )
where ψ ( x 1 ) is chosen such that
H d x ( x 1 * , x 2 * , x 3 * ) = 0 , H e s s ( H d ) | ( x 1 * , x 2 * , x 3 * ) > 0
hold. Noticing that
H d = H + H i n j = 1 2 x 2 2 + c d 2 e x 3 2 c x 3 cos x 1 a x 1 + ψ ( x 1 ) + 1 b ( a c x 3 sin x 1 ) x 2
we have
H d x 1 = c x 3 sin x 1 a + ψ x 1 c b x 2 x 3 cos x 1 H d x 2 = x 2 + 1 b ( a c x 3 sin x 1 ) H d x 3 = c d e x 3 c cos x 1 c b x 2 sin x 1
So, in order to guarantee that H d has a strict local minimum at the equilibrium point, we can choose
ψ ( x 1 ) x 1 = c o n s t = γ
Let ψ ( x 1 ) = γ x 1 , then we have
H e s s ( H d ) = c x 3 cos x 1 + γ c b x 2 x 3 sin x 1 c b x 3 cos x 1 c sin x 1 + c b cos x 1 c b x 3 cos x 1 1 c b sin x 1 c x 1 c b x 2 cos x 1 c b sin x 1 c d e
and
H e s s ( H d ( x * ) ) = c e d cos 2 x 1 * + γ c e b d cos 2 x 1 * c sin x 1 * + c b cos x 1 * c e b d cos 2 x 1 * 1 1 c sin x 1 * c b sin x 1 * 1
In order to guarantee a positive definite H d ( x * ) , we must choose γ such that
γ > c e d cos 2 x 1 * and γ > η c d e c 2 b 2 sin 2 x 1 *
hold simultaneously, where
η = c 3 e b 3 d cos 3 x 1 * sin x 1 * c sin 2 x 1 * c b cos x 1 * sin x 1 * c 3 e 2 b 3 d 2 cos 4 x 1 * sin x 1 *
According to Lemma 1, under the desired Hamiltonian realization the stabilization control scheme of the system can be chosen as
u = α c x 3 sin x 1 a + γ c b x 2 x 3 cos x 1 k c d e x 3 c sin x 1
where k > 0 is the feedback gain. The corresponding Hamiltonian function of the closed loop system is
H d = 1 2 x 2 2 + c d 2 e x 3 2 c x 3 cos x 1 + ( γ a ) x 1 + 1 b ( a c x 3 sin x 1 ) x 2
Remark 1
The internal transient voltage of the generator can be presented by measurable signals such as P e , Q e and V s . Noting that the active and reactive power of the generator can be written as P e = E q V s x d Σ sin δ and Q e = E q V s x d Σ cos δ V s 2 x d Σ respectively, we have E q = x d Σ V s 2 ( Q e V s 2 x d Σ ) 2 + P e 2 .

4. Simulation

To show the effectiveness of the proposed controller, we compare it with the conventional IDA-PBC controller [14] and PSS+AVR controller [16]. During the simulation, the effects of generator fields and rotor circuit dynamics are included in the model of power system [16]. The simulation is completed on a professional power system simulation platform developed by the China Electric Power Research Institute.
Figure 2. Response of the rotor angle under a circuit fault.
Figure 2. Response of the rotor angle under a circuit fault.
Information 04 00342 g002
During the simulation, a three-phase temporary short-circuit fault is assumed to occur at the transmission line, starting at time t = 1 s and cleared at t = 1.1 s . The parameters in the proposed feedback excitation controller are chosen as k = 2 and γ = 5 . The simulation results are shown in Figure 2 and Figure 3, where the solid black lines indicate the response under the proposed stabilization controller, the blue dashed lines corresponding to the response under conventional IDA-PBC controller and the red dotted lines correspond to PSS+AVR respectively.
Figure 3. Response of the terminal voltage of generator under a circuit fault.
Figure 3. Response of the terminal voltage of generator under a circuit fault.
Information 04 00342 g003
From the simulation results, it can be seen that, when the fault occurs, the proposed nonlinear stabilization control scheme makes the system responds much faster than the conventional passivity-based controller and PSS+AVR controller . It follows that the proposed controller outperforms the traditional controller in transient stability enhancement of the power system.

5. Conclusions

The excitation control of a single machine infinite bus power system based on the IDA-PBC method is investigated . We propose a novel excitation controller by choosing the desired interconnection and damping matrix, and solving the PDE to get the corresponding Hamiltonian function. An excitation controller is put forward based on the new Hamiltonian realization system, which not only injects damping into the electronic equation, but also reshapes the mechanical swing equation. Simulation results show that the proposed control scheme can effectively improve the transient stability of the power system.

Acknowledgements

This work was supported by the National Natural Science Foundation of China (No. 60974005 and 61104022), the Specialized Research Fund for the Doctoral Program of Higher Education (No. 20094101120008) and the Science and Technique Research Program of Henan Educational Committee (No. 13A520379).

Conflict of Interest

The authors declare no conflict of interest.

References

  1. Van der Schaft, A.J. Port-Hamiltonian Systems: Network Modeling and Control of Nonlinear Physical Systems. In Advanced Dynamics and Control of Structures and Machines; CISM Courses and Lectures, Volume 444; Springer: New York, NY, USA, 2004. [Google Scholar]
  2. Ortega, R.; van der Schaft, A.J.; Mareels, I.; Maschke, B. Putting energy back in control. IEEE Control Syst. Mag. 2001, 21, 18–33. [Google Scholar] [CrossRef]
  3. Ortega, R.; van der Schaft, A.J.; Maschke, B.; Escobar, G. Interconnection and damping assignment passivity-based control of port-controlled Hamiltonian systems. Automatica 2002, 38, 585–596. [Google Scholar] [CrossRef]
  4. Ortega, R.; García-Causeco, E. Interconnection and damping assignment passivity-based control: A survey. Eur. J. Control 2004, 10, 432–450. [Google Scholar] [CrossRef]
  5. Ortega, R.; Castaños, F.; Astolfi, A. Control by interconnection and standard passivity-based control of Port-Hamiltonian systems. IEEE Trans. Autom. Control 2008, 53, 2527–2542. [Google Scholar] [CrossRef]
  6. Fujimoto, K.; Sakurama, K.; Surgie, T. Trajectory Tracking Control of Portcontrolled Hamiltonian Systems and its Application to a Magnetic Levitation System. In Proceedings of the 40th IEEE Conference on Decision and Control, Orlando, FL, USA, 4–7 December 2001; pp. 3388–3393.
  7. Rodríguez, H.; Ortega, R. Stabilization of electromechanical systems via interconnection and damping assignment. Int. J. Robust Nonlinear Control 2003, 13, 1095–1111. [Google Scholar] [CrossRef]
  8. Batlle, C.; Dòria-Cerezo, A.; Fossas, E. Bidirectional power flow control of a power converter using passive Hamiltonian techniques. Int. J. Circuit Theory Appl., 2008, 36, 769–788. [Google Scholar] [CrossRef]
  9. Cheng, D.; Xi, Z.; Hong, Y.; Qin, H. Energy-based Stabilization in Power Systems. In Proceedings of the 14th IFAC World Congress, Beijing, China, 5–9 July 1999; pp. 297–303.
  10. Xi, Z.; Cheng, D.; Lu, Q.; Mei, S. Nonlinear decentralized controller design for multimachine power systems using Hamiltonian function method. Automatica 2002, 38, 527–534. [Google Scholar] [CrossRef]
  11. Wang, Y.; Cheng, D.; Hong, Y. Stabilization of synchronous generators with Hamiltonian function approach. Int. J. Syst. Sci. 2001, 32, 971–978. [Google Scholar] [CrossRef]
  12. Sun, Y.; Shen, T.; Ortega, R.; Lu, Q. Decentralized Controller Design for Multi-machine Power Systems on Hamiltonian Structure. In Proceedings of the 40th IEEE Conference on Decision and Control, Orlando, USA, 4–7 December 2001; pp. 3045–3050.
  13. Xi, Z.; Cheng, D.; Lu, Q.; Mei, S. Nonlinear decentralized controller design for multimachine power systems using Hamiltonian function method. Automatica 2002, 38, 527–534. [Google Scholar] [CrossRef]
  14. Galaz, M.; Ortega, R.; Bazanella, A.S.; Stankovic, A.M. An energy-shapping approach to the design of excitation control of synchronous generators. Automatica 2003, 39, 111–119. [Google Scholar] [CrossRef]
  15. Li, J.; Liu, Y.; Li, C.; Zhang, L. Interconnection and Damping Assignment Passivity-based Excitation Control of Power Systems. In Proceedings of IEEE Control and Decision Conference, Xuzhou, China, 26–28 May 2010; pp. 954–958.
  16. Kundur, P. Power System Stability and Control; McGraw-Hill: New York, NY, USA, 1994. [Google Scholar]

Share and Cite

MDPI and ACS Style

Li, J.; Liu, Y.; Li, C.; Chu, B. Passivity-Based Nonlinear Excitation Control of Power Systems with Structure Matrix Reassignment. Information 2013, 4, 342-350. https://doi.org/10.3390/info4030342

AMA Style

Li J, Liu Y, Li C, Chu B. Passivity-Based Nonlinear Excitation Control of Power Systems with Structure Matrix Reassignment. Information. 2013; 4(3):342-350. https://doi.org/10.3390/info4030342

Chicago/Turabian Style

Li, Jianyong, Yanhong Liu, Chunwen Li, and Bing Chu. 2013. "Passivity-Based Nonlinear Excitation Control of Power Systems with Structure Matrix Reassignment" Information 4, no. 3: 342-350. https://doi.org/10.3390/info4030342

Article Metrics

Back to TopTop