# Research on Reconfiguration of Distribution Network considering Three-Phase Unbalance

Authors: Xuejie Wang; Yanchao Ji; Jianze Wang; Yi Zhao; Peng Ye; Lei Qi; Shuo Yang; Siqi Liu
Journal: Wireless Communications and Mobile Computing
 (2022)

---

## Abstract

In the distribution network, the addition of distributed power sources can improve the voltage part of the distribution system, provide uninterrupted power supply, and reduce network losses. The location of the power supply and its output will also directly affect the system voltage and power loss. Therefore, determining the optimal installation location and capacity of distributed power generation is of great significance to the distribution network. However, in the actual system, due to the unbalanced load distribution and the three-phase unbalance of the transmission line, the power distribution system presents the characteristics of unbalanced distribution. As the number of electric vehicles increasing, the imbalance of the distribution network has become more prominent.


---

## Body

## 1. Introduction

Distribution network reconfiguration is of great significance to the planning and operation of the distribution network. Reconfiguration can determine the topology of the network to improve the economy and stability of the operation of the system. Distributed power has aggravated the three-phase imbalance of system. The goal of the reconfiguration program is to determine the optimal network topology path by controlling the switch state, to meet the system’s minimum network loss and load balance objective function, and to some given distribution network operation constraints [1, 2]. Due to the nonlinearity of power flow and the discreteness of switching states, the problem of distribution network reconstruction can be regarded as a mixed integer nonlinear problem in mathematics. In the previous theoretical research, heuristic algorithms have been studied in depth [3, 4]. In addition, many global optimization algorithms are applied, such as genetic algorithm [5], artificial neural networks [6], and particle swarm optimization [7]. The calculations of these methods are usually more complicated, and they cannot guarantee the global optimal solution. In mathematical optimization problems, reconstruction models are usually expressed as linear programming problems [8], quadratic programming problems [9], etc. Although a lot of results have been obtained in the algorithm for solving the reconstruction problem, most of the existing algorithms are obtained under the assumption that the distribution network is a three-phase balanced system. In the actual system, due to the unbalanced load and network structure, the distribution network operates in a three-phase unbalanced state. In order to solve the above imbalance problem, the literature [10] proposed a real-time distribution network reconstruction algorithm for three-phase unbalanced systems.


## 2. Three-Phase Unbalanced System Model

