Next Article in Journal
A Fast Lock All-Digital MDLL Using a Cyclic Vernier TDC for Burst-Mode Links
Previous Article in Journal
Multi-Objective Sizing Optimization of Hybrid Renewable Energy Microgrid in a Stand-Alone Marine Context
Previous Article in Special Issue
Investigation on Hybrid Energy Storage Systems and Their Application in Green Energy Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Mixed-Integer Conic Formulation for Optimal Placement and Dimensioning of DGs in DC Distribution Networks

by
Federico Molina-Martin
1,
Oscar Danilo Montoya
2,3,
Luis Fernando Grisales-Noreña
4 and
Jesus C. Hernández
5,*
1
Estudiante de Doctorado, Departamento de Ingeniería Eléctrica, Campus Lagunillas s/n, University of Jaén, Edificio A3, 23071 Jaén, Spain
2
Facultad de Ingeniería, Universidad Distrital Francisco José de Caldas, Bogotá 11021, Colombia
3
Laboratorio Inteligente de Energía, Universidad Tecnológica de Bolívar, Cartagena 131001, Colombia
4
Departamento de Electromecánica y Mecatrónica, Instituto Tecnológico Metropolitano, Medellín 050013, Colombia
5
Department of Electrical Engineering, Campus Lagunillas s/n, University of Jaén, Edificio A3, 23071 Jaén, Spain
*
Author to whom correspondence should be addressed.
Electronics 2021, 10(2), 176; https://doi.org/10.3390/electronics10020176
Submission received: 18 December 2020 / Revised: 12 January 2021 / Accepted: 12 January 2021 / Published: 14 January 2021
(This article belongs to the Special Issue Microgrids and Active Distribution Networks)

Abstract

:
The problem of the optimal placement and dimensioning of constant power sources (i.e., distributed generators) in electrical direct current (DC) distribution networks has been addressed in this research from the point of view of convex optimization. The original mixed-integer nonlinear programming (MINLP) model has been transformed into a mixed-integer conic equivalent via second-order cone programming, which produces a MI-SOCP approximation. The main advantage of the proposed MI-SOCP model is the possibility of ensuring global optimum finding using a combination of the branch and bound method to address the integer part of the problem (i.e., the location of the power sources) and the interior-point method to solve the dimensioning problem. Numerical results in the 21- and 69-node test feeders demonstrated its efficiency and robustness compared to an exact MINLP method available in GAMS: in the case of the 69-node test feeders, the exact MINLP solvers are stuck in local optimal solutions, while the proposed MI-SOCP model enables the finding of the global optimal solution. Additional simulations with daily load curves and photovoltaic sources confirmed the effectiveness of the proposed MI-SOCP methodology in locating and sizing distributed generators in DC grids; it also had low processing times since the location of three photovoltaic sources only requires 233.16 s, which is 3.7 times faster than the time required by the SOCP model in the absence of power sources.

1. Introduction

Direct current (DC) distribution networks have attracted much attention in recent years in specialized literature [1], since these networks have better voltage profiles [2] and low energy losses [3]. They are easily controllable since frequency or reactive power are nonexistent concepts in these electrical networks [4]. In the literature, multiple approaches regarding DC networks have been described, such as control methodologies for power electronic converters that interface renewable energies [5] and batteries [6], optimization approaches associated with power losses minimization [2] and analysis of convergence algorithms for power flow analysis [7].
Here, we explore the optimization approaches recently presented in scientific literature to analyze DC distribution networks under steady-state conditions regarding energy losses minimization. The authors in [3] performed a semidefinite programming relaxation for the optimal locating and dimensioning of distributed generators (DGs) in DC grids, combined with a heuristic method based on random hyperplanes. This approach demonstrated the possibility of finding better solutions than the different metaheuristic approaches. A similar approach was presented in [8] with a quadratic approximation of the optimal power flow problem; this was combined with a heuristic search based on hyperplanes to determine the best set of nodes for locating constant power sources. The authors of [9] had proposed the complete MINLP model for optimal placement and sizing of DGs in DC grids, taking into consideration the different levels of power penetration. This model was solved by using the GAMS software and comparing the different solvers available on it for addressing MINLP models. The authors of [10] proposed a mixed-integer quadratic transformation of the original MINLP model by solving both models in the GAMS software; they observed that the quadratic approximation was easier to solve and took tens of seconds, while the exact MINLP model took hundreds of seconds. The former also had the best numerical results for the minimization of the total grid losses. In [11] was proposed a multi-period MINLP model for the location and sizing of photovoltaic sources in DC distribution networks located in rural areas with the objective of minimizing the total greenhouse emissions (i.e., carbon dioxide CO 2 ) caused by diesel generators. The solution of the MINLP model was reached with the GAMS optimization package. The authors of [12] proposed multiple combinations of combinatorial optimization algorithms for optimal the placement and dimensioning of DGs in DC networks via sequential quadratic programming [13]. The optimal placement problem of generators was addressed using the population-based incremental learning (PBIL) optimization algorithm. Simultaneously, the sizing problem was solved using continuous metaheuristic algorithms such as particle swarm optimizer (PSO), constant genetic algorithm, and black hole optimizer. Numerical results showed that the best trade-off is found using the PBIL-PSO method. However, this approach has not yet been compared with exact techniques.
The DC distribution optimization problems were presented in [14,15], in which the issue of the optimal reconfiguration of DC feeders was addressed using convex approximations and exact MINLP models. The authors of [16] proposed a reformulation of the exact nonlinear programming (NLP) model for optimal power flow studies in DC grids using a conic representation of the power balance equations. This reformulation was compared with the exact one in different test feeders, and it was confirmed that convex approximation can find the global optimum, which coincides with the solution of the NLP model. The author of [17] proposed a convexification of the exact NLP model for optimal operation of battery energy storage systems using a semidefinite programming model, which was solved using the CVX package on the MATLAB environment. The semidefinite programming model ensured the global optimum finding due to the convex properties of the semidefinite matrices; in addition, this work has been further improved through the reduction of solution space by proposal of a second-order cone programming equivalent, as was presented in [18], that also ensures global optimum finding. The main advantage of this was the reduction of the processing times by more than 50% compared to the semidefinite programming approach.
Additional developments in the current literature regarding the inclusion of stochastic behaviors in expansion planning problems for distribution networks, transmission grids and generation sources have been proposed in [19,20,21], respectively. These works demonstrated the importance of having adequate representations of the electrical networks, including the related uncertainties in renewable generation, to make their expansion effective using efficient techniques that involve exact and metaheuristic optimization methods depending on the complexity of the problem under study.
Based on the importance of DC distribution networks in the literature, in this paper, we propose an exact mathematical optimization based on second-order cone programming (SOCP) optimization for optimal power flow analysis [16]. At the same time, integer optimization based on the branch and bound (B&B) method is used to solve the problem of the optimal placement and sizing of constant power sources in DC grids. This combination generates a mixed-integer second-order cone programming (MI-SOCP) reformulation of the exact MINLP model presented in [22].
The most important contributions of this paper can be condensed as follows:
The guarantee of reaching the global optimum for the problem of the optimal siting and dimensioning of constant power sources in DC grids via an MI-SOCP model solved through B&B and interior-point methods.
The verification of the complexity that the original MINLP model presents in finding the optimal solution using powerful optimization tools available in GAMS. Numerical results in the 21-bus test feeder are optimal based on the GAMS solutions. However, in the 69-bus test feeder case, most of the solvers available in GAMS are stuck in local optimum solutions. This is not the case with the proposed MI-SOCP model, in which global solutions are reached for both test systems.
The remainder of this document takes the following structure: Section 2 describes the exact MINLP formulation for the problem of the optimal siting and dimensioning of constant power sources in DC grids by remarking the main non-convexities of this model. Section 3 offers the reformulation of optimal power flow problems via SOCP. Additionally, the SOCP is combined with the model’s binary nature to produce an MI-SOCP reformulation for the exact MINLP optimization model. Section 4 describes the main properties of the MI-SOCP model and its solution via B&B and interior-point methods. Section 5 presents the main characteristics of the test feeders, which are composed of 21 and 69 nodes and have radial topology and the possibility of locating three distributed power sources. Section 6 shows the numerical performance of the proposed MI-SCOP and its comparison with multiple nonlinear solvers available in the GAMS optimization package. Finally, Section 7 lists the leading conclusions obtained from the proposed study as well as some possible future developments.

2. Exact MINLP Model

