Next Article in Journal
Towards Circular Economy—A Comparative Analysis of the Countries of the European Union
Previous Article in Journal
Safeguarding of Key Minerals Deposits as a Basis of Sustainable Development of Polish Economy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamic Reactive Power Compensation in Power Systems through the Optimal Siting and Sizing of Photovoltaic Sources

1
Ingeniería Eléctrica, Universidad Distrital Francisco José de Caldas, Bogotá 11021, Colombia
2
Facultad de Ingeniería, Universidad Distrital Francisco José de Caldas, Bogotá 11021, Colombia
3
Laboratorio Inteligente de Energía, Universidad Tecnológica de Bolívar, Cartagena 131001, Colombia
4
Facultad de Ingeniería, Institución Universitaria Pascual Bravo, Campus Robledo, Medellín 050036, Colombia
*
Author to whom correspondence should be addressed.
Submission received: 12 April 2021 / Revised: 28 April 2021 / Accepted: 4 May 2021 / Published: 11 May 2021

Abstract

:
The problem of the optimal placement and sizing of photovoltaic power plants in electrical power systems from high- to medium-voltage levels is addressed in this research from the point of view of the exact mathematical optimization. To represent this problem, a mixed-integer nonlinear programming model considering the daily demand and solar radiation curves was developed. The main advantage of the proposed optimization model corresponds to the usage of the reactive power capabilities of the power electronic converter that interfaces the photovoltaic sources with the power systems, which can work with lagging or leading power factors. To model the dynamic reactive power compensation, the η -coefficient was used as a function of the nominal apparent power converter transference rate. The General Algebraic Modeling System software with the BONMIN optimization package was used as a computational tool to solve the proposed optimization model. Two simulation cases composed of 14 and 27 nodes in transmission and distribution levels were considered to validate the proposed optimization model, taking into account the possibility of installing from one to four photovoltaic sources in each system. The results show that energy losses are reduced between 13% and 56% as photovoltaic generators are added with direct effects on the voltage profile improvement.

1. Introduction

Power systems are responsible for interconnecting power generation sources and consumers in high-voltage levels through transmission and sub-transmission systems since large power sources (i.e., hydraulic power plants) and loads are geographically separated by hundreds of kilometers [1,2]. The sector of electricity supply is composed of four main activities: (i) generation, (ii) transportation, (iii) distribution, and (iv) commercialization; these activities make possible the power supply from large and distributed power sources to all end-users from high- to low-voltage levels [3]. The main difference between the transmission and distribution system is associated with the grid topology since the former has multiple meshes to ensure reliability and security. At the same time, the latter is built typically with a radial structure to minimize investment costs and reduce the complexity in the coordination of the protective devices [4]. The electrical power system, in general, can be analyzed from the dynamical and static point of view. In the case of dynamic analysis, differential equations are used to determine the time-domain behavior of the electrical variables (i.e., frequency, active and reactive power flows, etc.) in the event of a large disturbance such as a short-circuit or a large load/generator disconnection [5]. While the static point of view is employed to determine the state variables of the network for a particular load and generation condition considering steady-state operative conditions [6].
In the static analysis of the power systems, one of the main preoccupations in the transmission and distribution sectors corresponds to energy losses caused during energy transportation. The Colombian system is not alien to this problem, since, in the case of high-voltage levels (i.e., ≥220 kV), the number of power losses oscillates between 1.5 % to 1.8 % of the total power generated; while in the case of the distribution grids, these percentages are between 6 % to 18 % [7]. To reduce the amount of energy wasted in electricity transmission and distribution activities, one of the main approaches studied by regulatory entities and academics corresponds to the inclusion of the renewable generation since it helps to: (i) reduce the harmful effects of greenhouse gas emissions in predominately thermal generation systems [8], (ii) improve the grid performance, i.e., voltage support and power loss minimization; and (iii) redistribute the power flow in lines, i.e., reduction of the power system stress, which helps make the system more robust regarding contingencies [9,10].
Different models and methodologies which have been developed to include renewable generation appropriately in conventional transmission and distribution networks are presented below.
The problem of the optimal placement and sizing of wind power sources in power systems considering variable power factor have recently been studied by Gil-González, et al. in [9], where the problem was formulated using a mixed-integer nonlinear programming (MINLP) model. The converter interfaces of the wind energy conversion system are used to provide reactive power to the grid when necessary via a nonlinear controller design [11]. The solution of the proposed model has been solved in the General Algebraic Modeling System (i.e., GAMS) software with excellent performance for radial distribution grids and meshed power systems.
For example, the authors of [12,13] have proposed differential evolution algorithms to address the problem of the optimal location, sizing, and power factor of multiple distributed generators (DGs) in order to minimize the power losses in a network. The bee swarm optimization method was employed for the sizing of a hybrid system composed of photovoltaic (PV) panels, wind turbines (WT), and fuel cells in [14]. In references [15,16], a genetic algorithm was implemented for optimal placement and sizing of DGs. This algorithm was also used by Ganguly and Samajpati in [17] to manage DGs’ optimal management when integrated into the network. In the specialized literature, hybrid algorithms have also been used, such as a mixed bee swarm and genetic algorithm for DG allocation as proposed by Mahmoud et al., in [18], in order to improve the voltage regulation and to reduce the active and reactive power losses in the distribution network. This method also reduced the number of iterations and improved the standard deviation values in all the case studies used when compared to other proposed algorithms in the specialized literature. Mathematical models of algebraic nature have also been proposed and solved in the GAMS software [19]. Moreover, it is possible to have adequate planning regarding the optimal location and dimensioning of DGs, as reported in [9]. An optimization methodology to size and locate the PV, WT, and capacitor banks considering the uncertainties of the demand caused by electric vehicles and the variations of irradiation, wind speed, and the usual load was presented in [20], where it was applied to a typical radial network. The author of [21] proposed a virtual power plant model that integrates the renewable energies available in the study area to maximize the profit of the operation by satisfying the given demand. The possible combinations of integrating DGs in a network involve multiple search engines based on hierarchical clustering, as proposed in [22], to reduce computation times.
The injection of reactive power to the grid is a significant study area since it allows correcting the power factor and avoiding voltage drops. For example, Quezada et al., in [23] explores the optimal location of capacitive banks through heuristic methods to improve the aspects mentioned above; using batteries, the same objective can also be achieved as observed in references [24,25]. Current technology allows the integration of photovoltaic generators via a three-phase inverter, which, in addition to complying with the basic requirements for their integration, can inject or absorb reactive power into the distribution system, as developed in [26] and improved by [27] using algorithms capable of searching for the maximum power point tracking.
In this paper, a model for the optimal location and sizing of PV plants in AC power systems considering reactive power injection was formulated. The proposed model generates the mixed-integer nonlinear programming (MINLP) model, which is solved using the GAMS software as a tool, as used in similar works [9]. Under the presented model, the PV can inject or absorb reactive power by employing the power electronics devices’ reactive power management supply capability for integration with the power grid, according to the grid requirements. Additionally, a loadability factor was integrated, which is used in the power electronic interface to limit the PV’s active power to the grid. Note that the main differences of this study regarding the reference [9] are: (i) the usage of real generation and demand curves in the Colombian power system context to characterize the demand and solar radiation behaviors; (ii) the inclusion of the sensitivity analysis in the optimization model regarding the numerical performance of the voltage profiles as a function of the number of PV plants integrated into the electrical grid; and (iii) the usage of the reactive power capabilities of the power electronic interface in photovoltaic applications which have not been previously explored in the relevant literature and was identified as an opportunity of research that this paper attempted to contribute to from the point of view of the exact mathematical optimization.
The remainder of this paper is structured as follows: Section 2 describes the mathematical optimization associated with the siting and sizing of PV sources considering the possibility of injecting reactive power; Section 3 presents the main aspects associated with the implementation of an optimization model in the GAMS optimization package; Section 4 shows the main characteristics of the test feeders composed of 14 and 27 nodes, respectively. Section 5 presents all the numerical simulations and their analysis and discussion; Section 6 presents the main concluding remarks derived from this work and some possible future developments.