In a three-phase unbalanced power distribution system, the relationship between voltage and current can be expressed in Figure1.
(1)ViaVibVic=VjaVjbVjc+zaazbazcazabzbbzcbzaczbczcc⋅IiaIibIic.Figure 1
Improved hierarchical decentralized reconstruction method process.The load model can be regarded as a negative sequence injection current in the system, and each load is assumed to be a constant power component, a linear combination of a constant impedance component, and a constant current component. Therefore, the three-phase injection current at thei bus can be expressed as
(2)IL,iabc=α1,i⋅IP,iabc+α2,i⋅IZ,iabc+α3,i⋅II,iabc,where IP,iabc is the constant power three-phase injection current component; IZ,iabc is the constant impedance three-phase injection current component; and II,iabc is the three-phase injection current component of constant current.The distributed power model at thei bus can be regarded as the positive sequence injection current as
(3)IDG,iabc=PDG,ia‐jQDG,iaVia∗PDG,ib‐jQDG,ibVib∗PDG,ic‐jQDG,icVic∗T,where PDG,ia, PDG,ib, and PDG,ic are, respectively, the three-phase active power of the distributed power supply at bus iand QDG,ia, QDG,ib, and QDG,ic are, respectively, the three-phase reactive power of the distributed power supply at bus i.When the distributed power supply operates in constant power operation mode, the equivalent injection current can be directly derived from Equation (3). When the distributed power supply operates in constant voltage mode, a double-loop calculation is required to obtain the equivalent injection current. The inner-loop calculation yields the amount of distributed power supply reactive power required to keep the bus voltage in a fixed range. The outer-loop calculation can obtain the injection current from the initial active power and the calculated reactive power. In a three-phase asymmetric system with two-phase or single-phase branches, the impedance of the missing phase in Equation (1) is zero, and the voltage and current of the missing phase can be excluded from the calculation results. Correspondingly, the current values of the missing phases in Equations (2) and (3) can be set to zero.Distribution network reconfiguration is generally viewed as feeder reconfiguration optimization problem, and it is usually assumed that only three-phase feeder branches are available for reconfiguration. Any line is equipped with three-phase switching equipment, and the switch states can be defined as(4)Sj=1,Switchjisclosed,consistentwiththeinitialdirection0,Switchjisoff−1,Switchjisclosed,oppositetotheinitialdirection,where the direction of reference is determined by the direction of current.The network connectivity of the distribution network is characterized by the node-branch correlation matrix. If the system is an ideal three-phase symmetric system, a single-phase equivalent circuit is usually adopted, as well as the node-branch correlation matrixAbalanced determined by the structure of system. The initial test state of the distribution network assumes that the switches are all closed. The initial node-branch correlation matrix A0 of the network is symmetric, and A0 is constant for a given network.whereaij0 is the matrix elements of A0.If the distribution network is a three-phase asymmetric system, the node-branch correlation matrix is three-phase and is obtained by multiplyingAbalanced by a unit matrix of order 3:
(5)Ai,j=aij000aij000aij=aij0000aij0000aij0⋅Sj.And(6)Ibus=A⋅Ybranch⋅AT⋅Vbus=Ybus⋅Vbus.Among them,WhereIbus is the node injection current vector; Vbus is the bus voltage vector; Ybus is the nodal conductivity matrix; and Ybranch is the tributary conductivity matrix.
(7)Ibus=IDG−IL,where IDG is the distributed power injection current vector and IL is the load injection current vector.According to Equations (6) and (7), the tidal current calculation can be expressed as follows: the initial voltage Vbusk−1 is given, the nodal injection current Ibusk can be solved by Equation (7), and Vbusk is calculated by substituting Ibusk into Equation (6). When the termination condition is satisfied, the iterative termination calculation in turn leads to the system nodal voltage and branch current values.Equation (7) can in turn be expressed as
(8)Ibus,il=∑k=1n∑p=actiklp⋅vkp,l=phasea,b,c.In the formula,tiklp=∑j=1maij0akj0yjlp⋅Sj2≜giklp+jbiklp.yjaa,yjab,yjac;yjba,yjbb,yjbc;yjca,yjcb,yjcc is the conjugate inverse matrix of the jth branch impedance matrix. Due to the existence of the relation of Equation (9):
(9)Ibus,il=Pinject,il−j⋅Qinject,ileil−j⋅fil,where Pinject,il,Qinject,il is the active and reactive injected power at node i of a phase.The tide calculation formula can be obtained as(10)Pinject,il=∑k=1n∑p=aceilgiklp⋅ekp−biklp⋅fkp+filgiklp⋅fkp−biklp⋅ekpQinject,il=∑k=1n∑p=acfilgiklp⋅ekp−biklp⋅fkp−eilgiklp⋅fkp+biklp⋅ekp.In power loss, in an unbalanced system, the result of power loss in the branch is the input power minus the output power of each phase, so we can obtain(11)AT⋅Vbus=Zbranch⋅Ibranch.Substituting Equations (1) and (5) into Equation (11), the total active power loss can be obtained as
(12)Ploss=Re∑j=1n∑k=acVjk∗⋅∑i=1n∑p=acVip⋅∑i=1mail0ajl0ylpk⋅sl2=∑j=1n∑k=ac∑i=1n∑p=acejk⋅eip⋅gijpk+fip⋅bijpk+fjk⋅fip⋅gijpk−eip⋅bijpk.The magnitude of the voltage depends on the network topology of the system and the output power of the distributed power supply, which are the two variables that satisfy the power loss minimization objective function. In addition to this, the operational constraints have to be satisfied:


### 2.1. Voltage Amplitude Constraint

(13)Vimin≤Via,Vib,Vic≤Vimax,i=1,2,⋯,N,where Vimin and Vimax are the allowable minimum and maximum voltages of node i, respectively, so that Vimin is 0.95 times the rated voltage and Vimax is 1.05 times the rated voltage.


### 2.2. Node Voltage Imbalance Constraint

Voltage imbalance can be obtained by dividing the maximum voltage deviation by the average voltage amplitude of the three phases. According to ANSI C84.1-2011, the voltage imbalance should be controlled within 3%. Therefore, the voltage imbalance constraint can be expressed as(14)Vip−avgiavgi≤3%,where,avgi=∑p=acVip/3,p=a,b,c.


### 2.3. Current Constraint

(15)Ibranch,ip≤Ii,max,where Ii,max is the load capacity of branch i.


### 2.4. Network Radiation-like Structural Constraints

Network constraints and balanced systems are systems where the distribution network remains radial in operation. The network of unbalanced system is radial network.(16)∑k=1MSk=N−d,where d is the number of all balancing nodes.All loads are supplied with uninterrupted power supply,rankA=N−d. At least one branch in each loop is open and closed.
(17)∑i=1MkSi≤Mk−1,where Mk is the number of branches of loop k.


## 2.1. Voltage Amplitude Constraint

(13)Vimin≤Via,Vib,Vic≤Vimax,i=1,2,⋯,N,where Vimin and Vimax are the allowable minimum and maximum voltages of node i, respectively, so that Vimin is 0.95 times the rated voltage and Vimax is 1.05 times the rated voltage.