The optimal siting and dimensioning of power sources in DC distribution networks is a complex optimization task in electrical engineering since it combines integer and continuous variables [10]. The continuous variables are related to power generations and voltage profiles, whereas the integer ones are associated with the possibility of placing (or not) in a particular node a distributed generator [9]. The exact MINLP formulation that describes the problem under study is shown below:
(1) Objective function : min p loss = { k m } L g k m v k v m 2 , (2) Subject to : p k s s + p k d g p k d = v k m N G k m v m , k N (3) p k m = g k m v k 2 v k v m , k m L (4) p k m min p k m p k m max , k m L (5) x k p k g d , min p k d g x k p k d g , max , k N (6) v min v k v max , k N (7) k N p k d g α k N p k d , (8) k N x k N g d a v a , (9) x k 0 , 1 , k N ,
where p loss is the objective function value associated with the total grid power losses; v k and v m are the values of the voltage variables associated with nodes k and m, respectively; g k m is the conductance parameter of the branch that connects nodes k and m; G k m corresponds to the element of the matrix of conductances that is associated with nodes k and m; p k s s , p k d g , and p k d are the power generation in slack nodes (voltage-controlled sources) and distributed generators, and power consumptions in load nodes, respectively; p k m is the flow of power in the line that connects nodes k and m; p k m min and p k m max are the lower and upper power flow bounds permitted at the branch k m ; p k g d , min and p k g d , max are the lower and upper power generation limits associated with a distributed generator connected at node k; x k is the binary variable associated with the placement ( x k = 1 ) or non-placement ( x k = 0 ) of a constant power source at node k; v min and v max are the minimum and maximum voltage regulation bounds allowed at each node; α is a positive constant between 0 and 1 that is related to the amount of distributed generation that is permitted in the DC grid; and N g d a v a refers to the number of power sources available for installation.
The complete interpretation of the general MINLP formulation for the problem of the optimal placement and sizing of DGs in DC grids defined from (1) to (9) is as follows: Equation (1) corresponds to the objective function of the optimization problem, which is associated with the total grid losses in all the branches of the network. This objective function is formulated using the voltage drop at each line and its conductance. The Equation (2) corresponds to the power balance constraint at each bus of the systems, which defines the hyperbolic relation between the voltages and powers generated or consumed. Please note that this set of equations defines the nonlinear quadratic constraints that make the continuous part of the problem non-convex [16]. Expression (3) determines the amount of power flow sent from node k to node m as a function of the voltage drop in the line, its conductance and the sending voltage. The inequality constraint (4) limits the amount of power flow that can be transported through the distribution line that connects nodes k and m. Please note that this restriction corresponds to the maximum thermal bound supported by the conduction written as a function of its power flow. The inequality constraint (5) determines the maximum amount of power that can be injected by a distributed generator connected at node k if the binary variable x k is activated, i.e., x k = 1 , if the DG is not connected at node k, then x k = 0 , which makes the amount of power injected be zero. Inequality expression (6) is known in the current literature as the voltage regulation constraint, since this limits the admissible voltage profiles in all the nodes of the distribution network. The upper and lower voltage bounds are defined by regulatory policies, and the operative consigns are defined by the grid operator. The inequality constraint (7) determines the maximum amount of power that can be injected via the distributed generators into the distribution network, which is also defined by regulatory policies in the electric distribution sector [10]. The inequality constraint (8) limits the maximum number of DGs that can be installed in the distribution network. This limitation is typically assigned by the grid operator as a function of the investment budget available. Finally, the binary nature of the decision variable, i.e., x k , is defined through (9), which defines the placement (or non-placement) of a constant power source in the DC grid.
Remark 1.
The expression (1), related to the objective function, is a quadratic function associated with the sums of the products between voltage variables; this is convex due to the properties associated with the matrix of conductance, G which is positive semidefinite [23,24].
Remark 2.
The set of nonlinear constraints associated with the balance of active power at each node and the power flow at each line, i.e., (2) and (3), are not convex as these generate quadratic equality constraints that are non-affine and non-convex [25]. However, as the unique non-linearity is the product between voltage variables, i.e., v k v m , these equations can be convexified using an SOCP approximation [18]. This has been described in the next section.

3. MI-SOCP Reformulation

The conic programming approach corresponds to a part of the convex optimization [18], which is a field of exact mathematical optimization that has attracted great attention in electrical engineering. The main advantage of it is that it can solve convex problems, guaranteeing a unique solution (global optimum) in a reliable and efficient manner [26]. The SOCP model works with a linear function over a convex region with the goal of minimizing it. In summary, the SOCP is an intersection of an affine linear space and conic constraints [27].
Here, we employ the hyperbolic equivalent of the product with two variables for rewriting Equations (1)–(3) using a conic representation [16]. Let us define the product between voltage variables as follows:
y k m = v k v m ,
in addition, if we pre-multiply (10) by v k v m on both sides, the following result is reached:
y k m 2 = v k 2 v m 2 .
Now, let us define the two auxiliary variables u k = v k 2 and u m = v m 2 ; they imply that (11) can be rewritten as follows:
y k m 2 = u k u m .
Please note that the product between the auxiliary variables u k and u m can be redefined by using its hyperbolic representation, i.e.,
y k m 2 = u k u m , = 1 4 u k + u m 2 1 4 u k u m 2 , 4 y k m 2 + u k u m 2 = u k + u m 2 , 2 y k m u k u m = u k + u m .
Remark 3.
Observe that Expression (13) continues being non-convex due to the equality restriction; however, to obtain a conic convex relaxation, the equality symbol can be replaced with the lower-equal symbol, as suggested in [28], which transforms this constraint into a second-order one.
Relaxing the equality constraint (13) allowed us to recover the original variables as y k k = u k and y m m = u m , which implies that (13) takes the following form:
2 y k m y k k y m m y k k + y m m .
According to the SOCP relaxation, as mentioned above, the exact MINLP formulation (1)–(9) can be rewritten as (15), which corresponds to an MI-SOCP reformulation for the problem of siting and dimensioning constant power sources in DC grids.
Objective function : min p loss = { k m } L g k m y k k 2 y k m + y m m , Subject to : p k s s + p k d g p k d = m N G k m y k m , k N p k m = g k m y k k y k m , k m L p k m min p k m p k m max , k m L x k p k g d , min p k d g x k p k d g , max , k N v min 2 y k k v max 2 , k N k N p k d g α k N p k d , k N x k N g d a v a , x k 0 , 1 , k N
Remark 4.
The main characteristic of the proposed MI-SOCP formulation defined in (15) is that it guarantees the reaching of the global optimum if it is solved through the B&B method since each explored node is convex, i.e., it has a unique solution. This has been briefly explained in the next section.

4. Strategy of Solution

