Next Article in Journal
Platelet Aggregation Inhibitors and Anticoagulants Delay Surgery for Hip Fractures
Next Article in Special Issue
Decentralized Voltage Optimization Based on the Auxiliary Problem Principle in Distribution Networks with DERs
Previous Article in Journal
Polyphenols Extracted from Chinese Hickory (Carya cathayensis) Promote Apoptosis and Inhibit Proliferation through the p53-Dependent Intrinsic and HIF-1α-VEGF Pathways in Ovarian Cancer Cells
Previous Article in Special Issue
Energy Commitment for a Power System Supplied by Multiple Energy Carriers System using Following Optimization Algorithm
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hybrid GA-SOCP Approach for Placement and Sizing of Distributed Generators in DC Networks

by
Oscar Danilo Montoya
1,2,*,
Walter Gil-González
3 and
Luis Fernando Grisales-Noreña
3
1
Facultad de Ingeniería, Universidad Distrital Francisco José de Caldas, Bogotá 11021, DC, Colombia
2
Laboratorio Inteligente de Energía, Universidad Tecnológica de Bolívar, Cartagena 131001, Colombia
3
Grupo GIIEN, Facultad de Ingeniería, Institución Universitaria Pascual Bravo, Campus Robledo, Medellín 050036, Colombia
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(23), 8616; https://doi.org/10.3390/app10238616
Submission received: 13 September 2020 / Revised: 16 October 2020 / Accepted: 26 October 2020 / Published: 2 December 2020
(This article belongs to the Special Issue Control, Optimization and Planning of Power Distribution Systems)

Abstract

:
This research addresses the problem of the optimal location and sizing distributed generators (DGs) in direct current (DC) distribution networks from the combinatorial optimization. It is proposed a master–slave optimization approach in order to solve the problems of placement and location of DGs, respectively. The master stage applies to the classical Chu & Beasley genetic algorithm (GA), while the slave stage resolves a second-order cone programming reformulation of the optimal power flow problem for DC grids. This master–slave approach generates a hybrid optimization approach, named GA-SOCP. The main advantage of optimal dimensioning of DGs via SOCP is that this method makes part of the exact mathematical optimization that guarantees the possibility of finding the global optimal solution due to the solution space’s convex structure, which is a clear improvement regarding classical metaheuristic optimization methodologies. Numerical comparisons with hybrid and exact optimization approaches reported in the literature demonstrate the proposed hybrid GA-SOCP approach’s effectiveness and robustness to achieve the global optimal solution. Two test feeders compose of 21 and 69 nodes that can locate three distributed generators are considered. All of the computational validations have been carried out in the MATLAB software and the CVX tool for convex optimization.

1. Introduction

Recently, the study and analysis of direct current (DC) networks have increased, due to the development of power electronics, advances in renewable generation, and the advantages over alternative current (AC) networks [1,2]. Initially, the DC networks began to operate at high-level voltages in the power transmission for long-distance between systems, since it allows for reducing transmission power losses, enhancing grid stability systems, easily incorporating off-shore wind systems, and managing power flow both direction simply [3,4,5]. Consecutively, the integration of the DC networks has been widened to medium- and low-voltage level applications, generating concepts, such as the microgrids [6,7].
The main advantages of DC systems as compared to conventional AC systems can be summarized, as follows: (i) fewer power converter devices, (ii) no requirement of synchronizing generators, (iii) better voltage profiles, (iv) less power losses, and (v) do not require frequency and reactive power analysis and controls [8,9,10].
Two kinds of analysis are used to study DC networks. The first analysis is dynamic, being performed under the time domain, and integrates concepts as primary and secondary controllers [11,12]. The second analysis is static, which computes all of the state variables under stationary conditions. This analysis is known as tertiary control [13], which addresses power flow analysis [2,14], optimal power flow [15], economic dispatch strategies [16], and voltage stability analysis [17,18]. Furthermore, these analyses are integrated with other problems, such as the placement and sizing of multiple distributed generators (DGs)  [19], as well as placement and operation of battery energy storage systems in DC networks [20].
In this study, we focus on the problem of the optimal placement and sizing of DGs in DC networks. Although, in the specialized literature, there are few approaches in this area. In [21], an analytical methodology for sizing energy storage systems in DC networks was proposed. This methodology evaluated the profiles of currents that were injected by energy storage systems in order to minimize power losses. In [22], a feasibility study for the integration of photovoltaic systems in a standalone building was carried out. This study considered strong weather scenarios for off-grid systems and variable electricity prices for the grid-connected system. In [23], a probabilistic optimization methodology for estimating uncertainties in loads, plug-in electric vehicles, and renewable generation in DC distribution networks was analyzed. In [19], a hybrid strategy for optimal location and dimensioning of DGs in DC grids was provided, which mixed a hyperplanes search method (HSM) with semidefinite programming (SDP) model. The HSM relaxed the binary variables, while the SDP computed the optimal power flow. In [24], a similar strategy using a sequential quadratic optimal power flow was presented. Both of the techniques ensure the optimal global solution only for a relaxed problem, where the complexity lies in the number of candidate solutions to recover the binary variables. In [25], a tutorial in the general algebraic modeling system (GAMS) for the optimal placement and dimensioning of DGs in DC grids were presented. In [26], the optimal placement and sizing of photovoltaic generators in DC networks were carried out, when considering an environmental-economic dispatch as an objective function. However, the approach presented in [26] was solved in GAMS, which does not guarantee the optimal global solution, since the problem is nonlinear non-convex mixed-integer. In [27], three metaheuristic methods for optimal location and dimensioning of DGs in DC systems were studied. A classical genetic algorithm was employed in the locations stage, which assigns the location of the DGs. While the black-hole, continuous genetic algorithm, and particle swarm metaheuristic methods were used to dimension the DGs. The main problem with these methods is that they have multiple tuning parameters that significantly affect their performance. Additionally, they do not ensure the optimal global solution to the problem. In [28], a methodology for optimal sizing of multiple devices, such as wind turbines, photovoltaic sources, and batteries in hybrid AC/DC grids was provided, which considered the life cycle cost and the life span of the renewable energy resources in the objective function.
After the revision mentioned above of the state-of-the-art, this research’s main contributions can be summarized, as follows:
🗸
The reformulation of the problem of the optimal placement and sizing of DGs in DC networks that allows transforming the original mixed-integer nonlinear programming (MINLP) model into a mixed-integer second-order cone programming (MI-SOCP) one.
🗸
The solution of the proposed MI-SOCP model with a hybrid master-slave optimization methodology, where the master stage is guided by a genetic algorithm to determine the best locations for the DGs, while the slave stage allows for solving the optimal power flow problem by using the SOCP formulation via interior point methods to determine the best sizes of the DGs.
🗸
The computational results demonstrate that the hybrid GA-SOCP approach has an efficiency higher than 70% to find the global optimal solution in both test feeders analyzed. In addition, the proposed approach presents better solutions than commercial solvers available in GAMS and classical metaheuristics such as genetic algorithms and particle swarm optimizer, and derived approaches
It is worth mentioning that the master–slave hybrid GA-SOCP approach has not previously been proposed in specialized literature for addressing the problem of the optimal placement and sizing of DGs in electrical networks, even if these are under AC or DC paradigms of operation. This situation was identified as a gap in the scientific literature that this research tries to fill.
The remainder of this document is organized, as follows: Section 2 presents the original MINLP formulation of the problem of the optimal placement and dimensioning of DGs in DC networks by highlighting its nonlinearities and non-convexities. Section 3 presents the proposed hybrid optimization method, which is composed of a conic reformulation problem of the optimal power flow problem in DC networks and the classical GA. Section 4 exposes the main characteristics of the test feeders, which are composed of 21 and 69 nodes with radial structures, where the possibility of installing three DGs is considered. Section 5 presents all of the computational validations of the proposed hybrid GA-SOCP approach and their comparisons with nonlinear solvers available in GAMS and master–slave optimization approaches that are based on genetic algorithms and particle swarm optimization methods and its variants. Section 6 shows the concluding remarks derived from this work and some possible future researches in the area of distribution system optimization.

