The example described in this section illustrates modeling of series compensation and related phenomena such as subsynchronous resonance in a transmission system.
The singleline diagram shown here represents a threephase, 60 Hz, 735 kV power system transmitting power from a power plant consisting of six 350 MVA generators to an equivalent system through a 600 km transmission line. The transmission line is split into two 300 km lines connected between buses B1, B2, and B3.
Series and Shunt Compensated Transmission System
To increase the transmission capacity, each line is series compensated by capacitors representing 40% of the line reactance. Both lines are also shunt compensated by a 330 Mvar shunt reactance. The shunt and series compensation equipment is located at the B2 substation where a 300 MVA735/230 kV transformer feeds a 230 kV250 MW load.
Each series compensation bank is protected by metaloxide varistors (MOV1 and MOV2). The two circuit breakers of line 1 are shown as CB1 and CB2.
This power system is available in the power_3phseriescomp
power_3phseriescomp
model.
Load this model and save it in your working directory as case1
to
allow further modifications to the original system.
Compare the SimPowerSystems™ circuit model (SeriesCompensated System (power_3phseriescomp)) with the schematic diagram above (Series and Shunt Compensated Transmission System). The generators are simulated with a Simplified Synchronous Machine block. A ThreePhase Transformer (Two Windings) block and a ThreePhase Transformer (Three Windings) block are used to model the two transformers. Saturation is implemented on the transformer connected at bus B2.
The B1, B2, and B3 blocks are ThreePhase VI Measurement blocks taken from the Measurements library. These blocks are reformatted and given a black background color to give them the appearance of bus bars. They output the three linetoground voltages and the three line currents. Open the dialog boxes of B1 and B2. Note how the blocks are programmed to output voltages in pu and currents in pu/100 MVA. Notice also that the voltage and current signals are sent to internal Goto blocks by specifying signal labels. The signals are picked up by the From blocks in the Data Acquisition subsystem.
The fault is applied on line 1, on the line side of the capacitor bank. Open the dialog boxes of the ThreePhase Fault block and of the ThreePhase Breaker blocks CB1 and CB2. See how the initial breaker status and switching times are specified. A linetoground fault is applied on phase A at t = 1 cycle. The two circuit breakers that are initially closed are 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 the line opening.
SeriesCompensated System (power_3phseriescomp)
Now, open the Series Compensation1 subsystem of the power_3phseriescomp
model.
The threephase module consists of three identical subsystems, one
for each phase. A note indicates how the capacitance value and the
MOV protection level are calculated. Open the Series Compensation1/Phase
A subsystem. You can see the details of the connections of the series
capacitor and the Surge Arrester block
(renamed MOV). The transmission line is 40% series compensated by
a 62.8 µF capacitor. The capacitor is protected by the MOV block.
If you open the dialog box of the MOV block, notice that it consists
of 60 columns and that its protection level (specified at a reference
current of 500 A/column or 30 kA total) is set at 298.7 kV. This voltage
corresponds to 2.5 times the nominal capacitor voltage obtained at
a nominal current of 2 kA RMS.
A gap is also connected in parallel with the MOV block. The gap is fired when the energy absorbed by the surge arrester exceeds a critical value of 30 MJ. To limit the rate of rise of capacitor current when the gap is fired, a damping RL circuit is connected in series. Open the Energy & Gap firing subsystem. It shows how you calculate the energy dissipated in the MOV 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.
Series Compensation Module
Series Compensation1/PhaseA Subsystem
Series Compensation1/PhaseA Subsystem/Energy and Gap Firing
Open the 300 MVA 735/230 kV Transformer dialog box and notice that the currentflux saturation characteristic is set at
[0 0 ; 0.0012 1.2; 1 1.45] in pu
These data are the current and flux values at points 1, 2, and 3 of the piecewise linear approximation to the flux linkage curve shown here.
Saturable Transformer Model
The fluxcurrent characteristic is approximated by the two segments shown in the graph here. The saturation knee point is 1.2 pu. The first segment corresponds to the magnetizing characteristic in the linear region (for fluxes below 1.2 pu). At 1 pu voltage, the inductive magnetizing current is 0.0010/1.0 = 0.001 pu, corresponding to 0.1% reactive power losses.
The iron core losses (active power losses) are specified by the magnetization resistance Rm = 1000 pu, corresponding to 0.1% losses at nominal voltage.
The slope of the saturation characteristic in the saturated region is 0.25 pu. Therefore, taking into account the primary leakage reactance (L1 = 0.15 pu), the air core reactance of the transformer seen from the primary winding is 0.4 pu/300 MVA.
Before performing transient tests, you must initialize your model for the desired load flow. Use the load flow utility of the Powergui to obtain an active power flow of 1500 MW out of the machine with a terminal voltage of 1 pu (13.8 kV).
Open the Powergui block and select Machine
Initialization. A new window appears. In the upper right
window you have the name of the only machine present in your system.
Its Bus type should be PV
Generator
and the desired Terminal
Voltage should already be set to the nominal voltage of 13800
V
. In the Active Power field,
enter 1500e6
as the desired output power. Click
the Compute and Apply button. Once
the load flow is solved, the phasors of AB and BC machine voltages
as well as currents flowing in phases A and B are updated in the left
window. The required mechanical power to drive the machine is displayed
in watts and in pu, and the required excitation voltage E is displayed
in pu.
Pmec 

