Optimal Voltage Control in MV Network with Distributed Generation

: The article presents the concept of voltage control in a medium-voltage network, using the classical control of a HV/MV (High Voltage/Medium Voltage) transformer and active participation of distributed generation sources. The proposed solution is based on the results of the optimization process. The objective function is considered as a single-criterion—the voltage quality indicator or value of power losses in the network, and optionally as two-criteria (voltage quality and losses combined, with appropriately selected weight factors). The analysis carried out for random selection of independent variables, using the original heuristic algorithm, indicated a very high efﬁciency of the proposed control process, compared to the traditional approach. A signiﬁcant improvement in the voltage quality index and reduction of losses was found, which justiﬁes the advisability of looking for new solutions in the ﬁeld of voltage control in MV taking into account active participation of distributed generation.


Introduction
Power networks are becoming increasingly saturated with distributed sources. Most often, they are wind and photovoltaic farms; in addition, the number of biogas plants is increasing. A significant number of these sources, especially in medium-and low-voltage networks, contribute to voltage problems. Grids, which until now, have worked as receiving networks, periodically become supply networks. Their basic supply from HV/MV transformers becomes supported by sources connected inside the MV network. A number of works have been created, in which problems resulting from a change in the characteristics of distribution networks, so far considered as typical (radial system of operation, unidirectional power flow), are considered in [1][2][3][4]. A change in the direction of power flow (in high-generation states-flow toward the transformer; in maximum load states-flow toward the grid) causes significant changes in the voltage values, which adversely affect the operating conditions of receiving devices. An overview of research works related to voltage problems in the distribution network can also be found in articles [5][6][7].
In line with the traditional approach, in MV distribution networks, the voltage is controlled by the HV/MV transformer equipped with an on-load tap changer (OLTC). The controlling operation maintains the voltage according to the preset value, selected in various ways. The simplest one so far has been to set it to 1.05 U nMV , which provided a voltage in the depths of the network equal to or slightly lower than the nominal one. The presence of distributed generation forces a change in approach to this problem. Due to the transmission of power toward the buses of the power supply station, the voltage in distant parts of the network may significantly exceed the nominal value. In this situation, it is natural to move toward a completely new concept of voltage control that takes into account the active contribution of distributed sources [1,6,8]. In order to realize this concept, in addition to the on-load tap changer, the possibility of data transmission from individual It is obvious that the proposed voltage control system is much more complicated than the traditional voltage controller that, using OLTC, maintains the set voltage value at a selected bus of the network. The aim of the article is to present the most serious problems related to the creation of such an innovative system. These problems are presented in the following sections: • search for an indicator of the quality of control (Section 2) • defining the objective function and task constraints (Section 3.1) • description of the optimization algorithm and justification for the choice of the heuristic method (Section 3. 2) The effectiveness of the proposed solution was confirmed by means of many computational cycles for the test network and randomly selected independent variables. The results are presented in the following sections: • MV test network and research methodology (Section 4) Figure 1. Illustration of the concept of the optimal voltage control system in an MV network with distributed generation (P l1 , P l2 . . . . P lj -loads; P g1 , P g2 , P g3 , P g4 -active power generation; Q g1 , Q g2 , Q g3 , Q g4 -reactive power generation/consumption by renewable energy sources (RES); U HV -high voltage; U MVs -medium voltage; ϑ Tr -HV/MV transformer ratio).
It is obvious that the proposed voltage control system is much more complicated than the traditional voltage controller that, using OLTC, maintains the set voltage value at a selected bus of the network. The aim of the article is to present the most serious problems related to the creation of such an innovative system. These problems are presented in the following sections: • search for an indicator of the quality of control (Section 2); • defining the objective function and task constraints (Section 3.1); • description of the optimization algorithm and justification for the choice of the heuristic method (Section 3.2).
The effectiveness of the proposed solution was confirmed by means of many computational cycles for the test network and randomly selected independent variables. The results are presented in the following sections: • MV test network and research methodology (Section 4); • results for three optional objective functions, finding the best one, verification of results for an alternative optimization algorithm (Section 5). The authors' intention is to convince the reader that the degree of complexity of the control system is justified due to the effects of regulation and meeting the expectations of smart grids.