## 2.2. Node Voltage Imbalance Constraint

Voltage imbalance can be obtained by dividing the maximum voltage deviation by the average voltage amplitude of the three phases. According to ANSI C84.1-2011, the voltage imbalance should be controlled within 3%. Therefore, the voltage imbalance constraint can be expressed as(14)Vip−avgiavgi≤3%,where,avgi=∑p=acVip/3,p=a,b,c.


## 2.3. Current Constraint

(15)Ibranch,ip≤Ii,max,where Ii,max is the load capacity of branch i.


## 2.4. Network Radiation-like Structural Constraints

Network constraints and balanced systems are systems where the distribution network remains radial in operation. The network of unbalanced system is radial network.(16)∑k=1MSk=N−d,where d is the number of all balancing nodes.All loads are supplied with uninterrupted power supply,rankA=N−d. At least one branch in each loop is open and closed.
(17)∑i=1MkSi≤Mk−1,where Mk is the number of branches of loop k.


## 3. Optimal Planning of DG Units

According to the literature [11], the results show that the factors affecting the system power loss include the network topology, the location, and the output power of DG units. The sensitivity of the injected active power loss per bus can be used to determine the optimal bus location for installing the distributed power units.


### 3.1. DG Optimal Location of the Unit

Since the integration of distributed power units increases the bus active power injection, the best installation location is the bus with the maximum negative active loss sensitivity [12] so that power loss will be effectively controlled. Define the power loss Equation (13) as
(18)Ploss=ge,f.The relationship between the voltage vector and the injected power is determined from Equation (10):
(19)he,f,Pinject,Qinject=0.If the small variationΔP,ΔQT is added to the injected power vector Pinject,QinjectT, the variation of the voltage vector can be solved by Equation (15):
(20)∂h∂e∂h∂fx0⋅△eΔf+∂h∂Pinject∂h∂Qinjectx0⋅△PinjectΔQinject=0.Therefore, the amount of change in power loss caused is shown in Equation (21):
(21)△Ploss=∂g∂e∂g∂f⋅△e△f.Then, the power loss sensitivity vector for each bus injected power is(22)Ms=‐∂g∂e∂g∂f⋅∂h∂e∂h∂f−1⋅∂h∂Pinject∂h∂Qinject.The sensitivity vector is solved to obtain the power loss sensitivity of the injected power of each phase of the busbar. It is considered that only the busbars of three phases can be used as candidate installation locations; then, the sensitivity of each busbar is taken as the average of the sensitivity of the three phases.


### 3.2. Optimal Capacity of DG Units

The optimal capacity of the DG unit minimizes the power loss in the unit in an unbalanced system that keeps the initial topology constant [13].This case can be considered an undesirable situation in reconfiguration studies, due to the fact that the DG unit needs to emit the maximum power and there is no other way to additionally support the reconfiguration. The influence of the optimal capacity of distributed generation is as follows: The optimal capacity of the DG unit is solved so that the power loss of the unit can be minimized in the unbalanced distribution system with initial topology. This scenario can be considered the worst case in reconfiguration research. Because DG units must generate maximum power into the grid, there is no need for additional support to reconfigure the network. After installing the DG unit and arranging its optimum capacity, the power loss of the system has been minimized under the initial system structure. Power optimization may be further reduced if the system is reconfigured. In addition, due to time-varying load, power losses are not always minimized in the case of a fixed network structure and constant DG output power. It is therefore necessary to reconfigure the network and reduce the DG power from time to time [14].The optimization problem of the optimal capacity can be expressed as(23)minuJ=Plossx,us.t.fx,u=0gx,u≤0.In the formula:x=e,fT, u=PDG,QDGT, fx,u is the tidal equation, and gx,u denotes the inequality constraint of Equations (13)–(17). The inequality constraint can be eliminated by introducing a penalty function in the objective function, and Equation (23) can be expressed as
(24)minJuc=Plossx,u+∑i=1Hϕiβi,gi,where H is the total number of inequality constraints.Each penalty function is defined as(25)ϕiβi,gi=0,ifgi≤0βi,gi2,ifgi>0andβi>0.The tide calculation is first performed to obtain the power loss such that the equation constraintfx,u is always satisfied. Thus, Equation (25) becomes an unconstrained optimization problem, and its minimum can be solved using the proposed Newtonian numerical method; the differentiation of Equation (26) is
(26)Fu=dJucx,udu=dPlossx,udu+ddu∑i=1Hϕiβi,gi.In the formula, the first part of the rightmost side of the equation is the same as the sensitivity matrix, whose matrix columns represent the bus of the selected DG unit. The second term can be obtained by differentiating the inequality constraint. Since the second-order partial derivative equation of the multivariate function is difficult to solve, the minimum of Equation (26) can be solved by the cutline method with positive definite cutline update.
(27)uk+1=uk−Hk−1⋅Fuk,sk=uk+1−uk,yk=Fuk+1−Fuk,Hk+1=Hk+yk⋅ykTykT⋅sk−Hk⋅sk⋅skT⋅HkskT⋅H⋅sk.


## 3.1. DG Optimal Location of the Unit

Since the integration of distributed power units increases the bus active power injection, the best installation location is the bus with the maximum negative active loss sensitivity [12] so that power loss will be effectively controlled. Define the power loss Equation (13) as
(18)Ploss=ge,f.The relationship between the voltage vector and the injected power is determined from Equation (10):
(19)he,f,Pinject,Qinject=0.If the small variationΔP,ΔQT is added to the injected power vector Pinject,QinjectT, the variation of the voltage vector can be solved by Equation (15):
(20)∂h∂e∂h∂fx0⋅△eΔf+∂h∂Pinject∂h∂Qinjectx0⋅△PinjectΔQinject=0.Therefore, the amount of change in power loss caused is shown in Equation (21):
(21)△Ploss=∂g∂e∂g∂f⋅△e△f.Then, the power loss sensitivity vector for each bus injected power is(22)Ms=‐∂g∂e∂g∂f⋅∂h∂e∂h∂f−1⋅∂h∂Pinject∂h∂Qinject.The sensitivity vector is solved to obtain the power loss sensitivity of the injected power of each phase of the busbar. It is considered that only the busbars of three phases can be used as candidate installation locations; then, the sensitivity of each busbar is taken as the average of the sensitivity of the three phases.


## 3.2. Optimal Capacity of DG Units

The optimal capacity of the DG unit minimizes the power loss in the unit in an unbalanced system that keeps the initial topology constant [13].This case can be considered an undesirable situation in reconfiguration studies, due to the fact that the DG unit needs to emit the maximum power and there is no other way to additionally support the reconfiguration. The influence of the optimal capacity of distributed generation is as follows: The optimal capacity of the DG unit is solved so that the power loss of the unit can be minimized in the unbalanced distribution system with initial topology. This scenario can be considered the worst case in reconfiguration research. Because DG units must generate maximum power into the grid, there is no need for additional support to reconfigure the network. After installing the DG unit and arranging its optimum capacity, the power loss of the system has been minimized under the initial system structure. Power optimization may be further reduced if the system is reconfigured. In addition, due to time-varying load, power losses are not always minimized in the case of a fixed network structure and constant DG output power. It is therefore necessary to reconfigure the network and reduce the DG power from time to time [14].The optimization problem of the optimal capacity can be expressed as(23)minuJ=Plossx,us.t.fx,u=0gx,u≤0.In the formula:x=e,fT, u=PDG,QDGT, fx,u is the tidal equation, and gx,u denotes the inequality constraint of Equations (13)–(17). The inequality constraint can be eliminated by introducing a penalty function in the objective function, and Equation (23) can be expressed as
(24)minJuc=Plossx,u+∑i=1Hϕiβi,gi,where H is the total number of inequality constraints.Each penalty function is defined as(25)ϕiβi,gi=0,ifgi≤0βi,gi2,ifgi>0andβi>0.The tide calculation is first performed to obtain the power loss such that the equation constraintfx,u is always satisfied. Thus, Equation (25) becomes an unconstrained optimization problem, and its minimum can be solved using the proposed Newtonian numerical method; the differentiation of Equation (26) is
(26)Fu=dJucx,udu=dPlossx,udu+ddu∑i=1Hϕiβi,gi.In the formula, the first part of the rightmost side of the equation is the same as the sensitivity matrix, whose matrix columns represent the bus of the selected DG unit. The second term can be obtained by differentiating the inequality constraint. Since the second-order partial derivative equation of the multivariate function is difficult to solve, the minimum of Equation (26) can be solved by the cutline method with positive definite cutline update.
(27)uk+1=uk−Hk−1⋅Fuk,sk=uk+1−uk,yk=Fuk+1−Fuk,Hk+1=Hk+yk⋅ykTykT⋅sk−Hk⋅sk⋅skT⋅HkskT⋅H⋅sk.


## 4. Network Reconfiguration of DG Units

