# An Optimization Method of Active Distribution Network Considering Time Variations in Load and Renewable Distributed Generation

Authors: Juan Wen; Xing Qu; Siyu Lin; Lin Ding; Lin Jiang
Journal: International Transactions on Electrical Energy Systems
 (2022)

---

## Abstract

Network optimization is one of an effective ways to enhance the performance of an active distribution network (ADN). Aiming to improve the operation and power quality of the ADN considering time variations in load and renewable distributed generation (RDG) power, a multi-time period optimization model and its dynamic solution method are proposed. Considering the real time load demand and power generation variation of RDG versus input parameters like wind speed and solar irradiance, the time variation models of load and RDG power output are developed. The minimum power loss and maximum absorption of RDG power are served as the optimization indexes to construct the dynamic muti-time period optimization model. A hybrid particle swarm optimization (HPSO) algorithm is presented based on integer coding and random coding technique, which can find the most satisfactory solutions for the proposed dynamic model. Considering the time variation of load and RDG power of ADN, the optimal network structure and RDG allocation scheme at any time interval are determined by analyzing the obtained solutions. Additionally, two ADNs with time variation in load and RDG are tested to verify the effectiveness and superiority of the proposed dynamic optimization model and HPSO algorithm. The simulation results show that the proposed method can improve the operation performance and RDG optimal utilization of the ADNs through muti-time period dynamic optimization.


---

## Body

## 1. Introduction

The application of RDG is a topic that attracts a great deal of interest in the electric power industry because it can reduce the dependency on fossil fuels and environmental pollution [1]. The integration of RDG in the ADN is growing rapidly in discussions about the future of distribution systems with the increased demand for electrical power and the requirement of environmental conservation [2, 3]. The share of RDG in primary energy supply would rise from 14% in 2015 to 63% in 2050 [4]. Recently, wind and solar photovoltaic (PV) RDG technologies are being integrated with ADN because they are easy to install, low operating costs, and mature technology [5, 6]. Taking wind power generation as an example, wind energy would share 15.7% of global electric power consumption by 2020 [6]. The integration of RDG would improve the operation performance and power supply shortage of ADN. Nevertheless, the high-level penetration of RDG brings new challenges affecting the dispatching operations of ADNs [7]. Moreover, The output power of RDGs exhibits fluctuations and intermittency because it depends on input parameters, which places high requirements on the adaptability and control of ADNs. For example, the loss of a single 400-kilovolt transmission line in the India blackout causes sizeable voltage deviations and around 300 million customers without power [8]. Various control actions have been performed to tackle the uncertainties in ADNs, such as network optimization, controllable loads, etc., [9, 10].Network optimization is a critical issue in distribution systems dispatching management. One effective optimization way is network reconfiguration which is implemented on ADN in the presence of RDG. The method is the process of finding the optimal topology to satisfy the operational objective and constraints [11]. There are some operational problems with the integration of RDG [12], such as economic dispatch, voltage oscillations, and harmonic distortion. Furthermore, operational conditions of ADN will become more severe when the characteristic of time-varying load conditions and the uncontrollability of RDG output power are both considered [7]. The operational conditions of ADN change frequently along with the time-varying characteristics of RDG and load. In the ADN, flexible and effective control measures are required to apply the frequently changing of operational conditions. We aim to analyze a potential network optimization scheme for improving the performance of the ADN.In [13], Merlin and Back first propose the concept of network optimization which is obtained the optimal topological configuration with the objective of minimum power loss. And then, the objective functions are extended as load balance [14], voltage deviations [15], power quality [16], and so on. The network optimization is defined as an optimization problem which should satisfy ADN operational objectives and constraints. The obtained optimal solutions of the optimization problem is a challenging task because the search space of solutions is typically large. Many intelligent evolutionary approaches are adopted to address the combinatorial problem such as genetic algorithm (GA) [17], improved decomposition based evolutionary algorithm (I-DBEA) [18], particle swarm optimization (PSO) [19], and so on.The integration of RDG brings many potential benefits to ADNs, such as minimizing power losses, improving voltage profile, enhancing system stability etc [20]. However, the fluctuation of RDG output in the ADN results in node overvoltage swells and reverse power flow of the system. For instance, reference [22] has shown the changes in ADN stability and economic indicators while the RDG sizing and siting are changed. For this reason, researchers have shown interest to obtain an optimal management scheme for RDG in ADNs using network optimization. In [23, 24], the optimization model with the objectives of RDG units sizing and sitting are established which allocates RDG to the optimal places of ADN. In [25], an energy management scheme of dispatchable wind and PV RDG is proposed to balance the power generation and load demand. It is noteworthy that the methods have been subjected to static network optimization theory and ignored the time variations in RDG. The static optimization model for non-variable RDGs cannot demonstrate real-time scenarios and cannot achieve the optimal RDG configuration scheme considering time scheduling of the ADNs.The RDG output power has strong uncertainty and time-varying because it is mostly influenced by wind speed and solar irradiance [26]. These characteristics bring new challenges to the calculation and evaluation of the ADN network optimization problem [27]. In [28], a forecast model is used to obtain the RDG variable power at future time interval. In [29], a planning model of wind and solar generations considering uncertainty is proposed to obtain an optimal mix scheme and sizing of various RDG in microgrids. However, The mentioned studies have been performed based on non-variable load models and addressed the network optimization during a predetermined time interval.Considering time scheduling of the load in ADN, reference [30] presents the optimization model of different time horizons with aiming to reduce switching costs. In [31], the optimally allocate RDG units have been investigated in ADN with time-varying load demand. The model considering load demand and RDG is proposed to solve the ADN optimization problem over a 24-hour time horizon [32]. An optimization model based on RDG fuzzy uncertainty is presented to obtain the optimal topology with objectives of power loss reduction and voltage stability improvement [33]. However, the obtained strategy is not suitable for online real-time network optimization due to the load and RDG based on a fixed prediction. In [34], the optimal management scheme of wind RDG in ADN is proposed on the basis of the probabilistic generation assumption. The solution for the optimization model in the presence of time variation in load and wind RDG is not evaluated markedly.Reference [35, 36] investigates the operation performance of the ADN based on the conditions of load variation and wind RDG penetration level. The obtained results have not been suitable for the ADN with the integration of multi-type RDG simultaneously. In [37, 38], the optimization methods considering multiple types of RDGs and load demand are proposed to improve the operational management of the ADN. However, the obtained solutions have the coordination between the load and RDG. During the process of ADN network optimization, operational and calculation conditions would become complicated with network size, the type and number of RDGs, and the length of the operation horizon. In addition, the above methods did not evaluate the operational indicators of the optimized ADN from the overall perspective.In an ADN with high penetration of RDGs, applying actual operation data to model time variation and fluctuation behaviours of the RDGs and loads is a practical solution for distribution operations. An ideal optimization model should consider dynamic optimization of network topology, optimal management scheme for RDGs, and satisfying demand in a reliable and stable way considering time variations in RDGs and loads. However, the operational and calculation conditions would become complicated since the characteristic of time variation in load and RDG during the network optimization of ADN. It is obviously that the mathematical methods are very hard to solve this optimization problem. Furthermore, the voltage oscillations and fault currents in the ADN may increase after RDG integration. The mentioned problems would be more severe in the process of optimization considering the time-varying load and RDG conditions. Although a countermeasure is presented to restrict the output power of RDG by installing control devices, it is not benefit to the effective utilization of RDG. Moreover, due to the characteristic of time variation in load and RDG, the optimization problem size is typically large which needs to identify the switch statuses of multiple time steps. Thus, the coordination mechanism between load and RDG is more complicated in the ADN optimization problem considering time variations feature. To address these challenges, we present the important theories of the network optimization problem for ADN [39]. Based on the theory in [39], a multi-period optimization model and its methods are proposed to obtain suitable optimal operation schemes for real-time ADN systems in this paper. The optimization model allows the topological network to change in each time interval of the optimization cycle. And the model can realize fast adjustment of operation scheme under real-time load and RDG output. Moreover, a dynamic optimization method is proposed to handle the multi-period optimization problem of conflicting objectives to find the optimal solution for each time interval. The key contributions of this paper are as follows:(1)
An ADN network model that is close to the actual situation of the project is established. In the network model, the maximum RDG power in any time interval is determined by combining network optimization, load demand, and RDG model.(2)
A multi-period optimization model for the ADN considering time variations in load and RDG is derived. The model aims to reduce the system power loss and improve RDG utilization by changing real-time topology at any time interval. Moreover, the overall evaluation indices of ADN operation performance are designed for a given optimization cycle.(3)
To balance the optimization objectives, the normalization and weight methods are used to establish a comprehensive objective function. Considering the dynamic characteristic of the optimization process, the normalized based values of the sub-objective functions fork time interval are the optimal value of the previous stage. Then the comprehensive objective function is derived by using the average weight method.(4)
A hybrid particle swarm optimization (HPSO) based on mixed code is used to find the optimal topology and RDG power for every time interval corresponding to conflicting objectives. Since the load demand and RDG power are time-varying, the proposed method is continuous guiding the search for the optimal solution in the next time interval according to the optimal solution of the previous period.(5)
Due to the differences in dimensions and representation between the topology and RDG power, the proposed method uses integer coding based on loop network and random number coding distribution to represent the variables. The globally optimal solution is obtained by simultaneously updating quantum position of mixed variables, which improves the search efficiency of the proposed method.The rest of this paper is organized as follows: Section II gives the time-varying of load and RDG model. Sections III presents the problem of mathematical formulation. Section IV provides the HPSO method to apply to the network optimization problem. Section V presents simulation results of two ADNs. Finally, section VI outlines conclusions.


## 2. Load and RDG model

### 2.1. Load model

In a practical scenario, the load of the ADN changes in real time because it depends on the actual demand of various type loads. Usually, the comprehensive load model of ADN is the sum of the demands of residential, commercial and industrial loads [37]. Fig 1 shows the actual hourly load for residential, commercial and industrial load patterns at a node in the ADN.Figure 1
Actual duration curves of three load patterns.For anyk time interval, the comprehensive load PLik of i node is represented as (1),(1)PLik=PRik+PIik+PCik,where PLik, PRik, PIik, and PCikare comprehensive power load, residential power load, industrial power load, and commercial power load at k time interval.


### 2.2. RDG model

#### 2.2.1. (1) wind RDG model

Wind RDG output power is often affected by some factors like weather conditions, temperature, and so on. In particular, wind speed is one of the influencing factors on the output power of wind RDG. Based on the statistical data and experience of the operators, the best expression for depicting the distribution of wind speed behavior is Weibull probability density function (PDF) [40]. The output power of wind RDG can be modeled by wind speed variable, in which the Weibull PDF is adopted to describe the distribution characteristics of wind speed. random feature and uncertainty of the wind speed. For k time interval, the actual wind RDG output behavior can be modeled as [39]:(2)Pwk=00≤vk≤vink∪voutk≤vk,Prkvk−vinkvrk−vinkvink≤vk≤vrk,Prkvrk≤vk≤voutk,where, for k time interval, Prk is the rated power of wind RDG, Pwk is the actual wind RDG output power, vk is average wind speed,vink is cut in wind speed, voutk is cut out wind speed, and vrk is rated wind speed.


#### 2.2.2. (2) PV RDG model

The PV RDG output is closely related to the factor of solar irradiance. The appropriate function to characterize solar irradiance is Beta PDF. Fork time interval, the obtained solar irradiance function by using Beta PDF is represented as follows [39]:(3)frk=Γα+βΓαΓβrkrmaxα−11−rkrmaxβ−1,where, for k time interval, Γ⋅ is the Gamma function, f(rk) is solar irradiance distribution function, rk and rmax represent actual value and maximum value of solar irradiance, α, β represent the shape variables of Gamma function.Similarly, the mathematical equation of PV RDG output is formulated as:(4)fPVwk=Γα+βΓαΓβPVwkEPα−11−PVwkEPβ−1,where PVwk is the actual output of PV RDG at k time interval, Ep is maximum power of PV RDG and its calculation equation is expressed as:(5)EP=rmaxAηpv,Where, A and ηpv are the total area of PV panels and the efficiency of PV conversion.Fork time interval, assuming that φ is a variable of power angle, the relational expression of active power PRDGk and reactive power QRDGk of RDG models is shown as follows [41].(6)QRDGk=PRDGktanφ.


### 2.3. Mathematical Formulation

#### 2.3.1. Objective Functions

