Next Article in Journal
Models of Teaching Science Communication
Next Article in Special Issue
Optimal Scheduling of AC–DC Hybrid Distribution Network Considering the Control Mode of a Converter Station
Previous Article in Journal
Shaping the Knowledge Base of Bioeconomy Sectors Development in Latin American and Caribbean Countries: A Bibliometric Analysis
Previous Article in Special Issue
Review on the Application of Photovoltaic Forecasting Using Machine Learning for Very Short- to Long-Term Forecasting
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Allocation and Sizing of Distributed Generation Using Interval Power Flow

by
Wallisson C. Nogueira
1,
Lina P. Garcés Negrete
1 and
Jesús M. López-Lezama
2,*
1
Electrical, Mechanical and Computer Engineering School, Federal University of Goiás, Av. Universitária No. 1488, Goiânia 74605-010, Brazil
2
Grupo en Manejo Eficiente de la Energía (GIMEL), Departamento de Ingeniería Eléctrica, Universidad de Antioquia (UdeA), Calle 70 No. 52-21, Medellin 050010, Colombia
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(6), 5171; https://doi.org/10.3390/su15065171
Submission received: 8 February 2023 / Revised: 4 March 2023 / Accepted: 9 March 2023 / Published: 14 March 2023
(This article belongs to the Special Issue Smart Grid Technologies and Applications)

Abstract

:
Modern distribution systems and microgrids must deal with high levels of uncertainty in their planning and operation. These uncertainties are mainly due to variations in loads and distributed generation (DG) introduced by new technologies. This scenario brings new challenges to planners and system operators that need new tools to perform more assertive analyses of the grid state. This paper presents an optimization methodology capable of considering uncertainties in the optimal allocation and sizing problem of DG in distribution networks. The proposed methodology uses an interval power flow (IPF) that adds uncertainties to the combinatorial optimization problem in charge of sizing and allocating DG units in the network. Two metaheuristics were implemented for comparative purposes, namely, symbiotic organism search (SOS) and particle swarm optimization (PSO). The proposed methodology was implemented in Python® using as benchmark distribution systems the IEEE 33-bus and IEEE 69-bus test distribution networks. The objective function consists of minimizing technical losses and regulating network voltage levels. The results obtained from the proposed IPF on the tested networks are compatible with those obtained by the PPF, thus evidencing the robustness and applicability of the proposed method. For the solution of the optimization problem, the SOS metaheuristic proved to be robust, since it was able to find the best solutions (lowest losses) while keeping voltage levels within the predetermined range. On the other hand, the PSO metaheuristic showed less satisfactory results, since for all test systems, the solutions found were of lower quality than the ones found by the SOS.

1. Introduction

Electric power systems (EPSs) play an important role in modern societies as they enable the use of technologies and supply electrical energy to industries and households. Nevertheless, their proper functioning can only be guaranteed if they have an adequate set of tools for their operation and planning. Among these studies, the most common analysis is the power flow (PF) in its AC and DC variants, which can be computed by using several widely known techniques [1,2]. Still, the biggest weakness of these studies is that the accuracy of their results is as good as the accuracy of the input data [3]. Power flow inaccuracy can be due to measurement errors, inaccurate forecasts, assumptions of limits on load, or scheduled outages [4].
The current political and economic landscape of the energy industry has increasingly encouraged the use of distributed energy resources (DERs), which include small-scale consumer-side generation, energy storage equipment, and demand response (DR) [5,6,7,8]. Thus, the role of consumers has been progressively changing, becoming more active. In this context, load models must adapt to characterize these new uncertainties [9]. Interval analysis can be seen as an effective tool for considering these uncertainties in PF problems. In this case, loads are modeled by a set of values varying linearly from the smallest to the largest possible value, and thus nodal voltages, power flows, and losses are calculated in an interval manner. This approach is known as interval power flow (IPF) [10,11,12,13].
Some mathematical tools are needed to calculate the IPF and ensure its convergence. This can be accomplished either by Newton’s interval method or Krawczyk’s method [14], which finds the interval solution of a system with a given tolerance. In [13], an optimal distribution planning model is proposed using IPF for uncertainty addition, where the interval results are used as merit functions in an optimization algorithm. Complementing IPF, interval mathematics are also used for uncertainty addition in problems such as state estimation and short-circuit calculation [15].
Another approach to adding uncertainty is the probabilistic power flow (PPF) [16,17,18]. In this case, loads are modeled by means of their density functions. Several research works have proposed mathematical models to solve this problem. Monte Carlo simulation (MCS) is the most widely used method for approaching a PPF. This method consists of solving several deterministic PF with random variations of their parameters based on their probability density functions, so if a large enough sample of results is obtained, a highly accurate result can be expected [19]. Probabilistic uncertainties can also be modeled by other methods, such as polynomial chaos expansion (PCE), which is able to work in a non-intrusive way and model random variables and stochastic processes with a high performance [20]. Table 1 presents the main characteristics of conventional power flow (CPF), IPF, and PPF.
IPF has been shown to be an effective tool to deal with uncertainties. In [21], the authors implemented an IPF analysis considering interval power injections from wind farms with the aim of guaranteeing interval values of the grid variables for safe system operation. In [22], an affine arithmetic-based PF is proposed for the consideration of unscheduled regional power fluctuations. The authors in [23] proposed an IPF based on a multi-stage affine arithmetic to consider the uncertainty of loads and DG in an unbalanced network. In [24], a hybrid between interval and probabilistic power flow is proposed using an analytical clustering method. In this case, the uncertainties of random and interval variables are considered simultaneously. The authors in [25] proposed an IPF to generate optimal scenarios. In this approach, the interval uncertainties are considered with variable bounds, and the objective function is programmed to compute such bound. In [26], the authors proposed an IPF based on the Taylor inclusion function, while in [27] an IPF is proposed using the correlation of the uncertainties of the power injections. Other variants of the IPF are also proposed in [28,29].
Even though power flow accuracy has improved, the distribution planning problem involving the integration of DG is still complex. Therefore, several studies have also been devoted to current and power control in distribution networks with DG [30,31,32]. Grid connections for DG are significantly different from traditional centralized generator connections. Several factors must be analyzed for planning the connection of DG, such as the number of generating sources, their locations, and connection types, among others [33,34,35].
The installation of DG in non-optimal locations brings along a series of problems to the network, including problems with protection coordination, over-voltages, and also an increment on technical losses, therefore, having the opposite effect of the desired one [36,37,38]. Thus, using a method capable of selecting the best locations and preferable sizing in a complex distribution network represents a useful tool for planners.
The problem of allocation and sizing of DG consists of two steps that are solved together, the selection of the best locations for the installation of the generation and the choice of the maximum capacity of each generating unit [39]. This problem has been widely studied using a series of different objectives. Such objectives can be minimizing losses [40], improving voltage profile [41], minimizing investments and maintenance costs [42] improvement of system reliability indices [43], and spinning reserve increment [44]; furthermore, there are proposals that combine several of these objectives [45,46,47].
The optimal location and sizing of DG is a complex optimization problem that can be tackled by metaheuristic techniques such as genetic and evolutionary algorithms [48], particle swarm optimization [49], symbiotic organism search (SOS) [50], gray wolf optimizer [51], crow search algorithm [52], artificial ecosystem-based optimization [53], and hybrid methodologies [54]. Recently, SOS and particle swarm optimization (PSO) [55,56] have demonstrated to be effective when solving complex optimization problems. PSO and SOS are optimization methodologies based on the behavior of swarms and symbiotic organisms, respectively.
This paper presents a methodology for the optimal allocation and sizing of DG in distribution networks that takes into account the existing uncertainties in loads and generators. The proposed methodology is based on an IPF for adding uncertainties; furthermore, PSO and SOS are used for the optimization of generation in the grid; nonetheless, any other metaheuristic technique might be applied. One of the strengths of IPF is that uncertainties are considered through the percentage variation of the variables of interest; then, it is possible to incorporate the expertise of the network operator to determine the most appropriate ranges of variation of these variables and therefore, obtain more accurate modeling of the operation. Furthermore, due to its nature, the IPF features a reduced computational effort when compared to probabilistic approaches that require thousands of simulations. In this sense, the main contributions of this paper are as follows: (i) it presents a methodology capable of adding and evaluating uncertainties in the process of optimal DG siting and sizing in distribution systems, (ii) two metaheuristic techniques are implemented for assessing the performance of the proposed approach; namely, PSO and SOS, and (iii) an analysis of the convergence of both metaheuristics operating with interval values is carried out along with an evaluation of their performance under two different metrics.

2. Uncertainty Modeling

