Next Article in Journal
A Novel Air-Cooled Thermal Management Approach towards High-Power Lithium-Ion Capacitor Module for Electric Vehicles
Previous Article in Journal
Identifying Key Factors in Determining Disparities in Energy Consumption in China: A Household Level Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Reconfiguration Strategy for DC Distribution Network Fault Recovery Based on Hybrid Particle Swarm Optimization

1
College of Computer and Electrical Engineering, Hunan University of Arts and Science, Changde 415000, China
2
Key Laboratory of Hunan Province for Control Technology of Distributed Electric Propulsion Air Vehicle, Changde 415000, China
3
College of Electrical and Information Engineering, Hunan University, Changsha 410082, China
*
Author to whom correspondence should be addressed.
Energies 2021, 14(21), 7145; https://doi.org/10.3390/en14217145
Submission received: 8 September 2021 / Revised: 10 October 2021 / Accepted: 25 October 2021 / Published: 1 November 2021
(This article belongs to the Topic Power Distribution Systems)

Abstract

:
DC distribution network faults seriously affect the reliability of system power supply. Therefore, this paper proposes a fault recovery reconfiguration strategy for DC distribution networks, based on hybrid particle swarm optimization. The original particle swarm algorithm is improved by simplifying the distribution network structure, introducing Lévy Flight, and designing an adaptive coding strategy. First, the distribution network structure is equivalently simplified to reduce the problem dimensionality. Further, the generated branch groups are ensured to satisfy the radial constraints based on the adaptive solution strategy. Subsequently, Lévy flight is introduced to achieve intra-group optimality search for each branch group. The method is simulated in several distribution systems and analyzed in comparison with the particle swarm algorithm, genetic algorithm, and cuckoo algorithm. Finally, the results validate the accuracy and efficiency of the proposed method.

Graphical Abstract

1. Introduction

With the increasing demand for electricity, the size and complexity of modern distribution networks have increased significantly [1]. In practice, most power outages are due to damage caused by distribution network failures, which can have serious consequences for both utilities and customers [2]. At the same time, medium-voltage direct current (MVDC) distribution networks are becoming increasingly important in view of the large-scale application of DC technology [3,4,5]. Compared with the traditional AC distribution network, the flexible DC distribution network based on voltage source converter (VSC) technology has the advantages of easy access to distributed energy and multiple loads, and flexible control [6]. Therefore, to improve power supply reliability and customer satisfaction, it is necessary to study the fault recovery reconfiguration of the DC distribution network.
Reconfiguration of the distribution network aims to find the optimal combination of all switches in the distribution, and achieve objects as follows: (1) maximizing demand, (2) minimizing the number of switching operations, (3) giving priority to automatic switching operations over manual operations, and (4) giving priority to special loads [7]. Distribution network reconfiguration (DNR) is essentially a multi-objective nonlinear hybrid optimization problem, which is mathematically an NP-hard problem [8].
There are several common reconfiguration methods, including heuristic algorithms, artificial intelligence (AI), and mathematical modeling. Heuristic algorithm is an intuitive analysis algorithm. It iterates step-by-step, according to certain principles, to obtain satisfactory results, and mainly includes branch-and-swap and optimal flow methods. In [9], the branch switching pair that caused the most network loss reduction is selected for branch switching in a simple way. In [10], all switches are closed, and the switches with the least optimal flow are subsequently opened sequentially until the reconfiguration solution is obtained. In addition, inspired by the fractal theory for solving optimization problems, implementation of a stochastic fractal search (SFS) algorithm is proposed to solve the distribution network reconstruction problem [11]. Although the method of [11] is an improvement on [9,10], it still cannot guarantee that the solution is the global optimal solution. In addition, ref [12] applies a heuristic algorithm based on Harris Hawks optimization to DNR, which can ensure global convergence performance. Ref [13] uses the discrete network reconfiguration of the data set method, which can significantly improve the effectiveness of the distribution network. This type of method has good convergence performance, and can obtain a global optimal solution for single-objective optimization problems.
AI algorithms for solving reconstruction problems are mainly based on swarm optimization methods. In [14], a genetic algorithm (GA) based on a standard benchmark problem is used to solve single-objective optimization problems. PSO has the advantages of simple principle, robustness, and easy implementation [15]. A binary particle swarm algorithm (BPSO) is proposed in [16] to solve the DNR problem. A new optimization algorithm called the Salp swarm algorithm can effectively solve the optimization problem. The beetle antennae search (BAS) algorithm is an intelligent optimization algorithm, which has the advantages of a simple principle, fewer parameters, and less calculation required [17]. The cuckoo algorithm (CA) can effectively solve the optimization problem of reconstruction by simulating the parasitic brooding of some species of cuckoo, and using the related Lévy flight search mechanism [18,19]. The above methods have been successfully applied to the field of DNR with better results than heuristic algorithms. However, a common problem is that they tend to fall into premature maturity and still cannot guarantee global optimality.
In addition, the reconfiguration problem can also be solved using mathematical modeling. In [20], a mixed-integer linear programming (MILP) model was developed to solve the distribution network fault recovery–reconstruction problem. Based on MILP, a mixed-integer second-order cone programming formulation for service restoration in distributed generation distribution networks was proposed, which relaxes the original non-convex tidal equation to quadratic form [21]. Meanwhile, the problem of uncertainty for renewable energy or load can be solved in the mathematical model by stochastic optimization [22] or robust optimization [23] to obtain the theoretically optimal solution. Mathematical modeling can be proved by formulas, but solving non-convex problems takes too long to solve. Therefore, linearized mathematical models and convex transformation of the original non-convex problems are further difficulties to be overcome for this kind of problem.
AI is a stochastic class of algorithms, and this class of methods has been successfully applied to solve DNR. Since traversing the feasible solution space takes a lot of computational time, the efficiency of the algorithm becomes crucial. The key to improving the efficiency of the algorithm lies in the encoding method, and the corresponding evolutionary approach. In traditional distribution networks, most distribution systems are designed as weak meshes, but usually operate in a radial topology to efficiently coordinate their protection systems [24,25]. In [26], a loop-based coding method was proposed to disconnect a branch in an independent loop to satisfy radial constraints, and increase the proportion of feasible solutions. The resulting solutions are all feasible, but the repetition rate is too high, reducing evolutionary efficiency. In [27], a new graph theory-based method is proposed for restoring the distribution system after multiple simultaneous faults due to extreme weather conditions. All of the above different encoding methods and evolutionary strategies improve the encoding efficiency to a certain extent, but still need to carry out the feasibility testing process of the solution, and the computation time is still long.
In summary, the reconstruction method based on the stochastic algorithm can solve the optimal solution more efficiently when solving large-scale distribution network problems. However, the stochastic algorithm generates many invalid solutions, and needs to judge whether the generated solutions are valid, which reduces the evolutionary efficiency. Meanwhile, the stochastic algorithm falls into local optimum easily. Therefore, a reconfiguration strategy based on HPSO is proposed. The main contributions of this paper are summarized as follows:
(1)
Considering the topological characteristics of DC distribution, the network structure is equivalently simplified to improve the search efficiency in the DNR process;
(2)
An adaptive coding strategy is designed to make the generated solutions satisfy the topological constraints. This strategy can improve the evolution efficiency of particles in the particle search process without judging infeasible solutions;
(3)
To avoid the algorithm from falling into premature, the idea of Lévy flight (LF) is introduced to improve the global search ability and convergence speed of the particles.
The paper is organized as follows. The mathematical modeling of DNR is in Section 2. In Section 3, the simplification and coding method of the distribution network is proposed. In Section 4, an HPSO combining adaptive coding strategy and improved discrete particle swarm optimization (DPSO) is proposed. In Section 5, the proposed algorithm is compared with other algorithms such as GA, PSO, and CA to verify the efficiency and accuracy of HPSO. Section 6 summarizes the paper.

