Next Article in Journal
Practices of Motivators in Adopting Agile Software Development at Large Scale Development Team from Management Perspective
Previous Article in Journal
Anonymous Sealed-Bid Auction on Ethereum
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Optimal Selection and Integration of Batteries in DC Grids through a Mixed-Integer Quadratic Convex Formulation

by
Federico Martin Serra
1,
Oscar Danilo Montoya
2,3,*,
Lázaro Alvarado-Barrios
4,
Cesar Álvarez-Arroyo
5 and
Harold R. Chamorro
6,*
1
Laboratorio de Control Automático (LCA), Facultad de Ingeniería y Ciencias Agropecuarias, Universidad Nacional de San Luis—CONICET, Villa Mercedes, San Luis 5730, Argentina
2
Facultad de Ingeniería, Universidad Distrital Francisco José de Caldas, Bogotá 110231, Colombia
3
Laboratorio Inteligente de Energía, Universidad Tecnológica de Bolívar, Cartagena 131001, Colombia
4
Department of Engineering, Universidad Loyola Andalucía, 41704 Sevilla, Spain
5
Department of Electrical Engineering, Universidad de Sevilla, 41092 Sevilla, Spain
6
Department of Electrical Engineering, KTH Royal Institute of Technology, SE-100 44 Stockholm, Sweden
*
Authors to whom correspondence should be addressed.
Electronics 2021, 10(19), 2339; https://doi.org/10.3390/electronics10192339
Submission received: 19 August 2021 / Revised: 10 September 2021 / Accepted: 22 September 2021 / Published: 24 September 2021
(This article belongs to the Special Issue AI for Cyber-Physical Power Systems Operation and Control)

Abstract

:
This paper deals with the problem of the optimal selection and location of batteries in DC distribution grids by proposing a new mixed-integer convex model. The exact mixed-integer nonlinear model is transformed into a mixed-integer quadratic convex model (MIQC) by approximating the product among voltages in the power balance equations as a hyperplane. The most important characteristic of our proposal is that the MIQC formulations ensure the global optimum reaching via branch & bound methods and quadratic programming since each combination of the binary variables generates a node with a convex optimization subproblem. The formulation of the objective function is associated with the minimization of the energy losses for a daily operation scenario considering high renewable energy penetration. Numerical simulations show the effectiveness of the proposed MIQC model to reach the global optimum of the optimization model when compared with the exact optimization model in a 21-node test feeder. All the validations are carried out in the GAMS optimization software.

1. Introduction

Electrical distribution networks have experienced important paradigm shifts associated with the large-scale insertion of renewable generation based on photovoltaic and wind sources in conjunction with energy storage technologies mainly focused on chemical storage [1,2,3,4]. All of these devices are interfaced with the distribution network through the implementation of power electronic converters with alternating current (AC) and direct current (DC) conversion stages as a function of the device interconnected and the technology of operation of the distribution grid, i.e., AC or DC network [5,6]. From the beginning to the electrical networks to our days, the predominant technology for building electrical networks in transmission levels has been the AC technology due to the most of the demands connected to the networks were simple, i.e., electrical rotate machines (induction motors), temperature conditioners, and illumination, among others; however, the nature of the loads have drastically changed with the appearance of computers, electric vehicles, household appliances, and small dispersed generating and storage systems, most of them operated with DC technologies [7,8,9,10]. Recent studies have demonstrated the advantages of having distribution networks in medium voltage levels operated with DC technologies due to important reductions in energy losses associated with the distribution activity since the reactive power is compensated directly in the point of connection of the load [11,12].
Recognizing the recent relevance that is having distribution networks operated with DC technologies, this research explores the impact of the integration of battery energy storage systems (BESSs) in these grids, considering the high penetration of renewable generation under an economic dispatch environment [13,14]. In the current literature, the problem of the optimal integration of BESSs has been explored in multiple pieces of research for AC and DC grids; here, we present some of these researches.
Authors of [15] have presented a mixed-integer nonlinear programming (MINLP) model to represent the problem of the optimal location of BESSs in AC distribution grids. To solve this model a decomposition method was proposed where the MINLP model is divided into two optimization sub-problems named planning and operation. In the planning stage is defined the optimal location of the BESSs, while the operation stage is entrusted with their optimal operation. For solving the planning stage the classical simulated annealing algorithm is implemented based on sensitivity indexes associated with the impedance matrix of the network. Numerical validations were carried out in two test feeders composed of 135 and 230 nodes; however, the authors did not provide comparisons with other optimization methodologies to confirm the effectiveness of the proposed approach. In Reference [16], the authors have implemented the problem of the optimal location-reallocation of batteries in DC microgrids by solving the exact MINLP model in the GAMS software. Numerical results were presented in a 21-bus system considering that the initial location of the batteries was heuristically defined by the utility company. The authors do not provide any comparative methodology to demonstrate the efficiency of the proposed optimization methodology since the research is presented in a tutorial style. Soroudi, in Reference [13] has presented different optimization models for the optimal operation of BESSs in electrical AC networks. Three models were presented, which correspond to (i) the economic dispatch model, (ii) the DC equivalent of the AC grid, (iii) and the complete AC model of the grid. All these optimization models were solved in the GAMS optimization package; nevertheless, no comparison with metaheuristic or approximated optimization models were provided, since the intention of the author is to provide optimization tools to introduce engineers in the power system optimization topics from the tutorial point of view. Authors of [17] have presented a master-slave optimization methodology to operate batteries in DC networks considering multiple loads curves and high renewable generation availability. In the master stage was proposed a particle swarm optimizer to define the optimal operation of the batteries during the day, while the slave stage was entrusted with solving the multi-period power flow problem. The objective persecuted by the authors corresponded to the minimization of the energy purchasing cost at the substation node. The proposed methodology was tested in a test feeder composed of 21 nodes, and the efficiency of the methodology was compared with different metaheuristic approaches such as black-hole optimizer and the genetic algorithm. In Reference [18] the authors have proposed the implementation of the genetic algorithm to select and operate BESSs in AC distribution networks. The genetic algorithms were entrusted with determining the size and the operation scheme of the batteries using an integer codification with three possible states to operate these batteries. Numerical results demonstrate that the total grid energy losses is reduced when batteries are installed by using the Baran & Wu test feeder composed of 69 nodes [19]; however, the authors did not provide comparisons with exact or metaheuristics optimizers to confirm the effectiveness of the proposed optimization approach.
Other authors have proposed multiple operative models to coordinate the daily operation of the batteries; some of these approaches are: mixed-integer linear programming [20,21,22]; second order cone optimization [23,24,25], semidefinite programming [26]; genetic algorithms [27,28,29], particle swarm optimization [30,31]; nonlinear programming [32,33,34,35,36], and reinforcement learning for energy system optimization [37,38]. The main characteristic of those researches is that the batteries are modeled through a linear relation between the state-of-charge and the amount of power injected/absorbed into the grid [11]; this linear representation allows solving efficiently the problem of the optimal dispatch of these batteries in AC and/or DC grids where these are previously located to the network.
To contribute to the research area associated with the optimal integration and operation of BESSs in electrical networks, here, we propose an efficient mixed-integer quadratic convex (MIQC) model to select and size batteries in DC networks. This corresponds to an improvement of the mixed-integer nonlinear programming (MINLP) model proposed in [16], with the main advantage that the optimal location and coordination of the batteries correspond to the global optimal solution of the problem since the MIQC model ensures the finding through the application of the Branch & Bound and interior-point methods [39]. To verify the efficiency of the proposed MIQC model is evaluated in a DC network composed of 21 nodes and its results are compared with the exact MINLP model implemented in the GAMS software.
To demonstrate the novelty of the proposed convex reformulation to select and locate batteries in DC grids, in Table 1 are summarized the main literature reports, by highlighting the type of mathematical models, solution techniques, and objective functions considered.
Note that the literature reports in Table 1 show that the most common objective functions are associated with the minimization of the grid operating costs and energy losses costs, which supports the selection of the objective function considered in this research to select and integrate BESSs in DC grids. In addition, the proposed MIQC reformulation for the exact MINLP model has not been previously proposed for AC and DC grids, which was identified in this research as a gap in the scientific literature that this investigation tries to fill.
Ahead, this document has the following organization: Section 2 presents the exact MINLP formulation of the problem of the optimal selection and operation of BESSs in DC grids considering as the objective function minimization of the energy losses costs during the period of operation; Section 3 presents the proposed MIQC reformulation by using a Taylor linearization of the product among voltages in the power balance equation; Section 4 presents the main characteristics of the 21-bus system used to validate the proposed optimization methodology. Section 5 presents the main numerical results when comparing the exact MINLP and the proposed MIQC models, including their analysis and discussion. Finally, Section 6 lists the main conclusions obtained from this research.

