Next Article in Journal
On the Error Estimation of the FEM for the Nikol’skij-Lizorkin Problem with Degeneracy in the Lebesgue Space
Next Article in Special Issue
Variable Matrix-Type Step-Size Affine Projection Sign Algorithm for System Identification in the Presence of Impulsive Noise
Previous Article in Journal
Analysis of Solutions to a Parabolic System with Absorption
Previous Article in Special Issue
Lagrangian Regularized Twin Extreme Learning Machine for Supervised and Semi-Supervised Classification
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Distributed Integrated Synthetic Adaptive Multi-Objective Reactive Power Optimization

1
Department of Mechanical and Electrical Engineering, Northeast Forestry University, Harbin 150040, China
2
State Grid HLJ Electric Power T&T Engineering Co., Ltd., Harbin 150040, China
*
Authors to whom correspondence should be addressed.
Symmetry 2022, 14(6), 1275; https://doi.org/10.3390/sym14061275
Submission received: 29 April 2022 / Revised: 6 June 2022 / Accepted: 17 June 2022 / Published: 20 June 2022
(This article belongs to the Special Issue Adaptive Filtering and Machine Learning)

Abstract

:
Reactive power is the core problem of voltage stability and economical operation in power systems. Aiming at the problem that multi-objective normalization reactive power optimization function is dependent on weight, an integrated synthesis of adaptive multi-objective particle swarm optimization (ISAMOPSO) is proposed to achieve weight adaptive. Through seven test functions and three algorithm comparison experiments, it is proved that the ISAMOPSO algorithm has stronger global search capability and better convergence. Considering the optimal access position and capacity of distributed generation (DG), the ISAMOPSO algorithm is used for three-objective reactive power optimization. Finally, the results indicate that the ISAMOPSO algorithm can not only provide a variety of optimization schemes to meet different needs, but also realize dynamic reactive power optimization, which further proves that the algorithm can provide effective technical support for solving reactive power optimization problems in practical engineering.

1. Introduction

With the proposal of the smart grid, it has become a developing trend that distributed generation (DG) gradually joins and becomes an essential part of the power distribution system [1,2]. Experts in the power industry generally agree that the adoption of DG technology is the main way to save money, reduce primary energy consumption, and improve the flexibility and stability of power systems [3,4,5,6].
However, the distribution network with the DG is a multi-power structure, and the capacity and location of the DG make the grid situation complicated when it is connected. The power flow and node voltage in the line will also be changed with the change of the power output of the DG [7,8], and even reverse power flow occurs. The distribution network is at risk of voltage overruns and voltage fluctuations [9]. Therefore, it is of great significance to study the control method of distributed generation in a distribution network to improve the reliability of the distribution network.
Since the multi-objective reactive power optimization problem of distribution network with DG is a very complex non-linear problem, domestic and foreign scholars have proposed various simplified mathematical models and optimization algorithms to solve it. Many heuristics have been applied to reactive power optimization problems, such as genetic algorithms (GA) [5,10], particle swarm algorithms (PSO) [6,11,12], and deep learning algorithms [13]. Meanwhile, many scholars have improved these algorithms, which have improved the speed and accuracy of the algorithms to solve reactive power optimization problems to different degrees [14,15]. Considering that the influence of distributed generation joining the distribution network is complex and multi-layered, some mathematical models of multi-objective optimization has been proposed. Taking Reference [16] and Reference [17] as examples, the optimization objectives are to reduce the system power consumption, improve the node voltage level and increase the voltage stability margin of the distribution network. After assigning a set of weights to several functions to be optimized and transforming them into a single objective function, the GA algorithm and PSO algorithm are used to solve the problem, respectively. The most common method in dealing with multi-objective problems is the simplification of multiple objectives, such as the weighting method, the ε-dominant method, and the affiliation method. The single solution is very dependent on the selection of weights, especially in the case of being unable to subjectively judge the relationship between various objectives. It is difficult to determine the weight of each function to be optimized. In contrast, multi-objective optimization can effectively solve the multi-objective optimization problem, multiple solutions can be obtained after the iteration of the optimization algorithm, and decision-makers can flexibly select schemes according to the actual situation [18,19,20]. In Reference [21], a static reactive power optimization model with minimum network loss and voltage deviation is established and solved by the multi-objective particle swarm optimization (MOPSO) algorithm. Based on Reference [21], Reference [22] also considered the maximization of static voltage stability margin, established a three-dimensional static reactive power optimization model and solved it by non-dominated sorting genetic algorithm II (NSGA-II). However, the power output of DG varies with time, so it is more practical to consider dynamic reactive power optimization. Reference [23] established a dynamic reactive power optimization model for the distribution networks and solved it using an improved particle swarm algorithm, which was accompanied by excessive iterations and premature convergence.
Aiming at the problem that single-objective reactive power optimization can only provide a single optimization scheme, an improved particle swarm optimization (PSO) algorithm is proposed to achieve a comprehensive and adaptive three-objective optimization of distribution networks with the DG. In addition, it can provide a variety of control strategies to achieve dynamic optimization.
The contributions are highlighted as follows:
  • A new optimization method (called ISAMOPSO) is proposed based on the MOPSO.
  • The proposed ISAMOPSO’s performance is evaluated by testing the selected benchmark functions.
  • The active power loss, voltage deviation, and reactive power compensation cost are taken as the optimization objectives, and the multi-objective reactive power optimization mathematical model of the distribution network of a distributed generation is established.
  • The new proposed ISAMOPSO approach is applied to solve the reactive power optimization in the distribution networks with DG.
The remaining paper is structured as follows: Section 2 presents the model description of the distribution networks with DG. Section 3 proposes the ISAMOPSO and tests its validation and performance. Section 4 presents two case simulations on the IEEE-33 node system with DG to verify the effectiveness of the proposed ISAMOPSO. Section 5 presents the conclusions.

2. Multi-Objective Reactive Power Optimization Model

2.1. Objective Function

The three most important objective functions in a power system are selected as the objective functions for reactive power optimization, which are active power loss, voltage deviation, and reactive power compensation cost.

2.1.1. Active Power Loss

The expression for active power loss is as follows:
min f P l o s s = k = 1 N b r a G k ( i , j ) [ V i 2 + V j 2 2 V i V j cos ( θ i θ j ) ]
where f P l o s s is active power loss, N b r a is the total branch number, G k ( i , j ) is the conductance of branch k connecting nodes i and j , voltage amplitude and voltage phase of the corresponding nodes V and θ , respectively.

2.1.2. Voltage Deviation

