IMODBO for Optimal Dynamic Reconﬁguration in Active Distribution Networks

: A dynamic reconﬁguration method based on the improved multi-objective dung beetle optimizer (IMODBO) is proposed to reduce the operating cost of the distribution network with distributed generation (DG) and ensure the quality of the power supply, while also minimizing the number of switch operations during dynamic reconﬁguration. First, a multi-objective model of distribution network dynamic reconﬁguration with the optimization goal of minimizing active power loss and voltage deviation is established. Secondly, the K-means++ clustering algorithm is used to divide the daily load of the distribution network into periods. Finally, using the IMODBO algorithm, the distribution network is reconstructed into a single period. The IMODBO algorithm uses the chaotic tent map to initialize the population, which increases the ergodicity of the initial population and solves the problem of insufﬁcient search space. The algorithm introduces an adaptive weight factor to solve the problem of the algorithm easily falling into a locally optimal solution in the early stage with weak searchability in the later stage. Levy ﬂight is introduced in the perturbation strategy, and a variable spiral search strategy improves the search range and convergence accuracy of the dung beetle optimizer. Reconﬁguration experiments on the proposed method were conducted using a standard distribution network system with distributed power generation. Multiple sets of comparative experiments were carried out on the IEEE 33-nodes and PG&E 69-nodes. The results demonstrated the effectiveness of the proposed method in addressing the multi-objective distribution network dynamic reconﬁguration problem.


Motivation
The prevalence of low-carbon power system networks is a trend that conforms to the development of new energy technologies. Such systems also represent a meaningful way to promote the high-quality development of clean energy [1][2][3]. Owing to the accessibility of distributed renewable energy [4], instability is associated with a number of non-traditional safety problems in the stable operation of the distribution network systems. By manipulating section switches and tie switches, distribution network reconfiguration can alter the network topology and optimize power transmission within the distribution network system.

Literature Review
Domestic and foreign scholars have been devoted to researching algorithms to solve the distribution network reconfiguration problem. Existing algorithms are classified into three categories: traditional mathematical optimization algorithms [5,6], heuristic algorithms [7,8], and intelligent optimization algorithms [9][10][11]. A reconfiguration model needs to be established for the distribution network before solving the problem using traditional mathematical optimization algorithms that only solve single-objective models. Given the increasing complexity of distribution network systems in recent years, the data involved in the distribution network have become more extensive, resulting in long computing times and difficulty in solving the distribution network reconfiguration problem. Heuristic algorithms establish some effective rules to solve the distribution network reconfiguration problem by analyzing the situation. These algorithms include branch exchange algorithms, optimal flow algorithms, etc. However, the optimization abilities of heuristic algorithms are poor.
Various optimization methods have been proposed to solve the distribution network reconfiguration problem. Li et al. [12] proposed a dynamic reconfiguration integrated optimization method based on a multi-objective sparrow search algorithm, optimizing multiple objectives, including power quality, economic efficiency, and energy loss, and comparing various algorithms, has achieved good results. Wang et al. [13] proposed a distribution network reconfiguration method that considers the network's reliability, economy, and environmental friendliness based on a parallel slime mold algorithm. When solving multi-objective optimization problems, using a weighted method to transform multiple objectives into a single objective may result in the inability to obtain the actual optimal value in a complex system. This is because the weights are preset. Gao et al. [14] proposed a multi-objective optimization method that effectively finds the optimal energy configuration scheme for intelligent communities using Levy flight and an improved chicken swarm algorithm, which enhances the algorithm's global and local search abilities. In addition, a comprehensive set of scenario simulations was conducted to test the model. The results demonstrate that the improved algorithm effectively overcomes the issue of being trapped in local optima. Mahdad et al. [15] proposed an optimal reconfiguration and reactive power planning method based on a fractal search algorithm, which can simultaneously consider multiple indicators and has high search efficiency and optimization accuracy. Kefayat et al. [16] proposed a hybrid approach that combines ant colony optimization and an artificial bee colony algorithm to optimize the placement and sizing of distributed energy resources with probabilistic constraints. The system dynamically switches between the two algorithms through an adaptive strategy to achieve better optimization performance. Jafari et al. [17] proposed a parallel implementation of the exchange market algorithm (EMA) and wild goats algorithm (WGA), which improved the computational speed and accuracy. However, their study on distribution network systems did not consider the presence of distributed energy resources. Parizad et al. [18] proposed a novel algorithm based on particle swarm optimization (PSO) for distribution feeder reconfiguration (DFR). The algorithm aims to minimize short-circuit levels in power distribution systems and is validated using the IEEE 83 bus power distribution system model, demonstrating its effectiveness. These optimization methods demonstrate the diverse strategies used to solve distribution network reconfiguration problems, and each method has its strengths and limitations depending on the problem being addressed.
Numerous researchers have proposed intelligent optimization algorithms by analyzing the habits and physical phenomena of various biological populations. Algorithm technology is a continuously evolving field that undergoes frequent updates and iterations. In 2022, Jiankai Xue introduced a new swarm intelligence optimization algorithm called the dung beetle optimizer (DBO) [19]. The DBO algorithm employs multiple search strategies and can perform local and global searches, demonstrating a fast convergence speed and high solution accuracy. The DBO algorithm has vast potential for application in various fields. The development of this algorithm represents a significant contribution to the advancement of optimization algorithms.