As mentioned in the previous section, the SOCP is a sub-field of the convex optimization that works with the optimization models composed of linear affine and conic constraints [27,29]. Multiple algorithms are available for efficiently solving these types of large SOCP models in only microseconds. To exemplify the structure of a conic constraint, let us define a general second-order cone as follows:
x z ,
where x R n and z R , and x are the Euclidean norm, i.e., norm-2, of the vector x. Figure 1 shows a second-order cone in R 3 , clearly a convex set [30]. For further details about optimization convex, see references [31,32].
The problem of the optimal siting and dimensioning of power sources in DC grids can be addressed with a combination of the SCOP formulation and the B&B algorithm. In general terms, an MI-SOCP takes the following structure:
A i x + b i α i x + β i z i + θ i ,
where decision variables are composed by a combination of continuous and binary ones (i.e., x and z); A i are real matrices; b i , α i and β i are real vectors; and θ i refers to the constants for each constraint i.
The solution of the MI-SOCP with the general structure (18) can be reached as follows: (i) each iteration of the B&B method defines the combination of the binary variables associated with the power sources’ placement (exploration node); (ii) for each one of the explored nodes is solved the SOCP relaxation, which deals with the problem of power sources’ sizing. This method benefits from the properties of the SOCP problems related to convexity and the fast convergence of the interior-point methods [33], and guarantees the finding of the optimal global solution at each node. Figure 2 illustrates the schematic solution of an MI-SOCP problem with two binary variables.
Please note that SOCPs have been widely studied in the literature, for example, the authors of [34] reported a general introduction of SOCP and an extensive list of applications. The authors of [27] presented an overview of the general properties, the duality theory and the interior-point methods applied to SOCP problems. Interior-point techniques have been the most popular methods for solving this class of problems, offering good theoretical convergence properties [35] and efficient computational performance in numerical validations, such as SeDuMi [36], and SDPT3 [37].
Remark 5.
To solve the MI-SOCP formulation reported in (15) to locate and size constant power sources in DC grids, the CVX optimization package available for MATLAB and the MOSEK solver has been selected here [38].
It is worth mentioning that in the case of the optimal power flow problems, we can ensure that the optimal solution of the original NLP model is the same one reached by conic approximations, as demonstrated in [39]; however, there exist multiple simple NLP programming with few variables, as in the case of the parametric estimation in photovoltaic modules [40], single-phase transformers [41] or induction motors [42], in which it is not possible to obtain conic equivalents. In this sense, the main limitation of the conic programming is its application to NLP models that only contain products between continuous variables.

5. Test Feeders

The computational testing of the proposed MI-SOCP reformulation was done in two conventional DC test feeders. These had 21 and 69 buses. These test feeders have been reported in the literature for addressing the studied problem with MINLP models [9] and heuristic methods [3]. The information regarding these distribution grids has been presented below.

5.1. 21-Bus Test Feeder

This is an electrical DC distribution grid conformed by 20 lines and 21 nodes in a radial connection (see Figure 3). The slack node (i.e., voltage-controlled source) is located at node 1 and supports a voltage magnitude of 1.00 pu [3]. The electrical parameters of this test system, i.e., line and load parameters, have been listed in Table 1 [43].
It is worth mentioning that the information reported in Table 1 was calculated with a voltage base of 1.0 kV, and a power base of 100 kW. Please note that with these values we calculated the base of the resistance to be 10 Ω .

5.2. 69-Bus Test Feeder

The 69-bus test system is an electrical distribution grid conformed by 68 lines and 69 nodes connected in a radial form and operated with 12.66 kV in the root node (see Figure 4) [44]. Here, we consider the DC adaptation of this test feeder as informed in [3,8,43]. The data about nodal and branch parameters has been listed in Table 2. Mark that the information reported in Table 1 was calculated with a voltage base of 12.66 kV and a power base of 100 kVA.

5.3. Implementation Characteristics of the Test Feeders

For the numerical implementation of both test feeders, the following characteristics were taken into account: (i) the possibility of installing three DGs in both test feeders, with single capacities of 1.5 pu for the 21-bus test system and 12 pu for the 69-bus test feeder; (ii) the admissible active power injection via distributed generation was fixed as 60% in both test feeders, where the percentage of penetration was related to the total power consumption of each test feeder.

6. Computational Implementation

The developed optimization model based on MI-SOCP formulation was programmed in the MATLAB 2020a with the help of the CVX package and the MOSEK solver in a personal 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 Single language. To validate our mixed-integer convex reformulation’s robustness and effectiveness for the optimal siting and dimensioning of constant power sources in DC grids, we compared our model with multiple MINLP solvers available in the GAMS optimization tool [10].

6.1. Solution under the Peak Load Condition

