Data Stream
Data Stream
DOCTOR OF PHILOSOPHY
2013
JINSHENG PENG
Contents
Contents ........................................................................................................................... 3
List of Figures .................................................................................................................. 7
List of Tables ................................................................................................................. 15
List of Abbreviations .................................................................................................... 17
Abstract .......................................................................................................................... 19
Declaration ..................................................................................................................... 21
Copyright Statement ..................................................................................................... 23
Acknowledgement ......................................................................................................... 25
Chapter 1 Introduction ............................................................................................. 27
1.1 Background ...................................................................................................... 27
1.1.1 Transformer energisation inrush phenomena ............................................... 27
1.1.2 Adverse effects of transformer energisation transients................................. 30
1.2 Objectives of research ...................................................................................... 33
1.3 Outline of the thesis .......................................................................................... 34
Chapter 2 Literature Review on Transformer Energisation Transients ............. 37
2.1 Approaches for calculating transformer inrush current .................................... 37
2.1.1 Simple analytical approaches for calculating inrush current ........................ 37
2.1.2 Numerical approaches for calculating inrush current ................................... 40
2.2 Modelling system components in EMTP for studying transformer energisation
transients ...................................................................................................................... 41
2.2.1 Transformer modelling ................................................................................. 41
2.2.2 Overhead line and cable modelling .............................................................. 47
2.2.3 Circuit breaker modelling ............................................................................. 53
2.2.4 Source and network equivalent modelling.................................................... 55
2.2.5 System load modelling.................................................................................. 56
2.3 Investigation case studies on transformer energisation transients .................... 56
2.3.1 Sympathetic interaction between transformers ............................................. 56
2.3.2 Mechanical forces induced by transformer inrush current ........................... 65
2.3.3 Energising transformers from a limited capacity generator.......................... 65
2.3.4 Harmonic incursion due to transformer energisation ................................... 67
2.3.5 Voltage dips caused by transformer energisation ......................................... 72
2.3.6 Statistical assessment of transformer energisation transients ....................... 82
2.4 Possible approaches for mitigating transformer inrush .................................... 85
3
Contents
4
Contents
5
Contents
6.5.4 Effect of transformer winding connections on voltage dips propagation ... 189
6.5.5 Voltage dips caused by consecutive energisation of wind turbine
transformers ........................................................................................................... 190
6.6 Sympathetic inrush between wind turbine transformers ................................ 192
6.6.1 Sympathetic inrush caused by energisation of multiple transformers ........ 192
6.6.2 Sympathetic inrush caused by independent energisation ........................... 195
6.7 Identification of energisation sequence resulting in less sympathetic inrush
between wind turbine transformers ........................................................................... 198
6.7.1 Sympathetic inrush level............................................................................. 198
6.7.2 Energisation sequence ................................................................................ 199
6.7.3 Study procedure correlating sympathetic inrush level and energisation
sequence ................................................................................................................. 199
6.7.4 Assessment of sympathetic inrush level under different energisation
sequences using deterministic approach ................................................................ 200
6.7.5 Assessment of sympathetic inrush level under different energisation
sequences using stochastic approach ..................................................................... 203
6.8 Summary ........................................................................................................ 206
Chapter 7 Conclusion and Future Work .............................................................. 209
7.1 Concluding remarks ....................................................................................... 209
7.2 Future work .................................................................................................... 213
References .................................................................................................................... 215
Appendix: List of Publications................................................................................... 225
6
List of Figures
List of Figures
Figure 1-1 Qualitative illustration of transformer core hysteresis loops and simplified
magnetization curve ............................................................................................ 28
Figure 1-2 Qualitative representation of voltage, flux and magnetizing current for a
transformer at steady state operation ................................................................... 28
Figure 1-3 Qualitative illustration of inrush phenomena and the effect of residual flux...... 29
Figure 1-4 Field measured long duration inrush current resulted from energising a 155
MVA GSU transformer [6] ................................................................................. 29
Figure 1-5 Measured RMS voltage dips caused by transformer energising at a 11 kV
distribution network [18] ..................................................................................... 32
Figure 1-6 Frequency of RMS voltage dip magnitudes out of 109 dip events measured at a
11 kV distribution network [18] .......................................................................... 32
Figure 1-7 RMS voltage dips caused by energising a 750/220/63 kV transformer (voltage
dips were measured on 220 kV side) [20] ........................................................... 32
Figure 2-1 Effect of circuit resistance during first cycle when switching in transformer at
the positive-going zero crossing of applied voltage [24] .................................... 38
Figure 2-2 Star-circuit representation of single-phase N-winding transformers [35] .......... 42
Figure 2-3 Connecting three two-winding STCs to represent a three-phase transformer .... 42
Figure 2-4 Schematic diagram of BCTRAN-based model for two-winding transformer, with
an externally connected core representation [35] ................................................ 44
Figure 2-5 Schematic diagram of the Hybrid transformer model [45] ................................. 44
Figure 2-6 Single-phase line with detail of a dx section ...................................................... 48
Figure 2-7 Pi-circuit model of a line [64] ............................................................................ 49
Figure 2-8 Equivalent two-port network for modelling a lossless line ................................ 51
Figure 2-9 Forming of Bergeron model based on two-port network model of lossless line 52
Figure 2-10 Statistical switching model involving closing time span among three phases [70]
............................................................................................................................. 54
Figure 2-11 Generic circuit for studying sympathetic interaction between transformers
connected in parallel............................................................................................ 57
Figure 2-12 Sympathetic inrush current waveforms simulated in [24] ................................ 58
Figure 2-13 One-line diagram of 20 kV converter test facility and recorded sympathetic
inrush current waveforms [84] ............................................................................ 63
Figure 2-14 Simplified electrical system circuit diagram [85] ............................................. 64
Figure 2-15 Measured voltage dips at 23 kV busbar [85] .................................................... 64
Figure 2-16 Simplified single-line diagram of wind farm collection grid during an emergent
islanded condition [2] .......................................................................................... 66
Figure 2-17 Variation of harmonic content of inrush current as a function of time [24] ..... 68
Figure 2-18 Field measured overvoltages caused by transformer energisation in HVDC
stations [13, 99] ................................................................................................... 70
Figure 2-19 System configuration, simulated harmonic resonant overvoltages and variation
of harmonic component [14] ............................................................................... 71
Figure 2-20 System configuration at the beginning of a restoration procedure and
overvoltage resulted from energising a transformer [96] .................................... 72
7
List of Figures
Figure 2-21 Size of voltage change against the time between each change [4] .................... 73
Figure 2-22 Simplified single-line diagram of a 138 kV BC Hydro system [19] ................. 74
Figure 2-23 Simplified diagram of a HV supply network in Australian system [51] ........... 75
Figure 2-24 Network configurations under comparison ....................................................... 77
Figure 2-25 Single line diagram of Jeju power system in Korea [101] ................................ 78
Figure 2-26 Voltage dip magnitudes resulted from different energisation angles when
residual flux is 28.3% and system loading is at its peak [101] ............................ 79
Figure 2-27 Schematic diagram of two wind farm configurations [75] ................................ 80
Figure 2-28 Wind farm topology and sequences for energising wind turbine transformers
[103] ..................................................................................................................... 81
Figure 2-29 Frequency of inrush current first peaks when residual flux and closing time
vary stochastically [104] ...................................................................................... 83
Figure 2-30 Network configuration studied in harmonic resonant overvoltages caused by
energising transformer during system restoration [43] ........................................ 83
Figure 3-1 South West Peninsula system as part of National Grid’s transmission system in
England & Wales ................................................................................................. 90
Figure 3-2 Schematic diagram of South West Peninsula system under detailed studies ...... 91
Figure 3-3 South West Peninsula voltage depression resulted from the first attempt........... 92
Figure 3-4 Three-phase inrush currents measured in the second attempt ............................. 93
Figure 3-5 RMS voltage dips measured at substation K in the second attempt .................... 93
Figure 3-6 Three-phase currents measured at power feeder 1 in Case E1 ............................ 95
Figure 3-7 Three-phase currents measured at the circuit I-K in Case E1 ............................. 96
Figure 3-8 Three-phase line-to-ground voltages measured at power feeder 1 in Case E1.... 97
Figure 3-9 Three-phase line-to-ground voltages measured at substation I in Case E1 ......... 97
Figure 3-10 Three-phase currents measured at power feeder 1 in Case E2 .......................... 98
Figure 3-11 Three-phase currents measured at circuit I-K in Case E2 ................................. 98
Figure 3-12 Three-phase line-to-ground voltages measured at power feeder 1 in Case E2.. 99
Figure 3-13 Three-phase line-to-ground voltages measured at substation I in Case E2 ....... 99
Figure 3-14 Three-phase currents measured at power feeder 1 in Case E3 (initial cycles) 100
Figure 3-15 Three-phase currents measured at power feeder 1 in Case E3 (long duration)
............................................................................................................................ 100
Figure 3-16 Three-phase currents measured at the circuit I-K in Case E3 ......................... 100
Figure 3-17 Three-phase line-to-ground voltages measured at power feeder 1 in Case E3 101
Figure 3-18 Three-phase RMS voltage dips measured at substation I in Case E3.............. 101
Figure 3-19 Sympathetic inrush currents measured at power feeder 1 in Case E4 ............. 102
Figure 3-20 RMS sympathetic inrush current measured at power feeder 1 in Case E4...... 103
Figure 3-21 RMS voltage dips measured at substation I in Case E4 .................................. 103
Figure 3-22 Currents measured at circuit I-K in Case E4 ................................................... 103
Figure 3-23 Basic tower structure used in South West system [113].................................. 105
Figure 3-24 Tower designs and transposing schemes associated with double circuit lines 106
Figure 3-25 Schematic diagram of SVC configuration ....................................................... 107
Figure 3-26 Procedure for generating firing pulses to control bi-directional thyristors...... 108
Figure 3-27 Conversion to derive saturation curve for type-96 nonlinear inductor ............ 110
Figure 3-28 lower half hysteresis curves for GSU transformers ........................................ 110
Figure 3-29 Open circuit test results deduced from GSU transformer model (415 MVA)
compared to manufacture test results ................................................................. 111
8
List of Figures
Figure 3-30 Open circuit test results deduced from GSU transformer model (345 MVA)
compared to manufacture test results ................................................................ 111
Figure 3-31 Comparison of inrush currents generated by Hybrid and BCTRAN+ for 345
MVA transformer (Energised at phase A voltage zero-crossing) ..................... 112
Figure 3-32 Comparison of inrush currents generated by Hybrid and BCTRAN+ for 415
MVA transformer (Energised at phase A voltage zero-crossing) ..................... 112
Figure 3-33 Comparison of inrush currents generated by Hybrid and BCTRAN+ for 415
MVA transformer (Energised at phase A voltage peak) ................................... 112
Figure 3-34 Conversion to derive saturation curve for type-93 nonlinear inductor ........... 113
Figure 3-35 Substation transformer saturation curves ........................................................ 114
Figure 3-36 Circuit diagram of CT model .......................................................................... 114
Figure 3-37 CT magnetization characteristic ..................................................................... 115
Figure 3-38 Simulated voltages at power feeder 1 compared to those measured in Case E1
........................................................................................................................... 116
Figure 3-39 Simulated currents at power feeder 1 compared to those measured in Case E1
........................................................................................................................... 116
Figure 3-40 Simulated currents at circuit I-K compared to those measured in Case E1 .... 117
Figure 3-41 Simulated voltages at substation I compared to those measured in Case E1.. 117
Figure 3-42 Simulated voltages at power feeder 1 compared to those measured in Case E2
........................................................................................................................... 117
Figure 3-43 Simulated currents at power feeder 1 compared to those measured in Case E2
........................................................................................................................... 118
Figure 3-44 Simulated currents at circuit I-K compared to those measured in Case E2 .... 118
Figure 3-45 Simulated voltages at substation I compared to those measured in Case E2.. 118
Figure 3-46 Simulated voltages at power feeder 1 compared to those measured in Case E3
........................................................................................................................... 119
Figure 3-47 Simulated currents at power feeder 1 compared to those measured in Case E3
........................................................................................................................... 119
Figure 3-48 Simulated currents at circuit I-K compared to those measured in Case E3 .... 119
Figure 3-49 Simulated RMS voltage variation at substation I compared to those measured in
Case E3 .............................................................................................................. 120
Figure 3-50 Comparison between measured and simulated inrush currents drawn by T2 and
T3 in Case E3 (simulated currents observed at the CT primary side) ............... 120
Figure 3-51 Comparison between measured and the simulated currents (observed at the CT
secondary side) .................................................................................................. 121
Figure 3-52 Simulated sympathetic inrush currents at power feeder 1 compared to those
measured in Case E4 ......................................................................................... 122
Figure 3-53 Simulated RMS sympathetic inrush currents at power feeder 1 compared to
those measured in Case E4 ................................................................................ 122
Figure 3-54 Simulated RMS voltage at power feeder 1 compared to those measured in Case
E4 ...................................................................................................................... 123
Figure 3-55 Simulated currents at circuit I-K compared to those measured in Case E4 .... 123
Figure 4-1 Inrush current observed at power feeder 2 (Case 5) ......................................... 128
Figure 4-2 Currents flowing through one of the circuits between substation I and K (Case 5)
........................................................................................................................... 128
9
List of Figures
Figure 4-3 Currents flowing through one of the circuits between substation J and K (Case 5)
............................................................................................................................ 129
Figure 4-4 Initiation of sympathetic inrush current observed at power feeder 1 (Case 5) .. 129
Figure 4-5 RMS sympathetic inrush current observed at power feeder 1 (Case 5) ............ 129
Figure 4-6 Voltage dips observed at Substation I (Case 5) ................................................. 129
Figure 4-7 Comparison between Case 5 and Case 4 regarding phase C voltage dip .......... 129
Figure 4-8 Inrush current observed at power feeder 1 (Case 10) ........................................ 130
Figure 4-9 Currents flowing through one of the circuits between I and K (Case 10) ......... 130
Figure 4-10 Currents flowing through one of the circuits between J and K (Case 10) ....... 131
Figure 4-11 Initiation of sympathetic inrush observed at power feeder 2 (Case 10) .......... 131
Figure 4-12 RMS sympathetic inrush current observed at power feeder 1 (Case 10) ........ 131
Figure 4-13 Comparison between Case 5 and Case 10 regarding phase C voltage dip ...... 131
Figure 4-14 Patterns of voltage dip magnitudes across all the network substations (voltage
dips observed at substation autotransformers’ 400 kV side versus 132 kV side)
............................................................................................................................ 132
Figure 4-15 Patterns of voltage dip duration across all the network substations (400 kV side
versus 132 kV side) ............................................................................................ 133
Figure 4-16 Voltage dip recovery traces observed at 400 and 132 kV busbars of substation I
............................................................................................................................ 133
Figure 4-17 Voltage dips caused by single-circuit outage between substation A and F ..... 134
Figure 4-18 Voltage dips caused by single-circuit outage between substation J and E ...... 134
Figure 4-19 Voltage dips caused by double-circuit outage between substation A and F.... 135
Figure 4-20 Voltage dips caused by double-circuit outage between substation B and C ... 136
Figure 4-21 Voltage dips caused by double-circuit outage between substation I and K .... 136
Figure 4-22 Voltage dips caused by double-circuit outage between substation J and E ..... 136
Figure 4-23 Voltage dips influenced by variation of key parameters ................................. 137
Figure 4-24 Example for illustrating two sensitivity factors Vd and Td............................. 138
Figure 4-25 Impacts of parameter variation on voltage dip magnitude .............................. 139
Figure 4-26 Impacts of parameter variation on voltage dip duration .................................. 139
Figure 4-27 Modified saturation curves for approximating maximum tap effect ............... 141
Figure 4-28 Voltage dip magnitudes observed in the case with GSU transformers set to their
maximum tap...................................................................................................... 142
Figure 4-29 Voltage dip duration observed in the case with GSU transformers set to their
maximum tap...................................................................................................... 142
Figure 4-30 Effect of SVC with different capacities on dip magnitude .............................. 143
Figure 4-31 Effects of SVC with different capacities on dip duration ................................ 143
Figure 4-32 Effects of SVC with different values of response time on dip duration .......... 144
Figure 4-33 Patterns of voltage dip duration at 400kV side for various SVC locations ..... 144
Figure 4-34 Effects of opening coupler CB1 on dip magnitude ......................................... 145
Figure 4-35 Effects of opening coupler CB1 on dip duration ............................................. 145
Figure 4-36 Dip magnitude pattern simulated under combined case .................................. 146
Figure 4-37 Dip duration pattern simulated under combined case...................................... 147
Figure 5-1 Procedure for generating stochastic circuit breaker closing time ...................... 153
Figure 5-2 Procedure for generating stochastic transformer core residual flux .................. 154
Figure 5-3 Single phase simulation circuit for preliminary Monte Carlo simulation ......... 155
10
List of Figures
Figure 5-4 Inrush current resulted from energising a single phase transformer under the
worst energisation condition ............................................................................. 155
Figure 5-5 Voltage dips resulted from energising a single phase transformer under the worst
energisation condition ....................................................................................... 156
Figure 5-6 Distribution of closing time in Case P1 ............................................................ 157
Figure 5-7 Distribution of residual flux in Case P1............................................................ 157
Figure 5-8 Relative voltage dip magnitudes plotted against relative inrush current peaks 158
Figure 5-9 Relative voltage dip magnitudes plotted against relative inrush current peaks 158
Figure 5-10 Frequency of voltage dips at different dip magnitude ranges ......................... 159
Figure 5-11 Frequency of voltage dips at different dip duration ranges ............................ 159
Figure 5-12 Relative inrush first peaks plotted against relative dip magnitudes ................ 160
Figure 5-13 Relative inrush first peaks plotted against relative dip durations ................... 160
Figure 5-14 Frequency plot of votlage dip magnitudes relative to the worst case dip
magnitude .......................................................................................................... 161
Figure 5-15 Frequency plot of votlage dip duration relative to the worst case dip magnitude
........................................................................................................................... 161
Figure 5-16 Distribution of offset closing time for three-phase poles in Case S1.............. 165
Figure 5-17 Distribution of residual flux in Case S1.......................................................... 165
Figure 5-18 Frequency of dip magnitude of each phase at substation I out of 1000 stochastic
runs .................................................................................................................... 165
Figure 5-19 Frequency of dip duration in each phase at substation I out of 1000 stochastic
runs .................................................................................................................... 166
Figure 5-20 Frequency of dip magnitude in each phase at substation I out of 5000 stochastic
runs .................................................................................................................... 166
Figure 5-21 Frequency of dip duration in each phase at substation I out of 5000 stochastic
runs .................................................................................................................... 166
Figure 5-22 Frequency of voltage dip magnitude in phase C at substation I under different
values of closing time span ............................................................................... 167
Figure 5-23 Uniform closing offset time distribution within ±2.5 ms range...................... 168
Figure 5-24 Exponential closing offset time distribution within ±2.5 ms range ................ 168
Figure 5-25 Frequency of voltage dip magnitudes in phase C at substation I for different
closing time span distributions .......................................................................... 169
Figure 5-26 Gaussian residual flux distribution ................................................................. 170
Figure 5-27 Exponential_1 residual flux distribution......................................................... 170
Figure 5-28 Exponential_2 residual flux distribution......................................................... 170
Figure 5-29 Frequency of voltage dips in phase C at substation I for different residual flux
distributions ....................................................................................................... 171
Figure 5-30 Frequency of voltage dips in phase C at substation I (comparing Case S9 with
S1) ..................................................................................................................... 172
Figure 5-31 Frequency of voltage dips in phase C at substation I of Case S10 contrasting
with that of Case S1........................................................................................... 173
Figure 5-32 Frequency of voltage dips in phase C at substation I of Case S10 contrasting
with that of Case S1 (two transformers with different residual flux) ................ 173
Figure 6-1 Layout of Nysted offshore wind farm collection grid and its connection with
onshore main grid .............................................................................................. 177
Figure 6-2 Measured three-phase voltages during energisation of feeder-A [128] ............ 178
11
List of Figures
Figure 6-3 Voltage waveforms around the energisation instants (at location P1) [128] ..... 178
Figure 6-4 Measured three-phase currents during energisation of feeder-A [128] ............. 179
Figure 6-5 Decay of phase A inrush current peaks measured at P1 .................................... 180
Figure 6-6 Complete network model of the Nysted wind farm collection grid and its
connection with the main grid ............................................................................ 181
Figure 6-7 Saturation curve of wind turbine transformers .................................................. 183
Figure 6-8 Comparison of voltage variation during energisation ....................................... 184
Figure 6-9 Comparison between measurement and simulation regarding the inrush currents
drawn by wind turbine transformer A9 .............................................................. 185
Figure 6-10 Comparison between measurement and simulation regarding the total inrush
currents drawn by nine wind turbine transformers in feeder-A ......................... 185
Figure 6-11 Decay trend comparison regarding feeder inrush currents .............................. 186
Figure 6-12 Decay trend comparisons regarding the inrush current drawn by transformer A1
and that drawn by A9 ......................................................................................... 186
Figure 6-13 Case studies of energising wind turbine transformers in one feeder ............... 188
Figure 6-14 Frequency of voltage dip magnitude in three phases at the point-of-common-
coupling under energising a feeder of wind turbine transformers ..................... 189
Figure 6-15 Frequency of voltage dip duration in three phases at the point-of-common-
coupling under energising a feeder of wind turbine transformers ..................... 189
Figure 6-16 Effects of transformer winding connections on voltage dip propagation ........ 190
Figure 6-17 Voltage dips caused by consecutive energisation under Case W4 .................. 191
Figure 6-18 Voltage dips caused by consecutive energisation under Case W5 .................. 191
Figure 6-19 Voltage dips caused by consecutive energisation under Case W6 .................. 192
Figure 6-20 Schematic diagram of two wind farm feeders connected at offshore platform
............................................................................................................................ 193
Figure 6-21 Sympathetic inrush currents drawn by the already connected feeder (1 km
electrical distance between two feeders) ............................................................ 193
Figure 6-22 Sympathetic inrush currents drawn by the already connected feeder (2 km
electrical distance between two feeders) ............................................................ 193
Figure 6-23 Sympathetic inrush currents drawn by the already connected feeder (3 km
electrical distance between two feeders) ............................................................ 193
Figure 6-24 Wind turbine transformers with different residual flux condition ................... 194
Figure 6-25 Sympathetic and inrush currents in the wind turbine transformers being
energised together .............................................................................................. 195
Figure 6-26 Sympathetic inrush currents observed in Case W6_1 simulation ................... 196
Figure 6-27 Sympathetic inrush currents observed in Case W6_2 simulation ................... 197
Figure 6-28 Definition of sympathetic inrush level ............................................................ 198
Figure 6-29 Procedure to correlate energisation sequence and sympathetic inrush level ... 200
Figure 6-30 Accumulated sympathetic inrush level on each wind turbine transformer
resulted from deterministic testing of S1 ........................................................... 202
Figure 6-31 Accumulated sympathetic inrush level on each wind turbine transformer
resulted from deterministic testing of S2 ........................................................... 202
Figure 6-32 Accumulated sympathetic inrush level on each wind turbine transformer
resulted from deterministic testing of S3 ........................................................... 203
Figure 6-33 Accumulated sympathetic inrush level on each wind turbine transformer
resulted from deterministic testing of S4 ........................................................... 203
12
List of Figures
Figure 6-34 Accumulated sympathetic inrush level of each wind turbine transformer
resulted from stochastic testing of S1................................................................ 204
Figure 6-35 Accumulated sympathetic inrush level of each wind turbine transformer
resulted from stochastic testing of S2................................................................ 205
Figure 6-36 Accumulated sympathetic inrush level of each wind turbine transformer
resulted from stochastic testing of S3................................................................ 205
Figure 6-37 Accumulated sympathetic inrush level of each wind turbine transformer
resulted from stochastic testing of S4................................................................ 205
13
List of Tables
List of Tables
Table 2-1 Origin of electrical transients and their associated frequency ranges [33]........... 41
Table 2-2 Summary of previous contributions on converting RMS V/I to λ/i curve ........... 46
Table 2-3 Guidelines for modelling circuit breaker [69] ...................................................... 53
Table 2-4 Summary of case studies carried out in [2] .......................................................... 67
Table 2-5 Estimated inrush current peaks, duration and voltage dip magnitudes resulted
from the worst case energisation under different network configurations [51] .. 77
Table 2-6 Energisation condition for simulation assessment [101]...................................... 78
Table 2-7 Summary of the influential parameters ................................................................ 82
Table 2-8 Cases studies of the influences of stochastic variables on the inrush current of a
single-phase transformer ..................................................................................... 82
Table 3-1 Plan of new generation installations at the South West Peninsula system [111] . 90
Table 3-2 Four energisation cases in the further field measurement .................................... 94
Table 3-3 Line dimension and conductor data [113] .......................................................... 105
Table 3-4 System loading data of the South West system.................................................. 106
Table 3-5 GSU transformer test report (T1&T2, 345 MVA) ............................................. 109
Table 3-6 GSU transformer test report (T3, 415 MVA) ..................................................... 109
Table 3-7 Comparison of three-phase voltage dip magnitudes .......................................... 123
Table 4-1 Voltage dips observed at substation I under different energisation conditions.. 126
Table 5-1 Case studies of stochastic estimation of voltage dips caused by energising a single
phase transformer .............................................................................................. 156
Table 5-2 List of case studies conducted in stochastic assessment .................................... 163
Table 5-3 Possible ranges and PDFs for random parameters ............................................. 172
Table 6-1 List of top ten operational offshore wind farms [125] ....................................... 176
Table 6-2 Inrush current first peaks resulted from energisation of feeder-A ..................... 180
Table 6-3 132 kV single core onshore cables [132] ........................................................... 182
Table 6-4 132 and 33 kV three-core offshore cables [133] ................................................ 182
Table 6-5 Main electrical information for modelling wind turbine transformers [128] ..... 183
Table 6-6 Network model parameter settings for simulating field measurement results ... 184
Table 6-7 Voltage dips resulted from energising a feeder of wind turbine transformer under
the strong source strength .................................................................................. 188
Table 6-8 Voltage dips resulted from energising a feeder of wind turbine transformer under
the weak source strength ................................................................................... 188
Table 6-9 Sequences for energising wind turbine transformers in a feeder ....................... 199
Table 6-10 Aggregation of sympathetic inrush levels resulted from each energisation
sequence ............................................................................................................ 202
15
List of Abbreviations
List of Abbreviations
CB Circuit Breaker
CT Current Transformer
DC Direct Current
FC Fixed Capacitor
HV High Voltage
LV Low Voltage
17
Abstract
Abstract
Transformers are essential components facilitating transmission and distribution of electric power.
Energisation of transformers, however, can cause core operating at deep saturation region and thereby
induce transient inrush currents of high magnitude and with rich harmonics. This can lead to undesirable
effects including potential damage to the transformer itself, relay mal-operation, harmonic resonant
overvoltages, and reduced power quality in the system (mainly in the form of voltage dips).
This thesis investigates voltage dips caused by energising generator step-up (GSU) transformers and two
types of generation connection are studied: one is a combine cycle gas turbine (CCGT) plant connected to
a 400 kV transmission grid and the other is a large offshore wind farm connected to a 132 kV distribution
grid. To carry out the investigation, detailed network models were developed in alternative transients
program/electromagnetic transients program (ATP/EMTP) and validated with the help of field
measurements.
For the connection of generation in the transmission grid, deterministic assessment was conducted to
comparatively analyse voltage dips caused by energising large GSU transformers under different
energisation conditions and different network conditions; special attention was paid to the energisation
cases involving sympathetic inrush between transformers by addressing its prolonging effects on voltage
dips, with sensitivity studies further carried out to identify the key influential parameters. In addition,
stochastic assessment was conducted by applying Monte Carlo method, which helps identify the dip
frequency pattern and the likelihood of reaching the dip magnitude resulted from the commonly agreed
worst case energisation condition; their sensitivities to the variation of circuit breaker closing time span,
transformer core residual flux, system condition and the number of transformers being energized together
were also investigated. Furthermore, possible cost-effective operational approaches to mitigate the
voltage dips were explored and compared. For the connection of large offshore wind farm, voltage dips
caused by energising wind turbine transformers under different scenarios were assessed; in particular,
sympathetic inrush between wind turbine transformers were studied, and the energisation sequence
resulting in less sympathetic inrush was deterministically identified and stochastically validated.
The simulation results of deterministic studies indicate that, when carrying out energisation of a large
GSU transformer in the transmission grid under the commonly agreed worst case energisation condition,
the dip magnitude can reach 9.6% and the duration 2.7 seconds; moreover, when coupled with
sympathetic inrush, the duration can be prolonged by 136%, lasting for 6.4 seconds. The sensitivity
studies show that transformer core saturation inductance is the key parameter determining dip magnitude
and transformer copper losses is the key parameter determining dip duration. Stochastic assessment of
voltage dips shows that, out of 1000 stochastic dip events, less than 0.5% of the dips can reach the worst
case dip magnitude and about 80% are of magnitudes less than 0.6 pu of the worst case dip magnitude;
the dip frequency pattern is found to be insensitive to the circuit breaker closing time variation but can be
considerably influenced by the residual flux distribution. In terms of mitigation measures, it was proven
that, by adjusting tap changer position, applying static var compensator and even opening coupler circuit
breaker in the substation, the degree of voltage dip especially the dip duration can be significantly reduced.
Contrasting to those observed in the transmission grid, voltage dips resulted from energising wind turbine
transformers in large offshore wind farms are of less concern; dip magnitudes are no more than 1% in the
case of energising a stand-alone wind turbine transformer. However, sympathetic inrush between wind
turbine transformers within one feeder was found to be significant and the energisation sequence resulting
in less sympathetic inrush is to separately energise the wind turbine transformer from the one closest to
the offshore platform to the one farthest away from the platform.
19
Declaration
Declaration
I declare that no portion of the work referred to in the thesis has been submitted in
support of an application for another degree or qualification of this or any other
university or other institute of learning.
21
Copyright Statement
Copyright Statement
(I). The author of this thesis (including any appendices and/or schedules to this thesis)
owns certain copyright or related rights in it (the “Copyright”) and s/he has given The
University of Manchester certain rights to use such Copyright, including for
administrative purposes.
(II). Copies of this thesis, either in full or in extracts and whether in hard or electronic
copy, may be made only in accordance with the Copyright, Designs and Patents Act
1988 (as amended) and regulations issued under it or, where appropriate, in accordance
with licensing agreements which the University has from time to time. This page must
form part of any such copies made.
(III). The ownership of certain Copyright, patents, designs, trademarks and other
intellectual property (the “Intellectual Property”) and any reproductions of copyright
works in the thesis, for example graphs and tables (“Reproductions”), which may be
described in this thesis, may not be owned by the author and may be owned by third
parties. Such Intellectual Property and Reproductions cannot and must not be made
available for use without the prior written permission of the owner(s) of the relevant
Intellectual Property and/or Reproductions.
(IV). Further information on the conditions under which disclosure, publication and
commercialisation of this thesis, the Copyright and any Intellectual Property and/or
Reproductions described in it may take place is available in the University IP Policy
(see http://www.campus.manchester.ac.uk/medialibrary/policies/intellectual-
property.pdf), in any relevant Thesis restriction declarations deposited in the University
Library, The University Library’s regulations (see
http://www.manchester.ac.uk/library/aboutus/regulations) and in The University’s
policy on presentation of Theses.
23
Acknowledgement
Acknowledgement
First and foremost I would like to thank my supervisors Dr. Haiyu Li and Prof.
Zhongdong Wang.
It is my privilege to be their student and this thesis would not have been possible
without their invaluable guidance, generous support and constant inspiration. I much
appreciate Dr. Haiyu Li for his generous supervision, forward-looking advices and
insightful discussions during my PhD studies. I am deeply grateful to Prof. Zhongdong
Wang for providing me the great opportunity to explore and endeavour the PhD journey
and the intensive technical guidance; she helped me build up confidence, cultivate
academic skills and enlighten my vision of future career.
I would like to thank Dr. Foroozan Ghassemi and Prof. Paul Jarman of National Grid,
for providing technical support.
Special thanks also to all my colleagues in the transformer research group, in Ferranti
Building and our school for their company and support; in particular, I would like to
express my gratitude to Dr. Swee Peng Ang for his kind help and interesting discussions
during my MSc and PhD studies in the University of Manchester.
With extreme appreciation, I wish to convey my sincere thanks to my family for their
continuous care and encouragement. Especially, I am deeply indebted to my beloved
wife Xiao Yi whose selfless love and devotion continuously motivate me to overcome
difficulties, make progress and achieve my best.
25
Chapter 1 Introduction
Chapter 1 Introduction
1.1 Background
Modern society critically relies on electric power as the key energy source and constant
efforts have been made by power system operators to maintain and operate
interconnected electrical systems as reliably as possible.
One of the challenges for the quality of power supply is the disturbances caused by
transformer energisation. Due to the nonlinearity of the magnetic characteristic of
transformer core, transformer energisation would result in inrush currents of high
magnitude and with rich harmonics, causing damage to transformer itself, and
influencing the system by harmonic resonant overvoltages, relay mal-operation and
reduced power quality mainly in terms of voltage dips [1]. Indeed, a modern high
voltage transmission grid normally consists of hundreds of power transformers and a
distribution grid may consist of thousands; the topology of some future network
configurations, such as in the case of wind farm grids, reflects the tendency that power
transformers would be connected more adjacent to one another and it is more likely to
simultaneously switch on a group of transformers [2, 3]. Therefore, more intensified
inrush transient interactions could occur, which would cause adverse impacts on power
system, hence affect industrial and commercial customers. This, without proper
management, could lead to significant economic losses; and consequently their
associated adverse effects should be carefully assessed to guide system operation and
planning so as to ensure the compliance with tightening standards that define secure and
high quality supply of electric power [4, 5].
27
Chapter 1 Introduction
obtained, as shown in Figure 1-1 (b). Normally, transformer core operates at the linear
region where the core magnetic permeability is high and the core magnetizing current is
low. As the voltage is increased, more and more flux is demanded (the flux is
proportional with the integral of the applied voltage) and the core would enter the
saturation region where a slight increase of flux would result in a significant increase of
magnetizing current.
Flux
Current
Figure 1-2 shows the relationship between voltage, magnetic flux, and magnetizing
current for a transformer under steady state operation. Due to the transformer core non-
linearity, the magnetizing current is non-sinusoidal. It follows the hysteresis loop
oscillating between ±im as the flux changes sinusoidally between ±Φm (the magnitude of
im is normally between 0.5% and 2% of the transformer rated current).
Figure 1-2 Qualitative representation of voltage, flux and magnetizing current for a transformer
at steady state operation
28
Chapter 1 Introduction
when the voltage peak is reached and continue to increase to 2Φm when the voltage
returns to zero again (as shown by the dash blue line in Figure 1-3). The excessive
demand of flux would saturate the core and result in sharp increase of the magnetizing
current. This sharply increased current is termed as inrush current (iinrush) and its
magnitude could be many times of the transformer nominal magnetizing current.
Usually, after a transformer is switched off, magnetizing current will follow a hysteresis
loop to zero and some residual flux Φr could be retained in the core. If this residual flux
is taken into account and suppose its polarity is in the direction of flux build-up, the
maximum flux resulted from the above-mentioned energisiation would become 2Φm+Φr
(as indicated by the solid blue line in Figure 1-3), resulting in even larger inrush current.
Time
Figure 1-3 Qualitative illustration of inrush phenomena and the effect of residual flux
Figure 1-4 shows a measured inrush current resulting from energising a 155 MVA
generator step-up (GSU) transformer [6]. The magnitude of the first peak of inrush
current is hundreds of times larger than the magnitude of the nominal magnetizing
current. The waveform of the inrush current looks like half cycle sinusoidal wave
superimposed by a DC component whose decay is largely influenced by system losses.
The inrush current is rich in harmonics especially even harmonics [7].
Figure 1-4 Field measured long duration inrush current resulted from energising a 155 MVA
GSU transformer [6]
29
Chapter 1 Introduction
The amplitude of inrush current can be equal to that of the short-circuit current [8], and
it lasts longer without enough damping in the system and consequently can seriously
damage the windings through excessive mechanical stresses. In fact, the axial forces due
to inrush current are always larger than those caused by short circuit current and the
radial force applied on transformer high voltage winding can be three times of the
corresponding force caused by short circuit condition [9]. These adverse effects on the
windings can result in pre-mature failure of a transformer.
Transformer inrush currents are rich in harmonics [10]. If one of the harmonic
components in the inrush current is close to the resonant frequency of the power system,
a sustained overvoltage might be produced [11]. This can be encountered in the
following scenarios:
• After system collapse, a black-start process is carried out by energising a remote
transformer against minimum generation and loading [12];
• In High-Voltage Direct Current (HVDC) scheme, the ac system can be in
resonance with the ac harmonic filters at particular harmonic frequencies, and
these resonances can be excited by inrush currents resulted from energising
HVDC converter transformers [13];
• In offshore electrical systems, such as offshore wind farms or oil production
facilities, the interconnection via a subsea cable introduces a significant shunt
capacitance to the source power system, resulting in low resonant frequency, and
therefore the resonant overvoltages may be excited by one of the harmonics of
the inrush currents caused by energising wind turbine transformers or offshore
platform transformers [14].
Harmonic resonant overvoltages may be amplified above the level the system
equipment can withstand and if these overvoltages last for a long period of time, they
may eventually damage the equipment [15].
30
Chapter 1 Introduction
Relay mal-operation
Transformer inrush may cause mal-operation of transformer relay protection [7]. During
transformer energisation, the inrush current only flows through the energised winding,
which has no equivalent currents from the other windings. This unbalance condition
may be treated by the transformer differential protection as transformer internal faults,
thereby tripping the circuit breaker immediately after transformer energisation.
Although the second harmonic of the inrush current can be used to identify the inrush
condition and restrain the relay operation, the restraining criteria for setting the relay
may affect the relay performance [16]; and in certain cases, ultra-saturation may be
induced by transformer energisation, which would inevitably cause tripping of a healthy
transformer [17].
Voltage dips
Due to the impedance between the supply source and the energised transformer, the
inrush currents may result in temporary voltage dips in the connected system. The dips
differ in the magnitudes among three phases and take a long time to recover [18]. If the
short circuit level at the transformer busbar is low, the resulted voltage dips can be
significant which may affect downstream customers having devices sensitive to power
quality variation [19].
Similar voltage dip events were observed in high voltage transmission networks. For
example, in [20], voltage dips caused by energising a no-load 750/220/63 kV auto-
31
Chapter 1 Introduction
transformer were measured on 220 kV side, which are shown in Figure 1-7. Compared
to those observed in distribution network shown in Figure 1-5, the dip pattern is similar
but the recovery is slower (the dips took about 750 ms to recover by 50%). This slower
recovery could be attributed to the larger L/R ratio in the transmission system, as the
L/R ratio determines the decaying time constant of the inrush current.
100
Number of events
50
0
0 5 10 15 20 25 30 35
Voltage dip (%)
Figure 1-6 Frequency of RMS voltage dip magnitudes out of 109 dip events measured at a 11
kV distribution network [18]
Figure 1-7 RMS voltage dips caused by energising a 750/220/63 kV transformer (voltage dips
were measured on 220 kV side) [20]
32
Chapter 1 Introduction
Recently, connections of wind farms are increasing and they are often located at remote
areas where the electrical network can be of relatively low fault level (i.e. the source
strength is weak). Voltage dips caused by energising wind turbine transformers have
been causing concerns.
Sympathetic inrush
Through the investigation, it is aimed to answer the following questions that have not
been addressed before:
Influence of sympathetic inrush on voltage dips caused by transformer
energisation;
Probability of encountering the worst case voltage dips;
Energisation sequence resulting in less sympathetic inrush between wind turbine
transformers.
To achieve the objectives, the scope of the work covers the following areas:
Use Alternative Transients Program/Electro-Magnetic Transients Program
(ATP/EMTP) to develop network models suitable for studying voltage dips and
sympathetic inrush caused by energising transformers;
Assess and compare voltage dips caused by energising GSU transformers in a
33
Chapter 1 Introduction
Chapter 1 Introduction
This chapter presents a general background about transformer inrush and its potential
adverse impacts on the power transformer itself and the power system. The objectives
and the scope of work of this thesis are presented.
This chapter reports the field measurements of inrush currents, sympathetic inrush
currents and voltage dips caused by energising 400 kV GSU transformers. The network
model development is described in detail and also its validation against field
measurement results.
34
Chapter 1 Introduction
This chapter describes the comprehensive assessment on voltage dips in the 400 kV
system caused by energising large GSU transformers, including: comparison of voltage
dips under different energisation conditions; the network-wide voltage dips under both
non-outage and outage scenarios; the influence of sympathetic inrush on voltage dips. It
also presents the work that has been done on identification of key influential parameters
on voltage dips and exploring operational approaches to cost-effectively reduce voltage
dips and sympathetic inrush.
This chapter attempts to extend the deterministic studies carried out in Chapter 4 by
taking into account stochastic variables. First, an ATP-EMTP interface to facilitate
stochastic assessment using Monte Carlo method is described. Second, the possible
stochastic variables are discussed and quantified. Then a preliminary stochastic
simulation based on a single phase circuit is presented. Finally, stochastic studies of the
three-phase system are presented, which include: calculating the distribution of voltage
dip magnitudes and durations; identifying the probability of the worst case voltage dips;
and testing the sensitivity of the results to various closing time span and residual flux
distributions.
This chapter applies the modelling and simulation methods proven in Chapter 4 and
Chapter 5 to assess transformer inrush transients during offshore wind farm connections.
Special attention is focusing on voltage dips and sympathetic inrush between wind
turbine transformers. It first describes the development of a wind farm collection grid
model and its validation against measurements, and then presents the studies of voltage
dips and sympathetic inrush caused by energising wind turbine transformers in a large
offshore wind farm collection grid. The studies estimate the possible voltage dips
caused by energising wind turbine transformers and assess energisation sequences with
the aim to reduce sympathetic inrush between wind turbine transformers.
35
Chapter 1 Introduction
This chapter summarizes the main finding of this thesis work. Future work is also
suggested for several aspects of the research on voltage dips and sympathetic inrush
caused by transformer energisation.
36
Chapter 2 Literature Review on Transformer Energisation Transients
Transients resulted from transformer energisation were first observed by Ferranti when
commissioning the Deptford to London 11 kV link in 1890 [23]. Afterwards, abundant
publications were devoted to the calculation of transformer inrush current, the
assessment of transformer energisation transients and the development of possible
mitigation measures.
37
Chapter 2 Literature Review on Transformer Energisation Transients
· 2
(2.1)
where
is the magnetic flux density outside the saturated core,
is the peak
nominal flux density,
and
are the residual and saturation flux densities, is the
length of the magnetic flux path in air, is the turn number of the energised winding,
is the cross-section of the space enclosed by the energised winding, is the cross-
section of the iron core and is the permeability of the air.
This formula, however, only considers one particular switching time (voltage zero
crossing) and assumes infinite short circuit capacity at the transformer terminal;
moreover, it neglects the resistance of the energisation circuit which contributes to the
decaying mechanism of inrush current (see Figure 2-1, the circuit resistance would
reduce the flux by an amount of per cycle, resulting in the attenuation of
!"
Figure 2-1 Effect of circuit resistance during first cycle when switching in transformer at the
positive-going zero crossing of applied voltage [24]
As one step forward, the work presented in [25] suggests an analytical formula which
takes into account the effect of switching angle and circuit resistance to predict the first
peak of inrush current. The formula was derived as follows:
#$
) *+,- 1/
% ! &'
( !
(2.2)
with additional parameters: #$ the magnitude of the applied voltage, & the angular
frequency, - the initial phase angle of the voltage source, the series resistance and
'
( the air-core inductance of the energised winding.
38
Chapter 2 Literature Review on Transformer Energisation Transients
Analytical equations proposed in [26], [27] and [28] extended the estimation of inrush
current peak from the first cycle to the following cycles. The work presented in [26]
estimates inrush current peaks via the following procedure:
Step 1 - calculate the saturation angle -
(i.e., the angle at which saturation occurs):
-
cos3 ) /
(2.3)
√2#$
1 cos -
&'
(
(2.4)
567
+
· · 2,5-
-
*+,-
&'
(
(2.5)
With the updated
, steps 1, 2 and 3 are repeated to calculate the inrush peaks of
subsequent cycles.
The method given in [26] neglects the integrated term for the first cycle and only
gives the peak values of the inrush current rather than the full inrush current waveform.
These limitations were overcome by including an exponential transient term in the
equations given by [27] and [28], which helps obtain the full current waveform as a
function of time, as shown by the following equation:
√2#$
% ! &'
( !
: C
· 8,5& 9 6 · ,5-
9G
D<E
(2.6)
;<=>?@A>B F
It should be noted that the above analytical approaches can only estimate the inrush
current peaks for a single-phase transformer. It might not be suitable to apply them for
estimating inrush transient of multi-winding and multi-phase transformers or assessing
transformer energisation transient with other network components involved.
39
Chapter 2 Literature Review on Transformer Energisation Transients
Some early attempts, such as those in [10, 29, 30], were focusing on numerical
prediction of inrush current in single-phase transformers. Basically, they adopted a time
stepping technique which can give successive discrete values of current at successive
chosen steps. The main purpose of these attempts was to investigate the effects of
varying point-on-wave switching, residual flux and transformer resistance on the
harmonic content variation of inrush current. However, they were not extended to study
the impacts of inrush transients on system operation in large-scale networks.
Inrush current calculation based on finite element method (FEM) was presented in [31].
In this method, both magnetic field and electrical circuit equations are solved
simultaneously. Similar contribution can be found in [32]. The use of FEM calculation
allows investigating mechanical stresses on the winding, internal flux distribution and
thermal condition. However, it is time consuming, computationally costly and not
suitable for studying inrush currents’ network impacts.
Up to date, the most frequently used numerical tool for calculating inrush current is the
Electromagnetic Transient Program (EMTP). In the EMTP-type simulation packages
(including ATP/EMTP and PSCAD/EMTDC), standard available transformer models
are provided together with non-linear elements describing transformer core hysteresis
and saturation features; models for other system components, including transmission
line, cable, circuit breaker, surge arrester, rotating machines and Flexible AC
Transmission Systems devices, are also provided. This enables the simulation of
complex networks and control systems of arbitrary structure, which make the EMTP a
preferable numerical tool for calculating inrush currents as well as simulating their
impacts on network operation. In the following section, modelling system components
in EMTP for transformer energisation transient studies are reviewed.
40
Chapter 2 Literature Review on Transformer Energisation Transients
In Table 2-1, various origins of transient and their associated frequency ranges are listed.
Table 2-1 Origin of electrical transients and their associated frequency ranges [33]
Origin of transients Frequency range
Transformer energisation
(DC) 0.1 Hz – 1 kHz
Ferroresonance
Load rejection 0.1 Hz – 3 kHz
Fault clearing 50/60 Hz – 3 kHz
Fault initiation 50/60 Hz – 20 kHz
Line energisation 50/60 Hz – 20 kHz
Line reclosing (DC) 50/60 Hz – 20 kHz
Transient recovery voltage
50/60 Hz – 20 kHz
Terminal faults
Short line faults 50/60 Hz – 100 kHz
Multiple re-strikes of circuit breaker 10 kHz – 1 MHz
Lightning surges, faults in substations 10 kHz – 3 MHz
Disconnector switching (single re-strike)
100 kHz – 50 MHz
Faults in GIS
Since the transformer energisation transients mainly range between DC to 1 kHz, when
modelling system components in EMTP to study transformer energisation transients, the
frequency range can be targeted on the range between DC and 1 kHz [34].
41
Chapter 2 Literature Review on Transformer Energisation Transients
L M N'O3 NPO N'O3 NONO
(2.7)
Core saturation and hysteresis effects are modelled by a nonlinear inductor connected at
the star point which can be located at the primary winding side.
i1 i2
in
The use of three single-phase two-winding STCs to model a three-phase transformer has
been applied in many studies [19, 36-39]. As shown in Figure 2-3, three single-phase
two-winding STCs were used to model a 315 MVA 138/21 kV, YNd connected GSU
transformer, with their primary side connected as grounded-star and secondary side
delta [19].
42
Chapter 2 Literature Review on Transformer Energisation Transients
These impedances can be calculated from positive and zero sequence impedances (Q3
and Q )
Q 2Q3
Q
3
(2.10)
Q Q3
Q$
3
(2.11)
being NO and N'O the real and imaginary part of the branch impedance matrix,
respectively.
In the case of a very low excitation current, the transformer should be described by an
admittance formulation
NRO NWON#O (2.13)
The BCTRAN model takes phase-to-phase coupling into account; it can model core
saturation but does not consider core topology. It is linear and is reasonably accurate for
frequencies below 1 kHz [41]. The effects of core saturation can be represented by a set
of externally connected non-linear inductances which are usually added at the terminals
43
Chapter 2 Literature Review on Transformer Energisation Transients
of the transformer winding nearest to the core, as shown in Figure 2-4. Application of
BCTRAN with external core representation has shown satisfactory performance for
studying transformer energisation transients in various networks [6, 42-44].
Figure 2-4 Schematic diagram of BCTRAN-based model for two-winding transformer, with an
externally connected core representation [35]
Capacitive
coupling
Winding
connection
Frequency dependent
winding resistance
Leakage equivalent
Zl Zl Zl
Core equivalent
L0 L0
Zy Zy
The equivalent circuit for the core modelling is derived from a simplified core magnetic
circuit via duality transformation (meshes in the magnetic circuit are transformed to
nodes in the electrical dual; reluctances are changed to inductances; sources of
44
Chapter 2 Literature Review on Transformer Energisation Transients
magnetomotive force become current sources) [46, 47]; in the circuit, a constant
resistance in parallel with a nonlinear inductance is used to represent each limb (Zl) and
yoke (Zy); L0 accounts for the zero-sequence flux paths. The leakage equivalent is an
inverse inductance matrix established based on BCTRAN approach.
Hybrid model has been implemented in ATPDraw [48]. Recent development of the
Hybrid model suggests the inclusion of type-96 non-linear inductors to produce residual
flux after de-energisation and accommodation of user-defined air-core inductance value
[49]. However, due to the requirement of core design data and the difficulty of manually
initializing residual flux, so far this model has not been widely used.
They are different from each other in terms of conversion approach, consideration of
core losses and handling of delta connection. Taking core losses into account improves
the accuracy of the conversion; handling of delta connection counts the fact that, in the
45
Chapter 2 Literature Review on Transformer Energisation Transients
tests carried out with closed delta-coupled windings, the triplen harmonics circulating in
the closed delta do not appear in the measured line currents [57]. The analytical
approach reported in [54] is the basis of the current main conversion routine (called
SATURA) used in ATP/EMTP [58].
Table 2-2 Summary of previous contributions on converting RMS V/I to λ/i curve
Contributors Conversion method Core losses Delta connection
Talukdar et al. [53] Numerical Not considered Not considered
Prusty and Rao [54] Analytical Not considered Not considered
Neves and Dommel [55] Analytical Considered Not considered
Neves and Dommel [57] Numerical Considered Considered
Chiesa and Høidalen [56] Analytical Considered Considered
In the open-circuit tests carried out by transformer manufacturer, the commonly applied
excitation levels are 90%, 100% and 110% of rated operating voltage, because, during
operation, system voltage variations are normally controlled to be within ±10% of the
nominal operation voltage. Hence, the converted λ/i characteristic usually consists of
only three points. These points can only form a very crude piecewise nonlinear core
saturation curve. Extension of the test report data is needed to form a more complete
core saturation curve.
There are two approaches for data extension: one is linear extrapolation and the other is
curve fitting. In linear extrapolation, the final segment of the crude piecewise nonlinear
curve is linearly extended to form a constant slope for representing λ/i characteristic in
transformer core deep saturation region. Although this method is simple, it may
severely underestimate the current resulted from any excitation level above the 110%
excitation level, because at the 110% excitation level the core has not reached the deep
saturation. Instead of simple linear extrapolation, curve fitting approach generates new
artificial points which form new segments to be added into the crude piecewise
nonlinear curve. One commonly used curve fitting function is a two term nth order
polynomial function:
X · Y Z · Y (2.15)
This function was first identified in [59] and used to develop an analytical approach to
evaluate ferroresnance. It was applied to curve fit non-linear saturation characteristics of
potential transformer and substation transformer in [60]. Other more complex functions
but rarely used for the curve fitting were documented in [61].
46
Chapter 2 Literature Review on Transformer Energisation Transients
The value of air-core inductance is important for the estimation of transformer core
saturation inductance which determines the final slope of core saturation curve [52]. In
BCTRAN and STC, the saturation inductance Lsat is deduced by taking into account
short circuit inductance:
'
'
( '^; (2.16)
In the case of Hybrid model, the leakage flux between the inner winding and the core
are taken into account by LLC; in this case, Lsat is calculated by:
'
'
( '^; ';_ (2.17)
47
Chapter 2 Literature Review on Transformer Energisation Transients
aR`, &
cf& d&g&e#`, & W&#`, &
a`
(2.19)
where #`, & and R`, & are the voltage and current of the line, respectively; &,
'& , f& and g& are line parameters expressed in per unit length and with
frequency dependent; Q& [=& '&] and W& [=f& g&] represent line
series impedance and shunt admittance in per unit length, respectively.
& d&'&
Q & h
f& d&g&
(2.20)
From the line equations above, one can further obtain the well-known relation between
the sending and receiving end [63]:
These equations are for describing steady-state conditions and form the basis for
deriving transient line models.
For modelling lines in time-domain simulation, two types of model are commonly used:
• Lumped-parameter models, usually known as pi-models, represent line by
48
Chapter 2 Literature Review on Transformer Energisation Transients
The exact pi-model is derived from equation 2.22 and can be described by the
equivalent circuit shown in Figure 2-7. In the circuit,
sinh i&
Q3 & Q& ·
i& (2.23)
tanh i&/2
W! & W& ·
i&/2 (2.24)
The scalar Q and W can be replaced by corresponding NQO and NWO matrixes to describe
N-phase transmission line. It is an exact representation of the line at a given frequency
and therefore named as exact-pi model. This model is suitable for steady-state or
harmonic analysis in which solutions are obtained for one frequency at a time. However,
they are difficult to be applied for time-domain transient analysis because the elements
in the exact pi-model are two-fold frequency dependent: one due to the line parameters
themselves (including Q and W) and second due to the propagation constant i . Even by
assuming the line parameters as constant, the elements in the exact pi-model are still
function of frequency.
49
Chapter 2 Literature Review on Transformer Energisation Transients
Nominal pi-model is described by the same equivalent circuit of exact pi-model (as
shown in Figure 2-7) but with the branch total series impedance and total shunt
admittance equal to Q& · and W& · , with the hyperbolic correction factors
neglected. The nominal pi-model is preferred over the exact pi-model because it can be
directly applied to time-domain transient simulations by assuming the line parameters as
constant. In addition, it approximates the performance of exact pi-model. In the cases of
short lines or low frequencies, the nominal pi-model is effectively identical to the exact
pi-model and it may be used for transient simulations in the proximity of the frequency
at which line parameter values are calculated. In the cases of long lines and high
frequencies, a number of cascaded short nominal pi-model sections, the so called
“cascaded nominal pi-model”, can be used to approximate the frequency dependent
effect of the propagation constant so as to mimic the performance of exact pi-model [12,
36]. However, the treatment of lumped parameters can give rise to spurious oscillations
and hence the pi-model is not preferable for representing the frequency dependent line
parameters.
Bergeron model is a simple, constant frequency model, based on travelling wave theory.
In general, it is a combination of lossless distributed parameter line and lumped series
resistances [65].
One important technique utilized to derive the Bergeron model (as well as the frequency
dependent model) is the decoupling between the line sending and receiving ends. To
understand this, a lossless distributed parameter line can be considered first. The general
solution for the wave propagation equations of a lossless line are:
'
Q h
g
(2.27)
1
p
√'g
(2.28)
50
Chapter 2 Literature Review on Transformer Energisation Transients
Clearly, when `– p is constant, r`, Q `, is constant. Given as the
length of the line, the travelling time for a constant wave to travel from the end k to the
other end m of the line is:
s √'g
p
(2.31)
Hence
1
$ P R$ s
Q $
(2.33)
1
where
R$ s P s $ s
Q
(2.34)
Similarly
1
$ P R s
Q
(2.35)
1
where
R s P s $ s
Q $
(2.36)
51
Chapter 2 Literature Review on Transformer Energisation Transients
As can be seen, the terminals of the two port network are topologically disconnected,
i.e., the sending and receiving ends of the line are effectively decoupled from each other
during the solution at time t. This is valid in time domain simulation provided that the
simulation time step, ∆t, is smaller than travelling time of the waves, τ [66].
The Bergeron model is formed by adding lumped series resistances into the lossless line
model, with the conductance to ground neglected. This is made by splitting the total line
resistance into three lumped parts and locating them at the middle and at the ends of the
line, as shown in Figure 2-9. The error incurred in lumping the series resistance as
compared to the distributed case is acceptable as long as · t Q [66].
Figure 2-9 Forming of Bergeron model based on two-port network model of lossless line
52
Chapter 2 Literature Review on Transformer Energisation Transients
Mechanical
Important Very important Negligible Negligible
pole spread
Closing
Prestrikes Negligible Important Important Very important
High current Important only for interruption
Negligible Negligible
interruption capability studies
Current Important only for interruption of
Negligible Negligible
chopping small inductive currents
Restrike
Opening Negligible Important only Very important Very important
characteristic
for interruption
High
of small
frequency
Negligible inductive Very important Very important
current
currents
intteruption
It can be seen that for low-frequency transients, such as those of transformer inrush
transients, high current and high frequency current interruption, current chopping and
re-strike characteristics of circuit breaker can be neglected in modelling opening
operation, and the prestrike can be neglected in modelling closing operation; the only
important feature that should be considered in detail is the mechanical pole spread, i.e.
closing time span (in general breaker poles do not close simultaneously, but with certain
time span). Therefore, in many previous studies of transformer inrush transients, such as
[6, 39, 70], each pole in a three-phase circuit breaker was modelled as an ideal time-
controlled switch: in opening operation, it opens at the first current zero crossing after
the ordered tripping instant (a current margin parameter can be included to approximate
current chopping); in closing operation, it behaves as an impedance changing
instantaneously from an infinite value to a zero value at the closing instant (the closing
instant can be at any part of a power cycle). Since three poles are represented by three
53
Chapter 2 Literature Review on Transformer Energisation Transients
separate ideal time-controlled switches, the closing time span between three poles can
be represented by closing time differences between the three time-controlled switches.
The circuit breaker closing time and closing time span between circuit breaker poles are
of stochastic nature; hence, several approaches have been proposed to construct
statistical switches with closing time and closing time span modelled by statistical
distributions. In [70], each pole of a circuit breaker was modelled by two contacts (one
is named as auxiliary contact and the other is named as main contact), as shown in
Figure 2-10.
Figure 2-10 Statistical switching model involving closing time span among three phases [70]
The closing time of the auxiliary contact Taux was used to represent the instant at which
the closing signal is ordered. It was considered to be the same for three phases and
follow a Uniform distribution with a typical range of one power frequency cycle (equal
to 20 ms for 50 Hz systems).
The closing times for three main contacts were defined as:
TAclose = Taux + τA ± TAr
TBclose = Taux + τB ± TBr (2.37)
TCclose = Taux + τC ± TCr
where TAr, TBr and TCr represent the closing offset time of three poles (each of them was
defined by a Gaussian distribution whose standard deviation was defined by one-sixth
of the maximum closing time span (MCTS)); τA, τB and τC represent the time delays
between the closing signal ordering and the actual closure of circuit breaker.
54
Chapter 2 Literature Review on Transformer Energisation Transients
Similar to modelling approach used in [70], the circuit breaker closing time span
modelled in [43] consists of four parameters:
• Common order time, torder. It was also characterized by a Uniform distribution
over a power frequency cycle.
• Random offset time for each pole (toffset,A, toffset,B and toffset,C). This offset time was
assumed to follow a Gaussian distribution, whose mean value is zero (assuming
that three poles tend to close simultaneously) and whose standard deviation is
MCTS/6.
The exact closing time for each pole was thus determined by:
TAclose = torder ± toffset,A
TBclose = torder ± toffset,B (2.38)
TCclose = torder ± toffset,C
It can be seen that in both modelling approaches, the MCTS determines the offset
closing time. However, MCTS is an uncertain value. According to [71], it is suggested
that the typical MCTS is between 3 and 5 ms. In [72], tests were carried out to study the
performance of a 110 kV circuit breakers (minimum oil circuit breaker and air-blast
circuit breaker) on energising transmission lines and it was shown that the MCTS is
normally smaller than 5 ms but could be as large as 10 ms. In [73], the performance of
400 kV and 220 kV circuit breakers for energising transmission lines in different
network topologies were experimentally investigated and it was shown that: for 400 kV
circuit breakers (minimum oil, air-blast or SF6 without switching resistance), the MCTS
was less than 9 ms; for 220 kV circuit breakers (minimum oil), the MCTS can
sometimes reach 16 ms.
In a large network, a proper boundary can be selected to reduce the network to a size
only covering the part of the network that is of interest for a specific study. This
boundary is normally set at the points where the system is very strong (i.e. large short
circuit level). These points can be the supply side of a substation transformer or a main
supply bus. The portion of network outside the boundary can be represented by a
55
Chapter 2 Literature Review on Transformer Energisation Transients
network equivalent. In many studies, such as [19, 36, 75, 76], the network equivalent is
modelled by an ideal voltage source together with a Thevenin equivalent impedance.
In the range of low frequency transients, loads are commonly modelled as a constant
impedance [77]. Naturally, this treatment also applies in the loads modelling for
studying transformer energisation transients; examples can be found in [2, 12, 78]. The
constant impedance model can be parallel-connected resistive and inductive elements
(loads vary with square of voltage magnitude) or can be of series-connected resistive
and inductive elements (loads vary with square of current magnitude). The power factor
of the load indicates the relative proportion of the resistive and inductive components in
the impedance.
56
Chapter 2 Literature Review on Transformer Energisation Transients
Figure 2-11 shows one commonly used generic circuit for studying sympathetic
interaction between two paralleled transformers. In the circuit, transformers TE1 and
TE2 are connected to an ideal voltage source through system resistance and
inductance ' ; 3 and ! are transformer winding resistances; '3u and '!u are
transformer leakage inductances; '3$ and '!$ are the magnetization inductances of
TE1 and TE2, respectively. # is equal to #$ ,5& -, where #$ is the amplitude of
source voltage and - is the energisation phase angle; # is the voltage of common
busbar; , 3 and ! are the currents flowing through the supply, TE1 and TE2,
respectively.
common S R2 L2σ
busbar
i2(t)
L2m
Rs Ls TE2
is(t)
R1 L1σ
Vs(t) Vc(t)
i1(t)
L1m
TE1
Figure 2-11 Generic circuit for studying sympathetic interaction between transformers
connected in parallel
Utilizing the coupled electromagnetic model proposed in [83], simulation was carried
out in [24] to study the sympathetic inrush interaction between the two parallel
connected transformers, in which case TE1 and TE2 are two identical 230/69 kV, 15
MVA single-phase transformers (referring to Figure 2-11). The simulated currents ! , 3
and are illustrated in Figure 2-12 (a), Figure 2-12 (b) and Figure 2-12 (c), respectively.
As can be seen, the inrush current ! reached maximum peak right after the energisation
of TE2 and then decayed gradually, while the sympathetic inrush current 3 built up in
TE1 gradually reached its maximum peak and then gradually decayed; the supply
current is the sum of the currents 3 and ! , showing the peaks of sympathetic inrush
current 3 and of the inrush current ! occur in direction opposite to each other, on
alternate half cycles.
57
Chapter 2 Literature Review on Transformer Energisation Transients
Inrush current i2
Supply current is
Y3
y ' 3 3 #$ ,5& -
w
Y3 Y! z
x
(2.39)
3 3 ! !
w
v 3 !
where Y3 , Y! are the flux-linkages of transformers TE1 and TE2, respectively, and
Y3 3 '3u '3$ and Y! ! '!u '!$ .
58
Chapter 2 Literature Review on Transformer Energisation Transients
' 1 : !:D
Y3 #$ sin& - p NY3 0 Y! 0O6 ; !;D
| }
Q 2
1 :
NY3 0 Y! 0O6
| }
(2.40)
;
2
' 1 : !:D
Y! #$ sin& - p NY3 0 Y! 0O6 ; !;D
| }
Q 2
1 :
NY3 0 Y! 0O6 ;
| }
(2.41)
2
where Q N 2 ' 2' ! O3/! , p X~*X5N&' 2' / 2 O ;
!
Y3 0 and Y! 0 are the initial flux of TE1 and the residual flux of TE2, respectively.
From equations 2.40 and 2.41, it can be seen that both Y3 and Y! consist of one
sinusoidal component and two exponential DC components. The AC component and the
first DC component are the same, but the second DC component in Y3 is opposite to that
in Y! , therefore 3 and ! are opposite to each other and appear alternately. Also,
because the DC components in Y! are negative, the maximum peak of ! would appear
right after the energisation of TE2, whilst the DC components in Y3 are of opposite
polarity and the time constant of the first DC component s3 N ' 2' / 2 O is
smaller than that of the second DC component s! N '/O, so 3 will gradually reach the
maximum peak, and gradually decay afterwards. The simplified analytical analysis
shows in a general way the variation of flux-linkages in TE1 and TE2 which depends on
the time constants formed by the inductances and resistances of the circuit branches. In
real situation, the core inductance is nonlinear and therefore the time constants cannot
be so readily determined.
In [22] and [24], the interactions between paralleled transformers were analysed using
the voltage drop across circuit resistances, with system and transformer winding
inductances neglected, which is summarized as follows (by referring to Figure 2-11).
59
Chapter 2 Literature Review on Transformer Energisation Transients
Before closing S, only the magnetizing current of the unloaded transformer TE1 flows
through the system; the source voltage # can be described by:
Y3
# 3 · 3
(2.42)
After closing S, saturation of transformer TE2 causes a transient inrush current ! which
flows through . Due to the unidirectional characteristic of the inrush current, each
cycle transformer T1 experiences an offset flux by an amount of:
Δ3 N 3 · 3 · ! O (2.45)
Meanwhile, an offset flux per cycle Δ! is produced in transformer TE2 by:
Δ! N ! · ! · 3 O (2.46)
At the initial stage, both Δ3 and Δ! are of the same polarity and mainly depend on
the voltage drop caused by the inrush current ! . The accumulation of Δ3 drives
transformer TE1 into saturation, while the effect of Δ! is to reduce the initial offset
flux in transformer TE2 so as to produce the decay of inrush current ! .
As the transformer TE1 becomes more and more saturated, a sympathetic inrush current
3 gradually increases from the steady state magnetizing current to a considerable
magnitude. Noted that as the transformer TE1 saturates with the polarity opposite to that
of transformer TE2, the peaks of the sympathetic inrush current 3 are with polarity
opposite to that of inrush current ! , on alternate half cycles. As a result, the voltage
asymmetry on transformer terminals caused by the inrush current ! during one half
cycle is reduced by the voltage drop produced by the sympathetic inrush current 3
during the subsequent half cycle. This decreases both Δ3 and Δ! , and therefore
60
Chapter 2 Literature Review on Transformer Energisation Transients
reduces the changing rate of the magnitude of both the increasing sympathetic inrush
current 3 and the decaying inrush current ! .
After a certain time, the increase of 3 and decay of ! can reach a point that:
61
Chapter 2 Literature Review on Transformer Energisation Transients
Sympathetic inrush has been encountered in many practical systems and caused
significant concerns. In [84], sympathetic interaction between transformers in a 20 kV
converter test facility was reported. The configuration of the test facility is shown in
Figure 2-13 (a). During converter testing, the active power is circulating between S1-
T1-TO-T2-S2-S1 and only the losses are compensated from a 20 kV grid (with 160
MVA short-circuit level). The two transformers T1 and T2 need to be energised on a
daily basis for carrying out tests. To reduce inrush current, a 100 Ohm short-circuit
current limiting resistor was connected, which is large enough to limit the inrush current
magnitude to values below 150 A and to damp out inrush current in less than 50 ms.
However, sustaining sympathetic inrush currents were encountered. As shown in Figure
2-13 (b), T1 was energised at 52.2 s and then T2 was energised at 52.26 s; although
inrush current caused by energising T1 was damped out in one cycle, long-duration
sympathetic inrush currents were induced after energising T2. The damping resistor,
62
Chapter 2 Literature Review on Transformer Energisation Transients
which helped reduce inrush current in the case of energising T1, caused voltage
asymmetry which resulted in sympathetic inrush when T2 was energised.
T1 Energised
T2 Energised
Voltage dips caused by sympathetic inrush between 100 MVA 220/23 kV transformers
were reported in [85]. The configuration of the electrical system subjected to
sympathetic inrush is shown in Figure 2-14. The substation is fed by two 220 kV
overhead lines with length of 178 km; it originally consisted of two 100 MVA, 220/23
kV transformers (T1 and T2) connected in parallel to supply power to mining facilities.
A new transformer T3 was added to meet increasing demands. When energising T3,
sympathetic inrush was induced in the two already connected transformers T1 and T2.
The energisation also resulted in high distortion of voltages and caused tripping of those
63
Chapter 2 Literature Review on Transformer Energisation Transients
(a) Measurement 1
(b) Measurement 2
Figure 2-15 Measured voltage dips at 23 kV busbar [85]
64
Chapter 2 Literature Review on Transformer Energisation Transients
In [8], how mechanical forces build up under inrush current was investigated, compared
to those occurring under short-circuit. The investigation was based on 2D and 3D
modelling of a 268 MVA, 525/17.75 kV three-legged step-up transformer. It was found
that inrush currents with peaks of more than 70% of the rated short-circuit current
magnitude would induce forces higher than those under short-circuit condition; these
forces summed up on the high voltage winding (normally the energised winding) can be
three times higher.
In [9], the radial and axial electromagnetic forces due to inrush currents were examined
for a three-phase, three-legged 66/11 kV, 40 MVA power transformer. The study shows
that the axial forces due to inrush current are always larger than those caused by short
circuit current and the radial force applied on High-Voltage (HV) winding is about three
times the corresponding force under short-circuit condition.
Even though inrush currents are normally smaller than short-circuit current, they are
with a much longer duration. In addition, the duration of inrush current can be further
prolonged under sympathetic inrush. This may cause winding damage or insulation
failures a certain time span after transformer energisation.
In [78], simulations were performed using EMTP to analyse the inrush transients
resulted from energising a 27 MVA transformer from an 8.3 MVA diesel generator
65
Chapter 2 Literature Review on Transformer Energisation Transients
under different operating conditions. The developed network model included generator,
generator control, transformer and loads; the governor was not included in the generator
control, because the governor time constant is longer than the analysed inrush transients
and the level of active power consumption during the transformer energisation is low.
The simulation results show that: although the inrush currents are of magnitudes lower
than three-phase short-circuit currents, they are much higher than normal operating
current and may generate high electromagnetic torch oscillations; these oscillation may
subject the shaft to high torsional stresses which could lead to fatigue failure if the
transformer energisation from the diesel generator is frequent.
Figure 2-16 Simplified single-line diagram of wind farm collection grid during an emergent
islanded condition [2]
As can be seen, it consists of eight cable feeders; each feeder contains five wind
turbines and each wind turbine connects a circuit breaker, a wind turbine transformer (4
MVA, 33/0.69 kV, Dyn) and a low-voltage auxiliary load (18 kVA); the emergency
diesel generator (1.6 MVA, 33/0.4 kV) is located at the offshore platform together with
two shunt reactors (1.5 MVar) which are to balance the reactive power generated by the
33 kV cables.
66
Chapter 2 Literature Review on Transformer Energisation Transients
The collection grid was modelled in PSCAD/EMTDC: the diesel generator was
modelled in a way similar to that used in [78]; the 33 kV submarine cable sections were
represented by nominal pi-sections (the length of the cable section between two adjacent
wind turbines is slightly above 1 km); the loads connected on the low-voltage side of
each wind turbine transformer were modelled by constant impedances; the wind turbine
transformers were modelled by the PSCAD classical model in which each phase of the
transformer is represented by a separate single-phase transformer model with no
coupling between phases; additional dc-current sources were connected to the wind
turbine transformer low-voltage side to simulate residual flux in the transformer.
In the paper, there were in total five case studies which are summarized in Table 2-4. As
can be seen, the case studies mainly considered the effects of residual flux, sympathetic
interaction between wind turbine transformers and the response of Automatic Voltage
Regulator (AVR). In all the cases, the shunt reactors and 33 kV cables were connected;
only the wind turbine transformer located farthest from the platform was energised and
the energisation instant was at the positive-going zero crossing of phase-to-ground
voltage. It was found that: the sympathetic inrush current induced in the already
connected wind turbine transformers imposes further reactive power demand on the
diesel generator; increase in the speed of AVR response from medium to high can result
in larger sympathetic inrush currents in the already connected wind turbine transformers
and hence higher reactive power demand from the diesel generator.
Table 2-4 Summary of case studies carried out in [2]
Energised wind turbine Adjacent wind turbine Residual AVR
Case
transformers transformers flux response
1 Not connected Zero
2 Not connected
Medium
The wind turbine Connected
3 speed
transformer farthest from (With saturation) 70%
the offshore platform Connected
4 (No saturation)
Connected High
5 (With saturation) speed
67
Chapter 2 Literature Review on Transformer Energisation Transients
In most previous contributions, the harmonic analysis of inrush current was performed
by looking at the variation of its harmonic content with time. The first contribution
showing the harmonic content variation was presented in [10]. In the paper, the
magnitude and phase shift of each harmonic component were obtained from a Fourier
analysis for each cycle of the inrush separately. This approach was also followed by
other contributions in [14, 24, 87, 88]. Typical harmonic analysis results given by [24]
are shown in Figure 2-17.
(a) Harmonic contents of the inrush current (b) Harmonic contents of combined inrush
drawn by the transformer being energised current and sympathetic inrush current
Figure 2-17 Variation of harmonic content of inrush current as a function of time [24]
Figure 2-17 (a) illustrates the harmonic components of inrush current alone (without
sympathetic inrush current involved). It can be seen that: the amplitude of any harmonic
component during one cycle is generally different from its amplitude during another
cycle; the second order harmonic is the dominant one; the higher the magnitude of the
inrush current at any one cycle, the higher the second order harmonic content of that
cycle; the higher the harmonic order, the smaller the magnitude of the corresponding
current component in the inrush current; for some harmonics, their highest amplitudes
do not appear at the first cycle after transformer energisation, such as the third and
fourth order harmonics; some harmonic components change their phase from negative
to positive, or vise versa, after their amplitude pass zero, such as the fourth and fifth
68
Chapter 2 Literature Review on Transformer Energisation Transients
order harmonics. Figure 2-17 (b) shows the harmonic components of the current
combining inrush current and sympathetic inrush current, in which case, the even order
harmonic components decay rather quickly, whereas the odd harmonics increase and
continue to stay for a considerable period of time.
During transformer energisation, the inrush current typically occurs in only one winding
of the transformer and thereby produce a differential current that may result in the
operation of transformer differential protection [89]. Since transformer inrush is not a
fault event, the differential protection must be restrained for this condition.
Experience shows that for most transformer application, the setting can effectively
differentiate the inrush events and internal fault events via harmonic restraint. However,
it should be aware of the fact that modern transformers may be characterized by lower
second harmonic ratios because of higher designed flux density and the use of step-lap
type joint [90]; in addition, in the case of transformer ultra-saturation, the percentage of
second harmonic can fall below 5%, inevitably leading to mal-operation [17, 91];
furthermore, in the case of current transformer (CT) saturation during internal faults, the
fault current transformed to the secondary may contain amounts of second harmonic
higher than the setting and thus cause incorrect restraining [16]. In view of the
limitations of the second harmonic restraint function, there are other methods proposed,
such as a complex second harmonic restraint [92], flux restraint [93], or the use of
artificial neutral network [94].
69
Chapter 2 Literature Review on Transformer Energisation Transients
with magnitude over twice the rated voltage for as long as 100 or more cycles, imposing
large risk of burning insulators, arresters and damaging transformer insulation. This
type of overvoltage has been identified in following cases:
• Energising the convertor transformers in HVDC substations consisting of ac
filter circuit [13, 95];
• Restoration of a bulk power supply system [12, 15, 76, 96];
• Energising transformer in systems with long length cables [14, 36];
• Energising transformer in some industrial distribution systems with installation
of power factor correction capacitors [97] or pulse-type loads [98].
In these cases, the systems consist of the following common characteristics: pronounced
parallel resonance points (such as in the systems with long transmission lines or reactive
components like filters and capacitor banks) and low degree of damping (the system is
light-loaded or non-loaded). Selected examples corresponding to some of the typical
cases are given below.
Normally, the HVDC station is directly fed by generators without local ac loads being
connected, i.e., low damping. The ac filter circuit connected at the HVDC stations can
form several parallel resonance points in the impedance-frequency characteristic of the
system. The inrush currents resulted from energising the converter transformer can
repeatedly shock the ac system – ac filter combination once per cycle, and due to the
slow decay of the inrush currents, result in overvoltages lasting many cycles, as shown
in Figure 2-18.
Figure 2-18 Field measured overvoltages caused by transformer energisation in HVDC stations
[13, 99]
The value of the overvoltages depends on the value of the harmonic current at which the
resonance occurs; the largest value of a certain harmonic might occur a long time after
energisation. Measurements show that the peak value of such overvoltages can be 1.7
pu [99], while simulations show that it can be over 2 pu [13, 95].
70
Chapter 2 Literature Review on Transformer Energisation Transients
It was also shown that severer overvoltages can be excited by energising the inter-
connector link (i.e., energising the step-up transformer, the sub-sea cable and the step-
down transformer simultaneously). The authors suggested that the overvoltage problem
can become less severe with increased levels of generation and load on the system.
71
Chapter 2 Literature Review on Transformer Energisation Transients
Harmonic resonant overvoltages are also likely to occur during system restoration,
because the networks, after a complete or partial collapse, are lightly loaded and with
low system resonant frequencies. An example configuration for a system restoration
was given in [96] and is shown here in Figure 2-20 (a). The basic procedure to restore
such a system would be: start up the generators; connect the local load at busbar B4;
energise the transmission line together with the shunt reactor connected to busbar B3;
finally, energise the unloaded transformer. The energisation of the unloaded transformer
through long transmission line resulted in significant resonant overvoltages at busbar B3,
as shown in Figure 2-20 (b). For harmonic resonant overvoltages during system
restoration, they can be controlled by several methods [15]: increase resistive loading,
bring additional generators on line or decrease the magnitude of generator terminal
voltage.
From the harmonic analysis of the system current combining inrush and sympathetic
inrush currents (see Figure 2-17 (b)), it can be seen that the sympathetic interaction can
probably reduce the severity of harmonic overvoltages for systems resonating at even
ordered harmonic frequencies; however, for systems resonating at odd ordered
harmonic frequencies, harmonic overvoltages are likely to be prolonged [24].
Figure 2-20 System configuration at the beginning of a restoration procedure and overvoltage
resulted from energising a transformer [96]
72
Chapter 2 Literature Review on Transformer Energisation Transients
drive, when the voltage drops below a critical level for a long duration, the drive might
function abnormally or even shut down.
Transformer energisation is a planned operation and may cause severe voltage dips. It
gains increasing attention in recent years due to the need to comply with tightened grid
code requirements. In UK, a 3% threshold is normally applied to the voltage dips
caused by transformer energisation. This threshold is derived from the Engineering
Recommendation P28 (ER-P28) which defines the curve describing tolerable dip
magnitude (in other words, the size of voltage change) against the interval between each
voltage change, as shown in Figure 2-21 [4].
The curve shows the tolerable size of voltage change increases with the time between
each change. For examples, if the time between each change is 1 second, the allowable
size of voltage change is 0.4%; if the time between each change is 200 seconds, the
limit will be 2%; and when the time between each change is equal to or more than 750
seconds, the maximum allowable size of voltage change is 3%.
Figure 2-21 Size of voltage change against the time between each change [4]
The voltage dips caused by energising large GSU transformers from HV transmission
grid were investigated in [19], [51] and [101].
In [19], the system under study comprises generating plants, long transmission lines and
power quality sensitive loads, which can be referred to the simplified single-line
diagram shown in Figure 2-22.
73
Chapter 2 Literature Review on Transformer Energisation Transients
3.7 km
89 km
An EMTP model of the 138 kV network was developed to carry out assessment of
voltage dips caused by energisation of the GSU transformer. In the network model:
• Positive and zero sequence Thevenin impedances were used to represent the
remaining network and the 500 kV connection to the main grid;
• Generating plant located at John Hart was represented by an ideal voltage source
connected to equivalent sub-transient impedance;
• Transmission lines were modelled by Bergeron model;
• GSU transformer was modelled by three single-phase two-winding STC
transformers (the three single-phase transformers are connected in grounded-star
on the 138 kV side, and in delta on the 21 kV side); type-96 nonlinear inductors
74
Chapter 2 Literature Review on Transformer Energisation Transients
Simulation assessment was carried out to estimate the voltage dips under the worst
energisation condition which was assumed in the study as:
• All three phases simultaneously switched at the zero crossing of phase A voltage;
• Maximum residual flux of negative polarity in phase A and the other two phases
with half of the maximum residual flux of positive polarity (the maximum
residual flux was assumed to be the flux retained at the instant when the
magnetizing current become zero following the core hysteresis curve).
In the simulation of the worst case energisation, it was estimated that the maximum dip
magnitude of the RMS voltage dips observed at the mill was about 0.27 pu.
The study carried out in [51] not only assessed voltage dips under the worst energisation
condition, but also investigated the influence of network configuration variation on
voltage dips. Figure 2-23 shows the network studied in [51].
A power station is located between a bulk supply point and a 132/33 kV transmission
substation. The bulk supply point contains three links to higher voltage grid through
three 330/132 kV step-down transformers. Network loads are connected to the bulk
supply point as well as the transmission substation. There are five feeders going into the
power station, two from the transmission substation and three from the bulk supply
point. The GSU transformer connected to the 400 MW generator needs to be energised
from the 132 kV grid. The nameplate data of the GSU transformer are 500 MVA,
75
Chapter 2 Literature Review on Transformer Energisation Transients
132/21 kV, 18% impedance, YNd11 winding connection with 132 kV star side being
solidly grounded and its tapping range is from +15% to -5% (each step is of 1.25%).
The network was modelled in PSCAD/EMTDC, consisting of the supply network and
the GSU transformer (the generator was not included in the model). The GSU
transformer was represented by three separate single-phase transformers. The effect of
residual flux in the GSU transformer was modelled by using adjustable DC current
sources in parallel with the transformer HV terminals.
In the assessment, the assumed worst energisation condition was that: energising at zero
crossing of one phase voltage and the maximum residual flux was assumed to be 90%
of the peak nominal flux. Based on the same energisation condition, voltage dips
resulted from energising GSU transformer under three different network configurations
were assessed. The three considered network configurations (C1, C2 and C3) are:
• C1, as shown in Figure 2-24 (a), all busbars are made solid; only two supply
sources are connected at the bulk supply point; all lines to the substation of the
generating plant are switched in;
• C2, as shown in Figure 2-24 (b), all busbars are made solid; all three supply
sources are connected at the bulk supply point; the generation plant is only
supplied by one single feeder from the bulk supply point;
• C3, as shown in Figure 2-24 (c), the 132 kV busbar at the bulk supply point are
split; all lines linking the generator substation and the transmission substation
are disconnected.
Under the three network configurations and the same worst case energisation condition,
the estimated inrush current (peak, duration) and voltage dip magnitudes observed at the
terminal of the GSU transformer and at the transmission substation 132 kV busbar are
summarized in Table 2-5. It was found that, by changing configuration from C1 to C3,
the RMS dip magnitude observed at the transmission substation 132 kV load bus can be
reduced from 16.8% to 6.3%. This is because the point of common coupling connecting
all loads and the energised GSU transformer has been effectively moved to 330 kV grid.
This increased the electrical distance between the TS and the energised transformer, so
the voltages at transmission substation 132 kV busbars become less sensitive to the
energisation of GSU transformer.
76
Chapter 2 Literature Review on Transformer Energisation Transients
(c) Configuration C3
Figure 2-24 Network configurations under comparison
Table 2-5 Estimated inrush current peaks, duration and voltage dip magnitudes resulted from
the worst case energisation under different network configurations [51]
77
Chapter 2 Literature Review on Transformer Energisation Transients
Figure 2-25 Single line diagram of Jeju power system in Korea [101]
The network of Jeju power system in Korea was modelled by following the approaches
below:
• Generator was modelled by using the SM model in EMTP;
• Transmission lines were represented by lumped parameter model;
• Main GSU transformers were modelled based on transformer test report data and
the core saturation curve was modelled to include hysteresis characteristic;
• Substation transformers were modelled but without considering hysteresis;
• Wind turbine and HVDC system were not considered in the network model.
Energisation of the main GSU transformer was simulated under the conditions shown in
Table 2-6. Voltage dip magnitudes at all the substations were observed, as shown in
Figure 2-26. It was found that the range of the voltage dip magnitudes is between 0.15
pu and 0.2 pu. The substations with closer distance to the energised transformer are with
relatively larger voltage dips. Regarding the effect of the parameter on voltage dip
magnitude, the assessment showed that voltage dip magnitude is sensitive to the circuit
breaker closing angle, core residual flux and load variation.
78
Chapter 2 Literature Review on Transformer Energisation Transients
Figure 2-26 Voltage dip magnitudes resulted from different energisation angles when residual
flux is 28.3% and system loading is at its peak [101]
In recent years, wind farms are becoming more prevalent. A large number of wind
farms have been connected to distribution or transmission grids and more are to be
designed and commissioned. A common requirement for the connection of wind farm in
the UK is that the energisation of transformers should not result in voltage dips
exceeding 3%, according to the P28 requirements. A number of contributions have been
devoted to transformer inrush studies for wind farm grid connections [50, 75, 102, 103],
all of which were performed using PSCAD/EMTDC.
In [102], simulations were conducted to identify the possible voltage dips resulted from
energising wind turbine transformers in two wind farm sites: one consists of 20 wind
turbines, each rated at 1.5 MVA; the other consists of 17 wind turbines, each rated at
1.0 MVA; both of them are directly connected to 33 kV distribution networks. The
simulation studies demonstrated that the voltage dips can be affected by circuit breaker
closing time, closing time scatter among three poles, the number of wind turbine
transformers to be energised simultaneously and the number of already connected
transformers. It was found that, to meet grid code requirements, up to two wind turbine
transformers can be energised simultaneously
Similar simulation studies were performed in [75]. Two wind farm sites (named here as
A and B) were investigated, as shown in Figure 2-27. Wind farm A contains 15 wind
turbine transformers (33/0.69 kV, 1.5 MVA, 6%) which are connected via two 33 kV
feeders to the main wind farm switchboard and the point-of-common-coupling (PCC) is
at 33 kV busbar; wind farm B contains 52 wind turbine transformers (33/0.69 kV, 2.6
MVA, 8.28%) which are connected via eight radial feeders back to the wind farm 33 kV
79
Chapter 2 Literature Review on Transformer Energisation Transients
Sympathetic inrush in wind farm B was also assessed. The scenario considered was
simultaneous energisation of nine 2.6 MVA wind turbine transformers under the worst
case energisation condition; the wind turbine transformers on other feeders were
assumed already connected, and, for each feeder with N wind turbine transformers, a
single equivalent transformer with N×2.6 MVA was used. It was found that the
sympathetic inrush current at one adjacent feeder can increase to a peak of 94 A about
12 cycles after the energisation.
80
Chapter 2 Literature Review on Transformer Energisation Transients
In [50], some curves were identified to correlate the system fault level with voltage dips
resulted from energising some typical wind turbine transformers. The main contribution
of these curves is to allow the P28 assessment to be preliminarily carried out for a
proposed wind farm connection before the detailed design information is available.
Based on the four contributions specifically targeted to assess voltage dips caused by
energisation of wind turbine transformers for wind farm connections, Table 2-7
summarized the influential parameters that have been considered. In addition, the
parameters studied by each paper are illustrated in Table 2-7 as well. As can be seen,
none of them could take into account all the possible influential parameters. The most
frequently concerned parameters are switching angle, the number of transformer being
energised simultaneously (aggregated energisation) and sympathetic inrush.
81
Chapter 2 Literature Review on Transformer Energisation Transients
References
Parameters
[102] [75] [50] [103]
Switching angle + + + +
Closing time span + - - -
Residual flux - - + +
System strength - - - +
Core saturation inductance + - - +
Aggregated energisation + + + +
Sympathetic inrush + + + -
Energisation sequence - + - -
Voltage dip propagation + - - -
+: parameter taken into account; -: parameter not considered
In [104], Monte-Carlo method was used to study the frequency of occurrence of inrush
current first peak when energising a single-phase transformer at stochastic circuit
breaker closing time and with stochastic transformer core residual flux. In the study, it
was assumed that: closing time varies in a range of one power frequency cycle; residual
flux varies in a range of ±1 pu of peak nominal flux; both ranges were characterized by
Uniform distribution. Three cases were simulated, see Table 2-8.
Table 2-8 Cases studies of the influences of stochastic variables on the inrush current of a
single-phase transformer
Circuit breaker closing time Transformer residual flux
Case
Range Distribution Range Distribution
1 One cycle (20 ms) Uniform Zero
Fixed at positive-going zero crossing of
2 ±1 pu of the peak
the applied voltage Uniform
nominal flux
3 One cycle (20 ms) Uniform
82
Chapter 2 Literature Review on Transformer Energisation Transients
The result of the Case 3 is shown in Figure 2-29 which suggests, under simultaneous
variation of closing time and residual flux, the inrush current first peaks obtained from
stochastic simulation runs seem to follow an asymmetrical exponential distribution.
Figure 2-29 Frequency of inrush current first peaks when residual flux and closing time vary
stochastically [104]
83
Chapter 2 Literature Review on Transformer Energisation Transients
In the system studied in [38], the source plant consists of three generators (total capacity
is 900 MW); the length of overhead line is 50 km and the auxiliary transformer is of
three-phase, shell-type, 96 MVA and 400/6.8 kV.
To represent the random closing times, it was assumed that: the first pole of the circuit
breaker is closed anytime within one cycle (20 ms for 50 Hz system); relative to the first
pole being closed, the second and third poles are closed within a standard deviation (one
cycle). Regarding the random residual flux, it was assumed that: the residual flux values
are in a range of ±0.8 pu of peak nominal flux; the sum of the three-phase residual flux
is recognized as zero. The possible range of the sub-transient reactance variation was
assumed to be ±15% of the nominal value; the range for the phase-to-earth line
capacitance was assumed to be ±5% of the base value. In the studies, 25 deterministic
combinations of sub-transient reactance and line capacitance were formed; for each
combination, 100 runs were carried out with stochastic residual flux and closing times.
Based on the same type of network, a similar study was performed in [43]. Compared to
the work shown in [38], the main differences made in [43] include:
• Adding transformer air-core reactance as a stochastic parameter and it was
described by a Uniform distribution covering ±20% of a base value;
• The stochastic variation of sub-transient reactance and phase-to-earth line
capacitances were described by a Uniform distribution covering their
corresponding ranges;
• Circuit breaker closing time were represented using the approach shown in
section 2.2.3;
• The residual flux was modelled by two parameters: the maximum residual flux
amplitude ( and a parameter identifying the flux pattern; twelve equally
probable residual flux patterns were considered, which includes:
, /2, /2
, , , 0
and , /2, /2
for all the possible
phase permutations.
84
Chapter 2 Literature Review on Transformer Energisation Transients
The oldest mitigation strategy to reduce inrush current is the use of the pre-insertion
resistor as an intermediate step in closing a switch [1, 105]. When a transformer draws a
large inrush current, the voltage drop across the pre-insertion resistor helps reduce the
voltage applied across the transformer and hence the flux in the transformer, quickly
reducing the inrush current. The effectiveness of this method is largely influenced by
the time of inserting the resistance and that of by-passing it. Although this method is
simple, it requires upgrading existing substation breakers, which is expensive; in
addition, if a pre-insertion resistor is not properly by-passed after the switching process,
the breaker could be damaged.
Reducing residual flux in the transformer core could be a possible strategy for reducing
inrush currents. In [106], it was shown that residual flux can be significantly reduced by
using a low, variable frequency voltage source to demagnetize the transformer. This
method has only been shown in simulation but has not yet been realized experimentally.
Modifying transformer design may help reduce inrush current. The modification
suggested by [107] is to change the distribution of the coil winding in a way to increase
the transient inductance of the primary coil. The one presented in [108] employs an
auxiliary winding to create an air gap in the transformer core so as to control the core
reluctance. These methods have only been tested in laboratory on small single-phase
85
Chapter 2 Literature Review on Transformer Energisation Transients
transformers. It is however challenging and costly to modify the design of large power
transformers.
Techniques involving the control of circuit breaker closing time are probably the most
widely discussed for mitigating inrush current. The most effective technique is to
optimize the circuit breaker closing time according to transformer core residual flux
[109]. The basic principle is that: measure the residual flux in the core during
transformer de-energisation; with the residual flux known, calculate the optimum
voltage angle resulting in the measured residual flux equal to the prospective flux; the
inrush current could be eliminated if the energisation is carried out at the optimum
voltage angle. Applying this technique requires installation of additional measurement
equipment on the transformer terminals, which is costly. In addition, its performance is
very sensitive to the accuracy of residual flux measurement and circuit breaker
operation.
The influence of tap position on the magnitudes of inrush current has been addressed in
[110]. It was considered that adjusting tap changer position to the lowest possible
energisation voltage could result in smaller inrush current. So far, this strategy is the
most commonly applied. The effectiveness of this technique has been proven in the
cases of energising generator step-up transformers in high voltage networks [51, 85].
2.5 Summary
This chapter reviewed the published work studying transformer energisation transients,
regarding the following aspects:
• Approaches for calculating transformer inrush current;
• Modelling system components in EMTP for transformer energisation studies;
• Investigation case studies of transformer energisation transients;
• Possible measures for mitigating transformer inrush current.
86
Chapter 2 Literature Review on Transformer Energisation Transients
87
Chapter 2 Literature Review on Transformer Energisation Transients
88
Chapter 3 Field Measurements, Network Model Development and Validation
89
Chapter 3 Field Measurements, Network Model Development and Validation
Figure 3-1 South West Peninsula system as part of National Grid’s transmission system in
England & Wales
In 2009, a new Combined Cycle Gas Turbine (CCGT) power plant, with total installed
capacity of 905 MW, was commissioned and connected to the Lagage substation, which
is to underpin the National Grid's reinforcement to the South West Peninsula system. In
the future, it is estimated that the local demand will modestly increase, yet significant
renewable generations are to be installed, driven by the commitments to connect more
low carbon generations to the transmission system. As shown in Table 3-1, substantial
amounts of generations have been approved and they will be connected to the South
West Peninsula in the near future.
Table 3-1 Plan of new generation installations at the South West Peninsula system [111]
90
Chapter 3 Field Measurements, Network Model Development and Validation
Figure 3-2 Schematic diagram of South West Peninsula system under detailed studies
91
Chapter 3 Field Measurements, Network Model Development and Validation
At the commissioning stage of the power plant, two attempts were made to
simultaneously energise GSU transformers T2 and T3. Prior to the first attempt:
• the synchronous compensator located at substation H, whose capacity is 140
MW, was not connected;
• the coupling circuit breaker CB1 at substation K was closed;
• all circuits in the 400 kV transmission grid were in service;
• GSU transformer T1 was already connected to the grid through power feeder 2,
but the generator unit was not synchronized.
When the first attempt was conducted by closing the CB2, voltage dips were observed
in the South West system, as shown in Figure 3-3 (voltages shown are line-to-line
voltage).
418 Pre-Energisation
414 Post-Energisation
Sub-A
418
Sub-G
414
414
Sub-F 414
Sub-B
418
413 416 414 416 413 413
399 401 388 Sub-J 410 413
Sub-D
Sub-H Sub-I Sub-K Sub-E Sub-C
Figure 3-3 South West Peninsula voltage depression resulted from the first attempt
Voltage dips in those substations in proximity to the power plant were relatively large;
the maximum voltage dip was observed at substation K, which was about 7.8%. As
consequences of the dips, there were responses from reactive compensation around the
South West system; low voltage alarms in the Integrated Energy Management System
were triggered and downstream distribution utilities were affected.
92
Chapter 3 Field Measurements, Network Model Development and Validation
Several days after the first attempt, the CB2 was opened for carrying out further
commissioning,, which led to the second attempt of energising T2 and T3. The initial
system conditions prior to the second attempt were similar to those in the first attempt
except that more proactive measures
measures were taken to reinforce reactive power support.
support
Dynamic System Monitoring equipment was set up at substation K to sensitively
sensitive
capture voltage and current variations.
variations During the second attempt, three-phase
t inrush
currents were captured by the current transformers located at power feeder 1, which are
shown in Figure 3-4.. As can be seen, the maximum peaks of phase A and phase B are
almost the same, both of which are about half that of phase C. In addition,
addition three-phase
RMS (root mean square) voltage dips were also measured at substation K, which are
shown in Figure 3-55. It can be seen that: the resulted three-phase
phase voltage dips were
unbalanced; the biggest voltage dip appeared at phase C, which was about 7%.
Figure 3-5 RMS voltage dips measured at substation K in the second attempt
93
Chapter 3 Field Measurements, Network Model Development and Validation
During the time of carrying out field measurements, there were four energisation cases,
which are illustrated in Table 3-2. The first two energisation cases (Case E1 and Case
E2) were trial energisation made by the power plant engineer by quickly pressing and
releasing the circuit breaker closing button. The trial energisation might be to reduce the
transformer residual flux. In the energisation Case E3 and Case E4, the circuit breaker
closing was completely conducted and the corresponding GSU transformers were
successfully energised.
For all the energisation cases, currents and voltages were simultaneously measured on a
couple of locations including:
• Three-phase currents and voltages at power feeder 1;
• Three-phase voltages at substation I;
• Three-phase currents at the circuit I-K linking substation I and K (the circuit
position is highlighted in Figure 3-2).
Measurement was not carried out at power feeder 2 because the monitor on that feeder
was not working at the time.
From the current waveforms shown in Figure 3-6, sharp pulses can be seen at the instant
94
Chapter 3 Field Measurements, Network Model Development and Validation
of energisation t0, which could be due to the charging of stray capacitance. The interval
between the energisation instant t0 and the inrush starting instant t1 is about 4 ms.
Because the transformer is of star-delta-connection, the inrush currents might be
qualitatively analysed by the helping effect theory explained in [112]. At t1, the inrush
started in phase B first; from t1 to t2, the inrush in phase B was aided by phases A and C
through the medium of the delta winding. At t2, phase A started to experience its own
inrush, as the current began to increase rapidly. This imposed demands on phase C from
both phase B and phase A. The two demands were in opposite direction and the
helping-current from phase C would flow in the direction of the maximum
instantaneous demand. From t2 to t3, the phase B current was larger than phase A, so the
phase C current remained in positive polarity but decreasing. At t3, the demand from
phase A and phase B were equal in magnitude but opposite in polarity, therefore the
helping-current from phase C became zero. After t3, the current magnitude of phase A
was larger than that of phase B, the helping-current in phase C therefore increased in
negative polarity. At t4, the inrush in phase B ceased. From t4 to t5, phase B and phase C
were helping phase A. This whole process could be repeated in the next cycle; however,
due to the circuit breaker opening at t5, the three-phase currents through power feeder 1
returned to zero.
From the circuit I-K currents shown in Figure 3-7, it can be seen that: before the
energisation, the feeder currents only contained steady state load currents flowing
95
Chapter 3 Field Measurements, Network Model Development and Validation
through circuit I-K (the peak magnitudes for phase A, phase B and phase C were 344 A,
346 A and 334 A, respectively); after the energisation, the load currents were
superimposed by a portion of inrush currents during the inrush period, which increased
current magnitudes of phase B and C and decreased the current magnitude of phase A.
After the inrush period, only the load currents flowed through the circuit I-K.
Figure 3-8 shows the voltage variation observed at power feeder 1. As can be seen,
three phases were almost simultaneously energised at the positive-going zero crossing
of phase B line-to-ground voltage. The build-up of voltages exhibited distortion due to
the influences of the inrush currents. The sharp dip at the voltage trace of phase C was
most probably due to the re-ignition when the circuit breaker pole connecting phase C
was trying to open at the current zero crossing appeared at the instant t3 (as shown in
Figure 3-6). CB2 was opened at the instant when all three-phase inrush currents
decreased to zero (i.e., at the instant t5 as shown in Figure 3-6). This instant was about
16.9 ms after the energisation instant. The opening of CB2 de-energised the GSU
transformers and initiated a ring-down process which involves exchange of energy
between electrical field in capacitances and magnetic field in the core of the GSU
transformers.
From Figure 3-9, it can be seen that: prior to the energisation, the peak magnitudes of
three-phase line-to-ground voltages at substation I for phase A, B and C were 333.8 kV,
333.5 kV and 331.0 kV, respectively. After the energisation, three-phase voltage
magnitudes were all reduced by the inrush currents during the inrush period; the most
affected phase was phase B whose peak magnitude was reduced by 36.1 kV.
96
Chapter 3 Field Measurements, Network Model Development and Validation
The measured inrush currents are shown in Figure 3-10. The closing time difference
among three phases can be interpreted by observing the zoom-in waveforms of the pulse
currents: the circuit breaker pole of phase C was the first one being closed, which was
followed by the closing of phase A and then the closing of phase B; between the closing
of phase C and phase A, the time interval was 0.3 ms; between the closing of phase C
and phase B, the time interval was 1.5 ms. The time interval between the energisation
instant of phase C and the inrush starting instant was about 4.5 ms, which is 12.5%
longer than that observed in Case E1. It can be also seen that the inrush started from
phase B first and, from the inrush starting instant to its peak, it was accompanied by the
helping-effect currents from phase A and phase C; after its peak, inrush also appeared in
phase A. The total inrush time only lasted for about 8.2 ms which is less than half cycle.
The inrush current peak magnitudes were significantly lower than those observed in
Case E1. As can be seen in Figure 3-11, during the inrush period, the small inrush
current being superimposed on the load currents flowing through circuit I-K, only
caused a slight distortion.
97
Chapter 3 Field Measurements, Network Model Development and Validation
The voltage waveforms measured at power feeder 1 are shown in Figure 3-12. It can be
seen that: phase A was energised almost at the negative-going zero crossing; however, it
was not the phase which experienced the biggest inrush peak. In fact, the peak
magnitude of phase B inrush current was the biggest, according to Figure 3-10. This
indicates that the forming of inrush might be influenced by residual flux. Once the
inrush currents decreased to zero, the transformers were de-energised. The de-
energisation time was about 13 ms after the energisation instant. Following the de-
energisation, the ring-down transient began. Due to the small inrush current magnitudes,
the voltages observed at substation I after the energisation were only slightly affected,
as shown in Figure 3-13.
98
Chapter 3 Field Measurements, Network Model Development and Validation
In Figure 3-15, the inrush currents with a longer time range is illustrated. As can be
seen, about seven cycles after the energisation, the inrush current waveforms of phase B
and phase C abruptly jumped. This abrupt change of inrush current waveform could be
due to the inaccurate measurement caused by the saturation of the current transformer
(CT). It will be further addressed in the section regarding CT modelling and simulation.
99
Chapter 3 Field Measurements, Network Model Development and Validation
Figure 3-14 Three-phase currents measured at power feeder 1 in Case E3 (initial cycles)
Figure 3-15 Three-phase currents measured at power feeder 1 in Case E3 (long duration)
Before and after energisation, three-phase voltages observed on power feeder 1 are
shown in Figure 3-17. The closing time, as one of the important initial conditions for
100
Chapter 3 Field Measurements, Network Model Development and Validation
inrush transients, can be identified. It can be interpreted that: the closing of circuit
breaker was made on phase C first at about 5 ms behind the negative-going zero
crossing of phase C line-to-ground voltage; the other two phases were energised both
with a closing time delay of about 1 ms. If the residual fluxes of the energised
transformers were zero, the closing time would have resulted in minimum peak of
inrush current on phase C. However, Figure 3-14 shows that the phase A was of
minimum inrush current peak, which could be due to the existence of residual flux in
the transformer core. Comparing the peak magnitudes of inrush currents with that
resulted from the second attempt presented in section 3.3, it can be seen the maximum
peak is twice as much as that shown in Figure 3-4.
The RMS variations of the voltages observed at substation I are shown in Figure 3-18
(Note that the RMS calculation for field measurement results was based on one power
frequency cycle window and refreshed half power frequency cycle). It can be seen that
the dip magnitudes were about 5.5 kV (2.3%), 14.2 kV (6%) and 9.35 kV (4%) for
phase A, B and C, respectively. In addition, it can be identified that the rate of the initial
voltage dip recovery was faster than that of the later stage, which will be shown that it
was due to the SVC compensation.
The three-phase currents observed at power feeder 1 are shown in Figure 3-19 in terms
of instantaneous current waveforms for a short-time range and in Figure 3-20 in terms
of RMS current waveforms for a long-time range. In addition, three-phase voltages
observed on power feeder 1 in terms of RMS value are shown in Figure 3-21; three-
phase currents observed at the circuit I-K are shown in Figure 3-22.
As can be seen in Figure 3-19, the magnetizing currents of GSU transformers T2 and T3
became gradually larger after transformer energisation, which indicates the initiation of
sympathetic inrush. The growth and the decay of the sympathetic inrush currents are
further illustrated in Figure 3-20 which shows that the sympathetic inrush currents took
about 2 seconds to reach their maximum magnitudes (the largest maximum magnitude
appeared in phase C, which is about 120 A in RMS) and the decay lasted more than 25
seconds. With the presence of the long duration sympathetic inrush, it can be seen in
Figure 3-21 that the system took a long time to achieve full recovery of three-phase
voltages.
102
Chapter 3 Field Measurements, Network Model Development and Validation
Figure 3-20 RMS sympathetic inrush current measured at power feeder 1 in Case E4
Similar to the one shown in Figure 3-16, the currents measured at circuit I-K shown in
Figure 3-22 can be divided into two parts: the first part is, before energisation, the
steady state load currents flowing through the circuit; the second part is, after
energisation, the steady state load currents superimposed by a portion of transformer
inrush currents.
The above field measurement results provide some insights of transformer inrush
current transients (in terms of initiation, waveform pattern, decay and peak magnitudes)
103
Chapter 3 Field Measurements, Network Model Development and Validation
and also show the possible influences of transformer inrush current on voltage dips (in
terms of the degrees of unbalance among three phases, dip magnitudes and recovery
trends). These results will be taken as the benchmarks for validating network models.
The lengths of the 400 kV double-circuit transmission lines in the South West system
can be referred to Figure 3-2. The longest one is between substation H and G, which is
about 97 km. They are based on the similar type of tower structure, as shown in Figure
3-23. Specifically, a, b and c define the horizontal distances between the phase
conductors and the tower central; g, d, e and f define the vertical heights of phase
conductor and earth wire relative to ground surface. They are all determined in
accordance to the tower designs shown in Table 3-3. In the table, some other
104
Chapter 3 Field Measurements, Network Model Development and Validation
Figure 3-23 Basic tower structure used in South West system [113]
105
Chapter 3 Field Measurements, Network Model Development and Validation
A
C B
L2/2 A A
G F B C
A B A B
L6/1 C C L6/1 C C
B A B A C B
L2/2 A A
B
B C C B
H L6 A A
B C
C B C B C B C B C B C B
L2 A A L2 A A L6/1 A A L6/1 A A L6 A A L6 A A
B C B C B C B C B C D B C
I K J E C
Figure 3-24 Tower designs and transposing schemes associated with double circuit lines
106
Chapter 3 Field Measurements, Network Model Development and Validation
RMS
voltage
FC TCR detector
VRMS
Vref
Regulator
Firing
pulse
generator
Control module
The control module, as shown in Figure 3-25, consists of three building blocks:
• RMS voltage detector
• Regulator
• Firing pulse generator
Continuously, the RMS voltage detector processes the phase voltages measured at the
400 kV busbar to provide the RMS value of three averaged line-to-line voltages (VRMS).
The VRMS is then compared with the Vref in the regulator to produce a benchmark. The
benchmark and the measured line voltages are the inputs for the firing pulse generator to
produce firing pulses for controlling the operation of bi-directional connected thyristors,
which is achieved by following the procedure shown in Figure 3-26.
First of all, the line voltage is compared with a zero threshold to produce two opposite
square waves with a 1 pu magnitude. In this way, the ranges between two adjacent
voltage zero crossings suitable to give firing pulses are identified. In the next step, both
square waves are processed by an integrator and a high pass filter, resulting in two
107
Chapter 3 Field Measurements, Network Model Development and Validation
triangle waves. Comparing these two triangle waves with a 0.005 threshold, the positive
parts of the two triangle waves are selected and then summed to produce another
triangle wave with a frequency of 100 Hz. This new triangle wave is then compared
with the benchmark obtained from the regulator: the parts of the triangle wave with
magnitude less than the benchmark result in zero output; the other parts result in outputs
of 0.5. The comparison gives a square wave that is processed by a high pass filter to
generate firing pulses. The variation of the benchmark will change the width of the
square wave so as to vary the timing of firing pulses. As a result, the reactive power
output of the TCR can be varied, which subtracts the output of the FC so as to vary the
output of the SVC.
Figure 3-26 Procedure for generating firing pulses to control bi-directional thyristors
3.5.5 Transformers
3.5.5.1 Modelling GSU transformers
The transformer modelling approach suggested and validated in [6] was used to model
GSU transformers. It consists of two parts: a linear BCTRAN object to represent
transformer short-circuit characteristics; a set of delta connected nonlinear inductors,
located at the Low-Voltage (LV) terminal of the BCTRAN object, to represent
transformer core saturation characteristics (this model is named here as BCTRAN+).
108
Chapter 3 Field Measurements, Network Model Development and Validation
The test report data of 345 MVA and 415 MVA GSU transformers are shown in Table
3-5 and Table 3-6, respectively. The BCTRAN object utilizes the nameplate data and
the short-circuit test data, including rated power, rated voltages, short-circuit
impedance, full load losses and winding connection, to model the transformers.
Open-circuit test data were used to derive saturation curves for the nonlinear inductors.
Type-96 nonlinear inductor (enables the setting of initial residual flux) was selected for
GSU transformers. The conversion procedure for obtaining hysteresis saturation curves
for type-96 nonlinear inductor is shown in Figure 3-27, which consists of four steps:
• Step 1: transform RMS voltage versus current data (Vrms-Irms) into peak flux-
current data (λpeak-ipeak) based on the algorithm presented in [57] (a MATLAB
program was specifically developed based on the algorithm for carrying out the
transformation);
• Step 2: curve fit the λpeak-ipeak data to obtain a piecewise nonlinear saturation
curve;
• Step 3: select a positive saturation point from the piecewise nonlinear curve for
ATP subroutine HYSDAT [58] to derive a hysteresis loop;
• Step 4: assign an additional point beyond the saturation point to set the
saturation inductance deduced from transformer air-core inductance, with the
effect of winding leakage inductance considered [52].
109
Chapter 3 Field Measurements, Network Model Development and Validation
Figure 3-27 Conversion to derive saturation curve for type-96 nonlinear inductor
Figure 3-28 shows the lower half hysteresis curves for the GSU transformers. The
saturation points used to define the corresponding major hysteretic loop are appointed;
beyond the saturation points, the final slopes of saturation curves are quantified by
∆λ/∆i.
It was assumed that the air-core inductance is equal to twice of transformer short-circuit
inductance. This assumption is reasonable, because, when an unloaded two winding
core type transformer is energised from the HV winding side (usually the outer
winding), the cross section area of the air-core cylinder enclosed by the HV winding
where the flux goes through under deep saturation is normally about twice the cross
section area of the gap between HV and LV windings where the flux goes through
during the short-circuit test.
To preliminarily check the accuracy of the models, simulation of open circuit test was
applied to the developed BCTRAN+ models for the GSU transformers. Results obtained
from testing 415 MVA and 345 MVA GSU transformer models are shown in Figure
3-29 and Figure 3-30, respectively. In both figures, the simulated no-load current and
110
Chapter 3 Field Measurements, Network Model Development and Validation
no-load losses as a function of magnetizing voltage (under 0.9 pu, 1 pu and 1.1 pu of
rated voltage) are compared with that given in the transformer test report, which show
good agreements between measured and simulation results.
Measurement Simulation
25 350
No-load current No-load losses
20 300
Current, A
Loss, kW
15 250
10 200
5 150
0 100
0.9 1.0 1.1 0.9 1.0 1.1
Voltage, p.u. Voltage, p.u.
Figure 3-29 Open circuit test results deduced from GSU transformer model (415 MVA)
compared to manufacture test results
Measurement Simulation
50 350
No-load current No-load losses
40 300
Current, A
Loss, kW
30 250
20 200
10 150
0 100
0.9 1.0 1.1 0.9 1.0 1.1
Voltage, p.u. Voltage, p.u.
Figure 3-30 Open circuit test results deduced from GSU transformer model (345 MVA)
compared to manufacture test results
Furthermore, the BCTRAN+ models were tested regarding inrush current calculation,
with results compared with those generated by the more advanced Hybrid Transformer
model (including 3-limb and 5-limb core models, as the core configuration of the
transformers was unclear). To do the comparison, the Y curves used in the
BCTRAN+ model (i.e., the curves shown in Figure 3-28) were also used to implement
the Hybrid Transformer models; the models were energised on the primary side at the
same switching angle against an ideal 400 kV voltage source, with residual flux
assumed to be zero. The energisation cases, including energising 345 MVA transformer
at phase A voltage zero-crossing, energising 415 MVA transformer at phase A voltage
zero-crossing and energising 415 MVA transformer at phase A voltage peak, were
simulated and the results comparison are shown in Figure 3-31, Figure 3-32 and Figure
3-33, respectively.
111
Chapter 3 Field Measurements, Network Model Development and Validation
Figure 3-31 Comparison of inrush currents generated by Hybrid and BCTRAN+ for 345 MVA
transformer (Energised at phase A voltage zero-crossing)
Figure 3-32 Comparison of inrush currents generated by Hybrid and BCTRAN+ for 415 MVA
transformer (Energised at phase A voltage zero-crossing)
Figure 3-33 Comparison of inrush currents generated by Hybrid and BCTRAN+ for 415 MVA
transformer (Energised at phase A voltage peak)
112
Chapter 3 Field Measurements, Network Model Development and Validation
As can be seen, the patterns and magnitudes of the three-phase inrush currents
calculated by the BCTRAN+ model can closely match those estimated by the more
advanced Hybrid Transformer model. Good agreements were achieved in the
energisation cases with different switching angles. This indicates that the relatively
simple core topology in the BCTRAN+ model is sufficient to give accurate simulation
results.
The way to model substation transformers is similar to that for modelling GSU
transformers, which is a BCTRAN object plus a set of delta connected nonlinear
inductors. The differences are that: type-93 nonlinear inductors were used for the core
representation of substation transformers (because there is no need to model their
residual flux condition since they are in service) and they were located at the tertiary
winding terminal. As the type-93 nonlinear inductor only requires single-value
saturation curve, the procedure to derive the saturation curve for type-93 nonlinear
inductor is simpler compared to that for type-96 nonlinear inductor. As shown in Figure
3-34, the basic procedure is the same with that shown in Figure 3-27 but with the step of
forming hysteresis loop skipped.
Figure 3-34 Conversion to derive saturation curve for type-93 nonlinear inductor
As shown in Figure 3-35, two saturation curves were derived from transformer test
report data: one is for conventional transformers (those manufactured before 1980); the
other is for modern transformers. This difference is due to the change of transformer
design, manufacture technique and core material around the time of 1980. The nominal
magnetizing current of modern transformer is much smaller than that of conventional
transformer, and therefore with lower core losses.
The air-core inductances of all the substation transformers were assumed to be four
times the transformer short-circuit inductance, as the substation transformers are
113
Chapter 3 Field Measurements, Network Model Development and Validation
1.25
1.00
Fluxlinkage, pu
0.75
0.50
0.25
Modern
Conventional
0.00
0 20 40 60 80 100 120 140
Current , A
In the field tests, the measurement of inrush currents can be influenced by the
performance of current transformer (CT). When the CT is saturated, the shape of inrush
currents given on the CT’s secondary side can be significantly different from those seen
on its primary side. A CT model was developed to simulate such a transient response.
In Figure 3-36, the circuit diagram of the CT model is shown. It is formed by three
single-phase transformer models. Each single-phase model consists of an ideal
transformer, a type-96 non-linear inductor and a lumped impedance. The ideal
transformer is with a ratio of 1200/5. The type-96 non-linear inductors were
characterized by the Y curve shown in Figure 3-37; the curve was derived based on
the procedure shown in Figure 3-27. The lumped impedance, whose value was assumed
to be 5 d0.62 Ohm, represented the relay burden and CT winding impedance.
114
Chapter 3 Field Measurements, Network Model Development and Validation
Flux-linkage, Wb.T
1
-1
-2
-3
-4
-0.5 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
Current, A
The simulation setup for each validation study consisted of two parts: the first part was
the initialization of network condition to define the values of system source voltages,
source strength impedance and system loading; the second part was the initialization of
energisation condition including the circuit breaker closing time and the transformer
residual flux.
Regarding the network condition, the network model was initialized via the following
procedure:
• The system source voltages were initialized by referencing the steady state
voltages (measured prior to energisation) at substation I: for simulating E1, E2
and E3, the applied source voltage (line-to-line) is 418 kV; for simulating E4,
the applied source voltage (line-to-line) is 414 kV;
• It was assumed that the synchronous compensator at substation H was not
connected;
• Equivalent source impedances for supply sources S1 and S2 were derived based
on the short-circuit levels: 7.1 GVA for S1 and 6.4 GVA for S2;
• Following the loading data provided by the network operator, the loading at each
substation was distributed in such a way that each substation transformer was
about half loaded: it should be noted that, due to loading variation, the loadings
115
Chapter 3 Field Measurements, Network Model Development and Validation
of substation H and G in the Case E1, Case E2 and Case E3 were set to be 10%
more than those in the Case E4.
As for the energisation condition, circuit breaker closing times were interpreted from
the measured voltage waveforms and transformer residual fluxes were all assumed.
Figure 3-38 Simulated voltages at power feeder 1 compared to those measured in Case E1
Figure 3-39 Simulated currents at power feeder 1 compared to those measured in Case E1
116
Chapter 3 Field Measurements, Network Model Development and Validation
Figure 3-40 Simulated currents at circuit I-K compared to those measured in Case E1
The voltage variation observed at power feeder 1, inrush current drawn by the T2 and
T3, current variation on circuit I-K and voltage variation at substation-I were all
simulated based on the network model, with the results plotted against the
measurements shown in Figure 3-42, Figure 3-43, Figure 3-44 and Figure 3-45,
respectively. As can be seen, the waveforms obtained from simulation closely agree
with test results.
Figure 3-42 Simulated voltages at power feeder 1 compared to those measured in Case E2
117
Chapter 3 Field Measurements, Network Model Development and Validation
Figure 3-43 Simulated currents at power feeder 1 compared to those measured in Case E2
Figure 3-44 Simulated currents at circuit I-K compared to those measured in Case E2
Circuit breaker closing time were interpreted from the voltage waveforms measured on
the power feeder 1: phase C was set to be energised first at 5 ms after the negative-
going zero crossing of phase C line-to-ground voltage; referencing the closing time of
phase C, the other two phases were energised both with a delay of 1.1 ms. Residual
fluxes in both T2 and T3 were initialized to -0.385 pu, 0.55 pu and -0.165 pu of peak
nominal flux for phase A, phase B and phase C, respectively.
The simulated currents and voltages on power feeder 1, currents flowing through circuit
I-K and also the RMS voltage variation at substation I are plotted against field
118
Chapter 3 Field Measurements, Network Model Development and Validation
measurement results in Figure 3-46, Figure 3-47, Figure 3-48 and Figure 3-49,
respectively. As can be seen, the simulated waveforms match the measurement results
well; specifically, in Figure 3-49, the part of voltage recovery with faster speed is also
replicated by considering the response of SVC in the network model.
Figure 3-46 Simulated voltages at power feeder 1 compared to those measured in Case E3
Figure 3-47 Simulated currents at power feeder 1 compared to those measured in Case E3
Figure 3-48 Simulated currents at circuit I-K compared to those measured in Case E3
119
Chapter 3 Field Measurements, Network Model Development and Validation
Figure 3-49 Simulated RMS voltage variation at substation I compared to those measured in
Case E3
During the validation against measurement E3, it was found that (as shown in Figure
3-50): good agreement between simulated and measured inrush currents (drawn by T2
and T3) can only be achieved in the first few cycles; after the sixth cycle, the simulated
inrush currents cannot replicate the abrupt changes of current peak and waveform
pattern of the measured inrush currents in phase B and phase C.
Figure 3-50 Comparison between measured and simulated inrush currents drawn by T2 and T3
in Case E3 (simulated currents observed at the CT primary side)
It is believed that the changes appeared in the measured phase B and phase C currents
were due to the CT saturation caused by the significant dc offset components presented
in both currents. No obvious change can be observed in the measured phase A current is
because its offset dc component was of relatively small magnitude, in which case, good
120
Chapter 3 Field Measurements, Network Model Development and Validation
agreement between simulated and measured currents can still be retained. To confirm
such reasoning, the developed CT model was incorporated in the network model to
reproduce the inrush currents altered by the CT saturation effects. The simulated inrush
currents shown in Figure 3-50 was taken as the input for the CT model, and the resulted
secondary side currents of the CT model were obtained and compared with the
measured inrush currents in Figure 3-51. As can be seen, the abrupt changes in phase B
and C can be very well matched, which indicates that the simulated inrush currents
shown in Figure 3-50 are accurate and confirms that the abrupt changes of the inrush
currents was due to CT saturation.
Figure 3-51 Comparison between measured and the simulated currents (observed at the CT
secondary side)
121
Chapter 3 Field Measurements, Network Model Development and Validation
Figure 3-52 Simulated sympathetic inrush currents at power feeder 1 compared to those
measured in Case E4
Figure 3-53 Simulated RMS sympathetic inrush currents at power feeder 1 compared to those
measured in Case E4
122
Chapter 3 Field Measurements, Network Model Development and Validation
Figure 3-54 illustrates the comparison regarding the long-duration RMS voltage dips,
with the dip magnitudes particularly compared in Table 3-7. As can be seen, the
simulated voltage dip recovery traces are similar to field measurement results; the
largest deviation of dip magnitude is approximately 2.8%.
The comparison in Figure 3-55 is the currents measured at circuit I-K. Good agreement
can be seen in the range of steady stage (i.e., prior to energisation), which confirms the
correct modelling of system loading. In the range of transient stage (i.e., after
energisation), simulation results also show good agreement with field test results. It
should be noted that the initial part of voltage recovery was affected by the response of
SVC, which has also been correctly replicated in the simulation results. Regarding this
SVC effect, more details will be described in Chapter 4.
Figure 3-54 Simulated RMS voltage at power feeder 1 compared to those measured in Case E4
Figure 3-55 Simulated currents at circuit I-K compared to those measured in Case E4
123
Chapter 3 Field Measurements, Network Model Development and Validation
3.7 Summary
In this chapter, voltage dip events caused by energising large GSU transformers is
reported, in conjunction with field measurements, network model development and its
validation. The system under study is a 400/132 kV transmission system featured by
long transmission lines between the supply source and the energised transformers.
In the voltage dip events, it was shown that the energisation of GSU transformers can
trigger a network-wide voltage dips; the recorded maximum RMS voltage dip was
about 7.85% of the initial voltage; the events triggered low voltage alarms and the dips
were noticed by the downstream distribution utilities.
Based on the system parameters provided by the network operator, a network model
was developed in ATP/EMTP by following the guidelines summarized in Chapter 2.
The network model was validated against the field measurement results, showing good
agreements. It thus confirms the accuracy of the following modelling approaches:
• The source network can be modelled by an ideal sine-wave source and a
Thevenin equivalent impedance of the part of the network not under study;
• The transmission network between the supply source and the energised
transformer should be represented in detail, taking into account the transmission
lines, system loading and reactive power compensation devices;
• The Bergeron model can be used to represent transmission lines, with line
dimension and transposing scheme considered;
• System loading can be represented by lumped constant resistance and
inductance connected in parallel;
• Transformers can be represented by the BCTRAN+ model (consisting of a
BCTRAN object with a set of delta-connected nonlinear inductors).
In the following chapters, computer simulation assessments of the voltage dips caused
by transformer energisation will be carried out based on the validated network.
124
Chapter 4 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Deterministic Approach
Being the first step of the assessment, different energisation conditions are analysed and
compared, considering circuit breaker closing time, transformer core residual flux, the
number of GSU transformers being energised and the number of GSU transformers
already connected. Based on the commonly agreed worst case energisation condition,
assessment of network-wide voltage dips is performed to investigate the voltage dip
patterns under non-outage, single-circuit outage and double-circuit outage situations.
Furthermore, sensitivity studies are carried out to identify the key influential
parameters. Finally, possible operational measures to control the voltage dips are
assessed and an optimized energisation procedure is proposed to effectively and
economically reduce the voltage dips caused by transformer energisation.
125
Chapter 4 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Deterministic Approach
regarded here (and also in the rest of the thesis) as the commonly agreed worst case
energisation condition.
It should be noted that, for quantifying the dip magnitude and duration, the following
two criteria are used in this chapter and also in the rest of the thesis:
• Dip magnitude: the dip magnitude is defined as the lowest RMS voltage among
three phases;
• Dip duration: the 3% dip was taken as the beginning and end threshold for
quantifying the voltage dip duration.
The results of the case studies are summarized in Table 4-1, with voltage dip
magnitudes and durations observed at the substation-I selected for the comparison,
because it is located closest to the substation-K which connects the generating plant.
Table 4-1 Voltage dips observed at substation I under different energisation conditions
GSU transformer closing Residual flux Voltage dips
Case time (pu) Magnitude (%) Duration (s)
Energised On-line
(ms) A B C A B C A B C
1 T1 -- 5 0 4.3 3.5 1.8 0.274 0.181 0
2 T1 -- 0 0 2.0 2.9 4.4 0 0 0.393
3 T1 T2&T3 0 0 2.0 2.9 4.4 0 0 0.631
4 T1 -- 0 -0.8 0 0.8 5.8 5 9.6 0.946 0.636 2.728
5 T1 T2&T3 0 -0.8 0 0.8 5.8 5 9.6 2.90 1.477 6.448
6 T2&T3 -- 5 0 7.8 6.4 3.3 0.906 0.601 0.006
7 T2&T3 -- 0 0 3.7 5.0 8.2 0.04 0.265 1.12
8 T2&T3 T1 0 0 3.7 5.0 8.2 0.04 0.279 2.614
9 T2&T3 -- 0 -0.8 0 0.8 10.7 9.2 18.4 1.576 1.199 3.564
10 T2&T3 T1 0 -0.8 0 0.8 10.7 9.2 18.4 4.0 2.68 7.919
Comparing Case 1 with Case 2, it can be seen that, energising at the peak of phase C
voltage only gives small voltage dip on phase C, but results in larger voltage dips on the
other two phases, which shows that, in the case of three-phase system, a closing time in
favour of one phase may act against the other two phases. This is also evidenced in the
comparisons between Case 6 and Case 7.
Comparing Case 2 and Case 4, it can be seen that, when the transformer is energised at
the commonly agreed worst case energisation condition, the resulted voltage dips are
much higher than those voltage dips resulted from energisation at zero-crossing
energisation instant but with zero residual flux. This can also been seen from other
comparisons, including the one between Case 3 and Case 5, the one between Case 7 and
Case 9 and the one between Case 8 and Case 10.
126
Chapter 4 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Deterministic Approach
Under the same closing time and residual flux condition, it can be seen that:
• Energising transformers T2 and T3 together would cause voltage dips with
longer dip duration (in most cases) and with magnitude almost twice that caused
by energising T1 alone, which, for instance, can be seen from the comparison
between Case 4 and Case 9;
• Comparison between Case 4 and Case 5 shows, with the presence of on-line
transformers, the voltage dip magnitudes are the same but the durations can be
prolonged by more than 100%; similar observation can be seen in the
comparison between Case 2 and Case 3, between Case 7 and Case 8, and
between Case 9 and Case 10.
The above case studies show that the sympathetic inrush does not affect the voltage dip
magnitude, which is obvious because at the initial stage of the energisation transient,
voltage dip is mainly determined by the inrush current in the energised transformers,
while at this time the sympathetic inrush drawn by the adjacent transformer has yet to
build up. With the increase of the sympathetic inrush current, the decay of the inrush
current in the energised transformer slows down and as a consequence the dip duration
is prolonged. The comparisons made between the case studies indicate that this
prolonging effect can be very significant. Comparing Case 4 with Case 5, the dip
duration of phase C is prolonged by 136%; similarly, comparing Case 9 with Case 10,
the dip duration of phase C is prolonged by 122%.
As far as the case studies of energising only one GSU transformer are concerned, Case 5
is the worst case energisation. Among all of the case studies, Case 10 is the worst
energisation case, as it results in the largest dip magnitude and the longest dip duration.
Both cases are further analysed in terms of the current and voltage variations.
The inrush current drawn by T1 after the energisation is shown in Figure 4-1. The
maximum current peak appears on phase C and it is about 2500 A. As shown in Figure
127
Chapter 4 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Deterministic Approach
4-2 and Figure 4-3, part of the inrush current flows through the double-circuit lines
between substation I and K and part of it flows through the lines between substation J
and K. The inrush current flowing through the lines between J and K (not observed in
the field measurement) is much larger than that flowing through the lines between I and
K, which is mainly due to the large impedance of the circuit between substation G and
H.
In Figure 4-6, the three-phase RMS voltage dips are illustrated. It shows that the voltage
dip recovery trace started to be altered once the sympathetic inrush currents started to
build up, which can be further evidenced by the comparison between Case 4 and Case 5
regarding the voltage dip recovery on phase C (see Figure 4-7).
Figure 4-2 Currents flowing through one of the circuits between substation I and K (Case 5)
128
Chapter 4 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Deterministic Approach
Figure 4-3 Currents flowing through one of the circuits between substation J and K (Case 5)
Figure 4-4 Initiation of sympathetic inrush current observed at power feeder 1 (Case 5)
Figure 4-5 RMS sympathetic inrush current observed at power feeder 1 (Case 5)
Figure 4-7 Comparison between Case 5 and Case 4 regarding phase C voltage dip
129
Chapter 4 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Deterministic Approach
The sympathetic inrush current drawn by transformer T1 alone (as shown in Figure 4-11
and Figure 4-12), is similar to that drawn by T2 and T3 together (see Figure 4-4 and
Figure 4-5 in the analysis of Case 5). This is because the total number of transformers
engaging sympathetic inrush in both cases is the same.
In Figure 4-13, the voltage dip of phase C observed in Case 10 is compared with that
observed in Case 5. As can be seen, the higher inrush current in Case10 results in
voltage dips with magnitude twice that observed in Case 5; the dip duration, however, is
of much smaller difference between the two cases, which can be attributed to the similar
level of sympathetic inrush for both cases.
Figure 4-9 Currents flowing through one of the circuits between I and K (Case 10)
130
Chapter 4 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Deterministic Approach
Figure 4-10 Currents flowing through one of the circuits between J and K (Case 10)
Figure 4-11 Initiation of sympathetic inrush observed at power feeder 2 (Case 10)
Figure 4-12 RMS sympathetic inrush current observed at power feeder 1 (Case 10)
Figure 4-13 Comparison between Case 5 and Case 10 regarding phase C voltage dip
131
Chapter 4 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Deterministic Approach
Figure 4-14 shows the patterns of voltage dip magnitudes of phase-C on both 400 kV
and 132 kV side. It can be seen that the two patterns are almost identical. This indicates
that, as far as dip magnitude is concerned, the voltage dips appear on the 132 kV side
are not much affected by the substation transformers in the system studied. The dip
magnitude observed at each substation is found to be related to the distance between the
substation and the supply source and also the distance between the substation and the
energised transformers. For those substations (including H, I and J) located in the
proximity of the energised transformer and relatively far away from the supply source,
the observed dip magnitudes are relatively large; for those substations (including D, E,
F and G) located relatively far away from the energised transformer and close to the
supply source, the observed dip magnitudes are relatively small.
400kV 132kV
20%
RMS voltage dip magnitude
18%
16%
14%
12%
10%
D E F G H I J
Substation
Figure 4-14 Patterns of voltage dip magnitudes across all the network substations (voltage dips
observed at substation autotransformers’ 400 kV side versus 132 kV side)
Figure 4-15 shows the patterns of voltage dip duration observed on 400 kV and 132 kV
side. It can be seen that the dip durations on 132 kV side are longer. The prolonged
voltage dip duration at 132 kV side is attributed to the sympathetic inrush of substation
132
Chapter 4 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Deterministic Approach
400kV 132kV
9
8
voltage dip duration, s
3
D E F G H I J
Substation
Figure 4-15 Patterns of voltage dip duration across all the network substations (400 kV side
versus 132 kV side)
Figure 4-16 Voltage dip recovery traces observed at 400 and 132 kV busbars of substation I
133
Chapter 4 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Deterministic Approach
To understand the impacts of single-circuit outage on the voltage dip pattern, various
single-circuit outage scenarios were simulated based on Case 10 energisation condition.
In the simulation studies, voltage dips in phase C were observed at substations D, E, F,
G, H, I and J and compared with those obtained under non-outage condition. The
relatively severe voltage dips scenarios have been found under the A-F and J-E single-
circuit outage, which are shown in Figure 4-17 and Figure 4-18. Note that in both
figures, the sub-figure on the left shows the comparison regarding the voltage dip
magnitude whilst the sub-figure on the right shows the comparison regarding the dip
duration. From the comparisons, it can be seen that: single-circuit outage only results in
slight increase of voltage dips; the maximum increase of dip magnitude and duration
both appear at substation J.
35% 16
Non-outage condition Non-outage condition
30% 14
Single-circuit outage between A and F Single-circuit outage between A and F
Voltage dip magnitude
12
25%
10
20%
8
15%
6
10%
4
5% 2
0% 0
D E F G H I J D E F G H I J
Substation Substation
(a) Dip magnitude (b) Dip duration
Figure 4-17 Voltage dips caused by single-circuit outage between substation A and F
35% 16
Non-outage condition Non-outage condition
30% 14
Single-circuit outage between J and E Single-circuit outage between J and E
Voltage dip duration, s
Voltage dip magnitude
12
25%
10
20%
8
15%
6
10%
4
5% 2
0% 0
D E F G H I J D E F G H I J
Substation Substation
(a) Dip magnitude (b) Dip duration
Figure 4-18 Voltage dips caused by single-circuit outage between substation J and E
134
Chapter 4 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Deterministic Approach
35% 16
Non-outage condition
Non-outage condition
14
30% Double-circuit outage between A and F
Double-circuit outage between A and F
Voltage dip duration, s
12
Voltage dip magnitude
25%
10
20%
8
15% 6
10% 4
5% 2
0% 0
D E F G H I J
D E F G H I J
Substation
Substation
(a) Dip magnitude (b) Dip duration
Figure 4-19 Voltage dips caused by double-circuit outage between substation A and F
135
Chapter 4 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Deterministic Approach
35% 16
Non-outage condition Non-outage condition
30% 14
Double-circuit outage between B and C Double-circuit outage between B and C
Voltage dip magnitude
0% 0
D E F G H I J D E F G H I J
Substation Substation
25%
10
20%
8
15%
6
10%
4
5% 2
0% 0
D E F G H I J D E F G H I J
Substation Substation
(a) Dip magnitude (b) Dip duration
Figure 4-21 Voltage dips caused by double-circuit outage between substation I and K
35% 16
Non-outage condition Non-outage condition
30% 14
Double-circuit outage between J and E Double-circuit outage between J and E
Voltage dip magnitude
12
25%
10
20%
8
15%
6
10%
4
5% 2
0% 0
D E F G H I J D E F G H I J
Substation Substation
136
Chapter 4 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Deterministic Approach
The energisation Case 10 was selected as the base case for the sensitivity assessment.
The sensitivity assessment addressed variations of the concerning parameters between
+100% and -50% of their values in the base case (i.e. up to 200% variation was applied).
The results of the sensitivity study are shown in Figure 4-23 in which voltage dips due
to variation of transformer load losses, core saturation inductance, source strength and
system loading are shown in Figure 4-23 (a), Figure 4-23 (b), Figure 4-23 (c) and Figure
4-23 (d), respectively.
In general, it can be seen that: both the core saturation inductance and the supply source
strength exhibit significant impact on dip magnitude; the dip duration is largely affected
by the variation of transformer load losses as well as the supply source strength; the
impact of system loading variation on voltage dip duration and magnitude is moderate.
(a) Effects of transformer load losses (b) Effects of transformer core saturation
variation inductance variation
(c) Effects of source strength variation (d) Effects of system loading variation
Figure 4-23 Voltage dips influenced by variation of key parameters
Furthermore, two sensitivity factors were used to quantify the effects of the parameter
variations on the voltage dip magnitude and duration: one is dip magnitude deviation Vd
137
Chapter 4 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Deterministic Approach
and the other is dip duration deviation Td. These two factors are presented in Figure
4-24 as an example which compares the phase C voltage dip between two cases: the
base case (plotting in solid line) and the case with 100% increase of source strength
(plotted in dotted line). The percentage of dip magnitude deviation is deduced from
100×(V1-Vbase)/Vbase; similarly, the percentage of dip duration deviation is calculated by
100×(T1-Tbase)/Tbase.
Tbase
T1 Td
Vbase
V1
Vd
Based on the results shown in Figure 4-23, the two sensitivity factors were analysed and
the obtained quantities of the deviations are detailed and compared in Figure 4-25 and
Figure 4-26, respectively. In both figures, each column represents the deviation caused
by variation of a certain parameter.
Figure 4-25 illustrates that change of core saturation inductance presents the biggest
impact on dip magnitude. Source strength is the second most influential one. In fact, in
the reduced cases, the impact of source strength variation on voltage dip magnitude is
comparable with that caused by variation of core saturation inductance. The variation of
system loading only results in minor impact on dip magnitude. The variation of
transformer copper losses barely impacts the dip magnitude.
In Figure 4-26, it can be observed that change of transformer load losses dominates the
impact on voltage dip duration (e.g., 50% reduction of transformer load losses could
increase dip duration by 88%); the cases of reduced full-load losses have more impact
than those of increased. The second most influential parameter is again the source
138
Chapter 4 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Deterministic Approach
impedance. 50% reduction of source strength will increase the duration by 36.5%.
Variation of system loading only shows minor impact on voltage dip duration, which is
similar to its impact on dip magnitude. The variation of core saturation inductance
shows very little impact on voltage dip duration, however, exhibits great impact on dip
magnitude. It should be noted that, although the above analysis focuses on phase C, the
other two phases also exhibit the same trend.
40%
Transfromer copper losses Core saturation inductance
Source strength System loading
30%
Dip magnitude deviation Vd
20%
10%
0%
-10%
-20%
-30%
-40%
+100%
100% +50%
50% +25%
25% -25%
-25% -50%
-50%
Percentage of parameter variation
Figure 4-25 Impacts of parameter variation on voltage dip magnitude
100%
Transfromer copper losses Core saturation inductance
Source strength System loading
80%
60%
Dip duration deviation Td
40%
20%
0%
-20%
-40%
-60%
+100%
100% +50%
50% +25%
25% -25%
-25% -50%
-50%
Percentage of parameter variation
Figure 4-26 Impacts of parameter variation on voltage dip duration
The minor impact of system loading variation suggests that the system loading
condition is of less concern when carrying out transformer energisation in the network
studied, which can be attributed to the low ratio of maximum substation loading to the
source strength (which is less than 0.09), i.e., the source strength is too strong relative to
139
Chapter 4 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Deterministic Approach
system loading. However, it is noteworthy that for those systems having a higher ratio
of system loading to source strength, the effect of system loading might become more
obvious. The network studied here is characterized by long transmission lines between
supply source and the transformers being energised, and therefore the effect of source
strength, to some extent, has been offset by the impedances of long transmission lines;
in those systems where the supply source is located closer to the transformers being
energised, the impact of source strength variation would be more significant. The
variation of core saturation inductance directly influences the magnitude of inrush
currents and thus produces great impact on the dip magnitude; however, most
transformers up to date are only tested at factory up to 110%; open circuit test at higher
voltages is required for more accurate estimation of core saturation inductance. On the
other hand, the saturation inductance only slightly affects the decay time constant of
inrush transients (due to the saturation inductance relatively small compared to network
impedance), and therefore shows little impact on the dip duration. The significant
impact of transformer load losses on dip duration is because the decay of sympathetic
inrush is highly determined by the losses of the GSU transformers and the losses of the
connection between them. Due to the short electrical distance between the GSU
transformers, the amount of load losses of the GSU transformers is the key contributor
to the damping of sympathetic inrush and therefore imposes the greatest impact on dip
duration.
140
Chapter 4 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Deterministic Approach
125 125
100 100
75 75
50 50
Flux-linkage, Wb. T
Flux-linkage, Wb. T
25 25
0 0
-25 -25
-125 -125
-150 -150
-20 0 20 40 60 80 100 120 140 160 180 200 -20 0 20 40 60 80 100 120 140 160 180 200
Current, A Current, A
141
Chapter 4 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Deterministic Approach
25%
Base case Maximum tap
20%
10%
5%
0%
D E F G H I J
Substation
Figure 4-28 Voltage dip magnitudes observed in the case with GSU transformers set to their
maximum tap
12
Base case Maximum tap
10
Voltage dip duration, s
0
D E F G H I J
Substation
Figure 4-29 Voltage dip duration observed in the case with GSU transformers set to their
maximum tap
In this study, it was assumed that the location of the SVC is at the 400 kV busbar of
substation K and its response time is fixed at 120 ms. With different capacities,
including 75 MVar, 150 MVar and 300 MVar, the effect of the SVC capacity on voltage
142
Chapter 4 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Deterministic Approach
dip performance was simulated. Regarding the effects on voltage dip magnitude, the
results obtained are shown in Figure 4-29; regarding the effects on voltage dip duration,
the results obtained are shown in Figure 4-31; both were compared with the patterns
observed in the base case. As can be seen, application of SVC device has negligible
impacts on dip magnitude but helps speed up the voltage dip recovery; the recovery can
be observed on all the substations; faster recovery can be achieved by the SVC with
larger capacity; applying 300 MVar SVC can result in 39% reduction of dip duration.
The negligible impacts on dip magnitude can be attributed to the fact that the SVC
normally needs several cycles to respond to the voltage dips within which maximum dip
magnitude has already been reached.
25%
Base case 75 MVar
20%
150 MVar 300 MVar
Voltage dip magnitude
15%
10%
5%
0%
D E F G H I J
Substation
12
Base case 75 MVar
10
150 MVar 300 MVar
Voltage dip duration, s
0
D E F G H I J
Substation
With the capacity of SVC fixed at 300 MVar and its location at the 400 kV busbar of
substation K, possible effect of response time variation on the SVC performance was
studied. The values of response time under study include 80 ms, 120 ms and 160 ms. It
143
Chapter 4 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Deterministic Approach
is shown in Figure 4-32 that the effect of SVC on voltage dip duration is not affected by
the variation of SVC response time.
With the capacity of SVC fixed at 300 MVar and its response time at 120 ms, further
studies were conducted to see whether the location of SVC may affect the dip duration.
Besides locating the SVC at substation K, two other locations were considered:
substations E and H. In Figure 4-33, the obtained patterns of dip duration are contrasted
to that observed in the base case. As can be seen, the substations nearest to the SVC
location will benefit the most; therefore it is desirable to locate SVC close to the
substation where the transformers are to be energised.
12
Base case 80 ms 120 ms 160 ms
10
Voltage dip duration, s
0
D E F G H I J
Substation
Figure 4-32 Effects of SVC with different values of response time on dip duration
12
Base case Locate SVC at H
10
Locate SVC at E Locate SVC at K
Voltage dip duration, s
0
D E F G H I J
Substation
Figure 4-33 Patterns of voltage dip duration at 400kV side for various SVC locations
144
Chapter 4 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Deterministic Approach
25%
Base case Opening coupler CB1
20%
Voltage dip magnitude
15%
10%
5%
0%
D E F G H I J
Substation
Figure 4-34 Effects of opening coupler CB1 on dip magnitude
12
Base case Opening coupler CB1
10
Voltage dip duration, s
0
D E F G H I J
Substation
Figure 4-35 Effects of opening coupler CB1 on dip duration
145
Chapter 4 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Deterministic Approach
As can be seen, opening the coupler CB1 moderately reduces the dip magnitudes (about
21%) but significantly reduces the dip duration (nearly 75%). This great impact on
duration is largely attributed to the fact that: the opening of CB1 adds long transmission
lines into the connection between the GSU transformer T1 and GSU transformers
T2&T3, which effectively increases the losses of the connection between the GSU
transformers. As the decay of sympathetic inrush is highly determined by the losses of
the connection between the GSU transformers and their own losses, the increase of
connection losses significantly reduces the dip duration.
25%
Base case Combined
20%
Voltage dip magnitude
15%
10%
5%
0%
D E F G H I J
Substation
Figure 4-36 Dip magnitude pattern simulated under combined case
146
Chapter 4 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Deterministic Approach
12
Base case Combined
10
0
D E F G H I J
Substation
Figure 4-37 Dip duration pattern simulated under combined case
4.5 Summary
In this chapter, the voltage dips due to energising large GSU transformers was
comprehensively studied based on a real 400/132 kV network, with special attention
paid to the influences of sympathetic inrush. The network model developed and
validated in Chapter 3 was used for carrying out the simulation studies.
The degrees of voltage dips under different energisation conditions were assessed. It
was found that, under the worst case energisation condition: with two GSU transformers
simultaneously energised, the maximum dip in the present system is at the substation
closest to the transformers being energised and its magnitude is about 18% and duration
about 3.5 seconds; with the presence of sympathetic interaction, the magnitude of
sympathetic inrush current can be more than twice the transformer rated current; the
sympathetic inrush can prolong the dip duration from 3.5 seconds to 7.9 seconds. The
voltage dips caused at 400 kV transmission system side can propagate into 132 kV
distribution system and the dip duration on 132 kV side is longer due to the sympathetic
inrush of substation transformers.
Network-wide voltage dips were assessed under non-outage network condition. The dip
magnitude observed at each substation was found to be related to the distance between
the substation and the supply source and also the distance between the substation and
the energised transformers. Those substations located in the proximity of the energised
transformer and relatively far away from supply source would experience larger dip
147
Chapter 4 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Deterministic Approach
magnitudes.
Furthermore, the network-wide voltage dips were assessed when the network under line
outage situations. It was found that: under single-circuit outage, the network voltage
dips performance is similar to that observed under non-outage network condition;
however, if there is double-circuit outage, both dip magnitude and duration can be
significantly exacerbated; in the present system under study, the most unfavorable
double-circuit outage can increase the dip magnitude to about 30% and the duration to
about 14 seconds.
Considering sympathetic inrush and voltage dips are controlled by multiple parameters,
sensitivity assessment was carried out to identify the most influential parameters on dip
magnitude and duration. It was found that transformer core saturation inductance has
the most profound impact on the voltage dip magnitude; however, this parameter is not
readily available from factory test report. The amount of transformer load losses has
been proven as the key influential parameter on determining the duration of sympathetic
inrush and voltage dip.
148
Chapter 5 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Stochastic Approach
However, since the circuit breaker closing time, transformer core residual flux as well
as system conditions are normally of stochastic nature, the assessments performed using
the deterministic approach as mentioned above could underestimate the outcome or at
least unable to give the probability distribution of the occurrence of inrush transients,
therefore, they cannot assist the realistic estimation of the adverse impacts of inrush
transients on the system. Hence, recently, a number of studies were devoted to studying
the impacts of parameter uncertainties on the calculation of transformer energization
transients [43, 104, 114]. The studies extended the few deterministically-defined studies
to many stochastically-defined case studies by considering the stochastic nature of the
key parameters.
In this chapter, Monte Carlo approach is used to stochastically evaluate voltage dips
caused by transformer energisation in the South West system. The 400 kV transmission
network model developed and validated in Chapter 3 is coordinated with MALAB
program to perform the Monte Carlo simulation. It is aimed to identify the dip
149
Chapter 5 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Stochastic Approach
frequency pattern and the likelihood of reaching the dip magnitude resulted from the
commonly agreed worst case energisation condition, with their sensitivities to the
variation of circuit breaker closing time span, transformer core residual flux, system
condition and the number of transformers being energized together also investigated.
Arbitrary number of simulation runs can be carried out. For each simulation run, the
coordination of MATLAB and ATP for Monte-Carlo simulation follows such a
procedure:
• First, the random number generators in MATLAB are called upon to generate
values for the stochastic variables such as closing time and residual flux;
• Second, MATLAB opens the ATP program and writes the generated values of
150
Chapter 5 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Stochastic Approach
Supposing detailed design data are available, transformer core saturation inductance can
be determined by the core and winding dimension; transformer copper losses can be
determined by the winding material, cross section area and length, which indicates they
are fixed by transformer design. Therefore, the two parameters can be regarded as
deterministic parameters.
As mentioned in Chapter 2, closing time and residual flux are normally of stochastic
nature. In general, when the closing of a circuit breaker is uncontrolled, the signal to
initiate breaker closing is stochastic with respect to the ac voltage wave; in addition, the
breaker poles of three-phase do not close simultaneously, but with some closing time
scatters which may vary with time and maintenance [116]. Hence, the circuit breaker
closing time is featured by the stochastic signalling time to order the closure of three
poles and the stochastic closing time span between the three poles.
Transformer residual flux is largely influenced by the ring down process initiated by
transformer de-energisation [117]. Resulted from this process, the residual flux in the
transformer core can be influenced by a number of factors including de-energisation
instants, circuit breaker chopping characteristics, transformer core material, winding
capacitances and other system capacitance connected to the transformer [118]. The
151
Chapter 5 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Stochastic Approach
uncertain properties of transformer residual flux can be categorized by two folds: one is
the magnitudes of residual flux in three phases and the other is their distribution among
three phases. Up to date, knowledge about these two aspects include: residual core flux
can be as high as 85% of peak nominal flux; residual flux in cores of three phase
transformers must inherently sum to zero [117, 118].
The amount of loads connected to the system could have daily and seasonal variations.
Similarly, the source strength may vary with system configurations and generation
connections. Therefore, to certain extend, both of these two parameters are of stochastic
feature.
Normally, Uniform distribution is used to describe the ranges of source strength and
system loading; Uniform distribution can also be used to describe ranges of switching
angle and residual flux in the case of single phase system. In the case of three-phase
system, two procedures shown in Figure 5-1 and Figure 5-2 are used for quantifying the
closing time and residual flux.
Utilizing the approach suggested in [43], modelling of the stochastic closing time
followed the procedures shown in Figure 5-1. First, a common order time T_cot, which
is the same for three poles of the circuit breaker, is defined by a Uniform distribution
ranging in one power frequency cycle. Second, the maximum closing time span
(MCTS), i.e. the time interval between the first pole and the last pole to close, is defined.
Third, the MCTS is used to define the range of the closing offset time for each pole
T_offset_i (i represents phase A, B or C); the range is assumed to be from –MCTS/2 to
+MCTS/2, referring to the T_cot; in addition, a probability distribution is assigned to
define the offset time range. Finally, the closing time of each pole can be determined by
the summation of T_cot and T_offset_i which are stochastically generated based on
their corresponding ranges and distributions.
152
Chapter 5 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Stochastic Approach
Figure 5-1 Procedure for generating stochastic circuit breaker closing time
For generating three-phase residual fluxes, it is commonly assumed that: the residual
flux of each phase is in a range whose absolute maximum value is no more than peak
nominal flux; the three-phase residual fluxes should sum to zero. With these
assumptions, the three-phase residual fluxes were generated through the procedure
shown in Figure 5-2. As can be seen, the first step is to define the maximum residual
flux Resi_max (normally in terms of the percentage of the peak nominal flux) for
determining the residual flux range ( i.e. from –Resi_max to +Resi_max); in the second
step, a probability distribution is assigned to characterize the residual flux range; in the
third step, based on the range and probability distribution, two residual flux values
Resi_1 and Resi_2 are stochastically generated; then, a check loop is called upon to
verify whether the absolute value of the sum of Resi_1 and Resi_2 is smaller than
Resi_max, which is to ensure that |Resi_3| = |-Resi_2-Resi_1| < Resi_max; if not, go
back to the second step; if yes, proceed to the fifth step to calculate Resi_3; finally,
Resi_1, Reis_2 and Resi_3 are stochastically assigned to phase A, B and C, respectively.
153
Chapter 5 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Stochastic Approach
Figure 5-2 Procedure for generating stochastic transformer core residual flux
154
Chapter 5 Assessment of Voltage
Voltage Dips Caused by Transformer Energisation Transients Using Stochastic Approach
Figure 5-3 Single phase simulation circuit for preliminary Monte Carlo simulation
Initially, a determinstic worst case simulation was carried out based on the commonly
agreed worst case energisation condition (i.e., the circuit breaker was closed at the
postive-going zero crossing of the applied voltage and the transformer core retained
retain 0.8
pu of peak nominal flux with positive polarity),
polarity) with the system condition set as: the
source fault level is of 9.5 GVA
G (X/R ratio is 10), the system loading is of 1800 MVA
(power factor is 0.95 inductive).
inductive) The resulted inrush currents and voltage dips observed
at the transformer primary terminal are shown in Figure 5-4 and Figure 5-5, respectively.
As can be seen, the
he first peak of the inrush current, the voltage dip magnitude and
duration were 4.63 kA,
k 12% and 1.02 s, respectively. These were taken as the base
values to scale the results obtained from the following case studies within this section.
Figure 5-4 Inrush current resulted from energising a single phase transformer under the worst
energisation condition
155
Chapter 5 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Stochastic Approach
Figure 5-5 Voltage dips resulted from energising a single phase transformer under the worst
energisation condition
Further, stochastic simulaiton was performed based on the single phase circuit. The
simulation consists of two cases, which are shown in Table 5-1.
Table 5-1 Case studies of stochastic estimation of voltage dips caused by energising a single
phase transformer
In Case P1, only the circuit breaker closing time and the transformer core residual flux
were considered as stochastic variables; both of them were defined by Uniform
distribution; the range for the closing time was one power frequency cycle and the range
for the residual flux was between -0.8 and +0.8 pu of peak nominal flux; the system
conditions were fixed at the level identical to those set in the worst case simultion (i.e.,
9.5 GVA source strength, 1800 MVA loading).
In Case P2, besides considering the variation of stochastic closing time and residual flux
(range and distirbution of the two variables were considered to be the same as in Case
P1), the variation of system condition (e.g., source strength and system loading) was
also considered and its range was considered to be ±25% (i.e., the source strength varied
in a range from 75 to 125% of 9.5 GVA and the system loading varied in a range 75 to
125% of 1800 GVA); both ranges were charaterized by Unifrom distribution. It should
be noted that the X/R ratio of the source impedance and the power factor for the system
loading were assumed to be constant.
156
Chapter 5 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Stochastic Approach
For both case studies, 1000 runs were conducted. The distributions of closing time and
residual flux obtained from Case P1 are shown in Figure 5-6 and Figure 5-7,
respectively.
157
Chapter 5 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Stochastic Approach
the inrush current peak, the larger the dip magntiude. There are two worst case
boundaries: the horizontal one is set by the worst case dip magnitude and the vertical
one is set by the highest inrush current first peak. None of the results obatined from
Case P1 study can exceed the two boudaries.
Figure 5-8 Relative voltage dip magnitudes plotted against relative inrush current peaks
Similar observations can be made in Figure 5-9 regarding the relation between the dip
duration and the inrush current peak; the section of zero dip duration is attributed to
those dips with magnitudes samller than the 3% dip threshold.
Figure 5-9 Relative voltage dip magnitudes plotted against relative inrush current peaks
Figure 5-10 illustrates the distribution of the relative voltage dip magntidutes. As can be
seen, about 50% of the dips are with magntiudes smaller than 0.2 pu of the worst case
dip magnitude; only about 6% of the dips are with magnitudes larger than 0.8 pu of the
worst case dip magnitude. The distribution of the relative voltage dip durations is shown
in Figure 5-11. Benchmarking the 3% dip threshold, about 55% of the dips are with
158
Chapter 5 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Stochastic Approach
duration equal to zero and less than 5% of the dips are with duration larger than 0.8 pu
of the worst case dip duration.
The resulted voltage dips of Case P2 stochastic simulation are illustrated as follows
(similar to Case P1, the obtained dip mangitudes, durations and inrush peaks are in per
unit values by scaling to 12%, 1.02 s and 4.63 kA, respectively):
• Scatter diagrams: relative dip magnitudes and relative dip durations plotted
against the relative inrush current peaks, which are shown in Figure 5-12 and
Figure 5-13, respectively;
• Bar charts: showing the frequency of different dip magntiudes and durations,
which are shown in Figure 5-14 and Figure 5-15, respectively.
In Figure 5-12, the scattered points indicate that the same inrush current peak may result
in a variaty of voltage dip magnitudes, mainly due to the vairation of system condition;
the width of the scatter range is proportional to the inrease of inrush current peak.
Benchmarking the scattered points with the worst case boundaries, most of the points
are capped inside; a few points are out of the voltage dip mangitude boundary in the
159
Chapter 5 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Stochastic Approach
cases of weak source strength; yet none of them can be of magnitudes exceeding both
boundaries.
Figure 5-12 Relative inrush first peaks plotted against relative dip magnitudes
Similarly, the relative inrush current first peaks are ploted against the realtive voltage
dip durations, as shown in Figure 5-13. It shows that: when the inrush current peak is
smaller than 0.3 pu of the worst case inrush current peak, the dip duration is zero, i.e.,
the voltage dip mangitude is smaller than the 3% threshold; when the inrush current
peak is bigger than 0.3 pu of the worst case, the voltage dip duration increases with the
inrush current peak magnitude; again, due to the vairation of system condition, the
possible voltage dip durations corresonponding to one inrush current peak is scattered.
The width of the scatter range is also propotional to the inrease of inrush current peak.
Figure 5-13 Relative inrush first peaks plotted against relative dip durations
Figure 5-14 shows the frequency of occurrence of different voltage dip mangitudes. As
can be seen, the proportion of dips with magnitudes less than 0.2 pu of the worst case
dip magnitude is about 65% which is larger than that observed in Case P1; only about 6%
of dips are of magnitudes larger than 0.8 pu of the worst case dip magnitude.
160
Chapter 5 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Stochastic Approach
Figure 5-15 shows the frequency of occurrence of different voltage dip durations. As
can be seen, about 70% of the dips are of zero duration; the proportion of dips with
duration reaching the worst case level is less than 1%.
The above findings suggest that the chance to encounter a voltage dip with a scale
matching that estimated by the commonly agreed worst case energisation condition is
very low. This can be attributed to the fact that, for any value of residual flux, there
exists a counteracting closing time that results in nearly zero inrush current, whereas the
probability is small for a maximum residual flux with a right polarity to meet a closing
time at voltage zero crossing.
Figure 5-14 Frequency plot of votlage dip magnitudes relative to the worst case dip magnitude
Figure 5-15 Frequency plot of votlage dip duration relative to the worst case dip magnitude
The preliminary study on the single phase circuit illustrates the influence of random
parameter variation on the outcomes of the inrush transient and the associated voltage
dips. It forms a basis for carrying out stochastic assessment for three-phase circuits.
161
Chapter 5 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Stochastic Approach
162
Chapter 5 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Stochastic Approach
transformer (T1) with the closing time and residual flux treated as stochastic
parameters and system condition (e.g., system loading and source strength)
fixed;
2) Influences of closing time on the dip frequency pattern are studied in Cases S2,
S3, S4 and S5; specifically, Cases S2 and S3 consider the influence of the
MCTS, Cases S4 and S5 study the influence of the closing offset time
distribution;
3) Influences of residual flux distribution on the dip frequency pattern are studied
in Cases S6, S7 and S8;
4) Influences of the variation of system condition including system loading and
source strength are considered in Case 9;
5) In Cases S10 and S11, the dip frequency patterns resulted from simultaneously
energizating two transformers are assessed and compared with that obtained in
the base case.
For all the case studies, the common order time is of one power frequency cycle range
defined by Uniform distribution and the residual flux is of a range between -0.8 pu and
+0.8 pu of peak nominal flux.
The dip magnitude observed in the case of energising T1 under the commonly agreed
worst case energisation condition (9.6%, named here as WCDM1) is selected as the
163
Chapter 5 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Stochastic Approach
base value to be referred to by all the dip magnitudes obtained from Case S1 to S9.
Similarly, the dip magnitude observed in the case of energising T2 and T3 under the
commonly agreed worst case energisation condition (18.4%, named here as WCDM2)
is used as the base value for scaling all the dip magnitudes obtained from Case S10 to
S11.
As can be seen, the dip frequency patterns of each phase are almost identical to each
other, indicating that the dip frequency pattern observed on one of the phases can
represent those observed on the other two phases. Therefore, analysis of the simulation
results can be focusing on one phase and the phase chosen here is phase C.
By observing the dip frequency of phase C shown in Figure 5-18, it can be seen that:
out of 1000 stochastic dip events, over 80% of the dips are with magnitudes less than
0.6 pu of WCDM1; no more than 6% of the dips are with magnitudes larger than 0.8 pu
of WCDM1; only about 0.2% of the dips are with magnitudes larger than the WCDM1
and their magnitudes are about 1.1 pu of WCDM1. This larger dip magnitude suggests
that calculation based on the commonly agreed worst case energisation condition
164
Chapter 5 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Stochastic Approach
(which assumes zero closing time span) may underestimate the worst case dip
magnitude by 0.1 pu.
Figure 5-16 Distribution of offset closing time for three-phase poles in Case S1
Figure 5-18 Frequency of dip magnitude of each phase at substation I out of 1000 stochastic
runs
165
Chapter 5 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Stochastic Approach
Figure 5-19 Frequency of dip duration in each phase at substation I out of 1000 stochastic runs
Figure 5-20 Frequency of dip magnitude in each phase at substation I out of 5000 stochastic
runs
Figure 5-21 Frequency of dip duration in each phase at substation I out of 5000 stochastic runs
166
Chapter 5 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Stochastic Approach
case studies were conducted: one with zero closing time span (Case S2) and the other
one with 10 ms MCTS (Case S3). For both cases, the distributions of closing offset time
and residual flux are the same as those used in the Case S1. Similar to Case S1, 1000
runs were made for Case S2 and Case S3 to predict the dip frequency patterns of phase
C at substation I; these two patterns are compared with that obtained from Case S1 in
Figure 5-22 (It should be noted that the comparison made on phase C can represent the
comparison made on the other two phases, because the dip frequency patterns of three
phases are almost identical to one another). As can be seen, there is not much difference
between any two of the dip frequency patterns. This indicates that the dip frequency
pattern is not sensitive to the variation of MCTS.
Figure 5-22 Frequency of voltage dip magnitude in phase C at substation I under different
values of closing time span
167
Chapter 5 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Stochastic Approach
5-23 and Figure 5-24, respectively. Obviously, compared with closing offset time with
Guassian distribution, the samples of larger offset time are largely increased in the case
of Exponential distribution.
Figure 5-23 Uniform closing offset time distribution within ±2.5 ms range
Figure 5-24 Exponential closing offset time distribution within ±2.5 ms range
Prediction of the dip frequency pattern of phase C at substation I with Uniform and
Exponential closing offset time was conducted in Case S4 and Case S5. The results are
compared with that obtained in Case S1 in Figure 5-25. As can be seen, the dip
frequency patterns are very similar to one another; the dip magnitudes are mostly
concentrated in the range between 0.2 and 0.4 pu of WCDM1; the frequency of dip
magnitudes larger than 0.8 pu of WCDM1 is less than 6% out of 1000 runs.
168
Chapter 5 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Stochastic Approach
Figure 5-25 Frequency of voltage dip magnitudes in phase C at substation I for different closing
time span distributions
169
Chapter 5 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Stochastic Approach
The dip frequency patterns observed on phase C at substation I were generated under
the above three residual flux distributions, which were compared with that obtained
under Case S1 in Figure 5-29. As can be seen, the dip frequency pattern resulted from
Gaussian residual flux distribution is almost the same with that given by Case S1; in the
case of Exponential_1, the frequency of dips with magnitudes less than 0.6 pu of
170
Chapter 5 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Stochastic Approach
WCDM1 is increased to about 90%, whilst the frequency of dips with magnitudes larger
than 0.8 pu of WCDM1 is reduced to less than 2%; residual flux with Exponential_2
distribution increases the frequency of dips between 0.8 and 1 pu of WCDM1 from 5%
to 11%, but it does not result in substantial increase of the dips with magnitudes
exceeding 1 pu of WCDM1. These findings suggest that the dip frequency pattern is
sensitive to the distribution of residual flux. For transformers prone to retain residual
flux of high magnitudes, the frequency of dips with magnitudes close to the worst case
dip magnitude is higher.
Figure 5-29 Frequency of voltage dips in phase C at substation I for different residual flux
distributions
The dip frequency pattern observed on phase C at substation I was calculated and
compared with that of Case S1, as shown in Figure 5-30. As can be seen, the two
171
Chapter 5 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Stochastic Approach
patterns are very similar to each other. This indicates that, the dip frequency pattern
estimated by Case S1 is not sensitive to ±25% variation of system condition.
Figure 5-30 Frequency of voltage dips in phase C at substation I (comparing Case S9 with S1)
The dip frequency patterns for both cases were obtained after 1000 runs. The dip
frequency pattern of phase C at substation I obtained from Case S10 simulation is
compared with that of Case S1 in Figure 5-31. Figure 5-32, similar comparison between
Case S11 and Case S1 is given. It can be seen that: if the transformers being energised
172
Chapter 5 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Stochastic Approach
simultaneously are of the same residual flux, the dip frequency pattern is identical to
that observed in the case of energising one transformer only; if the transformers are of
stochastically different residual flux, the frequency of dips with magnitudes between 0.2
and 0.6 pu of worst case dip magnitude is increased, whilst the frequency of dips with
other magnitudes decreased, which indicates that the likelihood of reaching the worst
case dip magnitude is reduced.
Figure 5-31 Frequency of voltage dips in phase C at substation I of Case S10 contrasting with
that of Case S1
Figure 5-32 Frequency of voltage dips in phase C at substation I of Case S10 contrasting with
that of Case S1 (two transformers with different residual flux)
5.10 Summary
In this chapter, Monte Carlo simulation has been conducted to stochastically assess the
voltage dips caused by transformer energisation in the South West Peninsula 400 kV
173
Chapter 5 Assessment of Voltage Dips Caused by Transformer Energisation Transients Using Stochastic Approach
grid, using the network model developed and validated in Chapter 3. The simulation
was automated by an ATP-MATLAB interfacing platform which uses ATP to handle
transient calculation and MATLAB to generate simulation inputs, control Monte Carlo
runs and process results.
A dip frequency pattern was produced over 1000 stochastic runs and it was found to be
sensitive to the distribution of residual flux but insensitive to the distribution of closing
offset time. This suggests that it is important to model the residual flux distribution in
transformer core while closing offset time distribution can be of less concern. In
addition, it was shown that the dip frequency pattern is insensitive to the system
condition when varying in a range of ±25% of the base case condition. The voltage dip
frequency pattern can be extended to cover the condition in which a number of
transformers are being energised simultaneously.
Furthermore, it was found that the probability of reaching the worst case dip magnitude
(estimated by the commonly agreed worst case energisation condition) is lower than
0.5%, indicating that the worst case scenario is unlikely to occur in a system; in fact,
about 80% of the dips are likely to be with magnitudes lower than 0.6 pu of the worst
case. Nevertheless, it was shown that there are dips with magnitudes exceeding the
worst case dip magnitude, indicating the inadequacy of the deterministic assessment
approach by using the commonly agreed worst case energisation condition.
174
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
To acquire more green energy, increasing installations of larger offshore wind farms are
being designed and commissioned. The total installed capacities of offshore wind farm
in counties such as Germany, UK and China are planned to be 10 GW, 18 GW and 30
GW by 2020, respectively [119-121]. The capacities of individual offshore wind farm
and wind turbine are also expanding: offshore wind farms with installed capacity
reaching 1 GW have been proposed [122]; wind turbines with rated capacity reaching
10 MW are commercially available [123].
Typical electrical system for an offshore wind farm involves a collection grid within the
wind farm and a transmission system to deliver the power to the onshore main grid. The
collection grid begins from the wind turbine transformers (usually at the base of the
wind turbine tower) which steps up the generation voltage from typically 690 V to a
medium voltage of 25-40 kV [124]. Depending on the wind farm total capacity and the
capacity of individual wind turbine, a large offshore wind farm may accommodate
dozens or even more than a hundred wind turbine transformers, as can be seen from the
top ten existing offshore wind farms (up to 2012) listed in Table 6-1 [125]. Usually, the
wind turbine transformers are distributed over a number of cable feeders; each feeder
may contain 5-10 wind turbines.
During the energisation of wind farm collection grid, there are two potential factors to
be considered: one is the possible voltage dip experienced at the point-of-common-
coupling between the electrical system of the wind farm and the utility company [75],
which concerns grid code compliance; the other is the sympathetic interaction between
wind turbine transformers [126], which on one hand may prolong the resulted voltage
dips and on the other hand prolong the mechanical and thermal stresses imposed on the
wind turbine transformers (according to the IEC 60076-16:2011 standard [127], due to
frequent energising wind turbine transformers during wind farm operation, wind turbine
transformers can be frequently exposed to mechanical and thermal stresses of inrush
175
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
currents). Both factors have not been thoroughly addressed for large offshore wind farm
connections.
In this chapter, modelling and simulation methodologies gained from the studies
conducted in Chapter 3 and Chapter 4 are used to develop network model of an existing
large offshore wind farm collection grid, based on which, deterministic and statistical
approaches presented in Chapter 4 and Chapter 5 are both employed to assess voltage
dips and sympathetic inrush caused by energising wind turbine transformers under
various scenarios. The assessment aims to provide guidance on planning and operating
offshore wind farms, especially focusing on identifying an optimum energisation
sequence to reduce sympathetic interaction between wind turbine transformers.
Table 6-1 List of top ten operational offshore wind farms [125]
176
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
small amount of auxiliary load. The cable connecting two adjacent wind turbine
transformers is 505 m long.
Figure 6-1 Layout of Nysted offshore wind farm collection grid and its connection with onshore
main grid
Three-phase voltages measured at P1, P2 and P3 are shown in Figure 6-2 (a), Figure 6-2
(b) and Figure 6-2 (c), respectively. In the figures, the base value of the measured
voltages is the phase-to-ground voltage (which is 33/√3 kV). As can be seen, right after
the energisation, there were high frequency oscillations of the voltages around the
energisation instants, which provoked the occurrence of overvoltages. The overvoltages
observed at P3 were slightly higher than those observed at the other two measurement
177
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
points. The high frequency oscillation of three-phase voltages decayed within about 1
ms. After the decay, certain level of voltage distortion can still be seen.
Figure 6-3 Voltage waveforms around the energisation instants (at location P1) [128]
Three-phase currents were also measured at P1, P2 and P3, which are shown in Figure
6-4 (a), Figure 6-4 (b) and Figure 6-4 (c), respectively. In the figures, the base value for
the currents measured at P1 is 420 A; for the currents measured at P2 and P3, the base
178
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
The first peak magnitudes of three-phase inrush currents are illustrated in Table 6-2. It
can be deduced that the inrush current peak magnitude (1191 A) measured at P1 was
close to nine times the peak magnitude of the inrush current drawn by transformer A1.
After the first inrush current peak, the inrush currents started to decay. In Figure 6-5, the
peak magnitudes of the first ten cycles of the phase A inrush current measured at P1 are
plotted, which shows that it decayed by about 60% (from 1191 A to 446 A) within 180
179
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
ms. For the same time span (180 ms), phase B inrush current decayed by 68% and phase
C by 62.5%.
Table 6-2 Inrush current first peaks resulted from energisation of feeder-A
1400
Peaks of phase A inrush current (A)
1200
1000
800
600
400
200
0
0 25 50 75 100 125 150 175 200
Time (ms)
Wind farm collection grid models developed in previous research, including those
presented in [36], [75], [102], [103] and [50], showed that: the auxiliary load on the
secondary side of the wind turbine transformer can be neglected or simply represented
by a constant resistance; the nominal pi model can be used to represent cables; an ideal
voltage source connected with an equivalent source impedance can be used to represent
the external grid. These modelling experiences together with those obtained in Chapter
3 were used to develop the collection grid model of the Nysted offshore wind farm
using ATP/EMTP.
The developed network model is shown in Figure 6-6. As can be seen, it consists of four
180
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
building blocks. Block 1 models the connection with the onshore main grid, which was
represented by an ideal voltage source connected in series with an equivalent source
impedance; the shunt reactor was modelled by linear inductor. In block 2, the 132 kV
onshore cable and offshore cable were represented by a number of pi sections. In block
3, the Hybrid Transformer model was used to represent the offshore platform
transformer. Block 4 is the main part of the wind farm collection grid model and it is
mainly formed by circuit breakers, cable sections and wind turbine transformers: the
circuit breakers were modelled by ideal time-controlled switch; the 33 kV cable sections
were represented by nominal pi model; the wind turbine transformers were modelled by
BCTRAN+. Further details of the building blocks are described in the following
subsections. This network model can be expanded to include more feeders.
Figure 6-6 Complete network model of the Nysted wind farm collection grid and its connection
with the main grid
181
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
The data used to define the geometry and material properties of the 132 kV single-core
onshore cable, 132 kV and 33 kV three-core submarine cables are given in Table 6-3
and Table 6-4, respectively (note that cables between A9 and A3 is of 95 mm2 cross
sectional area and cables between A3 and the offshore platform is of 185 mm2 cross
sectional area, which follows the design used in Lillgrund offshore wind farm [131]).
Table 6-3 132 kV single core onshore cables [132]
Cross-section of conductor [mm2] 1200
Laying depth [m] 1
Flat formation spacing [m] 0.25
Resistivity [Ohm.m] Al 3.6E-8
Radius of conductor [m] 0.0214
Thickness of conductor shielding [m] 0.0018
Thickness of insulation [m] 0.015
Thickness of insulation shielding [m] 0.0018
Thickness of lead sheath [m] 0.003
Thickness of anti-corrosion sheath [m] 0.004
Total radius [m] 0.04685
182
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
As mentioned before, the BCTRAN model can be derived based on transformer test
reports (mainly using short circuit and open circuit test data). The information obtained
from manufacturer for modelling the wind turbine transformer using BCTRAN is given
in Table 6-5. The core saturation effects were modelled by three wye-connected type-96
hysteretic inductors whose saturation characteristic is shown in Figure 6-7. The air-core
inductance of the transformer was assumed to be twice the transformer short-circuit
inductance.
Table 6-5 Main electrical information for modelling wind turbine transformers [128]
Parameters Value
Rated power [MVA] 2.5
Connection group Dyn7
Rated primary voltage [kV] 33
Rated secondary voltage [kV] 0.69
Copper losses [kW] 22
No-load losses [kW] 5.5
Short-circuit impedance [%] 8.3
183
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
Table 6-6 Network model parameter settings for simulating field measurement results
Parameters Setting Note
Source voltage level 107.7 kV Peak value of 132 kV (L-L).
Source impedance Rp=0.3 Ohm, Xp=3.454 Ohm Derived from 5 GW fault level
Ph. A 3.55 ms Closing time reference: positive-going
Circuit breaker closing
Ph. B 3.33 ms zero crossing of phase B voltage;
time
Ph. C 3.72 ms Only applied to feeder circuit breaker.
Ph. A 0
Applied to all wind turbine transformers;
Residual flux Ph. B +0.306 pu
Base value is peak nominal flux.
Ph. C –0.306 pu
In Figure 6-8, the simulated voltage waveforms at location P1 are compared with those
obtained from the measurements. The high frequency voltage oscillation appearing
around the energisation instant is not well replicated, due to the frequency response
range of the network model components is up to 1 kHz. Nevertheless, the simulated
voltages agree reasonably well with the measured voltages in terms of circuit breaker
closing times, peak magnitudes and waveform patterns.
Energisation Measured voltages
Va Vb Vc
The simulated three-phase inrush currents drawn by wind turbine transformer A9 are
compared with measurements in Figure 6-9. In addition, the total inrush currents drawn
184
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
by all the wind turbine transformers located at feeder-A are compared with
measurements in Figure 6-10. As can be seen from both comparisons, the peak
magnitudes and patterns of the simulated three-phase inrush currents are very similar to
those of the measured inrush currents. This indicates that the network model is not only
capable to estimate the inrush currents drawn by a stand-alone wind turbine transformer
but also the total inrush currents drawn by an array of wind turbine transformers.
Figure 6-9 Comparison between measurement and simulation regarding the inrush currents
drawn by wind turbine transformer A9
Figure 6-10 Comparison between measurement and simulation regarding the total inrush
currents drawn by nine wind turbine transformers in feeder-A
The inrush current decay trends were also compared between the simulated waveforms
and measured ones. The one shown in Figure 6-11 compares the inrush currents drawn
by all the wind turbine transformers connected at feeder-A. In Figure 6-12, two decay
trend comparisons are shown, one compares the phase A inrush current drawn by wind
turbine transformer A1 (in red curves) and the other compares the phase A inrush
current drawn by wind turbine transformer A9 (in blue curves). Each trend line indicates
the decay of inrush current peak magnitudes during the first twelve cycles. As can be
185
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
seen, the comparisons show that the decay trends of simulated inrush currents match
well with those measured decay trend lines. This indicates that the system parameters
that determining the decay time constants of the transformer inrush transients are
correctly represented by the network model.
Specifically, the comparisons in Figure 6-12 also shows the faster decay of inrush
currents in wind farm transformer A9 than that in A1 and this is also well replicated in
the simulation results, which indicates the losses associated with the cable sections
between A1 and A9 are accurately represented. This is important for studying
sympathetic inrush between wind turbine transformers, as the amount of cable section
losses would significantly influence sympathetic interaction between wind turbine
transformers.
The above validation comparisons indicate that the network model is capable to
reproduce field measurement results in terms of inrush current patterns, peak
magnitudes and decay trends, which confirms the accuracy of the developed network
model.
Figure 6-12 Decay trend comparisons regarding the inrush current drawn by transformer A1 and
that drawn by A9
186
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
187
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
under the weak source strength, up to three wind turbine transformers are allowed to be
energised together.
Figure 6-13 Case studies of energising wind turbine transformers in one feeder
Table 6-7 Voltage dips resulted from energising a feeder of wind turbine transformer under the
strong source strength
Voltage dips
Phase A Phase B Phase C
Case
Magnitude Duration Magnitude Duration Magnitude Duration
(%) (ms) (%) (ms) (%) (ms)
W1 0.55 0 1.24 0 1.2 0
W2 0.3 0 0.81 0 0.75 0
W3 0.3 0 0.75 0 0.72 0
W4 0.2 0 0.52 0 0.47 0
W5 0.13 0 0.36 0 0.32 0
W6 0.06 0 0.18 0 0.16 0
Table 6-8 Voltage dips resulted from energising a feeder of wind turbine transformer under the
weak source strength
Voltage dips
Phase A Phase B Phase C
Case
Magnitude Duration Magnitude Duration Magnitude Duration
(%) (ms) (%) (ms) (%) (ms)
W1 2.9 0 6.6 80 5.8 65
W2 1.8 0 4.4 51 3.7 30
W3 1.7 0 4.2 32 3.7 20
W4 1.1 0 2.89 0 2.4 0
W5 0.75 0 2.0 0 1.65 0
W6 0.35 0 1.0 0 0.85 0
188
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
breaker closing time was defined by one cycle Uniform distributed common-order-time
and a Gaussian distributed closing offset time (mean value equal to 0 and standard
deviation equal to 0.833 ms); the range of the residual flux was defined as ±0.8 pu and
of Uniform distribution (note that the residual flux in each wind turbine transformer is
independently assigned). The results obtained from 1000 stochastic runs are shown in
Figure 6-14 for the frequency of occurrence of dip magnitude and Figure 6-15 regarding
the frequency of occurrence of dip duration. The results indicate that: the frequency of
reaching the worst case voltage dip is actually negligible; over 95% of the dips are with
magnitudes less than 0.6 pu of the worst case.
Figure 6-14 Frequency of voltage dip magnitude in three phases at the point-of-common-
coupling under energising a feeder of wind turbine transformers
Figure 6-15 Frequency of voltage dip duration in three phases at the point-of-common-coupling
under energising a feeder of wind turbine transformers
189
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
As can be seen, on the secondary side of 132 kV transformer with YNyn0 connection,
the observed three-phase voltage dips are identical with those observed at the point-of-
common-coupling (this finding is similar to the voltage propagation through the
400/132 kV autotransformer shown in Figure 4-14); for the voltage dips observed on the
secondary side of 132 kV transformer with YNd1 and YNd11 connections, the biggest
voltage dip has been increased by 12.5%, comparing with those observed at the point-
of-common-coupling; in addition, the transformer with YNd1 connection changes the
phase with the biggest voltage dip from phase B to phase C.
The above analysis indicates that grid transformer with YNyn0 connection has no
impact on the propagation of three-phase voltage dips; the grid transformer with YNd1
or YNd11 connection would increase the dip magnitude seen by the end users.
Observation points
Figure 6-16 Effects of transformer winding connections on voltage dip propagation
190
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
to ER-P28, see Figure 2-21); for examples, the limit is about 1.3% if the interval is 1
minute; about 0.4% if the interval is 1 second.
Supposing the time interval between two consecutive energisations is 1 second, voltage
dips caused by Case W4, W5 and W6 were further tested (Case W2 and W3 were not
tested because the voltage dip magnitudes resulted from these cases already exceeded
the grid code limit). For all the cases, the energisation was conducted under the
commonly agreed worst case energisation condition against the weak source strength
and the energisation direction was from the end closest to the offshore platform to the
end farthest from the offshore platform.
Results are shown in Figure 6-17, Figure 6-18 and Figure 6-19, respectively. As can be
seen, the maximum dip magnitudes resulted from consecutive energisation under Case
W4, W5 and W6 are 2.95%, 2% and 1.2%, respectively. These dip magnitudes all
exceed the 0.4% limit stated in the grid code. This indicates that, under the assumed
energisation condition, the time interval between two consecutive energisations cannot
be as short as 1 second. By comparing the maximum dip magnitudes with those limits
given in Figure 2-21, it can be deduced that: to comply with the dip limits, the interval
for Case W4 should be more than 12.5 minutes; for Case W5, it should be more than 3.3
minutes; for Case W6, the time interval is the shortest, which is about 1 minute. This
shows that consecutively switching in only one transformer at a time is the
recommended way to energise a feeder of wind turbine transformers while keeping the
lowest risk of exceeding the grid code limit.
191
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
In Figure 6-20, the configuration of two feeders (feeder-A and feeder-B) connected to
offshore platform is schematically shown: feeder-A is to be energised and feeder-B is
already connected. The distance of the cable connection between transformer A1 and
B1 is defined as the electrical distance between the two feeders. With different electrical
distances, Case W1 energisation of feeder-A transformers were conducted to assess the
possible sympathetic inrush induced on feeder-B. Figure 6-21, Figure 6-22 and Figure
6-23 show the results corresponding to 1, 2 and 3 km electrical distance, respectively.
192
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
Figure 6-20 Schematic diagram of two wind farm feeders connected at offshore platform
As can be seen, in the case with 1 km electrical distance, the maximum instantaneous
peak of the total sympathetic inrush current drawn by feeder-B transformers is about
130 A and the duration of the sympathetic inrush is less than 0.9 second. Compared to
the sympathetic inrush currents drawn by large GSU transformers, the degree of
sympathetic inrush is much less significant. As the electrical distance between two
adjacent feeders increases, both the peak magnitude and the duration of sympathetic
inrush currents decrease considerably, because the submarine cables connected between
the two feeders provide resistive damping that suppresses the sympathetic inrush.
Figure 6-21 Sympathetic inrush currents drawn by the already connected feeder (1 km electrical
distance between two feeders)
Figure 6-22 Sympathetic inrush currents drawn by the already connected feeder (2 km electrical
distance between two feeders)
Figure 6-23 Sympathetic inrush currents drawn by the already connected feeder (3 km electrical
distance between two feeders)
193
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
Furthermore, similar studies were carried out for other energisation cases involving
energisation of multiple transformers. For Case W2 and W3 (with 1 km electrical
distance), it was found that the maximum peak of the total sympathetic inrush currents
is about 75 A and the duration of the sympathetic inrush is about 0.75 second. For Case
W4 and W5 (with 1 km electrical distance), the sympathetic inrush current is negligible,
which suggests that, in the present offshore wind farm grid, sympathetic interaction
between two feeders could be of little concern if no more than three wind turbine
transformers are being energised together.
When carrying out energisation of multiple transformers, the residual fluxes in wind
turbine transformers being energised may differ from each other. Figure 6-24 shows an
extreme condition that: the four transformers near to the offshore platform (A1-A4)
possess maximum residual flux opposite to the flux build up; the other five wind turbine
transformers (A5-A9) possess maximum residual flux in line with the flux build-up.
Figure 6-24 Wind turbine transformers with different residual flux condition
Simultaneous energisation of these nine wind turbine transformers was simulated under
the worst energisation condition against the weak source strength. The currents
observed at the primary sides of these wind turbine transformers are shown in Figure
6-25. Although the transformers were energised simultaneously, inrush took place in the
wind turbine transformer A5, A6, A7, A8 and A9, whilst sympathetic inrush took place
in the wind turbine transformer A1, A2, A3 and A4. This indicates that: sympathetic
interaction could even be induced between the energised transformers if their residual
flux conditions are different; the transformers which are relatively less saturated would
be forced to engage sympathetic inrush by the relatively more saturated transformers.
194
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
Figure 6-25 Sympathetic and inrush currents in the wind turbine transformers being energised
together
Under the Case W6_1 energisation, the resulted sympathetic inrush currents at the
already connected wind turbine transformers are shown in Figure 6-26. As can be seen,
sympathetic inrush currents are induced in A8, A7, A6, A5 and A4. A8 exhibits the
largest sympathetic inrush peak which is about 25 A. From A7 to A4, the magnitudes of
maximum peaks decrease from 24.5 to 4.6 A. Regarding the duration of sympathetic
195
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
inrush, the longest one is also seen at A8 and the shortest one is seen at A4. The
sympathetic inrush currents in transformers A1, A2 and A3 are of negligible level. Such
a distribution of sympathetic inrush current could be attributed to following reasons:
• When the high inrush currents drawn by A9 flow through the feeder cables, there is
voltage distortion being built up by the resistive elements of the cable sections,
especially on those 33 kV cable sections because of their higher resistivity; since A8
is the closest to A9, the resistive voltage distortion seen by A8 is the most significant
and therefore the largest sympathetic inrush is built up in A8;
• From A8 further away down to A4, smaller and smaller sympathetic inrush currents
could be induced, due to the increasing electrical distance from A9 causing bigger
losses and shorter feeder cable for building up resistive voltage distortion;
• The established sympathetic inrush currents would balance out the distortion caused
by inrush currents of A9 so that very little offset flux can be seen in A1, A2 and A3,
therefore they were of negligible sympathetic inrush footprints.
196
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
In Case W6_2, A1 is the wind turbine transformer being energised with other wind
turbine transformers already connected. The simulation results of this case are
illustrated in Figure 6-27. Following the energisation, the peaks of sympathetic inrush
currents induced in the already connected wind turbine transformers are relatively the
same. This is because the wind turbine transformers experience almost the same amount
of resistive voltage distortion built up on the cables connecting between the supply
source and the wind turbine transformer A1. The duration of the sympathetic inrush
currents, however, exhibits considerable differences. The duration of sympathetic inrush
in A2 is the largest, for it is located closest to A1; as the electrical distance between A1
and other already connected wind turbine transformers increases, the duration of
sympathetic inrush decreases. However, in contrasting to the results given by Case
W6_1, the initiation process of all the sympathetic inrush currents is much slower and
the maximum peaks are much smaller, which suggests that the resistive voltage
distortion across the system between supply source and A1 is much smaller than that in
Case W6_1.
197
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
From the above two case studies, it can be deduced that the degree of sympathetic
inrush induced on each wind turbine transformer is largely related to the location of the
wind turbine transformer being energised and the relative location of the other already
connected wind turbine transformers.
198
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
199
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
No Complete all
energization?
Yes
Figure 6-29 Procedure to correlate energisation sequence and sympathetic inrush level
200
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
deterministically assessed. In the assessment, the energisation condition for all the
energisation events was initialized to the commonly agreed worst case energisation
condition (same as in previous studies, i.e., three phases were simultaneously energised
at the positive-going zero crossing of phase C line-to-ground voltage; residual flux of
phase A, B and C was of -0.8 pu, 0 and +0.8 pu of peak nominal flux, respectively).
The simulation results given by assessing sequences S1, S2, S3 and S4 are shown in
Figure 6-30, Figure 6-31, Figure 6-32 and Figure 6-33, respectively. In all the figures,
the wind turbine transformer is indicated by the number shown on the horizontal axis;
each column group represents the sympathetic inrush level of three phases accumulated
on each wind turbine transformer under a particular energisation sequence.
As can be seen, for any column group, the sympathetic inrush level accumulated on
phase C is the highest among the three phases (this is due to the specific energisation
condition). The following analysis focuses on the sympathetic inrush level of phase C
only. Regarding the maximum accumulated sympathetic inrush level, the one in S1 is
about 130 A·s; in S2, the maximum level is about 295 A·s, which is almost identical to
that appear in S3 and is more than twice in S1; in S4, the maximum level is about 205
A·s. Comparing the severity of sympathetic inrush associated with each wind turbine
transformer, A8 and A9 would experience the most intensive sympathetic inrush if S2
or S3 is followed, whilst A5 and A6 would be more vulnerable if S1 or S4 is applied. It
can also be seen that the wind turbine transformers located close to the platform,
including A1, A2 and A3, experience rather small level of sympathetic inrush under all
the energisation sequences. In addition, aggregation of all the accumulated sympathetic
inrush levels under each energisation sequence was made (in total, four aggregations)
and they are compared in Table 6-10. It can be seen that the aggregation of the
accumulated sympathetic inrush levels resulted from S1 is the lowest, followed by S3,
S4 and S2.
The above results indicate that, when a feeder of wind turbine transformers are to be
sequentially energised, to result in less sympathetic inrush, wind turbine transformers
should be energised from the one closest to the offshore platform towards the one
farthest from the offshore platform, because it gives the smallest level of sympathetic
inrush on a stand-alone wind turbine transformer and also the least aggregated
sympathetic inrush effects compared to other energisation sequences.
201
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
Table 6-10 Aggregation of sympathetic inrush levels resulted from each energisation sequence
Energisation sequence S1 S2 S3 S4
Aggregation of sympathetic inrush level (A·s)
Accumulated sympathetic inrush level, A·s 1434 1796 1581 1590
A1 A2 A3 A4 A5 A6 A7 A8 A9
Wind turbine transformer
Figure 6-30 Accumulated sympathetic inrush level on each wind turbine transformer resulted
from deterministic testing of S1
Accumulated sympathetic inrush level, A·s
A1 A2 A3 A4 A5 A6 A7 A8 A9
Wind turbine transformer
Figure 6-31 Accumulated sympathetic inrush level on each wind turbine transformer resulted
from deterministic testing of S2
202
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
A1 A2 A3 A4 A5 A6 A7 A8 A9
Wind turbine transformer
Figure 6-32 Accumulated sympathetic inrush level on each wind turbine transformer resulted
from deterministic testing of S3
Accumulated sympathetic inrush level, A·s
A1 A2 A3 A4 A5 A6 A7 A8 A9
Wind turbine transformer
Figure 6-33 Accumulated sympathetic inrush level on each wind turbine transformer resulted
from deterministic testing of S4
The energisation condition was stochastically initialized as follows: the closing time of
each circuit breaker pole was defined by a common-order-time (same for three poles,
uniformly distributed over one power frequency cycle) plus a Gaussian distributed
closing offset time (with the mean value equals to zero and standard deviation equals to
0.833 ms) [43]; the range of the residual flux was defined as between -0.8 pu and +0.8
pu of the peak nominal flux and it was characterized by Uniform distribution (note that
the residual flux in each wind turbine transformer was independently assigned). For
each energisation sequence, 1000 stochastic simulation tests were carried out.
By using box plot, results obtained from stochastic testing energisation sequence S1, S2,
S3 and S4 are presented in Figure 6-34, Figure 6-35, Figure 6-36 and Figure 6-37,
respectively. In these figures, each column represents 1000 sympathetic inrush levels
accumulated on phase C of a specific wind turbine transformer (the results of other two
phases are similar to that of phase C); the central mark on each box is the median; the
lower and upper edges of the box are the 25th and 75th of the percentiles; the length of
the whisker is defined by w whose value is 1.5 (this value corresponds to 99.3%
coverage if the data are normally distributed); points are displayed using + as outliers if
they are larger than Q75+w·(Q75–Q25), where Q25 and Q75 are the 25th and 75th
percentiles, respectively.
Accumulated sympathetic inrush level, A·s
whisker
outliners
75 %
25 %
median
A1 A2 A3 A4 A5 A6 A7 A8 A9
Wind turbine transformer
Figure 6-34 Accumulated sympathetic inrush level of each wind turbine transformer resulted
from stochastic testing of S1
204
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
A1 A2 A3 A4 A5 A6 A7 A8 A9
Wind turbine transformer
Figure 6-35 Accumulated sympathetic inrush level of each wind turbine transformer resulted
from stochastic testing of S2
Accumulated sympathetic inrush level, A·s
A1 A2 A3 A4 A5 A6 A7 A8 A9
Wind turbine transformer
Figure 6-36 Accumulated sympathetic inrush level of each wind turbine transformer resulted
from stochastic testing of S3
Accumulated sympathetic inrush level, A·s
A1 A2 A3 A4 A5 A6 A7 A8 A9
Wind turbine transformer
Figure 6-37 Accumulated sympathetic inrush level of each wind turbine transformer resulted
from stochastic testing of S4
205
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
As can be seen from Figure 6-34 to Figure 6-37, the accumulated sympathetic inrush
level for each wind turbine transformer is scattered, due to the presence of stochastic
energisation conditions. The scatter ranges are relatively small for transformers A1, A2,
A3 and A4 and are relatively large for A6, A7, A8, which are evidenced in all
energisation sequences. The scatter range for A5 is relatively large in S1, S2 and S4,
while the scatter range for A9 is relatively large in the case of S2 and S4. The relatively
large scatter ranges for A5, A6, A7, A8 and A9 indicate that wind turbine transformers
located at these positions are likely to be affected by sympathetic inrush, which is also
evidenced by the findings obtained in previous deterministic studies.
As far as the median of each column is concerned, the lowest one is still found to be in
S1 and the profile of accumulated sympathetic inrush level formed by nine wind turbine
transformers is similar to that showed in Figure 6-30, which suggests that S1 is the
energisation sequence which would induce minimum sympathetic inrush between wind
turbine transformers. Therefore, the findings gained from the deterministic assessment
are further validated by the stochastic studies.
6.8 Summary
In this chapter, a network model of a wind farm collection grid was developed and
validated against field measurement results obtained in literature, which was then used
to study voltage dips and sympathetic inrush caused by energising wind turbine
transformers.
Regarding voltage dips, the above case studies show that, in the studied system,
energising one wind turbine transformer against the weak source strength can only
result in dip magnitude of no more than 1%, therefore causes no concern on breaching
grid code requirements. Concerns may be raised in the cases of simultaneously
energising multiple wind turbine transformers, because the resulted voltage dips may
reach 6.6%, even though stochastic estimation shows that the probability of reaching
such a dip magnitude is very low. The winding connection of 132/33 kV transformers
should be carefully considered, as it may change the phase with the biggest voltage dip
magnitude and result in larger dip magnitude seen by the end users. In addition, care
should be taken if consecutive energisation of wind turbine transformers (with time
206
Chapter 6 Assessment of Transformer Energisation Transients Due to Offshore Wind Farm Connection
interval shorter than 750 seconds) is carried out, because the voltage dip limit is further
tightened.
Regarding the sympathetic inrush caused by energising wind turbine transformers, the
performed studies consist of two parts: one is to on energisation of multiple wind
turbine transformers and the other is on energisation of a stand-alone wind turbine
transformer.
In the case of energisation of multiple transformers, it was found that: the degree of
sympathetic interaction between two adjacent feeders is rather minor, because the cable
connection between the feeders can contribute significant resistive losses; one should be
cautious that sympathetic inrush may occur between the transformers being energised
together, if they have a different residual flux.
In the case of energisation of a stand-alone transformer, it was found that the degree of
sympathetic inrush is largely related to the location of the wind turbine transformer
being energised and the relative location of the other already connected transformers.
Furthermore, the potential relationships between energisation sequence and sympathetic
inrush level were deterministically and stochastically evaluated, suggesting that the
energisation sequence that would result in less sympathetic inrush level between wind
turbine transformers is to energise wind turbine transformers from the one closest to the
offshore platform to the one farthest from the offshore platform.
207
Chapter 7 Conclusion and Future Work
209
Chapter 7 Conclusion and Future Work
Through analysing the field measurements which were carried out in the 400 kV
transmission grid, it was shown that the energisation of GSU transformers can trigger a
network-wide voltage dips (i.e., voltage dips not only appear at the substation connected
by the transformer being energised but also at other substations in the network); the
recorded maximum RMS voltage dip was about 7.85%; the measured sympathetic
inrush currents showed that the duration of sympathetic inrush lasted for more than 20
seconds and so did the full recovery of the resulted voltage dips.
Based on the system parameters provided by the network operator, a 400 kV grid
network model was developed in ATP/EMTP by following the modelling guidelines
summarized through the literature review. The successful validation of the network
model confirms that: the source network can be modelled by an ideal sine-wave source
and a Thevenin equivalent impedance of the part of the network not under study; the
transmission network between the supply source and the energised transformer should
be represented in detail, taking into account the transmission lines, system loading and
reactive power compensation devices; the constant frequency line model can be used to
210
Chapter 7 Conclusion and Future Work
represent transmission lines, with line dimension and transposing scheme considered;
system loading can be represented by lumped constant resistance and inductance
connected in parallel; transformers can be modelled by the BCTRAN routine with
transformer core externally represented by three delta-connected type-96 non-linear
inductors.
The validated network model was utilized to carry out a comprehensive study on
voltage dips using deterministic approach. Through assessing the degrees of voltage
dips under different energisation conditions, it was found that: in the present system
with two GSU transformers simultaneously energised under the worst energisation
condition, the maximum voltage dip, observed at the substation closest to the
transformers being energised, would be of magnitude about 18% and duration about 3.5
seconds; with the presence of sympathetic interaction, the dip duration can prolonged by
125% (increased from 3.5 seconds to 7.9 seconds); the voltage dips propagating to 132
kV side can be of longer dip duration due to the sympathetic inrush of substation
transformers.
Furthermore, assessing the network-wide voltage dips for the complete network under
non-outage condition suggests that: the dip magnitudes observed at each substation are
related to the distance between the substation and the supply source and also the
distance between the substation and the energised transformers; those substations
located in the proximity of the energised transformer and relatively far away from
supply source are subjected to larger dip magnitudes.
With line outage taken into account, it was found that: the network-wide voltage dip
outcome under single-circuit outage situation is similar to that observed under non-
outage condition; however, if there is double-circuit outage resulting in significant
network topology change, both dip magnitude and duration can be significantly
exacerbated; in the system under study, the most unfavorable double-circuit outage can
increase the dip magnitude from 18% to about 30%.
The sensitivity assessment shows that: transformer core saturation inductance has the
most profound impact on the voltage dip magnitude; the amount of transformer copper
losses is the most influential parameter on determining the voltage dip duration.
211
Chapter 7 Conclusion and Future Work
Possible operational measures to control the voltage dips were found to be adjusting
GSU tap changer to maximum tap, opening the coupler circuit breaker and applying
SVC. It was found that, if these operational measures are applied simultaneously, the
dip magnitude and duration resulted from worst case energisation can be reduced by
37% and 85%, respectively.
Monte Carlo simulation was conducted to extend the few deterministically-defined case
studies to many stochastically-defined case studies. A dip frequency pattern was
identified based on over 1000 stochastic runs and it was found to be sensitive to the
distribution of residual flux but insensitive to the distribution of closing offset time. In
addition, it was shown that the dip frequency pattern is insensitive to the system
condition variation in a range of ±25%. Furthermore, it was found that the probability of
reaching the worst case dip magnitude (estimated by the commonly agreed worst case
energisation condition) is lower than 0.5%, indicating that the worst case scenario is
unlikely to occur; in fact, about 80% of the dips are with magnitudes lower than 0.6 pu
of the worst case dip magnitude. Nevertheless, it was shown that there exist dips with
magnitudes exceeding the worst case dip magnitude, indicating the inadequacy of
deterministic assessment approach by using the commonly agreed worst case
energisation condition.
Voltage dips caused by energising wind turbine transformers were studied and
compared with those caused by energising large GSU transformers. It was found that, in
the present system under study, energising a stand-alone wind turbine transformer
against the weak source strength can only result in dip magnitude of no more than 1 %
(as the fault level of the source network is more than three hundred times larger than the
rating of the wind turbine transformer), therefore causes no concern on complying grid
code requirements. Concerns may be raised in the case of simultaneously energising a
feeder of wind turbine transformers or consecutive energisation of transformers with
short time interval.
212
Chapter 7 Conclusion and Future Work
significant resistive losses), one should be cautious that sympathetic inrush may occur
between the transformers being energised together, if they are of different residual flux.
The potential relationships between sympathetic inrush and energisation sequence were
both deterministically and stochastically evaluated and the main conclusion reached is
that the optimum energisation sequence to achieve less sympathetic inrush between
wind turbine transformers is to energise wind turbine transformers from the one closest
to the offshore platform to the one farthest away from the offshore platform.
2) The use of BCTRAN+ model in this thesis does not take into account the
influence of transformer core topology. The effect of core topology may be taken
into account by using Hybrid Transformer model to represent GSU transformers,
if more transformer design information is available. Then it would be interesting
to see if Hybrid Model can more accurately simulate sympathetic inrush
between GSU transformers.
1) In the deterministic case studies, it was assumed that the generators on the
secondary side of the GSU transformers were not connected yet. It is possible
that one of the generators is already connected while carrying out the
energisation of a GSU transformer. In this case, it would be interesting to
evaluate the influence of the generator’s AVR responses on the outcome of
sympathetic inrush.
213
Chapter 7 Conclusion and Future Work
For assessment of energising wind turbine transformer into offshore wind farm grid:
3) The degree of sympathetic inrush between wind turbine transformers might vary
with the layout of the wind farm collection grid. The studies carried out in this
thesis are focusing on radial layout only. Further studies can be carried out for
guiding the energisation of offshore wind farms with other layouts.
214
References
References
[1] L. F. Blume, G. Camilli, S. B. Farnham, and H. A. Peterson, "Transformer
magnetizing inrush currents and influence on system operation," Transactions of
the American Institute of Electrical Engineers, vol. 63, no. 6, pp. 366-375, 1944.
[2] I. Arana, A. Hernandez, G. Thumm, and J. Holboell, "Energization of wind
turbine transformers with an auxiliary generator in a large offshore wind farm
during islanded operation," IEEE Transactions on Power Delivery, vol. 26, no. 4,
pp. 2792-2800, 2011.
[3] N. Chiesa and H. K. Høidalen, "Novel approach for reducing transformer inrush
currents: laboratory measurements, analytical interpretation and simulation
studies," IEEE Transactions on Power Delivery, vol. 25, no. 4, pp. 2609-2616
2011.
[4] Engineering-Recommendation-P28, "Planning limits for voltage fluctuations
caused by industrial, commercial and domestic equipment in the UK," 1989.
[5] "National electricity transmission system security and quality of supply
standard," National Grid, Mar. 2012.
[6] L. Prikler, G. Bánfai, G. Bán, and P. Becker, "Reducing the magnetizing inrush
current by means of controlled energization and de-energization of large power
transformers," Electric Power Systems Research, vol. 76, no. 8, pp. 642-649,
2006.
[7] S. P. Patel, "Fundamentals of transformer inrush," in 64th Annual Conference
for Protective Relay Engineers, pp. 290-300, Texas, USA, Apr. 2011.
[8] M. Steurer and K. Frohlich, "The impact of inrush currents on the mechanical
stress of high voltage power transformer coils," IEEE Transactions on Power
Delivery, vol. 17, no. 1, pp. 155-160, 2002.
[9] J. Faiz, B. M. Ebrahimi, and T. Noori, "Three- and two-dimensional finite-
element computation of inrush current and short-circuit electromagnetic forces
on windings of a three-phase core-type power transformer," IEEE Transactions
on Magnetics, vol. 44, no. 5, pp. 590-597, 2008.
[10] R. Yacamini and A. Abu-Nasser, "Numerical calculation of inrush current in
single-phase transformers," IEE Proceedings B: Electric Power Applications, vol.
128, no. 6, pp. 327-334, 1981.
[11] R. A. Turner and K. S. Smith, "Transformer inrush currents," IEEE Industry
Applications Magazine, vol. 16, no. 5, pp. 14-19, 2010.
[12] A. Ketabi, A. M. Ranjbar, and R. Feuillet, "Analysis and control of temporary
overvoltages for automated restoration planning," IEEE Transactions on Power
Delivery, vol. 17, no. 4, pp. 1121-1127, 2002.
[13] R. Yacamini and A. Abu-Nasser, "Transformer inrush currents and their
associated overvoltages in HVDC schemes," IEE Proceedings C: Generation,
Transmission and Distribution, vol. 133, no. 6, pp. 353-358, 1986.
[14] R. A. Turner and K. S. Smith, "Resonance excited by transformer inrush current
in inter-connected offshore power systems," in IEEE Industry Applications
Society Annual Meeting, pp. 1-7, Edmonton, Canada, Oct. 2008.
215
References
216
References
217
References
219
References
220
References
221
References
[113] "Electrical parameters and impedance characteristics of plant, lines and cables,"
National Grid technical guidance note, 2002.
[114] M. Bolhasani, S. S. H. Kamangar, and S. Tavakoli, "Determination of
distribution function of inrush current in three phase transformer using Monte
Carlo method," in 47th International Universities Power Engineering Conference
(UPEC), pp. 1-7, Sep. 2012.
[115] J. A. Martinez and J. Martin-Arnedo, "Voltage sag stochastic prediction using an
electromagnetic transients program," IEEE Transactions on Power Delivery, vol.
19, no. 4, pp. 1975-1982, 2004.
[116] A. Greenwood, "Electrical transients in power systems," 2nd ed, New York:
Wiley, 1991.
[117] N. Chiesa, A. Avendano, H. K. Høidalen, B. A. Mork, D. Ishchenko, and A. P.
Kunze, "On the ringdown transient of transformers," in International Conference
on Power Systems Transients (IPST), Lyon, France, Jun. 2007.
[118] J. H. Brunke and K. J. Frohlich, "Elimination of transformer inrush currents by
controlled switching. I. Theoretical considerations," IEEE Transactions on
Power Delivery, vol. 16, no. 2, pp. 276-280, 2001.
[119] "UK renewable energy roadmap," available from
http://www.decc.gov.uk/assets/decc/11/meeting-energy-demand/renewable-
energy/2167-uk-renewable-energy-roadmap.pdf, accessed date: 25th Jan, 2013.
[120] "Germany's wind energy industry," available from
http://www.gtai.de/GTAI/Content/EN/Invest/_SharedDocs/Downloads/GTAI/Fa
ct-sheets/Energy-environmental/fact-sheet-wind-energy-in-germany.pdf,
accessed date: 25th Jan, 2013.
[121] "Pushing forward: the future of offshore wind energy," available from
http://www.rolandberger.de/media/pdf/Roland_Berger_taC_Offshore_20111206
.pdf, accessed date: 25th Jan, 2013.
[122] G. Quinonez-Varela, G. W. Ault, O. Anaya-Lara, and J. R. McDonald,
"Electrical collector system options for large offshore wind farms," IET
Renewable Power Generation, vol. 1, no. 2, pp. 107-114, 2007.
[123] "SeaTitan 10 MW wind turbine," available from http://www.amsc.com/library/,
accessed date: 25th Jan, 2013.
[124] A. B. J. Green, L. J. Fingersh and Y. Wan, "Electrical collection and
transmission systems for offshore wind power," in Offshore Technology
Conference, Houston, Texas, 2007.
[125] "List of offshore wind farms," available from
http://en.wikipedia.org/wiki/List_of_offshore_wind_farms, accessed date: 25th
Jan, 2013.
[126] M. Rioual and M. Sow, "Study of the sympathetic interactions when energizing
transformers for wind-farms: Description of the phenomena involved and
determination of the stresses during their energization," in 7th International
Workshop on Large-Scale Integration of Wind Power into Power Systems,
Madrid, 2008.
223
References
[127] IEC-60076-16, "Power transformers - part 16: Transformers for wind turbine
applications," 2011.
[128] I. Arana, "Modelling of switching transients in Nysted offshore wind farm and a
comparison with measurements," MSc Thesis, Technical University of Denmark,
2008.
[129] S. Das, N. Karnik, and S. Santoso, "Time-domain modeling of tower shadow
and wind shear in wind turbines," ISRN Renewable Energy, 2011.
[130] B. Gustavsen, J. A. Martinez, and D. Durbak, "Parameter determination for
modeling system transients-Part II: Insulated cables," IEEE Transactions on
Power Delivery, vol. 20, no. 3, pp. 2045-2050, 2005.
[131] J. Jeppsson, P. E. Larsen, and Å. Larsson, "Technical description Lillgrund wind
power plant," Vattenfall, Sep 2008.
[132] ABB, "XLPE land cable systems," 2010.
[133] ABB, "XLPE submarine cable systems," 2010.
[134] "National grid seven year statement," available from
http://www.nationalgrid.com/uk/Electricity/SYS/current/, accessed date: 15th
Nov, 2012.
[135] D. West, "Central networks grid transformer specification," E.ON Central
Networks 2008.
224
Appendix: List of Publications
225