According to various load demands, the comprehensive load of ADN would randomly decrease or increase in real time. And the output power of RDG changes in real time with the change of wind speed and solar radiation. The ADN network optimization is a dynamic process because of the time-varying features of load and RDG. Therefore, the optimization process in the model will be performed along the time intervals. Considering the engineering actual situation of the ADN, we seek the optimal operational scheme to enhance the operational performance of the ADN system. In this paper, the goal is to find the optimum configurations which have the minimum power losses and maximum RDG utilization simultaneously for all time intervals.(1)
Evaluation index of power lossThe active power loss on the line is an important index to evaluate the system loss. Fork time interval, the optimization problem is to find an optimal radial configuration that gives minimum active power loss. It is described as:(7)F1k=minf1k=min∑l=1NlzlkrlkPlk2+Qlk2Uik2,For a given time, the total active power lossF1total of the ADN is the sum of the power losses of all sub-time intervals.(8)F1total=∑k=1NHf1k=∑k=1NH∑l=1NlzlkrlkPlk2+Qlk2Uik2,where, for k time interval, f1k is the system power loss, l is branch label connecting node i and node j,Plk and Qlk are the active and reactive power of l branch, rlk is the l branch resistance, Uik is voltage amplitude at node i, zlk is a binary variable. If zlk represents the opening switch status, zlk=1, otherwise, zlk=0. Nl and NH represent the total number of branches and time sub-intervals.(2)
Evaluation index of RDG utilizationThe network operation mode of ADN is dynamic and variable because of time variations in load and RDG power. In the optimization method of ADN, the index of RDG utilization should be considered while improving the power quality and economy. Fork time interval, the sub-objective is to find the network configuration of ADN such that it has the maximum capability of RDG integration F2in.(9)F2ink=maxf2k=max∑h=1NGδhPRDGhk.To maintain consistency with the optimization method, the optimization objective can be transformed into (10).(10)F2k=min1f2k=min1∑h=1NGδhPRDGhk.RDG utilizationFRDG% is defined as an estimate of total absorption RDG power by the ADN during a given time.(11)FRDG%=∑k=1NH∑h=1NGδhPRDGhk∑k=1NH∑h=1NGδiPRDGhmaxk×100%,where h represents the index label of the RDG, NG is the number of RDG, δh is a binary variable representing the connection status between RDG and network, which takes δh=1 if the RDG is synchronized with the network and δh=0 otherwise. PRDGhk and PRDGhmaxk are actual output power and maximum output power of h RDG at k time interval, f2k is total output power of RDGs at k time interval.(3)
Comprehensive objective functionDue to the difference in the variation range and dimension of parameters, we should be normalized to eliminate the difference of sub-objective functions. Fork time interval, the topological structure of ADN optimization is based on the topology of the previous stage. Therefore, the normalized standard value of power loss of k time interval is obtained by combining the topological structure of k-1 time interval and load demand of k time interval. The normalized standard value of RDG output power is the sum of the upper limit RDG power in each period. The normalized standard expression of the sub-objective functions can be expressed as:(12)F1nork=∑l=1Nlzlk−1rlk−1Plk2+Qlk2Uik2.(13)F2nork=1∑h=1NGδhPRDGhmaxk,Where F1nork and F2nork represent the normalized standard values of power loss and RDG output power in k time interval, respectively.Conflict is certain existed in the objectives of RDG utilization and power loss. For example, the power loss of ADN will be very high if the penetration of RDG reaches a certain value. For coordinating the indices, we assign different significance weight coefficients to the sub-objective functions. Thus, a muti-period comprehensive optimization objective function is formulated by the weight coefficient method. The normalized comprehensive objective function fork time interval is defined as follows:(14)Fk=σ1F1kF1nork+σ2F2kF2norkwhere σ1+σ2=1 and σ1,σ2∈0,1.


#### 2.3.2. Constraints

The operating constraints include power flow constraint, voltage constraint, radial topology constraint, Branch apparent power flow constraint, and RDG constraint. These constraints are described in detail in [39], and only relevant expressions are listed here.(15)Psik+PRDGik−PLik=Uik∑j=1NzlkUjkglkcosθlk+blksinθlk,Qsik+QRDGik−QLik=Uik∑j=1NzlkUjkglkcosθlk−blksinθlk.(16)Uimink≤Uik≤Uimaxk.(17)Tgk∈TGr.(18)Slk≤Slmaxk,(19)PRDGhmink≤PRDGhk≤PRDGhmaxkh∈NG,QRDGhmink≤QRDGhk≤QRDGhmaxk.where, for k time interval, Psik and Qsik are the injected active and reactive power of the node i, PLik and QLik are the load demand of node i, PRDGik and QRDGik are the actual RDG output power on node i.glk, blk, θlk are conductance, susceptance, phase angle of branch l. Uik and Ujk are voltage amplitude of node i and j. N is the set of nodes of the ADN. Uik, Uimink, Uimaxk the real time voltage amplitude, lower limit amplitude and upper limit amplitude of node i at k time interval. Tgk is the reconfigured network topology at k time interval and TGr is the set of radial topologies. Slk, Slmaxk are apparent power flow and the apparent power limit of branch l at k time interval. PRDGhk and QRDGhk are the actual active and reactive power of i RDG at k time interval, PRDGhmink, QRDGhmink, PhRDGhmaxk, QRDGhmaxk are the minimum and maximum output power of RDGs for k time interval.


## 2.1. Load model

In a practical scenario, the load of the ADN changes in real time because it depends on the actual demand of various type loads. Usually, the comprehensive load model of ADN is the sum of the demands of residential, commercial and industrial loads [37]. Fig 1 shows the actual hourly load for residential, commercial and industrial load patterns at a node in the ADN.Figure 1
Actual duration curves of three load patterns.For anyk time interval, the comprehensive load PLik of i node is represented as (1),(1)PLik=PRik+PIik+PCik,where PLik, PRik, PIik, and PCikare comprehensive power load, residential power load, industrial power load, and commercial power load at k time interval.


## 2.2. RDG model

### 2.2.1. (1) wind RDG model

Wind RDG output power is often affected by some factors like weather conditions, temperature, and so on. In particular, wind speed is one of the influencing factors on the output power of wind RDG. Based on the statistical data and experience of the operators, the best expression for depicting the distribution of wind speed behavior is Weibull probability density function (PDF) [40]. The output power of wind RDG can be modeled by wind speed variable, in which the Weibull PDF is adopted to describe the distribution characteristics of wind speed. random feature and uncertainty of the wind speed. For k time interval, the actual wind RDG output behavior can be modeled as [39]:(2)Pwk=00≤vk≤vink∪voutk≤vk,Prkvk−vinkvrk−vinkvink≤vk≤vrk,Prkvrk≤vk≤voutk,where, for k time interval, Prk is the rated power of wind RDG, Pwk is the actual wind RDG output power, vk is average wind speed,vink is cut in wind speed, voutk is cut out wind speed, and vrk is rated wind speed.


### 2.2.2. (2) PV RDG model

The PV RDG output is closely related to the factor of solar irradiance. The appropriate function to characterize solar irradiance is Beta PDF. Fork time interval, the obtained solar irradiance function by using Beta PDF is represented as follows [39]:(3)frk=Γα+βΓαΓβrkrmaxα−11−rkrmaxβ−1,where, for k time interval, Γ⋅ is the Gamma function, f(rk) is solar irradiance distribution function, rk and rmax represent actual value and maximum value of solar irradiance, α, β represent the shape variables of Gamma function.Similarly, the mathematical equation of PV RDG output is formulated as:(4)fPVwk=Γα+βΓαΓβPVwkEPα−11−PVwkEPβ−1,where PVwk is the actual output of PV RDG at k time interval, Ep is maximum power of PV RDG and its calculation equation is expressed as:(5)EP=rmaxAηpv,Where, A and ηpv are the total area of PV panels and the efficiency of PV conversion.Fork time interval, assuming that φ is a variable of power angle, the relational expression of active power PRDGk and reactive power QRDGk of RDG models is shown as follows [41].(6)QRDGk=PRDGktanφ.


## 2.2.1. (1) wind RDG model

Wind RDG output power is often affected by some factors like weather conditions, temperature, and so on. In particular, wind speed is one of the influencing factors on the output power of wind RDG. Based on the statistical data and experience of the operators, the best expression for depicting the distribution of wind speed behavior is Weibull probability density function (PDF) [40]. The output power of wind RDG can be modeled by wind speed variable, in which the Weibull PDF is adopted to describe the distribution characteristics of wind speed. random feature and uncertainty of the wind speed. For k time interval, the actual wind RDG output behavior can be modeled as [39]:(2)Pwk=00≤vk≤vink∪voutk≤vk,Prkvk−vinkvrk−vinkvink≤vk≤vrk,Prkvrk≤vk≤voutk,where, for k time interval, Prk is the rated power of wind RDG, Pwk is the actual wind RDG output power, vk is average wind speed,vink is cut in wind speed, voutk is cut out wind speed, and vrk is rated wind speed.


## 2.2.2. (2) PV RDG model

The PV RDG output is closely related to the factor of solar irradiance. The appropriate function to characterize solar irradiance is Beta PDF. Fork time interval, the obtained solar irradiance function by using Beta PDF is represented as follows [39]:(3)frk=Γα+βΓαΓβrkrmaxα−11−rkrmaxβ−1,where, for k time interval, Γ⋅ is the Gamma function, f(rk) is solar irradiance distribution function, rk and rmax represent actual value and maximum value of solar irradiance, α, β represent the shape variables of Gamma function.Similarly, the mathematical equation of PV RDG output is formulated as:(4)fPVwk=Γα+βΓαΓβPVwkEPα−11−PVwkEPβ−1,where PVwk is the actual output of PV RDG at k time interval, Ep is maximum power of PV RDG and its calculation equation is expressed as:(5)EP=rmaxAηpv,Where, A and ηpv are the total area of PV panels and the efficiency of PV conversion.Fork time interval, assuming that φ is a variable of power angle, the relational expression of active power PRDGk and reactive power QRDGk of RDG models is shown as follows [41].(6)QRDGk=PRDGktanφ.


## 2.3. Mathematical Formulation

### 2.3.1. Objective Functions

According to various load demands, the comprehensive load of ADN would randomly decrease or increase in real time. And the output power of RDG changes in real time with the change of wind speed and solar radiation. The ADN network optimization is a dynamic process because of the time-varying features of load and RDG. Therefore, the optimization process in the model will be performed along the time intervals. Considering the engineering actual situation of the ADN, we seek the optimal operational scheme to enhance the operational performance of the ADN system. In this paper, the goal is to find the optimum configurations which have the minimum power losses and maximum RDG utilization simultaneously for all time intervals.(1)
Evaluation index of power lossThe active power loss on the line is an important index to evaluate the system loss. Fork time interval, the optimization problem is to find an optimal radial configuration that gives minimum active power loss. It is described as:(7)F1k=minf1k=min∑l=1NlzlkrlkPlk2+Qlk2Uik2,For a given time, the total active power lossF1total of the ADN is the sum of the power losses of all sub-time intervals.(8)F1total=∑k=1NHf1k=∑k=1NH∑l=1NlzlkrlkPlk2+Qlk2Uik2,where, for k time interval, f1k is the system power loss, l is branch label connecting node i and node j,Plk and Qlk are the active and reactive power of l branch, rlk is the l branch resistance, Uik is voltage amplitude at node i, zlk is a binary variable. If zlk represents the opening switch status, zlk=1, otherwise, zlk=0. Nl and NH represent the total number of branches and time sub-intervals.(2)
Evaluation index of RDG utilizationThe network operation mode of ADN is dynamic and variable because of time variations in load and RDG power. In the optimization method of ADN, the index of RDG utilization should be considered while improving the power quality and economy. Fork time interval, the sub-objective is to find the network configuration of ADN such that it has the maximum capability of RDG integration F2in.(9)F2ink=maxf2k=max∑h=1NGδhPRDGhk.To maintain consistency with the optimization method, the optimization objective can be transformed into (10).(10)F2k=min1f2k=min1∑h=1NGδhPRDGhk.RDG utilizationFRDG% is defined as an estimate of total absorption RDG power by the ADN during a given time.(11)FRDG%=∑k=1NH∑h=1NGδhPRDGhk∑k=1NH∑h=1NGδiPRDGhmaxk×100%,where h represents the index label of the RDG, NG is the number of RDG, δh is a binary variable representing the connection status between RDG and network, which takes δh=1 if the RDG is synchronized with the network and δh=0 otherwise. PRDGhk and PRDGhmaxk are actual output power and maximum output power of h RDG at k time interval, f2k is total output power of RDGs at k time interval.(3)
Comprehensive objective functionDue to the difference in the variation range and dimension of parameters, we should be normalized to eliminate the difference of sub-objective functions. Fork time interval, the topological structure of ADN optimization is based on the topology of the previous stage. Therefore, the normalized standard value of power loss of k time interval is obtained by combining the topological structure of k-1 time interval and load demand of k time interval. The normalized standard value of RDG output power is the sum of the upper limit RDG power in each period. The normalized standard expression of the sub-objective functions can be expressed as:(12)F1nork=∑l=1Nlzlk−1rlk−1Plk2+Qlk2Uik2.(13)F2nork=1∑h=1NGδhPRDGhmaxk,Where F1nork and F2nork represent the normalized standard values of power loss and RDG output power in k time interval, respectively.Conflict is certain existed in the objectives of RDG utilization and power loss. For example, the power loss of ADN will be very high if the penetration of RDG reaches a certain value. For coordinating the indices, we assign different significance weight coefficients to the sub-objective functions. Thus, a muti-period comprehensive optimization objective function is formulated by the weight coefficient method. The normalized comprehensive objective function fork time interval is defined as follows:(14)Fk=σ1F1kF1nork+σ2F2kF2norkwhere σ1+σ2=1 and σ1,σ2∈0,1.


### 2.3.2. Constraints

The operating constraints include power flow constraint, voltage constraint, radial topology constraint, Branch apparent power flow constraint, and RDG constraint. These constraints are described in detail in [39], and only relevant expressions are listed here.(15)Psik+PRDGik−PLik=Uik∑j=1NzlkUjkglkcosθlk+blksinθlk,Qsik+QRDGik−QLik=Uik∑j=1NzlkUjkglkcosθlk−blksinθlk.(16)Uimink≤Uik≤Uimaxk.(17)Tgk∈TGr.(18)Slk≤Slmaxk,(19)PRDGhmink≤PRDGhk≤PRDGhmaxkh∈NG,QRDGhmink≤QRDGhk≤QRDGhmaxk.where, for k time interval, Psik and Qsik are the injected active and reactive power of the node i, PLik and QLik are the load demand of node i, PRDGik and QRDGik are the actual RDG output power on node i.glk, blk, θlk are conductance, susceptance, phase angle of branch l. Uik and Ujk are voltage amplitude of node i and j. N is the set of nodes of the ADN. Uik, Uimink, Uimaxk the real time voltage amplitude, lower limit amplitude and upper limit amplitude of node i at k time interval. Tgk is the reconfigured network topology at k time interval and TGr is the set of radial topologies. Slk, Slmaxk are apparent power flow and the apparent power limit of branch l at k time interval. PRDGhk and QRDGhk are the actual active and reactive power of i RDG at k time interval, PRDGhmink, QRDGhmink, PhRDGhmaxk, QRDGhmaxk are the minimum and maximum output power of RDGs for k time interval.


## 2.3.1. Objective Functions