After installing the DG units and adjusting their optimal capacity as described above, the system network losses are at their lowest level under the initial network topology. A network reconfiguration is required to bring the power loss down further. In addition, due to the fluctuating load, the power loss is not always minimal with a constant topology and a constant DG [15]. Therefore, in each operational interim period, the objective function is defined as minimizing the total cost of power loss and limiting the output power of the DG set. The problem is described as
(28)minJ=w1⋅PlossS,PQDGact+w2∑i=1k∑p=acPDGmax,ip−PDGact,ipt+w3∑i=1k∑p=acQDGmax,ip−QDGact,ipt⋅△Ts.t.13∼17,andPDGact,ipt≤PDGmax,ip,QDGact,ipt≤QDGmax,ip,i=1∼k,p=a,b,c,t=1∼24.where T is the number of contact switches; K is the number of DG units; ΔT is the scheduled running time in hours; PDGmax,ip,QDGmax,ip is the optimal capacity of the i-th DG unit; and PDGact,ipt,QDGact,ipt is the p (p=a,b,c) phase active and reactive power output by the ith DG unit at time tof operation.Since the optimization problem of distribution network reconfiguration considering three-phase imbalance is a second-order partial derivative equation solution problem for multivariate functions, therefore, an improved hierarchical decentralized reconstruction method is adopted to simultaneously reconstruct the optimal topology of the distribution network considering the three-phase imbalance and regulate the DG output power. The flowchart is shown in Figure2:Figure 2
Improved hierarchical decentralized reconstruction method process.


### 4.1. Network Decomposition

Loops are more easily identified by the difference between closed contact switches and segmented opens for differentiation. In order to easily distinguish between tightly connected partitions and loosely connected partitions [16], the connectivity between the subdivisions is defined as follows:
(29)DA,B=NumberofpublicbusminTotalnumberofbusbarsinzoneA,TotalnumberofbusbarsinzoneB,IfDA,B=0, it means that partition A and partition B are independent of each other. Given a threshold δ, if, means that partition A and partition B are tightly connected and if DA,B≤δ means that partition A and partition B are loosely connected. However, if the threshold δ is set too small, the loosely connected partitions are not easily identified; if the threshold δ is set too large, the loosely connected partitions are easily expanded. Therefore, it is generally reasonable to choose a threshold δ of 0.2.The specific steps of network decomposition are as follows.The set of cut vertices of a connected graphG=V,E is a vertex set U⊆V satisfying the following conditions: G-U is not connected; G-K is connected, K⊂U; every vertex u of U in a connected graph G is connected to at least one vertex in every partition of G-U. A set of cut vertices is a cut vertex if there is only one vertex in the set.
(a)
The value of the modified adjacency matrixC is acquired. Let m be the number of loops, and the modified adjacency matrix is an m-dimensional square matrix(b)
If there are no partitions with a common bus other than the source node, these partitions are defined as basic partitions. Two basic partitions are connected to each other by a contact switch. The following steps are used for all basic partitions(c)
Draw the connectivity graph. In the connectivity graphG=V,E, the set of vertices V=v1,⋯,vnrepresents 1−n loops. The edge E connecting the vertices indicates that the loops are coupled to each other. The vertices are drawn in the studied region to represent the loops, and if the element C (i, j) is not zero, an edge is added between the vertices vi and vj. If loop i is loosely connected to loop j, then it is noted as L; otherwise, it is noted as T(d)
Be sure to check graph G for connectivity. If not, the isolated vertices need to be found, denoted asS. The corresponding loop of the isolated vertex contains the first member of the decomposition system(e)
Determine the cut vertex set of G-s. The search starts from the parent node, and the cut vertex set node is found from the child nodes associated with the parent node. If the cut vertex set is empty, the search moves to the vertex associated with the child node. If the cut vertex set is found, the cut vertex set is removed, and the associated edges define two separate components. Continue searching for each component’s cut vertex set until there are no more cut vertex sets. Each component represents the decomposed system, and the contact switch corresponding to the cut vertex set represents the decomposed system by the links between them(f)
In case two vertices are connected to the edge marked withL, their corresponding rings are decomposed into two subsystems. If two loosely connected vertices are tightly connected to another vertex, the contact switch corresponding to that vertex is an interconnection between the two subsystems. In the lowest level subsystem forming the basis of the entire system, which includes all busbars, loads, and DG units. Each higher-level subsystem consists of several subsystems representing the entire system


### 4.2. Applying a Multi-Agent Framework

Intelligent agents consisting of data units, computation units, and decision units need to be assigned to each subsystem for solving subproblems of the assigned subsystem and exchanging information with other agents. Since DG units may exist in the bottom most subsystem, the bottom most agent should be able to decide the optimal topology and the best DG output of the local system. The optimization problem for each bottom most agent layer is represented by Equation (29), and all variables are refined into local variables for its subsystem. Based on the decision plan solved by the lower layer agents, all upper layer agents ought to judge whether the common contact switch is off or not; hence, each embedded optimization problem is a pure reconstruction problem with the equation:
(30)minJ=PlossSsubsystem,Ssolved,PQsolved⋅ΔT.In the formula,Ssubsystem is the set of switching states to be solved; Ssolved is the set of switching states already solved by the lower layer agent; and PQsolved is the actual output power of the DG unit of the lowest layer agent.