The uncertainty modeling was carried out following the IPF implemented by the authors in [57] and is summarized below.
Step 1: Use a deterministic power flow to compute the bus voltages of the test system. In this case, the input variables may be fixed at their deterministic values. Once this is done, the percentage variations in demand and their respective interval values can be calculated using Equations (1) and (2). In this case,  P c k d  and  Q c k d  represent the deterministic values of the active and reactive loads at bus k, respectively.  P c k i  and  Q c k i  are the interval values of the active and reactive loads and  α  represents the percentage variation of the demand.
P c k i = P c k d 1 α , 1 + α
Q c k i = Q c k d 1 α , 1 + α
Step 2: The expected values of the interval active and reactive power are calculated using Equations (3) and (4). In this case,  P e x p k i  and  Q e x p k i  are the interval active and interval reactive powers expected at bus k, while  P e x p k d  and  Q e x p k d  are the expected active and reactive power used in the deterministic power flow.
P e x p k i = P e x p k d P c k d + P c k i
Q e x p k i = Q e x p k d Q c k d + Q c k i
The interval mismatches for initializing magnitudes and angles of bus voltages are calculated as the difference between the interval powers and their deterministic values calculated at each bus. Then, the interval increments of the modules and angles of the bus voltages are calculated using the Jacobian matrix of the last iteration of the deterministic power flow. Thus, the interval magnitudes and angles of the bus voltages can be calculated using their deterministic results and the interval increments as detailed in [57].
Step 3: The interval active and reactive powers are calculated using Equations (5) and (6), where  V k i ˙  and  V m i ˙  are the magnitudes of the interval voltages at buses k and m, respectively.  G i j  and  B i j  are the real and imaginary components of the system admittance matrix ( Y i j ) and  θ k i  and  θ m i  are the interval phase angles of the voltages at buses k and m, respectively. Finally,  P c a l c k i  and  Q c a l c k i  are the interval values of the active and reactive powers.
P c a l c k i = V k i ˙ 2 G k k + V k i ˙ m ϵ Ω k k V m i ˙ G k m c o s θ k i θ m i + B k m sin θ k i θ m i
Q c a l c k i = V k i ˙ 2 B k k + V k i ˙ m ϵ Ω k k V m i ˙ B k m c o s θ k i θ m i + G k m sin θ k i θ m i
Then, the interval mismatches are recalculated using Equations (7) and (8).
Δ P k i = P e x p k i P c a l c k i
Δ Q k i = Q e x p k i Q c a l c k i
Step 4: Apply the Krawczyk operator and update the interval values of voltage angles and magnitudes using Equation (9).
θ i V i h + 1 = θ i V i K x h , X h
Step 5: Check whether the largest variation in the radius of the variables between iterations is smaller than the specified tolerance, if yes, the interval power flow has converged, if not, return to step 4.
Step 6: Compute the desired interval electric output quantities,  L i ( V k , θ k , V m , θ m )  using Equations (10) and (11).
Δ L i = L d θ k X p + L d V k Y p + L d θ m Z p + L d V m W p [ Δ P i Δ Q i ] t
L i = L d + Δ L i
where  L i  and  L d  are, respectively, the interval and deterministic values of the output magnitude.  X p Y p Z p , and  W p  are the rows of the inverse Jacobian matrix evaluated after the convergence of the deterministic FP as shown in Equation (12).
[ Δ θ k i Δ V k i Δ θ m i Δ V m i ] t = [ X p Y p Z p W p ] t [ Δ P i Δ Q i ] t

3. Optimal Location and Sizing of DG

Consumers and investors are those who usually decide on the location of new DG units based on their needs and resources. Nonetheless, the utility may provide incentives for the location of DG in specific places of the network that would render other benefits. The approach presented in this paper considers that the utility may decide on the location and sizing of new DG units within the network. The objective function is the minimization of total active power losses with the sitting and sizing of DG as decision variables. Furthermore, an IPF is used to introduce uncertainty to the problem.

3.1. Mathematical Modeling

DG allocation in non-optimal locations may cause technical issues such as an increment in power losses and problems in protection coordination schemes. Therefore, it is important to have an approach capable of selecting the best locations and preferable sizing of DG in distribution networks [58].
Metaheuristic techniques are optimization tools better suited to handle combinatorial, non-convex optimization problems than classic approximation approaches [59]. They have been strongly developed since their introduction in the 1980s. Over time, these methods have added in their implementation strategies that helped them to escape from local optimal solutions in complex solution spaces, especially in methods that use one or more neighborhood structures as a metric for defining moves and transitions in solutions.
Algorithms based on metaheuristics evaluate several potential or candidate solutions and perform a series of operations so that new solutions are produced and re-evaluated. These algorithms are designed in such a way that candidate solutions can be stored and manipulated by the different operators used in the chosen technique. There are three types of metaheuristics according to the method they use to arrive at a solution, namely, constructive, local search, and population-based approaches. However, these are not mutually exclusive and, therefore, metaheuristics may combine different features to form hybrid methods [60,61,62]. In this paper, the SOS and PSO techniques were chosen to test the proposed approach; nonetheless, any other metaheuristic technique may be implemented. It is worth mentioning that the main purpose of the paper is not a comparison of metaheuristic techniques, but rather the inclusion of uncertainties within the optimal siting and sizing of DG in distribution networks by means of an IPF approach.

3.2. Symbiotic Organism Search

Living organisms rarely live in isolation, due to their dependence on other organisms for their sustenance for survival. This dependency-based relationship is called symbiosis. The most common symbiotic relationships between two organisms in nature are mutualism, commensalism, and parasitism. In mutualism both organisms benefit, in commensalism one organism is benefited and the other one is not affected, and in parasitism one of the organisms benefits, and the other one is actually harmed.
The symbiotic organism search (SOS) metaheuristic simulates symbiotic relationships with paired organisms, which are used to search for the fittest organism to solve the optimization problem. As in other population-based metaheuristics, the SOS algorithm iteratively uses a population of candidate solutions from promising areas to scan the solution space in search of the global optimum value. Initial organisms are randomly generated in the solution space. Each organism represents a candidate solution to the optimization problem with an associated fitness value, which reflects the degree to which the organism is adapted to the desired objective [50].
In the SOS algorithm, generating a new solution is a process governed by mimicking the biological interactions between two organisms in an ecosystem. Each organism interacts with the others in three phases that resemble the real world: mutualism, commensalism, and parasitism. With each generation, the organisms undergo interactions in each of the symbiotic phases until a predetermined stopping criterion is reached [50].

3.2.1. Mutualistic Phase

An example of mutualism that benefits each participating organism is the relationship between the woodpecker and the buffalo. The woodpecker feeds on insects present on the buffalo’s body, while the buffalo benefits by removing these parasites. In SOS,  X i  is the organism occupying the i-th position in the ecosystem. Another organism,  X j  is randomly chosen from the population to interact with  X i . New solution candidates  X i n e w  and  X j n e w  are calculated based on the mutualistic relationship which is modeled as indicated in Equations (13)–(15).
X i n e w = X i + r a n d ( 0 , 1 ) × ( X b e s t Mutual _ Vector × B F 1 )
X j n e w = X j + r a n d ( 0 , 1 ) × ( X b e s t Mutual _ Vector × B F 2 )
Mutual _ Vector = X i + X j 2
In Equations (13) and (14), the term  r a n d ( 0 , 1 )  is a random number between 0 and 1, while  B F 1  and  B F 2  represent how much the organisms will benefit from the mutual relationship. In this case,  B F 1  and  B F 2  vary randomly between 1 and 2.
The mutual vector indicated in Equation (15) represents the mutualism relationship between organisms  X i  and  X j . The expression  ( X b e s t Mutual _ Vector × B F )  reflects the mutual effort to achieve a good outcome and ensure survival. The term  X b e s t  represents the highest degree of fitness in the ecosystem, so using it can help increase the degree of fitness of each of the organisms. Initially, the organisms are updated only if their new fitness values are better than those prior to the organisms’ interaction.

3.2.2. Commensalism Phase

An example of commensalism is the relationship between the shark and the remora. The remora attaches itself to the shark and feeds on the leftover food that the shark leaves behind, while the shark is not affected by the presence and activities of the remora. In this phase, an organism denoted as  X j  is randomly chosen to relate to organism  X i . Thus, organism  X i  benefits from this interaction, while organism  X j  suffers no change. A new candidate  X i n e w  is calculated in this interaction which is mathematically modeled by Equation (16). Then, organism  X i  is updated only if the new organism is fitter than the previous one.
X i n e w = X i + r a n d ( 1 , 1 ) ( X b e s t X j )
In Equation (16), the term  r a n d ( 1 , 1 )  represents a random number between −1 and 1, while the expression  ( X b e s t X j )  reflects the benefit that organism  X i  is receiving from organism  X j , increasing its survival advantages in an ecosystem together with the current fittest organism  X b e s t .

3.2.3. Parasitism Phase