2. General Formulation

The study of the optimal siting and selection of BESS in DC microgrids considering high penetration of renewable sources corresponds to a mixed-integer nonlinear programming formulation that can be represented as a multi-period economic dispatch [15,41]. The nonlinear part of the optimization model is defined by the product among voltages in the power balance constraint [42]; while binary nature is associated with the variables that define the location or not of a BESS in an arbitrary node of the network [16].

2.1. Objective Function

Here, we present the exact MINLP for optimal locating-reallocating BESS in DC networks considering as the objective function the minimization of the costs of the daily energy losses.
min f 1 = C o E t H i B v i , t j B G i j v j , t Δ t ,
where f 1 is the value of the objective function related with the total costs of the daily energy losses, C o E t is the average cost of the energy in the spot market; v i , t and v j , t are the voltage variables at nodes i and j during the period of time t, respectively. G i j corresponds to the conductance parameter that relates nodes i and j; Δ t represents the length of the time period where the demand and generation are assumed constant; this parameter can be 1 h, 30 min, or 15 min, depending on the data resolution. It should be noted that H and B are the sets that contain all the periods of time and all the buses of the network, respectively.
Remark 1.
The main characteristic of the objective function f 1 is that it corresponds to an algebraic sum of quadratic terms with the main advantage that its combination is convex, since the conductance matrix, i.e., G , is a positive definite [43]; which implies that this objective function can be rewritten as follows
f 1 = C o E t H v t T G v t Δ t ,
where v t is the vector that contains all the nodal voltages at each period of time.
It is worth mentioning that the minimization of the objective function (2) involves the improvement of the voltage profiles in all the nodes of the networks since the total grid energy losses are a nonlinear function of the voltage profiles. In addition, due to the presence of the conductance matrix in this objective function, the usage of the renewable energies are in general maximized, since these helps with local power injection that reduces the current magnitudes provided by the slack source. For this reason, as confirmed by Table 1, in this research the minimization of the total costs of the energy losses is selected as the performance index for our proposed MIQC model.

2.2. Set of Constraints