2. MINLP Model

The problem of the optimal placement and sizing of DGs in DC networks corresponds to an MINLP model since: (i) the problem of sizing of these power sources is modeled with binary variables, and (ii) the problem of the location of these DGs is modeled with continuous variables (power generation and voltage profiles). The optimization objective in this MINLP formulation corresponds to the minimization of the grid power losses produced by the current circulation in all the networks’ branches. The mathematical structure of this objective function is presented in Equation (1).
min p loss = i N j N G i j v i v j ,
where p loss is the value of the objective function regarding power losses in all of the branches of the network, v i and v j are the voltage values at nodes i and j, respectively, and G i j is the conductance value taken from the admittance nodal matrix that defines the physical interconnection among all the nodes of the DC network. Observe that N is the set that contains all of the nodes of the network.
Regarding model constraints, these are related with operational conditions of the DC distribution networks. These constraints are presented below:
p i s + p i g d p i d = j N G i j v i v j , i N
p i j = g i j v i v i v j , i , j L
i N p i g d α i N p i d ,
i N y i N d g a v a ,
y i p g d , min p i g d y i p g d , max i N
p i j min p i j p i j max i , j L
v min v i v max i N
y i 0 , 1 i N
where p i s p i g d , p i d , and p i j represent the power generation in the voltage controlled source (slack nodes), power generation in distributed sources, power consumption in demand nodes, and power flow in distribution lines, respectively. α is a parameter that defines the percentage of power penetration of renewable generation as function of the total demand; y i is the binary variable associate with the location ( y i = 1 ) or no ( y i = 0 ) of a distributed generator at node i; N d g a v a is the maximum number of distributed generators available for location in the DC network, p g d , min and p g d , max are the minimum and maximum power generation capabilities allowed for each distributed generator connected at node i. p i j min and p i j max are the lower and upper bounds for the power flow in line that connects nodes i and j, respectively; and, v min and v max are the minimum and maximum voltage bounds in all nodes of the DC distribution network, respectively.
The complete interpretation of the MINLP model defined from (1) to (9) is as follows: Equation (1) defines the objective function of the problem regarding the minimization of the total grid power losses in the whole DC network; Equation (2) is the set of power balance equations applied at each node of the network; Equation (3) defines the amount of power that flows at each line of the DC grid. Expression (3) determines the maximum penetration of the distributed generation in the whole DC distribution network; Inequality constraint (4) bounds the number of generators available for location in the DC grid; inequalities (6), (7), and (8) present the upper and lower bounds for distributed generators, power flow in branches, and voltage variables, respectively. Finally, Expression (9) shows the binary nature of the decision variable regarding the location of distributed generators.
Remark 1.
In the problem of the optimal placement and sizing of distributed generators in DC distribution networks is complex due to two main reasons: (i) the binary nature of the location problem makes necessary integer programming methods to define the best possible location of the distributed generators (master optimization problem); and (ii) the optimal sizing of the distributed generators for each possible location provided by the master stage is a nonlinear non-convex optimization problem due to the product between voltage variables in Expressions (1) to (3).
To solve the MINLP model (1)–(9), this research proposes a hybrid optimization approach based on genetic algorithms to address the problem of the optimal location of distributed generators and the second-order cone programming model to solve the problem of the optimal sizing of these devices. The main advantage of using SCOP instead of metaheuristics for optimal power flow analysis is that for each possible combination of distributed generator locations, the solution is the global optima [29], which is not possible with any particle swarm derived method that corresponds to the most classical approach for solving this kind of problems in literature [27].

3. Proposed Hybrid GA-SOCP Optimization Approach

We proposed a hybrid GA-SOCP approach, which works by using a master–slave optimization strategy, where the master stage guides the discrete part of the problem (i.e., defines the location of the DGs), and the slave stage deals with their optimal sizing, in order to address the problem of the optimal location of distributed generators in DC networks in this research.

3.1. Slave Stage: Second-Order Cone Programming