According to various load demands, the comprehensive load of ADN would randomly decrease or increase in real time. And the output power of RDG changes in real time with the change of wind speed and solar radiation. The ADN network optimization is a dynamic process because of the time-varying features of load and RDG. Therefore, the optimization process in the model will be performed along the time intervals. Considering the engineering actual situation of the ADN, we seek the optimal operational scheme to enhance the operational performance of the ADN system. In this paper, the goal is to find the optimum configurations which have the minimum power losses and maximum RDG utilization simultaneously for all time intervals.(1)
Evaluation index of power lossThe active power loss on the line is an important index to evaluate the system loss. Fork time interval, the optimization problem is to find an optimal radial configuration that gives minimum active power loss. It is described as:(7)F1k=minf1k=min∑l=1NlzlkrlkPlk2+Qlk2Uik2,For a given time, the total active power lossF1total of the ADN is the sum of the power losses of all sub-time intervals.(8)F1total=∑k=1NHf1k=∑k=1NH∑l=1NlzlkrlkPlk2+Qlk2Uik2,where, for k time interval, f1k is the system power loss, l is branch label connecting node i and node j,Plk and Qlk are the active and reactive power of l branch, rlk is the l branch resistance, Uik is voltage amplitude at node i, zlk is a binary variable. If zlk represents the opening switch status, zlk=1, otherwise, zlk=0. Nl and NH represent the total number of branches and time sub-intervals.(2)
Evaluation index of RDG utilizationThe network operation mode of ADN is dynamic and variable because of time variations in load and RDG power. In the optimization method of ADN, the index of RDG utilization should be considered while improving the power quality and economy. Fork time interval, the sub-objective is to find the network configuration of ADN such that it has the maximum capability of RDG integration F2in.(9)F2ink=maxf2k=max∑h=1NGδhPRDGhk.To maintain consistency with the optimization method, the optimization objective can be transformed into (10).(10)F2k=min1f2k=min1∑h=1NGδhPRDGhk.RDG utilizationFRDG% is defined as an estimate of total absorption RDG power by the ADN during a given time.(11)FRDG%=∑k=1NH∑h=1NGδhPRDGhk∑k=1NH∑h=1NGδiPRDGhmaxk×100%,where h represents the index label of the RDG, NG is the number of RDG, δh is a binary variable representing the connection status between RDG and network, which takes δh=1 if the RDG is synchronized with the network and δh=0 otherwise. PRDGhk and PRDGhmaxk are actual output power and maximum output power of h RDG at k time interval, f2k is total output power of RDGs at k time interval.(3)
Comprehensive objective functionDue to the difference in the variation range and dimension of parameters, we should be normalized to eliminate the difference of sub-objective functions. Fork time interval, the topological structure of ADN optimization is based on the topology of the previous stage. Therefore, the normalized standard value of power loss of k time interval is obtained by combining the topological structure of k-1 time interval and load demand of k time interval. The normalized standard value of RDG output power is the sum of the upper limit RDG power in each period. The normalized standard expression of the sub-objective functions can be expressed as:(12)F1nork=∑l=1Nlzlk−1rlk−1Plk2+Qlk2Uik2.(13)F2nork=1∑h=1NGδhPRDGhmaxk,Where F1nork and F2nork represent the normalized standard values of power loss and RDG output power in k time interval, respectively.Conflict is certain existed in the objectives of RDG utilization and power loss. For example, the power loss of ADN will be very high if the penetration of RDG reaches a certain value. For coordinating the indices, we assign different significance weight coefficients to the sub-objective functions. Thus, a muti-period comprehensive optimization objective function is formulated by the weight coefficient method. The normalized comprehensive objective function fork time interval is defined as follows:(14)Fk=σ1F1kF1nork+σ2F2kF2norkwhere σ1+σ2=1 and σ1,σ2∈0,1.


## 2.3.2. Constraints

The operating constraints include power flow constraint, voltage constraint, radial topology constraint, Branch apparent power flow constraint, and RDG constraint. These constraints are described in detail in [39], and only relevant expressions are listed here.(15)Psik+PRDGik−PLik=Uik∑j=1NzlkUjkglkcosθlk+blksinθlk,Qsik+QRDGik−QLik=Uik∑j=1NzlkUjkglkcosθlk−blksinθlk.(16)Uimink≤Uik≤Uimaxk.(17)Tgk∈TGr.(18)Slk≤Slmaxk,(19)PRDGhmink≤PRDGhk≤PRDGhmaxkh∈NG,QRDGhmink≤QRDGhk≤QRDGhmaxk.where, for k time interval, Psik and Qsik are the injected active and reactive power of the node i, PLik and QLik are the load demand of node i, PRDGik and QRDGik are the actual RDG output power on node i.glk, blk, θlk are conductance, susceptance, phase angle of branch l. Uik and Ujk are voltage amplitude of node i and j. N is the set of nodes of the ADN. Uik, Uimink, Uimaxk the real time voltage amplitude, lower limit amplitude and upper limit amplitude of node i at k time interval. Tgk is the reconfigured network topology at k time interval and TGr is the set of radial topologies. Slk, Slmaxk are apparent power flow and the apparent power limit of branch l at k time interval. PRDGhk and QRDGhk are the actual active and reactive power of i RDG at k time interval, PRDGhmink, QRDGhmink, PhRDGhmaxk, QRDGhmaxk are the minimum and maximum output power of RDGs for k time interval.


## 3. Methodology

The power flow in ADN is changing frequently with the time variations of load demand and RDG outputs. This paper focuses on the multi-period optimization model which the studied time cycle divides into several sub-time intervals. Several requirements need to be considered in dealing with the multi-period optimization problem combining RDG integration, load demand and network topology. And the HPSO method is chosen to solve the coordination optimization problem.HPSO is a meta-heuristic optimization model that mimics the swarm characteristics of the foraging behaviour of bird flocks based on the traditional PSO algorithm [42]. The movement of each particle has a random speed that determines the direction and distance of its flight. Based on the speed update rule and location update rule, the algorithm brings about local search and global information exchange during the optimization process. In the traditional PSO, the speed and position of each particle are described as random vectors. In our study, a particle is composed of the opened switch number of a possible radial network configuration and potential RDG power. The opened switch number is usually expressed as an integer code, and the RDG power can be a random code. Moreover, the method should be efficient to address the multi-time period optimization problem. The concepts are considered within the traditional PSO to develop an HPSO method discovering a set of solutions of multi-time period and multi-objective optimization model. The procedure of HPSO is shown in Fig 2. Based on the [39], the detailed implementation process of the HPSO is as follows.Figure 2
The procedure of HPSO.Step 1.
Calculating the output power of RDGs for each time interval. The maximum outputs of RDGs for sub-time intervals are obtained by measuring the real-time input parameters of RDG such as wind speed and solar irradiance. The PSO is used to determine the base valuesPr and Ep of wind and PV RDGs output power based on the peak load point. We calculate time sequence RDG outputs based on wind and PV RDG models.Step 2.
Inputting the parameters like topology structure, RDG type and location, time sequence RDG outputs and load data, and so on.Step 3.
Designing the multi-period optimization model for a given time. The problem aims to find the best solution that satisfies the requirements of power loss and RDG utilization fixed by distribution manager. Fork time interval, the integration RDG number and type are not fixed. For example, there may be no PV RDG power since the solar irradiance is 0. Therefore, the type and number of RDGs connected to ADN for k time interval rely on the real time input parameters of RDG. After that, the multi-time period optimization model is designed as equations (7)-(19).Step 4.
Initializing the population. In the proposed HPSO method, the population is composed of a set of random particles. Fork time interval, the best solution to the optimization problem is a particle in the population. A particle is defined as a vector Xd=x1x2…xd which contains the radial topology and RDG power. The network topology is usually represented by the opened switches in loops that label the integer coded. And the RDG output power is generated randomly. Given an ADN with n loops and NG RDGs, a random population is represented as equation (20).(20)Xd=x1.1x1.2⋯x1.nx2.1x2.2⋯x2.n⋮⋮⋯⋮xnop.1xnop.2⋯xnop.n︷Swiches lablex1.n+1⋯x1.n+NGx2.n+1⋯x2.n+NG⋮⋯⋮xnop.n+1⋯xnop.n+NG︷RDGpower,Step 5.
Evaluating the fitness of the comprehensive objective functions. In the proposed method, a forward-backward sweep based load flow method is used to calculate the operation indices of the ADN. The weight coefficient of the objective functions isσ1=σ2=0.5. According to equation (14), the objective function values of each particle for k time interval are calculated to identify the best and worst solutions in the entire population. Assuming that the current and global best position of a D-dimensional particle are pmdk=pm1kpm2k...pmDk and pgdk=pg1kpg2k...pgDk. For t iteration, we obtain the average and optimal best position of M particles in a population are represented as mbestk(t) and pmdkt.(21)mbestkt=1M∑m=1Mpmdkt,pmdkt=φpmdkt+1−φpgdktφ∈0,1.Step 6.
Updating the position particles satisfying all the operating constraints. Fork time interval, the updated position of topology particle ixmdnkt+1 and RDG power particle i position xidNGkt+1 after t iteration are obtained.(22)xmdnkt+1=Pmdkt±λk∗mbestkt−xmdnkt∗ln1uxmdNGkt+1=roundPmdkt±λk∗mbestkt−xmdNGkt∗ln1u,λmaxk=tmaxk−tk×βmaxktmaxk.where, for k time interval, λk and λmaxk are the current value and upper limit of control coefficient of and convergence speed, u is the uniformly distributed random number, tmaxk is the maximum number iterationsStep 7.
Settingtk=tk+1 and updating the particle positions pmdk andpgdk.Step 8.
Checking the termination conditions. The optimization process will be continued until the current iteration is less than equal to the maximum iteration number. The algorithm repeat steps (4)-(7) until these termination conditions are satisfied.Step 9.
Checking the optimization completion for all time intervalsNH. The algorithm repeats steps(2)–(8) until k≥NH is satisfied.Step 10.
Output the best optimal results of all time intervals.


## 4. Case Study And Results

### 4.1. 33-bus ADN

The multi-time period model and proposed optimization method is tested on an ADN, presented in Fig3. The system has five loops which the tie-branches are e34-e37. The line impedances can be found in [43]. From Fig 3, the wind and PV RDGs are installed on the corresponding nodes. The voltage fluctuation range is 0.90-1.05 p.u. The comprehensive objective function is shown in (9) which the weight values are set σ1=σ2=0.5.Figure 3
The 33-bus ADN.The real-time load of one day is considered an optimization cycle. Based on the load and RDG models, Fig4 show the actual hourly load demand and output power of RDGs in a certain area.Figure 4
Hourly load demand and RDG power.


#### 4.1.1. Benefits of combining network optimization and RDG integration

To analyze the operation performance of ADN after RDG integration and network optimization, three cases are simulated to utilize the proposed method at peak load point.Case 1.
Initial system without RDGs and network optimization.Case 2.
Solving the optimal problem of RDG utilization maximization and power loss minimization without network optimization.Case 3.
Solving the problem of RDG utilization maximization and power loss minimization with network optimization.
In Fig4, the peak load point of actual hourly load demand is 7.4406MW+4.5314MVar at 19:00. The obtained results of the cases are presented in Table 1.
In case 2, the RDG is optimally integrated utilization without network optimization, providing an 80.2% reduction of system power loss (Ploss). However, the Ploss is reduced by 89.53% in case 3. Moreover, the absorption RDG power (PRDG) by the system has improved to 6.1836MW after optimization. The combination of network optimization and RDG integration is effective in reducing the system power loss of the ADN.
Fig5 shows the hourly voltage profiles of all the cases. The voltage amplitude is improved when combining network optimization with RDG integration than when integrating RDG only and base system. The minimum voltage (Vmin) amplitude of case 3 has increased to 0.9734p.u., which is above 0.9494p.u. of case 2 and 0.8104 p.u. of case 1. The results show that system efficiency in the optimization problem with adding the RDG is better than the base system and has resulted in more power loss reduction, more voltage profile change, and more RDG power utilization. The use of RDG along with network optimization leads to better operation conditions of the ADN due to the injection of RDG power.
To justify the superiority of the proposed method, the optimization problem in case 3 is also solved by the conventional algorithms, GA method [43] and moth-flame optimization (MFO) [44]. the obtained simulation results are compared with those obtained by the conventional methods. The methods are based on initial conditions and similar assumptions. Table 2 and Fig 6 illustrate the detailed comparison results.
The proposed method presents a superior performance to the GA and MFO. The proposed method provides an optimal solution with 99.9680kW of system power loss which is less than 101.3932kW in GA algorithm. According to Fig6, the proposed method could achieve optimal solution faster than other methods, i. e., 47 iterations against 95 in MFO, 47 iterations compared to 157 in GA. Consequently, the proposed method is an effective way for addressing ADN network optimization due to its high convergence ability and accuracy.Table 1
Benefits of combining network optimization and RDG integration.
CasesOpen switchesPloss/kWPRDG/MWVmin/p.u.1e33-e37954.864400.81042e33-e37189.03625.94230.94943e33, e3, e17, e35, e1399.96806.18360.9734Figure 5
The voltage profile for cases.Table 2
Comparison with last optimization methods.
ItemsOpen switchesPloss/kWPRDG/MWVmin/p.u.Proposed methode33, e3, e17, e35, e1399.96806.18360.9734GA [43]e33, e28, e17, e35, e13101.39326.18360.9662MFO [44]e33, e3, e17, e35, e1399.96806.18360.9734Fig 6
The convergence curve of the IEEE33 bus system.


#### 4.1.2. Experiment considering time variations in RDG and load