The problem of the optimal selection and operation BESSs in DC grids includes multiple constraints: power balance, state-of-charge (SoC) in batteries, devices’ capabilities, voltage regulation bounds, and the maximum number of BESSs that can be installed along the grid [17], among others. In this paper, the complete list of constraints considered is presented below.
p i , t + p i , t d g + b E p i , t b p i , t d = v i , t j B G i j v j , t , { i B & t H }
S o C i , t b = S o C i , t 1 b φ i b p i , t b Δ t , { b E , i B & t H }
S o C i , t 0 b = x i b S o C i b , i n i , { b E & i B }
S o C i , t f b = x i b S o C i b , f i n , { b E & i B }
p i , t min p i , t p i , t max , { i B & t H }
p i , t d g , min p i , t d g p i , t d g , max , { i B & t H }
x i b p i b , min p i , t b x i b p i b , max , { b E , i B & t H }
v i min v i , t v i max , { i B & t H }
x i b S o C i b , min S o C i , t b x i b S o C i b , max , { b E , i B & t H }
b E i B x i b = N b max ,
where p i , t , p i , t d g , p i , t b , and p i , t d are the power produced by the slack and disperse sources, the power absorbed/provided by the BESSs, and the constant power consumption at node i in the time period t, respectively; S o C i , t b is the BESS’ state-of-charge which is connected at node i in the period of time t; x i b represents a binary variable associated with the locating possibility of the BESS b at node i; S o C i b , i n i and S o C i b , f i n are the initial and final projected states of charge of the batteries, respectively, while S o C i b , min and S o C i b , max represent the minimum and maximum limits of the state-of-charge variables; p i , t min , p i , t max , p i , t d g , min , and p i , t d g , max correspond to the minimum and maximum limits of admissible generation for conventional and disperse generation sources located connected at node i in time period of time t, respectively, while p i b , min and p i b , max define the minimum and maximum charge/discharge capacities for the battery located at node i; v i min and v i max define the voltage regulation limits of the DC network. Finally, φ i b is the charge coefficient associated with the BESS connected at node i. Note that N b max is the maximum number of BESS available for being introduced in the DC grid, and E is the set that includes all the BESS technologies available.

2.3. Model Interpretation

The exact MINLP formulation (1) to (12) can be understood as follows: Equation (1) and its matricial form defined in (2) correspond to the objective function of the optimization problem that defines the daily cost of the energy losses of the network associated with the energy dissipation in all the distribution lines of the distribution system. Equality constraint (3) defines the power balance at each node of the network which results from the application of the nodal voltage method to the grid with constant power consumption; Equation (4) defines the linear relation between the state-of-charge of the battery and its power delivered/consumed [13]; Equations (5) and (6) define the operative characteristics for operating BEESs, which are the desired initial and final S o C s typically defined by the distribution company. Equations (7)–(9) define the upper and lower power limits for the slack and disperse generation sources as well as for BESSs, respectively. In Equation (10), the voltage regulation constraints are set. Inequality constraint (11) presents the upper and lower bounds for the S o C s in BESSs; Equation (12) defines the maximum number of BESSs available for installation in the DC distribution grid.
Remark 2.
The main complication of the MINLP model (1)–(12) corresponds to the power equilibrium constraints (see Equation (3)) since this is non-convex owing the products among continuous variables, i.e., voltages in all the nodes [42].
This set of equality constraints (power balance equations) will be treated with a Taylor-based linearization to become the MINLP model into a MIQC model as the main contribution of this work as will be presented in the next section.

3. Mixed-Integer Quadratic Reformulation

To deal with the MINLP model that represents the problem of the optimal selection and location of batteries in DC distribution grids, we present the proposed MIQC reformulation based on the linearization of the product among two continuous variables present in the power balance constraint (3). This linearization is based on the convex representation of the load flow problem for DC distribution grids proposed in [44]. The linearization of the product of two continuous variables is detailed as follows.
Let us suppose that a function to represent the product between two continuous positive variables defined as:
f ω 1 , ω 2 = ω 1 ω 2 ,
In addition, consider that the initial value assignable to these variables are ω 10 and ω 20 , which implies that after applied the Taylor’ series expansion to (13), the following representation for this product between two continuous variables is obtained,
f ω 1 , ω 2 = ω 20 ω 1 + ω 10 ω 2 ω 10 ω 20 + f H . O . T ω 1 , ω 2 , ω 10 , ω 20 ,
where f H . O . T · models the high-order-terms in the Taylor’s series expansion. To obtain a linear approximation of the product between two variables, the high-order-terms in f H . O . T · can be neglected due to their small contribution around the operating point ω 10 , ω 20 .
When the linear approximation defined in (14) is applied to electrical networks to transform the power balance set of constraints in (3), the following linear equivalent set of constraints is obtained [44].
p i , t + p i , t d g + b E p i , t b p i , t d = j B G i j v i , t v j 0 , t + v i 0 , t v j , t v i 0 , t v j 0 , t , { i B & t H }
where v i 0 , t and v j 0 , t represent the linearization points for the voltage profiles at each period of time; where if the per-unit representation is used, then these values are equal to 1.00 pu [44].
To demonstrate the effectiveness of the linearization of the product between voltages around the operative point v 10 , v 20 = 1 , 1 a graphic comparison between the nonlinearised and the linearised representation is considered for two small load systems. Considering a grid with two demand nodes with lower and upper voltage regulation bounds of 0.90 pu and 1.10 pu, the error between the nonlinear function in (13) and the linear representation (14) is presented in Figure 1.
Note that the percentage error depicted in Figure 1 shows that the linearization of the product between voltages has an estimation error of about 1 % in the extreme voltage points, which confirms that the linear representation in (14) is suitable to represent the product of the voltage variables as was redefined in Equation (15).
Remark 3.
Once the power balance equations are linearized as defined in (15), the complete structure of the proposed MIQC is composed by the objective function (1) or (2) and the restrictions (4)–(12) and (15).
To summarize the solution methodology for the proposed MIQC model and the exact MINLP formulation in Figure 2 is presented the flow diagram of both mathematical models in the GAMS software [13].

4. Analysis of a DC Network