2. Modeling of DC Distribution Network Restoration Reconfiguration

2.1. Mathematical Model of Converters

Before calculating the tidal distribution of DC distribution network, it is necessary to establish an equivalent circuit for the converters commonly found in DC distribution networks, and obtain the mathematical model of a DC distribution network. The voltage source converter (VSC) has obvious advantages in tidal current reversal and control, compared with the current source converter (CSC). In Figure 1, Rc + jXc is the equivalent impedance of the converter, Rc is the equivalent resistance of the transformer and the converter reactor (Ω), Xc is the equivalent reactance of the converter (Ω), Ps and Qs represent the active (W) and reactive power (Var) of the AC bus, Pc and Qc represent the active (W) and reactive power (Var) of the AC side, Us and δs represent the line voltage amplitude (V) and phase angle (rad) of the AC bus, Uc and δc represent the line voltage amplitude (V) and phase angle (rad) of the AC side, and Pdc and Qdc are the active power (W) and reactive power (Var) of the DC side for the converter, respectively.
The active power and reactive power on the AC side satisfy the following relationship:
{ P s = U s 2 G c U s U c [ G c cos ( δ s δ c ) + B c sin ( δ s δ c ) ] Q s = U s 2 B c U s U c [ G c sin ( δ s δ c ) B c cos ( δ s δ c ) ]
Equation (2) can be approximated to derive the internal losses of the converter P V S C l o s s , where A, B, and C are the coefficients of the fitted quadratic functions for converter current and active losses, respectively.
P l o s s V S C = a + b I c + c I c 2
a = 6.62 S N 600 S B ,   b = 1.8 V N d 600 V B d ,   c = 3 ( V N d ) 2 S B 600 ( V B d ) 2 S N
where SN and V N d denote the rated capacity (VA) of the converter and the rated voltage (V) of the DC side; and SB and V B d denote the reference capacity (VA) of the system and the reference voltage (V) of the DC network.
In Figure 1, the power relationship between the DC side of the converter and the AC side is shown in the following equation:
P c = P l o s s V S C + P d c

2.2. The Calculation of Power Flow in Medium Voltage DC Distribution Network

The DC distribution network does not have reactive power and voltage phase angle, and each node contains only two variables: power, and voltage. Therefore, only the value of one of the two variables is needed to solve the power flow. The Newton–Raphson method is the basic algorithm for solving the power flow. To solve the power flow, it is only necessary to obtain the voltage values of each node in the network. The nodal conductance matrix G needs to be established first in the tide calculation. Since the DC system is a purely resistive line and does not contain reactance, the conductance matrix is composed entirely of conductance, and the nodal current matrix I and the nodal voltage matrix U are related as follows:
I = G U
The amount of active power and power imbalance injected at node i on the DC side is shown as follows:
P i = U i I i = i = 1 n G i j U i U j
Δ P i = P i s P i = P i s j = 1 n G i j U i U j
where Ii is the node i injection current (A); Ui is the node i voltage (V); Gij is the conductance value (S) between node i and node j; and Pis denotes the fixed P node injection power (W).
The following equation is the formula for solving each element of the Jacobi matrix:
J i , j = Δ P i U j = { G i j U i , i j G i i U i k = 1 n G i k U k , i = j
The voltage adjustment value is obtained by using the Jacobi matrix, and the power flow of the distribution network can be obtained by iteratively correcting the voltage value.

2.3. Mathematical Model for Recovery Reconfiguration

DNR is a process of changing the switching state so that as much load as possible is transferred from the non-faulty zone to the normal feeder. DNR is not only a mixed-integer nonlinear optimization problem, but also a high-dimensional discrete non-convex optimization problem, making the optimization problem an NP-hard problem. The conventional solver takes too long to solve, and even the final solution cannot be obtained. Therefore, AI algorithms can be used to solve the problem. The main consideration is system economy, and the objective is to minimize network loss after reconfiguration.
The constraints that should be satisfied for fault recovery reconfiguration mainly include the tidal current constraint of the distribution network, the upper and lower voltage constraints of the nodes, the upper and lower branch current constraints, and the topological constraints of the distribution network structure:
o b j : min f = j = 1 N b k j I j 2 r j s . t .   A i = I U i , min   U i U i , max I i I i , max     g G
where Nb denotes the number of distribution network branches after reconfiguration; kj denotes the operation status of the first branch; 1 indicates normal operation and 0 indicates withdrawal; Ij is the current of the corresponding branch (A); rj is the impedance of the corresponding branch (Ω); A is the node-branch correlation matrix; i is the current vector of all branches; I is the current injection vector of all nodes; Ui is the actual voltage of node i; Ui,max and Ui,min denote the allowed upper (V) and lower voltage (V) of node i; Ii is the actual current (A) flowing through the branch; Ii,max is the maximum current (A); and g is the current network structure.
The decision variable is the branch switch number of the power distribution system, which is an integer decision variable. The number of tie switches in the distribution system is equal to the number of decision variables. If the IEEE33 node is used as an example, there are five decision variables for DNR. Among the constraints, active power balance constraints and reactive power balance constraints are equality constraints, and node voltage constraints and branch current constraints are inequality constraints.

3. Simplification and Coding of Power Distribution Networks

3.1. Simplification of Power Distribution Networks

A characteristic of smart algorithm solutions is randomness, which allows the algorithm to traverse the entire solution space, but it also encounters situations where the generated solutions are infeasible. Due to the existence of radial constraints in the distribution network, the algorithm may generate infeasible solutions during all iterations, and the decoded distribution network structure does not satisfy the radial constraints. The existence of infeasible solutions seriously reduces the search efficiency of distribution network reconstruction, so the generation of infeasible solutions should be avoided as much as possible during the iteration of the algorithm. Considering the large number of variables in the actual distribution network, the distribution network structure should also be simplified equivalently.
As shown in Figure 2, the following equivalent simplifications are made to the distribution network before fault recovery reconfiguration:
(1)
Some branches that already exist in radial branches need to be ignored. For example, the branch between node 1 and node 2, which is involved in splitting, will form an island that does not satisfy the radial constraint. Therefore, the first step should be to eliminate the branches that are not related to the ring network first;
(2)
Branches that have the same effect need to be disconnected. For example, disconnecting any two branches between node 3 and node 6 will have the same effect for disconnected loops, although the tide of the system will change accordingly. Therefore, in the second step, the concept of degree in graph theory is introduced to eliminate the nodes with 2°. The remaining equivalent branches are merged to form a set of branches with the two end nodes as branches.
As can be seen in the final simplified diagram generated in Figure 2, with all tie-switches closed, the original IEEE33 distribution system is reduced to an equivalent diagram containing 12 branched groups of 5 loops (L1–L5).

3.2. Coding Method

How the search space is defined is the key to determine the efficiency of the algorithm. For the simplified topology diagram generated at the end of Figure 2, a decision variable is used to represent a branch group, and the number of branch groups is the dimensionality of the particles. If the value is 0, it means that all the branches contained in the corresponding branch group are in the on state, and the number of removed branches is 0. If the value is a positive integer, the integer indicates the numbered branch corresponding to the removed branch. Therefore, the branch group generation process requires numbering of each branch. The upper value of each dimension is the number of branches in each branch group, and the lower limit is 0. In this coding method, the maximum number of branches broken in each branch group is guaranteed to be 1, which reduces the generation of infeasible solutions and improves the search efficiency.
The branch group matrix is defined with the branch group number as the row number and the internal branch number of each branch group as the internal element. For the initial state of the contact switch in Figure 2, the disconnected switch is numbered [3 3 3 4 3 5 3 6 3 7], and the decision variables of this network are coded [0 0 0 0 0 0 0 0 1 2 4 0 5 1].
After the initialization of the particle swarm is completed, the evolution of the particles is divided into two phases. The first stage uses an adaptive solution-based BPSO to ensure that the generated solution satisfies the topological constraints of the distribution network. The second stage uses the improved DPSO to select a disconnected switch in each branch group.

4. HPSO Algorithm

4.1. Branched Group Selection Optimization Based on Adaptive BPSO

The BPSO sets the individual pid and global pgd of the particle and the position vid of each dimension to 0 or 1, depending only on the magnitude of the particle velocity vid. If the velocity is larger, the chance that the corresponding particle dimension xid is 1 will be larger. Conversely, if the velocity is smaller, the chance that the current particle dimension xid is 0 will be larger. The update formula for the position update expression sigmod function of the BPSO algorithm is shown below:
{ v i d t + 1 = ω v i d t + c 1 r ( p i d x i d t ) + c 2 r ( p g d x i d t ) { x i d t + 1 = 1 ,                 r < s i g m o i d ( v i d t + 1 ) x i d t + 1 = 0 ,                 r s i g m o i d ( v i d t + 1 )
s i g m o i d ( x ) = 1 e x
where r is the random value on the interval [0, 1], ω is the inertia weight, c1 and c2 denote the two acceleration factors, vt+1 id denotes the velocity of particle i in the d-th dimension of the second iteration, and xt+1 id denotes the position of particle i in the d-th dimension of the second iteration.
The necessary condition for DNR is that the number of disconnected tributary switches must be the number of tie-switches, since it cannot be arbitrarily set to 0 or 1 in each dimension by the binary algorithm, and needs to be limited. According to the number of loops in the distribution network, disconnecting a tributary group number in each loop in turn can satisfy the necessary condition for DNR, but it is still not a sufficient condition. Here, according to the distribution network loops, a strategy to improve BPSO is proposed for adaptive solving to ensure that the generated solutions are all feasible solutions. Algorithm 1 shows the adaptive coding algorithm.
Algorithm 1: Adaptive coding algorithm.
  Input: loop branch group and nodal branch incidence matrix
  Output: s e l e c t e d _ g r o u p
1Initialization: c u r r e n t _ m e m b e r = 1 ; s e l e c t e d _ g r o u p = ; μ is the sigmoid parameter of BPSO;
2while s e l e c t e d _ b r a n c h n u m b e r _ l o o p s do
3  while c u r r e n t _ m e m b e r μ do
4    Pick a random value r uniformly from [ 0 , 1 ]
5    Set i = 1 ;
6    while R o u l e t t e _ p r o b a b i l i t y [ i ] r do
7     Set i = i + 1 ;
8    end
9    if n o d e _ d e g r e e [ i ] = 1 then
10      s e l e c t e d _ g r o u p = s e l e c t e d _ g r o u p c o n n e c t e d _ b r a n c h [ i ] ;
11    else
12     go to step 6;
13    end
14    Update n o d e _ d e g r e e [ j ] , j S n o d e ;
15    Set c u r r e n t _ m e m b e r = c u r r e n t _ m e m b e r + 1 ;
16  end
17  Derive s e l e c t e d _ g r o u p ;
18end
By choosing the disconnected branching group numbers sequentially through the loop number, the dimensionality of the decision variables can be guaranteed to be constant. The adaptive solution strategy for distribution network loops with branching degree can obtain feasible solutions that satisfy the radial constraint of DNR. For the DNR model, the stochastic algorithm can traverse the whole solution space. Based on the adaptive solution strategy, the solutions generated in are feasible solutions.

4.2. LF-Based Improved DPSO for Intra-Group Optimization

After determining the group of branches to be operated, an internal search for the best solution is required, at which point a feasible solution is generated using the encoding method. Since the decision variable at this point is an integer, it is necessary to use an algorithm that can handle decision variables of integer type. We define a random rounding method rand(N): with an integer N, if N ≥ 0, then a random integer value is selected from the interval [0, N]; if N ≤ 0, then a random integer value is selected from the interval [N, 0]. We can obtain the speed and position iteration formula of the DPSO as:
v i d t + 1 = v i d t + r a n d ( c 1 ( p i d x i d t ) ) + r a n d ( c 2 ( p g d x i d t ) )
x i d t + 1 = x i d t + v i d t
LF is a random wandering strategy obeying the Lévy distribution, and the position iterative update equations are:
x i d t + 1 = x i d t + α L e v y ( λ )
where is the point-to-point multiplication, α denotes the step control quantity, and L e v y ( λ ) is the random search path. In order to simplify the distribution function, the probability density function can be obtained by Fourier transform
L e v y u = t β , 1 β 3
The Mantegna algorithm, where the step size s can be calculated using two variables U and V that obey a Gaussian distribution:
s = U | V | 1 / λ
U ~ N ( 0 , σ 2 ) , V ~ N ( 0 , 1 )
σ 2 = [ Γ ( 1 + λ ) λ Γ ( ( 1 + λ ) / 2 ) · sin ( π λ / 2 ) 2 ( λ 1 ) / 2 ] 1 / λ
where U ~ N ( 0 , σ 2 ) indicates that the sample follows a Gaussian normal distribution with mean 0 and variance σ 2 ; V ~ N ( 0 , 1 ) indicates that the sample follows a standard normal distribution.
The improved DPSO combined with LF can broaden the spatial search capability of particles, and increase the diversity of particle search. Using the previously described encoding method to update the particles, the iterative formula of the improved discrete PSO based on LF for intra-group optimization can be obtained as shown below.
v i d t + 1 = { v i d t ,      ( x i d t = 0 )         v i d t ,      ( x i d t 0 ,   p i d t = p g d t = 0 )   v i d t + r a n d ( c 1 ( p i d t x i d t ) ) + r a n d ( c 2 ( p g d t x i d t ) ) + r o u n d ( l e v y ( x i d t ) ) , ( x i d t 0 , p i d t 0 , p g d t 0 )   v i d t + r a n d ( v max d ) + r a n d ( c 2 ( p g d t x i d t ) ) + r o u n d ( l e v y ( x i d t ) ) , ( x i d t 0 , p i d t = 0 , p g d t 0 )   v i d t + r a n d ( c 1 ( p i d x i d t ) ) + r a n d ( v max d ) + r o u n d ( l e v y ( x i d t ) ) , ( x i d t 0 , p i d t 0 , p g d t = 0 )  
x i d t + 1 = { r ( 1 , 2 , , L d ) , x i d t = 0 r ( 1 , 2 , , L d ) , x i d t 0 , p i d t = p g d t = 0 x i d t + v i d t + 1 ,     other
where, r ( 1   , 2   , , L d ) denotes a random value in the 1, 2, ⋯, Ld array, Ld is the number of branches contained in the d branch group, and the round symbol indicates rounding to the nearest integer.
If the d-th dimension of particle i is selected for disconnection in step 1, the optimization process within the branch group in step 2 is divided into the following cases.
(1)
If the particle i being optimized is set to 0, it means that the branch group in this dimension has no previously selected branches. Therefore, the individual and global extremes have no guiding effect on this particle and the velocity of this particle does not change in this dimension and retains the last value, see the first case of (19). However, since this iteration needs to find the internal optimization of the branching group, the position of this dimension should be determined by randomly selecting a branching path, see the first case of (20);
(2)
If the value of the optimized particle is not set to 0, it is further divided into the following cases:
If the individual extreme value and global extreme value velocity is 0, it cannot affect the particle velocity in this dimension. There is no guidance for non 0 particles, the velocity does not change, the position from the branch group randomly chooses a path, see the second case of (19) and (20);
If the individual extreme and the global extreme are all non 0, both are guided at this time, and can be updated according to the particle iteration formula after joining the LF, see the third case of (19) and (20);
If one of the individual and one of the global poles is 0 and one is not 0, the poles in the 0 state cannot be guided, and the poles in the non 0 state are guided to the velocity. The corresponding item is changed to a random value from the velocity space, LF is added to enhance the particle search ability, see the fourth and fifth cases of (19) and the 3rd case of (20).

5. Case Study

The structure of the original IEEE 33-node arithmetic example is shown in the original diagram in Figure 2. Since the original example is an AC system, a VSC converter is added between node 1 and node 2, and the reactive load of the node load and the reactance of the branch circuits are ignored to make it a DC distribution network. To prove the effectiveness of the algorithm, we take the same approach for the IEEE 69-node improvement as a test case. In the modified IEEE 33-node, branch 7 is set as a faulty branch; in the modified IEEE 69-node example, branch 14 is set as a faulty branch. Since the current study is a static reconstruction problem, the test is focused on a single time section.
Table 1 shows the parameters of different algorithms. This paper mainly considers the proposed HPSO in comparison with PSO, CA, GA, and BAS. The population size and maximum number of iterations are kept consistent for all algorithms. In addition to the variables that have already appeared, it is necessary to account for the added parameters. In GA, pc, pm, pw, and pn denote the crossover probability, the exchange variation probability, the inverse variation probability, and the addition probability, respectively. The η in BAS denotes the coefficient of variable step size.

5.1. Comparison of Differen Metaheuristic t Algorithms

The metaheuristic algorithm seeks the optimal solution mainly through iterations of pseudo-probabilities, leading to differences in the results of each run. Therefore, each algorithm is set to run 50× to compare the mean value of the optimal solution and running time. In addition, the number of successful convergences to the global optimal solution is used as a comparison metric.
Table 2 shows the comparison results of different algorithms. Among all algorithms, HPSO has the highest number of successful solutions, which results in the smallest mean value of all solutions. The convergence performance of HPSO is faster than other algorithms, because the adaptive solution of HPSO reduces the judgment process of feasible solutions, which further improves the operation efficiency. PSO and BAS have simpler algorithm principles and faster iterations, thus the running time is smaller than the other two algorithms. In addition, the network architecture of 69 nodes is more complex than that of 33 nodes, thus the algorithms take longer to solve 69 nodes and have fewer successful solutions.
Figure 3 shows the comparison of evolution curves of different algorithms. Among the results of the runs for the improved IEEE 69-node, the one with the fastest convergence is selected for comparison. In Figure 3, the maximum number of iterations for each algorithm is 50. In the enlarged local plot, HPSO converges to the global optimal solution most quickly after 9 iterations. Although the running time of GA is longer, the optimal solution can be obtained in fewer iterations than other methods, and it only took 12 iterations to converge successfully. Therefore, the fastest number of iterations and running time are a pair of contradictory metrics in the DNR, and HPSO is highly compatible with them.

5.2. Comparison of Different Algorithms

HPSO is divided into two main parts: branched group selection optimization based on adaptive solving, and the internal optimization of branched groups based on LF. In order to reflect the advantages of HPSO, the following three algorithms are set up separately for comparison:
Improved binary particle swarm optimization (IBPSO):The branch group selection optimization uses the improved BPSO algorithm based on adaptive coding strategy proposed in this paper, and the conventional DPSO is used for the internal optimization of the branch group.
Improved discrete particle swarm optimization (IDPSO): The internal optimization of the branch group adopts the improved DPSO based on LF proposed in this paper, and the optimization of branch group selection uses the conventional BPSO.
HPSO: The branch group selection optimization adopts the improved BPSO based on adaptive solution proposed in this paper, and the internal optimization of the branch group adopts the improved DPSO based on LF.
Table 3 shows the results of different algorithms. On the one hand, IBPSO has the fastest running speed, since LF is not included in the intra-branch group search, which saves computation time. On the other hand, the lack of LF also reduces the diversity of particles, which makes IBPSO the weakest in global search. IDPSO has the same global search ability asHPSO, which is due to the inclusion of the LF strategy in the intra-branch group search. The particle search space becomes wider, which improves the population’s merit-seeking ability. However, the average running time is longer than that ofHPSO. The lack of adaptive solution strategy causes IDPSO to require the inclusion of infeasible solution judgment, which increases the total running time of the intelligent algorithm. HPSO combines the advantages of both. It can not only ensure the global merit-seeking ability of the population by improving DPSO, but also improve the algorithm’s operation efficiency by adaptive solution method. In summary, HPSOhas a better global search capability than Algorithm 1, and also can converge faster than IDPSO. The test verifies the effectiveness of the proposed algorithm HPSO.
Figure 4 shows the comparison of evolution curves of different algorithms. Three algorithms obtained the global optimal solution. HPSO achieves convergence the fastest and obtains the global optimal solution at the ninth iteration. The convergence performance of IDPSO is basically the same as that of HPSO, but the convergence speed of IBPSO is much smaller than those of IDPSO andHPSO. Therefore, the convergence performance of Algorithm 1 is weaker than those of the other algorithms, which proves the necessity of the intra-branch group optimization search strategy.

5.3. Validity Analysis of Conclusions

Table 4 compares the network losses and minimum voltage of the distribution network at different periods. Figure 5 shows the voltage comparison of each node before and after DNR. To verify the effectiveness of the reconfiguration scheme, a comparative analysis is performed for the distribution network before and after reconfiguration.
Both the network loss and voltage excursion of the reconfigured system were significantly reduced. In the 33-node case, the network loss value is reduced by 34.4%, and the minimum voltage of the distribution network is increased by 3.1%, which improves the economic performance of the system. In the 66-node case, the network loss value is reduced by 55.8%, the system loss is reduced, and the minimum voltage of the distribution network is increased by 4.0%. In Figure 5, the voltage magnitudes after reconfiguration are all improved, which improves the reliability of the system. The above diagram proves the effectiveness of the HPSO, which can be successfully applied to practical engineering.

6. Conclusions

In this paper, a two-stage HPSO method for DC distribution networks is proposed. The first stage uses an improved BPSO based on adaptive solution strategy for inter-branch group search, and the second stage uses an improved DPSO based on LF idea for intra-branch group optimization. The combination of the two stages can effectively solve the DC distribution network fault recovery reconfiguration problem. The established model has the following advantages:
(1)
The improved BPSO based on the adaptive solution strategy can remove the judgment process of infeasible solutions and improve the efficiency of population evolution;
(2)
The improved DPSO based on LF idea can avoid the population from getting into local optimum, and find the global optimal solution faster;
(3)
The equivalent simplification of the distribution network reduces the problem dimension, and the special coding can further improve the computational efficiency.
Therefore, the method in this paper is more efficient and reliable than the existing methods, which is a guideline for the development and construction of DC distribution networks. However, the increasing penetration of new energy in the distribution network makes its impact on the distribution network increasingly obvious. Therefore, considering the connection of renewable energy sources and analyzing the uncertainty will be our further work.

Author Contributions

Conceptualization, methodology, validation, software, data curation, project administration, investigation, funding acquisition, formal analysis, and writing—original draft preparation, M.Y.; software, resources, visualization, J.L. (Jianqi Li); writing—review and editing, J.L. (Jianying Li), and X.Y.; project administration, funding acquisition, J.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Program of Natural Science Foundation of Hunan Province (Grant no. 2019jj40200), the Program of Natural Science Foundation of Hunan Province (Grant no. 2019jj6002), the Program of Natural Science Foundation of Hunan Province (Grant no. 2020jj6061), the Research Foundation of Education Bureau of Hunan Province, China (Grant no. 19B393), the Scientific And Technological Innovation And Development Project of Changde District (Grant no. 2020C083), and the Science And Technology Innovation Program of Hunan Province (Grant no. 2021GK2010).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ACAlternating current
AIArtificial intelligence
BPSOBinary particle swarm algorithm
BASBeetle Antennae search
CACuckoo algorithm
CSCCurrent source converter
DPSODiscrete particle swarm algorithm
DNRDistribution network reconfiguration
DCDirect current
GAGenetic algorithm
HPSOHybrid particle swarm algorithm
LFLévy flight
MILPMixed-integer linear programming
PSOParticle swarm algorithm
SFSStochastic fractal search
VSCVoltage source converter
Nomenclature
A. Mathematical model of converters
Ps, QsActive power and reactive power on the AC bus
Pc, QcActive power and reactive power on the AC side
Pdc, QdcActive power and reactive power on the DC side
Us, δsVoltage amplitude and phase angle of the AC bus
Uc, δcVoltage amplitude and phase angle of the AC side
P l o s s V S C Internal losses of the converter
a, b, cThe respective coefficients of the fitted quadratic functions for converter current and active losses
SN, V N d The rated capacity of the converter and the rated voltage of the DC side
SB, V B d Reference capacity and the reference voltage
I, G, UNodal current matrix, nodal conductance matrix, nodal voltage matrix
kjOperation status
A, iNode-branch correlation matrix and current vector
Ui,max, Ui,minUpper and lower voltage
Ii,maxUpper current
B. HPSO algorithm
Pid, PgdGlobal optimal solution and individual optimal solution
vid, xidParticle speed and particle position
c1, c2Inertial System of Evolutionary Algorithm
αStep control quantity
sStep size
U ~ N ( 0 , σ 2 ) Gaussian normal distribution
V ~ N ( 0 , 1 ) Standard normal distribution
C. Other Symbols
Lévy (·)The random search path
Point-to-point multiplication
rand (·)Random rounding method
round (·)Rounding to the nearest integer
r (·)Random value
Sigmoid (·)Sigmoid function

References

  1. Li, J.; Li, D.; Zheng, Y.; Yao, Y.; Tang, Y. Unified Modeling of Regionally Integrated Energy System and Application to Optimization. Int. J. Electr. Power Energy Syst. 2022, 134, 107377. [Google Scholar] [CrossRef]
  2. Kumar, Y.; Das, B.; Sharma, J. Multiobjective, Multiconstraint Service Restoration of Electric Power Distribution System with Priority Customers. IEEE Trans. Power Deliv. 2008, 23, 261–270. [Google Scholar] [CrossRef]
  3. Chen, T.-H.; Huang, W.-T.; Gu, J.-C.; Pu, G.-C.; Hsu, Y.-F.; Guo, T.-Y. Feasibility Study of Upgrading Primary Feeders from Radial and Open-Loop to Normally Closed-Loop Arrangement. IEEE Trans. Power Syst. 2004, 19, 1308–1316. [Google Scholar] [CrossRef] [Green Version]
  4. Li, Y.; He, L.; Liu, F.; Li, C.; Cao, Y.; Shahidehpour, M. Flexible Voltage Control Strategy Considering Distributed Energy Storages for DC Distribution Network. IEEE Trans. Smart Grid 2019, 10, 163–172. [Google Scholar] [CrossRef]
  5. Ma, J.; Zhu, M.; Cai, X.; Li, Y.W. DC Substation for DC Grid—Part I: Comparative Evaluation of DC Substation Configurations. IEEE Trans. Power Electron. 2019, 34, 9719–9731. [Google Scholar] [CrossRef]
  6. Li, X.; Zhu, M.; Li, Y.; Cai, X. Cascaded MVDC Integration Interface for Multiple DERs with Enhanced Wide-Range Operation Capability: Concepts and Small-Signal Analysis. IEEE Trans. Power Electron. 2020, 35, 1182–1188. [Google Scholar] [CrossRef]
  7. Debnath, S.; Qin, J.; Bahrani, B.; Saeedifard, M.; Barbosa, P. Operation, Control, and Applications of the Modular Multilevel Converter: A Review. IEEE Trans. Power Electron. 2015, 30, 37–53. [Google Scholar] [CrossRef]
  8. Romero, R.; Franco, J.F.; Leão, F.B.; Rider, M.J.; de Souza, E.S. A New Mathematical Model for the Restoration Problem in Balanced Radial Distribution Systems. IEEE Trans. Power Syst. 2016, 31, 1259–1268. [Google Scholar] [CrossRef] [Green Version]
  9. Khushalani, S.; Solanki, J.M.; Schulz, N.N. Development of Three-Phase Unbalanced Power Flow Using PV and PQ Models for Distributed Generation and Study of the Impact of DG Models. IEEE Trans. Power Syst. 2007, 22, 1019–1025. [Google Scholar] [CrossRef]
  10. Baran, M.E.; Wu, F.F. Optimal Capacitor Placement on Radial Distribution Systems. IEEE Trans. Power Deliv. 1989, 4, 725–734. [Google Scholar] [CrossRef]
  11. Shirmohammadi, D.; Hong, H.W. Reconfiguration of Electric Distribution Networks for Resistive Line Losses Reduction. IEEE Trans. Power Deliv. 1989, 4, 1492–1498. [Google Scholar] [CrossRef]
  12. Helmi, A.M.; Carli, R.; Dotoli, M.; Ramadan, H.S. Efficient and Sustainable Reconfiguration of Distribution Networks via Metaheuristic Optimization. IEEE Trans. Autom. Sci. Eng. 2021, 1–17. [Google Scholar] [CrossRef]
  13. Muhammad, M.A.; Mokhlis, H.; Naidu, K.; Amin, A.; Franco, J.F.; Othman, M. Distribution Network Planning Enhancement via Network Reconfiguration and DG Integration Using Dataset Approach and Water Cycle Algorithm. J. Mod. Power Syst. Clean Energy 2020, 8, 86–93. [Google Scholar] [CrossRef]
  14. Tran, T.T.; Truong, K.H.; Vo, D.N. Stochastic Fractal Search Algorithm for Reconfiguration of Distribution Networks with Distributed Generations. Ain Shams Eng. J. 2020, 11, 389–407. [Google Scholar] [CrossRef]
  15. Song, C.; Luo, Q.; Shi, F. Genetic Algorithms for Optimization of Complex Nonlinear System. In Proceedings of the 2008 International Conference on Computer Science and Software Engineering, Wuhan, China, 12–14 December 2008; Volume 1, pp. 378–381. [Google Scholar]
  16. Kennedy, J.; Eberhart, R. Particle Swarm Optimization. In Proceedings of the ICNN’95–International Conference on Neural Networks, Perth, Australia, 27 November–1 December 1995; Volume 4, pp. 1942–1948. [Google Scholar]
  17. Wang, J.; Chen, H. BSAS: Beetle Swarm Antennae Search Algorithm for Optimization Problems. arXiv 2018, arXiv:1807.10470. [Google Scholar]
  18. Rajabioun, R. Cuckoo Optimization Algorithm. Appl. Soft Comput. 2011, 11, 5508–5518. [Google Scholar] [CrossRef]
  19. Gao, S.; Gao, Y.; Zhang, Y.; Li, T. Adaptive Cuckoo Algorithm with Multiple Search Strategies. Appl. Soft Comput. 2021, 106, 107181. [Google Scholar] [CrossRef]
  20. Nguyen, T.T.; Nguyen, T.T. An Improved Cuckoo Search Algorithm for the Problem of Electric Distribution Network Reconfiguration. Appl. Soft Comput. 2019, 84, 105720. [Google Scholar] [CrossRef]
  21. Souza, E.S.; Romero, R.; Franco, J.F. Restoration of Electrical Distribution Systems Using a Relaxed Mathematical Model. J. Control. Autom. Electr. Syst. 2018, 29, 259–269. [Google Scholar] [CrossRef] [Green Version]
  22. Scarabaggio, P.; Grammatico, S.; Carli, R.; Dotoli, M. Distributed Demand Side Management with Stochastic Wind Power Forecasting. IEEE Trans. Control. Syst. Technol. 2021, 1–16. [Google Scholar] [CrossRef]
  23. Raposo, A.A.M.; Rodrigues, A.B.; da Guia da Silva, M. Robust Meter Placement for State Estimation Considering Distribution Network Reconfiguration for Annual Energy Loss Reduction. Electr. Power Syst. Res. 2020, 182, 106233. [Google Scholar] [CrossRef]
  24. Li, Y.; Xiao, J.; Chen, C.; Tan, Y.; Cao, Y. Service Restoration Model with Mixed-Integer Second-Order Cone Programming for Distribution Network with Distributed Generations. IEEE Trans. Smart Grid 2019, 10, 4138–4150. [Google Scholar] [CrossRef]
  25. De Oliveira, L.W.; Seta, F.D.S.; de Oliveira, E.J. Optimal Reconfiguration of Distribution Systems with Representation of Uncertainties through Interval Analysis. Int. J. Electr. Power Energy Syst. 2016, 83, 382–391. [Google Scholar] [CrossRef]
  26. Fathabadi, H. Power Distribution Network Reconfiguration for Power Loss Minimization Using Novel Dynamic Fuzzy C-Means (DFCM) Clustering Based ANN Approach. Int. J. Electr. Power Energy Syst. 2016, 78, 96–107. [Google Scholar] [CrossRef]
  27. De Macedo Braz, H.D.; de Souza, B.A. Distribution Network Reconfiguration Using Genetic Algorithms with Sequential Encoding: Subtractive and Additive Approaches. IEEE Trans. Power Syst. 2011, 26, 582–593. [Google Scholar] [CrossRef]
  28. Whitley, D. A Genetic Algorithm Tutorial. Stat. Comput. 1994, 4, 65–85. [Google Scholar] [CrossRef]
Figure 1. Equivalent model of VSC.
Figure 1. Equivalent model of VSC.
Energies 14 07145 g001
Figure 2. Simplified diagram of distribution network.
Figure 2. Simplified diagram of distribution network.
Energies 14 07145 g002
Figure 3. Comparison of evolution curves of different metaheuristic algorithms.
Figure 3. Comparison of evolution curves of different metaheuristic algorithms.
Energies 14 07145 g003
Figure 4. Comparison of evolution curves of different algorithms.
Figure 4. Comparison of evolution curves of different algorithms.
Energies 14 07145 g004
Figure 5. Voltage comparison of each node before and after DNR. (a) Improved IEEE 33 nodes; (b) improved IEEE 69 nodes.
Figure 5. Voltage comparison of each node before and after DNR. (a) Improved IEEE 33 nodes; (b) improved IEEE 69 nodes.
Energies 14 07145 g005
Table 1. Parameters of different algorithms.
Table 1. Parameters of different algorithms.
AlgorithmParameter Value
PSO M a x I t e r = 50 , S w a r m S i z e = 20 , c 1 = 2 , c 2 = 2 , v max = 4 , v min = 4
CA M a x I t e r = 50 , S w a r m S i z e = 20 , α = 0.75 , β = 1.5 , σ = 0.70
GA M a x I t e r = 50 , S w a r m S i z e = 20 , p c = 0.7 , p m = 0.8 , p w = 0.3 , p n = 0.3
BAS M a x I t e r = 50 , S w a r m S i z e = 20 , η = 0.95
HPSO M a x I t e r = 50 , S w a r m S i z e = 20 , c 1 = 2 , c 2 = 2 , v max = 4 , v min = 4
Table 2. Results of different metaheuristic algorithms.
Table 2. Results of different metaheuristic algorithms.
CasesItemPSO [16]CA [18]GA [28]BAS [17]HPSO
Modified IEEE 33-nodeAverage of all solutions/kW89.635290.611289.471589.311688.8219
Average running time/s8.273111.160513.15678.97318.0690
Number of successes3932414549
Modified IEEE 69-nodeAverage of all solutions/kW65.367268.019064.103765.478663.4467
Average running time/s11.635814.743218.762312.455210.9651
Number of successes3225343138
Table 3. Results of different algorithms.
Table 3. Results of different algorithms.
CasesItemIBPSO IDPSO HPSO
Modified IEEE 33-nodeAverage of all solutions/kW91.537489.521488.8219
Average running time/s6.75289.35668.0690
Number of successes414749
Modified IEEE 69-nodeAverage of all solutions/kW65.063763.861263.4467
Average running time/s9.374613.847310.9651
Number of successes273938
Table 4. Comparison before and after DNR.
Table 4. Comparison before and after DNR.
CasesItemClosed SwitchNetwork Loss/kWMinimum Voltage/pu
Modified IEEE 33-nodeNormal operation33, 34, 35, 36, 37135.27600.9338
After reconfiguration7, 9, 14, 32, 3788.81140.9629
Improved IEEE 69 caseNormal operation69, 70, 71, 72, 73143.42780.9320
After reconfiguration14, 55, 61, 69, 7063.42630.9695
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yang, M.; Li, J.; Li, J.; Yuan, X.; Xu, J. Reconfiguration Strategy for DC Distribution Network Fault Recovery Based on Hybrid Particle Swarm Optimization. Energies 2021, 14, 7145. https://doi.org/10.3390/en14217145

AMA Style

Yang M, Li J, Li J, Yuan X, Xu J. Reconfiguration Strategy for DC Distribution Network Fault Recovery Based on Hybrid Particle Swarm Optimization. Energies. 2021; 14(21):7145. https://doi.org/10.3390/en14217145

Chicago/Turabian Style

Yang, Minsheng, Jianqi Li, Jianying Li, Xiaofang Yuan, and Jiazhu Xu. 2021. "Reconfiguration Strategy for DC Distribution Network Fault Recovery Based on Hybrid Particle Swarm Optimization" Energies 14, no. 21: 7145. https://doi.org/10.3390/en14217145

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop