Scheduling of a Microgrid with High Penetration of Electric Vehicles Considering Congestion and Operations Costs

This paper reviews the impact that can be presented by the immersion of generation sources and electric vehicles into the distribution network, with a technical, operational and commercial approach, given by the energy transactions between customer and operator. This requires a mathematical arrangement to identify the balance between congestion and the operating cost of a microgrid when the operation scheduling of the system a day ahead of horizon time it is required. Thus, this research is directed to the solution, using heuristic algorithms, since they allow the non-convex constraints of the proposed mathematical problem. The optimization algorithm proposed for the analysis is given by the Multi-Object Particle Swarm Optimization (MOPSO) method, which provides a set of solutions that are known as Optimal Pareto. This algorithm is presented in an IEEE 141-bus system, which consists of a radial distribution network that considers 141 buses used by Matpower; this system was modified and included a series of renewable generation injections, systems that coordinate electric vehicles and battery storage, and the slack node was maintained and assumed to have (traditional generation). In the end it can be shown that the algorithm can provide solutions for network operation planning, test system robustness and verify some contingencies comparatively, always optimizing the balance between congestion and cost.


Introduction
In recent years, energy demand has increased significantly at the industrial and residential levels; this has led network operators (NOs) to increase their investments in terms of infrastructure and firm power. This means that the infrastructure should grow according to the needs of the demand. This is restrictive because the electric distribution network was designed based on a projected maximum demand and focused on a radial unidirectional flow, i.e., the energy supplied comes from a source to a load. In summary, this unidirectionality of demand leads the NO to frequently carry out studies to plan the network in order to ensure the continuous supply of energy.
At present, the low costs of the elements used in alternative energies, such as batteries, panels, and others, will make the unidirectionality of demand in the network unnecessary. This change is and will be caused by the installation of small energy sources by the customers in order to meet their energy needs. Given the variations in energy demand over time, there will be times when customers will have surplus energy and it will be profitable for them to sell it to the NO; consequently, energy supply will no longer depend on a single source but on several sources, either in close proximity or from the same customer or load. This immersion of more energy sources in the network is currently known as distributed generation (DG), which permits efficiently managing investments for the NO, meeting energy needs, and also a more active participation of customers by means of incentives at the regulatory, technical, and commercial level. On the technical side, it is important to mention that network connection generates some restrictions due to the Vehicles 2021, 3 579 technology (wind, solar, storage, etc.), impacts on infrastructure, electrical protections, and network congestion.
The immersion of these generation sources poses various challenges to the NO due to the change in the directionality of current flow, which should consider network congestion caused by voltage variations, overloads, increased reactive power generation, and problems in the delivery of efficient active power [1]. The immersion also manifests itself in the overload of lines and network transformers; this causes high current flows and heating in the distribution network conductors, the impact of which will result in equipment failures or sudden increased load or demand and in an undesired disconnection of service [2].
For congestion in the network, NOs are implementing technologies to better predict energy flow, and there is also more flexible equipment according to the needs of demand [3]. Congestion in a power system occurs because generation is trapped in the network, causing unscheduled load shedding and cascading failures. This is because networks are sometimes obliged to operate above their limits, due to the different restrictions that occur in the dispatch of energy [4].
In real-time operation, some practical methods have also been implemented, such as generation rescheduling, load shedding, load elimination, flexible alternating current transmission system (FACTS), and connection of storage systems. These methods involve large additional costs, as incentives should be created among generation and distribution companies to modify their preconceived schedules.
In order to optimize the operation and planning of distribution networks, it is necessary to carry out theoretical and forecast analyses aimed at making the best decisions, as well as to meet the real time operation, so that the network can be reliable and harmonized [2,5,6]. Approaches should lead the NO and the customer to optimal dispatch so as to minimize network congestion. For example, the different strategies for network congestion control defined by the actors include optimal power flow, nodal pricing, and structuring of contracts between the NO and the customer [2].
Optimal dispatch resulting from the analysis of several theoretical approaches includes genetic algorithms, numerical simulations, and heuristic optimization, among others. For any approach it is important to mention that there is a common and natural constraint: most radial-topology distribution networks have a high R/X ratio that affects the analysis of power flows due to their nonlinearity in the equations.
The results obtained from each formulation proposed lead to a more efficient way to decongest the network, to boost the electricity market, to increase the reliability of the network, to make investments, and other advantages.
To address the issue of network congestion is of great importance, since it is an essential factor to make the operation more effective and to measure the impact on the operation with the implementation of these new generation agents.
The proposal of this article aims to provide a clearer vision to NOs regarding the immersion of renewable sources and their effect on the network when congestion and operation costs are considered. The terms for this proposal arise from the current energy matrices of countries and the new changes in the operation of the NOs.