The expression for voltage deviation is as follows:
{ min f Δ V = i = 1 N L ( V i V i s p e c Δ V i max ) Δ V i max = V i max V i min V i s p e c = V i max + V i min 2
where N L is the number of load nodes, Δ V i max is the maximum allowable voltage deviation of node i , and V i s p e c is the expected voltage amplitude of node i .

2.1.3. Reactive Power Compensation Cost

The expression for reactive power compensation cost is as follows:
min f cos t = f cos t C + f cos t D G = s N C C C A P s | Q q s | + S D G ( r ( 1 + r ) n ( 1 + r ) n 1 C 1 + C 2 )
where s N C is the adjustable point of reactive power compensation equipment, Q q s and C C A P s are the actual input capacity and unit capacity investment of reactive power compensation equipment of node S , respectively. When C C A P s = 1 , the total investment of reactive power compensation equipment is to minimize the investment of reactive power compensation, S D G is the capacity of distributed power, r is the discount rate, n is the service life of distributed power, C 1 is the investment cost of distributed power, and C 2 is the operation cost of distributed power.

2.2. Constraint Condition

2.2.1. Equality Constraint

The mathematical expression of the power flow balance equation of active power and reactive power after DG access is shown below.
{ P G i + P D G i P L i = V i j = 1 N b V j ( G i j cos δ i j + B i j sin δ i j ) Q G i + Q D G i + Q C i Q L i = V i j = 1 N b V j ( G i j sin δ i j B i j cos δ i j )
where P G and Q G are the active power and reactive power injected into the generator, P L and Q L are the active power and reactive power consumed by the load, B is the admittance, G is the conductance, P G is the active output of DG, and Q D G and Q C are the reactive power capacity of DG and the reactive power capacity of reactive power compensation equipment.

2.2.2. Inequality Constraint

Where V G i min and V G i max are the minimum and maximum of the terminal voltage, and N G is the number of generator nodes. K j min and K j max are the lower and upper limits of the j adjustable transformer ratio, and N K is the number of adjustable transformers. C k min and C k max are the lower limit and upper limit of k compensation capacity of reactive power compensator, and N C is the number of nodes of reactive power compensation equipment. Q D G g min and Q D G g max are the lower and upper limits of the reactive power capacity of DG, and N D G is the number of nodes compensated by DG. Q G i min and Q G i max are the lower and upper limits of generator output reactive power, and N G is the number of generator nodes. V l min and V l max are the lower and upper limits of the voltage amplitude of the load node, and N B is the number of load nodes.
V G i min V G i V G i max , ( i = 1 , 2 , , N G )
K j min K j K j max , ( j = 1 , 2 , , N K )
C k min C k C k max , ( k = 1 , 2 , , N C )
Q D G g min Q D G g Q D G g max , ( g = 1 , 2 , , N D G )
Q G i min Q G i Q G i max , ( i = 1 , 2 , , N G )
V l min V l V l max , ( l = 1 , 2 , , N B )

3. ISAMOPSO Algorithm

To improve the global searching ability of PSO, the learning factor and inertia weight of PSO are improved adaptively.

3.1. The Basic Principle of the PSO Algorithm

The core operation of the PSO algorithm is the nearest individual away from the particle to avoid a collision and guides particles to move toward the center of the target and group. The core principle of PSO mainly includes the update of speed and position. The update process is expressed as:
v i k + 1 = ω v i k + c 1 r 1 ( x p i k x i k ) + c 2 r 2 ( x g i k x i k )
x i k + 1 = x i k + v i k + 1
where k is the number of iterations, i is the particle number, c 1 and c 2 are learning factors, and r 1 and r 2 are independent pseudorandom numbers, obeying uniform distribution on [ 0 , 1 ] . ω is the inertia weight. v i k is the search speed for a particle i at the kth iteration in the range [ V max , V max ] . x i k is the position of a particle i at the kth iteration, and its range is [ x max , x max ] . x p i k is the optimal position searched so far by a particle i in the kth iteration. x g i k is the optimal position in the kth iteration of the entire particle swarm.

3.2. Introducing Random Black Hole Strategy

The basic principle of introducing the black hole mechanism into PSO in [24] is to regard the particle population in the solution space as a planet in space, the fitness value of the particle as an attraction, and the real optimal solution as a black hole, and update the velocity and position of the particle according to the black hole mechanism. The principle is shown in Figure 1.
During the iteration process, for each dimension d of each particle, a black hole is randomly generated in the region with the dimension p g d corresponding to the global optimal solution as the center and r is the radius. A constant threshold p is set as the ability of black to capture the corresponding dimension x i d of the particle threshold p [ 0 , 1 ] . Set a random number l [ 0 , 1 ] for x i d . If the random number l p , x i d is captured by a black hole; if the random number l > p , after introducing the random black hole strategy, the particle update formula is:
x i d = { x i d + v i d , l > p p g d + 2 r ( r 3 0.5 ) , l p
where r 3 , l and p are random numbers with uniform distribution on [0,1]. l is the probability of a particle x i d appearing on its corresponding dimension, r is the radius of the black hole.

3.3. Adaptive Adjustment of Inertia Weight and Learning Factor

The ordinary MOPSO algorithm can easily fall into the local optimum. By changing the values of ω , c 1 and c 2 in the iterative process of the algorithm, it can effectively avoid the premature end of the iterative process without finding the optimal solution. The calculation method is:
ω = ω 0 + r 4 2 ( 1 ω 0 ) t T max
{ c 1 = 2 cos 2 [ π t 2 T max ] c 2 = 2 cos 2 [ π ( T max t ) 2 T max ]
where ω 0 takes the constant on [ 0 , 0.5 ] , r 4 takes the random number on [ 0 , 1 ] which obeys uniform distribution, and T max is the maximum iteration number. Additionally, t is the current iteration number.

3.4. The Diverse Selection of Leader Particles

The swarm induction strategy is introduced into PSO, and a disturbance is imposed on the population to increase the diversity of swarms. These disturbed particles as leading particles may lead the population to jump out of the local optimum to enhance its development ability. The individuals in the new sensing population are compared with the non-dominated individuals of the current generation, and the high-quality particles are selected for the next generation of particle swarm, and the formula is:
X t u r = X n e w + β V max s i g n ( 2 ( r a n d ( 1 , N ) 0.5 ) )
where X t u r is the particle position of the new sensing population, X n e w is the particle position of the next generation population, β is the disturbance amplitude obeying uniform distribution [ 0 , 1 ] , V max is the maximum flight velocity of the particle, and the size of the perceived group is N .
The particle with the largest fitness value was selected as the lead particle to guide the evolution of the next generation of the population. When a leading particle is selected, all particles in the non-dominant solution set are dynamically calculated according to the equations, which are expressed as:
{ f i t n e s s = 1 i = 1 M w i f i w i = λ i i = 1 M λ i , λ i ( 0 , 1 )
The number of objective functions is M , and the value of the ith objective function is f i , and the λ i obeys the uniformly distributed random number [ 0 , 1 ] .

3.5. Cyclic Elimination Strategy of Crowding Distance

Before sorting the crowding distance, the Rt population is sorted non-inferiorly and divided into different dominating levels. By calculating the crowding distance instead of one-time deletion, the crowding distance is gradually deleted from small to large using the idea of a dynamic cycle. The schematic diagram of the operation process is shown in Figure 2.
R t is divided into different classes according to the dominance relationship, where the F 1 class is the highest, and the class decreases in the following order. If the number of members of F 1 is less than N , then all the members of F 1 are selected for the population Pt+1, and the remaining members are selected from F 2 , F 3 ... until the size of N . The number of members in F 3 is greater than the number of remaining members in P t + 1 , so we need to perform crowding distance sorting at the same level of F 3 and select the number of members into P t + 1 . The crowded distance ranking is used to calculate the density of other solutions around a solution, which is calculated as shown in (18).
P i = m = 1 M P ( i + 1 ) · f m P ( i 1 ) · f m f max · m f min · m
Only one solution is deleted at each time, and the upper limit of deletion is determined immediately after deletion until the scale of P t + 1 is N . N d e l e t e is the number of eliminations, which is expressed as:
N d e l e t e = N 1 + N 2 + + N k N
where N k is the number of k-level non-inferior solutions, N is the size of the population P t + 1 , and the sum of the current number of k-level non-inferior solutions is just greater than N to achieve the elimination condition.

3.6. ISAMOPSO Algorithm Performance Verification and Analysis

The ISAMOPSO algorithm is used for the complex problem which is non-linear, multivariable, and under multi-constraint conditions. Therefore, the performance of the algorithm is evaluated by a comprehensive evaluation index that reflects the convergence and distribution uniformity of the algorithm, including the Inverted Generational Distance (IGD) and the Hypervolume (HV).
IGD is expressed as:
I G D = x * P * d ( x * , P ) n
where P is the real Pareto optimal solution set, P * is the non-dominated solution set solved by the algorithm, d ( x * , P ) represents the minimum Euclidean distance between the solution on the Pareto optimal surface and the individual in the non-dominated solution set solved by the algorithm, and n is the cardinal number of P * . The smaller the IGD value, the closer the non-dominated solution set is to the actual Pareto optimal solution set.
HV is expressed as:
H V = λ i = 1 | S | v i
where λ represents the Lebeg measure, v i represents the hypervolume composed of reference points and non-dominated individuals, and S represents the non-dominated solution set. The larger the HV value, the stronger the dominance of a single solution, and the closer it is to the Pareto boundary.
During the experiment, the performance of ISAMOPSO is compared with those of NSGA-II [25], MOPSO [26], and speed-constrained multi-objective particle swarm optimization (SMPSO) [27]; the swarm size is set to 100; the Maximum iteration is set to 10,000; and DTLZ1-DTLZ7 as test functions are run independently 30 times on the algorithm. The IGD performance index values and the mean values and standard deviations are listed in Table 1. It can be found that the mean value of IGD calculated by the ISAMOPSO algorithm is the smallest compared with the result of the other algorithm in the test functions DTLZ1, DTLZ2, DTLZ5, and DTLZ6, and the mean value of IGD calculated by the ISAMOPSO algorithm is slightly worse than that of SMPSO in the test function DTLZ3, and test functions, DTLZ4 and DTLZ7 are slightly worse than that of NSGAII. To more intuitively show the convergence and distribution of the optimal solution set obtained by various algorithms. The result with the best IGD value is selected to show the comparison between the approximate Pareto optimal frontier obtained by the four algorithms on all test functions and the real Pareto optimal front, as shown in Table 2. The results show that the non-dominated solution obtained by the ISAMOPSO algorithm is closer to the real Pareto optimal solution set than that obtained by other algorithms.
The HV performance index values and the mean values and standard deviations are listed in Table 3. The HV means calculated by the ISAMOPSO algorithm on the test function DTLZ1-DTLZ7 are greater than or approximately equal to the HV mean calculated by other algorithms, indicating that the dominant nature of the ISAMOPSO algorithm is stronger than that of its three algorithms, while the individual dominance of NSGA-II is stronger than the dominance of each individual obtained by the other two algorithms. Moreover, the Pareto optimal solution set obtained by the ISAMOPSO algorithm is closer to the real Pareto optimal front, indicating that the convergence and diversity of ISAMOPSO algorithm are better.

4. Applied ISAMOPSO for Reactive Power Optimization

To find the optimal capacity and location of the DG in the power distribution system and solve the multi-objective reactive power optimization problem of distribution networks with DG, this paper takes the IEEE-33 node distribution system as an application example [28]. There are 33 branches in the unimproved IEEE-33 node system. The active load is 3715 kW, the reactive load is 2300 kvar, the system reference voltage is 102.66 kV, and the reference capacity is 10 MVA. The calculation accuracy in the test is 10−4. The IEEE-33 node distribution system is shown in Figure 3.

4.1. Influence of Distributed Generation on Distribution Network

After the DGs are connected, the whole distribution network becomes a complex structure of multiple power sources, which is not only used to transmit power but also can generate power. The random power emitted by DG will cause the change in power flow, and some DGs with larger capacity will even cause the voltage instability at each node of the grid and the increase in network loss, which is closely related to the connection location and capacity of DGs. Therefore, it is necessary to consider the location and capacity of DGs connected to the distribution network.

4.1.1. Influence of Distributed Generation Capacity on Distribution Network

According to the structural characteristics of the IEEE-33 node distribution system, 14 groups of independent experiments were conducted on 14 representative nodes. The DG was incorporated into the 3rd, 6th, 9th, 12th, 15th, 18th, 19th, 21st, 23rd, 25th, 26th, 28th, 31st, and 33rd nodes to test the influence of the voltage and network loss of the distribution network system. The total load of the IEEE-33 node system is 3715 kW + j2300 kvar. The capacities of DG are 20%, 40%, 50%, 70%, 90% and 100% of the total load. The results of the experiments are shown below.
The per-unit values of the network loss when the different capacities of DG are incorporated into different 14 nodes are shown in Table 4.
Through the analysis of Figure 4a,g,i, it can be concluded when the DG is close to the power supply side, it has no obvious effect on the voltage rise no matter how much the capacity of the DG. Meanwhile, the voltage increases significantly and even exceeds the limit value when the DG is incorporated into the 15th, 18th, 31st, and 33rd nodes, which are far away from the power supply side. The higher the DG capacity, the higher the node voltage of the system, as can be seen in Figure 4e,f,l,m.
From Table 4, the network loss decreases with the increase in the capacity of the DG. The capacity of the DG should be between 20% and 40% of the total load. If the DG is incorporated near the power supply side, the capacity of the DG should be selected at 20% of the total load to reduce system power loss. Otherwise, the capacity should be 40% of the total load to increase the power flow transmission capacity and minimize the network loss so that the voltage distribution is more stable.

4.1.2. Influence of Distributed Generation Position on Distribution Network

Using the control variables method, when the capacity of the DG is certain, the location of the DG in the IEEE-33 node distribution system should be changed to study the influence of the node voltage distribution. To improve the experimental reliability and summarize the general rule, six groups of independent experiments with different capacities were set up. The capacities of DG were 20%, 40%, 50%, 70%, 90% and 100% of the total load. The voltage amplitudes of 14 nodes in each group of experiments are shown in the following (Figure 5):
The results of experiments show that no matter how much capacity of DG is incorporated into the 3rd, 6th, 19th, and 23rd nodes of the system, the voltage amplitude of 33 nodes is the same as the value there is no DG in the system. When the DG is incorporated into the 15th and 18th nodes, the voltage amplitude of 33 nodes is higher the higher the capacity of DG, but when the DG is incorporated into the 28th and 31st nodes, the voltage amplitude of 33 nodes is nearly no change. When the DG is incorporated into the 21st and 35th nodes, only the voltage of nearby nodes is affected. The closer the access location of DG supply is to the end of the line and the load center, the more obvious its effect on voltage lifting, so it should be chosen to merge into the middle and rear sections of the mainline. From the above experiments, it is most appropriate to select the 12th, 15th, and 18th nodes for the integration point, and the capacity of DG should be between 20% and 40% of the system load.

4.2. Multi-Objective Reactive Power Optimization Solution Process

The main steps of the ISAMOPSO for distribution grid reactive power optimization with DG are as follows:
  • Step 1: Read the system data, set the population size N and the maximum iteration number M t and the related parameters of the improvement strategy, and initialize the population P , personal best P t and global best P g .
  • Step 2: Update the personal best and leader for each particle. Carry out power flow computation and obtain the objective function values of each particle. The non-dominated solutions of the population P are confirmed according to the non-dominated ranking strategy, and the positions of the next generation of particles and the individual optimal positions are updated according to Equation (13).
  • Step 3: Generate P s off-spring population Q according to Equations (11) and (13) based on the current velocity and position of each particle in P , and then combine P and Q together and obtain population R whose size is 2 N .
  • Step 4: Identify the non-dominated solutions from R , and store them in N D .
  • Step 5: Generate the particle population P n e w for the next iteration and the evaluation criterion Formula (17) for the leader particles.
  • Step 6: Based on the state of velocity in the population, determine whether it constitutes a perturbation condition, and adjust the population again.
  • Step 7: Determine whether the maximum number of iterations M is reached; if so, proceed to the next step; otherwise, go to Step 3 for the next iteration.
  • Step 8: Output the non-dominated solutions as the final Pareto-optimal solutions.

4.3. Analysis and Discussion Results

According to the conclusions on the influence of DG capacity and location in the IEEE-33 node distribution system, the DG is, respectively, incorporated into the 15th node and 31st node. It is assumed that each DG can output active power or compensate for the reactive power of the distribution network.
In the example of the IEEE-33 node reactive power compensation distribution system [28], the active load is 3715 kW, the reactive load is 2300 kvar, the system reference voltage is 102.66 kV, and the reference capacity is 10 MVA. Based on the discussion of the capacity of DG in Section 4.1, the active output of the DG is all 1 MW, and the reactive power output range is −0.1~0.5 Mvar. Based on the load of the 33-node reactive power compensation distribution system, the shunt capacitor bank with a compensation capacity of 0.15 Mvar × 4 is incorporated into the 6th node, and the shunt capacitor bank with a compensation capacity of 0.15 Mvar × 7 is incorporated into the 24th node [29]. The on-load voltage regulating transformer with nine transformer taps is connected between 0 and 1 node which ratio range is 0.9–1.1. and the transformer tap adjustment step is 0.025. The range of node voltage variation is 0.85–1.15 per unit. The active network loss of the system before optimization is 0.1885 MW, and the node voltage deviation is 1.0435. The IEEE -33 node reactive power compensation distribution system is shown in Figure 6.

4.3.1. Static Reactive Power Optimization

During the process of static reactive power optimization process, the objective functions are the minimum distribution network loss, node voltage deviation, and reactive power compensation cost. The reactive power output of the DG, the capacity of the shunt capacitor, and the transformer tap of the on-load regulator transformer were selected as the control variables.
The optimal solution set in the objective space is shown in Figure 7. There is a certain non-linear relationship among the three objectives, which are in conflict with each other, and no optimization result can make each objective function reach the optimal at the same time.
The optimal solution set in the objective space is shown in Figure 7. The five minimum Euclidean distance compromise solutions were selected from the optimal solution set as 5 optimization schemes (W1~W5), and the control variables and objective function values were calculated, which are shown in Table 5. It was observed that, while the active power loss ( f P l o s s ) and voltage deviation ( f Δ V ) are all smallest in optimization scheme W3, but the smallest reactive power compensation ( f cos t ) is in W4 and the minimum Euclidean distance is W1. Because there is a certain non-linear relationship among the three objectives, which are in conflict with each other, and no optimization result can make each objective function reach the optimal at the same time. Therefore, it is necessary to provide multiple solutions to meet different requirements. The results of the minimum Euclidean distance of the ISAMOPSO algorithm were compared with the minimum Euclidean distance of the NSGA-II algorithm, as shown in Table 6, the network loss decreased by 11.4%, the voltage deviation decreased by 18.7%, and the compensation cost decreased by 6.75%, which proves that the ISAMOPSO algorithm is more applicable to the reactive power optimization problem than the NSGA-II algorithm.
Using the five static reactive power optimization schemes, the voltage amplitudes of the 33 nodes are different and compared with the value before static reactive power optimization, which is shown in Figure 8.
The results show that the voltage amplitude of each node optimized by the ISAMOPSO algorithm is within its upper and lower limits. Moreover, the voltage amplitude difference between nodes is reduced, and the voltage distribution after optimization is smoother. This shows that the distribution network is more stable. As known from the calculation, the active power network loss maximum decreased by 77.82%, and the maximum voltage deviation was reduced by 86% compared with the value before optimization. The effectiveness of the ISAMOPSO algorithm for multi-objective reactive power optimization is proved.

4.3.2. Dynamic Reactive Power Optimization

Affected by natural environmental factors such as wind speed, light intensity, and temperature, the output of DG will change over time, and the load of the actual distribution network is also changing. In this section, according to the natural environmental factors of a typical day in a certain region, the photovoltaic and wind power output at 24 different times of the day is predicted. This is incorporated into the IEEE-33 node distribution system as the access capacity of DG. The natural environmental factors for a typical day in a region are shown in Table 7. Based on the data given in Table 7, the output of PV and wind power was predicted, and the predicted results are shown in Figure 9.
During the process of the dynamic reactive power optimization, the objective functions are the minimum distribution network loss, node voltage deviation, and maximum total consumption ratio of DG.
Considering load time-varying over 24 h a day, the access capacity of DG, the capacity of the shunt capacitor, and the transformer tap of the on-load regulator transformer were selected as the control variables. The ISAMOPSO was used for dynamic reactive power optimization. The optimal solution set in the three objective spaces of dynamic reactive power optimization is shown in Figure 10.
The three objective functions have conflicting distributions in the objective space with non-linear relationships and present a spatial curve in the three-dimensional space. The uniform distribution and high diversity of the points on the curve further illustrate the effectiveness of the ISAMOPSO algorithm in solving the multi-objective reactive power optimization problem.
In this paper, a compromise solution with the minimum Euclidean distance from the ideal point (0, 0, 0) was selected from the optimal solution set as dynamic reactive power optimization schemes.
The power consumed by photovoltaic power generation and wind power generation at each time after grid connection is very close to their respective predicted power outputs. The specific situation is shown in Figure 11 and Figure 12.
The total daily consumption ratio of photovoltaic power generation is 85.61%, and that of wind power generation is 82.75%.
The total consumption and power output of DG can be seen in Figure 13. The total consumption ratio of the system after optimization is very high at 0.8383, which proves that the algorithm improves the ability of the distribution network to absorb DG.
The whole-day average network loss decreased by 21.08% compared with before optimization, which can be seen in Figure 14.
The whole-day average voltage deviation decreased by 4.18% compared with before optimization, which can be seen in Figure 15.
The comparison of the voltage distribution of different nodes at 24 different times in Figure 16 shows that the voltage at each node after optimization is smoother than before optimization and the system operation is more stable. This further proves the effectiveness of the ISAMOPSO algorithm to solve the multi-objective reactive power optimization problem and the rationality of the mathematical optimization model established.

5. Conclusions

To solve the reactive power compensation problem of distribution networks with DG for economical operation, reduction in network loss, and reduction in voltage deviation, this study investigates the optimal size and siting of DG, establishes a multi-objective optimization model for optimal regulation in the distribution system, and uses the ISAMOPSO algorithm as an optimization tool. The algorithm is compared with other algorithms and successfully applied to the IEEE-33 node system. In the static reactive power optimization experiment, the network loss decreases by 11.4%, the voltage deviation decreases by 18.7%, and the compensation cost decreases by 6.75% in the solution set obtained by the ISAMOPSO algorithm compared to the solution set obtained by the NSGA-II algorithm with the minimum Euclidean distance. In the dynamic reactive power optimization experiment, the multi-objective reactive power optimization strategy proposed in this paper takes into account the security and stability of the distribution network operation and the absorption capacity of DG. After optimization by the ISAMOPSO algorithm, the DG consumption rate reached 0.8383, the network loss for the whole day decreased by 21.08% compared with that before optimization and the voltage deviation for the whole day reduced by 4.18% compared with that before optimization. Thus, the ISAMOPSO algorithm proposed in this paper has good computational capability for solving non-convex, multi-objective, multi-constrained high-dimensional non-linear reactive power optimization problems in power systems, and can provide a set of Pareto-optimal solution sets as optimization schemes, each of which can optimize the operational performance of the system and provide a flexible reactive power optimization scheme for decision-makers.

Author Contributions

Conceptualization, J.S., H.Z., Q.Y. and Q.Z.; methodology, J.S.; software, Y.Z.; validation, J.S., C.L. and Q.M.; formal analysis, Q.M., Y.F. and Q.Y.; investigation, C.L.; resources, Q.Y.; data curation, Q.M., C.L. and Q.H.; writing—original draft preparation, J.S.; writing—review and editing, J.S.; visualization, Q.M., C.L.; supervision, J.S.; project administration, J.S.; funding acquisition, J.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by “The Fundamental Research Funds for the Central Universities”, grant number: 2572017CB13, and funded by “Heilongjiang Provincial Natural Science Foundation of China”, grant number: YQ2020C018.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Diahovchenko, I.; Viacheslav, Z. Optimal composition of alternative energy sources to minimize power losses and maximize profits in a distribution power network. In Proceedings of the 2020 IEEE 7th International Conference on Energy Smart Systems (ESS), Kyiv, Ukraine, 12–14 May 2020; pp. 247–252. [Google Scholar] [CrossRef]
  2. Ma, W.; Wang, W.; Chen, Z.; Wu, X.; Hu, R.; Tang, F.; Zhang, W. Voltage regulation methods for active distribution networks considering the reactive power optimization of substations. Appl. Energy 2021, 284, 116347. [Google Scholar] [CrossRef]
  3. Wang, Z.; Liu, H.; Li, J. Reactive Power Planning in Distribution Network Considering the Consumption Capacity of Distributed Generation. In Proceedings of the 2020 5th Asia Conference on Power and Electrical Engineering (ACPEE), Chengdu, China, 4–7 June 2020; pp. 1122–1128. [Google Scholar] [CrossRef]
  4. Side, D.O.; Picioroaga, I.I.; Tudose, A.M.; Bulac, C.; Tristiu, I. Multi-Objective Particle Swarm optimization Applied on the Optimal Reactive Power Dispatch in Electrical Distribution Systems. In Proceedings of the 2020 International Conference and Exposition on Electrical And Power Engineering (EPE), Iasi, Romania, 22–23 October 2020; pp. 413–418. [Google Scholar] [CrossRef]
  5. Cheng, X.; Chen, X.; Fan, Y.; Liu, Y.; Li, S.; Zhao, Q. Research on Reactive Power Optimization Control Strategy of Distribution Network with Photovoltaic Generation. In Proceedings of the 2019 Chinese Automation Congress (CAC), Hangzhou, China, 22–24 November 2019; pp. 4026–4030. [Google Scholar] [CrossRef]
  6. Bi, H. The significance and current situation of reactive power optimization. AIP Conf. Proc. 2019, 2122, 020025. [Google Scholar] [CrossRef]
  7. Chen, S.; Li, P.; Ji, H.; Yu, H.; Yan, J.; Wu, J.; Wang, C. Operational flexibility of active distribution networks with the potential from data centers. Appl. Energy 2021, 293, 116935. [Google Scholar] [CrossRef]
  8. Liu, J.; Lu, Y.; Lv, G.; Du, J. Research on Access Point and Capacity Optimization of Distributed Generation in Distribution Power Grid. In Proceedings of the 2021 China International Conference on Electricity Distribution (CICED), Shanghai, China, 7–9 April 2021; pp. 718–721. [Google Scholar] [CrossRef]
  9. Tonkoski, R.; Lopes, L.A.C.; El-Fouly, T.H.M. Coordinated Active Power Curtailment of Grid Connected PV Inverters for Overvoltage Prevention. IEEE Trans. Sustain. Energy 2021, 2, 139–147. [Google Scholar] [CrossRef]
  10. Siyuan, C.; Suwen, Z. Research on Reactive Power Optimization Algorithm of Power System Based on Improved Genetic Algorithm. In Proceedings of the 2017 International Conference on Industrial Informatics—Computing Technology, Intelligent Technology, Industrial Information Integration (ICIICII), Wuhan, China, 2–3 December 2017; pp. 21–24. [Google Scholar] [CrossRef]
  11. Liu, B.; Shi, L.; Yao, Z. Multi-objective optimal reactive power dispatch for distribution network considering pv generation uncertainty. In Proceedings of the 10th Renewable Power Generation Conference (RPG 2021), Online, 14–15 October 2021; pp. 503–509. [Google Scholar] [CrossRef]
  12. Soldevilla, F.R.C.; Huerta, F.A.C. Minimization of Losses in Power Systems by Reactive Power Dispatch using Particle Swarm Optimization. In Proceedings of the 2019 54th International Universities Power Engineering Conference (UPEC), Bucharest, Romania, 3–6 September 2019; pp. 1–5. [Google Scholar] [CrossRef]
  13. Srikanth, V.; Natarajan, V.; Jegajothi, B.; Arumugam, S.S.L.D.; Nageswari, D. Fruit fly Optimization with Deep Learning Based Reactive Power Optimization Model for Distributed Systems. In Proceedings of the 2022 International Conference on Electronics and Renewable Systems (ICEARS), Tuticorin, India, 16–18 March 2022; pp. 319–324. [Google Scholar] [CrossRef]
  14. Ali, M.H.; Kamel, S.; Hassan, M.H.; Tostado-Véliz, M.; Zawbaa, H.M. An improved wild horse optimization algorithm for reliability based optimal DG planning of radial distribution networks. Energy Rep. 2022, 8, 582–604. [Google Scholar] [CrossRef]
  15. Lian, L. Reactive power optimization based on adaptive multi-objective optimization artificial immune algorithm. Ain Shams Eng. J. 2022, 13, 101677. [Google Scholar] [CrossRef]
  16. Yuan, Z.; Dou, X.; Wang, L.; Dong, Y. Reactive Power Optimization of Distribution Network Based on BP Neural Network and Intelligent Scene Matching. In Proceedings of the 2021 IEEE Sustainable Power and Energy Conference (iSPEC), Nanjing, China, 23–25 December 2021; pp. 3048–3054. [Google Scholar] [CrossRef]
  17. Mendes, A.; Boland, N.; Guiney, P.; Riveros, C. (N-1) contingency planning in radial distribution networks using genetic algorithms. In Proceedings of the 2010 IEEE/PES Transmission and Distribution Conference and Exposition: Latin America (T&D-LA), Sao Paulo, Brazil, 8–10 November 2010; pp. 290–297. [Google Scholar] [CrossRef]
  18. Hui, Q.; Teng, Y.; Zuo, H.; Chen, Z. Reactive power multi-objective optimization for multi-terminal AC/DC interconnected power systems under wind power fluctuation. CSEE J. Power Energy Syst. 2020, 6, 630–637. [Google Scholar] [CrossRef]
  19. Bo, S.; Yujie, C.; Xiang, Z.; Yelin, W.; Ning, L. Multi-objective reactive power optimization in power system based on improved PSO-OSA algorithm. In Proceedings of the 2019 14th IEEE Conference on Industrial Electronics and Applications (ICIEA), Xi’an, China, 19–21 June 2019; pp. 1514–1519. [Google Scholar] [CrossRef]
  20. Zhang, C.; Chen, M.; Luo, C. A multi-objective optimization method for power system reactive power dispatch. In Proceedings of the 2010 8th World Congress on Intelligent Control and Automation, Jinan, China, 7–9 July 2010; pp. 6–10. [Google Scholar] [CrossRef]
  21. Dhifaoui, C.; Guesmi, T.; Welhazi, Y.; Abdallah, H.H. Application of multi-objective PSO algorithm for static reactive power dispatch. In Proceedings of the 2014 5th International Renewable Energy Congress (IREC), Hammamet, Tunisia, 25–27 March 2014; pp. 1–6. [Google Scholar] [CrossRef]
  22. Zhang, Z.; Wang, X. Optimal Reactive Optimization of Distribution Network with Wind Turbines Based on Improved NSGA-II. In Proceedings of the 2018 China International Conference on Electricity Distribution (CICED), Tianjin, China, 17–19 September 2018; pp. 2056–2060. [Google Scholar] [CrossRef]
  23. Li, S.; He, Y.; Zhou, Y.; Hu, Y. A New Reactive Power optimization Method For Distribution Network with High Participation of Renewable Distributed Generation. In Proceedings of the 2020 IEEE 4th Conference on Energy Internet and Energy System Integration (EI2), Wuhan, China, 30 October–1 November 2020; pp. 1298–1302. [Google Scholar] [CrossRef]
  24. Zhang, J.; Liu, K.; Tan, Y.; He, X. Random black hole particle swarm optimization and its application. In Proceedings of the 2008 International Conference on Neural Networks and Signal Processing, Nanjing, China, 7–11 June 2008; pp. 359–365. [Google Scholar] [CrossRef]
  25. Deb, K.; Pratap, A.; Agarwal, S.; Meyarivan, T. A fast and elitist multiobjective genetic algorithm: NSGA-II. IEEE Trans. Evol. Comput. 2002, 6, 182–197. [Google Scholar] [CrossRef] [Green Version]
  26. Coello, C.A.C.; Pulido, G.T.; Lechuga, M.S. Handling multiple objectives with particle swarm optimization. IEEE Trans. Evol. Comput. 2004, 8, 256–279. [Google Scholar] [CrossRef]
  27. Nebro, A.J.; Durillo, J.J.; Garcia-Nieto, J.; Coello, C.A.C.; Luna, F.; Alba, E. SMPSO: A new PSO-based metaheuristic for multi-objective optimization. In Proceedings of the 2009 IEEE Symposium on Computational Intelligence in Multi-Criteria Decision-Making(MCDM), Nashville, TN, USA, 30 March–2 April 2009; pp. 66–73. [Google Scholar] [CrossRef]
  28. Baran, M.E.; Wu, F.F. Network reconfiguration in distribution systems for loss reduction and load balancing. IEEE Trans. Power Deliv. 1989, 4, 1401–1407. [Google Scholar] [CrossRef]
  29. Cheng, S.; Chen, M.-Y. Multi-objective reactive power optimization strategy for distribution system with penetration of distributed generation. Int. J. Electr. Power Energy Syst. 2014, 62, 221–228. [Google Scholar] [CrossRef]
Figure 1. Particle position update process after introducing random black hole strategy.
Figure 1. Particle position update process after introducing random black hole strategy.
Symmetry 14 01275 g001
Figure 2. Operation diagram of crowding distance ranking strategy.
Figure 2. Operation diagram of crowding distance ranking strategy.
Symmetry 14 01275 g002
Figure 3. The IEEE-33 node distribution system.
Figure 3. The IEEE-33 node distribution system.
Symmetry 14 01275 g003
Figure 4. (a) The voltage amplitude of 33 nodes when the different capacities of DG are incorporated into the 3rd node; (b) The voltage amplitude of 33 nodes when the different capacities of DG are incorporated into the 6th node; (c) The voltage amplitude of 33 nodes when the different capacities of DG are incorporated into the 9th node; (d) The voltage amplitude of 33 nodes when the different capacities of DG are incorporated into the 12th node; (e) The voltage amplitude of 33 nodes when the different capacities of DG are incorporated into the 15th node; (f) The voltage amplitude of 33 nodes when the different capacities of DG are incorporated into the 18th node; (g) The voltage amplitude of 33 nodes when the different capacities of DG are incorporated into the 19tht node; (h) The voltage amplitude of 33 nodes when the different capacities of DG are incorporated into the 21st node; (i) The voltage amplitude of 33 nodes when the different capacities of DG are incorporated into the 23rd node; (j) The voltage amplitude of 33 nodes when the different capacities of DG are incorporated into the 25th node; (k) The voltage amplitude of 33 nodes when the different capacities of DG are incorporated into the 28th node; (l) The voltage amplitude of 33 nodes when the different capacities of DG are incorporated into the 31st node; (m) The voltage amplitude of 33 nodes when the different capacities of DG are incorporated into the 33rd node.
Figure 4. (a) The voltage amplitude of 33 nodes when the different capacities of DG are incorporated into the 3rd node; (b) The voltage amplitude of 33 nodes when the different capacities of DG are incorporated into the 6th node; (c) The voltage amplitude of 33 nodes when the different capacities of DG are incorporated into the 9th node; (d) The voltage amplitude of 33 nodes when the different capacities of DG are incorporated into the 12th node; (e) The voltage amplitude of 33 nodes when the different capacities of DG are incorporated into the 15th node; (f) The voltage amplitude of 33 nodes when the different capacities of DG are incorporated into the 18th node; (g) The voltage amplitude of 33 nodes when the different capacities of DG are incorporated into the 19tht node; (h) The voltage amplitude of 33 nodes when the different capacities of DG are incorporated into the 21st node; (i) The voltage amplitude of 33 nodes when the different capacities of DG are incorporated into the 23rd node; (j) The voltage amplitude of 33 nodes when the different capacities of DG are incorporated into the 25th node; (k) The voltage amplitude of 33 nodes when the different capacities of DG are incorporated into the 28th node; (l) The voltage amplitude of 33 nodes when the different capacities of DG are incorporated into the 31st node; (m) The voltage amplitude of 33 nodes when the different capacities of DG are incorporated into the 33rd node.
Symmetry 14 01275 g004aSymmetry 14 01275 g004bSymmetry 14 01275 g004c
Figure 5. (a) The voltage amplitude of 33 nodes when the capacity of DG is 20% of the total load, and it is incorporated into different 14 nodes; (b) The voltage amplitude of 33 nodes when the capacity of DG is 40% of the total load, and it is incorporated into different; (c) The voltage amplitude of 33 nodes when the capacity of DG is 50% of the total load, and it is incorporated into different 14 nodes; (d) The voltage amplitude of 33 nodes when the capacity of DG is 70% of the total load, and it is incorporated into different 14 nodes; (e) The voltage amplitude of 33 nodes when the capacity of DG is 90% of the total load, and it is incorporated into different 14 nodes; (f) The voltage amplitude of 33 nodes when the capacity of DG is 100% of the total load, and it is incorporated into different 14 nodes.
Figure 5. (a) The voltage amplitude of 33 nodes when the capacity of DG is 20% of the total load, and it is incorporated into different 14 nodes; (b) The voltage amplitude of 33 nodes when the capacity of DG is 40% of the total load, and it is incorporated into different; (c) The voltage amplitude of 33 nodes when the capacity of DG is 50% of the total load, and it is incorporated into different 14 nodes; (d) The voltage amplitude of 33 nodes when the capacity of DG is 70% of the total load, and it is incorporated into different 14 nodes; (e) The voltage amplitude of 33 nodes when the capacity of DG is 90% of the total load, and it is incorporated into different 14 nodes; (f) The voltage amplitude of 33 nodes when the capacity of DG is 100% of the total load, and it is incorporated into different 14 nodes.
Symmetry 14 01275 g005aSymmetry 14 01275 g005b
Figure 6. The IEEE-33 node reactive power compensation distribution system.
Figure 6. The IEEE-33 node reactive power compensation distribution system.
Symmetry 14 01275 g006
Figure 7. The optimal solution set in the three objective spaces of the static reactive power optimization.
Figure 7. The optimal solution set in the three objective spaces of the static reactive power optimization.
Symmetry 14 01275 g007
Figure 8. The voltage amplitude of 33 nodes using different static reactive power optimization schemes.
Figure 8. The voltage amplitude of 33 nodes using different static reactive power optimization schemes.
Symmetry 14 01275 g008
Figure 9. Prediction power output of DG.
Figure 9. Prediction power output of DG.
Symmetry 14 01275 g009
Figure 10. The optimal solution set in the three objective spaces of dynamic reactive power optimization.
Figure 10. The optimal solution set in the three objective spaces of dynamic reactive power optimization.
Symmetry 14 01275 g010
Figure 11. Photovoltaic power output and actual consumption.
Figure 11. Photovoltaic power output and actual consumption.
Symmetry 14 01275 g011
Figure 12. Wind power output and actual consumption.
Figure 12. Wind power output and actual consumption.
Symmetry 14 01275 g012
Figure 13. Comparison of total power output and total actual consumption.
Figure 13. Comparison of total power output and total actual consumption.
Symmetry 14 01275 g013
Figure 14. Comparison of network loss before and after optimization.
Figure 14. Comparison of network loss before and after optimization.
Symmetry 14 01275 g014
Figure 15. Comparison of voltage deviation before and after optimization.
Figure 15. Comparison of voltage deviation before and after optimization.
Symmetry 14 01275 g015
Figure 16. (a) Voltage distribution map before optimization; (b) voltage distribution map after optimization.
Figure 16. (a) Voltage distribution map before optimization; (b) voltage distribution map after optimization.
Symmetry 14 01275 g016
Table 1. The result of the mean and standard deviation of IGD using the test function in different algorithms.
Table 1. The result of the mean and standard deviation of IGD using the test function in different algorithms.
Test FunctionIndexAlgorithm
NSGAIIMOPSOSMPSOISAMOPSO
DTLZ1Mean2.6664 × 10−23.1767 × 10−23.0328 × 10−22.6595 × 10−2
Std1.28 × 10−31.75 × 10−21.20 × 10−21.09 × 10−3
DTLZ2Mean6.8392 × 10−27.7528 × 10−26.6816 × 10−26.3912 × 10−2
Std3.38 × 10−38.42 × 10−32.30 × 10−31.79 × 10−3
DTLZ3Mean6.8458 × 10−21.2056 × 10−16.8460 × 10−29.9426 × 10−2
Std2.19 × 10−37.13 × 10−22.63 × 10−34.33 × 10−2
DTLZ4Mean6.6723 × 10−21.7982 × 10−11.1460 × 10−11.5361 × 10−1
Std3.02 × 10−31.15 × 10−11.06 × 10−11.64 × 10−2
DTLZ5Mean5.7555 × 10−37.7760 × 10−35.3310 × 10−34.8635 × 10−3
Std3.61 × 10−46.66 × 10−44.20 × 10−41.35 × 10−4
DTLZ6Mean5.7767 × 10−38.5702 × 10−35.3689 × 10−34.6755 × 10−3
Std2.82 × 10−49.20 × 10−42.09 × 10−49.00 × 10−5
DTLZ7Mean7.9676 × 10−21.0997 × 10−11.0327 × 10−18.2313 × 10−2
Std6.02 × 10−39.98 × 10−23.81 × 10−26.18 × 10−3
Table 2. Pareto front comparison results of the test functions in different algorithms.
Table 2. Pareto front comparison results of the test functions in different algorithms.
FunctionNSGAⅡMOPSOSMPSOISAMOPSO
DTLZ1 Symmetry 14 01275 i001 Symmetry 14 01275 i002 Symmetry 14 01275 i003 Symmetry 14 01275 i004
DTLZ2 Symmetry 14 01275 i005 Symmetry 14 01275 i006 Symmetry 14 01275 i007 Symmetry 14 01275 i008
DTLZ3 Symmetry 14 01275 i009 Symmetry 14 01275 i010 Symmetry 14 01275 i011 Symmetry 14 01275 i012
DTLZ4 Symmetry 14 01275 i013 Symmetry 14 01275 i014 Symmetry 14 01275 i015 Symmetry 14 01275 i016
DTLZ5 Symmetry 14 01275 i017 Symmetry 14 01275 i018 Symmetry 14 01275 i019 Symmetry 14 01275 i020
DTLZ6 Symmetry 14 01275 i021 Symmetry 14 01275 i022 Symmetry 14 01275 i023 Symmetry 14 01275 i024
DTLZ7 Symmetry 14 01275 i025 Symmetry 14 01275 i026 Symmetry 14 01275 i027 Symmetry 14 01275 i028
Table 3. The result of the mean and standard deviation of HV using the test function in different algorithms.
Table 3. The result of the mean and standard deviation of HV using the test function in different algorithms.
Test FunctionIndexAlgorithm
NSGAIIMOPSOSMPSOISAMOPSO
DTLZ1Mean8.2429 × 10−17.9694 × 10−18.1700 × 10−18.1891 × 10−1
Std3.74 × 1034.61 × 10−25.11 × 10−33.13 × 10−2
DTLZ2Mean5.3293 × 10−15.0980 × 10−15.3286 × 10−15.3832 × 10−1
Std5.53 × 10−31.29 × 10−23.90 × 10−33.56 × 10−3
DTLZ3Mean5.0795 × 10−14.6962 × 10−15.2880 × 10−15.3488 × 10−1
Std5.30 × 10−26.38 × 10−26.18 × 10−34.64 × 10−3
DTLZ4Mean5.4192 × 10−14.9007 × 10−15.2322 × 10−15.2573 × 10−1
Std3.29 × 10−34.11 × 10−23.85 × 10−26.46 × 10−3
DTLZ5Mean1.9930 × 10−11.9691 × 10−11.9944 × 10−11.9965 × 10−1
Std1.70 × 10−42.61 × 10−32.09 × 10−41.25 × 10−4
DTLZ6Mean1.9958 × 10−11.9781 × 10−11.9962 × 10−12.0015 × 10−1
Std1.30 × 10−44.58 × 10−41.74 × 10−43.57 × 10−5
DTLZ7Mean2.7128 × 10−12.6964 × 10−12.6817 × 10−12.7167 × 10−1
Std1.68 × 10−31.17 × 10−22.67 × 10−32.57 × 10−3
Table 4. The per-unit values of the network loss when the different capacities of DG are incorporated into different 14 nodes.
Table 4. The per-unit values of the network loss when the different capacities of DG are incorporated into different 14 nodes.
Node020%40%50%70%90%100%
30.02930.02280.0220.02160.02110.02080.0207
60.02390.01590.0110.00970.00890.01040.012
90.02390.01410.01070.0110.01480.02250.0276
120.02390.01340.01180.01350.0210.03310.0405
150.02390.01340.0150.01890.03150.04910.0595
180.02390.01470.01980.02580.0430.06550.0785
190.02390.02320.02290.02290.02310.02370.0241
210.02390.02370.02640.02860.0350.04340.0484
230.02390.02240.02170.02170.02230.02360.0246
250.02390.02160.02260.02430.02960.03730.042
260.02390.01550.01070.00940.0090.01120.0132
280.02390.01350.0090.00870.01110.01730.0216
310.02390.01180.00950.0110.01850.03060.0382
330.02390.01210.01080.01310.0220.03590.0444
Table 5. Five static reactive power optimization schemes.
Table 5. Five static reactive power optimization schemes.
ResultDG1
(Mvar)
DG2
(Mvar)
C1
(Mvar)
C2
(Mvar)
K
(pu)
f P l o s s
(MW)
f Δ V f c o s t
(million)
W10.03350.29650.150.150.9750.06450.15810.0594
W20.13250.32140.150.150.9750.06630.15950.0569
W30.16180.29330.30.450.9750.05120.12870.0701
W40.31870.33960.450.7510.05470.27410.0421
W50.35710.41630.450.310.05820.28030.0486
Table 6. Comparison of the optimization results of the algorithms.
Table 6. Comparison of the optimization results of the algorithms.
Comparison ObjectBefore OptimizationNSGA-IIISAMOPSO
f P l o s s (MW)0.18850.07280.0645
f Δ V 1.04350.19450.1581
f cos t (million)-0.06370.0594
Table 7. Typical daily environmental factors and power output of DG.
Table 7. Typical daily environmental factors and power output of DG.
TimeLight Intensity (W/m2)Temperature (℃)Wind Speed (m/s)PV Output (kW)WT Output (kW) TimeLight Intensity (W/m2)Temperature (℃)Wind Speed (m/s)PV Output (kW)WT Output (kW)
0017.337.150219.412494.4427.465.90287.9134.3
1017.298.320313.613619.4430.106.72315.1188.1
2017.107.530248.414605.5629.265.82403.4129.5
3016.738.160299.915565.5629.185.18392.793.01
4016.3210.37051316522.2228.356.74365.2189.8
591.6716.069.25039917394.4426.305.65334.8119.2
6161.1116.118.8653.78361.61881.6723.936.78247.9192.6
7245.5616.669.5595.25428.319021.555.4849.25109.8
8333.3317.908.77146.8353.420019.506.220154.9
9341.6719.917.43202.124121018.035.960138.3
10391.1122.425.90208.8134.522017.096.110147.7
11457.2225.056.34242.4162.423016.526.810194.9
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Song, J.; Lu, C.; Ma, Q.; Zhou, H.; Yue, Q.; Zhu, Q.; Zhao, Y.; Fan, Y.; Huang, Q. Distributed Integrated Synthetic Adaptive Multi-Objective Reactive Power Optimization. Symmetry 2022, 14, 1275. https://doi.org/10.3390/sym14061275

AMA Style

Song J, Lu C, Ma Q, Zhou H, Yue Q, Zhu Q, Zhao Y, Fan Y, Huang Q. Distributed Integrated Synthetic Adaptive Multi-Objective Reactive Power Optimization. Symmetry. 2022; 14(6):1275. https://doi.org/10.3390/sym14061275

Chicago/Turabian Style

Song, Jiayin, Chao Lu, Qiang Ma, Hongwei Zhou, Qi Yue, Qinglin Zhu, Yue Zhao, Yiming Fan, and Qiqi Huang. 2022. "Distributed Integrated Synthetic Adaptive Multi-Objective Reactive Power Optimization" Symmetry 14, no. 6: 1275. https://doi.org/10.3390/sym14061275

APA Style

Song, J., Lu, C., Ma, Q., Zhou, H., Yue, Q., Zhu, Q., Zhao, Y., Fan, Y., & Huang, Q. (2022). Distributed Integrated Synthetic Adaptive Multi-Objective Reactive Power Optimization. Symmetry, 14(6), 1275. https://doi.org/10.3390/sym14061275

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