### 4.3. Optimization Problems for the Lowest-Level Agent

The optimization problem defined in the bottom-most agent is a mixed-integer nonlinear problem, switch states, and DG outputs selected as decision variables. The immune clonal selection differential evolution algorithm is utilized to solve the mixed integer nonlinear optimization problem defined in the bottom-most surrogate [17]. Due to the presence of DG in the decision variables, some changes have been made to the algorithm. A considerable number of antibody groups are generated from genes that determine antibodies, which are used as the initial antibody groups for the differential evolution algorithm of immune clone selection as shown in Table 1.Table 1
The structure of chromosomes.
OS1-OSTP1aQ1aP1bQ1bP1cQ1c…PKaQKaPKbQKbPKcQKcThe definition of antibody genes is shown in Table1, and there are T +6 K genes in total: The first T gene represents the open switch; the following 6 K genes are the active and reactive power generated by the Kth DG unit. Antibody clustering, clonal expansion, and high-frequency mutation randomly select a gene i from T +6 K genes; get high-frequency mutation to randomly change a selected gene, and introduce new information into offspring. The selected gene could represent the open switch. Before evaluating fitness values for new antibody clusters, all duplicate antibodies were removed, and the viability of each progeny was assessed by sequentially examining system structural constraints and voltage/current constraints. Then, the fitness values of all viable progeny are calculated, and the best antibody is selected. Finally, the optimal topology of all the lowest-level subsystems and the optimal output of the DG unit are solved. The agents coordinate when necessary, passing the final result to the upper agent to activate the computation. The network reconstruction problem is solved separately and decomposed into different subproblems. A hierarchical structure is used to assign independent agents to subproblems to realize parallel computing. Each agent is made up of three units. The data unit collects local information and communicates with other agents. The computational unit realizes the improved heuristic algorithm to solve the local reconstruction problem, and the decision unit completes the coordination control. The calculation result of the lower agent is transmitted to the data unit of the upper agent. The final optimal configuration is completed by the collaboration between multiple agents.


### 4.4. Upper Layer Proxy Startup Issues

A heuristic algorithm based on branch-and-switch and single-loop optimization is adopted to solve the optimal topology of the upper-level subsystem with the known switching states and DG outputs solved in the lower-level agent. Then, the optimal topology of the whole distribution system and the actual output power of all DG units are obtained based on real-time data within the current time window. At this point, the distribution reconfiguration is completed, and the DG units are in the optimal operating state, which will remain unchanged until the next time window when the scheduling plan for the next cycle is re-evaluated.


## 4.1. Network Decomposition

Loops are more easily identified by the difference between closed contact switches and segmented opens for differentiation. In order to easily distinguish between tightly connected partitions and loosely connected partitions [16], the connectivity between the subdivisions is defined as follows:
(29)DA,B=NumberofpublicbusminTotalnumberofbusbarsinzoneA,TotalnumberofbusbarsinzoneB,IfDA,B=0, it means that partition A and partition B are independent of each other. Given a threshold δ, if, means that partition A and partition B are tightly connected and if DA,B≤δ means that partition A and partition B are loosely connected. However, if the threshold δ is set too small, the loosely connected partitions are not easily identified; if the threshold δ is set too large, the loosely connected partitions are easily expanded. Therefore, it is generally reasonable to choose a threshold δ of 0.2.The specific steps of network decomposition are as follows.The set of cut vertices of a connected graphG=V,E is a vertex set U⊆V satisfying the following conditions: G-U is not connected; G-K is connected, K⊂U; every vertex u of U in a connected graph G is connected to at least one vertex in every partition of G-U. A set of cut vertices is a cut vertex if there is only one vertex in the set.
(a)
The value of the modified adjacency matrixC is acquired. Let m be the number of loops, and the modified adjacency matrix is an m-dimensional square matrix(b)
If there are no partitions with a common bus other than the source node, these partitions are defined as basic partitions. Two basic partitions are connected to each other by a contact switch. The following steps are used for all basic partitions(c)
Draw the connectivity graph. In the connectivity graphG=V,E, the set of vertices V=v1,⋯,vnrepresents 1−n loops. The edge E connecting the vertices indicates that the loops are coupled to each other. The vertices are drawn in the studied region to represent the loops, and if the element C (i, j) is not zero, an edge is added between the vertices vi and vj. If loop i is loosely connected to loop j, then it is noted as L; otherwise, it is noted as T(d)
Be sure to check graph G for connectivity. If not, the isolated vertices need to be found, denoted asS. The corresponding loop of the isolated vertex contains the first member of the decomposition system(e)
Determine the cut vertex set of G-s. The search starts from the parent node, and the cut vertex set node is found from the child nodes associated with the parent node. If the cut vertex set is empty, the search moves to the vertex associated with the child node. If the cut vertex set is found, the cut vertex set is removed, and the associated edges define two separate components. Continue searching for each component’s cut vertex set until there are no more cut vertex sets. Each component represents the decomposed system, and the contact switch corresponding to the cut vertex set represents the decomposed system by the links between them(f)
In case two vertices are connected to the edge marked withL, their corresponding rings are decomposed into two subsystems. If two loosely connected vertices are tightly connected to another vertex, the contact switch corresponding to that vertex is an interconnection between the two subsystems. In the lowest level subsystem forming the basis of the entire system, which includes all busbars, loads, and DG units. Each higher-level subsystem consists of several subsystems representing the entire system


