Digital Simulation of Power Systems and Power Electronics Using The MATLAB-Simulink Power System Blockset
Digital Simulation of Power Systems and Power Electronics Using The MATLAB-Simulink Power System Blockset
Digital Simulation of Power Systems and Power Electronics using the NIATLAB/Simulink Power System Blockset
by Gilbert Sybille Power System Simulation Laboratory IREQ - Hydro-Quebec Varennes, Quebec, CANADA and Hoang Le-Huy Dept. of Electrical and Computer Engineering Universite Lava1 Ste-Foy, Quebec, CANADA
I. DESCRIPTION
Digital simulation of power systems has been intensively studied during the last three decades because of increasing needs in the domain. Various approaches have been considered for the formulation and solution of the network equations. The main simulation methods that have been developed include state-variable analysis, nodal analysis, and modified nodal analysis. These methods have been applied to the development of different simulation packages. The most popular ones are EMTP for power networks (using nodal analysis with fixedstep integration) and SPICE for electronic circuits and power electronics (using modified nodal analysis with variable-step integration). The Power System Blockset (PSB) is a graphic tool that allows building schematics and simulation of power systems in the Simulink environment. The blockset uses the MATLAB/Simulink environment to represent common components and devices found in electrical power networks. Simulinks variable-step event-sensitive integration algorithms allow increased accuracy in zero-crossing detection of currents as compared with fixedstep algorithms. One important feature of the Power System Blockset which has been introduced with version 2.0 is its ability to simulate either with continuous variable time-step integration algorithms or with a discretized system. For small and medium size systems, variable time steps algorithms are usually faster because the number of steps will be less than with a fixed-time step method. However, for large systems which contain many states or many non-linear blocks such as power electronic switches, it is advantageous to discretize the electrical system. Both the continuous integration and fixed-time step methods will be illustrated. The graphic interface of Simulink provides a user-friendly environment where the power circuit and control system are represented in the same diagram. The results is displayed while the simulation is running. The processing power of MATLAB allows the designer to perform complex post-processing on simulation results. This special technical session presents the features of the Power System Blockset and its applications in the simulation of power systems and power electronics systems.
2973
[-qpqpq-jjjl~pq-Gq~~
Electrical Sources Elements Power Electronics Machines Connectors Measurements
Power System BIockset Copyright [c) 1997-1 998 by TEQSIM International, Inc., under sublicense from Hydro-Quebec, and The Mathworks, Inc. All Rights Resewed.
a
powergui
Block library
Elements in the PSB block library are classified in various groups according to their nature: Electrical Sources, Elements, Power Electronics, Machines, Connectors, and Measurements. The Elements group contains single-phase models RLC branches and loads, linear and saturable transformers, mutual inductances, E-section lines, MOV type surge arrester, circuit breaker and n-phase distributed-parameter line model. Using the masking facility of Simulink, the user can easily add more complex elements built from the basic PSB building blocks and associate a dialog box. This technique has been used to develop a three-phase library which is also provided. The Power Electronics group contains common semiconductor devices. Each element in this group (except Diode) has a Simulink gating control input and a Simulink output returning switch current and voltage. The Machines group contains simplified and detailed models of synchronous machine, an asynchronous machine, a permanent magnet synchronous machine, a model of hydraulic turbine governor, and an excitation system. Each machine block has a Simulink output returning measurements of internal variables.
This example illustrates phenomena related to subsynchronous resonance in a series-compensated AC transmission network. "Descriptionof the Transmission Network
6'350 MVA 13.8 kVl735 kV DNS MOVl
line 1; 300 km
MOV2
loo MW
- ... .. .
2975
I '
I I C
Ll
C7l
E-MOV
Energy WwMd
MOV
Latch
A gap is also connected in parallel with the MOV. The gap is fired when the energy absorbed by the surge arrester exceeds a critical value of 30 MJ. In order to limit the rate of rise of capacitor current when the gap is fired, a damping RL circuit is connected in series. The energy dissipated in the MOV is calculated by Energy & Gap firing subsystem. Energy is obtained by integrating the power (product of the MOV voltage and current). When the energy exceeds the 30 MJ threshold a closing order is sent to the breaker block simulating the gap.
Three-phase Fault on Line I Simulation results obtained for a three-phase-to-ground fault applied on line 1 are shown on Figure 5. Fault is applied at t=lcycle in substation B2 at the line side of the series compensation bank. The two line circuit breakers are initially closed and then open at t=5 cycles, simulating a fault detection and opening time of 4 cycles. The fault is eliminated at t=6 cycles, one cycle after line opening. The energy dissipated in the MOV (trace 6) builds up and reaches the 30 MJ threshold level after 3 cycles, one cycle before opening of the line breakers. As a result, the gap is fired and the capacitor voltage (trace 4)quickly discharges to zero through the damping circuit. Frequency Analysis One particular characteristic of series compensated systems is the existence of subsynchronous modes (poles and zeros of the system impedance below the fundamental frequency). Dangerous resonances can occur if the mechanical torsion modes of turbine/generator shafts is in the vicinity of the zeros of the system impedance. Also, high subsynchronous voltages due to impedances poles at subsynchronous frequencies will drive transformers into saturation. The transformer saturation due to subsynchronous voltages will be illustrated below. The PSB has the ability to compute the impedance vs. frequency between any pair of nodes. The positive-sequence impedance as function of frequency (magnitude and phase) seen from bus B2 is shown in Figure 6. Three main modes can be observed: 9 Hz, 175 Hz and 370 Hz.The 9 Hz mode is mainly due to a parallel resonance of the series capacitor with the shunt inductors. The 175 Hz and 370 Hz modes are due to the 600 km distributed parameter line. These three modes are likely to be excited at fault clearing. The system short circuit level at bus B2 is obtained from the impedance (58 Q) at 60 Hz. It corresponds to a three-phase short circuit power of (735kV) 2/58 = 9314 MVA
I
0.02
I
0.04
I
0.06 0.06
I
0.34
01 .
0.12
o.ia
o.le
1 0.2
Gap firing
2977
2000.
50
100
150
200
250
300
350
400
450
500
100
7 f ;
~
I
I I
^ _ _ _ _ _
I
I
- _ - I
I
I I I
~
50
0
I I
I I
I I
I
I
I I I
. _ - I
_ _ - _ L
_-__
_ - _ _
, - - _- _- _
I I I I
7 I I
_ _ _ _
I I
I
I
-100
50
100
150
200
250
Frequency (Hz)
300
350
400
450
500
Transientperformancefor a fault at bus B2 The simulation results obtained for a three-cycle, three-phase-to-ground fault at bus B2 as shown on Figure 7.
z
3
I I
I
I
I
I
I I
I
I
I
I
I I
I I
I
I I I
I
I
I
I
I
I
I
I
I I
I
I
I
I
-5000' 0
0.05
0.1
0.15
0.2
0.25
Time (s)
0.3
0.35
0.4
0.45
0.5
and capacitor voltage (trace 3). The 9 Hz voltage component appearing at bus B2 drives the transformer into saturation as shown on the transformer magnetizing current (trace 4). The flux in phase A of the transformer is plotted on trace 5. At fault application the voltage at transformer terminals drops to zero and the flux stays constant during the fault. At fault clearing, when the voltage recovers, the transformer is driven into saturation as a result of the flux offset created by the 60 Hz and 9Hz voltage components. The pulses of the transformer magnetizing current (trace 4) appear when the flux exceeds its saturation level of 1.2 p.u. The transformer magnetizing current contains a 60 Hz reactive component modulated at 9 Hz.
e1
Vab-inv
a+A
"-'
- -+r
A
-L
c p
d
c+c
cFihers
+c
-* pubes
IGBT INVERTER
25 k V / 2 kV
Vdc
kref
phi(deg)
ldq @U)
I
MVA,
Phase
Specify Ts-2c-6 s
lerel
Rate Limiter
CONTROL
Pulse Generator
3 kHZ
PWM
2979
Modulation index
VDC Regulator
Iq Regulator
.J3 x -Vdc 2 - h
= m x 0,612 x Vdc
Therefore, the minimum DC voltage required to generate 2000 V rms (rated secondary voltage) is 2000/0.612 = 3300 V rms. A DC bus voltage of 4000 V rms has been specified.
DSTATCOM start-up and step on reference current This system is discretized with a time step Ts=2e-6. In this test, a step is applied to the current reference input in order to observe the dynamic response when the DSTATCOM changes from full inductive to capacitive operation. When the simulation starts, the DC capacitor starts charging. This requires a Id component corresponding to the active power absorbed by the capacitor. When the DC voltage reaches its reference value (4000 V), the Id component drops to a value very close to zero and the Iq component stays at the 1 pu reference value (3 Mvar inductive). At t=O. 1 s the Iq reference current changes to + 1 pu. When the DSTATCOM changes from inductive to capacitive, a 180 degrees phase shift of the current Ia with respect to voltage Va is observed. The transition from inductive to capacitive is shown on Figure 10. It is noted that the modulation index increases (inverter
0-7803-5935-6/00/$10.00(c) 2000 IEEE 2980
'=a
I I
I I
I I
I
I
s-2'
0.05
0.06
0.07
0.08
0.09
0.1
0.11
0.12
0.13
0.14
0.15
5
v
I I
I
I I I
I
I I
I
I
I I I I
I
I
I
0I
I
1
I I
1
I
I I
g2 0.05
U 0 t
0.06
I
0.07
I
0.08
I
0.09
I
0.1
I
0.11
I
0.12
I
0.13
1
0.14
I
0.15
I
.c 0 .Q
U
0
= I
0.5
0.05
0.06
0.07
0.08
0.09
0.1
Time (s)
0.11
0.12
0.13
0.14
0.15
Figure 10 Step on operating point from inductive (Id=-1 pu) to capacitive (Id=+l pu)
Harmonic analysis In order to illustrate the Matlab post-processing capabilities, the frequency spectrum of the voltages at the inverter terminals as well as at the inverter output as also been computed. This is obtained from the standard Matlab fft function. Results are shown on Figure 11 for the STATCOM absorbing 3 Mvar (Iq =-1 pu). As expected, harmonics are clustered around multiples of the 3 kHz chopping frequency. All harmonics have been reduced to below 0.5% on the 25 kV bus by the filters. Also, the fundamental component of Va (indicated on Va spectrum bottom graph) drops to 0.88 pu, when the STATCOM absorbs its rated reactive current. This is the expected voltage drop for an inductive current of 1 pu/3 MVA on a 25 MVA network.
2981
5000
-5000
0 . 3 8 2
'
0.i84
0.186
0.188
0.1s
I
0.192
time (5) 0.194
0 1 6 . 9
0 1 8 . 9
0.2
0.202
60
Fundamental
2 1 i . 3 6 8 3 8
0 5.
0 -
-0.5
-i
0.5 0.4
.. -.
Fundam#%ntal0.88395 =
2 = & .
0 3-. .
I. .
Figure 11 Harmonic analysis of voltages top: Vab inverter voltage in volts, bottom: line-to-groundvoltage at 25 kV bus in pu
VI. CONCLUSION
This special technical session has presented the Power System Blockset as a unified framework for the simulation of power systems including power networks, power electronics and control systems. The simulation tool has been developed using state-variable approach and runs in MATLABNmulink environment. Application examples for power systems including power networks, power electronics and electric machines have been presented to illustrate the capability of the PSB to deal with various systems in power engineering. The simulation results obtained with PSB are comparable to those obtained by popular simulation packages (EMTP and PSpice). The user friendly graphical interfaces of Simulink and PSB allow faster development and closed-loop testing of control systems used in power systems and drives.
2982