An example of parasitism is the relationship between the chopper and other host birds. The chopper lays its eggs in the nest of another bird, which benefits the chopper because it transfers the job of caring for the chicks to the bird that had its nest invaded, and harms the host because the resources of its offspring will be consumed by the parasite, which can even kill the offspring of the host.
In the parasitism step, organism  X i  takes over the role of a parasite by creating an artificial organism called “Parasite Vector”. The parasite vector is created by duplicating vector  X i  and randomly changing one of its dimensions. An organism  X j  is randomly chosen and serves as a host for the parasite vector. The latter tries to replace  X j  organism, so both organisms are evaluated, and if the parasite vector is fitter, it will eliminate  X j  organism and assume its position in the ecosystem. The SOS metaheuristic was implemented according to the algorithm described below.
Step 1: Generate an ecosystem of size M, with N-dimensional organisms.
Step 2: Evaluate the fitness of each of the organisms in the ecosystem through IPF, where the fitness function  f i t i  of each of the organisms is the value of the interval technical losses obtained by IPF. Penalties are applied if a violation in voltage levels is observed. The penalty consists of adding a high value to the objective function that is wanted to minimize. Voltage level violations can be analyzed using Equation (17), where  V k i ˙  indicates the range of acceptable voltage values per unit. After that, the fittest organism,  X b e s t  is identified.
V k i ˙ [ 0.95 ; 1.05 ]
All organisms in the ecosystem will undergo all the interactions described in phases 1, 2, and 3 at each iteration. The organism undergoing the interaction will be called  X i .
Phase 1—Mutualism:
Step 3: Select an organism  X j  where  j i .
Step 4: The mutual vector is calculated by means of Equation (15).
Step 5: New organisms coming from mutualism are calculated using Equations (13) and (14), with the parameters  B F 1  and  B F 2  randomly drawn between 1 and 2.
Step 6: The new organisms are evaluated and if they are fitter than the original organisms, they replace them; otherwise, they are discarded.
Phase 2—Commensalism
Step 7: An organism  X j  is selected where  j i .
Step 8: A new organism is generated through commensalism according to Equation (16).
Step 9: The new organism is evaluated and if it is fitter than the original  X i  organism, it replaces it; otherwise, it is discarded.
Phase 3—Parasitism
Step 10: An organism  X j  is selected where  j i .
Step 11: A parasite vector is created by randomly changing one dimension of organism  X i .
Step 12: The new parasite organism is evaluated and if it is fitter than organism  X j  it replaces it; otherwise, it is discarded.
Step 13: It is checked if i is the last position in the population (if the whole population has passed through the three symbiotic stages), if yes, go to the next step, if no,  i = i + 1  and return to Step 2.
Step 14: It is checked if the stopping criterion has been reached, if not, return to Step 2. If yes, the search for the optimal solution is complete.
Figure 1 presents the flowchart of the implemented SOS algorithm.

3.3. Particle Swarm Optimization

PSO is a metaheuristic technique proposed in [63] and is based on the representation of the movement of organisms in a bird flock or fish school. PSO has shown to be effective for solving several electrical engineering problems such as optimal sitting and sizing of DG [55,56], ancillary service optimization [64], optimal sizing of photovoltaic systems [65], and optimal protection coordination [66], among many others. PSO has also been hybridized to improve its performance, as reported in [67,68,69]. According to the literature, the motion of birds in a flock or schools of fish, which are commonly named swarms, can be based on five principles [63]:
  • The proximity principle, where individuals in the population must be able to move around in a search space.
  • The quality principle, in which individuals must be able to respond to quality factors in the environment.
  • The principle of diverse responses, in which individuals should not be bound to a restricted path.
  • The stability principle, in which individuals should not change their behavior whenever environmental conditions change.
  • The adaptability principle, in which individuals must be able to change their behavior when it is no longer convenient.
The PSO algorithm starts from the basis of a swarm that is artificially generated in the search space of solutions. Each element of this swarm is called a particle and behaves as a point without mass and volume in an N-dimensional space. Each particle is a candidate solution to the optimization problem. The iterations of these particles are governed by the five principles mentioned above.
PSO is a search algorithm for optimal solutions based on the behavior of a swarm, in which each particle can memorize its own “optimal” value, the “optimal” value of all other particles in the swarm, and its own velocity, and has a fitness value that shows how close that particle is to the optimal position. Thus, in each generation, the particles’ information is combined to generate new velocities for each of the dimensions, which are used to generate a new particle [70]. The flowchart of the implemented PSO is presented in Figure 2.
To describe the PSO algorithm, a swarm of size M in an N-dimensional space is assumed, where the position of each particle i is  X i = ( x i 1 , x i 2 , , x i d , . , x i N ) , its velocity is  V i = ( v i 1 , v i 2 , , v i d , , v i N ) , and the optimal position that each particle has assumed, that is, the position at which that particle has assumed the best fitness value. This position is given by  P i = ( p i 1 , p i 2 , , p i d , , p i N ) . The optimal swarm position, which corresponds to the best of the optimal positions of the particles, is given by  P g = ( p g i 1 , p g i 2 , , p g i d , , p g i N ) .
Thus, new velocities and positions for the particles can be calculated as:
v i , t + 1 d = ω × v i , t d + c 1 × r a n d ( 0 , 1 ) × ( p i , t d x i , t d ) + c 2 × r a n d ( 0 , 1 ) × ( p g , t x i , t d )
x i , t + 1 d = x i , t d + v i , t + 1 d
From Equation (18) it can be seen that the first term  ω v i , t d  represents the influence of the particle’s prior velocity, so parameter  ω  is called the inertial weight. The second term,  c 1 × r a n d ( 0 , 1 ) × ( p i , t d x i , t d ) , represents the distance between the particle and its own optimal position. In this case, parameter  c 1  is called the cognitive learning factor. The third term,  c 2 × r a n d ( 0 , 1 ) × ( p g , t x i , t d ) , represents the distance between the particle and the optimal swarm position. It represents the fact that particles can share information and cooperate, moving according to the experience of the whole swarm. In this term, the parameter  c 2  is called the social learning factor.
The velocities and positions of each of the particles are iteratively updated. The local and global maxima,  P g  and  P i , are recalculated at each iteration until a predetermined stopping criterion is reached.
The PSO metaheuristic was implemented according to the algorithm described below.
Step 1: Generate a swarm of size M, with N-dimensional particles and set iteration zero with variable  i t e r = 0 .
Step 2: Evaluate each particle in the swarm via IPF, where the fitness function  f i t i  of each of the organisms is the value of interval technical losses obtained by IPF and penalties are applied if a violation in voltage levels is observed. Violations in voltage levels can be analyzed using Equation (17). After that, define the particle with the best position towards the solution  P g  in the swarm and also the best position of each of the particles  P i .
Step 3: A new velocity for the particles is calculated using Equation (18). In this case, learning factors c1 = 2 and c2 = 1.49445, which are commonly used, are adopted [70]. Furthermore, an inertia factor  ω  is used given by Equation (20).
ω = ω m a x ω m a x ω m i n m a x _ i t e r × i t e r
In Equation (20),  m a x _ i t e r  is the maximum number of iterations,  i t e r  is the value of the current iteration, and  ω m a x  = 0.9 and  ω m i n  = 0.4 are adopted as the maximum and minimum values for the inertial weight, respectively.
Step 4: Check if the stopping criterion  i t e r = m a x _ i t e r  has been reached, if not,  i t e r = i t e r + 1  and return to Step 2, if yes, the optimal solution has been found.

3.4. Metrics

From the information returned by the IPF intervals, it is necessary to define a criterion to evaluate which is the best value of the objective function; that is, the smallest value of the system power losses. In this work, the following comparison methodologies were implemented:

3.4.1. Comparison of the Midpoints

In this method, the interval  X = [ x 1 , x 2 ]  is considered to be smaller than the interval  Y = [ y 1 , y 2 ]  if their midpoints, obtained from Equation (1), follow the following constraint given by Equation (21).
X < Y M i d ( X ) < M i d ( Y )

3.4.2. Evaluation Using the Interval Measurement Function

In this method, the interval  X = [ x 1 , x 2 ]  is considered to be smaller than the interval  Y = [ y 1 , y 2 ]  if the interval measure function given by Equation (22) returns a value greater than zero, according to Equations (23)–(26) [71], where  m x  and  m y  represent the midpoint of the intervals X and Y, respectively, while  r x  and  r y  represent the radius of the intervals, respectively.
μ ( X , Y ) = m y m x + 2 · s g n ( m y m x ) , if r y + r x = 0 m y m x r y + r x + s g n ( m y m x ) , if m x m y and r y + r x 0 r y r x m a x { r y , r x } , if m x = m y and r y + r x 0
m x = M i d ( X )
r x = R a d i u s ( X )
m y = M i d ( Y )
r y = R a d i u s ( Y )
In this case,  s g n  is given by Equation (27).
s g n ( x ) = 1 if x < 0 0 if x = 0 1 if x > 0
Thus, X < Y can be mathematically described by Equation (28).
X < Y μ ( X , Y ) > 0

4. Tests and Results