2. Mathematical Formulation

The problem of the optimal location and sizing of PV systems in electrical grids is an MINLP problem [9], where the binary variables are associated with the location of the PV sources and the integer variables correspond to the power flow variables, i.e., voltages, currents, and powers, among others [9]. The complete mathematical formulation of this specific problem is presented below.

2.1. Objective Function

The objective function for the problem of the optimal placement and sizing of PV systems in electrical grids considered in this research corresponds to the amount of energy losses during a day (24 load/generation periods) of operation of the network. Equation (1) shows the objective function of the mathematical model.
min z = t Ω T i Ω N V i , t j Ω N V j , t Y i j cos ( θ i , t θ j , t ϕ i j ) Δ T ,
where z is the objective function value regarding with the daily energy losses in a power system, V i , t and V j , t are the magnitudes of the voltages at nodes i and j, respectively, Y i , j is the magnitude of the admittance between node i and j obtained from the matrix of admittances representing the power system, and θ i , t , θ j , t and ϕ i j are the angles of node i voltage, node j voltage and admittance between both, respectively. Δ T is the length of the period of time where the generations and loads are constant (typically one hour). Note that Ω N is the set that contains all the nodes of the grid and Ω T is the set that contains all the periods of time of the operation horizon.

2.2. Set of Constraints

The set of constraints associated with the problem of the optimal placement and sizing of PV systems in electrical grids includes the active and reactive power balance equations, the voltage regulation bounds, as well as PV generation bounds and the maximum number of PVs available for installation, among others. All the constraints considered in this study are listed below [28]
P i , t C G + y i P V P i , t P V n o m P i , t D = V i , t j Ω N V j , t Y i j cos ( θ i , t θ j , t ϕ i j ) , { i Ω N , t Ω T }
Q i , t C G + Q i , t P V Q i , t D = V i , t j Ω N V j , t Y i j sin ( θ i , t θ j , t ϕ i j ) , { i Ω N , t Ω T }
V i min V i , t V i max , { i Ω N , t Ω T }
0 y i P V P i P V max x i P V , { i Ω N }
y i P V η 1 η 2 ( P i , t P V n o m ) 2 Q i , t P V y i P V η 1 η 2 ( P i , t P V n o m ) 2 , { i Ω N , t Ω T }
i Ω N x i P V N G P V max ,
where Equations (2) and (3) define the balance of active and reactive power for each node of the system and for each period of time, P i , t C G and Q i , t C G being the active and reactive power delivered by the conventional generator located at node i; y i P V is the nominal power of the PV system; P i , t P V n o m is the active power curve of the PV system for node i at each period of time; Q i , t P V is the reactive capability of the PV system at node i; and P i , t D and Q i , t D are the active and reactive power demanded at node i in the period t. Equation (4) shows the voltage regulation constraint between a minimum ( V i min ) and a maximum ( V i max ) value allowed by regulatory entities. Equation (5) is the constraint on the power curve, where x i P V is a binary variable indicating whether a PV panel is located at node i; if x i P V = 1 , then the PV system is located at node i, and if x i P V = 0 , then it is not. Equation (6) shows the interval (area) where the reactive power delivered by the photovoltaic must be defined, and η is the chargeability factor of the PV system-inverter pair, and this value is between 0 and 1. Inequality constraint (7) shows the limit of generators to be placed in the power system, where the maximum quantity allowed is N G P V max . It observes that all the values of the electrical variables described here are per unit values.
Note that the derivation of the reactive power bounds presented by the nonlinear constraint (6) can be made using the power triangle applied to the power electronic converter that interface the PV system to the electrical grid [9]. The derivation of these active and reactive power bounds for the PV system is described below (the η -coefficient will be assumed in this research as 0.9 , as recommended in [9]). The relation between the active power and the nominal apparent power rate of the power converter is presented in (8).
max t Ω T { P i , t P V n o m } = η S i P V n o m .
Now, if we employ the reactive power as a function of the apparent and active power, we derive that it can be inductive or capacitive as
Q i , t P V = ± ( S i P V n o m ) 2 ( P i , t P V n o m ) 2 .
If we substitute (8) in (9), then, we have
Q i , t P V = ± 1 η max t Ω T { P i , t P V n o m } ) 2 η 2 ( P i , t P V n o m 2 .
Figure 1 shows the active and reactive power as a function of the coefficient η . This coefficient limits the quantity of active power that can be delivered and the maximum values of reactive power—either reactive or capacitive—that can take the system per hour. Note in Figure 1 that the active power is given as a function of a typical radiation curve.
Note that the solution of the optimization model defined from (1) to (10) was reached in this paper with the GAMS optimization package using the BONMIM solver that is capable of working with MINLP problems [9]. The main aspects of the GAMS optimizer are described in the next section.

3. Solution Methodology

The General Algebraic Modeling System (i.e., GAMS) is a programming environment that contains different optimization packages to solve significant problems compactly, as long as they can be represented in algebraic form [19]. This package is used in industry and academia because it contains an easy-to-implement environment, thus separating the proposed problem’s mathematical structure and the solution method. This nature of this software allows the user to focus mainly on the accurate modeling of the optimization problem being studied, rather than on the solution technique itself [19].
In the scientific literature, the GAMS optimization software has been employed for solving large-scale complex optimization problems; some of them are presented as follows: optimal planning and operation of power systems [19], optimal design of osmotic power plants [29]; optimal planning of water distribution systems [30]; optimal design of mechanical components [31]; optimal selection and location of batteries in distribution grids [7,25]; optimal location and sizing of distributed generators [32]; solution complementarity problems arising in applied economic analysis [33]; and analysis of the influence of reactive compensators and energy storage devices’ location on power systems [34], among others.
To understand the usage of the GAMS software to solve nonlinear programming problems, we present below the main aspects in the implementation of mathematical optimization models in the GAMS interface [35]:
  • Declaration of sets: Sets are a group of values of one or more elements, which in GAMS are traversed by subscripts. A set can also be made up of other sets.
  • Declaration of numerical parameters: The parameters are values not dependent on a set and will behave as constants throughout the resolution of the model. New constants can be reached by relating the mathematical operations of previously defined constants.
  • Declaration of Variables: The variables will depend on the sets and parameters defined above; accordingly, the type of the variable is declared. It is continuous or discrete, and its limits and initial values are defined. The user must define at least the variable that he wants to maximize or minimize.
  • Declaration of equations: At least one of the equations contained in the proposed model are declared, which is the objective function to minimize or maximize and which depends on the objective variable.
  • Declaration of the model and launch of the optimization program: A name is given to the program being carried out, and the type of program is determined (LP, MINLP, NLP, etc.), the maximization or minimization of the objective variable is carried out/determined, and it is also indicated if all the previously constructed equations are used.
  • Report of results: The user determines which variables he wants to print in the report of results. This report shows whether the optimal result was achieved, the value of the variables, and the number of iterations performed, among other information that the user may consider necessary.
Note that the GAMS optimization package was used in this paper to solve the mathematical model of the optimal location of PV systems in AC distribution networks considering reactive power injection. The flow diagram that summarizes the main aspects of the GAMS usage for addressing MINLP models is presented in Figure 2.

4. Test Systems

To verify the mathematical model of location and optimal sizing of photovoltaic generators in AC power systems considering reactive power injection, the variation of demand and solar power were considered as the conditions throughout this development. Additionally, two test feeders composed of 14 nodes (meshed power system topology) and 27 nodes (radial electric distribution grid) were considered for numerical validations in the software GAMS. All the information regarding the computational validations is described in the following subsections.

4.1. Demand and Solar Generation Curves

Table 1 shows the power demand and daily irradiation data, with the base powers of 8433 MW and 532.7 Wh/m 2 . These data are illustrated in Figure 3. It is worth mentioning that the renewable generation from photovoltaic technology can have important fluctuations during the year depending on the weather conditions in the geographical area that this will locate [36]. However, it is acceptable for planning purposes to assume an averaged generation profile to carry out studies about the location of PV sources as the case presented in this research, which is based on a typical Colombian PV curve.
The power demand curve is for the interconnected national system of Colombia as of 7 August 2020 [37]. The solar irradiation curve is the average power per hour for August in the measurement station of the National University of Colombia in Bogota. These data are available in [38]. Additionally, it is proposed that up to 4 PV generators will be installed to determine power losses and corresponding locations for each generator within each proposed system.

4.2. 14-Node Test System

The system of 14 nodes and three generators is used to validate the proposed model. Its parameters are shown in Table 2, and its configuration is presented in Figure 4. The base power is 20 MVA, and the base voltages are 69 kV and 13.6 kV. This system is modified according to the IEEE 14-node test system.

4.3. 27-Node Test System

Figure 5 depicts the 27-node test system, which is a system with 27 nodes and one generator. The parameters of this system are shown in Table 3. This is a typical AC distribution system, where the base power is 1 MVA, and the base voltage is 13.8 kV. This system is taken from [9].

5. Implementation and Results

The solution of the MINLP mathematical model defined in Section 2 uses the GAMS optimization package with the BONMIN solver on a laptop computer with an Intel(R) Core(TM) i5-3230M 2.6 Ghz processor and 6 GB of RAM running a 64-bit version of Windows 10 Pro.

5.1. 14-Node Test System

The optimal location and sizing of the PV systems in the 14-node test feeder is reported in Table 4. Observe that all the generators are located in demand nodes and not necessarily in those with higher demand values. It is important to note that the PV generators retain their locations and maintain similar power output as new PVs are added to the system. When PV generators are added to the system, the losses are reduced; however, this reduction becomes negligible from the third generator (see columns 3 and 4 in Table 4); this implies that the level of energy losses reduction goes to a saturation zone, where the number of the PV sources does not produce important effects on its reduction.
Figure 6 shows that regardless of the number of PV generators placed, the voltage of all nodes of the system for each hour is within limits set initially (these voltages fulfill the voltage regulation constraint (4)). Adding the PV generators flattens the minimum voltage curve compared to the base scenario and progressively increases the minimum voltage’s average value. An elevation of the maximum voltage curves of each scenario is noted, between times 6 and 18, without having a considerable difference between them; i.e., it could be considered that for this test system, only one maximum voltage curve is going to be reached, and variations are made in the minimum voltage curves in all the scenarios.
Figure 7 shows the reactive power injected by each generator in each scenario during the 24 h. In none of the scenarios, the defined power limits are surpassed. Note that in Figure 7, the PV generator located in node 6 has a slight increase in the night peak hour, where the reactive power injection follows the demand behavior curve. However, it is observed that from 2 PV, located at nodes 3, 10, and 14, the reactive power injected by them is almost negligible in each during entire period analyzed. This leads to the question of why the voltage curves have such similar data for all the placed PVs and do not have such an obvious flattening of the curves as the 14-node system. These aspects can be explained as due to the needs of the system itself to achieve the optimal value of losses.
To present the effect of the PV active and reactive power compensation in the 14-node test feeder, Figure 8 presents the daily performance of the conventional power sources connected at nodes 2, 9, and 14. From this figure, we can observe that: (i) all the conventional sources decrease their active power injections with respect to the benchmark case when PV plants start to provide active power to the grid from the period of time between 6 and 19 h; (ii) There is no period of time in which the PV plants completely replace the function of the conventional power sources, i.e., to minimize the daily energy losses, all the generators (PV and conventional power plants) provide active power to the grid; and (iii) all the conventional generators are provided during the entire period of time of reactive power to guarantee the voltage controlled output at these terminals and, at the same time, to help compensate the daily reactive power requirements of all the loads.

5.2. 27-Node Test System

Table 5 presents the power capacity and location of the various PV systems for the proposed scenarios. Similarly, as PV generators are added in previous 14-node test feeder, the power losses decrease. In this scenario, the PV generators are located at nodes distant from the conventional generator. This is expected, as the voltages at the nodes will drop as they move away from the slack node. For this reason, PV generators increase the voltages of the nodes. However, for this test system, the panels’ effect in improving the voltage is not as noticeable as in the previous cases. This is because the active and reactive power demands are lower than those in the previous cases.
It is observed that for the maximum voltage curve, the voltages do not vary considerably with respect to the base case; therefore, for this particular system, the effect of the PV generators is mainly on the minimum voltage curve and yet it does not change significantly, as seen in Figure 9.
Figure 10 shows the reactive power delivered by each PV generator in each of the suggested scenarios during 24 h of analysis. They do not exceed the established limits and have a slight elevation when in the evening peak hour of the demand curve. As a main characteristic, it is noticed that the reactive power injection tends to get considerably close to the upper power limit by the time it is at the peak value of the radiation curve.
To graphically present the effect of the PV plants in the performance of the 27-node distribution grid, Figure 11 presents the active and reactive power outputs in the slack node for all the periods under study. From the results in Figure 11, it is possible to note that: (i) the active power injections in the slack node are the same for all the periods between one to six hours and 19 to 24 h, which is expected behavior since in these periods there is no solar radiation, i.e., the possibility of injecting active power from PV sources is null, and (ii) for the cases of PV plants 3 and 4, when the solar radiation is high (periods between 11 to 14 h). The slack source does not provide active power to the grid, which implies that the PV sources supply all the active power consumption of the grid. However, the slack source provides the reactive power in all the periods to ensure the voltage-controlled output at its terminals.

5.3. Complementary Analysis

In all the proposed systems and scenarios, it is evidenced that the reactive power injection is associated with the needs of the system itself, almost following the power demand curve, with some particular exceptions. According to the results obtained from the voltage graphs in Figure 6 and Figure 9, as the number of PV generators increases, the average voltage of each test system increases, and the difference between the maximum and minimum voltage decreases, as can be seen in Figure 12 and Figure 13, respectively. In the specialized literature, this is referred to as stress profile enhancement. However, when there is no reactive power injection, this profile improvement does not appear; its values could be considered the same from one scenario to another.
The system that is less susceptible to the generators’ addition is the 27-node system since it does not improve its voltage behavior as the 14-node test feeder. This is not due to the number of nodes but the system’s topology (radial configuration). The 14-node system presents the lowest reduction of energy losses among the two systems analyzed in this research, with a maximum reduction of about 13% (Table 4). It is a test system where the reactive power injection is almost null as the PVs are added. This situation indicates the direct relationship between the reduction of energy losses and the injection of reactive power, along with the peak correction of the maximum and minimum voltage curves.
It is essential to highlight that the following aspects can influence the proposed study regarding the optimal location and sizing of PV plants in power systems: (i) the possible variations that can be experienced by the demand and solar radiation curves, (ii) the ability of the optimization tool to escape from local optimal solutions, (iii) the value assigned to the η -coefficient, and (iv) the economic and regulatory aspects that pertain to the integration of large-scale power plants in the conventional power systems. Note that these aspects can be used to improve the quality and effectiveness of the proposed MINLP model in future research.

6. Conclusions and Future Works

The problem of the optimal placement and sizing of PV plants in electrical networks from high- to medium-voltage levels was formulated in this research as an MINLP model, where the main characteristic corresponds to the possibility of using the PV converter to supply the reactive power to the power systems concurrently as the reactive power is injected. To make this apparent power compensation, the factor η was introduced in the optimization model to support reactive power from lagging or leading power factors as a function of the power system requirements. In addition, the radiation curve was characterized to ascertain the generation contributions of the various PV sources for 24 h in the context of the Colombian power system using the information as of 7 August 2020 to obtain a realistic simulation scenario of the analysis.
The solution of the proposed MINLP optimization model was addressed with the GAMS software by using the BONMIN optimization package. This software allowed the easy implementation of the model, and it was possible to separate it from the solution method.
According to the data obtained, the size was determined for each of the photovoltaic generators whose location was presented in the demand nodes of the system and, in those where the system was radial, in nodes far from the conventional generation source. It was evidenced that the reactive power injection goes hand in hand with the improvement of the voltage profile; in all the proposed systems and scenarios, it was possible to reduce energy losses between 13 and 56% according to the proposed system and the number of PV generators located. In addition, a methodology was proposed to calculate the average value and compare the improvement of the voltage profile for each system and proposed scenario, where the operational improvement of the inclusion of photovoltaic systems to electric power systems was graphically presented.
Future work is proposed to study the failure and maintenance scenarios in the PV systems and converters. The inclusion of other renewable sources within the model could be considered, such as wind generation or the contribution of electric vehicles and capacitive banks. Further, it is proposed to include economic variables to plan more appropriately for electric power systems.

Author Contributions

Conceptualization, A.F.B.-V., O.D.M., and W.G.-G.; methodology, A.F.B.-V., O.D.M., and W.G.-G.; investigation, A.F.B.-V., O.D.M., and W.G.-G.; writing—review and editing, A.F.B.-V., O.D.M., W.G.-G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

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

Acknowledgments

This work has been derived from the undergraduate project: “Localización y dimensionamiento óptimo de generadores fotovoltaicos en sistemas de potencia considerando inyección de potencia reactiva” presented by the student Andrés Felipe Buitrago Velandia to the Electrical Engineering Program of the Engineering Faculty at Universidad Distrital Francisco José de Caldas as a partial requirement for the Bachelor in Electrical Engineering.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yoon, M.; Lee, J.; Song, S.; Yoo, Y.; Jang, G.; Jung, S.; Hwang, S. Utilization of Energy Storage System for Frequency Regulation in Large-Scale Transmission System. Energies 2019, 12, 3898. [Google Scholar] [CrossRef] [Green Version]
  2. Bhatti, B.A.; Broadwater, R.; Dilek, M. Analyzing Impact of Distributed PV Generation on Integrated Transmission & Distribution System Voltage Stability—A Graph Trace Analysis Based Approach. Energies 2020, 13, 4526. [Google Scholar]
  3. Comisión de Regulación de Energía y Gas. CREG. Electrical Sector 2020. Available online: https://www.creg.gov.co/energia-electrica (accessed on 11 September 2020).
  4. Sorrentino, E.; Gupta, N. Summary of useful concepts about the coordination of directional overcurrent protections. CSEE J. Power Energy Syst. 2019, 5, 382–390. [Google Scholar] [CrossRef]
  5. Ebrahimzadeh, E.; Blaabjerg, F.; Lund, T.; Nielsen, J.G.; Kjær, P.C. Modelling and Stability Analysis of Wind Power Plants Connected to Weak Grids. Appl. Sci. 2019, 9, 4695. [Google Scholar] [CrossRef] [Green Version]
  6. Murty, P. Power Flow Studies. In Power Systems Analysis; Elsevier: Amsterdam, The Netherlands, 2017; pp. 205–276. [Google Scholar]
  7. Montoya, O.; Gil-González, W.; Hernández, J. Optimal Selection and Location of BESS Systems in Medium-Voltage Rural Distribution Networks for Minimizing Greenhouse Gas Emissions. Electronics 2020, 9, 2097. [Google Scholar] [CrossRef]
  8. Mundo-Hernández, J.; de Celis Alonso, B.; Hernández-Álvarez, J.; de Celis-Carrillo, B. An overview of solar photovoltaic energy in Mexico and Germany. Renew. Sustain. Energy Rev. 2014, 31, 639–649. [Google Scholar] [CrossRef]
  9. Gil-González, W.; Montoya, O.D.; Grisales-Noreña, L.F.; Perea-Moreno, A.J.; Hernandez-Escobedo, Q. Optimal Placement and Sizing of Wind Generators in AC Grids Considering Reactive Power Capability and Wind Speed Curves. Sustainability 2020, 12, 2983. [Google Scholar] [CrossRef] [Green Version]
  10. Streimikiene, D.; Klevas, V. Promotion of renewable energy in Baltic States. Renew. Sustain. Energy Rev. 2007, 11, 672–687. [Google Scholar] [CrossRef]
  11. Gil-González, W.; Montoya, O.; Escobar-Mejía, A.; Hernández, J. LQR-Based Adaptive Virtual Inertia for Grid Integration of Wind Energy Conversion System Based on Synchronverter Model. Electronics 2021, 9, 1022. [Google Scholar] [CrossRef]
  12. Arya, L.; Koshti, A.; Choube, S. Distributed generation planning using differential evolution accounting voltage stability consideration. Int. J. Electr. Power Energy Syst. 2012, 42, 196–207. [Google Scholar] [CrossRef]
  13. Huy, P.D.; Ramachandaramurthy, V.K.; Yong, J.Y.; Tan, K.M.; Ekanayake, J.B. Optimal placement, sizing and power factor of distributed generation: A comprehensive study spanning from the planning stage to the operation stage. Energy 2020, 195, 117011. [Google Scholar] [CrossRef]
  14. Maleki, A.; Askarzadeh, A. Artificial bee swarm optimization for optimum sizing of a stand-alone PV/WT/FC hybrid system considering LPSP concept. Sol. Energy 2014, 107, 227–235. [Google Scholar] [CrossRef]
  15. Abou El-Ela, A.; Allam, S.; Shatla, M. Maximal optimal benefits of distributed generation using genetic algorithms. Electr. Power Syst. Res. 2010, 80, 869–877. [Google Scholar] [CrossRef]
  16. Yang, N.C.; Chen, T.H. Evaluation of maximum allowable capacity of distributed generations connected to a distribution grid by dual genetic algorithm. Energy Build. 2011, 43, 3044–3052. [Google Scholar] [CrossRef]
  17. Ganguly, S.; Samajpati, D. Distributed generation allocation with on-load tap changer on radial distribution networks using adaptive genetic algorithm. Appl. Soft Comput. 2017, 59, 45–67. [Google Scholar] [CrossRef]
  18. Pesaran, H.A.M.; Nazari-Heris, M.; Mohammadi-Ivatloo, B.; Seyedi, H. A hybrid genetic particle swarm optimization for distributed generation allocation in power distribution networks. Energy 2020, 209, 118218. [Google Scholar] [CrossRef]
  19. Soroudi, A. Power System Optimization Modeling in GAMS; Springer International Publishing: Berlin/Heidelberg, Germany, 2017. [Google Scholar]
  20. Zeynali, S.; Rostami, N.; Feyzi, M. Multi-objective optimal short-term planning of renewable distributed generations and capacitor banks in power system considering different uncertainties including plug-in electric vehicles. Int. J. Electr. Power Energy Syst. 2020, 119, 105885. [Google Scholar] [CrossRef]
  21. Naval, N.; Sánchez, R.; Yusta, J.M. A virtual power plant optimal dispatch model with large and small-scale distributed renewable generation. Renew. Energy 2020, 151, 57–69. [Google Scholar] [CrossRef]
  22. Mena, R.; Hennebel, M.; Li, Y.F.; Zio, E. Self-adaptable hierarchical clustering analysis and differential evolution for optimal integration of renewable distributed generation. Appl. Energy 2014, 133, 388–402. [Google Scholar] [CrossRef] [Green Version]
  23. Quezada, C.; Torres, J.; Quizhpi, F. Optimal Location of Capacitor Banks by Implementing Heuristic Methods in Distribution Networks. In Proceedings of the 2019 IEEE CHILEAN Conference on Electrical, Electronics Engineering, Information and Communication Technologies (CHILECON), Valparaiso, Chile, 13–27 November 2019; pp. 1–6. [Google Scholar]
  24. Adewuyi, O.B.; Shigenobu, R.; Ooya, K.; Senjyu, T.; Howlader, A.M. Static voltage stability improvement with battery energy storage considering optimal control of active and reactive power injection. Electr. Power Syst. Res. 2019, 172, 303–312. [Google Scholar] [CrossRef]
  25. 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]
  26. Arellanes, A.; Rodríguez, E.; Orosco, R.; Perez, J.; Beristáin, J. Three-phase grid-tied photovoltaic inverter with reactive power compensation capability. In Proceedings of the 2017 IEEE International Autumn Meeting on Power, Electronics and Computing (ROPEC), Ixtapa, Mexico, 8–10 November 2017; pp. 1–6. [Google Scholar]
  27. Carrasco, M.; Mancilla-David, F. Maximum power point tracking algorithms for single-stage photovoltaic power plants under time-varying reactive power injection. Sol. Energy 2016, 132, 321–331. [Google Scholar] [CrossRef]
  28. Montoya, O.D.; Garcés-Ruíz, A.; Gil-González, W.; Escobar-Mejía, A. Compensación De Potencia Reactiva En Sistemas De Distribución: Un Enfoque Formal Basado en Optimización Matemática, 1st ed.; Ediciones UTB: Cartagena, Colombia, 2020. [Google Scholar]
  29. Naghiloo, A.; Abbaspour, M.; Mohammadi-Ivatloo, B.; Bakhtari, K. GAMS based approach for optimal design and sizing of a pressure retarded osmosis power plant in Bahmanshir river of Iran. Renew. Sustain. Energy Rev. 2015, 52, 1559–1565. [Google Scholar] [CrossRef]
  30. He, H.; Chen, A.; Yin, M.; Ma, Z.; You, J.; Xie, X.; Wang, Z.; An, Q. Optimal Allocation Model of Water Resources Based on the Prospect Theory. Water 2019, 11, 1289. [Google Scholar] [CrossRef] [Green Version]
  31. Andrei, N. Nonlinear Optimization Applications Using the GAMS Technology; Springer US: New York, NY, USA, 2013. [Google Scholar]
  32. Kaur, S.; Kumbhar, G.; Sharma, J. A MINLP technique for optimal placement of multiple DG units in distribution systems. Int. J. Electr. Power Energy Syst. 2014, 63, 609–617. [Google Scholar] [CrossRef]
  33. Rutherford, T.F. Extension of GAMS for complementarity problems arising in applied economic analysis. J. Econ. Dyn. Control 1995, 19, 1299–1324. [Google Scholar] [CrossRef]
  34. Ćalasan, M.; Kecojević, K.; Lukačević, O.; Ali, Z.M. Testing of influence of SVC and energy storage device’s location on power system using GAMS. In Uncertainties in Modern Power Systems; Elsevier: Amsterdam, The Netherlands, 2021; pp. 297–342. [Google Scholar]
  35. Montoya, O.D.; Gil-González, W.; Hernández, J.C.; Giral-Ramírez, D.A.; Medina-Quesada, A. A Mixed-Integer Nonlinear Programming Model for Optimal Reconfiguration of DC Distribution Feeders. Energies 2020, 13, 4440. [Google Scholar] [CrossRef]
  36. Yang, L.; Gao, X.; Li, Z.; Jia, D.; Jiang, J. Nowcasting of Surface Solar Irradiance Using FengYun-4 Satellite Observations over China. Remote Sens. 2019, 11, 1984. [Google Scholar] [CrossRef] [Green Version]
  37. XM. Real-Time Power Demand. 2020. Available online: https://www.xm.com.co.aspx (accessed on 7 March 2020).
  38. Benavides, H.; Simbaqueva, O.; Zapata, H. Atlas of Solar, Ultraviolet and Ozone Radiation of Colombia; Technical Report, Alianza entre IDEAM Y UPME; Instituto de Hidrología, Meteorología y Estudios Ambientales: Bogota, Colombia, 2018. (In Spanish) [Google Scholar]
Figure 1. Active and reactive power as a function of η during the day.
Figure 1. Active and reactive power as a function of η during the day.
Resources 10 00047 g001
Figure 2. Implementation of an MINLP model in GAMS.
Figure 2. Implementation of an MINLP model in GAMS.
Resources 10 00047 g002
Figure 3. Demand and solar radiation daily curves.
Figure 3. Demand and solar radiation daily curves.
Resources 10 00047 g003
Figure 4. The 14-node test system.
Figure 4. The 14-node test system.
Resources 10 00047 g004
Figure 5. The 27-node test system.
Figure 5. The 27-node test system.
Resources 10 00047 g005
Figure 6. Voltage limits by PV generator quantities for the 14-node system during the day.
Figure 6. Voltage limits by PV generator quantities for the 14-node system during the day.
Resources 10 00047 g006
Figure 7. Limits and reactive power injection by four PVs for the 14-node system during the day.
Figure 7. Limits and reactive power injection by four PVs for the 14-node system during the day.
Resources 10 00047 g007
Figure 8. Active and reactive power profiles: (a) generator at node 2, (b) generator at node 9, and (c) generator at node 14 during the day.
Figure 8. Active and reactive power profiles: (a) generator at node 2, (b) generator at node 9, and (c) generator at node 14 during the day.
Resources 10 00047 g008aResources 10 00047 g008b
Figure 9. Voltage limits by PV generator quantities for the 27-node system during the day.
Figure 9. Voltage limits by PV generator quantities for the 27-node system during the day.
Resources 10 00047 g009
Figure 10. Limits and reactive power injection by four PV’s for the 27-node system during the day.
Figure 10. Limits and reactive power injection by four PV’s for the 27-node system during the day.
Resources 10 00047 g010
Figure 11. Active and reactive power outputs in the slack node during the day.
Figure 11. Active and reactive power outputs in the slack node during the day.
Resources 10 00047 g011
Figure 12. Average voltage per number of PV generators.
Figure 12. Average voltage per number of PV generators.
Resources 10 00047 g012
Figure 13. Voltage difference by number of PV generators
Figure 13. Voltage difference by number of PV generators
Resources 10 00047 g013
Table 1. Demand and radiation curves for a measurement taken in Colombia.
Table 1. Demand and radiation curves for a measurement taken in Colombia.
Time (h)Demand Power (pu)Solar Rad. Power (pu)Time (h)Demand Power (pu)Solar Rad. Power (pu)
10.7650.000130.8791.000
20.7360.001140.8770.955
30.7140.001150.8680.823
40.7120.001160.8560.616
50.6980.000170.8450.352
60.6880.002180.8340.102
70.6770.102190.8340.001
80.7100.346200.9510.000
90.7530.570211.0000.000
100.7910.786220.9240.000
110.8290.920230.8690.001
120.8610.986240.7650.000
Table 2. The 14-node test system parameters.
Table 2. The 14-node test system parameters.
NodesGen. (pu)Parameters (pu)Power (pu)NodesGen. (pu)Parameters (pu)Power (pu)
i j V j R ij X ij P j Q j i j V j R ij X ij P j Q j
21-0.0040.0130.8000.100491.050.0070.0891.4750.113
421.050.0380.1171.9850.16379-0.0000.0460.0000.000
23-0.0350.1491.7100.125910-0.0100.0271.4500.129
34-0.0430.1110.7900.2196111.100.0720.1501.1750.109
15-0.0450.1891.3800.1081011-0.0600.1400.0000.000
25-0.0370.1150.0000.000612-0.1190.2490.2300.020
45-0.0000.0010.0000.000613-0.0320.0640.6360.129
56-0.0070.0891.1600.2751213-0.1680.1520.0000.000
47-0.0070.0890.9700.420914-0.1300.2780.5740.125
78-0.0000.1170.9000.1011314-0.2270.4620.0000.000
Table 3. The 27-node test system parameters.
Table 3. The 27-node test system parameters.
NodesGen. (pu)Parameters (pu)Power (pu)NodesGen. (pu)Parameters (pu)Power (pu)
i j V i R ij X ij P j Q j i j V i R ij X ij P j Q j
121.0000.000800.001040.0000.0001415-0.004600.002170.1060.066
23-0.003460.003140.0000.0001516-0.004600.002170.0250.158
34-0.001040.000940.2980.184317-0.004600.002170.2550.158
45-0.002300.001370.0000.0001718-0.002760.001300.1280.079
56-0.002560.001520.2550.1581819-0.004140.001950.2980.184
67-0.002530.001190.0000.0001920-0.004370.002060.3400.211
78-0.004600.002170.2130.1322021-0.004600.002170.0850.053
89-0.005750.002710.0000.000422-0.004600.002170.1060.066
910-0.004600.002170.2660.165523-0.004600.002170.0550.034
211-0.004600.002170.0850.053624-0.001840.000870.0700.043
1112-0.005660.002670.3400.211825-0.002760.001300.0640.040
1213-0.003450.001630.2980.1842526-0.003680.001740.1700.105
1314-0.002580.001220.1910.119827-0.002760.001300.2560.158
Table 4. Location and size of PV generators for the 14-node system.
Table 4. Location and size of PV generators for the 14-node system.
Quantity PVNodeGenerator Power MW (pu)Energy Losses (pu)Reduction (%)
0--17.187-
1669.480 (3.474)15.8547.755
2339.360 (1.968)15.38310.496
663.900 (3.195)
3339.360 (1.968)15.06712.334
658.500 (2.925)
1414.100 (0.705)
4339.360 (1.968)14.95113.009
658.480 (2.924)
1030.640 (1.532)
1414.100 (0.705)
Table 5. Location and size of PV generators for the 27-node system.
Table 5. Location and size of PV generators for the 27-node system.
Quantity PVNodeGenerator Power MW (pu)Energy Losses (pu)Reduction (%)
0--1.961-
181.740 (1.740)1.32332.534
281.479 (1.479)1.08644.620
201.013 (1.013)
381.441 (1.441)0.87555.380
140.925 (0.925)
200.982 (0.982)
481.437 (1.437)0.86056.145
130.833 (0.833)
160.188 (0.188)
200.979 (0.979)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Buitrago-Velandia, A.F.; Montoya, O.D.; Gil-González, W. Dynamic Reactive Power Compensation in Power Systems through the Optimal Siting and Sizing of Photovoltaic Sources. Resources 2021, 10, 47. https://0-doi-org.brum.beds.ac.uk/10.3390/resources10050047

AMA Style

Buitrago-Velandia AF, Montoya OD, Gil-González W. Dynamic Reactive Power Compensation in Power Systems through the Optimal Siting and Sizing of Photovoltaic Sources. Resources. 2021; 10(5):47. https://0-doi-org.brum.beds.ac.uk/10.3390/resources10050047

Chicago/Turabian Style

Buitrago-Velandia, Andrés Felipe, Oscar Danilo Montoya, and Walter Gil-González. 2021. "Dynamic Reactive Power Compensation in Power Systems through the Optimal Siting and Sizing of Photovoltaic Sources" Resources 10, no. 5: 47. https://0-doi-org.brum.beds.ac.uk/10.3390/resources10050047

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