To validate the MIQC model for the optimal selection and location of BESSs in DC networks using the GAMS software we use the 21-node test feeder reported in [16]. The complete information of this test feeder is presented below.
The 21-bus system corresponds to a radial DC network with 21 nodes and 20 branches with a radial configuration. A controlled voltage source is sited at node 1, which defines the voltage profile of the network in 1 kV. The connection among nodes in this test feeder is presented in Figure 3.
The information regarding constant power consumption, energy purchasing cost, variation of the demand along the day, and the location of the BESSs in this test feeder is reported in Table 2, Table 3 and Table 4.
Two distributed generators (DGs) are considered in the the 21-bus system, one wind power generator and one photovoltaic system. The wind power generator is located at node 12 with a maximum power rate of 221.52 kW. The photovoltaic source is located at node 21 with a maximum power rate of 281.58 kW. It is worth mentioning that the rated power of the DGs will multiply the normalized generation curves presented in Table 5.

5. Implementation and Results

The implementation of the proposed MIQC model and the exact MINLP model was made on a desktop computer running on INTEL(R) Core(TM) i7-7700, 3.60 GHz, 8 GB RAM with 64-bit Windows 10 Pro (Intel, Santa Clara, CA, USA). The optimization package corresponds to the GAMS software version 25.1.3 using the BONMIN solver.
To evaluate the performance of the BESSs in the DC network in all the simulations is considered that the initial and the final state-of-charge are assigned as 50 % , and along the day this variable can vary from 10 % to 90 % [36]. In addition, three simulation cases are studied as follows:
  • Case 1: The initial location of the BESSs is tested in the exact MINLP model and the MIQC model to determine the error introduced by the Taylor approximation in the daily cost of the energy losses. Note that this simulation case solves the power flow problem with multiple periods since the binary variables related with BESSs are fixed, i.e., the MINLP model becomes into a nonlinear programming model and the MIQC model becomes a quadratic convex model.
  • Case 2: The selection and location of the BESSs is determined by solving the exact MINLP model.
  • Case 3: The selection and location of the BESSs is defined by the solution of the MIQC model.
The main idea of the proposed simulation scenarios is to verify the effect of the proposed approximated MIQC model to select, locate and operate batteries in DC networks when compared with the exact MINLP model. For this reason, both models are implemented in the same optimization environment (i.e., GAMS interface). Regarding processing times, both optimization models take less than 5 min to be solved, which implies that for the 21-node test system the proposed optimization model allows evaluating multiple operative conditions (generation and demand combinations) to define the best operative scheme of the batteries, once these batteries have been installed, by ensuring that each solution will be optimal.

5.1. Comparative Results in the Case 1

In this simulation case, the initial BESS location is presented in Table 4 and Figure 3 are fixed into the exact MINLP and the MIQC model.
Once the exact MINLP model is solved in the GAMS package with the DICOPT solver, it is found a total daily cost of the energy losses of COP $ 52957.92 ; while the solution of the MIQC model with the same solver provides an optimal solution with a cost of COP $ 50 , 890.12 . The difference between both solutions is about 3.90 % , which corresponds to the approximation error between the exact and approximated power balance equations (see Equations (3) and (15)). In Figure 4 are reported the profiles of the state-of-charge of the batteries obtained with the exact and the convex approximation models.
The behavior of the state of the charge in the exact and convex model follows the same tendency with negligible errors, which confirms the effectiveness of the quadratic convex approximated model to operate batteries in DC networks. The main advantage of the quadratic approximation is that the existence of the global optimal solution is ensured via convexity theory, which is indeed more attractive for the grid operation since the solution with the same inputs will always be the same, which is not possible to ensure with nonlinear programming models or solution techniques based on metaheuristics. It is worth motioning that if the solution provided by the proposed quadratic model is evaluated in the exact nonlinear model, then, the error in the estimation of the daily energy losses costs is less than 0.5 % , which demonstrate the effectiveness of our proposal to determine the operative plan of BESSs under an economic dispatch environment for DC networks.
With respect to the behavior of the state of charge for all the BESSs depicted in Figure 4, we can observe that: (i) all the batteries begin and finish the day with the operative consign assigned, i.e., 50 % of the state of charge; (ii) the maximum value of the state of charge is about 77.20 % for the battery located at node 7, and the minimum value occurs for the same battery with a value of 48.90 % ; these values imply that the batteries maintain during all the day values between their minimum and maximum bounds, i.e., from 10 % to 90 % ; and (iii) the behavior of the state of charge at each battery is different since this depends on its location and the possibility of absorbing energy from the renewable energy source to return this energy to the grid when the load increases. Note for example that the BESS located at node system provides energy to the grid after period 21, while the remainder BESSs returns energy to the grid only after period 33, these behaviors confirm the complex relation among demand, power generation, and batteries to reach the minimization of the daily energy losses costs in the DC grid.

5.2. Comparative Results between Case 2 and Case 3

To determine the effectiveness and robustness of the proposed MIQC model to select and locate BESSs in DC networks, here, we compared the solution provided with the exact MINLP formulation and our proposed convex approximation. The solution of the exact MINLP and the MIQC models are obtained with the GAMS optimization package using the BONMIN solver. The objective function value for both models as well as the battery locations are reported in Table 6.
Results in Table 6 show that: (i) the BESSs’ location found by the exact MINLP model are the nodes 13, 20 and 21, with an objective function of COP $47,209.95, i.e., a reduction with respect to the base case of 10.85 % . The solution of the exact MINLP model reallocates all the batteries with respect to the benchmark case, which implies that the solution of the optimization model to select and locate the BESSs in the DC network is better than the heuristic approach reported in [36]; (ii) the proposed MIQC model finds a better optimal solution with an objective value of COP $41,627.34, which has an estimation error with respect to the evaluation of the BESSs’ location in the exact model about 3.49 % , which implies that the location of the batteries in nodes 5, 16, and 21 provides a daily energy losses cost of COP $43,134.59 (note that this objective function value corresponds to the evaluation of the batteries’ location provided by our proposed MIQC model into the exact MINLP approach to eliminate the estimation error introduced with the model linearization). This objective function shows that with respect to the benchmark case the effective reduction of the daily energy losses cost is about 18.55 % when the MIQC model is used to select and locate the BESSs in the DC grid; and (iii) the effective improvement of the proposed MIQC model with respect to the exact MINLP formulation is about 8.63 % .
Figure 5 presents the final location of the BESSs system obtained for the exact MINLP and the proposed MIQC model, where we can observe that only node 21 appears in both solutions, this is due to the presence of a distributed generator in this node that allows storing energy in the periods of high generation and low demand. To illustrate the behavior of the BESSs in both models, we also present the state-of-charge profile in this node for both solutions compared with the generation curve in this node.
Figure 6 shows that: (i) both optimization models present a similar state-of-charge behavior at node 21; however, the proposed MIQC model provides additional energy to the grid which helps with the total energy losses cost minimization (see periods of time between 3 to 33) and (ii) both batteries provide energy in the initial period of time (from 1 to 18) since this energy will recover in the period of times where the PV source increases its power injection (see periods of time between 15 to 33) to help with additional power injections in the periods of time that the PV source decreases (periods 33 to 48) and the demand increase, i.e., periods of time from 33 and beyond.