The proposed model and method is used to address the optimization problem of ADN considering the time variations in RDG and load. The hourly data in RDG and load of an average day as an example of study are illustrated in Fig4. In the 33-bus ADN, four cases are adopted to investigate the effect of time-varying characteristics on network optimization indicators.Case 4.
Load variation is considered without network optimization and RDG.Case 5.
Load variation and network optimization simulation are considered.Case 6.
Time variations in load and RDG are considered without network optimization.Case 7.
Simultaneous hourly network optimization and time variations in load and RDG are studied.
Using the proposed model and method, we obtained the selected open switches and RDG utilization for a given time certain. Tables3 and 4 present the simulation results and optimal solutions for each time interval for the cases. In each time interval, the open switches, Ploss, Vmin, and absorbed RDGs power are presented corresponding to each time interval for cases.
In cases 4 and 5, the simulation is carried out in the system with time variation of load demand. Table3 lists the results of active power loss, minimum node voltage, and open switches before and after optimization (case 4 and case 5). As seen in Table 3, the optimal solutions with opening switches e7, e14, e9, e32, and e28 are consistent at time intervals 0∼4th, 8∼11th, 13∼15th, 21th, and 23th. The opened switches in other time intervals are e7, e14, e9, e32, and e37 after optimization. The active power loss for each time interval reduces after optimization, i.e, 41.4271kW against 59.2882kW at 6:00, 458.5657kW compared to 707.9717kW at 10:00, 458.6379kW instead of 631.8764kW at 16:00, and so on. Further, it is noted that the minimum node voltage for each interval has increased in case 5 compared in case 4. For example, The percentage improvement in minimum node voltage value after optimization is 1.67%, 5.08%, and 6.40% at 6:00, 10:00, and 16:00, respectively. Fig 7 gives the total power loss and minimum voltage profile in a day of the cases. From Fig 7, the total daily power loss in case 4 is 10.0352MW. The daily total power loss after optimization has been reduced by 34.69%. The minimum voltage amplitude in a day improves to 0.8646p.u. in case 5, which implies an improvement of 6.69% compared to case 4. The results illustrate network optimization can effectively improve the system operation level. However, the daily minimum voltage amplitude 0.8646p.u. in case 5 is lower than the voltage limit of 0.90p.u. It is not conducive to the long-term stable operation of the active distribution system. The way of integration of RDG is considered to improve operating indices.
To perform the suggested method which requires time variations in RDG and load demand, different load profiles and RDG output power are added to the system at each time interval. Table4 depicts the variations of the operational indices of the system in cases 6 and 7. It is observed that the optimization effects of cases 6 and 7 are better than those of cases 4 and 5. In cases 6 and 7, the obtained optimal solutions provide the total daily power loss are 2.6758MW and 1.2290MW, respectively. The maximum power loss of the day in case 7 is 194.9442kw at 15:00, which implies a reduction of 762.7595kW, 263.0196kW, and 364.0683kW compared to case 4, case 5, and case 6. From Fig 7, the daily minimum voltage amplitude of case 7 is 0.9463p.u., which is higher than the minimum voltage limit 0.90p.u. It indicates that simultaneous network optimization with time-varying load and RDG can improve the performance of the ADN. Moreover, the day absorption RDG power for case 7 is 90.2614MW, which has increased by 12.42% with case 6. And the RDG utilization (RDG%) has been improved from 67.04% in case 6 to 75.37% in case 7.
Due to space reasons, the operational indices of the ADN at 11:00 and 21:00 are analyzed as examples to describe the operation characteristic at each time interval clearly. Fig8 shows the power loss of each branch at 11:00 and 21:00. In case 4, the power loss of ADN is relatively large without network optimization and integration of RDGs. In case 5, it is significantly reduced after optimization by changing the topology with opening switch combinations [e7, e14, e9, e32, e28]. The system power loss is reduced by 53.62% at 11:00, and 32.87% at 21:00, respectively. After integrating of RDGs, The system power loss for case 6 is reduced by 78.27% at 11:00, and 77.95% at 21:00, respectively. The benefits of simultaneously combining network optimization with the integration of RDG are investigated in case 7. The obtained final opening switch combinations at 11:00 and 21:00 are [e7, e28, e8, e10, e14] and [e33, e28, e17, e35, e13], respectively. It provides a reduction to 63.2703kW at 11:00 and 28.3297kW at 21:00 of system power loss, which is significantly lower than that of other cases.
Fig9 shows the voltage profile at two-time interval points. In case 4, voltage violations occur in most of the nodes which leads to an unsatisfactory low voltage profile. The network performance is improved by using network optimization and integration of RDGs simultaneously. It is observed that the system voltage volatility range is reduced in case 7. And the minimum voltage profile in cases 4-7 are 0.8449 p.u. for node 33 at 11:00, 0.8943 p.u. for node 32 at 11:00, 0.9539 p.u. for node 33 at 11:00, 0.9912 p.u. for node 24 at 11:00. The minimum voltage of ADN has raised dramatically with the network optimization and RDG integration.
As seen in Table4, the absorption RDG power at two-time intervals by the system of case 6 are 5.2497 MW and 3.4541MW, while the optimization system of case 7 reaches 6.3646MW and 3.5698MW. It turns out that the absorption power of RDGs in case 6 is increased by 21.24% at 11:00 and 3.4% at 21:00 in case 6. The RDG utilization of a day in case 7 has increased by 15.48% and 3.24% compared to case 6. Moreover, Fig 10 compares the absorption RDG power over 24 hours for cases 6 and 7. The minimum RDG utilization in case 7 is 19.68% at 17:00, which is higher 16.64% in case 6. In both cases, the maximum difference of RDG utilization has reached to 56.94%. Thus, considering the time variations in load and RDG, the operational performance of the ADN would improve by using the proposed dynamic optimization model and the HPSO method.Table 3
The simulation results for cases 4 and 5.
HourCase 4Case 5Ploss/kWVmin/p.u.Open switchesPloss/kWVmin/p.u.Open switches022.71690.9708e33, e34, e35, e36, e3715.42030.9800e7, e14, e9, e32, e28122.71680.970815.42030.9800216.30750.974610.45090.983636.05090.98494.09180.990046.05070.98494.09180.9900526.48390.967818.62090.9774e7, e14, e9, e32, e37659.28820.952141.42710.96757102.65280.938171.65230.95578454.68530.8711299.25720.9111e7, e14, e9, e32, e289634.78270.8465406.95870.895810707.97170.8399458.56570.890711667.97350.8449434.82380.894312450.93750.8707303.20610.9095e7, e14, e9, e32, e3713524.29160.8622353.28660.9055e7, e14, e9, e32, e2814655.23440.8467429.30500.894515957.70370.8154617.22440.873516972.44770.8126631.87640.8646e7, e14, e9, e32, e3717800.69860.8283523.44940.878618704.07660.8373460.81270.887019954.86440.8104612.00810.869220892.39200.8166574.12480.872621245.25010.9043164.64890.9364e7, e14, e9, e32, e2822107.07330.935974.21700.9554e7, e14, e9, e32, e372342.59500.960329.01500.9729e7, e14, e9, e32, e28Table 4
the simulation results for cases 6 and 7.
HourCase 6Case 7Ploss/kWVmin/puPRDG/MWOpen switchesPloss/kWVmin/puPRDG/MWOpen switches08.10670.99751.4218e33, e34, e35, e36, e376.03941.00001.6812e18, e26, e15, e35, e1318.10670.99751.42176.03941.00001.6812e18, e26, e15, e35, e1325.84690.99801.25324.84681.00001.4586e18, e26, e15, e35, e1333.29580.99890.84662.54421.00000.9858e18, e26, e15, e35, e1343.29580.99620.84662.54421.00000.9858e18, e26, e15, e35, e13515.40350.99261.504714.50191.00001.7885e18, e26, e15, e35, e13616.59980.99022.121210.66401.00002.5000e18, e26, e15, e35, e13720.76170.97912.860011.78041.00003.1844e33, e26, e34, e11, e13867.94960.90945.158332.08320.99615.5973e33, e27, e34, e10, e149260.68590.92503.739493.82670.96173.7394e20, e28, e7, e9, e1410255.55330.95394.905976.62670.97855.0553e19, e5, e34, e8, e1411145.13870.95395.249763.27030.99126.3646e7, e28, e8, e10, e1412135.43370.96034.255857.90870.98985.2516e3, e28, e8, e9, e1313123.24610.96254.732956.20660.99125.7489e7, e28, e8, e9, e1414126.85130.93735.364860.62880.99236.4603e7, e28, e8, e10, e1415559.01250.85131.5438194.94420.94633.5852e6, e37, e15, e8, e 1416174.54500.94715.485880.62160.97065.4860e20, e5, e15, e10, e 1417121.46350.96435.048072.54540.97065.0480e33, e4, e17, e11, e1418165.78090.92103.3588126.01270.94923.3588e33, e28, e32, e10, e1419194.79160.95166.071199.96800.97346.1836e33, e3, e17, e35, e1320166.38800.93825.163697.50130.96285.1636e33, e4, e17, e11, e132154.06690.97683.454128.32970.98593.5698e33, e28, e17, e35, e132230.28930.99082.655420.01710.99593.2057e18, e26, e15, e35, e132313.17780.99581.82979.57390.99942.1778e18, e26, e15, e35, e13Figure 7
the simulation results for all the cases in a day.Figure 8
Power loss of each branch at 11:00 and 21:00.Figure 9
Voltage profile of each node at 11:00 and 21:00.Fig 10
RDG utilization at of each time interval in case 6 and case 7.


#### 4.1.3. Comparing the results of HPSO and Other Methods

As a comparison, the optimization problem of ADN considering time variation feature is also addressed by the methods of BPSO [19], BAS (beetle antennae search algorithm) [37], GA [43] and MFO [44].Using all the methods, the list of final open switches for the best solutions during 24 hours are listed in Table5. The crucial constraint for the radial structure is satisfied for each topology. According to the optimal topologies, Fig 11 shows the power loss, minimum voltage standard amplitude, and RDG utilization based on case 7 using HPSO, BPSO, BAS, GA and MFO for the ADN. The obtained results show that the indexes of the proposed HPSO are better than other methods for ADN optimization considering time variations in load and RDG over 24 hours.Table 5
The open switches of four algorithms at each time interval.
Houropen switchesHPSOBPSOBASGAMFO0e18, e26, e15, e35, e13e7, e14, e9, e32, e28e3, e14, e11, e32, e28e18, e26, e15, e35, e13e18, e26, e15, e35, e131e18, e26, e15, e35, e13e7, e14, e9, e32, e28e33, e14, e11, e32, e3e18, e26, e15, e35, e13e18, e26, e15, e35, e132e18, e26, e15, e35, e13e7, e14, e9, e32, e28e33, e14, e10, e36, e3e18, e26, e16, e11, e12e18, e25, e16, e11, e123e18, e26, e15, e35, e13e7, e14, e9, e32, e28e33, e13, e10, e17, e3e18, e26, e16, e35, e13e18, e25, e16, e35, e134e18, e26, e15, e35, e13e7, e14, e9, e32, e28e33, e14, e11, e32, e4e18, e26, e16, e35, e13e18, e26, e16, e11, e125e18, e26, e15, e35, e13e7, e14, e9, e32, e28e33, e13, e11, e17, e3e18, e27, e16, e35, e13e18, e26, e16, e35, e136e18, e26, e15, e35, e13e7, e14, e9, e32, e37e33, e14, e11, e32, e3e19, e26, e15, e21, e13e19, e26, e15, e21, e137e33, e26, e34, e11, e13e7, e14, e9, e32, e28e33, e14, e11, e16, e4e33, e25, e15, e11, e13e33, e26, e15, e11, e138e33, e27, e34, e10, e14e7, e14, e9, e32, e28e33, e11, e21, e15, e4e33, e27, e34, e10, e14e33, e27, e34, e9, e139e20, e28, e7, e9, e14e7, e14, e9, e32, e28e6, e34, e9, e14, e28e20, e28, e7, e9, e14e20, e28, e6, e9, e1410e19, e5, e34, e8, e14e7, e14, e9, e32, e28e6, e14, e10, e15, e28e19, e5, e34, e8, e14e19, e 5, e8, e9, e1411e7, e28, e8, e10, e14e7, e14, e9, e32, e28e5, e14, e10, e15, e28e7, e28, e34, e8, e14e7, e14, e9, e32, e2812e3, e28, e8, e9, e13e7, e14, e9, e32, e37e33, e9, e6, e14, e28e7, e4, e8, e9, e14e7, e4, e8, e10, e1413e7, e28, e8, e9, e14e7, e14, e9, e32, e37e5, e13, e9, e15, e28e7, e28, e8, e10, e14e7, e28, e8, e10, e1414e7, e28, e8, e10, e14e7, e14, e9, e32, e28e33, e34, e9, e14, e5e7, e28, e8, e10, e14e7, e28, e8, e10, e1415e6, e37, e15, e8, e 14e7, e14, e9, e32, e28e6, e34, e9, e14, e37e6, e37, e15, e9, e 34e6, e37, e15, e9, e 3416e20, e5, e15, e10, e 14e7, e14, e9, e32, e37e5, e13, e9 15, e28e20, e5, e15, e9, e 14e20, e5, e15, e10, e 1417e33, e4, e17, e11, e14e7, e14, e9, e32, e28e33, e14, e11, e17, e4e33, e4, e32, e11, e14e33, e4, e17, e11, e1418e33, e28, e32, e10, e14e7, e14, e9, e32, e37e33, e11, e6, e14, e28e33, e28, e32, e10, e14e33, e28, e32, e10, e1419e33, e3, e17, e35, e13e7, e14, e9, e32, e37e33, e13, e11, e36 , e4e33, e3, e16, e35, e13e33, e3, e17, e35, e1320e33, e4, e17, e11, e13e7, e14, e9, e32, e28e33, e14, e11, e36, e4e33, e4, e17, e11, e13e33, e4, e17, e11, e1321e33, e28, e17, e35, e13e7, e14, e9, e32, e28e33, e14, e11, e36, e3e33, e28, e17, e35, e13e33, e28, e17, e35, e1322e18, e26, e15, e35, e13e7, e14, e9, e32, e37e33, e13, e10, e36, e3e18, e26, e15, e35, e13e18, e26, e15, e35, e1323e18, e26, e15, e35, e13e7, e14, e9, e32, e28e33, e14, e11, e36, e3e18, e2, e15, e35, e12e18, e26, e15, e 35, e13Table 6
Obtained results.
ItemsTotal Ploss/MWhRDG%Vmin/puHPSO1.229075.37%0.9463GA1.233575.22%0.9461MFO1.230275.28%0.9463BPSO2.303150.98%0.8848BAS1.344654.53%0.9409Figure 11
Comparison of power loss at each time interval.The power loss of the 33-bus ADN after optimization during 24 hours is shown graphically in Fig11. The maximum total power loss for the methods over 24 hours is 194.9442 kWh (HPSO), 389.0171 kWh (BPSO), 227.4633 kWh (BAS), 195.7518 kWh (GA), and 195.2001 kWh (MFO) at 15 p.m.Fig12 shows the RDG utilization after optimization at each time interval of a day. The results show that the minimum RDG utilization for each algorithm over 24 hours is 19.6811%, 9.0749%, 11.0089%, 19.6928%, and 19.6811%. The maximum RDG utilization for each algorithm over 24 hours is 100%, 92.9588%, 99.9970%, 99.9940%, and 100%. Moreover, the RDG utilization has reached 100% in eight-time intervals, which is obviously better than other methods.Figure 12
Comparison of RDG utilization at each time interval.The minimum voltage profiles of the ADN after optimization in a day are illustrated in Fig13. The day minimum voltage amplitude of each algorithm is 0.9463 p.u., 0.8848 p.u., 0.9409p.u., 0.9461 p.u., and 0.9463 p.u. at 15 p.m. The minimum node voltage for a day is almost the same as that of MFO, but 0.21%, 0.57%, and 6.95% higher than GA, BAS, and BPSO.Figure 13
Comparison of minimum voltage at each time interval.Table7 shows the total active power loss and RDG utilization in a day. The day power loss by the proposed HPSO is 1.2290MWh which is lower than 1.2302MWh (GA), 1.2335MWh (MFO), 1.3346MWh (BAS), and 2.3031MWh (BPSO). The RDG utilization for a day using HPSO, GA, MFO, BAS, and BPSO are 75.37%, 75.22%, 75.28%, 54.53%, and 50.98%, respectively. From all the results, we can see that the obtained results of HPSO algorithm are better to compare the GA, MFO, BAS, and BPSO. Thus, the proposed HPSO method is better in global optimization ability.Table 7
The simulation results of 69-bus ADN in a day.
ItemsTotal Ploss/MWTotal RDG/MWVmin/puCase 494.440400.4092Case 51.875200.8079Case 6180.967483.51320.5476Case 71.4379123.61390.9203