E/Vf 

Notice that Constant blocks containing these two values are
already connected to the Pm
and E
inputs
of the machine block. If you open the Machine block dialog box, you
see that the machine initial conditions (initial speed deviation dw
= 0; internal angle theta, current magnitudes, and phase angles) are
automatically transferred in the last line.
Once the load flow is performed, you can obtain the corresponding
voltage and current measurements at the different buses. In the Powergui
block, select Steady State Voltages and Currents.
You can observe, for example, the phasors for phase A voltages at
buses B1
, B2
, and B3
and
the current entering line 1 at bus B1
.
B1/Va 

B2/Va 

B3/Va 

B1/Ia  1560 A ; 30.50 degrees 
The active power flow for phase A entering line 1 is therefore
$${P}_{a}={V}_{a}\cdot {I}_{a}\cdot \mathrm{cos}\left({\phi}_{a}\right)=\frac{608.8\text{}kV}{\sqrt{2}}\cdot \frac{1.56\text{}kA}{\sqrt{2}}\cdot \mathrm{cos}\left(\mathrm{30.50.}18.22\right)=464\text{}MW$$
corresponding to a total of 464 * 3 = 1392 MW for the three phases.
To speed up the simulation, you need to discretize the power
system. The sample time is specified in the Powergui block
as a variable Ts
. The sample time Ts=50e6
has
already been defined in the Model Initialization function in the Callbacks
of the Model Properties. The sample time Ts
is
also used in the Discrete Integrator block of the MOV energy calculator
controlling the gap.
Ensure that the simulation parameters are set as follows.
Stop time 

Solver options type 

Fixed step size 

Ensure that the fault breaker is programmed for a linetoground fault on phase A. Start the simulation and observe the waveforms on the three scopes. These waveforms are shown here.
Simulation Results for a FourCycle LinetoGround Fault at the End of Line 1
The simulation starts in steady state. At the t = 1 cycle, a linetoground fault is applied and the fault current reaches 10 kA (a: trace 3). During the fault, the MOV conducts at every half cycle (b: trace 2) and the energy dissipated in the MOV (b: trace 3) builds up to 13 MJ. At t = 5 cycles the line protection relays open breakers CB1 and CB2 (see the three line currents on trace 2) and the energy stays constant at 13 MJ. As the maximum energy does not exceed the 30 MJ threshold level, the gap is not fired. At the breaker opening, the fault current drops to a small value and the line and series capacitance starts to discharge through the fault and the shunt reactance. The fault current extinguishes at the first zero crossing after the opening order given to the fault breaker (t = 6 cycles). Then the series capacitor stops discharging and its voltage oscillates around 220 kV (b: trace 1).
Doubleclick the ThreePhase Fault block to open the Block Parameters dialog box. Select the Phase B Fault and Phase C Fault check boxes, so that you now have a threephasetoground fault.
Restart the simulation. The resulting waveforms are shown.
Simulation Results for a FourCycle ThreePhasetoGround Fault at the End of Line 1
Note that during the fault the energy dissipated in the MOV (b: trace 3) builds up faster than in the case of a linetoground fault. The energy reaches the 30 MJ threshold level after three cycles, one cycle before the opening of the line breakers. As a result, the gap is fired and the capacitor voltage (b: trace 1) quickly discharges to zero through the damping circuit.
One particular characteristic of seriescompensated 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
are in the vicinity of the zeros of the system impedance. Also, high
subsynchronous voltages due to impedance poles at subsynchronous frequencies
drive transformers into saturation. The transformer saturation due
to subsynchronous voltages is illustrated at the end of this case
study. The torque amplification on a thermal machine is illustrated
in another example (see the power_thermal
power_thermal
model).
Now measure the positivesequence impedance versus frequency seen from bus B2.
The section Analyze a Simple Circuit explains how the Impedance Measurement block allows you
to compute the impedance of a linear system from its statespace model.
However, your case1
model contains several nonlinear
blocks (machine and saturation of transformers). If you connect the
Impedance Measurement block to your system, all nonlinear blocks are
ignored. This is correct for the transformer, but you get the impedance
of the system with the machine disconnected. Before measuring the
impedance, you must therefore replace the machine block with an equivalent
linear block having the same impedance.
Delete the Simplified Synchronous Machine block from your case1
model
and replace it with the ThreePhase Source block
from the Electrical Sources library. Open the block dialog box and
set the parameters as follows to get the same impedance value (L =
0.22 pu/ (6 * 350 MVA) Quality factor = 15).
Phasetophase rms voltage 

