1. Introduction
Rapid advances in power electronics, i.e., conversion technologies, have transformed the classical distribution networks, where power comes from large scale generators interconnected through transmission grids to the substations, into modern electrical networks where interactions between generators and end-consumers have a place in the same voltage level (e.g., medium- or low-voltage distribution) [
1,
2,
3]. These modern networks have been named as active distribution grids or smart grids, with these characteristics: (i) the large-scale integration of power electronic converters to interface renewable energy resources and battery energy storage systems [
4], (ii) the possibility to be designed in AC, DC or hybrid configurations [
5], (iii) the ability to work in microgrids that can be connected or disconnected from the main grid to manage its energy resources efficiently [
6].
In the case of new distribution networks or extensions of the existing ones, the DC paradigm is an attractive alternative since these grids are easily controllable due to the simplification of the number of variables, i.e., the objective of control in these grids is the voltage profile provided in all the nodes of the networks [
5,
7,
8]. This is because reactive power and frequency are concepts only applicable to AC distribution. In addition, DC networks present lower power losses and best voltage profiles (no inductive effect in steady-state conditions) since only active power is transported by the feeders, which is a remarkable advantage compared to classical AC networks [
9].
The DC networks can be designed and operated using multiple renewable energy resources, such as photovoltaic panels [
10], energy storage devices (batteries, supercapacitors or superconductors) [
11], home appliances [
12] and computer servers [
13], among others. Observe that in the case of applications where AC energy is required (i.e., motors), this energy can be provided by interfacing the DC distribution network with voltage source inverters that generate sinusoidal voltages and the required reactive power at the load point [
4]. The main advantage of using converters to interface AC loads is making them controllable, i.e., they can be added into the grid operation environment as active consumers that may or may not be dispatched inside of demand response plans [
14]. Nevertheless, the implementation of DC networks for power distribution has some difficulties, such as high costs in the application of protection schemes [
15,
16], power electronic devices such as voltage source converters required to generate the voltage output [
4] and high costs related to measuring devices [
17], among others. Notice that these aspects must be taken into account when DC technologies are preferred over classical AC technology, for the electrical networks expansion in urban and rural areas in order to guarantee the return of investment in the planning period [
17].
The DC paradigm distribution is definitely an option in modern power systems based on the advantages previously mentioned. In specialised literature during the last two decades, the DC networks have grabbed much attention, from large-scale to small-scale power applications (i.e., HVDC transmission networks and LVDC distribution applications) [
5,
18,
19]. To analyse these grids, new power flow methods and optimal power flows strategies have been recently proposed based on exact and/or metaheuristic optimisation for the analysis under steady-state conditions [
7,
20]. The main characteristic of these analyses is that the topology of the DC networks has been considered as a fixed input of the problem [
21]. Nevertheless, in real application cases, the DC grid configuration can be variable. Due to the high variation of the renewable energy resources, topology changes are needed to obtain the best performance of the network to efficiently meet the loads interconnected on it [
22]. Based on this classical consideration, when DC grids are analysed, this research is motivated in presenting an alternative of analysis for DC distribution feeders taking into account changes in the DC grid structure to improve its operative efficiency. For this purpose, a new mathematical model for optimal reconfiguration of DC distribution networks is presented, based on a mixed-integer nonlinear formulation that selects the best subset of distribution lines that must be connected to provide energy to the end-consumers with the lowest power losses.
The problem of the optimal reconfiguration of electrical distribution networks has been largely explored in AC electrical networks with multiple optimisation approaches, most of them, based on metaheuristics, such as genetic algorithms [
23], tabu search [
24], particle swarm optimisation [
25], ant colony [
26] or modified culture search algorithm [
27], and so on. Authors in [
28] have proposed a multiobjective optimisation approach to solve the problem of the optimal reconfiguration in distribution networks simultaneously considering the optimal location and sizing of renewable energy sources based on wind and photovoltaic plants. As a solution technique, the strength Pareto evolutionary algorithm II has been proposed, which allows the obtaining of the optimal Pareto from between annual operation costs and the number of greenhouse emissions. Fuzzy set theory is used to select the best compromise solution among acquired Pareto set. Numerical results demonstrate the applicability and effectiveness of the proposed approach, mainly in the electrical networks feed by fossil sources. In [
29], a metaheuristic approach named the salp swarm algorithm has been proposed to solve the problem of the optimal reconfiguration problem in distribution networks considering the optimal location and sizing of distribution generators simultaneously. Numerical results in two classical test feeders composed by 33 and 69 nodes demonstrate that the solution of this problem in the same stage is more efficient than the uncoupled approaches reported in the literature when the objective is to minimise the total grid losses and the voltage profile improvement. Authors of [
30] have presented an optimal reconfiguration approach for distribution networks considering realistic time-varying load models and renewable energy resources in a simultaneous manner. As a solution technique, a multiobjective genetic algorithm is employed to minimise the total power losses and maximise the annual cost savings. Numerical results in the 33-node test feeder demonstrate efficiency and robustness compared to literature reports based on metaheuristics regarding Pareto front and processing times. In [
31], an optimal reconfiguration approach has been proposed for active distribution networks. It is composed of multiple microgrids considering a daily economic dispatch environment as a function of the amount of renewable energy and energy storage availability. The main objective of the mathematical model presented is the minimisation of the total energy losses in the grid during its operation. The solution of this problem was to address a hybrid combination between conventional and selective particle swarm optimisation methods. Numerical results in the 33-node test feeder demonstrate the effectiveness and robustness of the proposed approach regarding processing times via parallel computing. In reference [
32] has recently proposed an optimal reconfiguration approach in distribution networks in order to improve reliability indices and decrease the operational cost of power systems under a microgrid environment by considering electric vehicles and renewable energy. As a solution technique, an improved version has been proposed of the classical genetic algorithm, and all the numerical validations were made in the 33-node test feeder where the optimal reconfiguration plan optimises the grid operation (i.e., energy losses minimisation) when a shortage in renewables or electric vehicles are redispatched. Authors of [
33] have presented a complete literature review about optimal reconfiguration approaches by offering the most typically objective functions and solution techniques with the main focus on the solution technique based on heuristic and metaheuristic optimisation approaches.
In the case of DC networks, the first approach reported in the literature was made in 2020 in [
22] where the authors propose the optimal reconfiguration of low-voltage DC networks to maintain voltage profiles in an acceptable range. To represent the problem, these authors reformulate the MINLP model into an integer second-order cone programming formulation, which is solved through a branch and bound approach. It is worth mentioning that the work reported in [
22] increases in a quadratic manner the number of variables required to represent the equivalent optimisation problem. This characteristic can increase the processing times with the additional need to use specialised software to achieve the optimal solution [
34]. The combination of convex optimisation models, i.e., second-order cone programming to solve the conventional power flow, with branch and bound methods to deal with binary variables, can guarantee the optimal global solution of the relaxed optimisation model [
35]. Nevertheless, in the final solution reported by the mixed-integer convex approximation, small estimation errors can be presented when these are evaluated in the exact MINLP model, mainly caused by the approximations employed in the optimal power flow problem [
36,
37]. Note that the main difference between our proposal and the work presented in [
22] are
(i) the use of power losses as an objective function by guaranteeing that all the voltage profiles are within allowable ranges and
(ii) the use of our approach to reconfigure the DC network or can even be used to make technical studies regarding the optimal expansion of the DC network to provide electrical service to new potential users. The complete list of contributions of our approach is listed below.
The main contributions of this research regarding DC distribution can be summarised as follows:
- ✓
The proposal of a new mathematical model that represents the problem of the optimal reconfiguration of DC distribution networks with an MINLP structure, taking into account the power balance equations as a function of the voltage profiles and the branch currents through an incidence branch-to-node incidence matrix.
- ✓
The GAMS package presentation as a solution tool in a tutorial form that introduces bachelor engineering students and researchers in the modelling of mathematical optimisation problems. This is specially important because GAMS interface allows the implementation of MINLP models with the same compact structure of the symbolic model that represents the problem under study, which helps to learn about it (GAMS) in an intuitive way applying to all the areas related with mathematical optimisation.
- ✓
The prospect to extend the proposed MINLP model to the optimal selection of routes in DC distribution networks in order to analyse future extensions of the existing DC networks or make projections for future distribution networks.
This research is framed within the active distribution network analysis, specifically in the DC distribution area, which presents a new MINLP model that represents the optimal reconfiguration problem in DC grids. In this sense, we concentrate our efforts on the representation of the technical operative aspects of the network, i.e., operational variables such as currents in lines, voltages at nodes and power losses. Thus, the proposed objective function is to minimise the total power losses caused by the distribution activity. Nevertheless, economic aspects, such as construction and maintenance of lines and cost of acquisition of power electronic interfaces are considered for future research. In addition, possible bipolar configurations of electrical DC distribution networks are not considered in this research. Nevertheless, as presented in [
38], it will be possible to reformulate the proposed MINLP model to introduce these possible grid configurations that can originate fascinating and attractive future works to the scientific community.
The remainder of this document is structured as follows:
Section 2 presents the mathematical formulation of the optimal reconfiguration problem in DC networks by highlighting the importance of the branch-to-node incidence matrix in this formulation.
Section 3 presents the solution strategy based on the implementation of the MINLP model in the GAMS interface, using a tutorial sense with a small DC test feeder composed of 6 nodes and 10 candidate lines (this example offers the optimal selection of routes in DC distribution networks).
Section 4 describes the main features of the test systems used to validate the proposed MINLP model, which are composed of 10 and 21 nodes. These are classically used in power flow and optimal power flow analysis for DC networks.
Section 5 presents the numerical results and their discussion by highlighting the main advantages of DC grids reconfiguration for power loss reduction and voltage profile improvements. Finally,
Section 6 presents the main conclusions derived from this research and some possible future investigation lines.
2. Mathematical Formulation
The problem of the optimal reconfiguration in direct current networks can be represented as a mixed-integer nonlinear programming model based on the following facts: (i) the presence of binary variables regarding the selection of the branch l, i.e., , and (ii) the nonlinear relation between voltage and currents related to the power balance at each node of the network when constant power consumption is presented.
Remark 1. The power balance equations in a dc network are typically formulated with nodal voltage variables when the interest is to solve optimal power flow problems or optimal location of shunt devices. Nevertheless, in the case of the optimal reconfiguration problem, the best representation of the problem is based on the currents flowing the branches. These currents allow to calculate the effect of the reconfiguration on the DC network, directly by modifying the resistive parameters in the distribution lines [39]. To develop the mathematical formulation of the reconfiguration problem in DC networks, let us define the branch-to-node incidence matrix as follows:
Definition 1 (Branch-to-node incidence matrix).
The branch-to-node incidence matrix is a rectangular matrix composed of n columns associated with the nodes of the network and b rows related to the all branches of the network [40,41]. This matrix is filled by based on the criteria presented in Equation (1)where is the unitary element that represents the line connection between vertices and on the graph that represents the electrical connection among nodes. The branch-to-node incidence matrix presented in Definition 1 relates nodal injected currents to branch currents, allowing the transformation of the classical power balance equations formulated with voltages into a formulation that combines voltages and branch currents.
The complete mathematical model representing the problem of the optimal reconfiguration of DC distribution feeders is described as follows [
39].
Objective function:
where
corresponds to the objective function value associated to the total power losses in all branches of the DC system,
is the resistive effect related to the branch
l (The term resistive effect makes a reference to the physical opposition to the current flow that is intrinsically associated with real conductors due to their atomic composition, i.e., corresponds to the resistance value of the conductor at ambient temperature), and
is the current that flows through this line.
Set of constraints:
where
and
represent the power generation and power consumption at node
k,
is the voltage value in the node
k,
is the component of the branch-to-node matrix that defines whether the node
k and line
l are connected,
is the binary variable associated with the selection of the branch
l in the reconfiguration plan,
is the conductance value associated with the branch
l, i.e.,
;
n represents the number of nodes and
corresponds to the number of voltage controlled sources (i.e., slack nodes);
and
are the minimum and maximum voltage values permitted at the node
k;
is the maximum current that is supported by the conductor associated to the branch
l. Note that
and
are the sets that contain all the branches and nodes, respectively.
The interpretation of the mathematical model (
2)–(
7) is as follows: Equation (
2) represents the objective function related to the minimisation of the power losses in all the conductors of the DC network. Expression (
3) defines the power balance equation at each node, i.e., the application of the first Kirchhoff’s law and the first Tellegen’s theorem. Equation (4) determines the magnitude of the current that flows through the branch
l, which is a function of the binary variable
. Expression (5) determines the number of conductors that can be activated in the final configuration of the problem, added to (
3) which is widely-known in specialised literature as the radial constraint [
39]. Equation (6) is a type-box constraint that determines the voltage regulation bounds allowed in the DC network. Finally, Expression (7) determines the thermal capabilities of the conductor connected to the branch
l.
Remark 2. The solution of the optimisation model (2)–(7) is not an easy task, since it corresponds to a MINLP formulation, which is indeed an NP-hard optimisation problem [39]. For this reason, we select the general algebraic modelling (GAMS) as a solution methodology due to its excellent results reported in literature [42]. It is important to mention that in modern power systems, the most used devices to implement reconfiguration plans are the reclosers, since they can be operated in remote manner trough the SCADA system [
43,
44]. This implies that real-time reconfiguration plans can be programmed as a function of the grid and load conditions [
22].
3. Solution Strategy
We employ the GAMS optimisation environment to solve the mathematical model associated with the optimal reconfiguration of DC distribution network. GAMS is focused on the mathematical model structure and not on the solution methodology. This feature makes GAMS an ideal alternative to introduce bachelor engineering students and researchers into the mathematical optimisation world.
We select the GAMS optimisation package as a solver since it has been employed mainly in mathematical optimisation, as reported in the scientific literature. Some of the problems solved with GAMS are optimal operation and allocation of batteries in AC and DC power systems [
45,
46], stability analysis in DC networks [
47], the optimal location of distributed generators in AC and DC networks [
8,
42], optimal design of water distribution networks [
48], optimal design of osmotic power plants [
49], optimal design of thermoacoustic engines [
50], optimal economic dispatch in thermal power systems [
51], and optimal design of a standing wave thermoacoustic refrigerator [
52], among others.
The main advantages of using GAMS software to solve optimisation problems are:
(i) the usage of sets to define subindexes of variables to make a compact formulation of the mathematical problem under study.
(ii) the organisation of the numerical information, i.e., constant parameters, vectors, and matrices, in tables as a function of the sets’ indexes.
(iii) the definition of all the variables of the problem, i.e., continuous, binary and/or, discrete variables, and their maximum and minimum bounds.
(iv) the definition of all the equations representing the mathematical structure of the problem in compact form, i.e., as a function of the sets previously defined.
(v) the selection of the solution approach, i.e., maximisation or minimisation and the nature of the problem (e.g., linear programming, nonlinear programming, quadratic programming, mixed-integer nonlinear programming, and so on).
Appendix A illustrates an example of the implementation of mixed-integer nonlinear programming model that represents the problem of the optimal reconfiguration of DC grids [
51]. In addition, in the flow chart presented in
Figure 1, the main steps using the GAMS environment is presented. Note that this implementation is based on the general description of GAMS software reported in [
46].
6. Conclusions and Future Works
A mixed-integer nonlinear programming model for optimal reconfiguration of DC feeders has been proposed in this research. The mathematical formulation was founded in the power balance equations rewritten as a function of the branch currents using the branch-to-node incidence matrix. Numerical results in five test feeders confirm that is possible to reach optimal reconfiguration plans by implementing the proposed MINLP model in the GAMS software with processing times lower than 30 s in all the cases.
A numerical validation in a medium-voltage test feeder composed of 23 nodes demonstrates that the proposed MINLP model can be adapted to reconfiguration and expansion plans, since its mathematical models are identical, with the main advantage that multiple one or multiple slack nodes can use it.
Simulations in the 33- and 69-node test feeders demonstrate that the proposed MINLP model is suitable for being applied to large scale medium-voltage DC grids with power losses minimisation about and , respectively. In addition, numerical validations in the 10-node test feeder show that the MINLP formulation combined to the GAMS software is suitable for implementing real-time reconfiguration plans with minimum computational efforts as a function of the power consumptions at the load points.
Some of the possible future works derived from this paper can be the following:
(i) extend the proposed MINLP model to the optimal selection of wires gauge for DC grids considering technical and economic aspects, as well as to propose a mixed-integer convex relaxation to address the problem studied in this contribution with optimality guarantee;
(ii) the inclusion of temporal variable in the MINLP formulation to address the issue of optimal reconfiguration of DC networks considering renewables and battery energy storage devices;
(iii) the implementation of metaheuristic optimisation techniques to solve the mathematical model (
2)–(7) in a sequential manner (master-slave optimisation) in order to develop free access software for students, researchers or even distribution companies; and
(iv) study the possibility of extending the proposed MINLP model to electrical networks that combines AC and DC feeders to offer optimal reconfiguration plans in modern hybrid distribution systems.
As a recommendation, we suggest free access software with the following tools: PowerModels (transmission, ac), PowerModelsDistribution (distribution, ac) and PowerModelsACDC (transmission, ac+dc+converters). These have been developed in Julia environment and can be used to implement power system optimisation problems, as an alternative to GAMS to solve the proposed MINLP model.