Problem Formulation
The main objective of this study is to propose a methodology for programming the operation of a test microgrid composed of renewable energies, energy storage, and electric vehicles. This study includes the minimization of congestion and operation costs by using heuristic optimization and by taking into account the own restrictions of the system.

Objective 1: Minimize Congestion
Based on the references presented in the introduction, the mathematical formulation that can represent this index is as follows: where Si is the MVA flow on line i in hour j Si max is the MVA capacity of line i in hour j The mathematical formulation described in Equation (1), and used in [7,10,14,16], was selected to carry out this study because it is not very complex, but is efficient and fast to implement and provides reliable results that will be discussed later. In addition, its congestion results are similar to those of the indicators analyzed in other studies. This is the first objective function selected to determine the level of congestion in the proposed problem.

Operation Costs of a Microgrid (Objective Function 2 of the Problem to Be Solved)
This section is based on the formulation developed in the study in [25]. When modeling a microgrid, the associated costs defined by the equipment (or sources) interacting with the system should be considered, in addition to the technical aspects. This article seeks to optimize the costs associated with the integration of alternative sources, such as solar panels, batteries, and vehicles, among others. The cost associated with each type of source varies according to its technology, and this generates a cost function. As can be seen, cost optimization will result in a more efficient operation of the system.
The costs or cost function that define the operation of the microgrid (CO) will be structured by the following generations costs: traditional generation F 1 (x), through electric mobility F 2 (x), battery storage F 3 (x), and uncertainty F 4 (x); these are represented in (2).
The inclusion of energy generation in the microgrid through fossil sources, such as diesel or gasoline, is a backup instrument when there are energy deficiencies in the network or demand from other sources such as the traditional one or renewable sources. To define the cost associated with this type of source, the generalized function in [26][27][28] is used, shown in Equation (3).

Cost of the Electric Vehicles F 2 (x)
In this section, plug-in electric vehicles (PEVs) are considered as part of the microgrid and have a parking aggregator (parking decks). Based on this consideration, vehicles are grouped so the dispatch occurs at the best time, defined by a cost-benefit analysis [29][30][31]. The main objective is to plan the charge of vehicles at the time of the day when the cost is lowest or optimal. If the PEV has not consumed all the energy stored, this excess can be returned to the batteries in the highest cost time slot [29].
The equation describing the planning, loading, operation, and transfer is based on [29,30,32]. The first step is defined through the aggregators when they calculate the maximum energy dispatch for each vehicle (PEV); the batteries are assumed to have charged all day. The start of charge is defined as SOC D,a ni and represented in Equation (4): PEVs have both upper and lower limits defined as t k , t k+H i −1 and the Equations (5) and (6). This group of equations helps define the expressions that denote the limits and the minimum energy that it can have in the period [k + j]. In the range, energy should be at least ρ·P max ·∆t and lower than the following interval [k + j + 1]. As it is a charging process, the energy obtained is incremental [29,32].
Vehicles (PEVs) have energy and power limits, determined by Equations (7) and (8). Naturally, the batteries have a limit in connection and load, described by Equations (9) and (10).
These equations define the energy and power limits for vehicles, in addition to Equations (11) and (12). Equation (13) defines the power limits of PEVs. , j = 0, . . . , Hi, ∀i ∈ I (12) The cost associated with the use of PEVs is determined through a function F 2 (x) defined in Equation (14). Equations (15)-(18) determine the restrictions. A cost equal to that generated by renewable generation equipment C bl is assumed in order to determine the storage cost (batteries) for this study [32][33][34][35][36][37]. Equation (19) determines this cost: where L loss refers to the useful life of the battery and C init−bat refers to the investment cost for the acquisition of batteries. The base case consists of the power being negative when the battery is charging and being positive when the battery is delivering power [32]. Due to the nature of batteries, there is a loss of useful life, represented as L loss . Equation (20) determines the cost associated with this loss; A c is the cumulative performance with each unit [Ah] effective in an interval of time and A total is the cumulative performance with [Ah] effective in the life cycle [36].
The variables that make up A c are the state of charge (SOC) and the actual performance Ah. Real Ac is represented in Equation (21): The parameter λ SOC is an effective weighted factor determined by Equation (22): The total cost will have the constraints established by the minimum and maximum state of charge of the batteries determined by (23). Equation (24) refers to the maximum power transfers at maximum discharge.