### 4.2. 69-bus ADN

The 69-bus ADN consists of 69 buses and 5 tie branches, as provided in Fig14 [44]. The wind RDGs are connected to nodes 11 and 50. The PV RDGs are connected to nodes 38 and 53. The open switches in the initial network are e69, e70, e71, e72, and e73. The proposed method is applied to obtain the optimal reconfigure to enable RDG penetration. The hourly load demand for the system is obtained by using the method in [43].Figure 14
69-bus ADN.For the sake of investigating the optimal problem of the ADN considering time variations in load and RDG, cases 4-7 will be implemented to the 69-bus ADN. Table7 provides simulated results of the test active distribution network for 24 hours using the proposed method.The original network considers the time-varying nature of load demand without the integration of RDGs and network optimization. From Table7, we observe that the total daily active power loss of 94.4404MW for case 4. After network optimization, the total power loss has been reduced to 1.8752MW in case 5. To discuss the RDG effect, consider the integration of RDGs as case 6. The total power loss of the case is calculated as 169.0020MW. When both network optimization and RDG integration are considered, the daily power loss is minimized to 1.4379MW. Fig 15 shows active power loss for the four cases during 24 hours. In all hours, the power loss has been reduced greatly using network optimization and integration of RDGs in case 7.Figure 15
Power loss of 69-bus ADN during 24 hours.The voltage profiles of the cases during 24 hours have been illustrated in Fig16. The hourly minimum voltages in cases are 0.4092p.u. at 15:00, 0.8079 p.u. at 21:00, 0.5476 p.u. at 8:00, and 0.9203 p.u. at 14:00, respectively. It is observed that the minimum voltage profile for cases 4-6 is lower than the lower voltage limit 0.9p.u. The network optimization based on time variation of RDGs and load demand brings much less voltage fluctuations with respect to case 7 than case 5 and case 6 because the scheme considers the network optimization and integration of RDGs.Figure 16
Voltage profiles of 69-bus ADN during 24 hours.Fig17 also compares the absorbed RDG power by the system at each interval before and after network optimization (case 6 and case 7). We observe that the total absorbed RDG power after optimization for a day is 123.6139MW, which has increased by 45.40% before the system. The RDG utilization has been improved from 52.93% to 76.96% after optimization. With the expansion of network structure, the proposed method is more effective for improving the operation indices of the ADN.Figure 17
Absorbed RDG power of 69-bus ADN during 24 hours.Table8 shows the open switches of the hourly optimization solution for the cases. The application of network optimization to the ADN has permitted to achieve a better operational conditions. As an example, the topological structure [e7, e18, e13, e50, e47] is optimal for the hour 8:00 of the day in case 7. It ensures a reduction of 93.37%, 49.97%, and 95.32% of power loss with other cases. In addition, The voltage profile is increased by 42.87%, 8.68% and 66.58% for case 4, case 5, and case 6, respectively. And the absorbed RDG power of case 7 is improved to 6.6721MW which is higher than 4.8958MW of case 6. These results show that the proposed model and method yield acceptable results in terms of power loss reduction, voltage improvement, and RDG utilization for both test systems under time variations of load demand and RDG. This demonstrates the ability of the proposed approach to find the optimal network topology along with absorbed RDG power and load demand for the ADN considering time characteristics.Table 8
Open switches of hourly optimization for 69-bus ADN.
HourCase 4Case 5Case 6Case 70e69, e70, e71, e72, e73e64, e14, e68, e44, e38e69, e70, e71, e72, e73e4, e18, e13, e53, e471e3, e18, e67, e50, e42e4, e17, e14, e26, e472e3, e18, e12, e47, e7e4, e18, e13, e53, e473e69, e19, e62, e50, e47e4, e18, e13, e53, e474e4, e16, e14, e53, e38e4, e18, e13, e25, e475e4, e17, e68, e52, e38e4, e18, e13, e72, e476e65, e17, e71, e52, e36e4, e18, e13, e24, e477e10, e15, e13, e12, e35e4, e16, e13, e20, e418e9, e15, e68, e20, e37e7, e18, e13, e50, e479e3, e16, e11, e13, e36e3, e15, e5, e12, e4510e3, e17, e 68, e11, e4e3, e17, e9, e12, e711e8, e70, e68, e11, e36e69, e13, e 67, e48, e712e9, e70, e12, e16 , e36e61, e16, e7, e12, e4113e9, e15, e13, e20, e36e10, e15, e68, e11, e514e3 e70, e 66, e11, e7e3, e15, e71, e50, e4115e10, e15, e13, e23, e37e69, e13, e5, e49, e3516e8, e16, e13, e11, e36e65, e16, e14, e20, e4717e8, e15, e71, e11, e35e3, e17, e68, e20, e718e8, e17, e11, e13, e36e60, e15, e9, e11, e3619e8, e15, e66, e12, e38e3, e19, e68, e21, e520e9, e15, e68, e22, e36e65, e16, e68, e53, e721e9, e18, e12, e14, e35e4, e18, e68, e20, e3622e10, e15, e68, e25, e36e4, e18, e13, e72, e4723e3, e19, e11, e26, e36e4, e15, e68, e25, e47To further illustrate the superiority of the HPSO, we compared and analyzed the 24 hours simulated results of case 7 by methods of HPSO, BPSO, BAS, GA, and MFO. Fig18–20 show the comparison of the proposed HPSO and the conventional approaches BPSO, BAS, GA, and MFO. And the comparison of the simulation results in a day is also shown in Fig 21–22. We can clearly the minimum voltage profiles in a day are 0.9203 p.u. for HPSO, 0.8487 p.u. for BPSO, 0.9048 p.u. for BAS, 0.9009 p.u. for GA, 0.9066 p.u. for MFO, respectively. From Fig 18, the total active power loss of a day 1.4379 MWh by HPSO, which is less 77.57%, 67.07%, 34.34%, and 21.09% than BPSO, BAS, GA, and MFO. The day RDG utilization with a penetration rate of 76.96% is the highest among all the methods. Thus, considering the overall operation economy in terms of power loss, voltage quality, and RDG utilization, the proposed HPSO is superior to the other methods.Figure 18
Comparison of power loss of 69-bus ADN.Figure19
Comparison of RDG utilization of 69-bus ADN.Figure 20
Comparison of minimum voltage of 69-bus ADN.Figure 21
Comparison total power loss and minimum voltage profile for all the methods in a day.Figure 22
Comparison RDG power and RDG utilization for all the methods in a day.


## 4.1. 33-bus ADN

The multi-time period model and proposed optimization method is tested on an ADN, presented in Fig3. The system has five loops which the tie-branches are e34-e37. The line impedances can be found in [43]. From Fig 3, the wind and PV RDGs are installed on the corresponding nodes. The voltage fluctuation range is 0.90-1.05 p.u. The comprehensive objective function is shown in (9) which the weight values are set σ1=σ2=0.5.Figure 3
The 33-bus ADN.The real-time load of one day is considered an optimization cycle. Based on the load and RDG models, Fig4 show the actual hourly load demand and output power of RDGs in a certain area.Figure 4
Hourly load demand and RDG power.


### 4.1.1. Benefits of combining network optimization and RDG integration

To analyze the operation performance of ADN after RDG integration and network optimization, three cases are simulated to utilize the proposed method at peak load point.Case 1.
Initial system without RDGs and network optimization.Case 2.
Solving the optimal problem of RDG utilization maximization and power loss minimization without network optimization.Case 3.
Solving the problem of RDG utilization maximization and power loss minimization with network optimization.
In Fig4, the peak load point of actual hourly load demand is 7.4406MW+4.5314MVar at 19:00. The obtained results of the cases are presented in Table 1.
In case 2, the RDG is optimally integrated utilization without network optimization, providing an 80.2% reduction of system power loss (Ploss). However, the Ploss is reduced by 89.53% in case 3. Moreover, the absorption RDG power (PRDG) by the system has improved to 6.1836MW after optimization. The combination of network optimization and RDG integration is effective in reducing the system power loss of the ADN.
Fig5 shows the hourly voltage profiles of all the cases. The voltage amplitude is improved when combining network optimization with RDG integration than when integrating RDG only and base system. The minimum voltage (Vmin) amplitude of case 3 has increased to 0.9734p.u., which is above 0.9494p.u. of case 2 and 0.8104 p.u. of case 1. The results show that system efficiency in the optimization problem with adding the RDG is better than the base system and has resulted in more power loss reduction, more voltage profile change, and more RDG power utilization. The use of RDG along with network optimization leads to better operation conditions of the ADN due to the injection of RDG power.
To justify the superiority of the proposed method, the optimization problem in case 3 is also solved by the conventional algorithms, GA method [43] and moth-flame optimization (MFO) [44]. the obtained simulation results are compared with those obtained by the conventional methods. The methods are based on initial conditions and similar assumptions. Table 2 and Fig 6 illustrate the detailed comparison results.
The proposed method presents a superior performance to the GA and MFO. The proposed method provides an optimal solution with 99.9680kW of system power loss which is less than 101.3932kW in GA algorithm. According to Fig6, the proposed method could achieve optimal solution faster than other methods, i. e., 47 iterations against 95 in MFO, 47 iterations compared to 157 in GA. Consequently, the proposed method is an effective way for addressing ADN network optimization due to its high convergence ability and accuracy.Table 1
Benefits of combining network optimization and RDG integration.
CasesOpen switchesPloss/kWPRDG/MWVmin/p.u.1e33-e37954.864400.81042e33-e37189.03625.94230.94943e33, e3, e17, e35, e1399.96806.18360.9734Figure 5
The voltage profile for cases.Table 2
Comparison with last optimization methods.
ItemsOpen switchesPloss/kWPRDG/MWVmin/p.u.Proposed methode33, e3, e17, e35, e1399.96806.18360.9734GA [43]e33, e28, e17, e35, e13101.39326.18360.9662MFO [44]e33, e3, e17, e35, e1399.96806.18360.9734Fig 6
The convergence curve of the IEEE33 bus system.


### 4.1.2. Experiment considering time variations in RDG and load