The SOCP is a branch of the convex optimization that allows transforming some nonlinear non-convex optimization problems into conic convex equivalents that guarantee the global optimum of the equivalent model [30]. In the case of the optimal power flow problem, the SOCP equivalent model has the same global optimum that the original nonlinear model, which is the main advantage of this technique in power system analysis [29]. Here, we present the SOCP version of the optimal power flow problem that addresses distributed generators’ optimal sizing. For this purpose, let us suppose that all of the binary variables regarding the distributed generators’ location, i.e., y i has been fixed by the master stage, which will be presented in this study. In this sense, the SOCP deals with the Equality constraints (2) and (3), which are the only continuous nonlinear non-convex constraints.
The convexification of the optimal power flow problem is made using the following auxiliary variable:
z i j = v i v j ,
note that, from (10), we know that z i i = v i 2 and z j j = v j , which implies that if we elevate both sides of (10) at the square, then we have:
z i j 2 = z i i z j j .
Now, in order to obtain an equivalent representation of (11), we use the hyperbolic relation of the product between two variables as follows:
z i j 2 = z i i z j j , = 1 4 z i i + z j j 2 1 4 z i i z j j 2 , ( 2 z i j ) 2 + z i i z j j 2 = z i i + z j j 2 , 2 z i j z i i z j j = z i i + z j j .
Note that Equation (12) is the conic equivalent of (11), since, to reach this result, only algebraic manipulations have been made; in addition, it is still non-convex due to the equality symbol; nevertheless, it can be relaxed with a lower-equal symbol as recommended in [13,31,32], which becomes (12) into a conic convex constraint, as presented below:
2 z i j z i i z j j z i i + z j j .
Remark 2.
To obtain a SOCP equivalent of the optimal power flow problem, we can recur to substitute (11) into Expressions (1) to (3), which produce equivalent mixed-integer second-order cone programming (i.e., MISOCP) model to address the problem of the optimal placement and sizing of distributed generators in DC distribution networks.
The complete MISOCP model is presented below:
min p loss = i N j N G i j z i j ,
p i s + p i g d p i d = j N G i j z i j , i N
p i j = g i j z i i z i j , i , j L
2 z i j z i i z j j z i i + z j j , i , j N
Expressions ( 4 ) to ( 9 ) ,
Remark 3.
Observe that the slave optimization stage solves the proposed MI-SOCP model once all the binary variables have been assigned on the master stage with the main advantage that the SOCP structure of the power flow model guarantees the global optimum finding [30].

3.2. Master Stage: Genetic Algorithm

The genetic algorithm is one of the most classical optimization methodologies from the family of combinatorial nature-inspired optimization methods [33]. The GA algorithm deals with complex optimization problems that involve integer (binary) variables, using three evolution rules: selection, recombination, and mutation [34]. Here, we select the GA optimizer to deal with the problem of the location of the distributed sources in DC networks using an integer combination that allows for minimizing the infeasibilities caused by classical binary codifications as recommended in [34]. The complete optimization structure of the proposed GA is presented below.

3.2.1. Initial Population