Phase angle of phase A 

Frequency (Hz) 

Internal connection Yg  Specify impedance using shortcircuit level 
3phase shortcircuit level 

Base voltage 

X/R ratio 

Save your modified model as case1Zf
.
Open the Fundamental Blocks/Measurements library and copy the Impedance Measurement block into your model. This block is used to perform the impedance measurement. Connect the two inputs of this block between phase A and phase B of the B2 bus. Measuring the impedance between two phases gives two times the positivesequence impedance. Therefore you must apply a factor of 1/2 to the impedance to obtain the correct impedance value. Open the dialog box and set the multiplication factor to 0.5.
In the Powergui block, select Impedance
Measurement. A new window opens, showing your Impedance
Measurement block name. Fill in the frequency range by entering 0:500
.
Select the linear scales to display Z magnitude vs. frequency plot.
Click the Save data to workspace button
and enter Zcase1
as the variable name to contain
the impedance vs. frequency. Click the Display button.
When the calculation is finished, the magnitude and phase as
a function of frequency are displayed in the two graphs on the window.
If you look in your workspace, you should have a variable named Zcase1
.
It is a twocolumn matrix containing frequency in column 1 and complex
impedance in column 2.
The impedance as a function of frequency (magnitude and phase) is shown here.
Impedance vs. Frequency Seen from Bus B2
You can observe three main modes: 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.
If you zoom in on the impedance in the 60 Hz region, you can find the system's shortcircuit level at bus B2. You should find a value of 58 Ω at 60 Hz, corresponding to a threephase shortcircuit power of (735 kV)^{2} / 58 = 9314 MVA.
The configuration of the substation circuit breakers normally
allows clearing a fault at the bus without losing the lines or the
transformers. You now modify your case1
model to
perform a threecycle, threephasetoground fault at bus B2:
Disconnect the ThreePhase Fault block and reconnect it so that the fault is now applied on bus B2.
Open the ThreePhase Fault block and make the following modifications in its dialog box:
Phase A, Phase B, Phase C, Ground Faults  All selected 
Transition times  [2/60 5/60] 
Transition status [1, 0, 1...]  (0/1) 
You have now programmed a threephasetoground fault applied at the t = 2 cycles.
Open the dialog boxes of circuit breakers CB1 and CB2 and make the following modifications:
Switching of Phase A  Not selected 
Switching of Phase B  Not selected 
Switching of Phase C  Not selected 
The circuit breakers are not switched anymore. They stay at their initial state (closed).
In the Data Acquisition subsystem, insert
a Selector block (from the Simulink^{®} Signals & Systems library)
in the Vabc_B2 output of bus B2 connected to the scope. Set the Elements parameter to 1
.
This allows you to see the phase A voltage clearly on the scope.
You now add blocks to read the flux and the magnetization current of the saturable transformer connected at bus B2.
Copy the Multimeter block
from the Measurements library into your case1
model.
Open the Transformer dialog box.
In the Measurements list, select Flux
and magnetization current
. Open the Multimeter block. Verify
that you have six signals available. Select flux and magnetization
current on phase A, and click OK.
You now have two signals available at the output of the Multimeter block. Use a Demux block to send these two signals on a twotrace scope.
In the Simulation
> Model Configuration Parameters dialog
box, change the stop time to 0.5
. This longer simulation
time allows you to observe the expected lowfrequency modes (9 Hz).
Start the simulation.
The resulting waveforms are plotted here.
Simulation Results for a ThreeCycle ThreePhasetoGround Fault at BusB2
The 9 Hz subsynchronous mode excited at fault clearing is clearly seen on the phase A voltage at bus B2 (trace 1) 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 5). The flux in phase A of the transformer is plotted on trace 4. 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 9 Hz voltage components. The pulses of the transformer magnetizing current appear when the flux exceeds its saturation level. This current contains a 60 Hz reactive component modulated at 9 Hz.