Contributions
In this study, we propose a dynamic reconfiguration method for distributed generation and distribution networks using the improved multi-objective dung beetle optimizer (IMODBO). This method establishes a multi-objective model with active power loss and voltage deviation as objective functions. According to the characteristics of the problem, the K-means++ clustering algorithm is used to divide the daily load into periods, and the IMODBO algorithm is used to transform the single-period distribution network. The IMODBO algorithm uses a tent chaotic map in the population initialization process to improve the quality and ergodicity of the initial population. The adaptive weight factor solves the problem of early convergence and weak searchability of the algorithm. Levy flight perturbation and variable spiral strategy enhance the global search ability of the algorithm. The validity of the proposed method was verified by a reconstruction experiment on the IEEE 33-node and the PG&E 69-node power distribution system. Against previous studies, the main contributions are summarized below.

•
The dynamic distribution network model with distributed power generation was established, aiming to minimize the active power loss and voltage deviation. • In response to the multi-objective optimization problem of distribution network reconfiguration, a improved multi-objective dung beetle optimizer (IMODBO) is proposed. Compared with various state-of-the-art single-objective and multi-objective algorithms, IMODBO demonstrates superior convergence speed and exploration capability. • The distribution network reconfiguration problem is solved using IMODBO combined with the K-means++ clustering algorithm. By comparing with existing algorithms and conducting multiple scenario tests on the IEEE 33-node and PG&E 69-node standard distribution network systems, significant improvements are achieved. • The experimental results demonstrate that the proposed method effectively reduces active power losses, stabilizes node voltage distribution, and significantly reduces the number of reconfigurations required.
The organization of this study is as follows: Section 2 introduces the dynamic reconfiguration optimization model for power distribution systems, including the constraints of distributed power generation. Section 3 presents the application of the K-means++ clustering algorithm, the improvement measures of IMODBO, and the optimization process for solving the distribution network. Section 4 validates the effectiveness of the improved algorithm through various test functions and verifies the proposed method in the different scenarios of two distribution network systems. Section 5 concludes the findings of this study, analyzes the limitations of the proposed method, and presents future research directions.

Problem Formulation
The proposed research methodology utilizes the IMODBO to optimize the active power loss and voltage deviation in distribution network systems. The algorithm is applied to determine the optimal configuration of the distribution network. The effectiveness of the proposed algorithm is validated using the IEEE 33-node system and the PG&E 69-node system.

