CC BY
44
UDC 621.3
doi: 10.20998/2074-272X.2016.4.06
D.V. Tugay
THREE-PHASE ENERGY SUPPLY SYSTEMS SIMULATION FOR THE TOTAL POWER LOSSES COMPONENTS ASSESSMENT
Purpose. The goal is to optimize a structure of Matlab-model of the three-phase energy supply system with power active filter. The mathematical model that describes the energy supply system modes of operation which contains additional losses is proposed. Methodology. We have applied concepts of the electrical circuits theory, mathematical modeling elements based on linear algebra and vector calculus, mathematical simulation in Matlab package. Results. We have developed two models of three-phase energy supply system. The first one is based on a vector representation, and the second one on the matrix representation of energy processes. Using these models we have solved the problem of maintaining unchanged the average useful power for 279 cases of energy supply system modes of operation. Originality. We have developed methods of mathematical analysis of a three-phase energy supply systems with polyharmonic voltages and currents in the symmetric and asymmetric modes. Practical value. We have created Matlab-model of a three-phase energy supply system with automated calculation of a correction factor. It allows reducing more than one order the time for energy processes elucidation in multiphase systems. References 11, tables 1, figures 13. Key words: energy supply system, power active filter, the minimum possible losses, total losses power, three-dimensional complex vector, Matlab-model of the three-phase energy supply system.
Цель. Целью статьи является оптимизация структуры Matlab-модели трехфазной системы электроснабжения с силовым активным фильтром с помощью математической модели, описывающей режимы работы системы электроснабжения, в которых возникают дополнительные потери. Методика. Для проведения исследований использовались положения теории электрических цепей, элементы математического моделирования, основанные на линейной алгебре и векторном исчислении, математическое моделирование в пакете Matlab. Результаты. Разработаны две модели трехфазной системы электроснабжения, первая, основанная на векторном представлении, а вторая на матричном представлении энергетических процессов, с помощью которых была решена проблема поддержания неизменной средней полезной мощности нагрузки для 279 случаев работы системы электроснабжения. Научная новизна. Получили дальнейшее развитие методы математического анализа режимов работы трехфазных систем электроснабжения с полигармоническими напряжениями и токами в симметричных и несимметричных режимах работы. Практическое значение. Создана Matlab-модель трехфазной системы электроснабжения с автоматизированным расчетом поправочного коэффициента, позволяющая более чем на порядок сократить время для исследования энергетических процессов в многофазных системах. Библ. 11, рис. 13, табл. 1.
Ключевые слова: система электроснабжения, силовой активный фильтр, минимально возможные потери, мощность суммарных потерь, трехмерный комплексный вектор, Matlab-модель трехфазной системы электроснабжения.
Introduction. Understanding the operation of three-phase energy supply systems (ESS) allows professionals in the field of power electronics to solve a number of problems related to the electromagnetic compatibility of electric power consumers from the industrial network as well as increasing the energy efficiency of such systems. Complex solution of problems presented in the current conditions based on the use of the active filtering means -power active filters (PAF). In the systems of energy supply of communal consumers, industrial plants, electric vehicles adaptation of PAF of parallel type is increased [1, 2]. Despite the simplicity of the power circuit which is a three-phase voltage inverter on the transistor modules, PAF are complex dynamic systems able to adapt to changes in the configuration and mode of operation of three-phase ESS. The feasibility of using PAF to reduce the losses in the ESS [3, 4] is due to the possibility of a preliminary assessment of the total power losses in the system before and after it is connected. The simplest way to perform such an assessment is a computer simulation of a three-phase ESS with PAF. This paper describes how to create a computer model to investigate diversity of power modes of the three-phase ESS with minimal time overhead.
A complicated extensive three-phase ESS can be replaced by a simple equivalent circuit shown in Fig. 1 [5].
Fig. 1. Equivalent circuit of the three-phase ESS with PAF
The scheme consists of three parts: a three-phase sinusoidal power source of sinusoidal (or nearly sinusoidal) voltages Source, a three-phase load Load and connecting cable line Line with resistance linear wire Rs and neutral wire resistance Rn. The three-phase load in the circuit of Fig. 1 may be any of: resistors, reactors, capacitors batteries, nonlinear consumers, current and voltage sources as well as the possible combinations of the elements listed. It should be noted that the inductances Ls and Ln of the cable line in the equivalent circuit are transferred in the load. Depending on the nature of the load and the operating mode of the ESS three power transmission modes are available: direct - energy is
© D.V. Tugay
transmitted from the source to the load, reverse - the energy is transferred from the load to the source, and mixed mode - in the period of recurrence it is possible to combine the first and the second modes. After the contactor SA switch on, the power active filter (block PAF) is connected parallel to the load.
In [6] a universal formula to determine the total losses power in the three-phase four-wire ESS through its components is proposed
APE* =
ap
E
P
= APm
= APm
usf
+ AP,
puls*
+ APq* + APn* +APm
+ AP
add *
■ q
Pusf = const,
mut*
(1)
=I 1 - _L
7max 2 + V 4
(2)
investigation of which confirmed the validity of (1) but was associated with significant time-consuming.
This work is dedicated to the creation of a universal Matlab-model of the three-phase ESS and optimization of its parameters for an accelerated assessment of the power losses components in any of the 288 possible ESS variants.
The structure of the three-phase ESS model.
Optimization of the structure of the Matlab-model of the three-phase ESS may be carried out at the transition from the circuit on Fig. 1 to the equivalent circuit of Fig. 2.
where APmin* is the relative minimum possible power losses determined by the absence of pulsations of the instantaneous active power and the vanishing of the instantaneous reactive power in three-phase ESS; APpub* is the relative components of additional power losses due to the transient component of the instantaneous power of the three-phase ESS; APq* is the relative component of additional losses power due to the instantaneous reactive power of the three-phase ESS; APn* is the relative power losses in the neutral conductor calculated in the period of repeatability caused by current flow in the neutral wire; hPmut* is the relative component of additional losses power due to the mutual influence of the electromagnetic processes in the phase conductors and the neutral wire of the three-phase ESS; APadd* is the relative additional losses power; Pusf is the average calculated in the period of repeatability useful load power the amount of which in the course of the experiment must be kept constant.
Calculation of the total power losses components is based on the concepts of modern theories of instantaneous active and reactive power [7-9] on the nature of the electromagnetic processes in three-phase ESS and algorithms of PAF control using matrix transformations of these theories. In the simplest case, if we consider the ideal active filter (with zero internal losses) then when it is connected to the three-phase ESS total losses in the system will be equal to the lowest possible, i.e. the ESS will operate with the highest possible efficiency [10]
ksc PsclPusf (3)
is the ratio of power of the three-phase resistive short circuit of the ESS to useful load power.
In above-mentioned publications [5, 6] it was shown that the validation of (1) is possible using the Matlab-model of three-phase four-wire ESS which has properties of universality and allows the calculation of power losses components at different modes of the ESS. Combinations of modes of three-phase voltage source, the nature of three-phase load, the neutral conductor resistance value and direction of the energy flow in the line wires allow to obtain 288 ESS variants, while in 279 ESS variants there are additional losses described by (1). The authors [5, 6] were able to check about 30 ESS variants, the
Fig. 2. Equivalent circuit of the three-phase ESS
The three-phase voltage source in this circuit consists of two series-connected sources - symmetrical sinusoidal voltage source usa, usb, usc and additional three-phase voltage source ua, ub, uc which is connected to the ESS by opening the bypass contactor SU. With the additional power source it is possible to specify amplitude or phase asymmetry as well as the addition of the higher harmonic components to the range of supply voltage. The three-phase load is modeled by two blocks - a symmetrical three-phase resistive load k¡-R¡ and regulated three-phase current source ja, jb, jc which is connected in parallel with the symmetrical resistive load after SJ contactor closing. Regulated current source provides in linear conductors the necessary current shape that will fit any linear or non-linear load with balanced or unbalanced phases load. The proportionality factor kl which is multiplied by the active resistances of all three phases is necessary to maintain a constant value of average active useful power of load Pusf = const as follows from (1) when making a particular disturbance in the system, and its definition is a separate problem.
Setting variable coefficients that determine the operation mode of the three-phase ESS. We return to the ESS equivalent circuits on Fig. 1 when the contactor SA is open and replace it with a simplified diagram shown in Fig. 3 explaining the cause of the power loss components in accordance with (1). Electrical quantities and parameters in the circuit of Fig. 3 are shown in the vector notation. The circuit consists of a voltage source the phase values of which are described by the
vector U = ku o us + uns o kun • su ;
resistance of the
line Rs - v , complex linear load resistance for the i-th harmonic component of the current which is described by
the vector zLi = k¡ • R¡ • ka + j • • kr + x!si • v j, and the
*
*
current source modeling the nonlinear load described by the vector kj ° J ■ Sj .
Fig. 3. Simplified circuit of the three-phase ESS for the parameters setting
In Fig. 3 the following notation is used: • the 3D complex vector of phase symmetrical sinusoidal voltages
" Us ■ ej
2rc
- j 'T
Us ' e 3
- j —
Us ' e 3
Us
U
U„
sb
(4)
where Us is the RMS phase voltage of the symmetric source;
• the 3D complex vector of phase voltage of the additional three phase source of non-sinusoidal voltage
]T V ' e}'^"ai J
i=2 \Ua
= —a Ub
i=2 Uc
¿ V ' ej i=2
(5)
J =
¿ (j, ' ej'v»J
i=2 \J ,
I (j-J = Jb
i=2 Jc
I (J-eJ' - J2 i=2
(6)
coefficient taking two values 0 or 1 depending on the switch SJ position in Fig. 2;
• the 3D complex vector of linear load resistances
zL, =
kl ' Rl ' kaa + j ' [xli ' kra + xs, kl ' Rl ' kab + j ' lxli ' krb + xs, kl ' Rl ' kac + j 'Vxli ' krc + Á
ZLai ' gJ ' ^ zai ZLa
ZLbi ej'<P zbi = ZLb
ZLci ej'<Pzci _Zlc _
(7)
where xli are the reactances for the i-th current harmonic; Xsi is the reactance of the cable line for the i-th current harmonic reduced to the load reactance; ZLah ZLbh ZLci are the modules of the impedance of linear load phases for the i-th current harmonic; ^za, are the phases of
corresponding impedances;
• the 3D vectors which store information about coefficients of amplitude asymmetry of the voltage source and non-linear load current, phase unbalance of active and reactive load as well as the value of the corresponding harmonic component of the input voltage range and a non-linear load ;urrent
ku ~ \kua kub kuc ] ,
k , k ]T
unb unc J
kj lT
kun \_ku
kj ~ l^ja
ka \kaa kr = \kr,
• the unit vector
"jb kab krb
kjc]T , k ]T
ac ,
krc] ;
(8) (9) (10)
(11) (12)
i
where U is the RMS of the i-th harmonic component of phase voltage of three-phase non-symmetric source, n is the number of higher harmonics in the voltage range; quai, tyuu, tyuci are the phases of corresponding i-th harmonic components of phase voltages; su is the coefficient taking two values 0 or 1 depending on the switch SU position in Fig. 2; ◦ is the operator of element-wise multiplication of elements of vectors and matrices (Hadamard product);
• the 3D complex vector of phase currents of the nonlinear load
v = [l 1 if. (13)
By setting the values of s„, Sj factors and elements of vectors (8) - (12) it is possible to simulate any of 279 modes of operation of the three-phase ESS and calculate the components of the total power losses. To summarize the results of the calculations we assume that the values of the elements of the vectors in (8), (11) are linked by a common regularity
kl + kl + £c2 = 3. (14)
The coefficient of the phase A ka remains unchanged and equal to 1, then
ka = const = 1
kb = var = 0...V2"
(15)
kc 2 - kl =V2...0
We assume that the values of the elements of the vectors in (9), (10), (12) subject to other regularity kb = var = 0...V!
ka = kc = kb
(16)
where J is the RMS of the i-th harmonic component of the phase current of the three-phase nonlinear load; 9iai, tyibt, 9ici are the phases of corresponding i-th harmonic components of phase currents; Sj is the
Regularities of change of vector elements values
(15) and (16) will be used later to summarize the simulation results for different operation modes of three-phase ESS. Note that these ratios may vary in any arbitrary manner, and the selected regularities (15) and
(16) represent the most characteristic changes influences of the three-phase ESS parameters on the total losses.
s
ns
Fulfillment of the condition Pusf = const at the three-phase energy supply system modeling.
Model representation in the vector form. Changes in the parameters of the elements of the equivalent circuit of three-phase ESS changes the value of the average active useful load power. Returning its initial value, as previously mentioned, is carried out by varying the coefficient kh Obtaining the analytical expressions for determination of kl is a difficult problem. It is easier to use the software for calculations. For automated calculation of the coefficient kl it is necessary to build the equation for the useful load power as shown in Fig. 3. There are two forms of this equation - in the vector and matrix form. We consider the first method.
We use the method of the current superposition. In the circuit of Fig. 3 we leave the voltage source. We define the vector of three-phase ESS conductivities for the i-th current harmonic
y, =-
1
¿Li
+Rs -v
(17)
The total conductivity of the three-phase four-wire ESS for the i-th current harmonic
ys, = y, - v+ R„l
(18)
We write the dependence for the neutral voltage displacement for the main harmonic of the three-phase voltage source
U 01 = y1 • ysl
•(us 0 ku ).
(19)
Vector of partial currents created by main voltage harmonic of the linear source in the ESS linear wires
ie1 = (us " U01 -V)o y1
(20)
If the higher harmonics (coefficient su = 1) are in the input voltage curve, the procedure for calculating the vectors of the partial currents created by the three-phase voltage source must be repeated for each harmonic component using (7), (17) - (20). The total vector of the partial currents created by the three-phase voltage source in the ESS line wires in this mode is described by the relation
1
((e1)2 + Z ( )2
i=2
(21)
„ Z ( J, • Yai • Rs eJ{(ai +9 yai )2
i=2
„ Z ( Jr • Ybi Rs ej(bi +9 ybi ))
i=2
„ Z„ Ji • Yc, • Rs ej(ci +9 yci ))
i=2
o kJ =
Z (J )
i=2
LaLj -bLj
■cLJ
where Yai, Ybi, Yci are the modules of total load phase conductivities for the i-th current harmonics; yyai, qybi, qyci are the initial phases of corresponding conductivities.
We find the vector of partial current created by the three-phase current source in ESS linear wires by the first Kirchhoff law
~ (23)
isJ =J iLJ
We determine the current vector in the linear load circuit
iL =$e )2 + )2 , (24)
and find the voltage vector on the load terminals
UL =
(i e1 0 zL1)2 + Z (iei 0 ZLi )2 + Z (i Lji 0 zLi f -(25) i=2 i=2
Average active useful power
Pusf = Rel
UL "J(e f+((sj f~(j ) = Re| UL )2 " 2 -(J 0 ^lJ)
(26)
Model representation in the matrix form. Vector form for three-phase models of ESS has some disadvantage manifests itself at polyharmonic shape of the supply voltage. With this it is difficult to describe the mutual influence of the higher voltage harmonic components of the three-phase source and higher current harmonics flowing in the non-linear load phases. In such a case it is convenient to use the matrix form of the model description.
Let's set the voltage at the three-phase source terminal by the matrix size «x3
where iei is the partial current vector of the i-th harmonic component in the circuit with three-phase voltage source.
In the circuit of Fig. 3 we leave the current source pre-shorting the voltage source, and define the vector of partial linear load current in the circuit with three-phase current source
Ua1 Ub1
U
c1
Ua„ Ub„ Uc„
(27)
Similarly, we write the matrix of the three-phase current source
j =
Ja2 Jb2
J.
c2
J,
Jb„ Jc
(28)
Separately, we define a matrix of complex load impedances for the three-phase voltage source currents and three-phase current source
u
-au1
-bul
zaun zbun zcun
z
-cul
z„
Zj =
Zaj2 zbj 2 zcj2
(29)
(30)
-ajn —bjn —cjn
We express the complex matrix of phase conductivities for three-phase voltage source currents
Mu
yu = -
Zu + Rs
(3l)
where Mu is the matrix of states size nx3 whose rows are taking the unit values in the case of presence of corresponding harmonic components in the range of the three-phase voltage source.
We write the vector of total complex impedances of the three-phase four-wire system
ys = yu 'V+ Rn"1. (32)
We represent the current in the neutral wire due to the higher harmonic components multiple three of the three-phase current source, as a vector of the length n
'3 j = j
Rs + 3 • Rn + Zj
(33)
The displacement voltage at the neutral at polyharmonic supply voltage and polyharmonic current is represented as a vector of the length n
(34)
- yu 0 U - r u0 = —--v"'3j
ys
Rn •
The line currents matrix in the circuit with three-phase voltage source
le = (u "
«0 • vT Io yu
(35)
Using matrix relations (28) and (30) we express the partial currents of the three-phase current source, respectively, in the load and linear wires
Rs
'Lj = j 0
Rs + Zj
j - 'Lj •
(36)
(37)
We determine total line current and total current in the linear load phases
(38)
'L
*Sj'
= 'e + 'Lj
And find the voltage on the load terminals
uL = ie ° zu + iLj ° zj .
Average active useful load power
I • mU1 ^
Pusf = ( o (( - |j|
(39)
(40)
(41)
Both described models allow to relatively simply calculate the value of the correction factor k returning the original value of the average load useful power even
when using standard MathCad tools (Given block and Find function). After that it is possible to determine components of the power losses and the comparison of calculation results at ESS different modes of operation.
As an example, we consider a three-phase ESS with a ratio of three-phase power of resistive short circuit to load active power ksc = 20. The RMS of the phase sinusoidal voltage of the three-phase symmetric source Us = 220 V. If the useful load power at the mode of minimum of losses equal to PusJ = 400.1 kW, then Rs = 0.01815 Q, R¡ = 0.3256 Q. According to (8) - (12) we consider separately 5 modes of ESS operation, each of which is caused by one factor influencing the occurrence of additional losses:
1. Voltage unbalance.
2. Resistive load unbalance.
3. Symmetric active-inductive load at yL = 15° (L = = 0.2777 mH).
4. Higher harmonics in the supply voltage curve (odd harmonics whose amplitudes vary inversely their serial number Um/n, where n=2i-1, i = 2, 3.. .19).
5. The symmetrical non-linear load resulting in the higher current harmonics (odd harmonics, the amplitudes of which vary inversely their serial number UJ(n-(Rs+R¡)), where n=2i-1, i = 2, 3.19).
Changing the ratio kb by (15) or (16) determine the alteration of the corresponding known power indicator: asymmetry coefficient by reverse sequence K2, the power factor cos^, THD of voltage THDu and current THDi. Fig. 4 illustrates the relationship of these power indicators with regularity of kb factor changes.
Fig. 5 shows the dependences of the correction factor kl on the kb factor which sets the pattern of changes of the asymmetrical operation mode of the three-phase three-wire ESS by (15). The values of the coefficient k¡ at
kb = 0 and at kb = V2 correspond to the three-phase ESS emergency operation mode at breakage of the line wire. The value of the coefficient kl at kb = 1 corresponds to the three-phase symmetrical operation at which the total power losses are equal to the lowest possible power losses.
Fig. 4. Dependence of power indicators vs coefficient kb
zu =
Z
ki
0.2 0.4 0,6 0.8 I 1.2 1,4 1.6 Fig. 5. Dependences kl = f kb) for three-phase ESS unbalanced modes: 1 - three-phase source unbalance;
2 - resistive load unbalance
Fig. 6 shows the dependence of the correction factor ki on the kb factor which sets the pattern of change of power factor and amplitudes of voltage and current higher harmonic of the three-phase four-wire ESS at Rn = Rs by (16).
Additional power losses determination. Deviation of the operation conditions of the three-phase ESS from the conditions in which the power losses correspond to the minimum possible value, leads to additional losses [10]. Developed mathematical models allow us to calculate the value of power the additional losses in any of 279 modes of operation of the three-phase ESS. Fig. 7 shows the dependences of the relative, as a fraction of the average load useful power, power of additional losses on the kb factor for five considered ESS modes of operation.
0.6 0.8 a
0 0.2 0.4 0,6 0.8 1 1.2 1.4 1.6 Fig. 6. Dependences ki = f(kb) at: 1 - active-inductive load; 2 -symmetric nonlinear load; 3 - three-phase voltage source higher harmonics
0.6 0.8
b
Fig. 7. Dependence of the relative additional power losses on the factor kb
Analysis of the dependences in Fig. 7 shows that additional power losses in three-phase ESS at ksc = 20 can range from a few tenths to a few percent of the average load useful power. The three-phase ESS with non-linear load have the smallest power efficiency. A combination of factors causing additional losses decreases energy efficiency of three-phase ESS.
Analytical data obtained on the developed model (27) - (41) for five of considered cases are presented in Table 1.
Table 1
Calculation of correcting factors and power indicators for five operation modes of the three-phase ESS
kb Mode 1 Mode 2 Mode 3 Mode 4 Mode 5
klu Ku APadd* kla K21 APadd* kir cos^ APadd* klun THDu APadd* kij THDi APadd*
0 0.8028 0.522 0.0137 1.6551 0.522 0.01612 1 1 0 1 0 0 1 0 0
0.1 0.8324 0.462 0.0112 1.4839 0.462 0.01251 0.9993 1 0.00004 1.0021 0.047 0.00005 0.9997 0.041 0.0003
0.2 0.8608 0.407 0.009 1.3576 0.407 0.0096 0.9971 0.999 0.00016 1.0084 0.094 0.0002 0.9989 0.082 0.0011
0.3 0.8877 0.354 0.0071 1.2615 0.354 0.0072 0.9935 0.997 0.00037 1.0188 0.141 0.0004 0.9975 0.123 0.0024
0.4 0.9129 0.304 0.0053 1.1868 0.304 0.0053 0.9883 0.994 0.00066 1.0336 0.188 0.0007 0.9957 0.164 0.0044
0.5 0.9359 0.254 0.0038 1.1282 0.254 0.0037 0.9816 0.991 0.00104 1.0526 0.235 0.00101 0.9932 0.205 0.0068
0.6 0.9564 0.206 0.0025 1.0824 0.206 0.0024 0.9733 0.987 0.00153 1.0759 0.282 0.00134 0.9902 0.246 0.0099
0.7 0.9738 0.157 0.0015 1.0473 0.157 0.0014 0.9634 0.983 0.00212 1.1037 0.329 0.00164 0.9866 0.287 0.0135
0.8 0.9875 0.107 0.0007 1.0218 0.107 0.00064 0.9516 0.978 0.00284 1.1359 0.367 0.0019 0.9825 0.328 0.0178
0.9 0.9966 0.055 0.0002 1.0058 0.055 0.00017 0.9378 0.972 0.0037 1.1726 0.423 0.0021 0.9778 0.368 0.0227
1.0 1 0 0 1 0 0 0.9219 0.966 0.0047 1.214 0.47 0.00213 0.9725 0.409 0.0283
1.1 0.9958 0.061 0.00023 1.0071 0.061 0.00021 0.9035 0.959 0.00595 1.26 0.517 0.0021 0.9665 0.45 0.0345
1.2 0.9808 0.133 0.0011 1.0338 0.133 0.001 0.8825 0.952 0.00743 1.3109 0.564 0.0019 0.96 0.491 0.0415
1.3 0.9479 0.227 0.0031 1.1008 0.227 0.0029 0.8581 0.944 0.00922 1.3665 0.611 0.00156 0.9525 0.532 0.0494
1.35 0.918 0.293 0.0049 1.1731 0.293 0.0049 0.8445 0.94 0.01027 1.3962 0.634 0.00135 0.949 0.553 0.0536
1.38 0.8901 0.35 0.0069 1.2538 0.35 0.007 0.8358 0.938 0.01095 1.4145 0.648 0.0012 0.9466 0.565 0.0562
1.4 0.8608 0.407 0.009 1.3576 0.407 0.0096 0.8298 0.936 0.01144 1.427 0.657 0.0011 0.945 0.573 0.058
1.414 0.8028 0.522 0.0137 1.6551 0.522 0.0161 0.8254 0.935 0.0118 1.436 0.664 0.00101 0.9438 0.579 0.0593
Three-phase ESS Matlab-model. In accordance with the equivalent circuit of Fig. 2 and matrix mathematical model (27) - (41) a Matlab-model of the three-phase ESS shown in Fig. 8 was created. It differs from the previously developed models [5, 6, 11] by the presence of optimized load structure which is divided into a symmetric resistive load and a adjustable symmetrical three-phase current source and a block of automated calculation of the correction factor k (kl Calculation). The model allows at setting the values of the elements of the vectors (8) - (12) without an operator to automatically calculate the coefficient k and calculate the values of the total power losses components taking into account the calculated ratio.
The Matlab-model consists of:
1. the power circuit configured by Fig. 2 and included:
• the three-phase source of symmetric sinusoidal voltages Usa, Usb, Usc;
• the three-phase source of non-sinusoidal voltages Ua, Ub, Uc;
• the symmetrical resistive load Rla, Rlb, Rlc;
• the three-phase current source Ja, Jb, Jc;
• the active resistances of the four-wire line Ra, Rb, Rc, Rn;
• the power active filter configured on three adjustable current sources SAF.
2. the current and voltage sensors connected to the model's power circuit.
3. the subsystems of the measured data processing:
• the subsystem of the total power losses components calculation by (1) Calculation 1 (Fig. 9);
• the Subsystem of the total power losses components
calculation in p-q-r coordinates Calculation 2 (Fig. 10);
• the subsystem of load useful power harmonic composition calculation Calculation 3;
• the subsystem of the total power losses components calculation by coefficients of direct, reverse and zero consequences Calculation 4 (Fig. 11);
• the subsystem of calculation and indication of measured information on instantaneous values of electrical quantities and three-phase ESS parameters Measurements 1;
• the subsystem of calculation and indication of measured information on average and RMS values of electrical quantities and three-phase ESS parameters Measurements 2 (Fig. 12).
4. the subsystems of the ESS state control and monitoring:
• the PAF control system SAF Controller (Fig. 13);
• the subsystem for setting control signals for the three-phase current source and three-phase non-sinusoidal voltage source Control System;
• the subsystem of the correcting factor kl calculation kl Calculation;
• the block for the PAF connection to the ESS PAF Connection.
5. the virtual measured devices.
The developed Matlab-model can be used to analyze the three-phase energy supply systems and for assessment the energy efficiency at to a power system of the active power filter at any of 288 possible modes of operation.
, -H I-Kiçç
Fig. 10. Subsystem of the total power losses components calculation in p-q-r coordinates Calculation 2
| qr )—► RMS
| qq >-► RMS
sqrt(u(1)*u(4))
-1 iq10
sqrt((u(1 )+u(3))/(sqrt(u(1 )*u(4)))-1 ) h
sqrt(u(2)/(sqrt(u(1 )*u(4)))*(1 +3*Rn/Rs))
iq5
u(1 )*(1 +u(2)A2+u(3)A2*(1 +3*Rn/Rs))
dP1*
►Œ)
1—1 km*
-KD
K2I*
Fig. 11. Subsystem of the total power losses components calculation by coefficients of direct, reverse and zero consequences Calculation 4
Fig. 12. Subsystem of calculation and indication of measured information on average and RMS values of electrical quantities and three-phase ESS parameters Measurements 2
Switch
Fig. 13. PAF control system SAF Controller
Conclusions.
1. We proposed two ways to create a mathematical model of the three-phase ESS, the first one is based on the vector form, and the second one is based on the matrix representation of power processes in multiphase systems. It allows solving the problem of maintaining of constant values of average active load power in the ESS operation at different modes.
2. Based on the analysis of the results of modeling of five operating modes of the three-phase ESS, at each of which a unique factor that causes the appearance of additional losses in the system is involved, it was concluded that the maximum reduction of power efficiency corresponds the ESS with non-linear load.
3. Using the developed mathematical models the operation of the MATLAB-model of the three-phase ESS with PAF is optimized. The possibility of the automated calculation of the correction factor k¡ is realized which is more than an order reduced the working time with model at investigation of the ESS operating modes at which additional losses arise
REFERENCES
1. Akagi H., Kanazawa Y., Nabae A. Instantaneous reactive power compensators comprising switching devices without energy storage components. IEEE Transactions on Industry Applications, 1984, vol.IA-20, no.3, pp. 625-630. doi: 10.1109/TIA.1984.4504460.
2. Wong M.-C., Dai N.-Y., Lam C.-S. «Active Power Filters», Parallel Power Electronics Filters in Three-Phase Four-Wire Systems. Springer Singapore, 2016. pp. 59-165. doi: 10.1007/978-981-10-1530-4_3.
3. Artemenko M.Yu., Batrak L.M., Mykhalskyi V.M., Polishchuk S.Y. Analysis of possibility to increase the efficiency of three-phase four-wire power system by means of shunt active filter. Tekhnichna elektrodynamika, 2015, no.6, pp. 12-18. (Ukr).
4. Artemenko M.Yu., Batrak L.M., Polishchuk S.Y., Mykhalskyi V.M., Shapoval I.A. The effect of load power factor on the efficiency of three-phase four-wire power system with shunt active filter. 2016 IEEE 36th International Conference on Electronics and Nanotechnology (ELNANO), IEEE, 2016. pp. 277-282. doi: 10.1109/ELNAN0.2016.7493067.
5. Zhemerov G.G., Tugay D.V. Physical meaning of the «reactive power» concept applied to three-phase energy supply systems with non-linear load. Electrical engineering & electromechanics, 2015, no.6, pp. 36-42. (Rus). doi: 10.20998/2074-272X.2015.6.06.
6. Zhemerov G.G., Tugay D.V. An universal formula clarification to determine the power losses in the three-phase energy supply systems. Bulletin of NTU«KhPI», 2015, no.12, pp. 339-343. (Rus).
7. Peng F.Z., Ott G.W., Adams D.J. Harmonic and reactive power compensation based on the generalized instantaneous reactive power theory for three-phase four-wire systems. IEEE Transactions on Power Electronics, 1998, vol.13, no.6, pp. 1174-1181. doi: 10.1109/63.728344.
8. Afonso J., Couto C., Martins J. Active filters with control based on p-q theory. IEEE Industrial Electronics Society Newsletter, 2000, vol.47, no.3, pp. 5-10.
9. Kim H.S., Akagi H. The instantaneous power theory on the rotating p-q-r reference frames. Proceedings of the IEEE 1999 International Conference on Power Electronics and Drive Systems. PEDS'99 (Cat. No.99TH8475), 1999, pp. 422-427. doi: 10.1109/PEDS.1999.794600.
10. G. Zhemerov, N. Ilina, D. Tugay. The Theorem of Minimum Energy Losses in Three-Phase Four-Wire Energy Supply System. 2016 2nd IEEE International Conference on
Intelligent Energy and Power Systems (IEPS-2016). June 07-11, 2016, Kyiv, Ukraine, pp. 52-54. doi: 10.1109/IEPS.2016.7521889.
11. Zhemerov G.G., Tugay D.V. Components of total electric energy losses power in pqr spatial coordinates. Electrical engineering & electromechanics, 2016, no.2, pp. 11-19. (Rus). doi: 10.20998/2074-272X.2016.2.02.
Received 01.06.2016
D.V. Tugay, Candidate of Technical Science, Associate Professor,
O.M. Beketov National University of Urban Economy in Kharkiv,
12, Revolution Str., Kharkiv, 61002, Ukraine. phone +38 057 7073111, e-mail: [email protected]
How to cite this article:
Tugay D.V. Three-phase energy supply systems simulation for the total power losses components assessment. Electrical engineering & electromechanics, 2016, no.4, pp. 43-53. doi: 10.20998/2074-272X.2016.4.06.