Cost of Operation Photovoltaic Generator and Wind Power Generator F 4 (x)
To establish the total cost of solar power generation (C PV ) in kW, two functions are used, as shown in Equation (25): the first function refers to the cost when the generator C PV,u,i is underestimated, and the second function, defined as C PV,o,i , refers to the opposite effect, i.e., when it is overestimated. The behavior of the two functions depends on both the power programmed in the photovoltaic power station (W PV,s,i ) and the defined available power (W PV,i ).
The costs underestimating the generator refer to the value that the microgrid does not receive when those kW are not sold as a traditional energy source [32,38]. The availability of solar energy power generates a variation in the cost with a probabilistic behavior due to the nature of the source itself, i.e., due to the variation of radiation [25,32,38].
Overestimation occurs when the photovoltaic generator cannot meet the expected power value and has to incur an additional cost as it has to resort to another generator [32,38].
The power obtained through wind energy also has probabilistic variability and causes overestimation or underestimation as occurs with solar energy [26,32,38]. Based on these two concepts, the total cost for this source can be estimated, as shown in (26).

583
Wind sources have the same underestimated behavior mentioned in the section on solar source: the microgrid operator does not need to incur a cost thanks to the dispatch of energy from a wind generator, while the overestimated cost refers to the operator's need to resort to other generators to meet the agreed demand.
Finally, the associated costs in F 4 (x) are denoted as the uncertainty probability cost of renewable agents; this is expressed by (25).

Optimization Methodology
A multi-objective problem is described as a vector of decision variables that meets constraints and optimizes a vector function whose elements represent the objective functions. These functions form a mathematical description of performance criteria that are often in conflict with each other. Therefore, the term "optimize" involves finding a solution that gives the values of all acceptable objective functions to the decision maker [39]. This type of problem can also maximize or minimize n functions or perform a combination of both together with a set of constraints of either equality or inequality.
Based on the previous definition taken from [40], a multi-objective problem is considered to have a set of solutions and not a single answer. This set of solutions is verified through Pareto optimal, which will be defined below.
Generally, most problems in engineering today involve the evaluation of several conflicting objectives, and feasible and optimal solutions are necessary to solve the proposed objectives.
Heuristic, metaheuristic, and genetic algorithms are widely used to develop multiobjective problems by using a good convergence speed to obtain the answer of all the possible spaces for the solution.

Multi-Objective Particle Optimization Algorithm
The Multi-Objective Particle Optimization Algorithm (MOPSO) was implemented by J. Moore and R. Chapman. It is a multi-objective optimization algorithm that belongs to evolutionary algorithms, and known to be competitive and efficient thanks to its convergence speed.
Its principle is based on the behavior of birds and a multidimensional search headed in the position taken by each individual. Therefore, each individual will always be affected by the individual with the best behavior, i.e., local or global individuals [41].
Its main parameters include the population to be evaluated and the unique memory used by the individuals to be considered, which is called "repository" [42]. The aim is to obtain a global repository in which each analyzed particle will store their flight experiences in each evaluated cycle [42]. In this way, the algorithm seeks to select the best positions by using a grid.
The grid method is used to establish a coordinate system and also check how many individuals are stored in the repository; the objective functions to be considered are established as the maximum and minimum fitness values of these individuals [41].
Equations (28) and (29) determine the grid limits (maximum and minimum fitness values).
The limits help define the neighborhoods to be considered in the multi-objective problem, which are known as grid index in the development of the problem.
In MOPSO, the identification of the set of solutions obtained in Pareto optimal is relevant to choose the best global particle (gbest). In multi-objective problems, this term may be different and have drawbacks for the selection of this one, so there is not a single solution but a set of solutions that make the solution of the problem possible.
In this algorithm optimal solutions are obtained through a defined search space; in each iteration new positions are found through the following equations [41,42]: In MOPSO, it is very important to choose the gbest, i.e., the best global particle from the set of non-dominated solutions obtained in Pareto optimal [43].
According to the above, in most cases multi-objective problems are addressed as a single-objective problem, considering the evaluation of solutions one by one (or objective by objective). When using these methods, it is evident that different simulations should be considered in order to solve certain problems. Consequently, addressing these multiobjective problems through genetic algorithms helps consider diverse solutions that are generated in Pareto optimal in a single algorithm: in most cases one solution improves one objective and in parallel worsens the next one.