Objective Function
To improve the transmission efficiency of the distribution network, the active network loss of the reconstructed distribution network is optimized to reduce the distribution network's power cost and enhance its reliability and stability. The minimum active power loss of the distribution network after reconstruction can be expressed as follows [20]: where T represents the reconfiguration interval; G represents the total number of branches in the distribution network system; k t,l represents the first state at the l-th branch during the t-th period, where k t,l = 1 means that the l-th branch during the t-th period is closed, and k t,l = 0 means that the l-th branch during the t-th period is disconnected; P t,l and Q t,l represent the active power and reactive power at the l-th branch during the t-th period, respectively; R l is the resistance of the l-th branch; and U t,l represents the node voltage at the end of the l-th branch during the t-th period.
A voltage offset model of the system is established to balance the voltage distribution of each node in the distribution network. Reducing the voltage offset can improve the stability and reliability of the distribution network, reduce energy waste, and reduce the incidence of equipment failure. The minimum voltage deviation of the distribution network after reconfiguration can be expressed as [21] where U t,i represents the actual voltage value at the i-th node during the t-th period, N represents the total number of nodes in the model, and U * i represents the rated voltage of the i-th node.

Constraint Conditions
In order to demonstrate the feasibility of the reconfiguration, the process must adhere to certain constraints, including power flow, nodal voltage, branch current, and distributed generation power.
(1) Equality constraint of the power flow equation [22]  where i and j represent the start and end nodes of the branch, respectively; P t,i , P t,DGi , and P t,loadi represent the power injected into the network, the distributed power output, and the active power at the i-th node during the t-th period, respectively; Q t,i , Q t,DGi , and Q t,loadi represent the network injected power, distributed power output, and reactive power at the i-th node during the t-th period, respectively; U t,i and U t,j represent the voltage of the i-th and j-th node during the t-th period, respectively; and δ t,ij , G ij , and B ij represent the phase angle difference, susceptance, and conductance between nodes i, j, respectively.
(2) Node voltage constraint [23] where U i,min and U i,max represent the lower and upper limits of the allowable voltage at the i-th node, respectively, and U min i and U max i take 0.9U t,i and 1.05U t,i , respectively. (3) Branch current constraints [24] I l ≤ I l,max where I l and I l,max represent the maximum value of the current flowing through the l-th branch.
(4) Distributed power generation constraint [25] P DGi,min ≤ P t,DGi ≤ P DGi,max Q DGi,min ≤ Q t,DGi ≤ Q DGi,max where P t,DGi and Q t,DGi represent the active and reactive power injected by DG at the i-th node during the t-th period, respectively; P DGi,max and P DGi,min represent the upper and lower limits of the active power that the DG can supply to the i-th node, respectively; and Q DGi,max and Q DGi,min represent the upper and lower limits of the reactive power that the DG can supply to the i-th node, respectively.

Distributed Power Generation Model
(1) Wind Power Output The wind power generation model is established according to [26], and the relationship between wind speed and output power at different times is shown in Equation (7).
where P r is the rated active power of the generator, v t is the actual wind speed on site, v ci is the minimum wind speed of the wind turbine, v co is the maximum wind speed of the wind turbine, and v r is the rated wind speed of the wind turbine.
(2) Photovoltaic Power Output According to reference [27], the power output of photovoltaic cells is mainly related to the intensity of solar irradiance, the area of solar photovoltaic panels, and the photoelectric conversion efficiency. As the solar irradiance intensity follows a Beta distribution, its probability density function can be described as: where p t represents the light intensity during period t, P max represents the upper limit of the power output capacity of the most comprehensive PV array in a given system, and ε and β represent parameters conforming to the Beta distribution. The photovoltaic output power can be expressed as where U is the total number of photovoltaic panels, A u represents the u-th cell the area of the plate, and η u represents the conversion efficiency of the battery plate.