Table 3 reports the optimal placement and sizing of the constant power sources in both DC distribution grids.
The results reported in Table 3 yield the following observations:
All the comparative solvers in the 21-bus test system reach the same numerical solution reported by the proposed MI-SOCP model, which corresponds to the final power losses of about 3.06 kW. This confirms that nodes 9, 12 and 16 are the ones nodes associated with the global solution of the problem of siting and dimensioning of DGs in DC grids.
For the 69-bus test feeder, we can highlight that the proposed MI-SOCP reformulation finds the global solution with 4.14 kW of final power losses. The BONMIN solver is the only solver that reaches this solution since the other optimization tools in GAMS are trapped in local solutions. The best places for locating constant power sources in this test feeder correspond to the nodes 17, 61 and 64, with a total power injection of about 2294.8 kW.
As for the power losses improvement in the 21-bus test system regarding the base case ( 0.2760 pu), this was about 88.91 % , while in the 69-bus test feeder (base case with 1.5385 pu) this reduction was about 97.31 % .
Even though the ANTIGONE solver reaches the global solution for the 21-bus test system, it fails in the 69-bus test system because it cannot find a combination of three nodes to minimize power losses; it only identifies one, leaving the other two options free. Therefore, this solver was not reported in the second test feeder.
Figure 5 shows the power losses enhancement found by each one of the solvers tested in GAMS and our proposed MI-SOCP for the 69-bus test system.
Please note that all the solution techniques presented in Figure 5 allow for improvements of higher than 90 % in power losses. Nevertheless, the best solution was found by the BONMIN and the MI-SOCP proposal, and the COUENNE solver had the worse result.
Remark 6.
The solutions reported by the GAMS solvers, as shown in Table 3, demonstrate that the MINLP model presented from (1) to (9) is complex to solve, and that it is impossible to ensure the global optimum achievement by using the GAMS optimization tools since these can be trapped in local solutions. However, the proposed MI-SOCP allows for the reaching of the global optimum by combining the B&B method and SOCP optimization, as presented in Table 3, for the 21- and 69-bus test feeders.
It is worth mentioning that in this research, we decided not to compare the proposed MI-SOCP model with metaheuristic optimization techniques as these are highly dependent on the programmer. These also require the tuning of multiple parameters, while this is not the case with the proposed model, which has exact mathematical methods. In addition, validation of metaheuristic methods requires statistical tests since their random procedures result in different solutions at each running. Thus, these are not recommended in problems where convex optimization is applied, as the one studied in this article.
Remark 7.
As for the comparison of the proposed MI-SOCP with the results reported in [3] and [8], where convex optimization methods based on semidefinite and sequential quadratic programming models were combined with a heuristic algorithm based on hyperplanes, we can affirm that our MI-SOCP approach allows for the reaching of the global optimum. It has better performance as regards the metaheuristics presented in those papers.
To demonstrate that the solutions found by the MI-SOCP formulation corresponded to the global optima, we evaluated all the possible combinations of three constant power sources in the 21- and 69-bus test feeder, which yielded about 1140 and 50,116 options, respectively. These exhaustive searches demonstrated that the nodes 9, 12 and 16 are the best places for the power sources in the 21-bus test feeder, and nodes 17, 61 and 64 are the best possible solution in 69-bus test feeder.

6.2. Solution Considering the Installation of Photovoltaic Generators

In this section, we discuss the ability of the proposed MI-SOCP formulation to site and size renewable distributed generators in DC grids, taking into consideration the daily load variations in the 69-bus test system (see mathematical model (16)). The curves analyzed have been presented in Figure 6. In this simulation, we considered one to three distributed generators available for installation without the limitations associated with the renewable energy penetration.
Objective function : min E loss = t T { k m } L g k m y k k , t 2 y k m , t + y m m , t , Subject to : p k , t s s + p k , t d g p k , t d = m N G k m y k m , t , k N , t T p k m , t = g k m y k k , t y k m , t , k m L , t T p k m min p k m , t p k m max , k m L , t T x k p k , t g d , min p k , t d g x k p k , t d g , max , k N , t T v min 2 y k k , t v max 2 , k N , t T k N max t T p k , t d g α k N max t T p k , t d k N x k N g d a v a , x k 0 , 1 , k N
Table 4 presents the numerical results associated with the installation of PV sources in DC grids using the proposed MI-SOCP reformulation.
The numerical results presented in Table 4 show the following: (i) in all the combinations of PV sources, node 61 is identified in the set of optimal solutions with power injections higher than 2000 kW; (ii) an injection of 2119.4980 kW at node 61 allows for the reduction of the daily energy losses by about 31.73 % with only one PV source; however, three PV sources in nodes 18, 49 and 61 inject about 3425.0456 , leading to a daily energy losses reduction of about 34.25 % . This difference corresponds to 46.2832 kWh/Day which no justifies an increment about 1000 kW in the PV capacity; (iii) the required processing time in the solution of the multi-period MI-SOCP increases as a function of the number of PV sources that will be located, since the dimension of the discrete space increases rapidly. In the case of one PV source, the dimension of the solution space is 68 nodal options, while in the case of three PVs the size of the solution space is 50,116 nodal combinations.
Please note that the numerical results reported in Table 4 confirm the robustness and effectiveness of the proposed MI-SOCP formulation in locating and sizing power sources in DC grids in one or over multiple periods of analysis; on the other hand, the software GAMS is trapped in local optimal solutions.

7. Conclusions and Future Works

This research paper proposed an MI-SOCP reformulation for the optimal siting and dimensioning of constant power sources in DC grids that allows for the reaching of global optimal solutions, which is not possible with the exact MINLP model. This solution of the proposed MI-SOCP model was reached by combining the classical branch and bound method with interior-point algorithms for SOCP problems. This combination guarantees the global optimum since each node explored by the branch and bound algorithm is, in fact, convex, i.e., it has a unique global solution.
Computational validations in the 21- and 69-bus test feeders showed that our proposed (i.e., the MI-SOCP approach) model has the best numerical performance when compared to multiple MINLP solvers in GAMS optimization since these reached the global solution only in the case of the 21-bus test feeder; they were stuck in locally optimal solutions in the 69-bus test feeder case due to the non-convex properties of the exact MINLP model.
The main result obtained after applying the proposed MI-SOCP model for optimal siting and dimensioning of constant power sources in the DC networks was the global optimization properties of this method. This demonstrated that this kind of MINLP problem can be addressed using exact mathematical optimization without going back to metaheuristic optimization techniques that cannot guarantee global optimal solution finding at every running. In addition, simulations considering daily load curves and renewable generation demonstrated this model’s applicability to distribution grids with high penetration of PV sources; another benefit of this model is that it requires low computational effort.
Future research can explore the following points: (i) the application of the proposed MI-SOCP model to solve the problem of the optimal location of capacitor banks and distributed generators in AC grids; (ii) the formulation of the economic dispatch problem, taking into consideration batteries and renewable generation in AC distribution networks, via SOCP representations; (iii) the study of the optimal location of batteries and renewables in DC grids with an MI-SOCP representation.