In this paper, the IPF approach described and validated in [57] was implemented. For the validation of the IPF, a deterministic power flow using the Newton–Rapshon method with sparse matrix treatment was developed. After running the deterministic power flow, the IPF is also processed as described in Section 2 (please see reference [57] for more details). Medium voltage distribution systems with radial topologies were simulated with load variations of  5 4  between iterations. All tests were carried out considering a single-phase equivalent of the distribution test systems. The stopping criterion considers a maximum variation in the radius of the variables of  10 4  between iterations. All computational simulations were implemented in Python®.
After the validation of the IPF results, the optimal allocation and sizing of DG in the test networks were developed in Python® based on the described metaheuristics: SOS [50] and PSO [63]. The DG units allocated by these algorithms were also simulated with 5% variations in their power. Regarding generation technologies, for the sake of simplicity, the proposed approach does not model the intrinsic generation variability of renewable energy resources; instead, an expected generation value is considered for the buses selected for DG allocation. The maximum numbers of iterations were 100 and 250 for the IEEE 33-bus and IEEE 69-bus test systems, respectively.

4.1. Results with the IEEE 33-Bus Test System

4.1.1. System Data

The data on the IEEE 33-bus test system can be found in [72]. This system has a total active load of 3715 kW and a total reactive load of 2300 kvar. The base values considered are 100 kVA and 12.66 kV. The maximum and minimum values of the voltage magnitude, obtained from the deterministic power flow, are 1.0 p.u and 0.909 p.u, respectively.
The optimal allocation and sizing of DG performed on the IEEE 33-bus test systems was carried out considering buses 7, 10, 13, 26, 31, and 33 as eligible for DG allocation. This consideration is based on the fact that in real distribution systems, not all network buses are eligible for DG allocation due to constraints related to the availability of resources as well as land use limitations. It was also considered the restriction that the maximum penetration rate of DG in this system is 30%, that is, the total generation cannot exceed 30% of the grid active load. Thus, for this system, the maximum generation to be allocated is 1114.5 kW.

4.1.2. SOS Applied to the IEEE 33-Bus Test System

The SOS algorithm was processed for the DG allocation, with 100 iterations in an ecosystem of 20 organisms, and had its active power losses evaluated under the two metrics presented in Section 3.4. The amounts of power allocated to each of the buses are presented in Table 2. To validate the results of the SOS algorithm, a convergence analysis was performed using the results of 10 optimization simulations using the midpoint comparison metric. Table 3 presents the results of these simulations. Note that all of them led to similar values, evidencing the convergence of the algorithm.
The proposed methodology only considers the results that meet voltage limits. Figure 3 allows us to verify this criterion. In this case, variables  V i n f _ m 1  and  V s u p _ m 1  are, respectively, the lower and upper bounds of the magnitude of the interval voltages on each of the buses in the solution obtained by the SOS algorithm using the midpoint comparison metric. On the other hand,  V i n f _ m 2  and  V s u p _ m 2  are, respectively, the lower and upper bounds of the magnitude of the interval voltage on each of the buses in the solution obtained by the SOS algorithm, using the evaluation metric via the interval measure function. Finally,  V i n f  and  V s u p  are, respectively, the lower and upper bounds of the magnitude of the interval voltage at each of the buses obtained by IPF without DG. Note that there is an important improvement in the voltage profile when DG is included in the network.
Table 4 presents the values of interval losses obtained by the SOS algorithm for the two metrics used. Figure 4 and Figure 5 show the convergence of the midpoints of the interval losses and their diameters for the two metrics. From Figure 4, it is observed that the convergence of the algorithm is obtained around iteration 55, which shows that the search process aiming at minimum losses is adequate. On the other hand, when the radius of the loss range of the best solution obtained at each iteration is evaluated, it is observed from Figure 5 that the diameter increases over the course of the simulation, i.e., for both metrics the uncertainty of the result increases over the course of the simulation. Also note that the simulation performed by the comparison metric using the interval measure function results in a smaller diameter.
Therefore, according to the results observed, it can be concluded that the SOS algorithm is robust and efficient in solving the DG allocation and sizing problem when both metrics are considered.

4.1.3. PSO Applied to the IEEE 33-Bus Test System

The PSO algorithm considers 100 iterations under a swarm of 20 particles. The active power losses were evaluated, during the process, by means of the two metrics presented in Section 3.4. The amount of power allocated to each bus is indicated in Table 5.
To validate the results of the PSO algorithm, a convergence analysis was performed using the results of 10 optimization simulations using the midpoint comparison metric. Table 6 presents the results of these simulations. Note that all of the values are similar, evidencing the convergence of the algorithm.
According to Equation (17), the methodology only considers the results that satisfy the pre-established voltage limits. Figure 6 evidences the enforcement of this criterion.
In Figure 6, the variables  V i n f _ m 1  and  V s u p _ m 1  are, respectively, the lower and upper bounds of the magnitude of the interval voltage at each bus of the system in the solution obtained by the PSO algorithm using the midpoint comparison metric. The variables  V i n f _ m 2  and  V s u p _ m 2  are, respectively, the lower and upper bounds of the magnitude of the interval voltage at each bus of the system in the solution obtained by the PSO algorithm using the evaluation metric via the interval measure function. Finally, the lines corresponding to  V i n f  and  V s u p  are, respectively, the lower and upper bounds of the magnitude of the interval voltage at each bus obtained by the IPF without DG.
Table 7 presents the values of interval losses obtained by the PSO algorithm for the two metrics used. Figure 7 and Figure 8 show the convergence of the midpoints of the interval losses and their diameters for the two metrics.
From Figure 7, it is observed that the convergence of the algorithm is obtained around iteration 9. This shows that the search process aiming at minimizing power losses is successful since the midpoint decreases with each iteration. However, it is observed that for both metrics the algorithm is stuck at a local minimum. On the other hand, when evaluating the radius of the loss range of the best solution obtained at each iteration, it is observed in Figure 8 that the radius increases over the course of the simulation, that is, for both metrics, the uncertainty of the result increases over the course of the simulation.
Therefore, regarding the convergence of the PSO algorithm, it can be concluded that it worked correctly when seeking to solve the DG allocation and sizing problem when considering the two metrics, but was not able to escape from local optimal solutions.

4.1.4. Power Loss Comparison

Figure 9 illustrates the values of losses considering midpoint comparison and evaluation with interval measurement function for both PSO and SOS algorithms. It is observed that SOS presents better performance than PSO, where the SOS midpoint is 75.42 MW while the PSO midpoint is 84.92 MW, approximately.
In summary, from the results presented in Section 4.1.2, Section 4.1.3 and Section 4.1.4 it is clear that the SOS algorithm performed better when compared to the PSO. From Figure 7 it is evident that the PSO algorithm was stuck at a local minimum for both metrics.

4.2. Results with the IEEE 69-Bus Test system

4.2.1. System Data

The data for this system can be found in [73]. The base values considered are 100 kVA and 12.66 kV. The 69-bus system has a total active load of 3802.19 kW and a total reactive load of 2694.60 kvar. The maximum and minimum values of the voltage magnitude, obtained from a deterministic power flow, are 1.0 p.u and 0.909 p.u, respectively.
The proposed methodology for optimal allocation and sizing of DG considers the following candidate buses: 10, 18, 27, 40, 49, 54, 62, and 68. It is also considered that the maximum DG penetration rate is 20%, which corresponds to 760.44 kW.

4.2.2. SOS Applied to the IEEE 69-Bus Test System

The SOS algorithm was simulated considering 250 iterations in an ecosystem of 20 organisms. The active power losses are evaluated under the two metrics presented in Section 3.4. The amounts of power allocated as a result of the implemented algorithm at each bus are presented in Table 8.
According to Equation (17), the methodology only considers the results that meet the pre-established voltage limits. Figure 10 allows us to verify this condition.
In Figure 10, the lines corresponding to  V i n f _ m 1  and  V s u p _ m 1  are, respectively, the lower and upper bounds of the magnitude of the interval voltage on the system buses in the solution obtained by the SOS algorithm, using the midpoint comparison metric. The lines of variables  V i n f _ m 2  and  V s u p _ m 2  are, respectively, the lower and upper bounds of the magnitude of the interval voltage on the system buses in the solution obtained by the SOS algorithm, using the evaluation metric via the interval measure function.  V i n f  and  V s u p  indicate, respectively, the lower and upper bounds of the magnitude of the interval voltage obtained by the IPF without DG.
The interval losses obtained by the SOS algorithm are presented in Table 9. Figure 11 and Figure 12 show the convergence of the midpoints of the interval losses and their diameters for the two metrics. From Figure 11, it can be inferred that the convergence of the algorithm is obtained around iteration 92. On the other hand, when the diameter of the loss range of the best solution obtained at each iteration is evaluated, it can be seen in Figure 12 that the radius varies during the simulation. Also note that simulation performed by the comparison metric using the interval measure function results in a smaller diameter.
Therefore, regarding the convergence of the SOS algorithm, it was observed that the algorithm is robust and efficient in solving the DG allocation and sizing problem when both metrics are considered.

