International Journal of Electrical and Computer Engineering (IJECE) Vol.
No.
February 2019, pp.
ISSN: 2088-8708.
DOI: 10.
11591/ijece.
Voltage sensitivity analysis to determine the optimal integration of distributed generation in distribution systems Katherine Cabana1.
John Candelo2.
Rafael Castillo3.
Emiro De-la-Hoz-Franco4 1,4Department of Computer Science and Electronics.
Universidad de la Costa.
Colombia 2Department of Electrical and Automatic Energy.
Universidad Nacional de Colombia Sede Medellin.
Colombia 3Department of Electrical and Electronic Engineering.
Universidad del Norte.
Colombia
Article Info
ABSTRACT
Article history:
This paper presents a voltage sensitivity analysis with respect to the real power injected with renewable energies to determine the optimal integration of distributed generation (DG) in distribution systems (DS).
The best nodes where the power injections improve voltages magnitudes complying with the constraints are determined.
As it is a combinatorial problem, particle swarm optimization (PSO) and simulated annealing (SA) were used to change injections from 10% to 60% of the total power load using solar and wind generators and find the candidate nodes for installing power sources.
The method was tested using the 33-node, 69-node and 118-node radial distribution networks.
The results showed that the best nodes for injecting real power with renewable energies were selected for the distribution network by using the voltage sensitivity analysis.
Algorithms found the best nodes for the three radial distribution networks with similar values in the maximum injection of real power, suggesting that this value maintains for all the power system cases.
The injections applied to the different nodes showed that voltage magnitudes increase significantly, especially when exceeding the maximum penetration of DG.
The test showed that some nodes support injections up to the limits, but the voltages increase considerably on all Received Feb 22, 2018 Revised Jul 19, 2018 Accepted Sep 8, 2018 Keywords:
Distributed generation Distribution networks Metaheuristic algorithms Sensitivity analysis Voltage magnitudes Copyright A 2019 Institute of Advanced Engineering and Science.
All rights reserved.
Corresponding Author:
Katherine Cabana.
Department of Computer Science and Electronics Universidad de la Costa.
CUC, 58 # 55 - 66.
Barranquilla.
Colombia.
Email: kcabana@cuc.
INTRODUCTION
The integration of DG in distribution network has increased rapidly .
, since the possibility of generating power close to consumers, reducing power losses, increasing voltages, improving reliability, and others .
, .
However, the large penetration of electric sources brings stability problems .
, .
For this reason, the voltage stability must be study as the increasing penetration of renewable-energy demand high levels of reactive power in the power grid and voltage support is a major challenge.
Ae.
According to .
, distribution networks have no stability problems in their normal design, but the integration of a large number of DG brings some stability issues.
Some previous studies have determine that the connection of DG may affect the voltage stability in distribution systems .
, .
In .
, the researchers presented the particle swarm optimization (PSO) as an alternative to optimize a power injection model from DG, with the aim of maximizing voltage stability.
In .
, a combination of evolutionary programming (EP) and PSO was presented to achieve convergence and accuracy of determining faster DG size and location.
Furthermore, the authors in .
have proposed the optimal placement, size and number of different types of DG units in distribution systems considering the voltage limits and the linesAo transfer capacities, using the Journal homepage: http://iaescore.
com/journals/index.
php/IJECE A ISSN: 2088-8708 genetic algorithm (GA) as an optimization technique and the backward/forward sweep method (BFS).
The authors in .
investigated the evaluation of static voltage stability on Ie 33-bus.
PG&E 69- bus and a real case with two stochastic DG units.
In the literature, the studies presented focus on determining how the integration affects the voltage stability of the distribution network and the maximum value of integrating the generators .
However, some of these studies have focused on the optimal location of distributed generation to minimize power loss or cost functions, not analysing the voltage sensitivities of the nodes.
And some studies uses metaheuristics applying the same objective functions to determine the best place for the power sources, but not other analysis such as voltage changes are considered to evaluate the possible integration of the renewable energy sources without affecting the operating conditions.
This paper focuses on determining the maximum integration of distributed generation in distribution networks by using the power flow iteration process integrated to the metaheuristic algorithms to search a better place for the power sources and the maximum power suported by the network without affecting the voltage conditions of the network.
Because this is a large combinatorial problem, we used particle swarm optimization (PSO) and simulated annealing (SA) as the optimization techniques to determine the place and the size of DG.
As the objective function we used the improvement of the voltage voltage profile in distribution networks.
The approach is applied on Ie 33-bus.
Ie 69-bus and Ie 118-bus radial distribution system.
To achieve this.
Section 2 describes the methodology of the research.
Section 3 presents the review and implementation of the methodology in a case study.
Section 4 presents a brief discussion about distributed generation effects.
Finally, the major contributions and conclusions of the papers are RESEARCH METHOD Figure 1 summarizes the method used in this research to determine the candidate nodes for installing DG.
The proposed method starts reading the input data from the models of the distribution network.
Then, an initial population is defined and evaluated.
Finally, the location and size of DG is determined during the iteration of the algorithms.
Start Parameters of the algorithms Input data Distribution system model Solar and wind models Initial population Probabilistic model of DG Fitness evaluation Save the best solution Find a new solution Evaluate the new solution Save the new best solution Yes Stoping criteria? End Location and size of DG Figure 1.
Flow chart of the methodology Int J Elec & Comp Eng.
Vol.
No.
February 2019 : 55 - 65
Int J Elec & Comp Eng
ISSN: 2088-8708
Nodes for injecting power Finding the node to inject power in the distribution network is important for this study, because it defines the best place to install the DG, maintaining good operating conditions of the system under study.
For this purpose, the following three stages were proposed in the analysis to determine the best nodes to inject power in the distribution networks .
Stage 1: distribution networks must be modelled to identify voltage sensitivities of different nodes when installing DG.
A static model of each power plant under analysis is considered .
The load is modelled as constant value maintaining the same power factor and the main source was considered as the slack node.
Stage 2: the operating state of the distribution network is found, determining the voltage sensitivities of nodes with respect to the real power injected .
Stage 3: candidate nodes are selected using power flows and analyzed with voltage sensitivities before and after injecting real power.
In this research, different scenarios of demand are used to identify weak nodes.
Nodes are classified according to the voltage variations when the injection of real power increase.
Scenarios for testing the method The following scenarios are considered to test the method and analyze the results:
Scenario 1.
In this first scenario, no DG is considered for the simulation, being a base case to compare other results of the integration of solar and wind energy generation.
Scenario 2.
In the second scenario, photovoltaic energy systems are integrated to the distribution network to determine the maximum power injection.
Scenario 3.
In the third scenario, wind energy systems are integrated to the distribution network to determine the maximum power injection.
Scenario 4.
In the fourth scenario, photovoltaic and wind energy systems are integrated to the distribution network to determine the maximum power injection.
Considerations For testing the method, the following considerations are stated:
More than one generator can be installed at each node.
All DG units operate with a unitary power factor, to avoid interference with voltage control devices connected to the system .
, .
In addition, the power factor remains constant for all tests.
Wind speed and solar radiation have the same values for all the points where the generators can be Some researchers recommend that DG penetration levels are equal to or less than 30% of the maximum load .
However, in this study it is up to 60%.
The increasing power steps are defined as 0, 10, 20, 30, 40, 50, and 60%.
The models of the elements proposed in .
are used in this research.
Load is modelled for the four seasons of the year, as considered in other studies .
The power demand of each distribution network is considered as the peak value.
Solar radiation and wind speed are modeled using the Beta and Weibull probability density functions, respectively, as considered in other studies .
, .
DG units are installed on a given node and the voltage changes are monitored.
Each generator supplies a constant power of 4,5MW, with unity power factor.
The minimum and maximum voltage values for all the distribution networks were defined as Vmin=0.
and Vmax=1.
1 p.
u, respectively.
Load model Table 1 presents the load profile, as percentages of the annual maximum load .
Annual maximum load demand is 16.
18 MVA.
Data is used to model solar radiation and wind speed with the Beta and Weibull probability functions, respectively.
Voltage sensitivity analysis to determine the optimal integration of distributed.
(Katherine Caban.
A ISSN: 2088-8708 Table 1.
Data Considered for the Loads of the Distribution Network .
Hours 12Ai1 am 1Ai2 2Ai3 3Ai4 4Ai5 5Ai6 6Ai7 7Ai8 8Ai9 9Ai10 10Ai11 11Ai2 pm 12Ai1 1Ai2 2Ai3 3Ai4 4Ai5 5Ai6 6Ai7 7Ai8 8Ai9 9Ai10 10Ai11 11Ai2 am Winter Spring Summer Fall Solar model Solar radiation was modeled using the beta probability density function .
, .
, as shown in .
Where yceyca .
cI) is the beta probability density function.
S is the solar radiation in kW/m2, considering that 0 O S O 1.
and are parameters of the distribution function, considering that > 0 and > 0.
u y.
cI) = { O ycI .
uOe.
O .
Oe ycI).
uOe.
Parameter yu can be calculated using .
Where is the mean distribution and yua is the standard deviation of the distribution function.
yu = .
Oe yuN) O ( yuN O .
yuN) Oe .
Parameter yua can be calculated using the average distribution parameter yuN and the parameter yu, as shown in .
yua= yuNyu 1OeyuN Wind model Wind speed variations can be described using the Weibull probability density function as shown in .
Where k is a shape parameter and c is a scale parameter.
When k is equal to 2, the probability density function is called Rayleigh yceyc .
as shown in .
Parameter yu was consider equal to 2.
02 and parameter equal to 9 .
, .
The scale parameter of the Rayleigh probability density function can be approximated as c=1.
128*Vm.
= yco yc ycoOe1
( )
yaycUycE [Oe [ ] ] yca yca yca 2yc yc2 yceyc .
= ( 2 ) yaycUycE [Oe [ ] ] yca yca Int J Elec & Comp Eng.
Vol.
No.
February 2019 : 55 - 65
Int J Elec & Comp Eng
ISSN: 2088-8708
Problem formulation The general model of the power system can be represented by the function yce.
cu, yuI), as shown in .
Where x represents the state variables and the load factor.
cu, yuI) = 0 When the load increases, the system can suffer variations in stress levels.
The variation in real and reactive power can be represented as shown in .
Where ycEycn and ycEycn represent the real and reactive power of the node i after changing yuI, respectively.
ycE0,ycn and ycE0,ycn are the initial real and reactive power of the node i, respectively.
yaycy,ycn and yayc,ycn are vectors that indicate the increasing power of node i, and OIyuI represents the variation of the load.
ycEycn = ycE0,ycn .
yaycy,ycn .
OIyuI) .
ycEycn = ycE0,ycn .
yayc,ycn .
OIyuI) .
The real and reactive power values of each node i, can be calculated as shown in .
Where n is the number of nodes, .
cOycn | represents the voltage magnitude of the node i, yuycn represents the voltage angle of the node i, .
cUycnyc | is the admittance magnitude of the element .
and yuEycnyc is the impedance angle of the element .
ycEycn = O.
ycOycn |.
cOyc |.
cUycnyc .
uEycnyc Oeyuycn yuyc ) yc=1 ycu ycEycn = O.
ycOycn |.
cOyc |.
cUycnyc .
uEycnyc Oeyuycn yuyc ) yc=1 The real power changes of the different generators can be modelled by the expression presented in .
Where ycEya0 the initial real power of each generation unit and represents the power changing The variation of is performed between zero .
harge nominal syste.
and the maximum value of convergence .
ycEya = .
yuI)ycEya0 The space of .
) variation in this research is limited between 0.
5 and 1.
5 times the load base.
The increase in the level of charge is carried out with same value of , for all nodes.
Objective function The objective function is defined to increase the generation at different nodes according to the voltage magnitudes of the network, as shown in .
Where ycOyco is the voltage at the selected node of the current scenario.
ycOycoycaycaycyce is the voltage of the load node in the previous scenario.
Where a high value of ycOycA indicates an excellent location of DG in terms of the voltage magnitudes.
At the maximum voltage values, the real or reactive powers are maximized .
ycAycaycuycnycoycnycyce ycOycA = ycOyco .
Oe .
cOyco Oe ycOycoycaycaycyce )] .
Electrical constraints The real power generated at node i, ycEycn , must be limited by the maximum and minimum value, as shown in .
Where ycEycn,ycoycnycu and ycEycn,ycoycaycu represent the maximum and minimum real power limits generated at node i, respectively.
ycEycn,ycoycnycu O ycEycn O ycEycn,ycoycaycu ycn = .
, .
A , ycu The reactive power generated at node i, ycEycn , must be restricted by the maximum and minimum value, as shown in .
Where ycEycn,ycoycnycu and ycEycn,ycoycaycu represents the maximum and minimum reactive power limits generated at node ycn, respectively.
Voltage sensitivity analysis to determine the optimal integration of distributed.
(Katherine Caban.
A ISSN: 2088-8708
ycEycn,ycoycnycu O ycEycn O ycEycn,ycoycaycu ycn = .
, .
A , ycu Voltage magnitude, ycOycn , of each node must be limited according to the maximum and minimum values, as presented in .
Where ycOycn,ycoycnycu ycaycuycc ycOycn,ycoycaycu represents the maximum and minimum voltage magnitude limits of the node ycn, respectively.
The slack node is assumed to have a voltage magnitude of 1 p.
ycOycn,ycoycnycu O ycOycn O ycOycn,ycoycnycu Test cases and simulations Three radial distribution systems were selected to test the method .
Ae.
The 33-node radial distribution network .
, .
has 33 nodes, 32 lines, 1 main source, and 32 loads.
The total load of the network is 3715 kW and 2300 kVAr and the total power supply of 3926 KW and 2443 kVAr.
The 69-node radial distribution network .
Ae.
has 69 nodes, 68 lines, 1 main source, and 49 loads.
The total load of the network is 4014 kW and 2845 kVAr and the total generation of 4265 kW and 2957 kVAr.
And finally, the 118-node radial distribution network .
has 118 nodes, 117 lines, 1 main source, and 117 loads.
The total load of the network is 22709 kW and 17041 kVAr and the total generation of 24000 kW and 18019 kVAr.
RESULTS AND ANALYSIS
1 Node selection Table 2 shows the results obtained when locating and sizing different generators in the distribution network with the objective function studied.
The first column is the distribution network test case, the second column is the node selected for installing DG, and the other columns correspond to the results obtained with the algorithm testing the four scenarios.
The higher power injection is obtained for the nodes away from the main source.
Additionally, from the table we can conclude that the algorithms find the same solutions, but SA have a large time to converge for the solution compared to the PSO.
The number of iterations of the SA is greater than the used with the PSO.
The voltage sensitivity analysis found with PSO and SA presented a correlation coefficient of 0.
9, indicating a strong and positive correlation between the data obtained by each Table 2.
Location and Size of DG with PSO and SA
Power System
Ie 33
Ie 69
Ie 118
Nodes PSO
PSO
Scenario (MW) PSO
PSO
Figure 2 shows the voltage sensitivity analysis when the power injected varies from 10% to 60% of the total load.
When the power injection exceeds the 30%, the voltage sensitivity values increase and separate from the initial values.
The voltages have an exponential increase and the power flow shows that some nodes are overloaded.
The simulations show a similar result obtained in previous results related to not exceed 30% of the power load .
, .
Int J Elec & Comp Eng.
Vol.
No.
February 2019 : 55 - 65
Int J Elec & Comp Eng
ISSN: 2088-8708
Figure 2.
Voltage sensitivity with different power injection 2 Sensitivity analysis of the 33-node test case Figure 3 presents the voltage sensitivity analysis of the 33-node radial distribution network.
The real power was injected changing from 10% to 60% of the total load in steps of 10%.
The axis y represents the change in voltage magnitudes with respect to the change in real power.
The axis x represents the node of the distribution network.
When the real power injection reaches 40%, the voltage changes The changes in voltage represents a large variation for all nodes of the distribution network.
This result confirms the maximum levels of DG penetration for the 30% of the maximum load .
Nodes selected by the algorithms .
, 30 and .
are the most sensitive of the network, changing voltages Furthermore, a node randomly selected .
shows that the maximum percentage has a different behavior in voltages of all nodes.
POWER INJECTION AT NODE 6
POWER INJECTION AT NODE 29
Voltage sensitivity POWER INJECTION AT NODE 30
POWER INJECTION AT NODE 32
10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33
Node Figure 3.
Voltage sensitivity analysis iV/iP for the 33-node radial distribution network Voltage sensitivity analysis to determine the optimal integration of distributed.
(Katherine Caban.
A ISSN: 2088-8708 3 Sensitivity analysis of the 69-node test case Figure 4 presents the voltage sensitivity analysis of the nodes in the 69-node radial distribution The real power was injected changing from 10% to 60% of the total load in steps of 10%.
The axis y represents the change in voltage magnitude with respect to the change in real power injection.
The nodes selected by the algorithms are 19, 25 and 68, and the node 7 was selected randomly to compare the voltage Voltage magnitude variations are confirmed, especially when the generation is larger than 30% of the load.
The voltage variations are similar for all nodes.
When the penetration is higher the voltages are largely increased.
Voltage sensitivity POWER INJECTION AT NODE 7
POWER INJECTION AT NODE 19
POWER INJECTION AT NODE 25
POWER INJECTION AT NODE 68
8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40 42 44 46 48 50 52 54 56 58 60 62 64 66 68
Node Figure 4.
Voltage sensitivity analysis iV/iP for the 69-node radial distribution network 4 Sensitivity analysis of the 118-node test case Figure 5 presents the voltage sensitivity analysis of the nodes in the 118-node radial distribution The real power was injected changing from 10% to 60% of the total load in steps of 10%.
The axis y represents the voltage magnitude change with respect to the real power change.
The axis x represents the node of the distribution network.
The PSO and SA selected the nodes 21, 22 and 117.
Int J Elec & Comp Eng.
Vol.
No.
February 2019 : 55 - 65
Int J Elec & Comp Eng
Voltage sensitivity POWER INJECTION AT NODE 8
POWER INJECTION AT NODE 21
POWER INJECTION AT NODE 22
POWER INJECTION AT NODE 117
ISSN: 2088-8708
Node Figure 5.
Voltage sensitivity analysis iV/iP for the 118-node radial distribution network Additionally, the node 8 was selected to compare the results with the best nodes for the power Similar to the previous results, the nodes selected show a high increase of voltage magnitudes, especially when the power injections are greater than 30% of the total load.
The node 8 presents a low voltage changes after the power injections and the variation do not represent significant changes on voltage magnitudes of all nodes.
CONCLUSION
This article presented the placement and size of DG in distribution systems using voltage sensitivity PSO and SA were used in this research to identify the nodes that accept the maximum real power The algorithms identified well the nodes for power injection, but PSO was faster than the SA.
Nodes away from the main source can increase more the voltage magnitudes and are more likely to be selected with the model applied.
The model presented in this paper showed how to place different types of DG in distribution systems to improve voltage profiles with a good percentage of success.
The results of the simulations show the location and size of the injections of power to impact positively on the system.
Voltage sensitivity analysis to determine the optimal integration of distributed.
(Katherine Caban.
A ISSN: 2088-8708 The model showed a good percentage of success when locating different types of DG to improve voltage magnitudes in the three distribution networks.
Location and size of power sources impact positively the radial network and after repeating all the test the results were similar, although in all the analysis performed the PSO is faster than the SA.
REFERENCES