6. Conclusions

This paper studied the problem of the optimal selection and location of the BESSs in DC grids by transforming the exact MINLP formulation into a MIQC equivalent. This transformation was applied to the power balance constraint using Taylor’s series expansion applied to the product among voltages at each node. The proposed obtained MIQC model has as the main advantage that it ensures the global optimal reaching via Branch & Bound and interior-point methods. Numerical results demonstrated that the error introduced by the Taylors’ approximation is less than 4.0 % for in the economic dispatch analysis when batteries are considered fixed.
When the exact MINLP and the proposed MIQC models are solved in the GAMS software with the BONMIM solver the numerical results demonstrated that the MINLP model obtains a reduction of 10.85 % in the daily energy losses cost when compared to the benchmark case by reallocating the BESSs in nodes 13 ( A ) , 20 ( B ) , 21 ( B ) ; however, the proposed MIQC model allowed a reduction of 18.55 % in the daily grid operation costs by reassigning batteries at nodes 5 ( A ) , 16 ( B ) , 21 ( B ) . The net difference between both models was about 7.70 % in favor of the proposed MIQC model when compared with the exact MINLP formulation, which confirmed the robustness and effectiveness of the proposed convex formulation to select and locate BESSs in DC networks.
It will be possible to develop the following future researches: (i) to include in the proposed MIQC model the possibility of doing the location of renewable energy sources and batteries for DC networks simultaneously; (ii) reformulate the exact MINLP model as a mixed-integer conic model to size and locate batteries in DC grids by minimizing the total cost estimation error; and (ii) improve the battery models with temperature and degradation factors using reinforcement learning techniques for being integrated into AC and DC networks considering different objective functions such as voltage maximization or renewables usage maximization.

Author Contributions