Review of Indicators for Assessing the Quality of Voltage Control in MV Networks
When considering advanced methods of voltage control or planning in the MV network, the basic step is to determine its quality indicator. As mentioned earlier, when the consumers are supplied with a voltage other than the nominal one, their efficiency and reliability are deteriorated. Therefore, one of the most important objectives of voltage control in the distribution network is to keep the voltage as close as possible to the nominal value in all its nodes. Works [11,12] used the indicator (wsk U ) described by the formula: where U i is the voltage at the i-th node, U n is the nominal voltage of the network, and N is the number of network nodes.
In the works mentioned above, the value of indicator (1) was optimized by adjusting only the transformer tap changer. Similar considerations were conducted in works [1,13,14].
Studies [15,16] considered an index of a modified form: where U oi is the expected voltage value in the i-th network node. In this case, the control process covered both the transformer transmission and the reactive power of the sources, but only in a few operating states. A linear grid model was developed, which reflected the operation of the real grid. The applied linear optimization method was locally convergent. Similar considerations were conducted, among others, in [2,[17][18][19], where other methods of voltage control were additionally considered. In article [20], three objective functions were analyzed, taking into account the deviations in voltages and reactive power in generator nodes from the set values, with appropriate weighting factors. Similar considerations were conducted in paper [21]. An interesting function was considered in [22]: where VD abc is the positive three-phase voltage deviation function. The authors of the work [22] minimized the voltage deviation in the network taking into account individual phases, thus taking into account the asymmetry of phase voltages.
The second problem of voltage control, which was also mentioned earlier, is the economic aspects of distribution network operation, on which part of the research work devoted to these issues focuses. In paper [15], the following 'cost' objective functions can be found: • minimum cost of economic losses for consumers where K G is the cost of losses incurred by costumers due to voltage deviation from the nominal value; • the minimum cost of power and energy losses in the network to the network operator where K ∆P is the cost of power losses and K ∆A is the cost of energy losses in the grid; • minimum relative energy losses in the grid where ∆A is the energy losses in the network during optimization period T, A is the energy fed into the network during period T; • the maximum profit achieved by the network operator from the sale of electricity where D S are the revenues generated from the sale of electricity and K Z are costs of purchasing power and energy; • minimum total costs (both on the part of the operator and the costumers) where K S are the costs of power and energy losses in the grid, K G are costs of the economic losses incurred by costumers, • minimum costs of the distribution network operator where K U is the cost of the discounts granted to consumers in connection with exceeding the permissible voltage deviations.
Minimization of losses of active power in the grid (by changing transformer ratio and reactive power in sources) is considered in article [23]. In paper [24], the authors apply the objective function also taking into account the costs related to the flow of reactive power: (10) where C p is the cost of the active power, ∆P loss is active power losses, and Q cost are the costs related to reactive power flow.
Tengku H. et al. [25] analyzes a two-criteria objective function of the form: where U dev is the voltage deviation, U i is the actual voltage at bus i, and α, β are weighting factors. The voltage deviation in this case refers to the current voltage values at individual nodes. When searching for the voltage control quality indicator to be used in the real-time process (Figure 1), it was decided to use relations (1) and (11) in a compromise manner.

The Formulation of the Objective Function and Task Constraints
As mentioned earlier, the traditional setpoint voltage of the on-load tap changer control is generally set to a value slightly higher than the nominal mains voltage due to the expected voltage drops. In Poland, two settings are often used: Day and night, taking into account the load differences during the day. The values of set voltages for 15 kV networks range from 15.5 to 16.3 kV. Taking into account the above, in the conducted analyses, the setting of the voltage controller was assumed to be 15.75 kV (1.05 U nMV ).
The objectives of MV network voltage control discussed previously, qualitative and economic, can be mathematically mapped in different ways, as described in point 2. The problem is treated as an objective of nonlinear optimization with limitations; control variables can be continuous or discrete. After analyzing a number of quality indicators described in the literature, the paper considers two single-criterion objective functions and one two-criteria function. They were expressed in the following relationships: 1.
single-criterion objective function-voltage quality indicator: where U o is the expected voltage value in individual network nodes (the analyses assumed that U o = 1.05 · U n ), 2.
single-criterion objective function-grid active power losses (∆P loss ): 3. two-criteria objective function-sum of the F 1 (x, y, z) function and modified F 2 (x, y, z) function with weighting factors w 1 , w 2 : whereby: .. Q Gp -vector of control variables formed by the transformer ratio (ϑ-discrete variable) and reactive power of p sources connected to the MV network; y = U HV , P L1 ... P Lm , Q L1 ... Q Lm , P G1 ... P Gp -vector of independent variables, formed by: HV network supply voltage, active and reactive power received in m nodes, and power generated in p sources, not subject to change during optimization calculations; z = U 1 ... U j , δ 1 ... δ j -the vector of state variables containing nodal voltages and their arguments (total number of network nodes j = p + m).
The determination of the state vector z, for the known values of x and y vectors, is one of the basic computational problems included in the computer analysis of power systems. This problem, known as load flow analysis (LF), is described in textbooks such as [26,27].
Many algorithms and computer programs have been developed to solve this problem. Some of them are available as the Toolbox of the Matlab package (Matpower), but for the solutions provided in the proposed control system, they must be in the form of the dedicated 'calculation engine'. Generally speaking, the components of the unknown vector z (state variables) are determined from a system of nonlinear equations of the form: Each equation is formulated for a node of the analyzed network so that the power balance for this node is equal to zero.
On the basis of nodal voltages (state variables), in view of the knowledge of the impedance parameters of the analyzed network, the power losses are determined based on the formula: where m, n are terminal nodes of the network branch with impedance Z mn = R mn + jX mn ; all network branches form the set {G}. Considering the above-described objective functions together with the vectors x, y, z, the optimization task can be formulated in the following way: • objective function-the vector of control variables x providing the minimum function described (optionally) as (12), (13) or (14) is sought: F(x, y, z) → min (17) • vector of inequality constraints h, ensuring that the elements of the vector of state variables and elements of the vector of control variables are kept within the range specified by the technical requirements: Adopted as inequality restrictions were: • minimum and maximum values for transformer ratio (ϑ) • minimum and maximum reactive power values of each renewable source (Q Gk ) • minimum and maximum voltage values for all network nodes (U i ) • current-capacity limit values for power lines (I lmax ) • the permissible load of the HV/MV transformer (S nT )

Description of the Heuristic Optimization Algorithm
After defining the objective function, one should choose the method of solving the optimization task. Over the past half-century, the optimal power flow has gained great attention due to its importance in the power system operation of HV grids [28][29][30]. Optimal power flow is considered an important tool for efficient planning and enhancing the operation of electric power systems. The main task of optimal power flow is to determine the best or the most secure operating point (control variables) of certain objective functions while satisfying the system equality and inequality constraints.
Given the small number of nodes in the distribution network and its simple (radial) structure, it is possible to apply real-time Optimal Power Flow (OPF) calculations into the voltage control concept, which is seen in Figure 1. In order to do so, it is necessary to create a network calculation model for each optimization cycle. This is possible thanks to the process of estimating its condition, which can be carried out thanks to measurements of voltage and power values- Figure 1. Further in this case, the radial structure and small size of the network allow this task to be performed in a much simpler mode than for multi-node meshed networks [31][32][33][34][35][36]. Being aware of the fundamental importance of the problem of network condition estimation for the effectiveness of the proposed voltage control solution, the authors decided to discuss this problem further in subsequent publications.
Several classical (deterministic) and recent (nondeterministic) heuristic optimization techniques have been proposed to find the solution to the optimal power flow problem. It should be pointed out that the classical methods may be trapped in a local optimum due to the nonlinearity of optimal power flow problems [37][38][39]; therefore, the metaheuristic optimization techniques are widely employed for solving the optimal power flow problems [38,40,41]. In part of the presented works, various heuristic methods were used to solve the optimization tasks [22,25,[42][43][44]. Compared to classical methods, heuristic methods are characterized by the fact that they do not require knowledge of the derivative form of the objective function; they are resistant to discontinuities of this function and to 'getting stuck' in the calculation process at a local minimum. In the case of a divergent iterative process, data can be reloaded. The best solution recently found is not lost, as the solution vector is remembered at each stage of the calculation. The classical method would then have to be interrupted and the calculation restarted. The optimization method used in the research is called the Algorithm of the Innovative Gunner (AIG) [45]. It differs from more widely known heuristic methods [46][47][48], the way of generating new (successive) solutions. Generally speaking, in most heuristic methods (nearly 100 such methods are defined and described [37,[49][50][51][52][53]), the iterative process of determining optimal values of the vector x components is similar, because most often, in the next iterative step, which can be considered as additive correction, the following action is performed: (24) and the correction of the l-th component of the vector, in the i-th iterative step, is determined as a function value: ∆x whereby function f l and vector ξ are characteristic of the heuristic method used [40,46,48].
In the case of the AIG algorithm, relation (24) is replaced by the expression: Thus, the originality of the AIG method lies in the fact that in subsequent steps, the components of the x vector are changed by multiplying, not by adding. As functions g l1 (ξ 1 ), g l2 (ξ 2 ), assuming that ξ 1 = α, ξ 2 = β, the following ones are respectively adopted: where α and β are 'correction angles', drawn from a range (−α max , α max ) and (−β max , β max ) using a normal or uniform probability distribution. The process of solving optimization tasks for the objective function and the vector of decision variables x is as follows [45]: a. drawing of initial values of the components of the decision vector x 0 , b.
finding the objective function F value for the start vector x 0 , c.
setting the iteration counter to t = 1, d.
determining the scope of α max and β max , e.
correction of the decision vector x components by means of multiplicative correction factors g(α), g(β), f. determination of a new state of the analyzed system (in the case of the OPF task, solution of the system of Equation (15)), g.
determining the value of objective function F for subsequent solutions, checking if the objective function for the next solution is better than the solution from the previous iteration (if so, replace this solution with the current-better-solution, and if not-leave it), identification of the q solution for which the function reaches the minimum value, h.
checking whether for the identified solution q, the relation F x i+1 q < F best x i best occurs; if so, then accepting it as the next decision vector in the iterative process, i.e., verification of the stopping criterion (t + 1 = t max ), end or subsequent iteration (t = t + 1), and return to point d, Considering a population of new solutions in the AIG algorithm, it is characterized by a large exploration of the search space. The subsequent solutions differ (sometimes quite significantly) from the previous ones, which allows the avoiding of the local optimum. The effectiveness of the AIG algorithm has been confirmed for many test tasks in the area of mathematics and technology [45]. In power engineering applications (also OPF, [41]), the AIG algorithm works on the basis of an autonomous 'calculation engine', performing load flow (LF) calculations in subsequent steps.
The authors wish to emphasize that the use of the AIG algorithm results from the fact it was invented by them, and new confirmations are sought that, despite its simplicity, it can be effectively applied in many optimization tasks. In order to confirm the correctness of the obtained results, calculations were also carried out with the use of the CS (Cuckoo Search [48]) algorithm, one of the best-known metaheuristic algorithms.

Description of the Test Network and Methodology of the Tests Carried Out
In order to verify the effectiveness of the proposed voltage regulation concept, a multi-variant test calculation was performed. For the analyses, a network model consisting of a 15 kV overhead power line with branches, with a total length of 33.5 km, made with ACSR 70 mm 2 conductors, was used. It corresponds to the actual line located in the central part of Poland in Figure 2. There is a transformer with a power of 10 MVA and a ratio ϑ = 115/16.5 kV/kV ± 10% (taps every 1.11%). Four sources with rated capacities of P nG1 = 1 MW, P nG2 = 2 MW, P nG3 = 0.5 MW, and P nG4 = 0.5 MW are connected to the network. Each source has the capacity to generate reactive power in the range (−0.4, 0.4)P nG , and this power is remotely controlled in the system proposed in Figure 1.     A single solution to the optimization task would not allow us to draw reliable conclusions about the effectiveness of the proposed voltage control system using the HV/MV transformer tap changer and reactive power of generators connected to the MV network. A single solution to the optimization task would not allow us to draw reliable conclusions about the effectiveness of the proposed voltage control system using the HV/MV transformer tap changer and reactive power of generators connected to the MV network. That is why the search for the solution was repeated many times-in the initial stage, the number was one hundred. The values of the control variables were drawn by means of a random number generator with a rectangular distribution, from the following ranges: U HV voltage from the range (0.9, 1.1)U nHV • load at the node i 1 3 , 1 P Lmaxi • generation in node k including 1 3 , 1 P nGk Figure 4a shows the selected values of U HV voltage for a cycle of 100 draws, Figure 4b shows the drawn values of the power generated in sources, and Figure 4c shows the total value of the power received in the receiving nodes resulting from draws. In the case of the receiving nodes, the random load of reactive power determined by the coefficient tgφ = 0.4 was also modeled. That is why the search for the solution was repeated many times-in the initial stage, the number was one hundred. The values of the control variables were drawn by means of a random number generator with a rectangular distribution, from the following ranges: UHV voltage from the range ( ) nHV Figure 4a shows the selected values of UHV voltage for a cycle of 100 draws, Figure  4b shows the drawn values of the power generated in sources, and Figure 4c shows the total value of the power received in the receiving nodes resulting from draws. In the case of the receiving nodes, the random load of reactive power determined by the coefficient tgϕ = 0.4 was also modeled.
There is no correlation between the drawn values of UHV voltage, loads in nodes, and active power generated in sources. These values do not represent a specific time of day or weather condition. The intention of the authors was to check the effectiveness of optimization calculations for a large number of completely random independent variable values.
The results of these analyses should be regarded as positive.   There is no correlation between the drawn values of U HV voltage, loads in nodes, and active power generated in sources. These values do not represent a specific time of day or weather condition. The intention of the authors was to check the effectiveness of optimization calculations for a large number of completely random independent variable values. The results of these analyses should be regarded as positive.

The Way of Presenting the Research Results
In order to justify the effectiveness of the proposed method of voltage control in the MV network, the results of the optimization process repeated for 100 calculation cycles were presented. For each of them, a set of independent variables (y vector) was randomly generated. In order to justify the final choice of the objective function, the effects of the six considered methods of control have been compared graphically below: a.
the effects of traditional control, i.e., maintaining a constant voltage value of 1.05 U nMV , on the MV busbars, using OLTC; b.
control effects, the aim of which is to minimize the single-criterion objective function (12), i.e., to minimize the square of the relative voltage deviations from the value of 1.05 U nMV determined for all network nodes, using OLTC and appropriately selected reactive power values of distributed sources, as for all the next cases; c. effects of regulation, the aim of which is to minimize the single-criteria objective function (13), i.e., to minimize power losses in the MV grid; d.
the effects of minimizing the two-criteria objective function (15), where the voltage deviation component has the weighting factor w1 = 0.1, and the relative power loss indicator w2 = 0.9; e. the effects of minimizing the two-criteria objective function (15), where the voltage deviation component has the weighting factor w1 = 0.9, and the relative power loss indicator w2 = 0.1; f. the effects of minimizing the two-criteria objective function (15), where the voltage deviation component has the weighting factor w1 = 0.7, and the relative power loss indicator w2 = 0.3. Figure 5 shows the voltage values for all network nodes, obtained for all 100 calculation cycles. The charts form characteristic 'bands' that are used for the qualitative assessment of the effects of control. The larger the band area, the worse the concentration of the voltages around 1.05 U nMV is, and therefore, the control is less effective. Undoubtedly, the greatest approximation of the voltage value to the expected value 1.05 U nMV is shown in Figure 5b,e for which the voltage criterion has the only or dominant significance. Minimizing the power losses (Figure 5c,d) noticeably worsens the voltage quality. For the two-criteria optimization and the weighting factors w1 = 0.7 and w2 = 0.3, the band area slightly widens, compared to the single-criteria optimization according to (12). ment of the effects of control. The larger the band area, the worse the concentration of the voltages around 1.05 UnMV is, and therefore, the control is less effective. Undoubtedly, the greatest approximation of the voltage value to the expected value 1.05 UnMV is shown in Figures 5b,e for which the voltage criterion has the only or dominant significance. Minimizing the power losses (Figure 5c,d) noticeably worsens the voltage quality. For the twocriteria optimization and the weighting factors w1 = 0.7 and w2 = 0.3, the band area slightly widens, compared to the single-criteria optimization according to (12).

Quantitative Assessment of Voltage Control
The qualitative evaluation of the voltage for various control methods was supplemented with the quantitative evaluation with the use of the indicator value (12), as shown in Figure 6. The results for the OLTC regulation case, shown in Figure 6a,b,d (blue), were taken as reference. The average value of the indicator (for all calculation cycles) was F1av = 0.0262, while the worst value was F1 = 0.1052. Optimization according to function (12)- Figure 6a-reduces the value of the voltage quality indicator so dramatically that while keeping the y-axis scale the same as for the traditional regulation, its values obtained for the subsequent calculation cycles are almost invisible. The worst case is F1 = 0.0037, mean value F1av = 0.0010-about 26 times less than in the case of the traditional solution. Minimizing power losses causes a very slight reduction in the voltage quality indicator ( Figure  6b). The same indicator is very sensitive to the optimization process, because even with the weighting factor w1 = 0.1, its value significantly improves compared to the reference value (Figure 6c). For the weighting factor w1 = 0.3, the two-criteria optimization is almost as effective as for the single-criteria approach (12) (Figure 6d).

Quantitative Assessment of Voltage Control
The qualitative evaluation of the voltage for various control methods was supplemented with the quantitative evaluation with the use of the indicator value (12), as shown in Figure 6. The results for the OLTC regulation case, shown in Figure 6a,b,d (blue), were taken as reference. The average value of the indicator (for all calculation cycles) was F 1av = 0.0262, while the worst value was F 1 = 0.1052. Optimization according to function (12)- Figure 6a-reduces the value of the voltage quality indicator so dramatically that while keeping the y-axis scale the same as for the traditional regulation, its values obtained for the subsequent calculation cycles are almost invisible. The worst case is F 1 = 0.0037, mean value F 1av = 0.0010-about 26 times less than in the case of the traditional solution. Minimizing power losses causes a very slight reduction in the voltage quality indicator (Figure 6b). The same indicator is very sensitive to the optimization process, because even with the weighting factor w1 = 0.1, its value significantly improves compared to the reference value (Figure 6c). For the weighting factor w1 = 0.3, the two-criteria optimization is almost as effective as for the single-criteria approach (12) (Figure 6d).

Power Losses, Selection of Weighting Factors
Good voltage quality, although justified from the consumers' point of view, does not significantly reduce the value of power losses. This is due to the fact that reactive power is used in the control process, which increases the value of losses. This is shown in Figure 7. For the traditional control, the mean losses ∆P loss = 34.1 kW, and after optimization (according to voltage quality index), it changed very slightly to ∆P loss = 28.6 kW- Figure 7a. On the other hand, the value of the objective function (level of power losses) has clearly decreased- Figure 7b. The average value of F 2 (power losses (13)) after optimization decreased significantly, to ∆P loss = 19.3 kW.
The qualitative evaluation of the voltage for various control methods was supplemented with the quantitative evaluation with the use of the indicator value (12), as shown in Figure 6. The results for the OLTC regulation case, shown in Figure 6a,b,d (blue), were taken as reference. The average value of the indicator (for all calculation cycles) was F1av = 0.0262, while the worst value was F1 = 0.1052. Optimization according to function (12)- Figure 6a-reduces the value of the voltage quality indicator so dramatically that while keeping the y-axis scale the same as for the traditional regulation, its values obtained for the subsequent calculation cycles are almost invisible. The worst case is F1 = 0.0037, mean value F1av = 0.0010-about 26 times less than in the case of the traditional solution. Minimizing power losses causes a very slight reduction in the voltage quality indicator ( Figure  6b). The same indicator is very sensitive to the optimization process, because even with the weighting factor w1 = 0.1, its value significantly improves compared to the reference value (Figure 6c). For the weighting factor w1 = 0.3, the two-criteria optimization is almost as effective as for the single-criteria approach (12) (Figure 6d).

Power Losses, Selection of Weighting Factors
Good voltage quality, although justified from the consumers' point of view, does not significantly reduce the value of power losses. This is due to the fact that reactive power is used in the control process, which increases the value of losses. This is shown in Figure  7. For the traditional control, the mean losses ∆Ploss = 34.1 kW, and after optimization (according to voltage quality index), it changed very slightly to ∆Ploss = 28.6 kW- Figure 7a. On the other hand, the value of the objective function (level of power losses) has clearly decreased- Figure 7b. The average value of F2 (power losses (13)) after optimization decreased significantly, to ∆Ploss = 19.3 kW.
The analyses carried out for the one-criteria objective functions F1 and F2 indicate that effective optimization of the voltage indicator does not decrease power losses, and minimizing losses has a negative impact on the value of the voltage indicator. Therefore, an attempt to apply polyoptimization is justified, according to the proposed function (15). The question of selecting the 1 2 ( , ) w w weighting factors remains. Below, Table 1 presents the results of two-criteria optimization for five pairs of these factors. As expected, the use of a pair of weighting factors (w1 = 0.1, w2 = 0.9) gives better results in reducing power losses than (w1 = 0.9, w2 = 0.1)- Figure 7c. Both Table 1 and Figures 5f, 6d, and 7d confirm that the selection of the weighting factors (w1 = 0.7, w2 =

The Values of Decision Variables
The following figures show the values of the decision variables (vector x) occurring in successive variants of the optimal voltage control process. Figure 8 shows the tap The analyses carried out for the one-criteria objective functions F 1 and F 2 indicate that effective optimization of the voltage indicator does not decrease power losses, and minimizing losses has a negative impact on the value of the voltage indicator. Therefore, an attempt to apply polyoptimization is justified, according to the proposed function (15). The question of selecting the w1, w2 weighting factors remains. Below, Table 1 presents the results of two-criteria optimization for five pairs of these factors. Table 1. Summary of results of two-criteria optimization (15) for different weighting factors (w1-responsible for voltage quality, w2-responsible for power losses), one hundred calculation cycles.

Weight w1
Weight As expected, the use of a pair of weighting factors (w1 = 0.1, w2 = 0.9) gives better results in reducing power losses than (w1 = 0.9, w2 = 0.1)- Figure 7c. Both Table 1 and Figures 5f, 6d and 7d confirm that the selection of the weighting factors (w1 = 0.7, w2 = 0.3) gives both a significant reduction in the voltage quality indicator and a reduction in power losses by an average of 40%, hence the authors' recommendation for these factor values.

The Values of Decision Variables
The following figures show the values of the decision variables (vector x) occurring in successive variants of the optimal voltage control process. Figure 8 shows the tap changer numbers of the OLTC determined for voltage quality optimization (a), power loss minimization (b), and two-criteria optimization for the three pairs of weighting factors (c, d, f) listed in Table 1. Figure 9 shows the values of the reactive power of the G1, G2, G3, G4 generators corresponding to the above-mentioned optimization processes. It is difficult to formulate general conclusions based on the qualitative observation of changes in the presented values, but it can be noticed that the optimization process adjusts the voltage values to the form of the objective function. Thus, it is not justified to make simplified statements that distributed sources overestimate the voltage values in the network; these voltages vary within wide limits, and only the comprehensive interaction of control variables effectively reduces this variability and equalizes the voltage in the entire network.

The Effectiveness of the AIG Optimization Algorithm and Verification of the Calculation Results
In order to verify the correctness of the proprietary AIG algorithm, analogous calculations were performed according to the very popular metaheuristic algorithm-Cuckoo Search (CS). The results are shown in Figure 10 and in Tables 2-4. In all cases, the optimization process is very fast; about 50 iterations lead practically to the minimum value of the objective function.

The Values of Decision Variables
The following figures show the values of the decision variables (vector x) occurring in successive variants of the optimal voltage control process. Figure 8 shows the tap changer numbers of the OLTC determined for voltage quality optimization (a), power loss minimization (b), and two-criteria optimization for the three pairs of weighting factors (c, d, f) listed in Table 1. Figure 9 shows the values of the reactive power of the G1, G2, G3, G4 generators corresponding to the above-mentioned optimization processes. It is difficult to formulate general conclusions based on the qualitative observation of changes in the presented values, but it can be noticed that the optimization process adjusts the voltage values to the form of the objective function. Thus, it is not justified to make simplified statements that distributed sources overestimate the voltage values in the network; these voltages vary within wide limits, and only the comprehensive interaction of control variables effectively reduces this variability and equalizes the voltage in the entire network.

The Effectiveness of the AIG Optimization Algorithm and Verification of the Calculation Results
In order to verify the correctness of the proprietary AIG algorithm, analogous calculations were performed according to the very popular metaheuristic algorithm-Cuckoo Search (CS). The results are shown in Figure 10 and in Tables 2-4. In all cases, the optimization process is very fast; about 50 iterations lead practically to the minimum value of the objective function.  (12), (b) optimization according to (13), (c) polyoptimization according to (15) for W1 = 0,7, W2 = 0,3.
Below, Table 2 presents the optimal voltage values in network nodes for two heuristic algorithms (AIG and CS) and three objective functions, for one of the calculation cycles. The table also includes the voltage values before optimization and the load power. The random voltage value on the HV side was UHV = 1.0514 p.u.      Table 2 presents the optimal voltage values in network nodes for two heuristic algorithms (AIG and CS) and three objective functions, for one of the calculation cycles. The table also includes the voltage values before optimization and the load power. The random voltage value on the HV side was U HV = 1.0514 p.u.
Below, Table 3 presents the optimal values of the reactive power of sources for two heuristic algorithms (AIG and CS) and three objective functions for one sample of the calculation cycles (the same as in Table 2). Table 4 contains the optimal values of the HV/MV transformer tap changer and optimal value of the objective function, for two heuristic algorithms (AIG and CS) and three objective functions for one sample of the calculation cycles (the same as in Tables 2 and 3).
Full compliance of the results obtained with the proprietary AIG algorithm and the well-known CS algorithm, with a noticeably better convergence and speed of AIG, confirms its good properties demonstrated in other applications.

Discussion and Conclusions
The voltage control system proposed in Figure 1 for MV networks with distributed sources allows the achievement of favorable voltage level, both in terms of values and in terms of minimizing power losses. This is possible through the traditional application of the on-load tap changer of the HV/MV transformer and the innovative, appropriate use of the reactive power control capabilities of distributed sources. The key to achieving a control quality significantly better than that of traditional fixed-value control is the use of an algorithm to minimize the objective function, with the best results being achieved through a two-component function-based polyoptimization (14).
The concept presented contains a number of assumptions simplifying and idealizing the optimization process, performed as a heuristic optimization using the proprietary AIG algorithm. This process can take place provided that the network model is properly estimated and iterative load flow calculations, which are an essential element of the OPF algorithm, are performed on it. This issue was not discussed in detail in the article, but the technological development of network data transmission and processing systems makes it possible to look realistically at obtaining such possibilities on-line at the control system level. The concept is realistic because the dimension of the calculation task is small, and there are no particularly high requirements in view of the speed of the control process. In this way, the characteristics of smart grids will be achieved. It will be possible to introduce such operation at the MV level as the answer of control variables to the load changes and generation variability of distributed sources.