Time Division Method Based on the K-Means++ Clustering Algorithm
To address the issue of frequent reconfigurations in dynamic distribution networks, we employs the K-means++ [28,29] clustering algorithm to divide the daily load into periods, effectively reducing the active power loss and voltage offset and minimizing the number of reconfigurations.
Since the load status of each node in the distribution network changes dynamically with time, the load within 24 h is considered a dataset, with each dimension corresponding to the total number of nodes in the system. Using 1 h as a period, the node load of each period is treated as a sample set (O) for cluster analysis, representing the apparent power of each node in the distribution network system within that period. During dynamic reconfiguration optimization, reconfiguration is carried out continuously, and interperiod reconfigurations must meet actual needs. Thus, distance segment constraints need to be added to cluster analysis, which is achieved as follows: Step 1: Determine the collection dataset (O) using the elbow method to determine the optimal range of clustering, i.e., the K value.
Step 2: Randomly select a value from the sample dataset (O) as the initial value center and set the initial value of K to 1.
Step 3: Compute the Euclidean distance from each data object to the center of K. Using the roulette method, select the next cluster center.
Step 4: Repeat Step 3 until K initial value centers are generated.
Step 5: Divide the 24 moments in the sample dataset (O) into the nearest initial value center to form a cluster. Calculate the mean value of each dataset as the new initial value center.
Step 6: Repeat Step 5 until the sum of the squares of the distances from each data point to the corresponding cluster remains unchanged. Then, output the period in each group.

Dung Beetle Optimizer
The dung beetle optimizer takes inspiration from the biological behavior of dung beetles. The optimal solution is achieved by continuously updating the position updating formulas of four dung beetles, each with its position renewal formula, as follows: The dung beetles that roll balls are updated and can be expressed as Encountering obstacles, being unable to move forward, and repositioned dancing can be expressed as The position-updating formulas for female dung beetles can be expressed as: The position-updating formulas for small dung beetles can be expressed as: The position-updating formulas for thief dung beetles can be expressed as: The values of k 1 and k 2 are fixed, with k 1 set to 0.1 and k 2 set to 0.3. The flexure angle (tan θ) belongs to the interval [0, π] and is used in the position-updating formulas of the four dung beetles, where X * and X b represent the current best position and the global best position, respectively. The oviposition and reproduction areas have upper and lower limits represented by Ub * , Lb * , Ub b , and Lb b . The constants b 1 and b 2 represent independent random vectors of size 1 × dim. The variable C 1 is a random number drawn from a normal distribution, while C 2 is a random vector with values ranging from 0 to 1. Additionally, a random vector of size 1 × dim that follows a normal distribution is represented by g, and S represents a constant value. The four dung beetle position-updating formulas are constantly updated until the optimal solution is found.

IMODBO Based on Mixed Strategy
After partitioning the data into time intervals, IMODBO is used to solve the distribution network reconstruction problem for each interval.

Population Initialization Based on a Tent Map
The standard dung beetle optimizer uses random initialization, which cannot guarantee the diversity of the population, affecting the quality of the global optimal solution. Using chaotic mapping to generate individuals randomly, the algorithm can effectively improve the speed and accuracy of the search for the optimal solution. As mentioned in [30], through many comparative experiments, chaotic tent mapping has better bias adaptation. The iteration equation of chaotic tent mapping is expressed as follows: Equation (15) shows that the chaotic tent mapping equation has fewer parameters, making the application convenient and straightforward. The population size is set to 500, the population dimension is to 2, and the parameter s is set to 0.5. The initial population distribution generated by chaotic tent mapping in the interval [−100, 100] is shown in Figure 1. The initial population generated by chaotic tent mapping is distributed uniformly in space. Using these chaotic sequences, the initial population pop of the dung beetle algorithm can be obtained according to Equation (16).
In Equation (16), ub and lb are the upper and lower limits of the population.