Conceptualization, F.M.S., O.D.M., L.A.-B., C.Á.-A. and H.R.C.; Methodology, F.M.S., O.D.M., L.A.-B., C.Á.-A. and H.R.C.; Investigation, F.M.S., O.D.M., L.A.-B., C.Á.-A. and H.R.C.; Writing—review and editing, F.M.S., O.D.M., L.A.-B., C.Á.-A. and H.R.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partially supported in part by the Laboratorio de Simulación Hardware-in-the-loop para Sistemas Ciberfísicos de la Universidad Loyola Andalucía.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

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. Guerrero, J.; Blaabjerg, F.; Zhelev, T.; Hemmes, K.; Monmasson, E.; Jemei, S.; Comech, M.; Granadino, R.; Frau, J. Distributed Generation: Toward a New Energy Paradigm. IEEE Ind. Electron. Mag. 2010, 4, 52–64. [Google Scholar] [CrossRef]
  2. Saberi, H.; Nazaripouya, H.; Mehraeen, S. Implementation of a Stable Solar-Powered Microgrid Testbed for Remote Applications. Sustainability 2021, 13, 2707. [Google Scholar] [CrossRef]
  3. Akinyele, D.; Rayudu, R. Review of energy storage technologies for sustainable power networks. Sustain. Energy Technol. Assess. 2014, 8, 74–91. [Google Scholar] [CrossRef]
  4. Pandey, G.; Singh, S.N.; Rajpurohit, B.S.; Gonzalez-Longatt, F.M. Smart DC Grid for Autonomous Zero Net Electric Energy of Cluster of Buildings. IFAC-PapersOnLine 2015, 48, 108–113. [Google Scholar] [CrossRef]
  5. Blaabjerg, F.; Yang, Y.; Ma, K.; Wang, X. Advanced Grid Integration of Renewables Enabled by Power Electronics Technology. In Nachhaltige Energieversorgung und Integration von Speichern; Springer Fachmedien: Wiesbaden, Germany, 2015; pp. 3–9. [Google Scholar] [CrossRef] [Green Version]
  6. Carrasco, J.; Franquelo, L.; Bialasiewicz, J.; Galvan, E.; PortilloGuisado, R.; Prats, M.; Leon, J.; Moreno-Alfonso, N. Power-Electronic Systems for the Grid Integration of Renewable Energy Sources: A Survey. IEEE Trans. Ind. Electron. 2006, 53, 1002–1016. [Google Scholar] [CrossRef]
  7. Wehbring, N.; Bleilevens, R.; Tepasse, B.; Priebe, J.; Moser, A. Strategies to convert AC into DC Medium Voltage Grids. In Proceedings of the 2018 53rd International Universities Power Engineering Conference (UPEC), Glasgow, UK, 4–7 September 2018. [Google Scholar] [CrossRef]
  8. Shaqsi, A.Z.A.; Sopian, K.; Al-Hinai, A. Review of energy storage services, applications, limitations, and benefits. Energy Rep. 2020, 6, 288–306. [Google Scholar] [CrossRef]
  9. Gonzalez-Longatt, F.; Sanchez, F.; Singh, S.N. On the topology for a smart direct current microgrid for a cluster of zero-net energy buildings. In Distributed Energy Resources in Microgrids; Elsevier: Amsterdam, The Netherlands, 2019; pp. 455–481. [Google Scholar] [CrossRef]
  10. Wei, Z.; Quan, Z.; Wu, J.; Li, Y.; Pou, J.; Zhong, H. Deep Deterministic Policy Gradient-DRL Enabled Multiphysics-Constrained Fast Charging of Lithium-Ion Battery. IEEE Trans. Ind. Electron. 2021, in press. [Google Scholar] [CrossRef]
  11. 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]
  12. Stieneker, M.; Doncker, R.W.D. Medium-voltage DC distribution grids in urban areas. In Proceedings of the 2016 IEEE 7th International Symposium on Power Electronics for Distributed Generation Systems (PEDG), Vancouver, BC, Canada, 27–30 June 2016. [Google Scholar] [CrossRef]
  13. Soroudi, A. Power System Optimization Modeling in GAMS; Springer: Berlin/Heidelberg, Germany, 2017. [Google Scholar] [CrossRef]
  14. Zocca, A.; Zwart, B. Minimizing heat loss in DC networks using batteries. In Proceedings of the 2016 54th Annual Allerton Conference on Communication, Control, and Computing (Allerton), Monticello, IL, USA, 27–30 September 2016. [Google Scholar] [CrossRef] [Green Version]
  15. Valencia, A.; Hincapie, R.A.; Gallego, R.A. Optimal location, selection, and operation of battery energy storage systems and renewable distributed generation in medium–low voltage distribution networks. J. Energy Storage 2021, 34, 102158. [Google Scholar] [CrossRef]
  16. Montoya, O.D.; Gil-González, W.; Rivas-Trujillo, E. Optimal Location-Reallocation of Battery Energy Storage Systems in DC Microgrids. Energies 2020, 13, 2289. [Google Scholar] [CrossRef]
  17. Grisales-Noreña, L.F.; 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]
  18. Grisales-Noreña, L.F.; Montoya, O.D.; Gil-González, W. Integration of energy storage systems in AC distribution networks: Optimal location, selecting, and operation approach based on genetic algorithms. J. Energy Storage 2019, 25, 100891. [Google Scholar] [CrossRef]
  19. Vai, V.; Suk, S.; Lorm, R.; Chhlonh, C.; Eng, S.; Bun, L. Optimal Reconfiguration in Distribution Systems with Distributed Generations Based on Modified Sequential Switch Opening and Exchange. Appl. Sci. 2021, 11, 2146. [Google Scholar] [CrossRef]
  20. Cardoso, G.; Stadler, M.; Siddiqui, A.; Marnay, C.; DeForest, N.; Barbosa-Póvoa, A.; Ferrão, P. Microgrid reliability modeling and battery scheduling using stochastic linear programming. Electr. Power Syst. Res. 2013, 103, 61–69. [Google Scholar] [CrossRef] [Green Version]
  21. Zhang, Y.; Jia, Q.S. Operational Optimization for Microgrid of Buildings with Distributed Solar Power and Battery. Asian J. Control 2017, 19, 996–1008. [Google Scholar] [CrossRef]
  22. Luna, A.C.; Diaz, N.L.; Graells, M.; Vasquez, J.C.; Guerrero, J.M. Mixed-Integer-Linear-Programming-Based Energy Management System for Hybrid PV-Wind-Battery Microgrids: Modeling, Design, and Experimental Verification. IEEE Trans. Power Electron. 2017, 32, 2769–2783. [Google Scholar] [CrossRef] [Green Version]
  23. Home-Ortiz, J.M.; Pourakbari-Kasmaei, M.; Lehtonen, M.; Mantovani, J.R.S. Optimal location-allocation of storage devices and renewable-based DG in distribution systems. J. Energy Storage 2019, 172, 11–21. [Google Scholar] [CrossRef]
  24. Macedo, L.H.; Ortega-Vazquez, M.A.; Romero, R. Optimal operation of storage systems in distribution networks considering battery degradation. In Proceedings of the 2018 Simposio Brasileiro de Sistemas Eletricos (SBSE), Niteroi, Brazil, 12–16 May 2018. [Google Scholar] [CrossRef]
  25. 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]
  26. 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, 1–8. [Google Scholar] [CrossRef]
  27. Mansuwan, K.; Jirapong, P.; Burana, S.; Thararak, P. Optimal Planning and Operation of Battery Energy Storage Systems in Smart Grids Using Improved Genetic Algorithm Based Intelligent Optimization Tool. In Proceedings of the 2018 International Conference and Utility Exhibition on Green Energy for Sustainable Development (ICUE), Phuket, Thailand, 24–26 October 2018. [Google Scholar] [CrossRef]
  28. Magnor, D.; Sauer, D.U. Optimization of PV Battery Systems Using Genetic Algorithms. Energy Procedia 2016, 99, 332–340. [Google Scholar] [CrossRef] [Green Version]
  29. Boonluk, P.; Siritaratiwat, A.; Fuangfoo, P.; Khunkitti, S. Optimal Siting and Sizing of Battery Energy Storage Systems for Distribution Network of Distribution System Operators. Batteries 2020, 6, 56. [Google Scholar] [CrossRef]
  30. Ikeda, S.; Ooka, R. Metaheuristic optimization methods for a comprehensive operating schedule of battery, thermal energy storage, and heat source in a building energy system. Appl. Energy 2015, 151, 192–205. [Google Scholar] [CrossRef]
  31. Kai, H.; Yong-Fang, G.; Zhi-Gang, L.; Hsiung-Cheng, L.; Ling-Ling, L. Development of Accurate Lithium-Ion Battery Model Based on Adaptive Random Disturbance PSO Algorithm. Math. Probl. Eng. 2018, 2018, 3793492. [Google Scholar] [CrossRef]
  32. Grisales, L.F.; Grajales, A.; Montoya, O.D.; Hincapie, R.A.; Granada, M.; Castro, C.A. Optimal location, sizing and operation of energy storage in distribution systems using multi-objective approach. IEEE Lat. Am. Trans. 2017, 15, 1084–1090. [Google Scholar] [CrossRef]
  33. Aaslid, P.; Geth, F.; Korpås, M.; Belsnes, M.M.; Fosso, O.B. Non-linear charge-based battery storage optimization model with bi-variate cubic spline constraints. J. Energy Storage 2020, 32, 101979. [Google Scholar] [CrossRef]
  34. Montoya, O.D.; Gil-González, W. Dynamic active and reactive power compensation in distribution networks with batteries: A day-ahead economic dispatch approach. Comput. Electr. Eng. 2020, 85, 106710. [Google Scholar] [CrossRef]
  35. Berglund, F.; Zaferanlouei, S.; Korpås, M.; Uhlen, K. Optimal Operation of Battery Storage for a Subscribed Capacity-Based Power Tariff Prosumer—A Norwegian Case Study. Energies 2019, 12, 4450. [Google Scholar] [CrossRef] [Green Version]
  36. Montoya, O.D.; Gil-González, W.; Grisales-Noreña, L.; Orozco-Henao, C.; Serra, F. Economic Dispatch of BESS and Renewable Generators in DC Microgrids Using Voltage-Dependent Load Models. Energies 2019, 12, 4494. [Google Scholar] [CrossRef] [Green Version]
  37. Wu, J.; Wei, Z.; Liu, K.; Quan, Z.; Li, Y. Battery-Involved Energy Management for Hybrid Electric Bus Based on Expert-Assistance Deep Deterministic Policy Gradient Algorithm. IEEE Trans. Veh. Technol. 2020, 69, 12786–12796. [Google Scholar] [CrossRef]
  38. Wu, J.; Wei, Z.; Li, W.; Wang, Y.; Li, Y.; Sauer, D.U. Battery Thermal- and Health-Constrained Energy Management for Hybrid Electric Bus Based on Soft Actor-Critic DRL Algorithm. IEEE Trans. Ind. Inform. 2021, 17, 3751–3761. [Google Scholar] [CrossRef]
  39. Benson, H.Y.; Sağlam, Ü. Mixed-Integer Second-Order Cone Programming: A Survey. In Theory Driven by Influential Applications; INFORMS: Holland, MI, USA, 2013; pp. 13–36. [Google Scholar] [CrossRef] [Green Version]
  40. Molina-Martin, F.; Montoya, O.D.; Grisales-Noreña, L.F.; Hernández, J.C.; Ramírez-Vanegas, C.A. Simultaneous Minimization of Energy Losses and Greenhouse Gas Emissions in AC Distribution Networks Using BESS. Electronics 2021, 10, 1002. [Google Scholar] [CrossRef]
  41. Dos Santos, C.; Cavalheiro, E.; Bartmeyer, P.; Lyra, C. A MINLP Model to Optimize Battery Placement and Operation in Smart Grids. In Proceedings of the 2020 IEEE Power & Energy Society Innovative Smart Grid Technologies Conference (ISGT), Washington, DC, USA, 17–20 February 2020. [Google Scholar] [CrossRef]
  42. 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]
  43. 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]
  44. 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]