Results
This study was based on the selection of a test system chosen for the validation of the proposed congestion method, which corresponds to a radial distribution system. The case consisting of 141 nodes of Khodr, Olsina, De Jesus, and Yusta from the IEE was selected, corresponding to a portion of the system of the metropolitan area of Caracas (Venezuela); the data associated can be found in [44] (Tables 1-4). This radial network comprises 141 voltage nodes, 84 load nodes, 1 generator, 3 controlled voltage nodes, and 140 distribution elements (lines).   The base case of the IEEE 141 system was chosen. Some of its parameters were modified, such as the capacities of the lines. Two new generators were also considered, which were located in nodes 109 and 110 with powers of 50 and 10 MW, respectively; these nodes are supposed to have more solar and wind resources. These values were subsequently modified by joining them in the functions used in the algorithm. The application of distributed generators was also considered with the help of the MOPSO algorithm. Changes in some nodes were expected in the real and reactive power demands, given the different types of energies from electric vehicles and battery storage, to verify the proposed method of congestion.
There was generation only in some bars according to the modification made. This was modified to verify the impact of the inclusion of renewable agents on the chosen test network. Some characteristics of the case study are shown below.
Based on the study in [25], the base case was chosen to check the congestion index of the system so as to weaken the original network, which is originally very robust.
The MOPSO optimization algorithm was adjusted and applied to the formulation of the problem, whose parameters are shown below: With respect to the values obtained shown in Figure 1, the optimal solution can be discarded when the congestion value is 8.28%, since it means that the lines of the test system will be overloaded. This leads to its possible exclusion, in spite of its cheaper operation cost, which may cause instability, generate load shedding, and result in economic and social damage. Consequently, the set of optimal solutions ensuring the system's safe operation should range between 0.65% and 1.03% of congestion. This prevents any failure from leading the study system to a blackout, in order to avoid any damage. In this case, the operator will consider the scenarios with a minimal congestion index to ensure a robust system ( Figure 2):  Figure 2 shows that maintaining a low congestion level involves a significant increase in the operation cost. The study therefore aims to optimize the operation cost corresponding to the decrease of network congestion, as described below.    Consequently, the set of optimal solutions ensuring the system's safe operation should range between 0.65% and 1.03% of congestion. This prevents any failure from leading the study system to a blackout, in order to avoid any damage. In this case, the operator will consider the scenarios with a minimal congestion index to ensure a robust system ( Figure 2 Consequently, the set of optimal solutions ensuring the system's safe operation should range between 0.65% and 1.03% of congestion. This prevents any failure from leading the study system to a blackout, in order to avoid any damage. In this case, the operator will consider the scenarios with a minimal congestion index to ensure a robust system ( Figure 2):  Figure 2 shows that maintaining a low congestion level involves a significant increase in the operation cost. The study therefore aims to optimize the operation cost corresponding to the decrease of network congestion, as described below.   Figure 2 shows that maintaining a low congestion level involves a significant increase in the operation cost. The study therefore aims to optimize the operation cost corresponding to the decrease of network congestion, as described below.

Variation of the Decision Variables at the Pareto Optimal Points
Based on the results obtained, Figure 3 shows the dispatch required to reach the minimal operation cost with this type of resource, in which the storage system has three charging cycles. During the storage system's charging times through materials, the traditional electric generation dispatches energy to the load and allows the storage system to recharge (Figure 4).

Variation of the Decision Variables at the Pareto Optimal Points
Based on the results obtained, Figure 3 shows the dispatch required to reach the minimal operation cost with this type of resource, in which the storage system has three charging cycles. During the storage system's charging times through materials, the traditional electric generation dispatches energy to the load and allows the storage system to recharge (Figure 4). The dispatch of solar energy takes place only during daylight hours: this energy source is neither available before dawn nor after dusk. Consequently, Figure 4 shows that the energy dispatch is lower between 11:00 a.m. and 3:00 p.m. than the available energy; likewise, from 12:00 to 2:00 p.m., a portion of energy is stored. In doing so, it is possible to supply more energy than that available when necessary or to dispatch it in the evening when the operation cost is more profitable.

Variation of the Decision Variables at the Pareto Optimal Points
Based on the results obtained, Figure 3 shows the dispatch required to reach the minimal operation cost with this type of resource, in which the storage system has three charging cycles. During the storage system's charging times through materials, the traditional electric generation dispatches energy to the load and allows the storage system to recharge (Figure 4). The dispatch of solar energy takes place only during daylight hours: this energy source is neither available before dawn nor after dusk. Consequently, Figure 4 shows that the energy dispatch is lower between 11:00 a.m. and 3:00 p.m. than the available energy; likewise, from 12:00 to 2:00 p.m., a portion of energy is stored. In doing so, it is possible to supply more energy than that available when necessary or to dispatch it in the evening when the operation cost is more profitable.  The dispatch of solar energy takes place only during daylight hours: this energy source is neither available before dawn nor after dusk. Consequently, Figure 4 shows that the energy dispatch is lower between 11:00 a.m. and 3:00 p.m. than the available energy; likewise, from 12:00 to 2:00 p.m., a portion of energy is stored. In doing so, it is possible to supply more energy than that available when necessary or to dispatch it in the evening when the operation cost is more profitable.
Similarly, Figure 5 shows that the energy dispatch surpassed the energy available in an 8 h period, and in the remaining 16 h of the day that amount of energy remained stored to be dispatched at convenient times, based on the costs of underestimation and overestimation. The actual dispatch performed during the daytime comprises a significant increase regarding the wind resource, given the operation costs involved when this sort of energy is dispatched to the grid.
Vehicles 2021, 3, FOR PEER REVIEW 11 Similarly, Figure 5 shows that the energy dispatch surpassed the energy available in an 8 h period, and in the remaining 16 h of the day that amount of energy remained stored to be dispatched at convenient times, based on the costs of underestimation and overestimation. The actual dispatch performed during the daytime comprises a significant increase regarding the wind resource, given the operation costs involved when this sort of energy is dispatched to the grid. Finally, Figure 6 shows that the highest power occurs during early morning hours; aggregator 2 presents a high power peak at around midnight.  Finally, Figure 6 shows that the highest power occurs during early morning hours; aggregator 2 presents a high power peak at around midnight.
Vehicles 2021, 3, FOR PEER REVIEW 11 Similarly, Figure 5 shows that the energy dispatch surpassed the energy available in an 8 h period, and in the remaining 16 h of the day that amount of energy remained stored to be dispatched at convenient times, based on the costs of underestimation and overestimation. The actual dispatch performed during the daytime comprises a significant increase regarding the wind resource, given the operation costs involved when this sort of energy is dispatched to the grid. Finally, Figure 6 shows that the highest power occurs during early morning hours; aggregator 2 presents a high power peak at around midnight.

Sensibility Analysis of Line Capacity in the Modified IEEE Case No. 141
This section describes the simulations run for the IEEE case study No. 141, which used the MOPSO genetic algorithm to verify the variation of the maximum power limit in the lines of the analyzed system. Based on this, the system underwent a number of modifications intended to determine the impact of congestion and operation costs, based on the assumption that the power capacity limit might change within five years.
This limit was modified for values above and below 45 MW, since the system had a better performance in that range (Figure 7). Vehicles 2021, 3, FOR PEER REVIEW 12

Sensibility Analysis of Line Capacity in the Modified IEEE Case No. 141
This section describes the simulations run for the IEEE case study No. 141, which used the MOPSO genetic algorithm to verify the variation of the maximum power limit in the lines of the analyzed system. Based on this, the system underwent a number of modifications intended to determine the impact of congestion and operation costs, based on the assumption that the power capacity limit might change within five years.
This limit was modified for values above and below 45 MW, since the system had a better performance in that range (Figure 7). According to Figure 8, the congestion index corresponding to a minimal operation cost is not directly correlated to the system's power carrying capacity.  Figure 9 shows that the network does not become congested based on the operation cost structure depicted in Figure 10. The analysis was focused on carrying capacities of more than 50 MW.  According to Figure 8, the congestion index corresponding to a minimal operation cost is not directly correlated to the system's power carrying capacity.

Sensibility Analysis of Line Capacity in the Modified IEEE Case No. 141
This section describes the simulations run for the IEEE case study No. 141, which used the MOPSO genetic algorithm to verify the variation of the maximum power limit in the lines of the analyzed system. Based on this, the system underwent a number of modifications intended to determine the impact of congestion and operation costs, based on the assumption that the power capacity limit might change within five years.
This limit was modified for values above and below 45 MW, since the system had a better performance in that range (Figure 7). According to Figure 8, the congestion index corresponding to a minimal operation cost is not directly correlated to the system's power carrying capacity.  Figure 9 shows that the network does not become congested based on the operation cost structure depicted in Figure 10. The analysis was focused on carrying capacities of more than 50 MW.  Figure 9 shows that the network does not become congested based on the operation cost structure depicted in Figure 10. The analysis was focused on carrying capacities of more than 50 MW.          As demonstrated in Figure 11, however much the operation cost increases, it is not possible to decongest the network for the test system; a capacity lower than 45 MW is not optimal. As demonstrated in Figure 11, however much the operation cost increases, it is not possible to decongest the network for the test system; a capacity lower than 45 MW is not optimal.

Discussion and Conclusions
The multi-objective solution made it possible to evaluate the congestion resulting from the immersion by means of different energy sources in a test microgrid. The study found a very significant relation between improvement and decrease in grid operation, based on operation costs and the congestion index of the system subject to its own restrictions.
The MOPSO algorithm, being a population-based optimization model, does not converge at the same optimal point, but does provide a suboptimal solution quite close to that point. The execution of the optimization model provided a set of Pareto optimal solutions for any congestion index; this allowed the system operator to choose from several options to make better operational and commercial decisions.
Evidence shows that there are solutions that oppose each other, i.e., obtaining a minimal congestion index involves a higher operation cost. In other words, the improvement of one objective-or n objectives-leads to a decrease in the other ones. The load curve used presented most congestion problems in the test system considered.
A system with ample carrying capacity helps reduce congestion in the network. This study aimed at optimizing the congestion index of the microgrid analyzed, since it is an objective of the NOs; ensuring that the final user's demand will be met; providing continuity of service; and obtaining greater system reliability.
This study provides a tool to run simulations that may contribute to network planning and to ensure strategies for sufficient generation and carrying capacity. The equivalent model makes it possible to approach the system's initial state in a close and simplified way. The analysis of the incursion of generation projects soon to be installed will offer an idea of the system variables and their behavior to define the strategies to be chosen, such as the repowering of congested lines.
The evaluation of the decision variables at each Pareto point shows that, when dispatching the system, it is sometimes more favorable for energy to be stored when there is exceeding solar and wind energy. It was also found that diesel-based generation is always used because the dispatch involves a lower cost at given hours and the supply is necessary at more critical times.

Discussion and Conclusions
The multi-objective solution made it possible to evaluate the congestion resulting from the immersion by means of different energy sources in a test microgrid. The study found a very significant relation between improvement and decrease in grid operation, based on operation costs and the congestion index of the system subject to its own restrictions.
The MOPSO algorithm, being a population-based optimization model, does not converge at the same optimal point, but does provide a suboptimal solution quite close to that point. The execution of the optimization model provided a set of Pareto optimal solutions for any congestion index; this allowed the system operator to choose from several options to make better operational and commercial decisions.
Evidence shows that there are solutions that oppose each other, i.e., obtaining a minimal congestion index involves a higher operation cost. In other words, the improvement of one objective-or n objectives-leads to a decrease in the other ones. The load curve used presented most congestion problems in the test system considered.
A system with ample carrying capacity helps reduce congestion in the network. This study aimed at optimizing the congestion index of the microgrid analyzed, since it is an objective of the NOs; ensuring that the final user's demand will be met; providing continuity of service; and obtaining greater system reliability.
This study provides a tool to run simulations that may contribute to network planning and to ensure strategies for sufficient generation and carrying capacity. The equivalent model makes it possible to approach the system's initial state in a close and simplified way. The analysis of the incursion of generation projects soon to be installed will offer an idea of the system variables and their behavior to define the strategies to be chosen, such as the repowering of congested lines.
The evaluation of the decision variables at each Pareto point shows that, when dispatching the system, it is sometimes more favorable for energy to be stored when there is exceeding solar and wind energy. It was also found that diesel-based generation is always used because the dispatch involves a lower cost at given hours and the supply is necessary at more critical times.
The modeled system helps propose several congestion index optimization strategies, such as the optimal location of renewable energy injection or controllable loads taking into account the generation capacity necessary to obtain a minimal congestion index and a lower operation cost. The verification of possible implications on the congestion index and on the decision variables evaluated can be performed through a sensibility analysis of load curves in the studied system. This study is a tool that can be modeled in the distribution systems of NOs to optimize congestion and potential costs derived from a widespread distributed generation.