Adaptive Weight Factor
According to the position-updating formula of the standard dung beetle optimizer, the four dung beetles strongly depend on the best position in the position-updating phase. Therefore, an adaptive weight factor is introduced in during the position updating, the formula of which is expressed as follows [31]: where W max is the maximum value of the adaptive weight, W min is the minimum value of the adaptive weight, T max represents the maximum number of iterations of the algorithm, and t * indicates the current iteration number. After testing, the best algorithm performance is achieved when W max = 0.9 and W min = 0.4. As the number of iterations increases, the adaptive weight decreases from 0.9 to 0.4. At the beginning of the iteration, a more considerable adaptive weight improves the exploration ability of the algorithm. As the number of iterations increases, a smaller adaptive weight in the later stage improves the searchability of the algorithm. This prevents falling into a local optimum early, and the convergence in the last stage can gradually approach the optimal value. The position-updating equations are obtained by adding adaptive weighting factors into Equations (10)- (14), which can be expressed as follows: The dung beetles that roll balls are updated and can be expressed as: The position-updating formulas for female dung beetles can be expressed as: The position-updating formulas for small dung beetles are expressed as follows: The position-updating formulas for thief dung beetles are expressed as follows:

Levy Flight Disturbance Strategy and Variable Spiral Search Strategy
According to the characteristics of the ball-rolling dung beetle algorithm, the ballrolling dung beetles can expand their search range by "dancing". However, the search range of the three remaining dung beetles is relatively simple, so the Levy flight perturbation strategy and the variable spiral search strategy should be introduced to increase the global optimal search range. The Levy flight perturbation strategy is introduced, which enables the algorithm to search randomly at different distances, maximizes the diversification of the search domain, and enhances the optimal global searchability. Levy flight not only satisfies the small-scale refined search, ensuring that each target can be searched, but also satisfies the large-scale rough search, avoiding the limitations of local search. The distribution density function of the Levy flight step size change is expressed as [32]: where f is the motion step of Levy flight, which can be expressed by Equation (24) L v = µ |ν| 1/∂ (24) where µ and v are random numbers that conform to a normal distribution, and σ µ and σ r are obtained by Equation (25): where the value range of parameter ∂ is 0 < ∂ < 2, and generally, ∂ = 1.5.
The variable helical search strategy allows for the development of various positionupdating search paths, balancing the global and local searches of the algorithm. In the position-updating process, the parameter z is designed as an adaptive variable to solve the monotony of the spiral parameter in the search method in case of falling into a locally optimal solution. It is used to dynamically adjust the spiral shape of the three dung beetle searches, thereby broadening the ability of three dung beetles to explore unknown areas and improving the algorithm's search efficiency and global search performance. The formula for the variable helical search strategy is as follows [33]: where Q is the spiral factor, and h is the random number of (−1, 1), where z = e k·cos((t * /T max )·π) , and k = 5.
In summary, the position updates of the four improved dung beetles are expressed as follows: The position-updating formulas for the dung beetles that roll balls are updated and can be expressed as: The position-updating formulas for female dung beetles can be expressed as: The position-updating formulas for small dung beetles are expressed as follows: The position-updating formulas for thief dung beetles are expressed as follows:

Pareto Dominance Theory
The traditional dung beetle optimizer selection strategy is suitable for single-objective optimization problems but not for multi-objective distribution network optimization and reconstruction problems. Therefore, in this paper, we use the Pareto dominance theory to evaluate the solutions generated by the distribution network optimization and reconstruction problem. Individuals with high degrees of optimization and extensive coverage ranges are selected to form the population for the next iteration. According to the Pareto theory, the objective function of this iteration is ranked by non-dominant levels. Every time a fitness calculation is performed, it is compared with the solution set. The solution with the highest non-dominant level is stored in the external archive as the current optimal solution, and solutions with lower non-dominant levels are removed. This process is repeated until the fitness calculation is completed.
Step 1 involves reading the distribution network data and determining the loop and the branch switches of each loop participating in the coding.
Step 2 sets the initial parameters, including the number of populations, the maximum number of iterations (T), the population dimension (dim), and the upper and lower limits of each dimension component of IMODBO based on the number of branch switches of each loop participating in the encoding.
Step 3 generates the population using the chaotic tent map according to Equation (16).
Step 4 calculates the fitness of each individual in the population, performs Pareto screening, selects the solution with the highest non-domination level, and stores it in an external file.
Step 5 divides the population into four types of dung beetles and updates the positions of the four dung beetles according to Equations (27)-(31).
Step 6 recalculates the fitness of the updated positions, compares them with the external file using greedy selection, and updates the external file.
Step 7 checks whether the current iteration number of the algorithm reaches T. If not, the algorithm proceeds to Step 8; otherwise, it returns to Step 4.
Step 8 outputs the optimal solution from the external file. In summary, the steps of the improved multi-objective dung beetle algorithm for solving distribution network optimization and reconstruction problems are shown in the Figure 2.