4.2.3. PSO Applied to the IEEE 69-Bus Test System

The PSO algorithm was run considering 250 iterations on a swarm of 20 particles. The active power losses are evaluated under the two metrics presented in Section 3.4. The amounts of power allocated to the candidate buses are presented in Table 10.
As already mentioned, the proposed approach only considers results that meet voltage limits. This can be verified in Figure 13, where  V i n f _ m 1  and  V s u p _ m 1  are the lower and upper bounds of the magnitude of the interval voltage on the system buses in the solution obtained by the PSO algorithm, using the midpoint comparison metric. The lines of variables  V i n f _ m 2  and  V s u p _ m 2  are, respectively, the lower and upper bounds of the magnitude of the interval voltage on the system buses in the solution obtained by the PSO algorithm, using the evaluation metric via the interval measure function. Finally, the variables  V i n f  and  V s u p  are, respectively, the lower and upper bounds of the magnitude of the interval voltage on the system buses obtained by IPF without DG.
Table 11 shows the interval loss values obtained by the PSO algorithm for the two metrics used. Figure 14 and Figure 15 show the convergence of the midpoints of the interval losses and their diameters for the two metrics. From Figure 14, it is observed that the convergence of the algorithm is obtained around iteration 9 and, it shows that the search process aiming at minimum losses is adequate since the midpoint decreases with each iteration. However, it is clearly observed that for both metrics the algorithm is stuck in a local minimum. On the other hand, when evaluating the radius of the loss range of the best solution obtained at each iteration, it is observed from Figure 15 that the diameter increases over the course of the simulation, i.e., for both metrics the uncertainty of the result increases over the course of the simulation. Also note that the simulation performed by the comparison metric using the interval measure function results in a smaller diameter.
Therefore, regarding the convergence of the PSO algorithm, it was observed that the algorithm worked correctly in trying to solve the DG allocation and sizing when considering the two metrics, but it was not able to escape from local optimal solutions. Furthermore, the solution at which the algorithm was stuck is not even able to regulate the feeder voltage within the desired limits.

4.2.4. Power Loss Comparison

Figure 16 presents the power losses considering midpoint comparison and evaluation with interval measurement function for the metaheuristics under study. Note that the midpoint of lower power losses obtained with the SOS and PSO approaches approximately correspond to 104.57 MW and 169.63 MW, respectively.
In summary, from the results presented in Section 4.2.2, Section 4.2.3 and Section 4.2.4 it is observed that the SOS algorithm performed better when compared to the PSO. Additionally, according to Figure 14, it is evident that the PSO algorithm was stuck at a local minimum for both metrics.

5. Conclusions

This paper presented an optimal location and sizing approach of DG in distribution networks using PSO and SOS metaheuristics; furthermore, the proposed approach implements an interval power flow that adds uncertainty to the optimization problem. The results obtained in the two distribution test systems evidenced the superiority of the SOS methodology over PSO. With the parameters and modeling used in this work, the SOS metaheuristic obtained better results in all simulations. One of the advantages of SOS over PSO lies in the fact that SOS does not require configuration or fine-tuning of parameters. The implemented methodologies were evaluated using two different metrics, the midpoint comparison and the interval evaluation function. From the results obtained, both metrics present equivalent results; however, the interval evaluation function results show a smaller diameter.
In particular, for the IEEE 69-bus network, it was observed that the PSO algorithm did not perform well, not even being able to propose a solution capable of meeting the constraints regarding the system voltage levels. Nonetheless, in cases where the optimal solutions involved the allocation of DG in all candidate buses, the PSO algorithm showed reasonable results. The PSO algorithm did not have good convergence when considering the simultaneity of the DG allocation and sizing problems; however, according to the results obtained, adequate convergence can be expected when only the sizing problem is addressed.
The SOS algorithm using IPF showed to be a robust methodology capable of adding uncertainty to the optimal location and sizing problem of DG units in distribution systems, which is becoming increasingly common in modern power distribution systems.

Author Contributions

Conceptualization, W.C.N., L.P.G.N. and J.M.L.-L.; Data curation, W.C.N. and L.P.G.N.; Formal analysis, W.C.N. and L.P.G.N.; Funding acquisition, J.M.L.-L.; Investigation, W.C.N., L.P.G.N. and J.M.L.-L.; Methodology, W.C.N.; Project administration, W.C.N., L.P.G.N. and J.M.L.-L.; Resources, W.C.N., L.P.G.N. and J.M.L.-L.; Software, W.C.N.; Supervision, L.P.G.N.; Validation, L.P.G.N.; Visualization, L.P.G.N.; Writing—original draft, W.C.N.; Writing—review and editing, W.C.N., L.P.G.N. and J.M.L.-L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Colombia Scientific Program within the framework of the call Ecosistema Cientifíco (Contract no. FP44842-218-2018). The authors also acknowledge the financial support of the Brazilian National Council for Scientific and Technological Development (CNPq), grant 408898/2021-6, FAPEMIG-APQ-03609-17, CAPES-Brazil under Grant 001, and INERGE.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is available via authors by email.

Acknowledgments