## 4.2. Applying a Multi-Agent Framework

Intelligent agents consisting of data units, computation units, and decision units need to be assigned to each subsystem for solving subproblems of the assigned subsystem and exchanging information with other agents. Since DG units may exist in the bottom most subsystem, the bottom most agent should be able to decide the optimal topology and the best DG output of the local system. The optimization problem for each bottom most agent layer is represented by Equation (29), and all variables are refined into local variables for its subsystem. Based on the decision plan solved by the lower layer agents, all upper layer agents ought to judge whether the common contact switch is off or not; hence, each embedded optimization problem is a pure reconstruction problem with the equation:
(30)minJ=PlossSsubsystem,Ssolved,PQsolved⋅ΔT.In the formula,Ssubsystem is the set of switching states to be solved; Ssolved is the set of switching states already solved by the lower layer agent; and PQsolved is the actual output power of the DG unit of the lowest layer agent.


## 4.3. Optimization Problems for the Lowest-Level Agent

The optimization problem defined in the bottom-most agent is a mixed-integer nonlinear problem, switch states, and DG outputs selected as decision variables. The immune clonal selection differential evolution algorithm is utilized to solve the mixed integer nonlinear optimization problem defined in the bottom-most surrogate [17]. Due to the presence of DG in the decision variables, some changes have been made to the algorithm. A considerable number of antibody groups are generated from genes that determine antibodies, which are used as the initial antibody groups for the differential evolution algorithm of immune clone selection as shown in Table 1.Table 1
The structure of chromosomes.
OS1-OSTP1aQ1aP1bQ1bP1cQ1c…PKaQKaPKbQKbPKcQKcThe definition of antibody genes is shown in Table1, and there are T +6 K genes in total: The first T gene represents the open switch; the following 6 K genes are the active and reactive power generated by the Kth DG unit. Antibody clustering, clonal expansion, and high-frequency mutation randomly select a gene i from T +6 K genes; get high-frequency mutation to randomly change a selected gene, and introduce new information into offspring. The selected gene could represent the open switch. Before evaluating fitness values for new antibody clusters, all duplicate antibodies were removed, and the viability of each progeny was assessed by sequentially examining system structural constraints and voltage/current constraints. Then, the fitness values of all viable progeny are calculated, and the best antibody is selected. Finally, the optimal topology of all the lowest-level subsystems and the optimal output of the DG unit are solved. The agents coordinate when necessary, passing the final result to the upper agent to activate the computation. The network reconstruction problem is solved separately and decomposed into different subproblems. A hierarchical structure is used to assign independent agents to subproblems to realize parallel computing. Each agent is made up of three units. The data unit collects local information and communicates with other agents. The computational unit realizes the improved heuristic algorithm to solve the local reconstruction problem, and the decision unit completes the coordination control. The calculation result of the lower agent is transmitted to the data unit of the upper agent. The final optimal configuration is completed by the collaboration between multiple agents.


## 4.4. Upper Layer Proxy Startup Issues

A heuristic algorithm based on branch-and-switch and single-loop optimization is adopted to solve the optimal topology of the upper-level subsystem with the known switching states and DG outputs solved in the lower-level agent. Then, the optimal topology of the whole distribution system and the actual output power of all DG units are obtained based on real-time data within the current time window. At this point, the distribution reconfiguration is completed, and the DG units are in the optimal operating state, which will remain unchanged until the next time window when the scheduling plan for the next cycle is re-evaluated.


## 5. Algorithm Simulation