Result and Discussion
In this study, we conducted performance testing on the algorithm using both singleobjective and multi-objective test functions. A comparison was made between IMODBO and existing algorithms. Furthermore, a comparative analysis was performed on two distribution network systems through multiple experimental comparisons.

Name
Function Interval

Test Functions and Performance Evaluation
The multi-objective intelligent algorithm was evaluated with three evaluation indices. The specific calculation methods for each index are as follows: (1) Inverted Generational Distance Indicator The inverted generational distance (IGD) indicator is commonly used to evaluate the convergence and diversity of the non-dominated solutions obtained by the algorithm. It measures the average Euclidean distance between all solutions in the Pareto front and the non-dominated solutions. A low IGD value indicates the better convergence and diversity of the algorithm. The calculation formula of IGD is as follows [45]: |pop| (32) where P is the distribution of a real set of points on the Pareto frontier, |pop| represents the number of set points pop, P * for the algorithm to obtain the Pareto optimal solution set, and d(h, P * ) represents the minimum Euclidean distance from a point in the set.
(2) Hypervolume Indicator The hypervolume (HV) metric is used to evaluate the quality of the non-dominated solution set generated by the algorithm concerning a reference point. A high HV value indicates better comprehensive algorithm performance in developing well-distributed and diverse non-dominated solutions. The calculation formula of HV is as follows [46]: where λ represents the Lebesgue measure, and v i represents the super volume formed by the reference point.
(3) Spacing To evaluate the diversity of the solution set, the standard deviation of the minimum distance from each solution to the other solutions is calculated, which is commonly referred to as the spacing metric. A low spacing value indicates a more uniform distribution of solutions in the objective space. The calculation formula of spacing is as follows [47]: where d i and d represent the minimum distance from the i-th solution to other solutions and the average value of d i , respectively. The single-objective test function was calculated, and the results of the experiment are presented in Table 2. Additionally, the iteration curve is depicted in Figure 3.

Analysis of Test Results
According to the data in Table 2, the performance of IMODBO on different test functions far exceeds those of other algorithms, and according to Figure 3, IMODBO only required a small number of iterations to find the optimal solution. The results show that the introduction of different strategies can greatly improve the optimization accuracy of the algorithm and improve the performance of the algorithm.
The calculation results of the three indices are shown in Tables 3-5. The Pareto front solutions of seven different multi-objective algorithms on each test function and the accurate Pareto front distribution of the test function are shown in Figure 4. Tables 3-5 show three comprehensive indices of each algorithm under different test functions.According to the IGD index, IMODBO performs well on all four test functions, with only slightly inferior performance to that of MOWOA on the ZDT4 function, indicating that IMODBO has a strong searching ability. According to the HV index, the results obtained by IMODBO are similar to those obtained by MOGWO and MOSMA but better than the results obtained by the other algorithms, indicating that IMODBO has a rich solution diversity and the best all-around performance. From the perspective of the SP index, IMODBO has the best solution distribution space, widest solution set distribution, and best distribution performance. Based on the performance indices obtained from the test function, the IMODBO algorithm has apparent advantages in terms of distribution and convergence.