Figure 1. Percentage of error when compared the convex approximation and the exact product of two positive variables around the operative point ( x 0 , y 0 ) = ( 1 , 1 ) .
Figure 1. Percentage of error when compared the convex approximation and the exact product of two positive variables around the operative point ( x 0 , y 0 ) = ( 1 , 1 ) .
Electronics 10 02339 g001
Figure 2. Procedure to solve an optimization model in GAMS.
Figure 2. Procedure to solve an optimization model in GAMS.
Electronics 10 02339 g002
Figure 3. Electrical topology for the 21-bus system.
Figure 3. Electrical topology for the 21-bus system.
Electronics 10 02339 g003
Figure 4. Profile of the state of charge in all the batteries with the exact and convex quadratic models: (a) S o C profile at node 7, (b) S o C profile at node 10, and (c) S o C profile at node 15.
Figure 4. Profile of the state of charge in all the batteries with the exact and convex quadratic models: (a) S o C profile at node 7, (b) S o C profile at node 10, and (c) S o C profile at node 15.
Electronics 10 02339 g004
Figure 5. BESSs’ locations obtained with the exact MINLP and the proposed MIQC model.
Figure 5. BESSs’ locations obtained with the exact MINLP and the proposed MIQC model.
Electronics 10 02339 g005
Figure 6. Behavior of the state-of-charge at node 21 for the exact MINLP and the proposed MIQC model compared with the PV generation in the same node.
Figure 6. Behavior of the state-of-charge at node 21 for the exact MINLP and the proposed MIQC model compared with the PV generation in the same node.
Electronics 10 02339 g006
Table 1. Main approaches reported in the specialized literature.
Table 1. Main approaches reported in the specialized literature.
Math. ModelObjective FunctionSolution MethodRef.
MINLPMinimization of the grid generation costsGeneral algebraic modeling system[34]
MINLPMinimization of energy losses costs and investment costsSensitivity index combined with simulated annealing[15]
MINLPMinimization of energy losses costsGeneral algebraic modeling system[16]
MINLPMinimization of energy losses costsGenetic algorithms and multiperiod power flow[18]
NLPSimultaneous minimization of energy losses costs and greenhouse gas emissionsGeneral algebraic modeling system[40]
NLPMinimization of grid generation costsGeneral algebraic modeling system[13]
NLPMinimization of energy losses costsParticle swarm optimization and multiperiod power flow[17]
LPMinimization of grid operation costs and greenhouse gas emissionsStochastic linear programming[20]
MILPMinimization of the operating costs reduction by promoting self-consumptionGeneral algebraic modeling system[22]
MILPMinimization of operative costs in microgridsSimulation scenarios in the CPLEX solver[21]
MICPMinimization of the grid expansion planning costsCPLEX solver in the AMPL software[23]
SDPMinimization of the grid generation costsConvex solvers in the CVX environment for MATLAB[26]
SOCPMinimization of the grid generation costsConvex solvers in the CVX environment for MATLAB[25]
Table 2. Data for the 21-bus system.
Table 2. Data for the 21-bus system.
From iTo j R i j (pu) P j (pu)From iTo j R i j (pu) P j (pu)From iTo j R i j (pu) P j (pu)
1 (slack)20.00530.70790.00720.8015160.00640.23
130.00540.003100.00530.0016170.00740.43
340.00540.3610110.00380.4516180.00810.34
450.00630.0411120.00790.6814190.00780.09
460.00510.3611130.00780.1019200.00840.21
370.00370.0010140.00830.0019210.00810.21
780.00790.3214150.00650.22
Parameters in per-unit representation using the following bases P base = 100 kW and V base = 1 kV.
Table 3. Hourly demand behavior and energy buying costs.
Table 3. Hourly demand behavior and energy buying costs.
Time (h) C o E (pu)Dem. Var. (%)Time (h) C o E (pu)Dem. Var. (%)Time (h) C o E (pu)Dem. Var. (%)
0.50.8105348.50.92636216.50.973790
1.00.7789289.00.94216817.0190
1.50.7474229.50.95797217.50.994790
2.00.73682210.00.95797818.00.989590
2.50.72632210.50.95798418.50.973786
3.00.73162011.00.95798619.00.957984
3.50.73681811.50.95799019.50.952692
4.00.74741812.00.95269220.00.9474100
4.50.75791812.50.94749420.50.921198
5.00.80002013.00.94749421.00.894794
5.50.84212213.50.94219021.50.868490
6.00.87892614.00.93688422.00.842184
6.50.91582814.50.94218622.50.794776
7.00.93683415.00.94749023.00.747468
7.50.95794015.50.94749023.50.721158
8.00.94215016.00.94749024.00.694750
The cost of the energy is assumed as COP$/kWh 479.3389 (CODENSA utility in Colombia [16].
Table 4. BESSs information.
Table 4. BESSs information.
NodeType φ b p b , max p b , min
7A0.06254−3.2
10B0.08133.2−2.4616
15B0.08133.2−2.4616
Table 5. Normalized Curve of generation: Wind power generator and PV system.
Table 5. Normalized Curve of generation: Wind power generator and PV system.
Time (h) P W T (pu) P P V (pu)Time (h) P W T (pu) P P V (pu)Time (h) P W T (pu) P P V (pu)
0.50.630308.50.82710.040316.50.98920.4193
1.00.619409.00.85230.134417.00.96520.2784
1.50.609809.50.87880.271017.50.92440.1373
2.00.6050010.00.90640.367318.00.86070.0374
2.50.6122010.50.93280.458418.50.77430.0007
3.00.6411011.00.95200.612519.00.72510
3.50.6927011.50.96400.813419.50.71670
4.00.7395012.00.97000.912220.00.71670
4.50.7779012.50.97480.963320.50.72510
5.00.7887013.00.97841.000021.00.72630
5.50.7671013.50.98320.958221.50.71790
6.00.7479014.00.98800.879122.00.70950
6.50.7287014.50.99400.730822.50.69870
7.00.7371015.00.99880.764523.00.69150
7.50.7731015.51.00000.686623.50.68670
8.00.80310.001616.00.99640.589324.00.68310
Table 6. Optimal selection and location of the BESSs with the exact MINLP and the proposed MIQC models.
Table 6. Optimal selection and location of the BESSs with the exact MINLP and the proposed MIQC models.
Model TypeLocation and Type of the BESSLosses Cost (COP$)Error MINLP (%)
MINLP 13 ( A ) , 20 ( B ) , 21 ( B ) 47209.950.00
MIQC 5 ( A ) , 16 ( B ) , 21 ( B ) 41627.343.49
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Serra, F.M.; Montoya, O.D.; Alvarado-Barrios, L.; Álvarez-Arroyo, C.; Chamorro, H.R. On the Optimal Selection and Integration of Batteries in DC Grids through a Mixed-Integer Quadratic Convex Formulation. Electronics 2021, 10, 2339. https://doi.org/10.3390/electronics10192339

AMA Style

Serra FM, Montoya OD, Alvarado-Barrios L, Álvarez-Arroyo C, Chamorro HR. On the Optimal Selection and Integration of Batteries in DC Grids through a Mixed-Integer Quadratic Convex Formulation. Electronics. 2021; 10(19):2339. https://doi.org/10.3390/electronics10192339

Chicago/Turabian Style

Serra, Federico Martin, Oscar Danilo Montoya, Lázaro Alvarado-Barrios, Cesar Álvarez-Arroyo, and Harold R. Chamorro. 2021. "On the Optimal Selection and Integration of Batteries in DC Grids through a Mixed-Integer Quadratic Convex Formulation" Electronics 10, no. 19: 2339. https://doi.org/10.3390/electronics10192339

APA Style

Serra, F. M., Montoya, O. D., Alvarado-Barrios, L., Álvarez-Arroyo, C., & Chamorro, H. R. (2021). On the Optimal Selection and Integration of Batteries in DC Grids through a Mixed-Integer Quadratic Convex Formulation. Electronics, 10(19), 2339. https://doi.org/10.3390/electronics10192339

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