CC BY
38
UDC 621.318.001.4
TOPOLOGICAL TRANSIENT MODELS OF THREE-PHASE FIVE-LIMB
TRANSFORMER
ZIRKA S. E. MOROZ Y. I. ARTURI C. M. BONNMAN D.
Professor of the Department of Physics and Technology of the Dnipro National University, Dnipro, Ukraine, e-mail: [email protected];
Ph.D., Associate professor of the Department of Physics and Technology of the Dnipro National University, Dnipro, Ukraine, e-mail: [email protected];
Professor of the Department of Electronics, Information and Bioengineering of the Politechnico di Milano, Milan, Italy, e-mail: [email protected];
Ph.D., ABB AG, Bad Honnef, Germany, e-mail: [email protected].
Purpose. To show capabilities of topological models of three-phase, five-limb transformer to correctly represent transformer operation in regimes with high flux densities in the core. As a practically important example, time domain response of transformer subjected to geomagnetically induced currents (GIC) is analyzed and compared with results of a comprehensive field experiment.
Methodology. Transformer magnetic model, which takes into account geometry of the core and windings, is transformed in a dual electric equivalent scheme, whose transient is calculated by using EMTP-ATP. The results obtained demonstrate the importance of incorporating the positive and zero sequence impedances ofpower network.
Findings. A simple and reliable model of five-limb transformer is proposed. It was found that the presence of the transformer tank can be effectively accounted for by linear inductances representing the paths of the off-core fluxes from yoke to yoke. The modeling of GIC events represented in the paper is the most accurate ever obtained for three-phase, five-leg transformers. The model is validated by close agreement of the predicted values and waveforms of the phase currents and reactive power with those measured in tests performed on two 400 MVA transformers connected back-to-back and to a 400 kV power network.
Originality. It is shown that a simplified non-hysteresis model developed portrays the behavior of the of five-limb transformer under GIC condition with the same good accuracy as its hysteresis model. Both the transformer models are well grounded. So, the paper dispels some misconceptions about the influence of the hysteresis properties of the core material and tank in transient modeling of five-leg transformers.
Practical value. The practical value and significance of the paper is caused by the fact that the model proposed is a simple and reliable tool for power system studies. The paper warns of using unnecessary complicated models whose parameters are difficult to be estimated by experiments or calculations.
Key words: five-limb transformer, topological transient models, transients, off-core magnetic fluxes, current waveforms, reactive power, experimental validation, back-to-back transformer test.
for large transformers remained open up to now, although
I. INTRODUCTION
Over the last decades, a common practice in transient modeling of large power transformers is the use of their topological models [1] - [6]. The advantage of these low-frequency models is their ability to reproduce individually dynamic processes in the core limbs and yokes, as well as the off-core magnetic fluxes outside the windings. The latter feature becomes especially important when the core approaches saturation and some part of the flux leaves the core and closes through the air (oil) and transformer tank. The role of the off-core fluxes and thus the necessity for a detailed tank model is different for different transformer constructions. The off-core flux is significant in three-limb three-phase transformers, so a transient model of the tank was proposed in [4]. The lesser role of the tank in five-limb transformers is caused by the presence of the end (lateral) core limbs providing the paths for the unbalanced (zero-sequence) flux.
Perhaps the first topological model of the five-limb transformer, which covers the range of high flux densities, was proposed in [1], in which a discussion was raised concerning the necessity to model transformer tank. Because of the lack of experimental data, the “tank problem”
systematic measurements on five-leg transformers with saturated cores were reported in 2002 [7]. These measurements were carried out on two 400 MVA units connected back-to-back and to the Fingrid power network to simulate conditions experienced by transformers in the presence of geomagnetically-induced currents (GIC). The distinguishing feature of the work described in [7], is a large (up to 200 A) DC current injected in the transformers’ neutrals and a considerable reactive power (55 Mvar) consumed by each of the transformers in this extreme regime.
Although the mentioning of the comprehensive experiment in [7] has become customary in the literature, there was no successful attempt to reproduce it with a model. An effort to repeat experimental results of [7] has been recently undertaken in [8]. However it will be shown in Section V that the outwardly plausible results represented in [8] were obtained due to deviations from the actual experimental conditions and hence cannot validate the “advanced” model [9], [10].
The aim of the present paper is to show that the measurements in [7] can be portrayed accurately by topological models developed in the framework of the conven-
© Zirka S. E., Moroz Y. I., Arturi C. M., Bonnman D., 2017
DOI 10.15588/1607-6761-2017-2-2
tional approach [1]-[5] if the network impedances are taken into account.
II. THE BACK-TO-BACK TRANSFORMER TEST
The experimental setup in fig. 1 consists of two similar YNyn0d11 400/400/125 MVA full transformers, T1 and T2, with rated voltages 410/120/21 kV. The medium-voltage (MV) 120-kV windings are nearest to the core, the high-voltage (HV) 410-kV windings are in the middle, and the outermost are the low-voltage (LV) 21-kV delta-connected windings. The turn numbers NM, NH, and NL in these windings of T2 are 224, 766, and 68.
The welding generator G serves to inject the dc biasing current Idc into the HV neutrals of both T1 and T2, thus providing the DC currents Idc/3 in each of the HV windings.
Figure 1. Experimental setup in [7].
Important elements of both the experimental setup and the corresponding model are shunt capacitors C, which form the AC earthing of the transformers and isolate the generator G from earth. Two current transformers are to monitor the variables designated in fig. 1. A summary of other transformer data and measurements can be found in [7].
III. TOPOLOGICAL TRANSFORMER MODEL
It was already noted in the literature that a lot of topological models exist for a given transformer configuration [11]. Some of them are unnecessary complicated and contain parameters difficult to be estimated by experiments or calculations. An example is the division of the yokes into sections [9], [11]; another one is the planar model in which the cover and bottom of the tank are represented by one-dimensional elements that results in the “tank”, which consists of only the left and right walls [8]. As explained in [12], questionable is also the lumped character of the tank elements, and that each wound leg of the core has its own “tank” [10].
The tractable model employed in this paper is based on the magnetic circuit in fig. 2, which is a development of that described in [5]. To save space, only the elements related to the leg A are mainly shown in fig. 2. In particular, the sources of magnetomotive force (MMF) Fm, Fh, and Fl represent the MV, HV, and LV windings. The magnetic flux paths in the core are shown as solid rectangles, which represent the legs (elements Ra, Rb, Rc), yokes (Rab and Rbc) and the end limbs (Rd and Rf). The innermost channel Leg-MV and the equivalent leakage channels
MV-HV and HV-LV are characterized by the linear reluctances R01, R12, and R23 respectively. The negative (fictitious) reluctances Rp are added to match all three binary short-circuit inductances [2].
R04 R04
Figure 2. Magnetic (reluctive) transformer model.
Linear reluctances R03 represent the flux paths from yoke to yoke, which are outside the windings. Reluctances R04 are for fringing flux paths in parallel to the yokes. Reluctances Rg = g/(^o^ieg) take into account the air gaps g at the core joints, which are related to the legs with cross-section Sleg.
The electrical (duality-derived) equivalent of the model is shown in fig. 3 between the nine ideal transformers (ITs). The linear inductances L in fig. 3 are indexed with the same symbols as the linear reluctances R in fig. 2, and L = NH /R . The seven hysteretic elements are the ATPDraw implementations of the dynamic hysteresis model (DHM) [13]. The 1:1 turn ratio of three ITs at HV terminals points out that the model parameters are referred to NH turns. So, the turn ratios n of the ITs at MV and LV terminals are NM/NH and NH/NL respectively.
As in many practical cases, it was not possible to determine the material of the core, so it was first assumed
that the cores of both transformers are assembled from grain-oriented steel 27ZDKH85. Then it was observed that the same results (the ones presented in the paper) are obtained when using AK steel H1 in transformer T1 and steel 27ZDKH85 in the T2 unit (the DHM-inductors of both these steels can be taken at [14] or in the current version of EMTP-ATP [15]).
The lengths of the legs, yokes, and lateral limbs of T2 are 3.340 m, 2.421 m, and 6.257 m respectively, their cross-sections are 0.8309 m2, 0.4959 m2, and 0.4155 m2. Transformer T1 has 11 percent less cross sections and, respectively, larger turn numbers.
The percentage short circuit reactances of T2 provided in [7] yield the following binary short-circuit inductances referred to the HV side: LS12 = 263.5 mH, LS23 = 540.4 mH, and LS13 = 877.5 mH. The winding resistances rb r2, and r3 are brought outside the inductive part of the model. Their values in T2 are 0.34 Q, 28 mO, and 16 mQ respectively. Resistances of T1 are 30% greater. Resistances Rinf = 109 Q make the MV windings effectively open-circuited and the delta LV windings unloaded, according to [7].
To match LS12, LS23, and LS13, the conventional star-connected inductances L12, L23, and Lp are used in the model of fig. 3. The value of the negative inductance Lp is
(LS12 + LS23 - LS13)/2 [4] then L12 = LS12 - Lp, 3nd
L23 = LS23 - Lp.
In the absence of winding design data, inductance L01 of the innermost channel Leg-MV was related to the short-circuit inductance LS12 of the next channel MV-HV. Following [3], the ratio k01 = L01/ LS12 can be evaluated by the insulating clearances of these channels, which are assumed proportional to the voltages across the channels. Taking into account the zero potential of the core, k01 = VM/(VH - VM) = 0.413, which is rounded to 0.42.
The role of inductance L04 was explained in [4]. Since the modeled five-limb transformers have reinforced yokes (Ayoke = 0.6 A[eg), the variation of L04 in the range [0, LSJ2] does not affect the model behavior, and L04 was set equal to LS12.
So, the only fitting parameter of the model is the value of L03. Similarly to L01 and L04, it is related to LS12 by the ratio k03 = L03 /LS12. In general, all three inductances L03 can be chosen individually, but for simplicity, they are assumed to be the same. The iterative search of L03 was resulted in k03 = 14 for transformer T2 and k03 = 10 for T1, which will be used in all calculations below.
IV. MODELING TECHNIQUE
The initial stage of the modeling is the fit of the transformer models to the measured no-load losses of T1 (172 kW) and T2 (100 kW). As proposed in [5], the model fitting to these losses is carried out by choosing the coefficient Kbss of the DHM. The reactive and hence the apparent power in the rated no-load regime can be fitted by changing the core air gaps g. We do not dwell on these details, because, as shall be shown below, the no-load losses do not practically influence the transformer behav-
ior under large GlCs.
As will be shown in the next section, it is important to take into account the positive (Zi) and zero-sequence (Z0) impedances of the power network. According to Fin-grid evaluations at the time of the test [7], the per phase impedance Z1 at f =50 Hz is formed by R1 = 10.54 Q and L1 = 181.3 mH. The impedance Z0 was evaluated using the frequency content of the current in the earthling conductor of the equivalent network source in fig. 1. It was found by a Fourier analysis of the calculated earth current (current iearth in fig. 1) that the third (150-Hz) harmonic distinctly dominates in its waveform. This is because, in the configuration considered, the current returning to the remote generator is the sum of currents in the neutrals of T1 and T2, and both of them are sums of three strongly asymmetrical (phase) currents shifted by 120 degrees that causes three positive and three negative peaks in their waveforms. This means that the network impedance Z0(3) = 748.7 exp(j66.3°) = R0(3) + j2nf3 L0(3) calculated by Fin-grid for f3 = 150 Hz should be used.
If the solver employed has a standard n-equivalent of the network, the values of R1, L1, R0(3), and L0(3) are entered directly into the window of this component. An alternative way of accounting for Z0 is the use of the network equivalent in fig. 4 in which the boxes T1 and T2 contain transformer model shown in fig. 3.
V —1=1—
(\ p* isum iT1 Г 1
. Network 1
ln 2 equivalent 1
rn [
T1
iT2
T2
:r©-^4
Edc Ri
Figure 4. Model configuration employed.
Using the values R0(3) = 300.7 Q, L0(3) = 727.5 mH, and the generic formula for ground-return impedance, ZN = (Z0 - Z1)/3 [16], [17] the following network parameters should be used in fig. 4: RN = 96.7 Q and LN = 182.1 mH. The elevated source voltage (V = 416 kV) is to compensate for the voltage drops over the network and provide the voltage (~ 404 kV) observed across the modeled transformers at 200-A GIC. The value of the shunt capacitors C has not been documented, but this is not an influencing parameter if C > 100 qF. The current Idc is determined by the voltage Edc of the generator G, its internal resistance Ri, and resistances r11 and r21 of the HV windings of T1 and T2: Idc = Edc / (Ri + r11/3 + r21/3).
In all the simulations below, a three-phase ramped voltage of the network is initially used in the model to establish symmetrical flux densities in the legs of both the transformers. To observe dynamics of GIC events, it is supposed that the GIC considered is initiated by a step voltage of generator G in fig. 4, which is switched on at t = 2 s. The following transient is characterized by a growth in the currents and the reactive power consumed by T1 and T2. The calculated rms values of the currents in phas-
es B and C of T2 as well as the reactive power Q2 consumed by this transformer are shown in fig. 5. In steady state, the rms currents in fig. 5 reach 108 and 124 A, which are close to the measured values of 110 and 125 A (for phases B and C respectively). The steady state reactive power Q2 is 53.1 Mvar that also is close to the measured 55 Mvar.
The waveforms of all three phase currents in the steady state (calculated in the presence of 200-A DC current) are shown in fig. 6. They agree quite well with the measured currents in fig. 7. The same good agreement was obtained for the current waveforms in the network (three-phase currents isum in fig. 4) so we do not show them in the paper.
Figure 5. Calculated RMS currents in phases B and C of transformer T2 and the reactive power Q2 consumed by T2 during the transient.
Figs. 8, 9, and 10 show calculated flux densities in the core leg, yoke, and the end limbs respectively.
Figure 8. Flux density in the leg A.
a
Он
й
3
JJ
3
и
50
0
-50
-100
-150
-200
-250
-300
-350
Д rr\. Д\ f? л\ . /А ґ\ /;
M \ |ii. \^/ \ / \Jf
M I .. I , і.. і M... w r ...M | і м
1 ! і 1 1 ' ........ 1 i і і 1
Ill 1 ! \ ...i 1 ' ' / > / і і \
1 'I і/. \. 1 1 1 \.../... ті \
і ...U v l.B......! \ і \
1 \J A i 1 V Ci У I 1
0.00
0.01
0.02
0.03 t [s] 0.04
Figure 6. Calculated phase currents in T2.
50
^ 0
g -50 Ь
3 -100
<D
§ -150
£ -200
<D
S -250
сл
I -300 -350
0.00 0.01 0.02 0.03 t [s] 0.04
Figure 7. Measured phase currents in T2.
Figure 9. Flux density in the yoke AB.
0 20 40 l Ы
59.94 59.96 59.98 t [sj 60.01
Figure 10. Flux densities in the outer limbs D and F during (a) transient and (b) steady state.
These demonstrate their different waveforms during the transient and steady state. As can be seen in fig. 10(b), the end limbs D and F are saturated non-simultaneously, and there is no point of time at which both BD and BF drop below -2 T. This shows that the zero-sequence flux closes mainly within the core and explains why there is no need for a detailed tank model when the five-leg transformer is considered as network element.
V. UNREALISTIC MODELING IN [8]
It should be recalled at this point that the influence of the network parameters (impedances Z1 and Z0) was not immediately evident, and considerable unsuccessful efforts were made to obtain the realistic waveforms in fig. 6 using the initial model configuration in fig. 1. The power network is represented in fig. 1 by an ideal voltage source, which was also used in [8]. The current waveforms in transformer T2 calculated for the initial configuration are shown in fig. 11.
< о
1 -100
2
0
1 -200
f -300 ”3
-400
0.00 0.01 0.02 0.03 t [s] 0.04
Figure 11. Phase currents in T2 calculated with the ideal network.
The currents in fig. 11 are markedly different from the measured currents in fig. 7. All our attempts to improve the waveforms in fig. 11 by changing the structure and parameters of the transformer model were in vain, pointing out that the reason lies outside the transformer representation.
<:
о
50 0 -50 -100 -150 -200 -250 -300 -350
0.00 0.01 0.02 0.03 t [s] 0.04
Figure 12. Phase currents in T2 calculated with ideal network and in the absence of earthing capacitors used in the measurements.
A pragmatic measure to circumvent this difficulty
was observed by one of the co-authors (D.B.) and then used independently in [8]. This artificial measure is to substitute capacitors C in the scheme of fig. 1 with large resistances. The current waveforms calculated in the absence of capacitors C and for the idealized 400-kV network are shown in fig. 12.
The currents in fig. 12 and similar waveforms in [8] outwardly resemble the measured currents in fig. 7, but we should stress again that they were obtained using the model configuration different from that in the actual test
[7] and hence cannot validate the “advanced” model [9], [10].
It should be noted that plausible currents in the HV windings of transformer T2 can also be simulated in the model configuration in fig. 13, which contains only one this transformer. In accordance with [8], an ideal voltage source V is used in this configuration (Zi=0), and there are no grounding capacitances C. The latter simplification means that the internal inductance L of the welding generator G can be included in series with its resistance R1. To obtain the needed current Idc in the neutral, the e.m.f. Edc in fig. 13 is set equal to (rn/3+Rj) Idc.
iT2
Figure 13. Model configuration with one transformer T2.
Figure 14. Phase currents of T2 calculated with the model configuration in fig. 13.
Steady-state currents in T2 calculated at V = 400 kV and Lj = 400 mH are shown in fig. 14. Although, in comparison to the waveforms in fig. 11, they are closer to the measured currents in fig. 7, these curves have flat segments and “inverted convexities”, which are absent in the measured waveforms. The presence of similar segments
ISSN 1607-6761 (Print)
ISSN 2521-6244 (Online)
and convexities in the currents shown in [8, fig. 7] allows one to suppose that the analysis in [8] could be carried out for a model configuration similar to that in fig. 13, but not for the two-transformer arrangement in fig. 1.
Although the model configuration in fig. 13 does not reflect the actual test conditions in [7] where the impact of Li is cancelled by capacitors C, the favorable influence of Lj in fig. 13 should be discussed further. It can be seen that all the elements of the circuit in fig. 13 are connected in series, so inductance Li can be moved into the grounding conductor of the voltage source. In this case, inductance Lj takes the place and role of inductance LN in the circuit of fig. 4. This explains the need for accounting for the zero-sequence impedance of the network implemented in the model configuration in fig. 4.
In concluding this section, we mention a separate study conducted to estimate the influence of hysteretic properties of the core material. With this purpose, each of the DHM-inductors in the model of fig. 3 was replaced by a lossless non-hysteretic inductance. The flux-current relationships of these inductances (different for legs and yokes) were obtained from the “middle” (anhysteretic) curves of the steel employed. It was found by numerical experiments that all the results calculated with nonhysteresis and hysteresis models practically almost coincide. It is remarkable that the coincidence was also obtained in portraying the “hysteresis loop” shown in fig. 10 of [7]. This shows that this loop is not caused by hysteresis properties of the core, but is a kind of Lissajous figure occurred due to the phase shift between the flux density of the leg and corresponding terminal current.
When the modeling of GIC events had been completed, it turned out that the model also predicts accurately the zero sequence impedances of the transformers measured from both HV and MV sides (they are 40.1% and 65.2% respectively).
VI. CONCLUSION
In addition to accuracy, an important feature of the transformer transient model intended for network studies is its simplicity and ease of implementation. It is always desirable to understand the capabilities of the models developed within the well-established framework, before proposing new, often unnecessary complicated tools. With this idea in mind, we have considered the capabilities of the conventional topological model of five-leg transformer in simulating transformer behavior at high flux densities in the core legs and yokes. In fitting the model, we used field test results from two 400 MVA transformers with 200-A DC currents entering their neutrals. It was found that accurate modeling of the test in [7] requires the network impedances to be taken into account, as well as the actual back-to-back connection of the transformers considered.
It was found that the presence of the tank can be effectively accounted for by linear inductances representing the off-core fluxes from yoke to yoke. This finding is not unexpected because the same conclusion was made in
studying inrush current events in a three-limb Yyn transformer [5] in which an accurate tank model was anticipated due to the core design and higher flux densities in the core.
Having started with hysteretic transformer model, we have ended up with a simplified non-hysteretic model, which reproduces transformer currents and reactive power with the same good accuracy. This means that in its an-hysteretic form, the model can be implemented in any computing environment, and the data provided in the paper is sufficient to repeat the calculated waveforms with the use of any circuit simulator.
The model is mainly a tool for power system studies, in which the application of finite-element models is inappropriate because of the computational cost and the presence of several transformers in the studied network.
REFERENCES
[1] . Arturi, C. M. (1991). Transient simulation of a three
phase five limb step-up transformer following an out-of-phase synchronization. IEEE Trans. Power Delivery, 6, 1, 196 - 207. DOI: 10.1109/61.103738.
[2] . Chen, X., Venkata, S.S. (1997). A three-phase three-
winding core-type transformer model for low-frequency transient studies. IEEE Trans. Power Delivery, 12, 2, 775 - 782. DOI: 0885-
8977/97/$10.00.
[3] . Mork, B. A., Gonzalez, F., Ishchenko, D., Stuehm, D.
L., Mitra J. (2007). Hybrid transformer model for transient simulation - Part I: Development and parameters. IEEE Trans. Power Delivery, 22, 1, 248 -255. DOI: 10.1109/TPWRD.2006.883000.
[4] . Zirka, S. E., Moroz, Y. I., Arturi, C. M. (2014).
Accounting for the influence of the tank walls in the zero-sequence topological model of a three-phase, three-limb transformer. IEEE Trans. Power Delivery, 29, 5, 2172-2179. DOI: 10.1109/
TPWRD.2014.2307117.
[5] . Zirka, S. E., Moroz, Y.I., Hoidalen, H. Kr., Lotfi, A.,
Chiesa, N., Arturi, C. M. (2017). Practical experience in using a topological model of a core-type three-phase transformer - No-load and inrush conditions.
IEEE Trans. Power Delivery, 32, 4, 2081-2090. DOI: 10.1109/TPWRD. 2016.2618900.
[6] . Tikhovod, S.M. (2014). Modelirovanie perehodnyh
processov v transformatorah na osnove magnitojelek-tricheskih shem zameshhenija. [Modeling transformer transients using magnitoelectric equivalent schemes]. Elektrotehnika i elektroenergetika, 2, 59-68, (in Russian).
[7] . Lahtinen, M., Elovaara, J. (2002). GIC occurrences
and GIC test for 400 kV system transformer. IEEE Trans. Power Delivery, 17, 2, 555-561. DOI: 0885-8977(02)02750-4.
[8] . Rezaei-Zare, A. Marti, L., Narang, A., Yan, A.
(2016). Analysis of three-phase transformer response due to GIC using an advanced duality-based model.
IEEE Trans. Power Delivery, 31, 5, 2342-2350.
DOI: 10.1109/TPWRD. 2015.2505499.
[9] . Rezaei-Zare, A. (2015). Enhanced transformer model
for low- and mid-frequency transients-Part I: Model development. IEEE Trans. Power Delivery, 30, 1, 307-315. DOI: 10.1109/TPWRD.2014.2347930.
[10] . Rezaei-Zare, A. (2015). Enhanced transformer
model for low- and mid-frequency transients-Part II: Validation and simulation results. IEEE Trans. Power Delivery, 30, 1, 316-325. DOI:
10.1109/TPWRD.2014.2347934.
[11] . Lambert, M., Mahseredjian, J. (2013) Electromagnetic transient type transformer models for geomagnetically-induced current (GIC) studies. EPRI Report 3002000832.
[12] . Zirka, S. E., Moroz, Y. I., Rahimpour, E. (2017). Towards a transformer transient model as a lumped-distributed parameter system. Compel, 36, 3, 741750. DOI: 10.1108/COMPEL-09-2016-0389.
[13] . Zirka, S. E., Moroz, Y. I., Chiesa, N., Harrison, R.
G., Heidalen, H. Kr. (2015). Implementation of inverse hysteresis model into EMTP - Part II: Dynamic model. IEEE Trans. Power Delivery, 30, 5, 2233-2241. DOI: 10.1109/TPWRD. 2015.2416199.
[14] . Moroz, Y. I., Zirka, S. E. (2014). Inverse models of magnetic hysteresis, [Online]. Available: https://sites.google.com/site/inverse-hysteresismodel.
[15] . Alternative Transients Program, ATP-EMTP, (2016). [Online]. Available: http://www.eeug.org.
[16] . Tleis, N. D. (2008). Power systems modelling and fault analysis: Theory and practice. New York: New-nes/Elsevier, 625.
[17] . Evdokunin, G. A. (2016). Jelektricheskie sistemy i seti [Electrical systems and networks]. Saint Petersburg: Rodnaya Ladoga, 384, (in Russian).
Manuscript received 01.10.2017
ТОПОЛОГІЧНІ МОДЕЛІ ТРИФАЗНОГО П’ЯТИСТРИЖНЕВОГО
ТРАНСФОРМАТОРА
ЗІРКА С. Є.
МОРОЗ Ю. І. АРТУРІ Ц. М. БОННМАН Д.
д-р техн. наук, професор, професор кафедри систем автоматизованого управління Дніпровського національного університету, Дніпро, Україна, e-mail: [email protected];
канд. техн. наук, доцент, доцент кафедри систем автоматизованого управління Дніпровського національного університету, Дніпро, Україна, e-mail: [email protected]; професор факультету електроніки, інформації і біоінженерії Міланського Політехнічного інституту, Мілан, Італія, e-mail: [email protected]; д-р техн. наук, ABB, Бад-Гоннеф, Німеччина, e-mail: [email protected].
Мета роботи. Показати, що на основі топологічних моделей трифазного п ’ятистрижневого трансформатора можливо коректно відтворювати його роботу в режимах великих індукцій в осерді. В якості практично важливого прикладу, аналізується часовий відгук трансформатора на дію геомагнітно-індукованих струмів (ГІТ). Результати моделювання порівнюються з результатами натурного експерименту.
Методи дослідження. Магнітна модель трансформатора, що враховує геометрію осердя і обмоток, перетворюється в дуальну електричну схему заміщення, перехідний процес в якій розраховується в середовищі EMTP-ATP. Отримані результати демонструють необхідність урахування опорів прямої і нульової послідовності енергосистеми.
Отримані результати. Розроблено адекватну, просту і надійну модель п ’ятистрижневого трансформатора. Встановлено, що присутність бака трансформатора може бути врахована за допомогою лінійних індуктивностей, що представляють шляхи магнітних потоків поза осердям (потоків від ярма до ярма). Точність представленого моделювання процесів в трансформаторі в присутності ГІТ перевищує точність відомих моделей трифазних п ’ятистрижневих трансформаторів. Адекватність моделі підтверджується близькістю прогнозованих діючих значень і кривих фазних струмів, а також споживаної реактивної потужності, до відповідних величин, виміряних в експерименті, проведеному на двох 400 MBA трансформаторах, які були під ’єднані паралельно до енергосистеми напругою 410 кВ і послідовно відносно джерела постійної напруги.
Наукова новизна. Розроблено спрощену безгістерезисну модель п ’ятистрижневого трансформатора, яка відтворює поведінку трансформатора з такою ж високою точністю, що і гістерезисна модель, яка відрізняється тем, що вибір безгістерезисної моделі є науково обґрунтованим.
Практична цінність. Практична цінність і значимість статті обумовлено тим, що запропонована модель трансформатора являє собою простий і надійний інструмент для дослідження електричних мереж. Стаття застерігає від використання надмірно ускладнених моделей, параметри яких не можуть бути визначені в експерименті або знайдені шляхом розрахунків.
Ключові слова: п ’ятистрижневий трансформатор, топологічні моделі, перехідний режим, магнітні потоки поза осердям, форми струмів, реактивна потужність, експериментальна перевірка, послідовно-паралельне включення.
ТОПОЛОГИЧЕСКИЕ МОДЕЛИ ТРЕХФАЗНОГО ПЯТИСТЕРЖНЕВОГО
ТРАНСФОРМАТОРА
ЗИРКА С. Е.
МОРОЗ Ю. И.
АРТУРИ Ц. М. БОННМАН Д.
д-р техн. наук, профессор, профессор кафедры систем автоматизированного управления Днепровского национального университета, Днепро, Украина, e-mail: [email protected];
канд. техн. наук, доцент, доцент кафедры систем автоматизированного управления Днепровского национального университета, Днепро, Украина, e-mail:
профессор факультета элэктроники, информации и биоинженерии Миланского Политехнического института, Милан, Италия, e-mail: [email protected]; д-р техн. наук, ABB, Бад-Гоннеф, Германия, e-mail: [email protected].
Цель работы. Показать, что на основе топологических моделей трехфазного пятистержневого трансформатора возможно корректно воспроизводить его работу в режимах больших индукций в сердечнике. В качестве практически важного примера, анализируется временной отклик трансформатора на воздействие геомагнитно-индуцированных токов (ГИТ). Результаты моделирования сравниваются с результатами натурного эксперимента.
Методы исследования. Магнитная модель трансформатора, учитывающая геометрию сердечника и обмоток, преобразуется в дуальную электрическую схему замещения, переходный процесс в который рассчитывается в среде EMTP-ATP. Полученные результаты демонстрируют необходимость учета сопротивлений прямой и нулевой последовательности энергосистемы.
Полученные результаты. Разработана адекватная простая и надежная модель пятистержневого трансформатора. Установлено, что присутствие бака трансформатора может быть учтено посредством линейных индуктивностей, представляющих пути магнитных потоков вне сердечника (потоков от ярма к ярму). Точность представленного моделирования процессов в трансформаторе при наличии ГИТ превышает точность известных моделей трехфазных пятистержневых трансформаторов. Адекватность модели подтверждается близостью предсказанных действующих значений и кривых фазных токов, а также потребляемой реактивной мощности, к соответствующим величинам, измененным в эксперименте, проведенном на двух 400 MBA трансформаторах, которые соединены параллельно к энергосистеме напряжением 410 кВ и последовательно по отношению к источнику постоянного напряжения.
Научная новизна. Разработана упрощенная безгистерезисная модель пятистержневого трансформатора, которая воспроизводит поведение трансформатора с той же высокой точностью, что и гистерезисная модель, отличающаяся тем, что выбор безгистерезисной модели научно обоснован.
Практическая ценность. Практическая ценность и значение статьи обусловлено тем, что предложенная модель трансформатора является простым и надежным инструментом для исследования электрических сетей. Статья предостерегает от использования излишне усложненных моделей, параметры которых не могут быть определены в эксперименте или посредством вычислений.
Ключевые слова: пятистержневой трансформатор, топологические модели, переходный режим, магнитные потоки вне сердечника, формы токов, реактивная мощность, экспериментальная проверка, последовательно-параллельное включение трансформаторов.