Reconstruction of the IEEE 33 and PG&E 69 Power Distribution Systems
We conducted simulation verification tests using MATLAB R2022a on the IEEE 33 [48] and PG&E 69 [49] power distribution systems. The proposed method was validated on the IEEE 33-node power distribution system to demonstrate its effectiveness and was further applied to the PG&E 69 distribution system to showcase its universality. Two sets of comparative experiments were conducted on each distribution network system as follows: • The K-means++ clustering algorithm was employed to divide the distribution network load dataset into different time periods, facilitating comparative analyses under various scenarios. • A comprehensive test was conducted during a single time period to compare the performance of IMODBO with various other algorithms.

IEEE 33-Node System Reconstruction
As shown in Figure 5, the IEEE 33-node system contains five contact switches, 37 branches, and 33 nodes. Wind power generation was connected to nodes 8 and 18, and photovoltaic cells were connected to node 32. The load added to DG was calculated according to Tables 6 and 7. Table 6. Photovoltaic cell parameters.

Solar Cell
A/m 2 η p max /(W/m 2 ) PV 4500 15 20 The dynamic reconfiguration analysis of the distribution network was carried out, and the load distribution of the distributed power supply was added, as shown in Figure 6. The load is different at each moment and each node, which includes commercial, residential, and industrial loads. A load of each period was calculated according to Equation (35) to simulate the dynamic operation of the distribution network. The daily load variation curve of the IEEE 33-node distribution network system is shown in Figure 7.
where i is the node number, t is the time point, h is the load type, L i max is the peak load of the i-th node, H i,h is the proportion of the class h load in the i-th node, and C h,t is the proportion of the class h load during the t-th period. We take 110% of the standard load peak value of IEEE 33-nodes. The period was divided using the elbow method and the K-means++ clustering algorithm. The inflection point was determined to be 4, and the time period was divided into four class clusters, as shown in Figure 8.    Three scenarios were set up for comparative analysis: Scenario 1: Maintain the original system structure without reconfiguration; Scenario 2: Use IMODBO to reconfigure once every period; Scenario 3: Dynamic reconfiguration via IMODBO. The final operation effect is shown in Table 8. voltage offset are relatively high. Dynamic reconfiguration using the IMODBO proposed in this paper can effectively reduce the functional network loss and voltage offset. Before reconfiguration, the active network loss is 1492.14 kW; after reconfiguration, the active network loss is 1105.74 kW. The voltage offset before reconfiguration is 17.3396 p.u, and after reconfiguration, the voltage offset is 12.3715 p.u, which is 30.2% lower than without reconfiguration. The effect is remarkable. A comparison of Scenario 2 with Scenario 3 shows that the daily load in Scenario 2 is divided into 24 periods and reconstructed once each period, leading to a slower reconfiguration speed and more switching actions. In Scenario 3, the active power network loss is 15.34 kW lower than that in Scenario 2. The voltage offset decreases by 0.22 p.u from 12.3715 p.u to 12.1026 p.u.
The above results show that the method proposed in this paper reduces the number of switching actions and ensures calculation accuracy. Figure 9 shows that Scenario 3 can effectively minimize the work packet loss in most periods. Figures 10-12 show that the node voltage distribution in Scenario 3 is more uniform.

PG&E 69-Node System Reconstruction
The PG&E 69-node distribution system is shown in Figure 13. Wind power generation is connected to nodes 33 and 64, and photovoltaic cells are connected to node 14. The load of the connected distributed power supply is the same as previously described. The daily load curve of the PG&E 69-node distribution system is shown in Figure 14.
The time period is divided according to the algorithm described in Section 2, setting the same scene as the IEEE 33-nodes. Table 9 shows the final running effect.  According to Table 9, in Scenario 3, compared with Scenario 1 and Scenario 2, the active power network loss decreases by 1194.37 kW and 13.29 kW, respectively, and the voltage offset decreases by 12.4101 p.u and 1.343 p.u. After reconfiguration, the effect is significantly improved. According to the comparison of the network loss distribution of different scenarios in different periods in Figure 15, Scenario 3 can effectively reduce the network loss value in most periods compared with Scenario 2 and Scenario 1. Figures 16-18 show that the voltage distribution of each node at each time in Scenario 3 is more uniform than that in Scenario 1 and Scenario 2.

