Sizing Power Transformers Using Thermal Limits
Sizing Power Transformers Using Thermal Limits
Article history: It has already been shown in the literature that power transformers may be more accurately sized by their
- thermal limits than by their rated power limit. In practice however, thermal limits are usually considered
only in operations; but not at planning stage, where the more usual notion of rated power is used. This
Keywords: paper proposes a novel method to take into account (and benefit from) thermal limits directly at planning
Power systems planning stage. This is made possible by quantifying separately the impact of each generator and load on the
Power transformer temperature of the distribution transformer. Decoupling the effect of individual generators and loads is
Smart grid achieved by linearizing and rewriting the analytical expression of the hot-spot transformer temperature.
Dynamic thermal rating The practical value of the method is assessed using a real-world dataset, by estimating the increase in the
hosting capacity of the considered transformer for additional generators and loads. Significant gains are
obtained when the transformer is sized by generation, in particular when photovoltaic (PV) generators are
involved.
∆
Top-oil temperature °C
evolution of load and generation.
∆
Hot-spot-to-top-oil gradient °C
After loading scenarios have been defined, one of the following two
Top-oil temperature rise °C
categories of methods may be used [2]:
Parameters • The first and simpler one is to simply define physical limits, typically
Thermal model constant on instantaneous active power but possibly on other criteria, and to size
Thermal model constant the equipment so that these limits are not violated in any of the stress-
Thermal model constant tests.
• And the second, more elaborate one, is to perform a so-called “lower
x
Ratio of load losses at rated current to no-load losses
cost optimization” that aims at finding an economical trade-off between
y
Oil exponent
a smaller transformer that will be heavily loaded and undergo
Winding exponent
α
accelerated aging and increased losses, and a larger one that will cost
Transformer constant stating the relative influence of
more but not suffer from increased aging and losses [3-4], [15].
γ
and on
The second method involves estimating the loss of life of the
Benefit from using the thermal criterion %
transformer depending on its loading conditions — see for example [5]
Average oil time constant min
where the impact of high PV penetration on transformer lifetime is
∆
Winding time constant min
studied, and [6-7] where the impact of Electric Vehicles (EVs) on
∆
Hot-spot-to-top-oil gradient at rated load °C
transformer lifetime is investigated. However, predicting with great
Top-oil temperature rise at rated load °C
accuracy the useful lifetime of a transformer based on its loading
condition is currently considered difficult if not impossible [3], and as a
Distribution transformers and power (HV/MV) transformers are
consequence this paper will focus on the first method, namely, sizing
predominant assets in a distribution grid, and sizing them well is a crucial
transformers using physical limits.
task for distribution system operators (DSOs). When replacing,
___________
* Corresponding author at: Univ. Grenoble Alpes, CNRS, Grenoble INP, G2Elab, 38000 Grenoble, France
E-mail address: [Link]@[Link]
© 2019 published by Elsevier. This manuscript is made available under the CC BY NC user license
[Link]
are a current limit and a hot-spot temperature limit (2).
≤ 1.3
The simplest criterion is to never exceed the rated current (or power). It
(
is the method most commonly used by DSOs, and it is used in several
ℎ ≤ 120°$
(2)
recent research papers on the topic of distribution network planning [8-9].
≤ 1.5
used here for power transformers becomes
short duration above the threshold would be forbidden – although such a
for distribution transformers and medium power transformers.
peak would actually have little impact on transformer temperature, and
thus on aging and deterioration. An improvement over this criterion is 2.2. Presentation of the Thermal Model
thus to use a criterion with thermal limits [3-4], [10] (most importantly on
In a planning study, in order to compute the hot-spot temperature ,
the hot-spot temperature), which makes it possible to benefit from the
one may use the model from the standard IEC 60076-7, presented on
thermal inertia of the transformer, in particularly the thermal inertia of its
Fig. 1. In this paper, we take this well-known model for granted; that is to
oil. This criterion is the core of dynamic thermal rating (DTR) and real-
say, it is assumed that the model is accurate and it is used as a black-box,
time thermal rating (RTTR), two technologies that are studied in [10-13],
without discussing the physical meaning of the equations and parameters.
[16-17] and that allow to operate a transformer closer to its thermal limit
The important question of assessing the value of the various parameters
and to extend its useful life. These two approaches have much in common
contained in this model is discussed in [7].
with the method proposed in this paper, with the crucial difference that
Three temperatures are involved in the model: the ambient temperature
they do not address planning stage but only operations. The already
1+ - 1
/
distributed generation (DG) owners should take with respect to the DSO
& =& ℒ (, . 0
1+ 1+
in order to enforce the thermal limits. It will be seen that the thermal (4)
−1
criterion needs to be linear; so, the model from IEC 60067-7 will be
& =& ℒ1 - 23 4 − 6
1+ 1+ /
linearized in Section 4.
ℎ ℎ (5)
= 42°$ , which is
Winding hot-spot The hot-spot temperature is computed using the IEC model. The
120 120 120
temperature (°C)
would reach 120°$ when = 1 for a long time.
ambient temperature is assumed to be constant,
Note that 40°$ is a standard maximal ambient temperature [4], [14]. One
Top-oil chosen such that
105 105 105
temperature (°C)
may notice that the temperature is “smoothed” with respect to current, as
Three limits are thus defined; one of them is however always less current peaks are damped by the thermal inertia of the transformer.
& ≥ 15°$ for every type of transformers (see Appendix, Table 2). The
and transformer is fully loaded in all three test-cases, and the remaining
hosting capacity is zero. However, for the load and the PV curves, and
winding hot-spot temperature limit is thus always reached before the top-
oil temperature limit, and there is no need to study the latter. contrary to the wind curve, some hosting capacity margin does remain
So the limits for a large power transformer under normal cyclic loading according to the thermal criterion (2): the red temperature curve does not
3
reach the 120°$ limit. We conclude that for a transformer sized by a generators whose power output exhibits a clear plateau during sunny days
typical load or by a PV, criterion (1) may lead to slightly oversizing the (the power output seems to saturate, possibly because the power converter
transformer. However, for a wind farm (and by extension all plants except was slightly undersized with respect to the peak power of photovoltaic
PV) there is nothing to gain as the power may remain close to its modules). Also, Fig. 3 shows that the typical value of the hosting capacity
maximum for a long period. increase, for a pure mix of “normal” PV generators that exhibit no power
output saturation, would be in the range of 10 to 40 % for distribution
transformers, and 2 to 20 % for ONAN power transformers.
In a real distribution system however, PV generators would probably
be mixed with other types of generators and loads, solar irradiance would
be different from one place to another, etc. These impediments probably
explain why the hosting capacity gain associated with thermal limits
remains currently untapped: calculating the exact factor by which a given
transformer may be overloaded in order to account for its thermal inertia
would depend on many local details and prove difficult for planning
engineers. The sequel of the paper is thus motivated by the will to find a
practical way to include the thermal inertia of transformers in planning
studies.
maximal current using (2) is = 1.1 then the ratio ; is 1.1, and using
allowed value. The histogram from Fig. 3 gives an idea of the possible
benefit for a transformer loaded by solar panels. Provided that the
the thermal criterion yields a hosting capacity increase of 10 %. technical details of the transformer and PV plant are accurately known,
the only caveat of this approach is that the ambient temperature might be
difficult to forecast. Techniques are presented in [4] to tackle this
problem.
Frequency of occurence
the transformer over its entire lifetime, with sufficient accuracy to directly
compute the temperature using IEC model [3]. As observed above, this is
probably the reason why the hosting capacity gains that may be unlocked
by using thermal limits at planning stage are currently untapped: DSOs
generally do not use the thermal criterion (2) today, but use a limit on
rated current or power.
We argue that the current impossibility to use a thermal criterion at
Fig. 3. Histrogram of the benefit of thermal criterion γ for 518 PV.
planning stage fundamentally results from the following fact: the
The histogram of γ spreads over a wide range of values. Looking into expression of hot-spot temperature from the IEC model [3], described by
the details of the simulations that yielded extreme values of γ, we Fig. 1, is not proportional to the current. Indeed, with a linear expression,
observed that the lowest values of γ were obtained for photovoltaic it would be possible to study separately the impact of each generator and
4
1
specific to the considered transformer, defined by (9).
ℎ: = +
wind generators on Fig. 2.
1+ 1+
(7) • Secondly, the oil-time constant, or more exactly the product : a
:=
larger time constant makes smoother and is thus beneficial.
1+
(8) • A third crucial parameter is α. It is also determined by the physics of the
2
& : −& = −1
considered transformer; the higher its value, the more is weighted in
:= 1+
& =
ℎ (9) , which is beneficial.
ℎ 4.3. The Linear Thermal Criterion in Practice
Parameter α expresses the relative weight of both variables and Most importantly, the new thermal criterion is a linear expression, an
with weights 1⁄ 1 +
(see typical values of in Table 2). Equation (7) expresses that is the indispensable property to decouple the study of various loads and
DE 89: ≤ 1.3
B ?
?
CE 89:F ℎ,? G ≤1
(11)
B
A ?
For DSOs the latter formulation (11) is even easier to use than (10):
89: ? and a maximum thermal load factor 89:F ,? G, that only have
to each generator or load, we simply associate a maximum current
it is impractical to contract with all customers as a whole and require them calculated by linearizing in the vicinity of its nominal value. The
to collectively ensure that constraint (10) will be satisfied, it is much more method is dependent on the IEC and IEEE models and its accuracy. The
practical to contract with each generator individually, which is what linearity of the proposed criterion is a crucial characteristic in practice, as
formulation (11) makes possible. it allows the DSO to contract with each generator and load individually. A
case study gives an example of how to use this linear criterion in practice,
In practice, computing ,? would be straightforward and could be a
and the expected benefits. The proposed criterion can also be used for
new feature embedded inside a standard smart meter without requiring
existing transformer in order to size a new PV installation. The paper also
any thermal sensor. The meter could be set to open its internal circuit
provides several figures on the hosting capacity gains which can be
breaker whenever constraint (11) is reached, allowing for a practical
expected by sizing transformers with thermal limits. Benefits are
enforcement of the thermal limit. Observe that this scheme would allow
substantial when the transformer is sized by generation including PV,
enforcing the global constraint — not overheating the transformer — by
which will be more and more common in the future. In the end, using this
coordinating the actions of individual customers without communication.
criterion would facilitate the installation of low-carbon technologies like
This is in contrast with the technology of real-time thermal rating [12-13]
PV in locations dominated by generators, as it would reduce the cost
where thermal constraints must be managed in real-time; implementing
induced by transformers.
such a scheme requires that some “active resources”, loads and/or
generators, are available to the DSO upon request, and thus calls for a
6. Appendices
communication link with these resources.
6.1. Simplification of the thermal model
4.4. Case Study: Transformer Sized by Generation
According to the value of and from Table 2, for every type of
A HV/MV power transformer has several feeders. The first is medium and large power transformer, equation (12) holds.
⁄ =
connected to a 10 MW plant with wind generators, the second to a 10 MW
plant with PV generators, and all the other feeders are only connected to (12)
loads, with a maximal consumption of 15 MW, and an estimated minimal
For a distribution transformer, (12) no longer holds but, as is equal
consumption of 3 MW. To simplify, a constant voltage will be assumed so to 1, the fraction (in the upper-right block on Fig. 1) using N⁄ 22
K can be assimilated to a power. Then, the two sizing methods can be vanishes. It can be deduced that N⁄ 22 may be replaced by 11 N for all
compared: types of transformers.
• The current method, based on power limits. Both plants would commit Table 2 Example characteristics related to the load ability of transformers
x
commit to a 10 MW capacity on both and , but the PV plant
0,8 0,8 0,8 1 1
y
would rather commit to a 10 MW capacity on and to a 9 MW capacity
Recommended Values
on . The capacity would be computed beforehand, based on the 1,6 1,3 1,3 1,3 2
time-profile of the power generated by the plant, and on the transformer 1 0,5 0,5 1 1
parameters (given by the DSO). The value of 9 MW is an example
1 2 2 1,3 1
∑ = 10 + 10 − 3 = 17 JK and ∑ = 10 + 9 − 3 = 16 JK. So,
value consistent with Fig. 3. Then the power balance gives:
2 2 2 1 1
min
a transformer with a nominal power of at least 16 MW would be enough
180 210 150 90 90
min
to satisfy criterion (11).
4 10 7 7 7
R
In the end, 1 MW are gained on the sizing of the transformer by
considering its thermal capacity instead of its power capacity. In this 5 6 6 6 6
∆ °$
example, gains are only linked to the PV generation. Further gains might
Example
23 26 26 22 29
∆ °$
be achieved, but they are expected to be either negligible, either too
complicated: 55 52 52 56 49
• The of the load might also be taken into account and would α 2 0.55 0.55 2.4 1.5
potentially be strictly lower (e.g. -3.1 MW instead of -3 MW).
• Criterion (10) would give higher gains than (11), but it would be too Two time constants still need to be managed: and . It can
complicated, if not impossible, to use in planning studies — this is be noted that for all transformer types, the winding time constant is
because as stated above, it does not decouple individual loads and much lower (a few minutes, see Table 2) than . In addition, for the
generators. purpose of sizing transformers, we consider that load variations over a
timeframe of a few minutes are too short regarding the dynamics involved
5. Conclusion (see e.g. Fig. 2); its impact would be negligible. For these reasons,
1
is neglect by approximating
≅ 1,
This paper proposes a method to size a power transformer using its
1+
thermal limit directly at planning stage. The method relies on the linear (13)
criterion (11) which is based on so-called thermal load factor . is
and the expression of (14) is obtained.
an image of the hot-spot temperature and its expression has been
6
−1
≅& ℒ1 - 23 4 − 6
1+
• to further simplify the mathematical expression of the constraint,
ℎ ℎ
• and to take the ambient temperature into account only implicitly,
1+ - 1
/
+& ℒ (, . 0
1+ 1+
through the fact that the rating of the transformer is assumed well
(14)
1 −1
chosen for the geographical area.
+ & = ⋅ ℒ1 − 13 4 − 6
1+ 1+
2 1
ℎ
+& : ℒ1 − 13 ≤0
(18)
1+ 1+
6.2. Linearization
ℒ1 − 13
ℒ1 − 13 + ≤0
As a second step, (14) is further simplified by linearizing it. Since, in
1+
the context of sizing a transformer, it is important to get a good behaviour (19)
- + - ≤1
−1 1+ 1+
(20)
≅& ℒ11 + = - − 1 34 − 6
1+ - ≤1
2 1
ℎ ℎ
+& ℒ T1 + : - −1 U
ℎ (21)
1+ 1+
(15)
1
References
+
1+ [1] ‘Enedis-PRO-RES_06E’, Enedis Technical Document [Online]
[Link]
RES_06E.pdf, assessed 08 August 2018
[2] Willis HL. Power Distribution Planning Reference Book, Second
Edition. CRC Press; 2004.
[3] Power Transformers Part 7: Loading Guide for Oil-Immersed
Power Transformers, International standard IEC 60076-7, 2005.
[4] IEEE Guide for Loading Mineral-Oil-Immersed Transformers and
Step-Voltage Regulators - Redline. IEEE Std C5791-2011
(Revision of IEEE Std C5791-1995) - Redline 2012:1–172.
[5] Pezeshki H, Wolfs P. Impact of high PV penetration on
distribution transformer life time. 2013 IEEE Power and Energy
Society General Meeting (PES), 2013, p. 1–5.
doi:10.1109/PESMG.2013.6672074.
[6] Turker H, Bacha S, Chatroux D, Hably A. Low-Voltage
Transformer Loss-of-Life Assessments for a High Penetration of
Plug-In Hybrid Electric Vehicles (PHEVs). IEEE Transactions on
Power Delivery 2012;27:1323–31.
doi:10.1109/TPWRD.2012.2193423.
[7] Gao Y, Patel B, Liu Q, Wang Z, Bryson G. Methodology to assess
distribution transformer thermal capacity for uptake of low carbon
denoted by :
(and that still depends on through a 1st order differential equation) is [9] Dias FM, Canizes B, Khodr H, Cordeiro M. Distribution networks
planning using decomposition optimisation technique. IET
- := & - +& - + ℒ -1 4 6
Generation, Transmission & Distribution 2015;9:1409–20.
1+
(16) doi:10.1049/iet-gtd.2014.0860.
ℎ ℎ
[10] Lachman MF, Griffin PJ, Walter W, Wilson A. Real-time
dynamic loading and thermal diagnostic of power transformers.
Then, (17) is obtained by combining (15) and (16). IEEE Transactions on Power Delivery 2003;18:142–8.
≅ doi:10.1109/TPWRD.2002.803724.
−1
+& = ⋅ ℒ1 − 13 4 − 6
ℎ ℎ
[11] Douglass DA, Edris AA. Real-time monitoring and dynamic
1+
thermal rating of power transmission circuits. IEEE Transactions
2 1
ℎ
(17) on Power Delivery 1996;11:1407–18. doi:10.1109/61.517499.
+& : ℒ1 − 13
1+ 1+
[12] Degefa MZ, Humayun M, Safdarian A, Koivisto M, Millar RJ,
Lehtonen M. Unlocking distribution network capacity through
real-time thermal rating for high penetration of DGs. Electric
temperature is the constant value 120°$. Assuming that the rating of the
We started from constraint (2) in which the maximum hot-spot Power Systems Research 2014;117:36–46.
doi:10.1016/[Link].2014.07.032.
transformer was well chosen, this hot-spot temperature limit should not be [13] Greenwood DM, Gentle JP, Myers KS, Davison PJ, West IJ, Bush
JW, et al. A Comparison of Real-Time Thermal Rating Systems in
≤ 120°$.
exceeded for any plausible ambient temperature curve in the considered the U.S. and the U.K. IEEE Transactions on Power Delivery
≤
Thus, constraint [14] ‘Power and Distribution Transformers’, in Low Voltage and
by (18) or more compactly (21). This approximation allows Medium Voltage Products, Schneider Electric, Catalogue 2015.
[15] Biçen Y, Aras F, Kirkici H. Lifetime estimation and monitoring of
both: power transformer considering annual load factors. IEEE
7
Transactions on Dielectrics and Electrical Insulation [17] Yip HT, An C, Lloyd JG, Taylor P, Michiorri A, Jupe S, et al.
2014;21:1360–7. doi:10.1109/TDEI.2014.6832284. Dynamic thermal rating and active control for improved
[16] Degefa MZ, Humayun M, Safdarian A, Koivisto M, Millar RJ, distribution network utilisation. Managing the Change, 10th IET
Lehtonen M. Unlocking distribution network capacity through International Conference on Developments in Power System
real-time thermal rating for high penetration of DGs. Electric Protection (DPSP 2010), 2010, p. 1–5. doi:10.1049/cp.2010.0213.
Power Systems Research 2014;117:36–46.
doi:10.1016/[Link].2014.07.032.