CC BY
10Digital Control of Electric Power Flows in Unbalanced Distribution Networks as Part of The Automated Metering and Control System
T.T.Omorov1*' B.K.Takyrbashev1 K.E.Zakiriaev2, T.Zh.Koibagarov1
1 National Academy of Sciences of the Kyrgyz Republic, Institute of Mechanical Engineering and Automation, Bishkek
2 Issyk-Kul State University named after K. Tynystanov, Karakol, Kyrgyz Republic
Abstract — This paper aims to address the problem of controlling the electricity flows in power distribution networks (PDN) operating under current and voltage unbalance. As is known, the unbalance factor is responsible for significant losses of active power and, therefore, is detrimental to the PDN efficiency and technical and economic performance. The purpose of control is to minimize technical power losses in the distribution network. This is to be achieved by building an information control system (ICS) for balancing a three-phase network as part of the automated metering and control system (AMCS). The latter is currently being widely adopted to automate information processes in PDNs. However, the AMCS does not include technologies for solving the problem in question. We propose an algorithm of the digital controller operation for the ICS. Its primary function is to maintain phase power at a given level in real-time. The algorithm concept is based on the idea of required redistribution of electricity flows between the phases of the distribution network by appropriately switching single-phase loads of consumers (users) to ensure a minimum spread of phase powers relative to their required level. To achieve the goal of control, we construct criterion functions that determine qualitative indices of the ICS operation and develop computational schemes for their minimization. Control actions to be generated by the digital controller and performed on the facility represent a digital code that contains data on the coordinates of three-phase network loads to be switched to another phase.
Index Terms, distribution network, unbalanced mode, control of electricity flows, digital controller, control algorithm.
* Corresponding author. E-mail: omorovtt@mail.ru
http://dx.doi.org/10.38028/esr.2021.01.0005 Received March 26, 2021. Revised April 03, 2021. Accepted April 13, 2021. Available online May 26, 2021.
This is an open access article under a Creative Commons Attribution-NonCommercial 4.0 International License.
© 2021 ESI SB RAS and authors. All rights reserved.
I. Introduction
At present, power distribution networks (PDN) are widely adopting advanced digital technologies in the form of automated metering and control systems (AMCSs) to automate their processes and make them data-driven [1]. These systems are elements of the Smart Grid technology [2]. An analysis of the functional structure of existing (conventional) AMCSs shows that they mainly perform metering and fiscal metering of electricity and belong to the class of information and measurement systems. The latter lack tools designed to solve the problem of balancing three-phase distribution networks. At the same time, most of the 0.4 kV power distribution networks operate under current and voltage unbalance [3-5], which causes nonoptimal operating conditions. Unbalance of PDN results from unbalance of loads and their uneven distribution over phases of a three-phase network. According to the published research data, in networks with distributed load, which have a deviation (unbalance) of phase currents from their average values, technical losses increase significantly [6-8]. The results of experimental studies of power supply systems indicate that technical losses in power lines and transformer substations (TS) supplying them amount to more than 6% of the total electricity supplied to the inputs of distribution networks [9]. At the same time, the quality of electricity deteriorates, and the likelihood of failure of household appliances and industrial installations also increases [8]. Although many methods and devices were proposed [1013] for balancing distribution networks, their practical use proves difficult. This is due to the technical challenges of creating automated balancing systems on their basis for three-phase distribution networks with a large number of electricity consumers. For the technical implementation of these solutions, at points of load connection to the network phases, it is necessary to create separate local subsystems with individual units as their parts to perform functions of parameter measurement, control, and data processing. However, building such subsystems involves rather elaborate processing algorithms based on measurements of a group of information parameters. For example, the method outlined in [10] requires measurements of the active
power of a zero wire of line, voltage and current signals of the zero wire at the points of load connection to the network, while the reactive power control system [13] uses positive sequence, negative sequence, and zero sequence currents of loads. In modern electricity meters [16], which are part of the AMCS, these measurement functions are not available. In this regard, the use of the above technologies as part of modern AMCSs seems to be problematic. In [14, 15], one of the possible approaches proposed to solve the given problem is based on a digital system for automatic control of the network balancing process as part of the data concentrator of the AMCS. As is known, a data concentrator in conventional automated systems collects data from a group of electricity meters installed at the side of network users and performs computing operations related to fiscal metering of electricity. According to this approach, the data concentrator is assigned with an additional function of digital control, which provides the required redistribution of electricity flows between the network phases by switching single-phase loads from a more loaded phase to a less loaded one. At the same time, electricity meters must have dedicated actuating elements, i.e., phase current switches (PCS), as their part. The creation of such PCSs is a feasible task since modern single-phase and three-phase meters [16, 17] have relays with powerful contacts up to 100 A, which can be used to switch power consumers from one phase to another [15].
The present study aims to develop this approach for solving problems of control of electricity flows in unbalanced distribution networks.
II. Problem statement Let us consider a three-phase four-wire distribution network to which a power supply transformer substation (TS) and a group of n electricity meters {Emvk} installed at singlephase consumers are connected. Its simplified diagram is shown in Fig. 1, where k, v are index variables denoting, respectively, the numbers of phases A, B, C (k = 1,3) and loads (electrical consumers) of network users (v = 1, n ); Ik, Uk are sinusoidal instantaneous current and voltage at the input of the k-th phase; Zvk is the load with the coordinate (v, k); Ivh Uvk are instantaneous current and voltage across the load Zvk; J is instantaneous current in the zero wire; Pk is complex power of the k-th phase; Pvk is complex power consumed by the consumer load with the coordinate (v, k); PCSvk is a switch of phase currents of the meter Emvk. Next, we make the following assumptions:
1. the three-phase network operates under current and voltage unbalance;
2. the phase current switches (PCSvk) are a part of the electricity meters (Emvk);
3. the concentrator at discrete moments of time t = ^ (4 = 1, 2, ...) periodically polls the head-end three-phase electricity meter (EmTS) and groups of consumers' meters connected to phases A, B, C, and writes to the database AMCS values of the following quantities:
• active prk and reactive qrk power consumed by loads of network users with coordinates (v, A);
• active pk and reactive qk powers consumed by the corresponding phases at the input of the three-phase network.
TS
X Em22
PCS22
Z2 2 P2 2
¡22 U22
Fig. 1. Simplified diagram of the distribution network with Emvk and PCSvk
Technical losses of electricity in the transformer substation and the network are largely determined by the value of current J flowing at the head-end of the zero wire. Moreover, the lower the value of this current, the higher the level of the network balance and the lower the power loss there. Therefore, the quality and efficiency of the distribution network can be assessed by the following objective function [15]:
E = J. (1)
Minimizing the efficiency index of system E is equivalent to optimizing the operating condition of an unbalanced distribution network, which provides a minimum unbalance in the values of the total power consumed by each phase of the network. To this end, we propose creating an appropriate information control system (ICS) as part of the modern (conventional) AMCS. The development of the ICS is to be based on the principles of building automatic systems [18] and digital technologies, the functional structure of which is shown in Fig. 2. In this case, the group of loads {Zvk} of network consumers to which the electricity meters {Emvk} are connected, together with the
ICS actuating elements, is a control object. The key functional block of the system is a digital controller (DC) based on a microprocessor controller [19, 20]. The DC block generates control actions u to be performed on the facility based on a dedicated algorithm (control rule). The DC operation program is built by the automatic system initial state identifier (ISI) as setting action p*. In this case, the control signal u* is a digital code of the command, which is formed as vector F = [Fj, F2, p], where F1, F2 are numbers (names) of a pair of phases in which it is necessary to switch network consumers with a more loaded phase (/•',) to a less loaded one (I' -,): p is a vector composed of addresses (numbers) of users of phase /•', to be switched. This control signal u is transmitted via the communication channel (CC) to the system actuating elements, which are used as phase current switches (PCS) designed to perform the required switching of the network users' loads from
one phase to another [15]. Conventional AMCSs employ various data transmission technologies (PLC, GSM, and others) as the CC.
The task is to synthesize an algorithm for the operation of a digital controller that optimizes the functioning of an unbalanced distribution network.
III. The control problem
The above control problem-solving process includes the following main steps:
1. Define the input data of the control problem.
2. Conduct a case study of the facility.
3. Construct criterion functions.
4. Build an algorithm for minimizing criterion functions.
5. Generate control actions to be performed on the facility.
A. Define the input data of the problem
Based on the data obtained from the main three-phase meter and consumer meters, the complex powers Pk and Pvk, respectively, at the network input and consumer loads, are represented as [21]:
Pk = Pk + Mk, k =1,3 ,
Pvk = Pvk + M*, v = 1,n , k = 1,3 , j =V-T,
Where j is an imaginary number. In this case, absolute values of these complex quantities are determined by the expressions:
Pk = Jfl + j , k = 1,3, (2)
P = Jplk + q2 , v = M ,k = U. (3)
It is worth noting that to solve the above-stated control problem, it is sufficient to poll only a small group of single-phase consumer electricity meters, except for meters installed at two-phase consumers. Further, we denote by ;?2, and ;?3 the total number of simplified consumer meters connected to phases .4, 5, and C, respectively. For brevity of notation and convenience of further transformations, the data on absolute values of complex powers of network consumers, determined as per
Controlled facility
Emrs
{Emvk}
CC *
Fig. 2. ICS functional structure.
formulas (3), are provisionally presented in the form of the following vectors:
pA = [P, p,..., pa ],
P* = [P*, PP* ],
pc = [ pc, p2c p ],
and the coordinates of their elements are set by vectors:
pf = [pf ,pf ,...,pf],
pf = [pf ,pf ,...,pf ], pC =[pf ,PC.....Pi; ],
where superscripts A, B, and C indicate that the introduced vectors and their elements refer to the corresponding phases of the three-phase network. When forming vectors PA, PB, PC, their components can be positioned arbitrarily. In this case, the set of components of indicated vectors make up the sets Aa, Ab, Ac, respectively:
A, = {Pt,' = Kt , Ac = {PL,
B. Conduct the case study of the facility Such a study is performed to determine the structure of phase switching, i.e., to identify names of phases (Fl, F2), in which corresponding loads of network consumers are to be switched. To this end, the average value P0 of the phase power absolute values is provisionally calculated based on the formulas (2):
Po = P + P2 + Ps)/3. (4)
Similarly, one can determine the average value of effective phase currents I1, I2, and I3 at the network input:
Io = I + I2 + Is)/3.
To assess the level of unbalance (deviations) of phase powers and currents from their average values, we introduce the following differences (residuals) into consideration:
ek = Pk - Po, k = 1^3,
— (5)
Ek = Ik - I0, k = 1, 3 .
Note that the values e;, e2, and e3, as determined by formulas (5), are essentially control errors (residuals) that should tend to zero during control. In this case, current absolute values of phase powers at the input of the distribution network are determined by the formulas Pk = Po + ek, k = U.
Hence one can_ see that at ek ^ 0 phase powers are Pk ^ P0, k = 1,3 . Thus, the problem of optimizing the operating condition of the distribution network in each interval [t^, t^+1] is reduced to minimizing control errors e1, e2, and e3. In this case, control actions u , when performed on the facility, also lead to desired changes in effective phase currents I1, I2, and I3 at the network inputs. Indeed, by definition, complex powers Pk at inputs of corresponding phases are determined by the formulas [21]:
Pk = Ujl = I2k zk, k = U,
where Uk, Ik are complex voltage and current at the input of the k-th phase; I* is a complex conjugate value corresponding to current ik; zk is complex resistance of the k-th phase. Therefore, the following expressions can be written for absolute values:
P
12 Z
k = 1,3
where Zk is the absolute value of complex impedance zk. Therefore, for effective phase currents Ik, the following formulas are valid:
a)
ek
b)
ek
i 0 i e i 1 2 e, ^ t- i ei e2 1 3 L ei 2 3
1 Ï2 3 0 1 2 e: 1 e ej
Fig. 3. Illustration of possible variants of control errors et
I„
k = 1,3
It is worth noting here that the connection of a new load to the selected phase leads to an increase in the magnitude of the absolute value of power Pk and a decrease in its impedance Zk. Therefore, in this case, the ratio Pk /Zk increases, and when the load is disconnected, the value of Pk /Zk decreases. Thus, in the process of control, an increase or decrease in the absolute values of power Pj, P2, and P3 leads to the same changes in the dynamics of operating phase currents Ij, I2, and I3.
In this case, corresponding absolute values of residuals for phase currents |ej|, |e2|, and |e3| decrease so that the minimum of their unbalance is ensured. As a result, along with the equalization (balancing) of phases in terms of power, control of electricity flows in the distribution network also makes phase currents get balanced, which leads to the minimization of effective current J flowing in the head-end of the zero wire of the three-phase network.
Based on the results of calculating the residuals ek, k = 1,3 , the following cases can be distinguished:
Case 1. One of the components of the control error vector has a positive value e = [ej, e2, e3], one of them is negative, and one of them is equal to zero. For example, ej > 0, e2 < 0, e3 = 0, which is shown in Fig. 3a. In this case, phase A is overloaded, therefore, it is necessary to switch its consumers to phase B, which we denote further by A ^ B (Fj = A, F2 = B). In this case, phase C should retain the previous state since e3 = 0.
Case 2. Two components of the residual vector e are positive, and the third one has a negative value. For example, ej > 0, e2 > 0, e3 < 0 (Fig. 3b). In this case, phases A and B are overloaded. It is, therefore, necessary to switch A ^ C and B ^ C.
Case 3. One of the elements of vector e is positive, and the rest are negative. For example, ej > 0, e2 < 0, e3 < 0 (Fig. 3c). In this case, it is necessary to switch A ^ B and A ^ C.
Thus, to identify the required switchings of the distribution network phases, it is necessary to perform a case study of the signs of values of ek, k = J, 3 . In this case, it is necessary to switch some consumers from phases with positive values of errors ek to those with negative signs of control errors.
As a result, the case study of the facility provides for the following operations:
J. Poll the head-end three-phase electricity meter (EmTS) and a given set of meters installed at network consumers connected to phases A, B, C, and define input data of the control problem.
2. Form vector P = [PJ, P2, P3], consisting of the absolute values of complex powers consumed by phases A, B, and C, according to formulas (2).
3. Calculate average value P0 of the phase power absolute values, as determined by formula (4).
4. Determine values and signs of the components of the vector of control errors e = [eJ, e2, e3] as per formulas (5).
5. Identify names of phases (FJ, F2), in which load switching operations are necessary, based on the analysis of signs of control error vector components
e [ej, e2, e3].
C. Construct the criterion function
It is easy to show that the initial state of the facility as it manifests itself by the components of the control error vector e = [eJ, e2, e3] satisfies the following condition:
eJ + e2 + e3 = 0. (6)
Next, without loss of generality, consider Case 1 at the facility. In the case under consideration, the synthesized controller must output the initial mismatch p* = p* through the control error, which is formed by the ISI block. Its value can be determined from condition (6). Since following the considered Case the residual e3 = 0, condition (6) will have the form
ei = |e2|.
Therefore, we arrive at p* = |e2| . Next, we form subset consisting of a set of elements of the set Aa, i.e., ^ e Aa. In particular, can consist of only one element. Let ^ denote subset of indices (numbers) of vector components PA: = {1, 2, ..., nj}.
Now let us introduce the following function:
F (X ) = |S (X,) - p* |, (7)
where S1(X1) is a discrete function determined by the sum of elements of subset :
51 () =
leR
The function F1(X1) defined by formula (7) can be used to indirectly estimate the value of the objective function E defined by expression (1). The analysis shows that when F1(X1) ^ 0, the control errors e1 and e2 also tend to zero when switching consumers with more loads connected to the first phase (phase A).
D. The algorithm for minimizing criterion functions
Now it can be seen that in the case under consideration,
coordinates of loads of phase A users can be identified by solving the following extremal problem:
min F (>•, ) = F *), (8)
where X* is the optimal set of absolute values of complex power consumed by phase A loads (k =1) to be switched
to phase B. The set X* found this way contains complete information about the coordinates of consumer loads to be switched.
The minimization problem (8) is discrete and can be solved by enumerating options, which leads to time-consuming computational operations and inadequate efficiency of the procedure for generating control actions to be performed on the facility. The analysis shows that a
simpler and more efficient scheme for solving the extremal problem (8) can be constructed by transforming the input data of the problem in a certain way and constructing criterion functions so that in the process of solving it, a purposeful movement towards the optimum is provided by analogy with the gradient method [22]. For this purpose, we propose a computational procedure. First, we preliminarily transform input data of the problem, represented by power vectors PA, PB, and PC. Such a transformation involves sorting the latter to arrange their components in ascending order. In this case, new vectors PA , PB , and PC with the same dimensions as the original ones are formed:
PA = [pA, Pp,..., pa PB = [P* , p,..., P*
where
p c = [ pf, p2 p
pA <■ pA pB <- pB p c ^ p c rv-l— ry ' -1 — ' n '
v = 2, n , t = 2, n2, n = 2, n3.
In this case, new coordinates of the components of the indicated vectors are formed:
ß r =
ßA,ß?,
,ß n
ß * = [ ßf ,ßf ,... ,ß* ],
ßc =
ßf ,ß C,-,ß C
A simple analysis shows that to define functions F1(X1), depending on a situation at the facility, it is more expedient to use components of new vectors PA , PB , and PC . Now let us consider Case 1 of the facility, residuals are e1 > 0, e2 < 0, and error is e3 = 0. In this case, the criterion function F1(X1) must be specified in the form
F () = F (r ) = ¿P4, (9)
i=1
where r is some integer (r > 1).
As a result, problem (8) will be written as
mm F (r)= F (O . (10)
Thus, the original problem is reduced to determining integer r = r* at which the optimum criterion function F1(X1) is achieved.
By solving (10), we identify vector p* that consists of addresses of loads of phase A users to be switched to phase B. Vector p* found in this way is a truncated part of original vector p4 = [p„p2,..Jj, that is, p* = [p„p2,..„p,.] . In this case, the desired control vector u* has the following form:
u* = [a, b, P*] = [1,2, p„ p2,...,pr.].
The above algorithm for solving the control problem can also be used for cases where the facility operation corresponds to Case 2 and Case 3. Consider the features of constructing algorithms of the functioning of a digital controller for these cases.
Suppose that Case 2 is observed in the controlled facility. Without loss of generality, we assume that the control errors are e1 > 0, e2 > 0, and e3 < 0. In this case, condition (6) will be written in the form
e1 + e2 = |e3|.
In this case, it is necessary to carry out two switchings A ^ C and B ^ C. Therefore, based on the data received from the consumer electricity meters, we will form two vectors of the input data:
pa = [ pa, pa,..., p
pb =
[ PB
pB pB
<r2 m
and two vectors of their coordinates:
pa=[pf ,pA ], pb=[pf ,pf ,...,pm ],
where m1, m2 are dimensions of corresponding vectors.
Next, components of vectors PA and PB are sorted similarly to the previous case, i.e., new vectors PA and PA are formed with their components arranged in ascending order:
pa = [ pxa, pfp
PA = [ p
A pA 2 '
Now, as previously, subsets R2 and R3 are introduced for the indices of vectors pA and pB, respectively, i.e., R2 = {1, 2, ..., m1} and R3 = {1, 2, ..., m2}.
Following condition (6), the levels of mismatches p*2 and p3 for the absolute values of corresponding complex powers of switched phases A and B are set as follows:
p*2 = e, p3 = e2.
Now, to find the required coordinates of network consumer loads, it is necessary to carry out two types of switchings, A ^ C and B ^ C. In this case, the following two criterion functions are constructed:
F2 (^2 )= \S2 (l2 ) - P2| >
f3 () = \s3 () -P*| >
Where X2, X3 are subsets composed of components of vectors PA and PB , respectively, i.e., e Aa, X3 e AB. In this case, functions S2(X2) and S3(X3) are determined using the components of vectors PA and PB similarly to the previous case (Case 1) according to (9):
^2 (* 2 )= ^2 (n )= Yff,
1=1
S3 (^) = S3 (0 )=±p: .
Y=1
Further, to find the desired coordinates of the loads connected to phases A and B, two extremal problems are solved:
min F (n) = F2 (n*), (11)
ruin F3 (0 ) = F3 (0* ) (12)
where n* and 0* are the optimal quantities of absolute values of complex powers consumed by loads of consumers of phases A and B, respectively.
Extremal problems (11) and (12) are solved using the above algorithm for minimizing the criterion function F1(X1). Since in the considered case it is necessary to carry out two types of switchings A ^ C and B ^ C, then solving the last extremal problems, we obtain vector u* = [u1, u2], where u. = [A, C, p*], u2 = [B, C, P*2],
P* = [pf, PA,...,pJA], p; = [ Pf, Pf,..Jf.].
When Case 3 is observed at the facility and, for example, errors are e1 > 0, e2 < 0, e3 < 0, condition (6) takes the form
e1 = N + N.
In this case, it is necessary to carry out switchings A ^ B and A ^ C. Therefore, two vectors are formed from the set Aa:
pa = [pa,papa ],
pA=[ P A, P A p;A ],
as well as vectors PA and (3A, respectively, for their coordinates. Moreover, vectors P4 and PA must have different components. Respectively, the components of these vectors form sets that do not intersect in this case. Next, we form the subsets e Aa and X5 e Aa. These input data are used in the above computational scheme to determine the desired control vector u*. Firstly, by sorting components of vectors PA and PA, new vectors PA and 13A are formed with their components arranged in ascending order. Then, criterion functions F4(X4) and F5(X5) are determined:
F4 () = () - P4|' f5 () = () - P5| >
where p*4 = \e21 and p5 = |e31, and the functions S4(X4) and S5(X5) are defined similarly to the previous case. As in previous cases, the desired control vector u* is determined based on minimizing criterion functions F4(X4) and F5(X5). In this case, the algorithm for solving the extremal problem (10) is used.
To improve the ICS reliability and extend the service life of the phase current switch (PCS), it is advisable to limit the number of switchings of consumer electricity meters. To this end, a current sensor installed at the transformer substation can be used, the output of which is fed to a digital controller (DC). Based on the value of current J of zero wire, DC checks the balance (quasi-optimality) condition of the system in each control interval. This condition, in particular, can be specified by the following inequality:
J < A, (13)
where A is the maximum permissible value of current J, at which permissible level of the distribution network
balance is ensured. To limit the number of switchings of consumer meters, control actions are generated and implemented only when condition (13) is not met.
The ICS constructed in this way with a digital controller can be considered as a system of a variable structure [23] since this will ensure the adaptation of the three-phase distribution network to the conditions that arise during uncontrolled random changes in its loads by automatic self-adjustment of its structure in real-time.
Based on the results obtained above, one can formulate a generalized algorithm of the ICS operation, which includes the following basic operations:
1. Poll the main three-phase electricity meter (EmTS) and a given set of single-phase consumer meters connected to phases A, B, and C, except for meters installed at two-phase consumers.
2. Measure effective current J of zero wire and verify condition (13). If it is not met, go to step 3. Otherwise, the DC block does not generate control commands for the facility.
3. Formulate the input data of the control problem determined by vectors P, P4, PB, PC, PA, PB, PC.
4. Estimate components of control error vector e = [e1, e2, e3] by formulas (5).
5. Conduct a case study of the facility based on the signs of components of the control error vector e = [e1, e2, e3] and identify an interphase switching structure determined by vector F = [F1, F2].
6. Determine the initial mismatch p* by the ICS unit and construct corresponding criterion functions {Fi (X)} depending on the identified situation at the facility. For example, for Case 1 at the facility p* = p* = \e2\
, and F1(X1) is constructed as a criterion function.
7. Minimize selected criterion functions {F,- (X)}. In this case, the DC block solves the extremal problems (10) or (11), (12), depending on the situation at the facility.
8. Generate digital control commands u and transfer them to corresponding phase current switches to implement the required switching of the selected loads (meters) of the distribution network users.
IV. Conclusion
Most three-phase distribution networks operate under unbalanced currents and voltages due to the unbalance of loads and their uneven distribution over phases. The unbalance factor causes an increase in active power losses in networks and transformer substations, which significantly diminishes the technical and economic performance of distribution companies.
At the same time, modern AMCS systems do not solve the problem of electricity flow control in unbalanced distribution networks as they focus on optimization of their operating conditions. We have proposed a novel approach to the construction of an information control system (ICS) aimed at solving this problem as part of a conventional AMCS system. The main idea, in this case, is the distribution
of electricity flows between the phases of the network by switching its consumers from a more loaded phase to a less loaded one. To assess the quality of control processes, we have introduced criterion functions that characterize technical losses of electricity in the transformer substation and distribution network. We have proposed computational schemes for constructing them in terms of the situation at the facility. We have developed algorithms of operation of ICS of digital controller and minimization of criterion functions, which ensure minimization of unbalance of phase powers and currents at the network input, and that of technical electricity losses in network and supply transformer substation. The obtained findings are aimed at creating more efficient information control systems by improving conventional AMCS systems that belong to the category of information and measurement systems.
References
[1] Yakushkov K.V. "Automated system of the fiscal electricity metering for the retail market," Informatizatsiya i sistemy upravleniya v promyshlennosti, 2009, no. 3 (23), pp. 9-13 (in Russian).
[2] El-Hawary, M.E. "The Smart Grid—State-of-the-art and Future Trends," Electr. Power Compon. Syst.
2014, 42, pp. 239-250.
[3] Kosoukhob F.D., Vasilev N.V, Filipov A.O. "Reducing losses from current unbalance and improving electric-energy quality in 0.38-kV networks with municipal loads", Russian electrical engineering, Allerton Press, Inc. (New York), vol. 85, no. 6, 2014, pp. 350-353.
[4] Omorov T.T. "AMCS-based assessment of the current and voltage unbalance effect on power losses in the distribution grid," Elektrichestvo, 2017, no. 9, pp. 1723 (in Russian).
[5] Czarnecki L.S., Bhattarai P.D. "Currents' physical components (CPC) in three-phase systems with asymmetrical voltage," Przegl^d Elektrotechniczny,
2015, no.6, pp. 40-47. DOI: 10.15199/48.2015.06.06.
[6] Ponomarenko O. I., Holodinov I.H. "Effect of unbalanced conditions on power losses in power grids of distributed power supply systems," Energetik, 2015, no.12, pp. 6-8 (in Russian).
[7] Mahmoud, K.; Yorino, N.; Ahmed, A. "Optimal Distributed Generation Allocation in Distribution Systems for Loss Minimization," IEEE Trans. Power Syst. 2016, 31, pp. 960-969.
[8] Harutyunyan A.G. "On the calculation of additional power losses in three-phase four-wire grids," Elektrichestvo, 2015, no.10, pp. 55-58 (in Russian).
[9] Averbukh M.A., Zhilin E.V. On the loss of electricity in power supply systems for individual housing construction. Energetik. 2016. no. 6. pp. 54-56 (in Russian).
[10] Samokish VV. Patent No. 2548656 (Russian Federation). "Method of balancing phase currents of three-phase four-wire line and device for
its implementation," Bull. no.11. 20.04.2015 (in Russian).
[11] G.A.Bolshanin. "Patent No. 2249286 (Russian Federation). Method for computer-aided online measurement of current and voltage unbalance level, Bull. no.9. 28.03.2005. (in Russian)
[12] I.VNaumov, D.A.Ivanov, S.V.Podjachikh. Patent No. 2490768 (Russian Federation). Balancer for three-phase networks with neutral wire. Bull. no. 23. 21.08.2013 (in Russian).
[13] Kiselyov M. G., Lepanov M.G. "Balancing of currents in electrical grids by a power electronic controller of non-active power." Elektrotekhnika, 2018, no.11, pp. 63-70 (in Russian).
[14] Omorov T.T. "Balancing of a power distribution network by digital control," Mekhatronika, Avtomatizatsiya, Upravlenie, 2018, vol. 19, no.3, pp. 194-200 (In Russian).
[15] Omorov T.T. "Synthesis of a digital controller for balancing the power distribution network," Instruments and Systems: Monitoring, Control, and Diagnostics, 2019, no. 11, pp.51-56 (in Russian).
[16] https://mir-omsk.ru/products/oborudovanie-dlya-tsifrovykh-setey/umnye-schetchiki-elektroenergii/ .
[17] Tanveer A. "Non-technical loss analysis and prevention using smart meters," Renewable and Sustainable Energy Reviews. 2017. Vol. 72. pp. 573-589.
[18] "Methods of the classical and modern theory of automatic control: in 5 volumes", Ed. by K.A. Pupkova, N. D. Egupova. M.: Publishing House of the Bauman Moscow State Technical University, 2004. (in Russian)
[19] Yiu, Joseph. The Definitive Guide to the Arm Cortex-M3 / Joseph Yiu; translated by A. V. Evstifeev. Moscow: Dodeka XXI, 2012. 552 p.
[20] Developing a New Generation of Smart Automation with ARM Cortex-M3 (https://www.radiolocman.com/ review/article.html?di=134036).
[21] Demirchyan K.S., Nejman L.R, Korovkin N.V et al. Theoretical foundations of electrical engineering, vol.1. St.Petersburg: Piter, 2009, 512 p. (in Russian).
[22] Aoki M. Introduction to optimization methods. Moscow: Nauka, 1977, 334 p. (in Russian).
[23] Emeljanov S.I., Utkin V I., Taran V. A. The theory of systems of a variable structure. Moscow: Nauka, 1970, 592 p. (in Russian).
Turatbek Omorov is D.Sc. in engineering, Head of the Laboratory "Adaptive and Intelligent Systems" of the National Academy of Sciences of the Kyrgyz Republic (NAS KR), Bishkek. He is a Corresponding Member of the NAS KR. He graduated from the Faculty of Automation and Electrical Engineering of Leningrad Electrotechnical Institute in 1975. T. Omorov received his Ph.D. degree in engineering from the Bauman Moscow State Technical University in 1981 and a D.Sc. degree in engineering from the Scientific and Production Association "Cybernetics" of the Academy of Sciences of the Republic of Uzbekistan in 1997. He worked as the Director of the Institute of Automation (20002008), and Vice-President of the National Academy of Sciences of the Kyrgyz Republic (2008-2013). Key research areas are automatic control, informatization and optimization of control processes, and automation of power systems.
Taalaibek Koibagarov is a Ph.D. student at NAS KR, Bishkek. He graduated from the Faculty of Energy of the Kyrgyz Technical University in 2002. He is the Head of the Metrological Service of JSC "Severelectro." Key research areas are automatic control and automation of power systems.
Takyrbashev Beishenaly is a Ph.D. in engineering, a Senior Researcher at the "Adaptive and Intelligent Systems" Laboratory at the NAS KR, Bishkek. He graduated from Frunze Polytechnic Institute in 1973 with a major in automation and remote control. He was the Head of the Power Systems Protection Service at JSC «Severelectro».
Key research areas are automatic control and automation of power systems.
Kubanychbek Zakiriaev is a senior lecturer at Issyk-Kul State University, Karakol, Kyrgyz Republic. He graduated from the Faculty of Automation and Computer Engineering of Tomsk Polytechnic University in 1992. Key research areas are information technology and automation of power systems