The proposed model and method is used to address the optimization problem of ADN considering the time variations in RDG and load. The hourly data in RDG and load of an average day as an example of study are illustrated in Fig4. In the 33-bus ADN, four cases are adopted to investigate the effect of time-varying characteristics on network optimization indicators.Case 4.
Load variation is considered without network optimization and RDG.Case 5.
Load variation and network optimization simulation are considered.Case 6.
Time variations in load and RDG are considered without network optimization.Case 7.
Simultaneous hourly network optimization and time variations in load and RDG are studied.
Using the proposed model and method, we obtained the selected open switches and RDG utilization for a given time certain. Tables3 and 4 present the simulation results and optimal solutions for each time interval for the cases. In each time interval, the open switches, Ploss, Vmin, and absorbed RDGs power are presented corresponding to each time interval for cases.
In cases 4 and 5, the simulation is carried out in the system with time variation of load demand. Table3 lists the results of active power loss, minimum node voltage, and open switches before and after optimization (case 4 and case 5). As seen in Table 3, the optimal solutions with opening switches e7, e14, e9, e32, and e28 are consistent at time intervals 0∼4th, 8∼11th, 13∼15th, 21th, and 23th. The opened switches in other time intervals are e7, e14, e9, e32, and e37 after optimization. The active power loss for each time interval reduces after optimization, i.e, 41.4271kW against 59.2882kW at 6:00, 458.5657kW compared to 707.9717kW at 10:00, 458.6379kW instead of 631.8764kW at 16:00, and so on. Further, it is noted that the minimum node voltage for each interval has increased in case 5 compared in case 4. For example, The percentage improvement in minimum node voltage value after optimization is 1.67%, 5.08%, and 6.40% at 6:00, 10:00, and 16:00, respectively. Fig 7 gives the total power loss and minimum voltage profile in a day of the cases. From Fig 7, the total daily power loss in case 4 is 10.0352MW. The daily total power loss after optimization has been reduced by 34.69%. The minimum voltage amplitude in a day improves to 0.8646p.u. in case 5, which implies an improvement of 6.69% compared to case 4. The results illustrate network optimization can effectively improve the system operation level. However, the daily minimum voltage amplitude 0.8646p.u. in case 5 is lower than the voltage limit of 0.90p.u. It is not conducive to the long-term stable operation of the active distribution system. The way of integration of RDG is considered to improve operating indices.
To perform the suggested method which requires time variations in RDG and load demand, different load profiles and RDG output power are added to the system at each time interval. Table4 depicts the variations of the operational indices of the system in cases 6 and 7. It is observed that the optimization effects of cases 6 and 7 are better than those of cases 4 and 5. In cases 6 and 7, the obtained optimal solutions provide the total daily power loss are 2.6758MW and 1.2290MW, respectively. The maximum power loss of the day in case 7 is 194.9442kw at 15:00, which implies a reduction of 762.7595kW, 263.0196kW, and 364.0683kW compared to case 4, case 5, and case 6. From Fig 7, the daily minimum voltage amplitude of case 7 is 0.9463p.u., which is higher than the minimum voltage limit 0.90p.u. It indicates that simultaneous network optimization with time-varying load and RDG can improve the performance of the ADN. Moreover, the day absorption RDG power for case 7 is 90.2614MW, which has increased by 12.42% with case 6. And the RDG utilization (RDG%) has been improved from 67.04% in case 6 to 75.37% in case 7.
Due to space reasons, the operational indices of the ADN at 11:00 and 21:00 are analyzed as examples to describe the operation characteristic at each time interval clearly. Fig8 shows the power loss of each branch at 11:00 and 21:00. In case 4, the power loss of ADN is relatively large without network optimization and integration of RDGs. In case 5, it is significantly reduced after optimization by changing the topology with opening switch combinations [e7, e14, e9, e32, e28]. The system power loss is reduced by 53.62% at 11:00, and 32.87% at 21:00, respectively. After integrating of RDGs, The system power loss for case 6 is reduced by 78.27% at 11:00, and 77.95% at 21:00, respectively. The benefits of simultaneously combining network optimization with the integration of RDG are investigated in case 7. The obtained final opening switch combinations at 11:00 and 21:00 are [e7, e28, e8, e10, e14] and [e33, e28, e17, e35, e13], respectively. It provides a reduction to 63.2703kW at 11:00 and 28.3297kW at 21:00 of system power loss, which is significantly lower than that of other cases.
Fig9 shows the voltage profile at two-time interval points. In case 4, voltage violations occur in most of the nodes which leads to an unsatisfactory low voltage profile. The network performance is improved by using network optimization and integration of RDGs simultaneously. It is observed that the system voltage volatility range is reduced in case 7. And the minimum voltage profile in cases 4-7 are 0.8449 p.u. for node 33 at 11:00, 0.8943 p.u. for node 32 at 11:00, 0.9539 p.u. for node 33 at 11:00, 0.9912 p.u. for node 24 at 11:00. The minimum voltage of ADN has raised dramatically with the network optimization and RDG integration.
As seen in Table4, the absorption RDG power at two-time intervals by the system of case 6 are 5.2497 MW and 3.4541MW, while the optimization system of case 7 reaches 6.3646MW and 3.5698MW. It turns out that the absorption power of RDGs in case 6 is increased by 21.24% at 11:00 and 3.4% at 21:00 in case 6. The RDG utilization of a day in case 7 has increased by 15.48% and 3.24% compared to case 6. Moreover, Fig 10 compares the absorption RDG power over 24 hours for cases 6 and 7. The minimum RDG utilization in case 7 is 19.68% at 17:00, which is higher 16.64% in case 6. In both cases, the maximum difference of RDG utilization has reached to 56.94%. Thus, considering the time variations in load and RDG, the operational performance of the ADN would improve by using the proposed dynamic optimization model and the HPSO method.Table 3
The simulation results for cases 4 and 5.
HourCase 4Case 5Ploss/kWVmin/p.u.Open switchesPloss/kWVmin/p.u.Open switches022.71690.9708e33, e34, e35, e36, e3715.42030.9800e7, e14, e9, e32, e28122.71680.970815.42030.9800216.30750.974610.45090.983636.05090.98494.09180.990046.05070.98494.09180.9900526.48390.967818.62090.9774e7, e14, e9, e32, e37659.28820.952141.42710.96757102.65280.938171.65230.95578454.68530.8711299.25720.9111e7, e14, e9, e32, e289634.78270.8465406.95870.895810707.97170.8399458.56570.890711667.97350.8449434.82380.894312450.93750.8707303.20610.9095e7, e14, e9, e32, e3713524.29160.8622353.28660.9055e7, e14, e9, e32, e2814655.23440.8467429.30500.894515957.70370.8154617.22440.873516972.44770.8126631.87640.8646e7, e14, e9, e32, e3717800.69860.8283523.44940.878618704.07660.8373460.81270.887019954.86440.8104612.00810.869220892.39200.8166574.12480.872621245.25010.9043164.64890.9364e7, e14, e9, e32, e2822107.07330.935974.21700.9554e7, e14, e9, e32, e372342.59500.960329.01500.9729e7, e14, e9, e32, e28Table 4
the simulation results for cases 6 and 7.
HourCase 6Case 7Ploss/kWVmin/puPRDG/MWOpen switchesPloss/kWVmin/puPRDG/MWOpen switches08.10670.99751.4218e33, e34, e35, e36, e376.03941.00001.6812e18, e26, e15, e35, e1318.10670.99751.42176.03941.00001.6812e18, e26, e15, e35, e1325.84690.99801.25324.84681.00001.4586e18, e26, e15, e35, e1333.29580.99890.84662.54421.00000.9858e18, e26, e15, e35, e1343.29580.99620.84662.54421.00000.9858e18, e26, e15, e35, e13515.40350.99261.504714.50191.00001.7885e18, e26, e15, e35, e13616.59980.99022.121210.66401.00002.5000e18, e26, e15, e35, e13720.76170.97912.860011.78041.00003.1844e33, e26, e34, e11, e13867.94960.90945.158332.08320.99615.5973e33, e27, e34, e10, e149260.68590.92503.739493.82670.96173.7394e20, e28, e7, e9, e1410255.55330.95394.905976.62670.97855.0553e19, e5, e34, e8, e1411145.13870.95395.249763.27030.99126.3646e7, e28, e8, e10, e1412135.43370.96034.255857.90870.98985.2516e3, e28, e8, e9, e1313123.24610.96254.732956.20660.99125.7489e7, e28, e8, e9, e1414126.85130.93735.364860.62880.99236.4603e7, e28, e8, e10, e1415559.01250.85131.5438194.94420.94633.5852e6, e37, e15, e8, e 1416174.54500.94715.485880.62160.97065.4860e20, e5, e15, e10, e 1417121.46350.96435.048072.54540.97065.0480e33, e4, e17, e11, e1418165.78090.92103.3588126.01270.94923.3588e33, e28, e32, e10, e1419194.79160.95166.071199.96800.97346.1836e33, e3, e17, e35, e1320166.38800.93825.163697.50130.96285.1636e33, e4, e17, e11, e132154.06690.97683.454128.32970.98593.5698e33, e28, e17, e35, e132230.28930.99082.655420.01710.99593.2057e18, e26, e15, e35, e132313.17780.99581.82979.57390.99942.1778e18, e26, e15, e35, e13Figure 7
the simulation results for all the cases in a day.Figure 8
Power loss of each branch at 11:00 and 21:00.Figure 9
Voltage profile of each node at 11:00 and 21:00.Fig 10
RDG utilization at of each time interval in case 6 and case 7.


### 4.1.3. Comparing the results of HPSO and Other Methods

As a comparison, the optimization problem of ADN considering time variation feature is also addressed by the methods of BPSO [19], BAS (beetle antennae search algorithm) [37], GA [43] and MFO [44].Using all the methods, the list of final open switches for the best solutions during 24 hours are listed in Table5. The crucial constraint for the radial structure is satisfied for each topology. According to the optimal topologies, Fig 11 shows the power loss, minimum voltage standard amplitude, and RDG utilization based on case 7 using HPSO, BPSO, BAS, GA and MFO for the ADN. The obtained results show that the indexes of the proposed HPSO are better than other methods for ADN optimization considering time variations in load and RDG over 24 hours.Table 5
The open switches of four algorithms at each time interval.
Houropen switchesHPSOBPSOBASGAMFO0e18, e26, e15, e35, e13e7, e14, e9, e32, e28e3, e14, e11, e32, e28e18, e26, e15, e35, e13e18, e26, e15, e35, e131e18, e26, e15, e35, e13e7, e14, e9, e32, e28e33, e14, e11, e32, e3e18, e26, e15, e35, e13e18, e26, e15, e35, e132e18, e26, e15, e35, e13e7, e14, e9, e32, e28e33, e14, e10, e36, e3e18, e26, e16, e11, e12e18, e25, e16, e11, e123e18, e26, e15, e35, e13e7, e14, e9, e32, e28e33, e13, e10, e17, e3e18, e26, e16, e35, e13e18, e25, e16, e35, e134e18, e26, e15, e35, e13e7, e14, e9, e32, e28e33, e14, e11, e32, e4e18, e26, e16, e35, e13e18, e26, e16, e11, e125e18, e26, e15, e35, e13e7, e14, e9, e32, e28e33, e13, e11, e17, e3e18, e27, e16, e35, e13e18, e26, e16, e35, e136e18, e26, e15, e35, e13e7, e14, e9, e32, e37e33, e14, e11, e32, e3e19, e26, e15, e21, e13e19, e26, e15, e21, e137e33, e26, e34, e11, e13e7, e14, e9, e32, e28e33, e14, e11, e16, e4e33, e25, e15, e11, e13e33, e26, e15, e11, e138e33, e27, e34, e10, e14e7, e14, e9, e32, e28e33, e11, e21, e15, e4e33, e27, e34, e10, e14e33, e27, e34, e9, e139e20, e28, e7, e9, e14e7, e14, e9, e32, e28e6, e34, e9, e14, e28e20, e28, e7, e9, e14e20, e28, e6, e9, e1410e19, e5, e34, e8, e14e7, e14, e9, e32, e28e6, e14, e10, e15, e28e19, e5, e34, e8, e14e19, e 5, e8, e9, e1411e7, e28, e8, e10, e14e7, e14, e9, e32, e28e5, e14, e10, e15, e28e7, e28, e34, e8, e14e7, e14, e9, e32, e2812e3, e28, e8, e9, e13e7, e14, e9, e32, e37e33, e9, e6, e14, e28e7, e4, e8, e9, e14e7, e4, e8, e10, e1413e7, e28, e8, e9, e14e7, e14, e9, e32, e37e5, e13, e9, e15, e28e7, e28, e8, e10, e14e7, e28, e8, e10, e1414e7, e28, e8, e10, e14e7, e14, e9, e32, e28e33, e34, e9, e14, e5e7, e28, e8, e10, e14e7, e28, e8, e10, e1415e6, e37, e15, e8, e 14e7, e14, e9, e32, e28e6, e34, e9, e14, e37e6, e37, e15, e9, e 34e6, e37, e15, e9, e 3416e20, e5, e15, e10, e 14e7, e14, e9, e32, e37e5, e13, e9 15, e28e20, e5, e15, e9, e 14e20, e5, e15, e10, e 1417e33, e4, e17, e11, e14e7, e14, e9, e32, e28e33, e14, e11, e17, e4e33, e4, e32, e11, e14e33, e4, e17, e11, e1418e33, e28, e32, e10, e14e7, e14, e9, e32, e37e33, e11, e6, e14, e28e33, e28, e32, e10, e14e33, e28, e32, e10, e1419e33, e3, e17, e35, e13e7, e14, e9, e32, e37e33, e13, e11, e36 , e4e33, e3, e16, e35, e13e33, e3, e17, e35, e1320e33, e4, e17, e11, e13e7, e14, e9, e32, e28e33, e14, e11, e36, e4e33, e4, e17, e11, e13e33, e4, e17, e11, e1321e33, e28, e17, e35, e13e7, e14, e9, e32, e28e33, e14, e11, e36, e3e33, e28, e17, e35, e13e33, e28, e17, e35, e1322e18, e26, e15, e35, e13e7, e14, e9, e32, e37e33, e13, e10, e36, e3e18, e26, e15, e35, e13e18, e26, e15, e35, e1323e18, e26, e15, e35, e13e7, e14, e9, e32, e28e33, e14, e11, e36, e3e18, e2, e15, e35, e12e18, e26, e15, e 35, e13Table 6
Obtained results.
ItemsTotal Ploss/MWhRDG%Vmin/puHPSO1.229075.37%0.9463GA1.233575.22%0.9461MFO1.230275.28%0.9463BPSO2.303150.98%0.8848BAS1.344654.53%0.9409Figure 11
Comparison of power loss at each time interval.The power loss of the 33-bus ADN after optimization during 24 hours is shown graphically in Fig11. The maximum total power loss for the methods over 24 hours is 194.9442 kWh (HPSO), 389.0171 kWh (BPSO), 227.4633 kWh (BAS), 195.7518 kWh (GA), and 195.2001 kWh (MFO) at 15 p.m.Fig12 shows the RDG utilization after optimization at each time interval of a day. The results show that the minimum RDG utilization for each algorithm over 24 hours is 19.6811%, 9.0749%, 11.0089%, 19.6928%, and 19.6811%. The maximum RDG utilization for each algorithm over 24 hours is 100%, 92.9588%, 99.9970%, 99.9940%, and 100%. Moreover, the RDG utilization has reached 100% in eight-time intervals, which is obviously better than other methods.Figure 12
Comparison of RDG utilization at each time interval.The minimum voltage profiles of the ADN after optimization in a day are illustrated in Fig13. The day minimum voltage amplitude of each algorithm is 0.9463 p.u., 0.8848 p.u., 0.9409p.u., 0.9461 p.u., and 0.9463 p.u. at 15 p.m. The minimum node voltage for a day is almost the same as that of MFO, but 0.21%, 0.57%, and 6.95% higher than GA, BAS, and BPSO.Figure 13
Comparison of minimum voltage at each time interval.Table7 shows the total active power loss and RDG utilization in a day. The day power loss by the proposed HPSO is 1.2290MWh which is lower than 1.2302MWh (GA), 1.2335MWh (MFO), 1.3346MWh (BAS), and 2.3031MWh (BPSO). The RDG utilization for a day using HPSO, GA, MFO, BAS, and BPSO are 75.37%, 75.22%, 75.28%, 54.53%, and 50.98%, respectively. From all the results, we can see that the obtained results of HPSO algorithm are better to compare the GA, MFO, BAS, and BPSO. Thus, the proposed HPSO method is better in global optimization ability.Table 7
The simulation results of 69-bus ADN in a day.
ItemsTotal Ploss/MWTotal RDG/MWVmin/puCase 494.440400.4092Case 51.875200.8079Case 6180.967483.51320.5476Case 71.4379123.61390.9203