Author Contributions

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

Funding

This research was funded by the Agencia Estatal de Investigación, Spain (AEI) and the Fondo Europeo de Desarrollo Regional (FEDER) aimed at the Challenges of Society (grant No. ENE 2017-83860-R “Nuevos servicios de red para microredes renovables inteligentes. Contribución a la generación distribuida residencial”).

Acknowledgments

This work was supported in part by the Centro de Investigación y Desarrollo Científico de la Universidad Distrital Francisco José de Caldas under grant 1643-12-2020 associated with the project: “Desarrollo de una metodología de optimización para la gestión óptima de recursos energéticos distribuidos en redes de distribución de energía eléctrica.” and in part by the Dirección de Investigaciones de la Universidad Tecnológica de Bolívar under grant PS2020002 associated with the project: “Ubicación óptima de bancos de capacitores de paso fijo en redes eléctricas de distribución para reducción de costos y pérdidas de energía: Aplicación de métodos exactos y metaheurísticos.”

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lotfi, H.; Khodaei, A. AC Versus DC Microgrid Planning. IEEE Trans. Smart Grid 2017, 8, 296–304. [Google Scholar] [CrossRef]
  2. 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]
  3. 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]
  4. Murillo-Yarce, D.; Garcés-Ruiz, A.; Escobar-Mejía, A. Passivity-Based Control for DC-Microgrids with Constant Power Terminals in Island Mode Operation. Rev. Fac. Ing. Univ. Antioq. 2018, 32–39. [Google Scholar] [CrossRef]
  5. Kumar, J.; Agarwal, A.; Agarwal, V. A review on overall control of DC microgrids. J. Energy Storage 2019, 21, 113–138. [Google Scholar] [CrossRef]
  6. Gao, F.; Kang, R.; Cao, J.; Yang, T. Primary and secondary control in DC microgrids: A review. J. Mod. Power Syst. Clean Energy 2018, 7, 227–242. [Google Scholar] [CrossRef] [Green Version]
  7. 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]
  8. 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]
  9. Montoya, O.D.; Garrido, V.M.; Grisales-Norena, L.F.; Gil-Gonzalez, 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. [Google Scholar] [CrossRef]
  10. Montoya, O.D.; Gil-González, W. A MIQP model for optimal location and sizing of dispatchable DGs in DC networks. Energy Syst. 2020. [Google Scholar] [CrossRef]
  11. 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]
  12. 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 Communications in Computer and Information Science; Springer: Berlin/Heidelberg, Germany, 2019; pp. 214–225. [Google Scholar] [CrossRef]
  13. 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. [Google Scholar] [CrossRef]
  14. Altun, T.; Madani, R.; Yadav, A.P.; Nasir, A.; Davoudi, A. Optimal Reconfiguration of DC Networks. IEEE Trans. Power Syst. 2020, 35, 4272–4284. [Google Scholar] [CrossRef]
  15. Montoya, O.D.; Gil-González, W.; Hernández, J.C.; Giral-Ramírez, D.A.; Medina-Quesada, A. A Mixed-Integer Nonlinear Programming Model for Optimal Reconfiguration of DC Distribution Feeders. Energies 2020, 13, 4440. [Google Scholar] [CrossRef]
  16. 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]
  17. Gil-González, W.; Montoya, O.D.; Holguín, E.; Garces, A.; Grisales-Noreña, L.F. Economic dispatch of energy storage systems in dc microgrids employing a semidefinite programming model. J. Energy Storage 2019, 21. [Google Scholar] [CrossRef]
  18. Gil-González, W.; Montoya, O.D.; Grisales-Noreña, L.F.; Cruz-Peragón, F.; Alcalá, G. Economic Dispatch of Renewable Generators and BESS in DC Microgrids Using Second-Order Cone Optimization. Energies 2020, 13, 1703. [Google Scholar] [CrossRef] [Green Version]
  19. Jannesar, M.R.; Sedighi, A.; Savaghebi, M.; Anvari-Moghaddam, A.; Guerrero, J.M. Optimal probabilistic planning of passive harmonic filters in distribution networks with high penetration of photovoltaic generation. Int. J. Electr. Power Energy Syst. 2019, 110, 332–348. [Google Scholar] [CrossRef]
  20. Navidi, M.; Tafreshi, S.M.M.; Anvari-Moghaddam, A. A game theoretical approach for sub-transmission and generation expansion planning utilizing multi-regional energy systems. Int. J. Electr. Power Energy Syst. 2020, 118, 105758. [Google Scholar] [CrossRef]
  21. Khaligh, V.; Anvari-Moghaddam, A. Stochastic expansion planning of gas and electricity networks: A decentralized-based approach. Energy 2019, 186, 115889. [Google Scholar] [CrossRef]
  22. Montoya, O.D.; Serra, F.M.; Angelo, C.H.D. On the Efficiency in Electrical Networks with AC and DC Operation Technologies: A Comparative Study at the Distribution Stage. Electronics 2020, 9, 1352. [Google Scholar] [CrossRef]
  23. Garcés, A.; Montoya, O.D. A Potential Function for the Power Flow in DC Microgrids: An Analysis of the Uniqueness and Existence of the Solution and Convergence of the Algorithms. J. Control. Autom. Electr. Syst. 2019, 30, 794–801. [Google Scholar] [CrossRef]
  24. Gil-González, W.; Molina-Cabrera, A.; Montoya, O.D.; Grisales-Noreña, L.F. An MI-SDP Model for Optimal Location and Sizing of Distributed Generators in DC Grids That Guarantees the Global Optimum. Appl. Sci. 2020, 10, 7681. [Google Scholar] [CrossRef]
  25. Chen, Y.; Xiang, J.; Li, Y. SOCP Relaxations of Optimal Power Flow Problem Considering Current Margins in Radial Networks. Energies 2018, 11, 3164. [Google Scholar] [CrossRef] [Green Version]
  26. Hindi, H. A tutorial on convex optimization. In Proceedings of the 2004 American Control Conference, Boston, MA, USA, 30 June–2 July 2004. [Google Scholar] [CrossRef]
  27. Alizadeh, F.; Goldfarb, D. Second-order cone programming. Math. Program. 2003, 95, 3–51. [Google Scholar] [CrossRef]
  28. Jeyakumar, V.; Li, G. Exact Conic Programming Relaxations for a Class of Convex Polynomial Cone Programs. J. Optim. Theory Appl. 2016, 172, 156–178. [Google Scholar] [CrossRef]
  29. Montoya, O.D.; Molina-Cabrera, A.; Chamorro, H.R.; Alvarado-Barrios, L.; Rivas-Trujillo, E. A Hybrid Approach Based on SOCP and the Discrete Version of the SCA for Optimal Placement and Sizing DGs in AC Distribution Networks. Electronics 2020, 10, 26. [Google Scholar] [CrossRef]
  30. Gil-González, W.; Garces, A.; Montoya, O.D.; Hernández, J.C. A Mixed-Integer Convex Model for the Optimal Placement and Sizing of Distributed Generators in Power Distribution Networks. Appl. Sci. 2021, 11, 627. [Google Scholar] [CrossRef]
  31. Nesterov, Y. Lectures on Convex Optimization; Springer: Berlin/Heidelberg, Germany, 2018. [Google Scholar] [CrossRef]
  32. Tuy, H. Convex Analysis and Global Optimization; Springer: Berlin/Heidelberg, Germany, 2016. [Google Scholar] [CrossRef]
  33. Benson, H.Y.; Ümit, S. 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]
  34. Lobo, M.S.; Vandenberghe, L.; Boyd, S.; Lebret, H. Applications of second-order cone programming. Linear Algebra Appl. 1998, 284, 193–228. [Google Scholar] [CrossRef] [Green Version]
  35. Nesterov, Y.; Nemirovskii, A. Interior-Point Polynomial Algorithms in Convex Programming; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 1994. [Google Scholar] [CrossRef]
  36. Sturm, J.F. Using SeDuMi 1.02, A Matlab toolbox for optimization over symmetric cones. Optim. Methods Softw. 1999, 11, 625–653. [Google Scholar] [CrossRef]
  37. Toh, K.C.; Todd, M.J.; Tütüncü, R.H. SDPT3—A Matlab software package for semidefinite programming, Version 1.3. Optim. Methods Softw. 1999, 11, 545–581. [Google Scholar] [CrossRef]
  38. Grant, M.; Boyd, S. CVX: Matlab Software for Disciplined Convex Programming, Version 2.1. 2014. Available online: http://cvxr.com/cvx (accessed on 9 November 2020).
  39. Lavaei, J.; Low, S.H. Zero Duality Gap in Optimal Power Flow Problem. IEEE Trans. Power Syst. 2012, 27, 92–107. [Google Scholar] [CrossRef] [Green Version]
  40. Ridha, H.M.; Gomes, C.; Hizam, H. Estimation of photovoltaic module model’s parameters using an improved electromagnetic-like algorithm. Neural Comput. Appl. 2020, 32, 12627–12642. [Google Scholar] [CrossRef]
  41. Ćalasan, M.; Mujičić, D.; Rubežić, V.; Radulović, M. Estimation of Equivalent Circuit Parameters of Single-Phase Transformer by Using Chaotic Optimization Approach. Energies 2019, 12, 1697. [Google Scholar] [CrossRef] [Green Version]
  42. Ćalasan, M.; Micev, M.; Ali, Z.M.; Zobaa, A.F.; Aleem, S.H.E.A. Parameter Estimation of Induction Machine Single-Cage and Double-Cage Models Using a Hybrid Simulated Annealing–Evaporation Rate Water Cycle Algorithm. Mathematics 2020, 8, 1024. [Google Scholar] [CrossRef]
  43. Montoya, O.D.; Gil-González, W.; Garces, A. Power flow approximation for DC networks with constant power loads via logarithmic transform of voltage magnitudes. Electr. Power Syst. Res. 2019, 175, 105887. [Google Scholar] [CrossRef]
  44. 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]