Comparison Test of Different Algorithms
To verify the effectiveness of the IMODBO in the dynamic reconfiguration of a distributed power distribution network, the NSGA-II, MOALO, MOPSO algorithms were used in this study to compare the solution results of the IEEE 33-node and PG69-node system models and distributed power supply.
Because the real Pareto solution set in the system cannot be known, the four algorithms were run independently ten times, and their Pareto solution sets were taken for nondominated sorting to simulate the natural Pareto frontier of the IEEE 33-node distribution network. Taking the distribution network load in the eighth hour as an example, the Pareto solution set obtained after the four-stage algorithm was run, as shown in Figure 19. To further demonstrate the effectiveness of the approach, the Pareto solution set for the PG&E 69-node distribution network during the fourteenth hour is also presented. Moreover, two multi-objective performance indicators are calculated, as shown in Table 10, through which the performances of different algorithms can be judged.  In terms of the IGD index, IMODBO demonstrates the lowest value among the four algorithms, indicating its superior convergence and diversity. This signifies that IMODBO is capable of finding solutions that are closer to the true Pareto front. Additionally, the HV index also shows that IMODBO achieves the highest values, outperforming the other algorithms by a significant margin. This confirms the exceptional performance of IMODBO in generating solutions that cover a wide range of objectives. In summary, the multiobjective distribution network dynamic reconfiguration method based on the IMODBO proposed in this paper exhibits significant advantages in terms of convergence and diversity. It showcases superior performance compared to MOALO, NSGA-II, and MOPSO.

Conclusions
This study proposes an optimization method based on the improved multi-objective dung beetle optimizer (IMODBO) to address the dynamic reconfiguration optimization problem in distribution networks with distributed power sources. The IMODBO method, in conjunction with the K-means++ clustering algorithm, is proposed and validated for this problem. Furthermore, multiple scenarios are designed in standard distribution network systems to verify the effectiveness of the proposed method. The main contributions of this study are as follows: • Through testing on both single-objective and multi-objective test functions, the IMODBO algorithm outperforms other optimization algorithms such as DBO, ISSA, NGO, WOA, and GWO in terms of convergence rate, convergence accuracy, and overall performance in single-objective optimization. In multi-objective optimization, the IMODBO algorithm exhibits better performance than MOALO, NSGA-II, MODA, MOSMA, MOPSO, and MOWOA in terms of convergence, solution set distribution, and comprehensive performance. • An optimization scheme is proposed for dynamic reconfiguration in distribution networks with distributed power sources, utilizing the IMODBO algorithm and K-means++ clustering algorithm. The objective is to reduce the active power loss, stabilize node voltages, and minimize switch operations. • In the IEEE-33 nodes and PG69 nodes test systems, the proposed method achieves the minimum active power loss and voltage deviation among multiple scenarios. In a single-period comparison experiment, the Pareto solutions obtained by the IMODBO algorithm dominate over the solutions obtained by NSGA-II, MOALO, and MOPSO algorithms. The proposed method also exhibits good performance in terms of multiobjective performance metrics. The Pareto solutions obtained by the IMODBO algorithm have the minimum IGD and maximum HV values in both test systems. This indicates that the proposed method can simultaneously seek solutions that are close to the optimal Pareto front and provide diverse options for decision makers.
This study has significant implications for achieving sustainable energy goals. The proposed IMODBO method demonstrates superior optimization performance and provides an effective solution to reduce energy loss and improve power quality in distribution networks. However, there are limitations to this study. For instance, it does not consider the randomness of photovoltaic power generation and wind power generation. Future research should take into account more factors to make the results more applicable to real-world projects. Additionally, testing in larger and more complex systems with a wider range of distributed power sources can be conducted to further increase the penetration rate of renewable energy and achieve carbon reduction and neutrality goals.