## 4.1.1. Benefits of combining network optimization and RDG integration

To analyze the operation performance of ADN after RDG integration and network optimization, three cases are simulated to utilize the proposed method at peak load point.Case 1.
Initial system without RDGs and network optimization.Case 2.
Solving the optimal problem of RDG utilization maximization and power loss minimization without network optimization.Case 3.
Solving the problem of RDG utilization maximization and power loss minimization with network optimization.
In Fig4, the peak load point of actual hourly load demand is 7.4406MW+4.5314MVar at 19:00. The obtained results of the cases are presented in Table 1.
In case 2, the RDG is optimally integrated utilization without network optimization, providing an 80.2% reduction of system power loss (Ploss). However, the Ploss is reduced by 89.53% in case 3. Moreover, the absorption RDG power (PRDG) by the system has improved to 6.1836MW after optimization. The combination of network optimization and RDG integration is effective in reducing the system power loss of the ADN.
Fig5 shows the hourly voltage profiles of all the cases. The voltage amplitude is improved when combining network optimization with RDG integration than when integrating RDG only and base system. The minimum voltage (Vmin) amplitude of case 3 has increased to 0.9734p.u., which is above 0.9494p.u. of case 2 and 0.8104 p.u. of case 1. The results show that system efficiency in the optimization problem with adding the RDG is better than the base system and has resulted in more power loss reduction, more voltage profile change, and more RDG power utilization. The use of RDG along with network optimization leads to better operation conditions of the ADN due to the injection of RDG power.
To justify the superiority of the proposed method, the optimization problem in case 3 is also solved by the conventional algorithms, GA method [43] and moth-flame optimization (MFO) [44]. the obtained simulation results are compared with those obtained by the conventional methods. The methods are based on initial conditions and similar assumptions. Table 2 and Fig 6 illustrate the detailed comparison results.
The proposed method presents a superior performance to the GA and MFO. The proposed method provides an optimal solution with 99.9680kW of system power loss which is less than 101.3932kW in GA algorithm. According to Fig6, the proposed method could achieve optimal solution faster than other methods, i. e., 47 iterations against 95 in MFO, 47 iterations compared to 157 in GA. Consequently, the proposed method is an effective way for addressing ADN network optimization due to its high convergence ability and accuracy.Table 1
Benefits of combining network optimization and RDG integration.
CasesOpen switchesPloss/kWPRDG/MWVmin/p.u.1e33-e37954.864400.81042e33-e37189.03625.94230.94943e33, e3, e17, e35, e1399.96806.18360.9734Figure 5
The voltage profile for cases.Table 2
Comparison with last optimization methods.
ItemsOpen switchesPloss/kWPRDG/MWVmin/p.u.Proposed methode33, e3, e17, e35, e1399.96806.18360.9734GA [43]e33, e28, e17, e35, e13101.39326.18360.9662MFO [44]e33, e3, e17, e35, e1399.96806.18360.9734Fig 6
The convergence curve of the IEEE33 bus system.


## 4.1.2. Experiment considering time variations in RDG and load

The proposed model and method is used to address the optimization problem of ADN considering the time variations in RDG and load. The hourly data in RDG and load of an average day as an example of study are illustrated in Fig4. In the 33-bus ADN, four cases are adopted to investigate the effect of time-varying characteristics on network optimization indicators.Case 4.
Load variation is considered without network optimization and RDG.Case 5.
Load variation and network optimization simulation are considered.Case 6.
Time variations in load and RDG are considered without network optimization.Case 7.
Simultaneous hourly network optimization and time variations in load and RDG are studied.
Using the proposed model and method, we obtained the selected open switches and RDG utilization for a given time certain. Tables3 and 4 present the simulation results and optimal solutions for each time interval for the cases. In each time interval, the open switches, Ploss, Vmin, and absorbed RDGs power are presented corresponding to each time interval for cases.
In cases 4 and 5, the simulation is carried out in the system with time variation of load demand. Table3 lists the results of active power loss, minimum node voltage, and open switches before and after optimization (case 4 and case 5). As seen in Table 3, the optimal solutions with opening switches e7, e14, e9, e32, and e28 are consistent at time intervals 0∼4th, 8∼11th, 13∼15th, 21th, and 23th. The opened switches in other time intervals are e7, e14, e9, e32, and e37 after optimization. The active power loss for each time interval reduces after optimization, i.e, 41.4271kW against 59.2882kW at 6:00, 458.5657kW compared to 707.9717kW at 10:00, 458.6379kW instead of 631.8764kW at 16:00, and so on. Further, it is noted that the minimum node voltage for each interval has increased in case 5 compared in case 4. For example, The percentage improvement in minimum node voltage value after optimization is 1.67%, 5.08%, and 6.40% at 6:00, 10:00, and 16:00, respectively. Fig 7 gives the total power loss and minimum voltage profile in a day of the cases. From Fig 7, the total daily power loss in case 4 is 10.0352MW. The daily total power loss after optimization has been reduced by 34.69%. The minimum voltage amplitude in a day improves to 0.8646p.u. in case 5, which implies an improvement of 6.69% compared to case 4. The results illustrate network optimization can effectively improve the system operation level. However, the daily minimum voltage amplitude 0.8646p.u. in case 5 is lower than the voltage limit of 0.90p.u. It is not conducive to the long-term stable operation of the active distribution system. The way of integration of RDG is considered to improve operating indices.
To perform the suggested method which requires time variations in RDG and load demand, different load profiles and RDG output power are added to the system at each time interval. Table4 depicts the variations of the operational indices of the system in cases 6 and 7. It is observed that the optimization effects of cases 6 and 7 are better than those of cases 4 and 5. In cases 6 and 7, the obtained optimal solutions provide the total daily power loss are 2.6758MW and 1.2290MW, respectively. The maximum power loss of the day in case 7 is 194.9442kw at 15:00, which implies a reduction of 762.7595kW, 263.0196kW, and 364.0683kW compared to case 4, case 5, and case 6. From Fig 7, the daily minimum voltage amplitude of case 7 is 0.9463p.u., which is higher than the minimum voltage limit 0.90p.u. It indicates that simultaneous network optimization with time-varying load and RDG can improve the performance of the ADN. Moreover, the day absorption RDG power for case 7 is 90.2614MW, which has increased by 12.42% with case 6. And the RDG utilization (RDG%) has been improved from 67.04% in case 6 to 75.37% in case 7.
Due to space reasons, the operational indices of the ADN at 11:00 and 21:00 are analyzed as examples to describe the operation characteristic at each time interval clearly. Fig8 shows the power loss of each branch at 11:00 and 21:00. In case 4, the power loss of ADN is relatively large without network optimization and integration of RDGs. In case 5, it is significantly reduced after optimization by changing the topology with opening switch combinations [e7, e14, e9, e32, e28]. The system power loss is reduced by 53.62% at 11:00, and 32.87% at 21:00, respectively. After integrating of RDGs, The system power loss for case 6 is reduced by 78.27% at 11:00, and 77.95% at 21:00, respectively. The benefits of simultaneously combining network optimization with the integration of RDG are investigated in case 7. The obtained final opening switch combinations at 11:00 and 21:00 are [e7, e28, e8, e10, e14] and [e33, e28, e17, e35, e13], respectively. It provides a reduction to 63.2703kW at 11:00 and 28.3297kW at 21:00 of system power loss, which is significantly lower than that of other cases.
Fig9 shows the voltage profile at two-time interval points. In case 4, voltage violations occur in most of the nodes which leads to an unsatisfactory low voltage profile. The network performance is improved by using network optimization and integration of RDGs simultaneously. It is observed that the system voltage volatility range is reduced in case 7. And the minimum voltage profile in cases 4-7 are 0.8449 p.u. for node 33 at 11:00, 0.8943 p.u. for node 32 at 11:00, 0.9539 p.u. for node 33 at 11:00, 0.9912 p.u. for node 24 at 11:00. The minimum voltage of ADN has raised dramatically with the network optimization and RDG integration.
As seen in Table4, the absorption RDG power at two-time intervals by the system of case 6 are 5.2497 MW and 3.4541MW, while the optimization system of case 7 reaches 6.3646MW and 3.5698MW. It turns out that the absorption power of RDGs in case 6 is increased by 21.24% at 11:00 and 3.4% at 21:00 in case 6. The RDG utilization of a day in case 7 has increased by 15.48% and 3.24% compared to case 6. Moreover, Fig 10 compares the absorption RDG power over 24 hours for cases 6 and 7. The minimum RDG utilization in case 7 is 19.68% at 17:00, which is higher 16.64% in case 6. In both cases, the maximum difference of RDG utilization has reached to 56.94%. Thus, considering the time variations in load and RDG, the operational performance of the ADN would improve by using the proposed dynamic optimization model and the HPSO method.Table 3
The simulation results for cases 4 and 5.
HourCase 4Case 5Ploss/kWVmin/p.u.Open switchesPloss/kWVmin/p.u.Open switches022.71690.9708e33, e34, e35, e36, e3715.42030.9800e7, e14, e9, e32, e28122.71680.970815.42030.9800216.30750.974610.45090.983636.05090.98494.09180.990046.05070.98494.09180.9900526.48390.967818.62090.9774e7, e14, e9, e32, e37659.28820.952141.42710.96757102.65280.938171.65230.95578454.68530.8711299.25720.9111e7, e14, e9, e32, e289634.78270.8465406.95870.895810707.97170.8399458.56570.890711667.97350.8449434.82380.894312450.93750.8707303.20610.9095e7, e14, e9, e32, e3713524.29160.8622353.28660.9055e7, e14, e9, e32, e2814655.23440.8467429.30500.894515957.70370.8154617.22440.873516972.44770.8126631.87640.8646e7, e14, e9, e32, e3717800.69860.8283523.44940.878618704.07660.8373460.81270.887019954.86440.8104612.00810.869220892.39200.8166574.12480.872621245.25010.9043164.64890.9364e7, e14, e9, e32, e2822107.07330.935974.21700.9554e7, e14, e9, e32, e372342.59500.960329.01500.9729e7, e14, e9, e32, e28Table 4
the simulation results for cases 6 and 7.
HourCase 6Case 7Ploss/kWVmin/puPRDG/MWOpen switchesPloss/kWVmin/puPRDG/MWOpen switches08.10670.99751.4218e33, e34, e35, e36, e376.03941.00001.6812e18, e26, e15, e35, e1318.10670.99751.42176.03941.00001.6812e18, e26, e15, e35, e1325.84690.99801.25324.84681.00001.4586e18, e26, e15, e35, e1333.29580.99890.84662.54421.00000.9858e18, e26, e15, e35, e1343.29580.99620.84662.54421.00000.9858e18, e26, e15, e35, e13515.40350.99261.504714.50191.00001.7885e18, e26, e15, e35, e13616.59980.99022.121210.66401.00002.5000e18, e26, e15, e35, e13720.76170.97912.860011.78041.00003.1844e33, e26, e34, e11, e13867.94960.90945.158332.08320.99615.5973e33, e27, e34, e10, e149260.68590.92503.739493.82670.96173.7394e20, e28, e7, e9, e1410255.55330.95394.905976.62670.97855.0553e19, e5, e34, e8, e1411145.13870.95395.249763.27030.99126.3646e7, e28, e8, e10, e1412135.43370.96034.255857.90870.98985.2516e3, e28, e8, e9, e1313123.24610.96254.732956.20660.99125.7489e7, e28, e8, e9, e1414126.85130.93735.364860.62880.99236.4603e7, e28, e8, e10, e1415559.01250.85131.5438194.94420.94633.5852e6, e37, e15, e8, e 1416174.54500.94715.485880.62160.97065.4860e20, e5, e15, e10, e 1417121.46350.96435.048072.54540.97065.0480e33, e4, e17, e11, e1418165.78090.92103.3588126.01270.94923.3588e33, e28, e32, e10, e1419194.79160.95166.071199.96800.97346.1836e33, e3, e17, e35, e1320166.38800.93825.163697.50130.96285.1636e33, e4, e17, e11, e132154.06690.97683.454128.32970.98593.5698e33, e28, e17, e35, e132230.28930.99082.655420.01710.99593.2057e18, e26, e15, e35, e132313.17780.99581.82979.57390.99942.1778e18, e26, e15, e35, e13Figure 7
the simulation results for all the cases in a day.Figure 8
Power loss of each branch at 11:00 and 21:00.Figure 9
Voltage profile of each node at 11:00 and 21:00.Fig 10
RDG utilization at of each time interval in case 6 and case 7.


## 4.1.3. Comparing the results of HPSO and Other Methods