The authors gratefully acknowledge the financial support provided by the Colombia Scientific Program within the framework of the call Ecosistema Cientifíco (Contract no. FP44842-218-2018). The authors also acknowledge the support of the Brazilian National Council for Scientific and Technological Development (CNPq), grant 408898/2021-6 and FAPEMIG-APQ-03609-17, CAPES-Brazil under Grant 001, and INERGE.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bai, X.; Qu, L.; Qiao, W. Robust AC Optimal Power Flow for Power Networks With Wind Power Generation. IEEE Trans. Power Syst. 2016, 31, 4163–4164. [Google Scholar] [CrossRef]
  2. Lopez-lezama, J.M.; Murillo-Sanchez, C.; Zuluaga, L.; Gutierrez-Gomez, J. A Contingency-Based Security-Constrained Optimal Power Flow Model for Revealing the Marginal Cost of a Blackout Risk-Equalizing Policy in the Colombian Electricity Market. In Proceedings of the 2006 IEEE/PES Transmission Distribution Conference and Exposition: Latin America, Caracas, Venezuela, 15–18 August 2006; pp. 1–6. [Google Scholar] [CrossRef]
  3. Bian, J.; Wang, H.; Wang, L.; Li, G.; Wang, Z. Probabilistic optimal power flow of an AC/DC system with a multiport current flow controller. CSEE J. Power Energy Syst. 2021, 7, 744–752. [Google Scholar] [CrossRef]
  4. Amiri, H.; Mohammadi, M.; Karimi, M.; Rastegar, M.; Rostami, A. Probabilistic Load Flow Based on Parameterized Probability-Boxes for Systems with Insufficient Information. IEEE Access 2021, 9, 161038–161045. [Google Scholar] [CrossRef]
  5. Villada Duque, F.; López Lezama, J.M.; Muñoz Galeano, N. Effects of Incentives for Renewable Energy in Colombia. Ing. Univ. 2017, 21, 257–272. [Google Scholar] [CrossRef] [Green Version]
  6. Strezoski, L.; Padullaparti, H.; Ding, F.; Baggu, M. Integration of Utility Distributed Energy Resource Management System and Aggregators for Evolving Distribution System Operators. J. Mod. Power Syst. Clean Energy 2022, 10, 277–285. [Google Scholar] [CrossRef]
  7. Ryu, J.; Kim, J. Virtual Power Plant Operation Strategy Under Uncertainty with Demand Response Resources in Electricity Markets. IEEE Access 2022, 10, 62763–62771. [Google Scholar] [CrossRef]
  8. Rezaei, E.; Dagdougui, H.; Ojand, K. Hierarchical Distributed Energy Management Framework for Multiple Greenhouses Considering Demand Response. IEEE Trans. Sustain. Energy 2023, 14, 453–464. [Google Scholar] [CrossRef]
  9. Primadianto, A.; Lu, C. A Review on Distribution System State Estimation. IEEE Trans. Power Syst. 2017, 32, 3875–3883. [Google Scholar] [CrossRef]
  10. de Oliveira, L.W.; Seta, F.d.S.; de Oliveira, E.J. Optimal reconfiguration of distribution systems with representation of uncertainties through interval analysis. Int. J. Electr. Power Energy Syst. 2016, 83, 382–391. [Google Scholar] [CrossRef]
  11. Pereira, L.; da Costa, V.; Rosa, A. Interval arithmetic in current injection power flow analysis. Int. J. Electr. Power Energy Syst. 2012, 43, 1106–1113. [Google Scholar] [CrossRef]
  12. Araujo, B.M.C.; da Costa, V.M. New Developments in the Interval Current Injection Power Flow Formulation. IEEE Lat. Am. Trans. 2018, 16, 1969–1976. [Google Scholar] [CrossRef]
  13. Seta, F.D.S.; de Oliveira, L.W.; de Oliveira, E.J. Distribution System Planning with Representation of Uncertainties Based on Interval Analysis. J. Control Autom. Electr. Syst. 2020, 31, 494–510. [Google Scholar] [CrossRef]
  14. Mori, H.; Yuihara, A. Calculation of multiple power flow solutions with the Krawczyk method. In Proceedings of the 1999 IEEE International Symposium on Circuits and Systems (ISCAS), Orlando, FL, USA, 30 May–2 June 1999; Volume 5, pp. 94–97. [Google Scholar] [CrossRef]
  15. Rasmussen, T.B.; Yang, G.; Nielsen, A.H. Interval estimation of voltage magnitude in radial distribution feeder with minimal data acquisition requirements. Int. J. Electr. Power Energy Syst. 2019, 113, 281–287. [Google Scholar] [CrossRef]
  16. Lin, X.; Shu, T.; Tang, J.; Ponci, F.; Monti, A.; Li, W. Application of Joint Raw Moments-Based Probabilistic Power Flow Analysis for Hybrid Power Systems. IEEE Trans. Power Syst. 2022, 37, 1399–1412. [Google Scholar] [CrossRef]
  17. Li, Y.; Wan, C.; Chen, D.; Song, Y. Nonparametric Probabilistic Optimal Power Flow. IEEE Trans. Power Syst. 2022, 37, 2758–2770. [Google Scholar] [CrossRef]
  18. Shu, T.; Lin, X.; Peng, S.; Du, X.; Chen, H.; Li, F.; Tang, J.; Li, W. Probabilistic Power Flow Analysis for Hybrid HVAC and LCC-VSC HVDC System. IEEE Access 2019, 7, 142038–142052. [Google Scholar] [CrossRef]
  19. Peng, S.; Lin, X.; Tang, J.; Xie, K.; Ponci, F.; Monti, A.; Li, W. Probabilistic Power Flow of AC and DC Hybrid Grids With Addressing Boundary Issue of Correlated Uncertainty Sources. IEEE Trans. Sustain. Energy 2022, 13, 1607–1619. [Google Scholar] [CrossRef]
  20. Fu, C.; Xu, Y.; Yang, Y.; Lu, K.; Gu, F.; Ball, A. Response analysis of an accelerating unbalanced rotating system with both random and interval variables. J. Sound Vib. 2020, 466, 115047. [Google Scholar] [CrossRef]
  21. Cheng, W.; Cheng, R.; Shi, J.; Zhang, C.; Sun, G.; Hua, D. Interval Power Flow Analysis Considering Interval Output of Wind Farms through Affine Arithmetic and Optimizing-Scenarios Method. Energies 2018, 11, 3176. [Google Scholar] [CrossRef] [Green Version]
  22. Luo, L.; Gu, W.; Wang, Y.; Chen, C. An Affine Arithmetic-Based Power Flow Algorithm Considering the Regional Control of Unscheduled Power Fluctuation. Energies 2017, 10, 1794. [Google Scholar] [CrossRef] [Green Version]
  23. Wang, Y.; Wu, Z.; Dou, X.; Hu, M.; Xu, Y. Interval power flow analysis via multi-stage affine arithmetic for unbalanced distribution network. Electr. Power Syst. Res. 2017, 142, 1–8. [Google Scholar] [CrossRef]
  24. Wang, C.; Liu, D.; Tang, F.; Liu, C. A clustering-based analytical method for hybrid probabilistic and interval power flow. Int. J. Electr. Power Energy Syst. 2021, 126, 106605. [Google Scholar] [CrossRef]
  25. Zhang, C.; Chen, H.; Shi, K.; Qiu, M.; Hua, D.; Ngan, H. An Interval Power Flow Analysis Through Optimizing-Scenarios Method. IEEE Trans. Smart Grid 2018, 9, 5217–5226. [Google Scholar] [CrossRef]
  26. Liao, X.; Liu, K.; Zhang, Y.; Wang, K.; Qin, L. Interval method for uncertain power flow analysis based on Taylor inclusion function. Transm. Distrib. IET Gener. 2017, 11, 1270–1278. [Google Scholar] [CrossRef]
  27. Liu, B.; Huang, Q.; Zhao, J.; Hu, W. A Computational Attractive Interval Power Flow Approach With Correlated Uncertain Power Injections. IEEE Trans. Power Syst. 2020, 35, 825–828. [Google Scholar] [CrossRef]
  28. Zhang, C.; Chen, H.; Ngan, H.; Yang, P.; Hua, D. A Mixed Interval Power Flow Analysis Under Rectangular and Polar Coordinate System. IEEE Trans. Power Syst. 2017, 32, 1422–1429. [Google Scholar] [CrossRef]
  29. Tang, K.; Dong, S.; Zhu, C.; Song, Y. Affine Arithmetic-Based Coordinated Interval Power Flow of Integrated Transmission and Distribution Networks. IEEE Trans. Smart Grid 2020, 11, 4116–4132. [Google Scholar] [CrossRef]
  30. Çelik, D.; Meral, M.E. Current control based power management strategy for distributed power generation system. Control Eng. Pract. 2019, 82, 72–85. [Google Scholar] [CrossRef]
  31. Chen, G.; Lewis, F.L.; Feng, E.N.; Song, Y. Distributed Optimal Active Power Control of Multiple Generation Systems. IEEE Trans. Ind. Electron. 2015, 62, 7079–7090. [Google Scholar] [CrossRef]
  32. Machado, L.F.M.; Santo, S.G.D.; Junior, G.M.; Itiki, R.; Manjrekar, M.D. Multi-Source Distributed Energy Resources Management System Based on Pattern Search Optimal Solution Using Nonlinearized Power Flow Constraints. IEEE Access 2021, 9, 30374–30385. [Google Scholar] [CrossRef]
  33. Ackermann, T.; Andersson, G.; Söder, L. Distributed generation: A definition. Electr. Power Syst. Res. 2001, 57, 195–204. [Google Scholar] [CrossRef]
  34. Saldarriaga-Zuluaga, S.D.; López-Lezama, J.M.; Muñoz-Galeano, N. An Approach for Optimal Coordination of Over-Current Relays in Microgrids with Distributed Generation. Electronics 2020, 9, 1740. [Google Scholar] [CrossRef]
  35. Rajagopalan, A.; Swaminathan, D.; Alharbi, M.; Sengan, S.; Montoya, O.D.; El-Shafai, W.; Fouda, M.M.; Aly, M.H. Modernized Planning of Smart Grid Based on Distributed Power Generations and Energy Storage Systems Using Soft Computing Methods. Energies 2022, 15, 8889. [Google Scholar] [CrossRef]
  36. Gallego Pareja, L.A.; López-Lezama, J.M.; Gómez Carmona, O. A Mixed-Integer Linear Programming Model for the Simultaneous Optimal Distribution Network Reconfiguration and Optimal Placement of Distributed Generation. Energies 2022, 15, 3063. [Google Scholar] [CrossRef]
  37. Pérez Posada, A.F.; Villegas, J.G.; López-Lezama, J.M. A Scatter Search Heuristic for the Optimal Location, Sizing and Contract Pricing of Distributed Generation in Electric Distribution Systems. Energies 2017, 10, 1449. [Google Scholar] [CrossRef] [Green Version]
  38. Zhang, C.; Li, J.; Zhang, Y.J.; Xu, Z. Optimal Location Planning of Renewable Distributed Generation Units in Distribution Networks: An Analytical Approach. IEEE Trans. Power Syst. 2018, 33, 2742–2753. [Google Scholar] [CrossRef]
  39. Afraz, A.; Malekinezhad, F.; Shenava, S.; Jlili, A. Optimal sizing and sitting in radial standard system using pso. Am. J. Sci. Res. 2012, 67, 50–58. [Google Scholar]
  40. Mareddy, P.; Reddy, V.; Vyza, U. Optimal DG placement for minimum real power loss in radial distribution systems using PSO. J. Theor. Appl. Inf. Technol. 2010, 13, 107–116. [Google Scholar]
  41. Sedighizadeh, M.; Rezazadeh, A. Using Genetic Algorithm for Distributed Generation Allocation to Reduce Losses and Improve Voltage Profile. Int. J. Comput. Syst. Eng. 2008, 2, 50–55. [Google Scholar]
  42. Liu, K.Y.; Sheng, W.; Liu, Y.; Meng, X.; Liu, Y. Optimal sitting and sizing of DGs in distribution system considering time sequence characteristics of loads and DGs. Int. J. Electr. Power Energy Syst. 2015, 69, 430–440. [Google Scholar] [CrossRef]
  43. Bhadoria, V.S.; Pal, N.S.; Shrivastava, V.; Jaiswal, S.P. Reliability Improvement of Distribution System by Optimal Sitting and Sizing of Disperse Generation. Int. J. Reliab. Qual. Saf. Eng. 2017, 24, 1740006. [Google Scholar] [CrossRef]
  44. Abou El-Ela, A.; Allam, S.; Shatla, M. Maximal optimal benefits of distributed generation using genetic algorithms. Electr. Power Syst. Res. 2010, 80, 869–877. [Google Scholar] [CrossRef]
  45. Moradi, M.; Abedini, M. A novel method for optimal DG units capacity and location in Microgrids. Int. J. Electr. Power Energy Syst. 2016, 75, 236–244. [Google Scholar] [CrossRef]
  46. Moradi, M.H.; Abedini, M.; Hosseinian, S.M. A Combination of Evolutionary Algorithm and Game Theory for Optimal Location and Operation of DG from DG Owner Standpoints. IEEE Trans. Smart Grid 2016, 7, 608–616. [Google Scholar] [CrossRef]
  47. Ameli, A.; Farrokhifard, M.R.; Davari-nejad, E.; Oraee, H.; Haghifam, M.R. Profit-Based DG Planning Considering Environmental and Operational Issues: A Multiobjective Approach. IEEE Syst. J. 2017, 11, 1959–1970. [Google Scholar] [CrossRef]
  48. Ali, A.; Keerio, M.U.; Laghari, J.A. Optimal Site and Size of Distributed Generation Allocation in Radial Distribution Network Using Multi-objective Optimization. J. Mod. Power Syst. Clean Energy 2021, 9, 404–415. [Google Scholar] [CrossRef]
  49. Purlu, M.; Turkay, B.E. Optimal Allocation of Renewable Distributed Generations Using Heuristic Methods to Minimize Annual Energy Losses and Voltage Deviation Index. IEEE Access 2022, 10, 21455–21474. [Google Scholar] [CrossRef]
  50. Cheng, M.Y.; Prayogo, D. Symbiotic Organisms Search: A new metaheuristic optimization algorithm. Comput. Struct. 2014, 139, 98–112. [Google Scholar] [CrossRef]
  51. Sanjay, R.; Jayabarathi, T.; Raghunathan, T.; Ramesh, V.; Mithulananthan, N. Optimal Allocation of Distributed Generation Using Hybrid Grey Wolf Optimizer. IEEE Access 2017, 5, 14807–14818. [Google Scholar] [CrossRef]
  52. Farh, H.M.H.; Al-Shaalan, A.M.; Eltamaly, A.M.; Al-Shamma’A, A.A. A Novel Crow Search Algorithm Auto-Drive PSO for Optimal Allocation and Sizing of Renewable Distributed Generation. IEEE Access 2020, 8, 27807–27820. [Google Scholar] [CrossRef]
  53. Eid, A.; Kamel, S.; Korashy, A.; Khurshaid, T. An Enhanced Artificial Ecosystem-Based Optimization for Optimal Allocation of Multiple Distributed Generations. IEEE Access 2020, 8, 178493–178513. [Google Scholar] [CrossRef]
  54. Nayeripour, M.; Mahboubi-Moghaddam, E.; Aghaei, J.; Azizi-Vahed, A. Multi-objective placement and sizing of DGs in distribution networks ensuring transient stability using hybrid evolutionary algorithm. Renew. Sustain. Energy Rev. 2013, 25, 759–767. [Google Scholar] [CrossRef]
  55. da Rosa, W.; Gerez, C.; Belati, E. Optimal Distributed Generation Allocating Using Particle Swarm Optimization and Linearized AC Load Flow. IEEE Lat. Am. Trans. 2018, 16, 2665–2670. [Google Scholar] [CrossRef]
  56. Arif, S.M.; Hussain, A.; Lie, T.T.; Ahsan, S.M.; Khan, H.A. Analytical Hybrid Particle Swarm Optimization Algorithm for Optimal Siting and Sizing of Distributed Generation in Smart Grid. J. Mod. Power Syst. Clean Energy 2020, 8, 1221–1230. [Google Scholar] [CrossRef]
  57. Nogueira, W.C.; Garcés Negrete, L.P.; López-Lezama, J.M. Interval Load Flow for Uncertainty Consideration in Power Systems Analysis. Energies 2021, 14, 642. [Google Scholar] [CrossRef]
  58. Alinejad-Beromi, Y.; Sedighizadeh, M.; Sadighi, M. A particle swarm optimization for sitting and sizing of Distributed Generation in distribution network to improve voltage profile and reduce THD and losses. In Proceedings of the 2008 43rd International Universities Power Engineering Conference, Padova, Italy, 1–4 September 2008; pp. 1–5. [Google Scholar] [CrossRef]
  59. Saldarriaga-Zuluaga, S.D.; López-Lezama, J.M.; Muñoz-Galeano, N. Adaptive protection coordination scheme in microgrids using directional over-current relays with non-standard characteristics. Heliyon 2021, 7, e06665. [Google Scholar] [CrossRef] [PubMed]
  60. Gendreau, M.; Potvin, J.Y. Handbook of Metaheuristics; Springer: New York, NY, USA, 2010; Volume 2. [Google Scholar]
  61. Radosavljević, J.; Arsić, N.; Milovanović, M.; Ktena, A. Optimal Placement and Sizing of Renewable Distributed Generation Using Hybrid Metaheuristic Algorithm. J. Mod. Power Syst. Clean Energy 2020, 8, 499–510. [Google Scholar] [CrossRef]
  62. Ellahi, M.; Abbas, G. A Hybrid Metaheuristic Approach for the Solution of Renewables-Incorporated Economic Dispatch Problems. IEEE Access 2020, 8, 127608–127621. [Google Scholar] [CrossRef]
  63. Kennedy, J.; Eberhart, R. Particle swarm optimization. In Proceedings of the ICNN’95—International Conference on Neural Networks, Perth, WA, Australia, 27 November–1 December 1995; Volume 4, pp. 1942–1948. [Google Scholar] [CrossRef]
  64. Dos Santos Alonso, A.M.; Pereira Junior, B.R.; Brandao, D.I.; Marafao, F.P. Optimized Exploitation of Ancillary Services: Compensation of Reactive, Unbalance and Harmonic Currents Based on Particle Swarm Optimization. IEEE Lat. Am. Trans. 2021, 19, 314–325. [Google Scholar] [CrossRef]
  65. Silva de Souza, J.; Percy Molina, Y.; Silva de Araujo, C.; Pereira de Farias, W.; Santos de Araujo, I. Modified Particle Swarm Optimization Algorithm for Sizing Photovoltaic System. IEEE Lat. Am. Trans. 2017, 15, 283–289. [Google Scholar] [CrossRef]
  66. Ghanbari, M.; Gandomkar, M.; Nikoukar, J. Protection Coordination of Bidirectional Overcurrent Relays Using Developed Particle Swarm Optimization Approach Considering Distribution Generation Penetration and Fault Current Limiter Placement. IEEE Can. J. Electr. Comput. Eng. 2021, 44, 143–155. [Google Scholar] [CrossRef]
  67. Prashant; Sarwar, M.; Siddiqui, A.S.; Ghoneim, S.S.M.; Mahmoud, K.; Darwish, M.M.F. Effective Transmission Congestion Management via Optimal DG Capacity Using Hybrid Swarm Optimization for Contemporary Power System Operations. IEEE Access 2022, 10, 71091–71106. [Google Scholar] [CrossRef]
  68. Zhen, L.; Liu, Y.; Dongsheng, W.; Wei, Z. Parameter Estimation of Software Reliability Model and Prediction Based on Hybrid Wolf Pack Algorithm and Particle Swarm Optimization. IEEE Access 2020, 8, 29354–29369. [Google Scholar] [CrossRef]
  69. Dabhi, D.; Pandya, K. Enhanced Velocity Differential Evolutionary Particle Swarm Optimization for Optimal Scheduling of a Distributed Energy Resources With Uncertain Scenarios. IEEE Access 2020, 8, 27001–27017. [Google Scholar] [CrossRef]
  70. Wang, D.; Tan, D.; Liu, L. Particle swarm optimization algorithm: An overview. Soft Comput. 2018, 22, 387–408. [Google Scholar] [CrossRef]
  71. Alolyan, I. A new method for comparing closed intervals. Aust. J. Math. Anal. Appl. 2011, 8, 1–6. [Google Scholar]
  72. Alam, A.; Gupta, A.; Bindal, P.; Siddiqui, A.; Zaid, M. Power Loss Minimization in a Radial Distribution System with Distributed Generation. In Proceedings of the 2018 International Conference on Power, Energy, Control and Transmission Systems (ICPECTS), Chennai, India, 22–23 February 2018; pp. 21–25. [Google Scholar] [CrossRef]
  73. Baran, M.E.; Wu, F.F. Optimal capacitor placement on radial distribution systems. IEEE Trans. Power Deliv. 1989, 4, 725–734. [Google Scholar] [CrossRef]