Figure 1. Schematic representation of an SOCP constraint Ω = x z . with x R 2 and z R [30].
Figure 1. Schematic representation of an SOCP constraint Ω = x z . with x R 2 and z R [30].
Electronics 10 00176 g001
Figure 2. Exploration of the solution space with a branch and bound method to solve an MI-SOCP formulation with two binary variables.
Figure 2. Exploration of the solution space with a branch and bound method to solve an MI-SOCP formulation with two binary variables.
Electronics 10 00176 g002
Figure 3. Grid configuration of the 21-bus test system.
Figure 3. Grid configuration of the 21-bus test system.
Electronics 10 00176 g003
Figure 4. Grid configuration of the 69-bus test system.
Figure 4. Grid configuration of the 69-bus test system.
Electronics 10 00176 g004
Figure 5. Improvement of the power losses using GAMS solvers and the proposed MI-SOCP approach.
Figure 5. Improvement of the power losses using GAMS solvers and the proposed MI-SOCP approach.
Electronics 10 00176 g005
Figure 6. Daily behavior of the load and renewable generation.
Figure 6. Daily behavior of the load and renewable generation.
Electronics 10 00176 g006
Table 1. Branch and load information of the 21-bus test system.
Table 1. Branch and load information of the 21-bus test system.
Bus iBus j R ij (pu) P j (pu)Bus i (pu)Bus j (pu) R ij (pu) P j (pu)
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. Branch and load information of the 69-bus test system.
Table 2. Branch and load information of the 69-bus test system.
Bus iBus j R ij (pu) P j (pu)Bus i (pu)Bus j (pu) R ij (pu) P j (pu)
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. Comparative results between the proposed MI-SOCP approach and the GAMS package.
Table 3. Comparative results between the proposed MI-SOCP approach and the GAMS package.
21-Bus Test Feeder
SolverLocation (Nodes)Size(pu) p loss (pu)
ALPHAECP 9 , 12 , 16 0.8441 , 1.0254 , 1.4544 0.0306
ANTIGONE 9 , 12 , 16 0.8441 , 1.0254 , 1.4544 0.0306
BARON 9 , 12 , 16 0.8441 , 1.0254 , 1.4544 0.0306
BONMIN 9 , 12 , 16 0.8441 , 1.0254 , 1.4544 0.0306
DICOPT 9 , 12 , 16 0.8441 , 1.0254 , 1.4544 0.0306
SBB 9 , 12 , 16 0.8441 , 1.0254 , 1.4544 0.0306
MI-SOCP 9 , 12 , 16 0.8441 , 1.0254 , 1.4544 0.0306
69-Bus Test Feeder
SolverLocation(Nodes)Size (pu) p loss (pu)
ALPHAECP 16 , 61 , 64 5.0537 , 12.0000 , 5.7763 0.0425
BARON 11 , 61 , 69 8.0194 , 12.0000 , 4.2392 0.0812
BONMIN 17 , 61 , 64 4.9245 , 12.0000 , 5.7944 0.0414
COUENNE 52 , 60 , 61 5.8706 , 6.2416 , 12.0000 0.1216
DICOPT 17 , 61 , 67 5.0311 , 12.0000 , 5.0797 0.0524
SBB 18 , 61 , 63 4.8851 , 12.0000 , 6.0630 0.0458
MI-SOCP 17 , 61 , 64 4.9245 , 12.0000 , 5.7944 0.0414
Table 4. Siting and dimensioning of PV sources in DC grids using the proposed MI-SOCP.
Table 4. Siting and dimensioning of PV sources in DC grids using the proposed MI-SOCP.
No. of PVsLocation (Nodes)Size (kW)En. Losses (kWh/Day)Proc. Time (s)
0 1838.521863.73
1 61 2119.4980 1255.082775.70
2 17 , 61 535.8994 , 2025.81399 1215.216495.22
3 18 , 49 , 61 535.3449 , 864.2630 , 2025.4377 1208.7995233.16
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Molina-Martin, F.; Montoya, O.D.; Grisales-Noreña, L.F.; Hernández, J.C. A Mixed-Integer Conic Formulation for Optimal Placement and Dimensioning of DGs in DC Distribution Networks. Electronics 2021, 10, 176. https://doi.org/10.3390/electronics10020176

AMA Style

Molina-Martin F, Montoya OD, Grisales-Noreña LF, Hernández JC. A Mixed-Integer Conic Formulation for Optimal Placement and Dimensioning of DGs in DC Distribution Networks. Electronics. 2021; 10(2):176. https://doi.org/10.3390/electronics10020176

Chicago/Turabian Style

Molina-Martin, Federico, Oscar Danilo Montoya, Luis Fernando Grisales-Noreña, and Jesus C. Hernández. 2021. "A Mixed-Integer Conic Formulation for Optimal Placement and Dimensioning of DGs in DC Distribution Networks" Electronics 10, no. 2: 176. https://doi.org/10.3390/electronics10020176

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