The most of nature-inspired optimization algorithms work with a list of potential solutions named “population”. This list has the complete information regarding the solution of the problem under study. Here, we consider an integer population to determine the location of the DGs in the Dc network. This population W t (being t the iterative counter) takes the following form:
W t = w 11 w 12 w 1 a w 21 w 22 w 2 a w b 1 w b 2 w b a
where a is the number of variables of interest and b is the number of the individuals in the population.
Note that the component w k m in the initial population matrix is calculated, as follows:
w k m = round 1 + rand ( · n 1 , k = 1 , 2 , , b , m = 1 , 2 , , n .
where round ( · ) determines the integer part of the argument and rand ( · ) is a random number generated with a normal Gaussian distribution in the interval 0 , 1 . It is worth remembering that n is the number of nodes and it corresponds to the cardinality of the set N , i.e., n = | N | .
Remark 4.
An essential stage during the creation of the initial population is that it is mandatory to guarantee that each potential individual W i t is different from each other W j t ; j i . This condition is known in genetic algorithms terminology as the diversity verification.

3.2.2. Fitness Function Evaluation

It is needed to define the fitness function in order to evaluate the initial population. The fitness function in metaheuristic optimization is the “adaptation” of the objective function to deal with possible infeasibilities in the solution space [35,36]. However, in the proposed MI-SOCP model, we guarantee that each individual is feasible, since the SOCP model is correctly constrained by the optimization model (14)–(18). In this sense, our proposed methodology’s fitness function is equal to the objective function that is defined in Equation (14).

3.2.3. Selection

The selection (tournament) is the procedure where a subset of parents is chosen from the initial population to obtain a set of potential solutions (offsprings). Here, we select the most straightforward selection method that consists of selecting two individuals (a different one to each other) from the population, named as the parents.

3.2.4. Recombination

The recombination is a heuristic procedure where the genetic information of the parents is crossed to generate offsprings. To present this procedure, let us consider two parents that are provided by the selection stage as follows:
W i t = 1 0 1 1 0 1 1 0 0 0
W j t = 0 0 1 0 0 0 1 0 1 1 .
Note that, in this example, the each parent has 10 genes ( n = 10 ). To make the recombination it is generate a random number between 1 and n 1 . Consider that this number is 6. With this number is make the recombination procedure by adding the first 6 genes of the individual W i t with the last 4 genes of W j t ; in addition, the first 6 genes of the individual W j t are added to the last 4 genes of W i t . The resulting offsprings ( Os i t and Os j t ) are presented below.
Os i t = 1 0 1 1 0 1 1 0 1 1 ,
Os j t = 0 0 1 0 0 0 1 0 0 0 .

3.2.5. Mutation

The mutation is a procedure in which a subset of genes is changed. This procedure modifies the initial information of the parents after the recombination procedure. To make this procedure, we use the mutation with a unique point for each offspring, which implies that two random number between 1 and n 1 is generated. Suppose that, for Os i t , the random number is 3 and for Os j t this number is 9, which produces the following mutations in both offsprings.
Os i t = 1 0 0 1 0 1 1 0 1 1 ,
Os j t = 0 0 1 0 0 0 1 0 1 0 .
Note that both of the offsprings are different from their parents, since both have new information added by the mutation procedure.

3.2.6. Fitness Function Evaluations and Replacement of the Population

Once the offsprings have been created by the selection, recombination, and mutation stages, these are evaluated in the fitness function, i.e., optimization model (14)–(18). Let us suppose that the fitness function of the offspring Os i t is z Os i t and the fitness function of the offspring Os j t is z Os j t . Subsequently, the potential offspring that will be part of the population is selected as the minimum of both fitness function values: (i) select Os i t if z Os i t z Os j t , or (ii) select Os j t otherwise.
In order to decide whether the winner offspring is added to the population, two aspects must be verified: (i) is the fitness function of the winner offspring better than the worst individual in the population? If the answer is positive, then, (ii) is the winner offspring different from each of the individuals in the population? If the answer is also positive, then the winner offspring replace the worst individual in the population.
Note that, if one of the previous questions has a negative response, then a new iteration is started without modification in the population.

3.2.7. Stopping Criteria

In order to finalize the searching procedure of the proposed genetic algorithm, the following aspects are considered:
🗸
If the maximum number of iteration t max is reached, the optimization procedure ends by reporting the best solution obtained in the current population.
🗸
If during τ max consecutive iterations, the best objective function has not improved, then the optimization procedure ends by reporting the best solution obtained in the current population.
Remark 5.
The genetic algorithm that is presented in the subsections mentioned above is known in the specialized literature as the Chu & Beasley genetic algorithm [34]. It does not replace all the population individuals in the same iteration, as is made by the classical genetic algorithm.

3.3. Pseudocode of the Proposed Hybrid Optimization Algorithm

Algorithm 1 summarized the main aspects of the proposed hybrid GA-SOCP method for optimal placement and sizing of distributed generators in DC distribution networks.
Algorithm 1: Proposed master-slave optimization approach for the optimal location and sizing of distributed generators (DGs) in DC distribution networks
1:Inputs:
2:Read the data of the DC network;
3:Define the number of generators available, i.e., N d g a v a
4:Define the percentage of penetration of distributed generation i.e., α
5:Generate the initial population W 0
6:Evaluate the fitness function for each individual in the population (see the slave stage); 
7:Make t = 0
8:while  t t max  do
9: Select two parents from the population; 
10: Make recombinations for generating two offspring individuals;
11: Apply the mutation operator for each offspring population;
12: Evaluate the fitness function of the offspring individuals, i.e., slave stage (see optimization model (14)–(18)); 
13: Select the offspring with the minimum fitness function; 
14: Replace the worst individual in the population with the selected offspring if the offspring has the best objective function and fulfills the diversity criterion;
15:if  τ τ max  then
16:  Select as the solution of the problem, the individual with the best objective function in the current population W t
17:  Return the optimal solution concerning the location and sizing of DGs (see the slave stage);
18:  break;
19:end if 
20: t = t + 1
21:end while 
22:Output:  
23:The best solution is found for the MINLP model;

4. Test Systems

In this section, the main characteristics of the 21- and 69-node test feeders are presented. These grids are mainly used in power system analysis for DC networks, including voltage stability, power flow analysis, and distributed generators’ optimal location [37]. Both of the test feeders are described below.

4.1. 21-Node Test Feeder

The 21-node test feeder is a DC distribution network that is composed of 21 nodes and 20 branches that generates a radial configuration where the voltage-controlled source is connected at node 1. This source supports a voltage profile of 1.00 pu [19]. The complete information regarding constant power consumption and branches are presented in Figure 1 and Table 1, respectively [19].
Note that all of the values in Table 1 are calculated when considering 1 kV and 100 kW as the voltage and power bases, respectively.

4.2. 69-Node Test Feeder

This test feeder is a classical AC distribution network that is employed in the literature, such as shunt devices’ optimal location for power losses minimization [37]. Notwithstanding, in this study, we used its DC adaptation reported in [37]. Figure 2 depicts the 69-node test feeder configuration and Table 2 reports the branch parameters and load information. Furthermore, for simulation purposes, 12.66 kV and 100 kVA are used as the test system’s voltage and power bases.

5. Computational Implementation

The proposed hybrid optimization model, i.e., the GA-SOCP method, has been implemented in the MATLAB software version 2019b on a desktop computer with an INTEL(R) Core(TM) i7-7700 2.8-GHz processor and 16.0 GB of RAM running on a 64-bit version of Microsoft Windows 10 Home. To validate our proposed MI-SOCP formulation solved with the proposed hybrid GA-SOCP approach, we employed different MINLP solvers that were available in the GAMS software (i.e., SBB and CONOPT solvers), and some metaheuristic approaches reported in [37]. In [37], the genetic algorithm has been combined with three continuous methods, as follows: black-hole optimizer (GA-BHO), particle swarm optimization (GA-PSO), and continuous genetic algorithm (GA-CGA). In addition, the following simulating conditions are considered:
🗸
Three distributed generators are available for installation in the DC test feeder; wherein, the case of the 21-node test system, their maximum individual allowed size is 1.5 pu and, in the case of the 69-node test feeder, this bound is 12 pu.
🗸
The maximum power penetration in the 21-node test feeder is 40 % and, in the case of the 69-node test feeder, the maximum penetration is 40 % regarding the total power demand of each one of the DC networks.
Remark 6.
All the genetic algorithms employed as the master stage in the proposed and comparison methods have adjusted using a trial and error procedure, which have been produced as parameters: 10 individuals in the population, 100 iterations, recombination, and mutation rates 50%. In addition, to evaluate statistically speaking these methods, 100 consecutive iterations were made for each method.

5.1. 21-Node Test Feeder

Table 3 reports the numerical comparisons between the proposed hybrid optimization method and the literature approaches for the 21-node test feeder. Note that this problem’s optimal solution corresponds to nodes 9, 12, and 16, with power injections of 0.8350 pu, 1.0258 pu, and 1.4632 pu, respectively.
From Table 3, we can observe that:
🗸
All of the comparison methods, as well as the proposed approach, identifies the best optimal location in nodes 9, 12, and 16. However, due to the aleatory nature of the continuous metaheuristics (i.e., BHO, PSO, and CGA), both of them are stuck local optimums, and only the GA-PSO approach can deal with the global optimal solution.
🗸
The mixed-integer quadratic programming (MIQP) model and the SBB solver (for the exact MINLP model) reach the global optimal solution with the lowest computational effort, which is attributable to the efficiency of the GAMS interface.
🗸
The standard deviation of the proposed GA-SOCP programming is the lowest when compared with metaheuristics; this situation occurs because the proposed approach works with an exact optimization method based on second-order cone programming that has a unique solution for each combination of DGs provided by the GA in the master stage.
Figure 3 presents the histogram after the 100 consecutive evaluations of the proposed methodology for optimal placement and sizing distributed generators in DC distribution networks in order to present the efficiency of the proposed GA-SOCP approach.
Note that the histogram presented in Figure 3 shows that more than 90 solutions (93 exactly) reached by the proposed hybrid GA-SOCP method has the capability of finding the global optimal solution (i.e., 0.0306 pu) regarding power losses, which confirms the efficiency and robustness of the studied methodology for optimal placement and sizing DGs in DC distribution networks.
Remark 7.
We can ensure that the solution finds by the proposed GA-SOCP approach is the optimal one, since all the possibilities (i.e., 1140 options) for the location of DGs in the 21-node test feeder have been evaluated comprehensively via nested loops, which have confirmed that nodes 9, 12 and 16 are the best possible location for these DGs reaching a final power losses of 0.0306 pu. It is worth mentioning that this exhaustive evaluation takes about 800 s.
Figure 4, which presents the percentage of improvement reached by the comparative methods in contrast to the proposed hybrid optimization approach, shows the total grid improvement about the power losses minimization.
Figure 4 shows the positive impact that has the inclusion of distributed generation in DC distribution networks, since all of the comparative methods and the proposed hybrid approach present improvements that are higher than 88 % with a total penetration of 60% in a distributed generation regarding total power consumption. In addition, these results confirm the global capabilities of the proposed hybrid GA-SOCP method for dealing with the problem of the optimal location and sizing DGs in DC feeders.

5.2. 69-Node Test Feeder

Table 4 presents the numerical results of the proposed and comparative methods. These results allow observing that:
🗸
All of the comparative methods are stuck in locally optimal solutions, including the MIQP approach and the MINLP model solved in GAMS through the SBB solver since these finds final power losses about 0.1660 pu and 0.1574 pu, respectively. However, the best optimal solution is reached by the GA-SOCP method with 0.1573 pu.
🗸
Some of the metaheuristic approaches (i.e., GA-CGA and GA-PSO) find the best positions for installing the GDs being these nodes 21, 61, and 64; nevertheless, due to the random procedures in the slave stages, these find local solutions.
🗸
Even if standard deviations of the GA-PSO and GA-CGA are lower than the proposed GA-SOCP, these values imply that these methods are always trapped in the local solutions, and they do not have the capacity of scape from these solutions to explore other promissory regions of the solution space.
Note that the best optimal solution reported by the GA-SOCP proposed approach connects DGs in nodes 21, 61, and 64 with generation capabilities of about 1.4140 pu, 10.2630 pu, and 3.8803 , respectively. Once again, we can confirm that this solution is the global optimum, since, after the exhaustive evaluation of the solution space (i.e., 50116 possible locations) via nested loops, the same set of nodes and power generations have been identified as the best global solution. It is important to highlight that this exhaustive evaluation takes about 7000 s to revise the whole solution space, which implies that 66.91 s reported by the proposed GA-SOCP approach demonstrates its efficiency and robustness global finding capabilities.
Figure 5 presents the histogram with the optimization behavior of the proposed GA-SOCP method for locating and sizing DGs in DC distribution networks.
From Figure 5, it can be noted that more than 70 % of the solutions that were found by the proposed GA-SOCP method find the global optimal solution, and 24 % of the additional solutions are better than the mean solutions reported in Table 4 (see μ column). These results confirm the proposed GA-SOCP optimization approach’s efficiency and robustness for solving the optimal placement and sizing of DGs in DC distribution networks with global optimum finding capabilities.
Figure 6 presents the percentage of power losses improvement reached by the comparative and the proposed optimization approaches.
Note that the percentage of reduction of power losses reported in plot Figure 6 confirms the global optimization capabilities of the proposed GA-SOCP approach, since it is the only one methodology that reports the highest percentage of reduction with 89.78 % , only followed by the SBB methodology with 89.77 % of reduction.

5.3. General Commentaries

After applying the comparative and the proposed optimization methodology for optimal locating and sizing DGs in DC distribution networks, we can affirm that:
🗸
The processing times in both test systems have been lower than 70 s, which implies that for planning purposes, our approach is perfectly applicable since, with a fraction of the time taken by the exhaustive search method, this can deal with the global optimal solution. In addition, the processing times of the comparative metaheuristic methods increase rapidly as a function of the test feeder size (see the last columns of Table 3 and Table 4) with the main disadvantage that due to the random behavior of the slave stages, i.e., BHO, CGA, and PSO, it is not possible to find the optimal solution of the studied optimization problem.
🗸
The methodologies implemented in GAMS, i.e., the MIQP and the MINLP model, both solved with the SBB tool available in GAMS, have speedy processing times. However, the linearization of the power flow equations in the MIQP approach and the nonlinear and non-convex structure of the exact MINLP model make that these get stuck in optimal locations as the solution space grows.
🗸
The main advantages of the proposed hybrid GA-SOCP master-slave optimization method are: (i) the possibility of using a compact codification with integer numbers (nodes) that allows for guaranteeing feasibility in all the stages of the genetic algorithm procedures, and (ii) the guarantee of the global optimal solution for each combination of distributed generators thanks to the second-order cone structure of the optimal power flow problem in the slave stage.

6. Conclusions and Future Works

A hybrid master-slave optimization approach for optimal placement and sizing distributed generators in DC distribution networks while using a combination of the classical Chu & Beasley genetic algorithm in the master stage and a second-order cone programming model in the slave stage, i.e., GA-SOCP. This algorithm has an important advantage regarding classical metaheuristics, since solving the optimal power flow problem for optimal dimensioning the distributed generators provided by the master stage with a SCOP model guarantees the existence and uniqueness of the optimal solution, which is not possible with combination methods, such as BHO, CGA, and PSO, where random procedures affect the optimization process. Numerical comparisons with classical metaheuristics and exact methodologies demonstrate that, for the studied test feeders, the proposed GA-SOCP approach deals with the global optimal solution. Simultaneously, it is not the case of the comparative methods in the 69-node test feeder scenario. The effectiveness of the proposed GA-SOCP method to find the optimal solution after 100 consecutive evaluations is about 73 % for the 21-node test feeder, and 93 % in the case of the 69-node test system, which implies the high efficiency of the studied hybrid optimization approach to deal with MINLP problems.
Regarding future works, it will be possible to make the following researches: (i) extend the proposed optimization methodology to the optimal location and sizing of distributed generators and capacitor banks in AC distribution networks, (ii) to apply the proposed GA-SOCP approach to the problem of the optimal placement and sizing battery energy storage systems in AC and DC networks, and (iii) to develop an economic dispatch approach to optimal locating renewable generation at isolated distribution networks when considering variations in the load curve and the primary energy resources.

Author Contributions

Conceptualization, O.D.M. and W.G.-G.; Methodology, O.D.M. and W.G.-G.; Investigation, O.D.M., W.G.-G., and L.F.G.-N.; Writing—review and editing, O.D.M., W.G.-G., and L.F.G.-N. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

This work was partially supported by the National Scholarship Program Doctorates of the Administrative Department of Science, Technology, and Innovation of Colombia (COLCIENCIAS), by calling contest 727-2015.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
Acronyms
ACAlternating Current
DGDistributed Generator
GAChu & Beasley genetic algorithm
HSMHyperplanes search method
NLPNonlinear Programming
CVXMatlab Software for Disciplined Convex Programming
GAMSGeneral Algebraic Modeling System
MIMixed-integer
MINLPMixed-integer nonlinear programming
SOCPSecond-order cone programming
SDPSemidefinite programming
Sets and subscripts
L Set of lines
N Set of nodes
dDemand
g d Distributed generator
sSlack node
i or jNode
Parameters
GConductance matrix
G i j Component of the conductance matrix that relates nodes i and j
N d g a v a Maximum number of distributed generators available
g k m Conductance between nodes k and m
p d Power demanded by loads
p s Power generated by generators
p g d Power generated by distributed source
round ( · ) Integer part of the argument
rand ( · ) Random number
Os i t Offsprings of population i in period t
Variables
v i Voltage profile at node i
z i j cross-product of voltage i with voltage j
w k m Population position k and m
y i binary variable associate with the location at node i
WPopulation matrix
Control variables
p s Power generated by conventional generator
p g d Power generated by renewable energies
p k m Power flow between nodes k and m
y i Binary variable for location of DG
P s Power generated vector by conventional generator
P g d Power generated vector by renewable energies
Limits
N g d max Maximum number of DG
v min , v max Minimum and maximum voltage profile
V min , V max Minimum and maximum voltage profile squared
p s , min , p s , max Minimum and maximum by conventional generator
p g d , min , p g d , max Minimum and maximum by generated source
p k m min , p k m max Minimum and maximum power flow between nodes k and m

References

  1. Kwon, M.; Choi, S. Control scheme for autonomous and smooth mode switching of bidirectional DC–DC converters in a DC microgrid. IEEE Trans. Power Electron. 2017, 33, 7094–7104. [Google Scholar] [CrossRef]
  2. Garcés, A. Stability Analysis of DC-Microgrids: A Gradient Formulation. J. Control Autom. Electr. Syst. 2019, 30, 985–993. [Google Scholar] [CrossRef]
  3. Wang, L.; Wang, K.H.; Lee, W.J.; Chen, Z. Power-flow control and stability enhancement of four parallel-operated offshore wind farms using a line-commutated HVDC link. In Wind Energy Conversion Systems; Springer: Berlin/Heidelberg, Germany, 2012; pp. 385–414. [Google Scholar]
  4. Gavriluta, C.; Candela, I.; Citro, C.; Luna, A.; Rodriguez, P. Design considerations for primary control in multi-terminal VSC-HVDC grids. Electr. Power Syst. Res. 2015, 122, 33–41. [Google Scholar] [CrossRef]
  5. Simiyu, P.; Xin, A.; Wang, K.; Adwek, G.; Salman, S. Multiterminal Medium Voltage DC Distribution Network Hierarchical Control. Electronics 2020, 9, 506. [Google Scholar] [CrossRef] [Green Version]
  6. Justo, J.J.; Mwasilu, F.; Lee, J.; Jung, J.W. AC-microgrids versus DC-microgrids with distributed energy resources: A review. Renew. Sustain. Energy Rev. 2013, 24, 387–405. [Google Scholar] [CrossRef]
  7. Parhizi, S.; Lotfi, H.; Khodaei, A.; Bahramirad, S. State of the Art in Research on Microgrids: A Review. IEEE Access 2015, 3, 890–925. [Google Scholar] [CrossRef]
  8. Kumar, D.; Zare, F.; Ghosh, A. DC microgrid technology: System architectures, AC grid interfaces, grounding schemes, power quality, communication networks, applications, and standardizations aspects. IEEE Access 2017, 5, 12230–12256. [Google Scholar] [CrossRef]
  9. Lotfi, H.; Khodaei, A. AC versus DC microgrid planning. IEEE Trans. Smart Grid 2015, 8, 296–304. [Google Scholar] [CrossRef]
  10. Gelani, H.; Dastgeer, F.; Siraj, K.; Nasir, M.; Niazi, K.; Yang, Y. Efficiency Comparison of AC and DC Distribution Networks for Modern Residential Localities. Appl. Sci. 2019, 9, 582. [Google Scholar] [CrossRef] [Green Version]
  11. Rouzbehi, K.; Miranian, A.; Candela, J.I.; Luna, A.; Rodriguez, P. A Generalized Voltage Droop Strategy for Control of Multiterminal DC Grids. IEEE Trans. Ind. Appl. 2015, 51, 607–618. [Google Scholar] [CrossRef]
  12. Rouzbehi, K.; Miranian, A.; Luna, A.; Rodriguez, P. DC Voltage Control and Power Sharing in Multiterminal DC Grids Based on Optimal DC Power Flow and Voltage-Droop Strategy. IEEE J. Emerg. Sel. Top. Power Electron. 2014, 2, 1171–1180. [Google Scholar] [CrossRef]
  13. Garcés, A. Convex Optimization for the Optimal Power Flow on DC Distribution Systems. In Handbook of Optimization in Electric Power Distribution Systems; Springer: Berlin/Heidelberg, Germany, 2020; pp. 121–137. [Google Scholar]
  14. Garcés, A. On the Convergence of Newton’s Method in Power Flow Studies for DC Microgrids. IEEE Trans. Power Syst. 2018, 33, 5770–5777. [Google Scholar] [CrossRef] [Green Version]
  15. Garces, A. Uniqueness of the power flow solutions in low voltage direct current grids. Electr. Power Syst. Res. 2017, 151, 149–153. [Google Scholar] [CrossRef]
  16. Hamad, A.A.; El-Saadany, E.F. Multi-agent supervisory control for optimal economic dispatch in DC microgrids. Sustain. Cities Soc. 2016, 27, 129–136. [Google Scholar] [CrossRef]
  17. Barabanov, N.; Ortega, R.; Griñó, R.; Polyak, B. On Existence and Stability of Equilibria of Linear Time-Invariant Systems With Constant Power Loads. IEEE Trans. Circuits Syst. I 2016, 63, 114–121. [Google Scholar] [CrossRef] [Green Version]
  18. Chauhan, R.K.; Chauhan, K.; Guerrero, J.M. Controller design and stability analysis of grid connected DC microgrid. J. Renew. Sustain. Energy 2018, 10, 035101. [Google Scholar] [CrossRef] [Green Version]
  19. Montoya, O.D.; Gil-González, W.; Grisales-Noreña, L. Relaxed convex model for optimal location and sizing of DGs in DC grids using sequential quadratic programming and random hyperplane approaches. Int. J. Electr. Power Energy Syst. 2020, 115, 105442. [Google Scholar] [CrossRef]
  20. Grisales-Noreña, L.; Montoya, O.D.; Ramos-Paja, C.A. An energy management system for optimal operation of BSS in DC distributed generation environments based on a parallel PSO algorithm. J. Energy Storage 2020, 29, 101488. [Google Scholar] [CrossRef]
  21. Fantauzzi, M.; Lauria, D.; Mottola, F.; Scalfati, A. Sizing energy storage systems in DC networks: A general methodology based upon power losses minimization. Appl. Energy 2017, 187, 862–872. [Google Scholar] [CrossRef]
  22. Castillo-Calzadilla, T.; Macarulla, A.M.; Kamara-Esteban, O.; Borges, C.E. A case study comparison between photovoltaic and fossil generation based on direct current hybrid microgrids to power a service building. J. Clean. Prod. 2020, 244, 118870. [Google Scholar] [CrossRef]
  23. Chiodo, E.; Fantauzzi, M.; Lauria, D.; Mottola, F. A Probabilistic Approach for the Optimal Sizing of Storage Devices to Increase the Penetration of Plug-in Electric Vehicles in Direct Current Networks. Energies 2018, 11, 1238. [Google Scholar] [CrossRef] [Green Version]
  24. Montoya, O.D. A convex OPF approximation for selecting the best candidate nodes for optimal location of power sources on DC resistive networks. Eng. Sci. Technol. Int. J. 2020, 23, 527–533. [Google Scholar] [CrossRef]
  25. Montoya, O.D.; Garrido, V.M.; Grisales-Noreña, L.F.; Gil-González, W.; Garces, A.; Ramos-Paja, C.A. Optimal Location of DGs in DC Power Grids Using a MINLP Model Implemented in GAMS. In Proceedings of the 2018 IEEE 9th Power, Instrumentation and Measurement Meeting (EPIM), Salto, Uruguay, 14–16 November 2018; pp. 1–5. [Google Scholar]
  26. Montoya, O.D.; Grisales-Noreña, L.F.; Gil-González, W.; Alcalá, G.; Hernandez-Escobedo, Q. Optimal Location and Sizing of PV Sources in DC Networks for Minimizing Greenhouse Emissions in Diesel Generators. Symmetry 2020, 12, 322. [Google Scholar] [CrossRef] [Green Version]
  27. Grisales-Noreña, L.F.; Garzon-Rivera, O.D.; Montoya, O.D.; Ramos-Paja, C.A. Hybrid Metaheuristic Optimization Methods for Optimal Location and Sizing DGs in DC Networks. In Proceedings of the Workshop on Engineering Applications, Santa Marta, Colombia, 16–18 October 2019; Springer: Berlin/Heidelberg, Germany, 2019; pp. 214–225. [Google Scholar]
  28. Huang, L.; Chen, Z.; Cui, Q.; Zhang, J.; Wang, H.; Shu, J. Optimal planning of renewable energy source and energy storage in a medium- and low-voltage distributed AC/DC system in China. J. Eng. 2019, 2019, 2354–2361. [Google Scholar] [CrossRef]
  29. Li, J.; Liu, F.; Wang, Z.; Low, S.H.; Mei, S. Optimal Power Flow in Stand-Alone DC Microgrids. IEEE Trans. Power Syst. 2018, 33, 5496–5506. [Google Scholar] [CrossRef] [Green Version]
  30. Benson, H.Y.; Sağlam, Ü. Mixed-Integer Second-Order Cone Programming: A Survey. In Theory Driven by Influential Applications; INFORMS: Catonsville, MD, USA, 2013; pp. 13–36. [Google Scholar] [CrossRef] [Green Version]
  31. Li, Q.; Vittal, V. Convex Hull of the Quadratic Branch AC Power Flow Equations and Its Application in Radial Distribution Networks. IEEE Trans. Power Syst. 2018, 33, 839–850. [Google Scholar] [CrossRef] [Green Version]
  32. Farivar, M.; Low, S.H. Branch Flow Model: Relaxations and Convexification-Part I. IEEE Trans. Power Syst. 2013, 28, 2554–2564. [Google Scholar] [CrossRef]
  33. Lambora, A.; Gupta, K.; Chopra, K. Genetic Algorithm- A Literature Review. In Proceedings of the 2019 International Conference on Machine Learning, Big Data, Cloud and Parallel Computing (COMITCon), Faridabad, India, 14–16 February 2019; IEEE: Piscataway, NJ, USA, 2019. [Google Scholar] [CrossRef]
  34. Montoya, O.D.; Gil-González, W.; Orozco-Henao, C. Vortex search and Chu-Beasley genetic algorithms for optimal location and sizing of distributed generators in distribution networks: A novel hybrid approach. Eng. Sci. Technol. Int. J. 2020. [Google Scholar] [CrossRef]
  35. Sahin, O.; Akay, B. Comparisons of metaheuristic algorithms and fitness functions on software test data generation. Appl. Soft Comput. 2016, 49, 1202–1214. [Google Scholar] [CrossRef]
  36. Badem, H.; Basturk, A.; Caliskan, A.; Yuksel, M.E. A new hybrid optimization method combining artificial bee colony and limited-memory BFGS algorithms for efficient numerical optimization. Appl. Soft Comput. 2018, 70, 826–844. [Google Scholar] [CrossRef]
  37. Montoya, O.D.; Gil-González, W. A MIQP model for optimal location and sizing of dispatchable DGs in DC networks. Energy Syst. 2020, 1–22. [Google Scholar] [CrossRef]
Figure 1. Electrical connection between nodes in the 21-node test feeder.
Figure 1. Electrical connection between nodes in the 21-node test feeder.
Applsci 10 08616 g001
Figure 2. Schematic connection among nodes for the 69-node test feeder.
Figure 2. Schematic connection among nodes for the 69-node test feeder.
Applsci 10 08616 g002
Figure 3. Behavior of the proposed hybrid GA-PSO after 100 consecutive evaluations in the 21-node test feeder.
Figure 3. Behavior of the proposed hybrid GA-PSO after 100 consecutive evaluations in the 21-node test feeder.
Applsci 10 08616 g003
Figure 4. Percentage of improvement regarding power losses minimization in the 21-node test feeder.
Figure 4. Percentage of improvement regarding power losses minimization in the 21-node test feeder.
Applsci 10 08616 g004
Figure 5. Behavior of the proposed hybrid GA-PSO after 100 consecutive evaluations in the 69-node test feeder.
Figure 5. Behavior of the proposed hybrid GA-PSO after 100 consecutive evaluations in the 69-node test feeder.
Applsci 10 08616 g005
Figure 6. Percentage of improvement regarding power losses minimization in the 69-node test feeder.
Figure 6. Percentage of improvement regarding power losses minimization in the 69-node test feeder.
Applsci 10 08616 g006
Table 1. Parametric information of the 21-node test system.
Table 1. Parametric information of the 21-node test system.
Node iNode j R ij [p.u] P j [p.u]Node iNode j R ij [p.u] P j [p.u]
120.00530.7011120.00790.68
130.00540.0011130.00780.10
340.00540.3610140.00830.00
450.00630.0414150.00650.22
460.00510.3615160.00640.23
370.00370.0016170.00740.43
780.00790.3216180.00810.34
790.00720.8014190.00780.09
3100.00530.0019200.00840.21
10110.00380.4519210.00820.21
Table 2. Parametric information of the 69-node test system.
Table 2. Parametric information of the 69-node test system.
Node iNode j R ij [ Ω ] P j [kW]Node iNode j R ij [ Ω ] P j [kW]
120.000503360.004426
230.0005036370.064026
340.0015037380.10530
450.0215038390.030424
560.36602.639400.001824
670.381040.440410.7283102
780.09227541420.31000
890.04933042430.04106
9100.81902843440.00920
10110.187214544450.108939.22
11120.711414545460.000939.22
12131.030084470.00340
13141.0440847480.085179
14151.0580048490.2898384.7
15160.19664549500.0822384.7
16170.3744608510.092840.5
17180.00476051520.33193.6
18190.327609530.17404.35
19200.2106153540.203026.4
20210.341611454550.284224
21220.0140555560.28130
22230.1591056571.59000
23240.34632857580.78370
24250.7488058590.3042100
25260.30891459600.38610
26270.17321460610.50751244
3280.00442661620.097432
28290.06402662630.14500
29300.3978063640.7105227
30310.0702064651.041059
31320.3510065660.201218
32330.83901066670.004718
33341.70801467680.739428
34351.4740468690.004728
Table 3. Comparison of different approaches in the 21-node test feeder.
Table 3. Comparison of different approaches in the 21-node test feeder.
MethodNodes σ [p.u.] μ [p.u.]min [p.u.] t ave [s]
GA-BHO{9,12,16}2.2761 × 10−030.03680.0318111.9440
GA-CGA{9,12,16}1.3537 × 10−030.03290.031133.1974
GA-PSO{9,12,16}1.8437 × 10−030.03190.030658.0934
MIQP{9,12,16} 0.03060.03068.9688
SSB{9,12,16} 0.03060.03065.3653
GA-SOCP{9,12,16}1.1050 × 10−030.03090.030652.146
Table 4. Comparison of different approaches in the 69-node test feeder.
Table 4. Comparison of different approaches in the 69-node test feeder.
MethodNodes σ [p.u.] μ [p.u.]min [p.u.] t ave [s]
GA-BHO{23,61,67}2.5207 × 10−030.16330.1593713.7193
GA-CGA{21,61,64}3.4801 × 10−040.16480.1603218.0169
GA-PSO{21,61,64}5.4023 × 10−040.16890.1588378.4731
MIQP{22,61,64} 0.16600.166024.1868
SBB{22,61,65} 0.15740.157412.4270
GA-SOCP{21,61,64}3.1295 × 10−030.15910.157366.9100
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Montoya, O.D.; Gil-González, W.; Grisales-Noreña, L.F. Hybrid GA-SOCP Approach for Placement and Sizing of Distributed Generators in DC Networks. Appl. Sci. 2020, 10, 8616. https://doi.org/10.3390/app10238616

AMA Style

Montoya OD, Gil-González W, Grisales-Noreña LF. Hybrid GA-SOCP Approach for Placement and Sizing of Distributed Generators in DC Networks. Applied Sciences. 2020; 10(23):8616. https://doi.org/10.3390/app10238616

Chicago/Turabian Style

Montoya, Oscar Danilo, Walter Gil-González, and Luis Fernando Grisales-Noreña. 2020. "Hybrid GA-SOCP Approach for Placement and Sizing of Distributed Generators in DC Networks" Applied Sciences 10, no. 23: 8616. https://doi.org/10.3390/app10238616

APA Style

Montoya, O. D., Gil-González, W., & Grisales-Noreña, L. F. (2020). Hybrid GA-SOCP Approach for Placement and Sizing of Distributed Generators in DC Networks. Applied Sciences, 10(23), 8616. https://doi.org/10.3390/app10238616

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