As a comparison, the optimization problem of ADN considering time variation feature is also addressed by the methods of BPSO [19], BAS (beetle antennae search algorithm) [37], GA [43] and MFO [44].Using all the methods, the list of final open switches for the best solutions during 24 hours are listed in Table5. The crucial constraint for the radial structure is satisfied for each topology. According to the optimal topologies, Fig 11 shows the power loss, minimum voltage standard amplitude, and RDG utilization based on case 7 using HPSO, BPSO, BAS, GA and MFO for the ADN. The obtained results show that the indexes of the proposed HPSO are better than other methods for ADN optimization considering time variations in load and RDG over 24 hours.Table 5
The open switches of four algorithms at each time interval.
Houropen switchesHPSOBPSOBASGAMFO0e18, e26, e15, e35, e13e7, e14, e9, e32, e28e3, e14, e11, e32, e28e18, e26, e15, e35, e13e18, e26, e15, e35, e131e18, e26, e15, e35, e13e7, e14, e9, e32, e28e33, e14, e11, e32, e3e18, e26, e15, e35, e13e18, e26, e15, e35, e132e18, e26, e15, e35, e13e7, e14, e9, e32, e28e33, e14, e10, e36, e3e18, e26, e16, e11, e12e18, e25, e16, e11, e123e18, e26, e15, e35, e13e7, e14, e9, e32, e28e33, e13, e10, e17, e3e18, e26, e16, e35, e13e18, e25, e16, e35, e134e18, e26, e15, e35, e13e7, e14, e9, e32, e28e33, e14, e11, e32, e4e18, e26, e16, e35, e13e18, e26, e16, e11, e125e18, e26, e15, e35, e13e7, e14, e9, e32, e28e33, e13, e11, e17, e3e18, e27, e16, e35, e13e18, e26, e16, e35, e136e18, e26, e15, e35, e13e7, e14, e9, e32, e37e33, e14, e11, e32, e3e19, e26, e15, e21, e13e19, e26, e15, e21, e137e33, e26, e34, e11, e13e7, e14, e9, e32, e28e33, e14, e11, e16, e4e33, e25, e15, e11, e13e33, e26, e15, e11, e138e33, e27, e34, e10, e14e7, e14, e9, e32, e28e33, e11, e21, e15, e4e33, e27, e34, e10, e14e33, e27, e34, e9, e139e20, e28, e7, e9, e14e7, e14, e9, e32, e28e6, e34, e9, e14, e28e20, e28, e7, e9, e14e20, e28, e6, e9, e1410e19, e5, e34, e8, e14e7, e14, e9, e32, e28e6, e14, e10, e15, e28e19, e5, e34, e8, e14e19, e 5, e8, e9, e1411e7, e28, e8, e10, e14e7, e14, e9, e32, e28e5, e14, e10, e15, e28e7, e28, e34, e8, e14e7, e14, e9, e32, e2812e3, e28, e8, e9, e13e7, e14, e9, e32, e37e33, e9, e6, e14, e28e7, e4, e8, e9, e14e7, e4, e8, e10, e1413e7, e28, e8, e9, e14e7, e14, e9, e32, e37e5, e13, e9, e15, e28e7, e28, e8, e10, e14e7, e28, e8, e10, e1414e7, e28, e8, e10, e14e7, e14, e9, e32, e28e33, e34, e9, e14, e5e7, e28, e8, e10, e14e7, e28, e8, e10, e1415e6, e37, e15, e8, e 14e7, e14, e9, e32, e28e6, e34, e9, e14, e37e6, e37, e15, e9, e 34e6, e37, e15, e9, e 3416e20, e5, e15, e10, e 14e7, e14, e9, e32, e37e5, e13, e9 15, e28e20, e5, e15, e9, e 14e20, e5, e15, e10, e 1417e33, e4, e17, e11, e14e7, e14, e9, e32, e28e33, e14, e11, e17, e4e33, e4, e32, e11, e14e33, e4, e17, e11, e1418e33, e28, e32, e10, e14e7, e14, e9, e32, e37e33, e11, e6, e14, e28e33, e28, e32, e10, e14e33, e28, e32, e10, e1419e33, e3, e17, e35, e13e7, e14, e9, e32, e37e33, e13, e11, e36 , e4e33, e3, e16, e35, e13e33, e3, e17, e35, e1320e33, e4, e17, e11, e13e7, e14, e9, e32, e28e33, e14, e11, e36, e4e33, e4, e17, e11, e13e33, e4, e17, e11, e1321e33, e28, e17, e35, e13e7, e14, e9, e32, e28e33, e14, e11, e36, e3e33, e28, e17, e35, e13e33, e28, e17, e35, e1322e18, e26, e15, e35, e13e7, e14, e9, e32, e37e33, e13, e10, e36, e3e18, e26, e15, e35, e13e18, e26, e15, e35, e1323e18, e26, e15, e35, e13e7, e14, e9, e32, e28e33, e14, e11, e36, e3e18, e2, e15, e35, e12e18, e26, e15, e 35, e13Table 6
Obtained results.
ItemsTotal Ploss/MWhRDG%Vmin/puHPSO1.229075.37%0.9463GA1.233575.22%0.9461MFO1.230275.28%0.9463BPSO2.303150.98%0.8848BAS1.344654.53%0.9409Figure 11
Comparison of power loss at each time interval.The power loss of the 33-bus ADN after optimization during 24 hours is shown graphically in Fig11. The maximum total power loss for the methods over 24 hours is 194.9442 kWh (HPSO), 389.0171 kWh (BPSO), 227.4633 kWh (BAS), 195.7518 kWh (GA), and 195.2001 kWh (MFO) at 15 p.m.Fig12 shows the RDG utilization after optimization at each time interval of a day. The results show that the minimum RDG utilization for each algorithm over 24 hours is 19.6811%, 9.0749%, 11.0089%, 19.6928%, and 19.6811%. The maximum RDG utilization for each algorithm over 24 hours is 100%, 92.9588%, 99.9970%, 99.9940%, and 100%. Moreover, the RDG utilization has reached 100% in eight-time intervals, which is obviously better than other methods.Figure 12
Comparison of RDG utilization at each time interval.The minimum voltage profiles of the ADN after optimization in a day are illustrated in Fig13. The day minimum voltage amplitude of each algorithm is 0.9463 p.u., 0.8848 p.u., 0.9409p.u., 0.9461 p.u., and 0.9463 p.u. at 15 p.m. The minimum node voltage for a day is almost the same as that of MFO, but 0.21%, 0.57%, and 6.95% higher than GA, BAS, and BPSO.Figure 13
Comparison of minimum voltage at each time interval.Table7 shows the total active power loss and RDG utilization in a day. The day power loss by the proposed HPSO is 1.2290MWh which is lower than 1.2302MWh (GA), 1.2335MWh (MFO), 1.3346MWh (BAS), and 2.3031MWh (BPSO). The RDG utilization for a day using HPSO, GA, MFO, BAS, and BPSO are 75.37%, 75.22%, 75.28%, 54.53%, and 50.98%, respectively. From all the results, we can see that the obtained results of HPSO algorithm are better to compare the GA, MFO, BAS, and BPSO. Thus, the proposed HPSO method is better in global optimization ability.Table 7
The simulation results of 69-bus ADN in a day.
ItemsTotal Ploss/MWTotal RDG/MWVmin/puCase 494.440400.4092Case 51.875200.8079Case 6180.967483.51320.5476Case 71.4379123.61390.9203


## 4.2. 69-bus ADN

The 69-bus ADN consists of 69 buses and 5 tie branches, as provided in Fig14 [44]. The wind RDGs are connected to nodes 11 and 50. The PV RDGs are connected to nodes 38 and 53. The open switches in the initial network are e69, e70, e71, e72, and e73. The proposed method is applied to obtain the optimal reconfigure to enable RDG penetration. The hourly load demand for the system is obtained by using the method in [43].Figure 14
69-bus ADN.For the sake of investigating the optimal problem of the ADN considering time variations in load and RDG, cases 4-7 will be implemented to the 69-bus ADN. Table7 provides simulated results of the test active distribution network for 24 hours using the proposed method.The original network considers the time-varying nature of load demand without the integration of RDGs and network optimization. From Table7, we observe that the total daily active power loss of 94.4404MW for case 4. After network optimization, the total power loss has been reduced to 1.8752MW in case 5. To discuss the RDG effect, consider the integration of RDGs as case 6. The total power loss of the case is calculated as 169.0020MW. When both network optimization and RDG integration are considered, the daily power loss is minimized to 1.4379MW. Fig 15 shows active power loss for the four cases during 24 hours. In all hours, the power loss has been reduced greatly using network optimization and integration of RDGs in case 7.Figure 15
Power loss of 69-bus ADN during 24 hours.The voltage profiles of the cases during 24 hours have been illustrated in Fig16. The hourly minimum voltages in cases are 0.4092p.u. at 15:00, 0.8079 p.u. at 21:00, 0.5476 p.u. at 8:00, and 0.9203 p.u. at 14:00, respectively. It is observed that the minimum voltage profile for cases 4-6 is lower than the lower voltage limit 0.9p.u. The network optimization based on time variation of RDGs and load demand brings much less voltage fluctuations with respect to case 7 than case 5 and case 6 because the scheme considers the network optimization and integration of RDGs.Figure 16
Voltage profiles of 69-bus ADN during 24 hours.Fig17 also compares the absorbed RDG power by the system at each interval before and after network optimization (case 6 and case 7). We observe that the total absorbed RDG power after optimization for a day is 123.6139MW, which has increased by 45.40% before the system. The RDG utilization has been improved from 52.93% to 76.96% after optimization. With the expansion of network structure, the proposed method is more effective for improving the operation indices of the ADN.Figure 17
Absorbed RDG power of 69-bus ADN during 24 hours.Table8 shows the open switches of the hourly optimization solution for the cases. The application of network optimization to the ADN has permitted to achieve a better operational conditions. As an example, the topological structure [e7, e18, e13, e50, e47] is optimal for the hour 8:00 of the day in case 7. It ensures a reduction of 93.37%, 49.97%, and 95.32% of power loss with other cases. In addition, The voltage profile is increased by 42.87%, 8.68% and 66.58% for case 4, case 5, and case 6, respectively. And the absorbed RDG power of case 7 is improved to 6.6721MW which is higher than 4.8958MW of case 6. These results show that the proposed model and method yield acceptable results in terms of power loss reduction, voltage improvement, and RDG utilization for both test systems under time variations of load demand and RDG. This demonstrates the ability of the proposed approach to find the optimal network topology along with absorbed RDG power and load demand for the ADN considering time characteristics.Table 8
Open switches of hourly optimization for 69-bus ADN.
HourCase 4Case 5Case 6Case 70e69, e70, e71, e72, e73e64, e14, e68, e44, e38e69, e70, e71, e72, e73e4, e18, e13, e53, e471e3, e18, e67, e50, e42e4, e17, e14, e26, e472e3, e18, e12, e47, e7e4, e18, e13, e53, e473e69, e19, e62, e50, e47e4, e18, e13, e53, e474e4, e16, e14, e53, e38e4, e18, e13, e25, e475e4, e17, e68, e52, e38e4, e18, e13, e72, e476e65, e17, e71, e52, e36e4, e18, e13, e24, e477e10, e15, e13, e12, e35e4, e16, e13, e20, e418e9, e15, e68, e20, e37e7, e18, e13, e50, e479e3, e16, e11, e13, e36e3, e15, e5, e12, e4510e3, e17, e 68, e11, e4e3, e17, e9, e12, e711e8, e70, e68, e11, e36e69, e13, e 67, e48, e712e9, e70, e12, e16 , e36e61, e16, e7, e12, e4113e9, e15, e13, e20, e36e10, e15, e68, e11, e514e3 e70, e 66, e11, e7e3, e15, e71, e50, e4115e10, e15, e13, e23, e37e69, e13, e5, e49, e3516e8, e16, e13, e11, e36e65, e16, e14, e20, e4717e8, e15, e71, e11, e35e3, e17, e68, e20, e718e8, e17, e11, e13, e36e60, e15, e9, e11, e3619e8, e15, e66, e12, e38e3, e19, e68, e21, e520e9, e15, e68, e22, e36e65, e16, e68, e53, e721e9, e18, e12, e14, e35e4, e18, e68, e20, e3622e10, e15, e68, e25, e36e4, e18, e13, e72, e4723e3, e19, e11, e26, e36e4, e15, e68, e25, e47To further illustrate the superiority of the HPSO, we compared and analyzed the 24 hours simulated results of case 7 by methods of HPSO, BPSO, BAS, GA, and MFO. Fig18–20 show the comparison of the proposed HPSO and the conventional approaches BPSO, BAS, GA, and MFO. And the comparison of the simulation results in a day is also shown in Fig 21–22. We can clearly the minimum voltage profiles in a day are 0.9203 p.u. for HPSO, 0.8487 p.u. for BPSO, 0.9048 p.u. for BAS, 0.9009 p.u. for GA, 0.9066 p.u. for MFO, respectively. From Fig 18, the total active power loss of a day 1.4379 MWh by HPSO, which is less 77.57%, 67.07%, 34.34%, and 21.09% than BPSO, BAS, GA, and MFO. The day RDG utilization with a penetration rate of 76.96% is the highest among all the methods. Thus, considering the overall operation economy in terms of power loss, voltage quality, and RDG utilization, the proposed HPSO is superior to the other methods.Figure 18
Comparison of power loss of 69-bus ADN.Figure19
Comparison of RDG utilization of 69-bus ADN.Figure 20
Comparison of minimum voltage of 69-bus ADN.Figure 21
Comparison total power loss and minimum voltage profile for all the methods in a day.Figure 22
Comparison RDG power and RDG utilization for all the methods in a day.


## 5. Conclusions

In this paper, a successful multi-period optimization strategy of the ADN considering time variations in load demand and RDG has been presented. The objective function of the dynamic optimization model includes the power loss and the RDG utilization for a given time to ensure optimal operation of the ADN. Additionally, the average weigh factors approach is used to balance the relationship of the objectives in multi-period optimization model. Furthermore, the proposed HPSO algorithm and the optimization model are combined to determine the optimal topology and RDG resource allocation scheme for each time interval in the ADN.The optimization strategy considers different possibilities of network optimization, time variations of load demand and RDG. The proposed optimization model and HPSO method test on the 33-bus and 69-bus ADNs to verify the effectiveness. The simulation results show that the proposed strategy obtains reasonable and high-quality schedules for switching and the active power values of RDGs in multi-period optimization objective frameworks. Additionally, considering the network optimization helps significantly increase RDG utilization and improve the voltage profile in a day. Case studies are conducted on 33-bus and 69-bus ADNs under different operating conditions. Simulation results show that the proposed strategy can find the corresponding optimal solution in each time interval. Thus, the performance of the ADN can be improved by integrating RDG and network optimization. Moreover, the superiority of the proposed optimization algorithm considering time variations of load and RDG is proved by comparing conventional methods. This paper presents a management strategy for ADN considering time variations in load and RDG, which is an important requirement to mostly explore the benefits without extra investment costs or system risks. However, the switches in the optimization process would be changed constantly according to time intervals. The influence of frequent switch action is ignored in this paper. The author suggestion for future study would be to modify the method according to the switch constraint. [21, 45].



---

*Source: 5771094-2022-12-07.xml*