The simulated system of the distribution network used in this paper is a 25-node unbalanced system [18] with unbalanced line impedance and load distribution as shown in Figure 3.The system has three tie switches, 25, 26 and 27, respectively. Under normal operation, the interconnection switch is open. The initial three-phase power loss is 450.38 kW, and the minimum voltage is 0.93 pu.Figure 3
The single-line diagram of the 25-bus unbalanced distribution system.In this paper, the load is modeled with a Gaussian mixture model (GMM) proposed in the literature [19], and the probability density function of the GMM is determined by the following equation:
(31)fz=∑i=1AMwiNμi,σi.In the formula, AM is the number of components in the mixture model,μi is the mean, σi is the standard deviation, and wi is the weight.A GMM is applied randomly to each load. So the actual load power is obtained by multiplying the initial value by the unit value generated by the GMM. Monte Carlo simulations were performed, and the results showed that buses 12, 13, 11, 10, 15, 17, 14, and 9 were consistently the most sensitive 8 buses out of all 200 samples. To demonstrate the effectiveness of power loss reduction by installing DG units on the most sensitive buses, a system with 8 DG units on each of the 8 buses mentioned above is installed. After the DG unit locations are determined, the optimal capacity of each DG unit can be solved. Three scenarios are tested with a capacity limit of 500 kW/500kVA per phase DG, and the proposed Newton algorithm converges quickly after 5 to 20 iterations for different scenarios. The optimal capacities of DG units for the three scenarios are given in Table2:
(1)
Install a DG unit at the most sensitive bus 12(2)
Install a DG unit at bus 12 and 13, respectively(3)
Install a DG unit at bus 12, 13 and 11, respectivelyTable 2
Optimal capacity of DG units for three scenarios.
ScenarioOptimal capacity regionInstallation sitePhase APhase BPhase CActive(kw)Reactive(KVar)Active(kw)Reactive(KVar)Active(kw)Reactive(KVar)112356.2362.4420.5368.3258.4189.3Three−phaseminimumvoltage=0.966,0.965,0.956puPower loss =224.2 kW; maximum voltage imbalance =2.236%212220.5158.2210.8172.4214.8165.913214.8156.7208.6163.2212.4162.5Three−phaseminimumvoltage=0.964,0.965,0.963puPowerloss=192.8kW; maximumvoltageimbalance=0.562%312136.5103.6145.2106.3148.6108.213135.2104.5143.8102.6146.3105.811152.6116.2160.5116.5162.6115.6Three−phaseminimumvoltage=0.963,0.965,0.965puPowerloss=186.2kW; maximumvoltageimbalance=0.358%Tests were conducted on several scenarios, and the results are shown in Table2.The results show that the optimal capacity of the DG units can be successfully solved while satisfying all system constraints. The integration of DG units helps to reduce power losses, increase voltage, and reduce voltage imbalance. The results of scenarios 2 and 3 are significantly better than those of scenario 1, indicating that the integration of multiple decentralized small-capacity DG units is more helpful than the integration of one large-capacity DG unit. Comparing the results from scenarios 2 and 3, the integration of two DG units already reduces the losses significantly, and adding a third DG unit does not help much. Therefore, the final decision was made to install two DG units at buses 12 and 13.The operation period is set to 1 h, and the starting point is 0 : 00. 4 sets of 24-hour load curves are input to the system as real-time load data as shown in Figure4.Figure 4
Four groups of load shapes.The optimal switching scheme is shown in Table3:Table 3
The optimal switching plan.
Time windowOpen switchesTime windowOpen switchesTime windowOpen switches0-1 h15, 17, 229-10 h15, 17, 2217-18 h15, 17, 222-3 h15, 17, 2210-11 h17, 24, 618-19 h15, 17, 223-4 h15, 17, 2211-12 h15, 17, 2219-20 h15, 17, 224-5 h15, 17, 2212-13 h17, 24, 2220-21 h15, 17, 225-6 h15, 17, 2213-14 h15, 17, 2221-22 h15, 17, 226-7 h15, 17, 2214-15 h17, 24, 2222-23 h17, 24, 67-8 h17, 24, 615-16 h17, 24, 623-24 h15, 17, 228-9 h15, 17, 2216-17 h15, 17, 22The following graphs give the actual output of the DG unit for 24 hours as shown in Figures5–9.Figure 5
Active output power of the DG unit at Bus-12.Figure 6
Reactive output power of the DG unit at Bus-12.Figure 7
Active output power of the DG unit at Bus-13.Figure 8
Reactive output power of the DG unit at Bus-13.Figure 9
System power loss within 24 hours.


## 6. Conclusion

A hierarchical decentralized distribution network reconfiguration method is proposed based on the characteristics of three-phase unbalanced distribution network reconfiguration. Network decomposition and multi-agent architecture are used to obtain the optimal reconfiguration scheme. In addition, although multiple agents are set up, the necessary information exchange between them is only the switching status of their subsystems, so the information transfer is low. The network reconfiguration algorithm has been applied in an unbalanced 25 node system. The simulation results show that the algorithm can reduce the power loss and the imbalance degree of the system.



---

*Source: 9906100-2022-08-18.xml*