Figure 1. Flowchart of the implemented SOS algorithm.
Figure 1. Flowchart of the implemented SOS algorithm.
Sustainability 15 05171 g001
Figure 2. Flowchart of the implemented PSO.
Figure 2. Flowchart of the implemented PSO.
Sustainability 15 05171 g002
Figure 3. Voltage profile of the IEEE 33-bus test system—SOS.
Figure 3. Voltage profile of the IEEE 33-bus test system—SOS.
Sustainability 15 05171 g003
Figure 4. Convergence of SOS algorithm for the IEEE 33-bus test system—midpoints.
Figure 4. Convergence of SOS algorithm for the IEEE 33-bus test system—midpoints.
Sustainability 15 05171 g004
Figure 5. Convergence of SOS algorithm for the IEEE 33-bus test system—diameters.
Figure 5. Convergence of SOS algorithm for the IEEE 33-bus test system—diameters.
Sustainability 15 05171 g005
Figure 6. Voltage profile of the IEEE 33-bus test system—PSO.
Figure 6. Voltage profile of the IEEE 33-bus test system—PSO.
Sustainability 15 05171 g006
Figure 7. Convergence of the PSO algorithm for the IEEE 33-bus test system—midpoints.
Figure 7. Convergence of the PSO algorithm for the IEEE 33-bus test system—midpoints.
Sustainability 15 05171 g007
Figure 8. Convergence of the PSO algorithm for the IEEE 33-bus test system—diameters.
Figure 8. Convergence of the PSO algorithm for the IEEE 33-bus test system—diameters.
Sustainability 15 05171 g008
Figure 9. Convergence of PSO and SOS—33-bus test system.
Figure 9. Convergence of PSO and SOS—33-bus test system.
Sustainability 15 05171 g009
Figure 10. Voltage profile of the IEEE 69-bus test system—SOS.
Figure 10. Voltage profile of the IEEE 69-bus test system—SOS.
Sustainability 15 05171 g010
Figure 11. Convergence of the SOS algorithm for the IEEE 69-bus—midpoints.
Figure 11. Convergence of the SOS algorithm for the IEEE 69-bus—midpoints.
Sustainability 15 05171 g011
Figure 12. Convergence of the SOS algorithm for the IEEE 69-bus—diameters.
Figure 12. Convergence of the SOS algorithm for the IEEE 69-bus—diameters.
Sustainability 15 05171 g012
Figure 13. Voltage profile of the IEEE 69-bus test system—PSO.
Figure 13. Voltage profile of the IEEE 69-bus test system—PSO.
Sustainability 15 05171 g013
Figure 14. Convergence of the PSO algorithm for the IEEE 69-bus test system—midpoints.
Figure 14. Convergence of the PSO algorithm for the IEEE 69-bus test system—midpoints.
Sustainability 15 05171 g014
Figure 15. Convergence of the PSO algorithm for the IEEE 69-bus test system—diameters.
Figure 15. Convergence of the PSO algorithm for the IEEE 69-bus test system—diameters.
Sustainability 15 05171 g015
Figure 16. Power loss comparison between PSO and SOS—69-bus test system.
Figure 16. Power loss comparison between PSO and SOS—69-bus test system.
Sustainability 15 05171 g016
Table 1. Main characteristics of CPF, IPF, and PPF.
Table 1. Main characteristics of CPF, IPF, and PPF.
Conventional PFInterval PFProbabilistic PF
Nature of input variablesDeterministic variablesInterval variablesProbability distribution
functions
Mathematical modelPower balance equations
with conventional mathematical
operations
Power balance equations
with interval mathematical
operations
Conventional power balance
equations considering the
realizations of random variables
Solution approachNewton–RaphsonKrawczyk methodA pre-determined number
of simulations, each
running a Newton–Raphson
Obtained solutionUnic solutionInterval form (interval
ranges of the output data
are not necessarily the
same as those of the
input data)
Probability distribution
function (not necessarily
the same as the one of the
input data)
Table 2. Active power of DG allocated per bus by the SOS algorithm in the IEEE 33-bus test network.
Table 2. Active power of DG allocated per bus by the SOS algorithm in the IEEE 33-bus test network.
Active Power of DG Allocated (kW)
Bus71013263133Total
Midpoint Comparison00528.20304.8281.31114.4
Evaluation with Interval Measurement Function0145.6410.470527.330.631114.3
Table 3. Convergence analysis of the SOS algorithm on the IEEE 33-bus test network.
Table 3. Convergence analysis of the SOS algorithm on the IEEE 33-bus test network.
SOS
SimulationLower Limit of Losses (kW)Upper Limit of Losses (kW)
172.9577.88
272.9177.84
371.1576.88
473.2077.73
572.9677.69
673.0878.01
772.8377.86
874.0078.83
972.9377.76
1072.9177.93
Table 4. Interval losses obtained by the SOS algorithm for the IEEE 33-bus network.
Table 4. Interval losses obtained by the SOS algorithm for the IEEE 33-bus network.
MetricLower Limit of Losses
(kW)
Upper Limit of Losses
(kW)
Midpoint comparison72.9577.88
Evaluation with Interval Measurement Function73.1477.91
Table 5. Active power of DG allocated per bus by the PSO algorithm in the IEEE 33-bus test network.
Table 5. Active power of DG allocated per bus by the PSO algorithm in the IEEE 33-bus test network.
Active Power of DG Allocated (kW)
Bus71013263133Total
Midpoint Comparison169.6186.4172.9184.1203.0191.61107.8
Evaluation with Interval Measurement Function159.1185.5169.2183.0198.6214.71110.3
Table 6. Convergence analysis of the PSO algorithm on the IEEE 33-bus network.
Table 6. Convergence analysis of the PSO algorithm on the IEEE 33-bus network.
PSO
SimulationLower Limit of Losses (kW)Upper Limit of Losses (kW)
182.7287.10
282.5287.42
381.7487.16
482.9288.54
582.5987.11
683.0887.79
782.5188.22
881.6187.52
984.6289.34
1082.7788.58
Table 7. Interval losses obtained by the PSO algorithm for the IEEE 33-bus network.
Table 7. Interval losses obtained by the PSO algorithm for the IEEE 33-bus network.
MetricLower Limit of Losses
(kW)
Upper Limit of Losses
(kW)
Midpoint Comparison82.7287.10
Evaluation with Interval Measurement Function82.2586.65
Table 8. Active power of DG allocated per bus by the SOS algorithm in the IEEE 69-bus network.
Table 8. Active power of DG allocated per bus by the SOS algorithm in the IEEE 69-bus network.
Allocated Distributed Generation Active Power (kW)
Bus1018274049546368Total
Midpoint Comparison0018.72021.320720.40760.44
Evaluation with Interval Measurement Function13.211.102.9162.630670.60760.44
Table 9. Interval losses obtained by the SOS algorithm for the IEEE 69-bus test network.
Table 9. Interval losses obtained by the SOS algorithm for the IEEE 69-bus test network.
MetricLower Limit of Losses
(kW)
Upper Limit of Losses
(kW)
Midpoint Comparison79.7129.5
Evaluation Using the Interval Measurement Function85.7135.3
Table 10. Active power of DG allocated per bar by the SOS algorithm in the IEEE 69-bus network.
Table 10. Active power of DG allocated per bar by the SOS algorithm in the IEEE 69-bus network.
Allocated Distributed Generation Active Power (kW)
Bus1018274049546368Total
Midpoint Comparison69.9118.195.4752.986102.8115.188.5728.7
Evaluation with Interval Measurement Function92.997.582.3100.575.797.3114.185.3745.7
Table 11. Interval losses obtained by the PSO algorithm for the IEEE 69-bus network.
Table 11. Interval losses obtained by the PSO algorithm for the IEEE 69-bus network.
MetricLower Limit of Losses
(kW)
Upper Limit of Losses
(kW)
Midpoint Comparison145.56193.70
Evaluation with Interval Measurement Function148.77196.30
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Nogueira, W.C.; Garcés Negrete, L.P.; López-Lezama, J.M. Optimal Allocation and Sizing of Distributed Generation Using Interval Power Flow. Sustainability 2023, 15, 5171. https://doi.org/10.3390/su15065171

AMA Style

Nogueira WC, Garcés Negrete LP, López-Lezama JM. Optimal Allocation and Sizing of Distributed Generation Using Interval Power Flow. Sustainability. 2023; 15(6):5171. https://doi.org/10.3390/su15065171

Chicago/Turabian Style

Nogueira, Wallisson C., Lina P. Garcés Negrete, and Jesús M. López-Lezama. 2023. "Optimal Allocation and Sizing of Distributed Generation Using Interval Power Flow" Sustainability 15, no. 6: 5171. https://doi.org/10.3390/su15065171

APA Style

Nogueira, W. C., Garcés Negrete, L. P., & López-Lezama, J. M. (2023). Optimal Allocation and Sizing of Distributed Generation Using Interval Power Flow. Sustainability, 15(6), 5171. https://doi.org/10.3390/su15065171

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

Article Metrics

Back to TopTop