Received 27 November 2015; accepted 25 January 2016; published 28 January 2016
In the past few years, great progress has been achieved in wind energy (WE) field as distributed generation source (DGS) because wind energy has a lot of merits such as environment, economic and technical. Environmentally it’s clean energy source and meets the requirement of reducing global warming. Technically the single wind turbine capable of generating power in few MW at present and in near future 2020 the rated power of WT may reach 10 - 20 MW  . Economically it’s the cheapest renewable energy source   , these merits make wind energy the most suitable source for future energy. Therefore, the penetration of wind turbines in electrical power systems will increase and they may begin to influence overall power system behavior, making it impossible to run a power system by only controlling large scale power plants. So it is important to study the behavior of wind turbines in an electrical power system. Refs  - were concerned in impact of connecting distributed generation (DG) to the grid on voltage stability, to study the optimal penetration level of wind farms without compromising the system voltage stability in transmission level using energy storage systems studied in ref  , while ref.  presented the impact of fixed speed WT on grid voltage stability in transmission level.
Distribution system (DS) is an important part of power system, and it represents the link between the bulk system and customers. It is estimated that as much as 13% of the total power generation is lost in the distribution networks. Approximately 70% of the total electric power system real power losses are associated with the distribution level  . Wind energy farms usually located at distribution area, so it’s essential to study the behavior associated with connecting wind turbine generation units to distribution level.
Voltage stability and power losses problems are annoying problems in DS, voltage stability may spread to transmission system and cause blackout for the whole system. Extensive studies are needed to determine the best connection of WT to DSs with acceptable voltage level.
This paper is organized as follows. Section 2 describes load flow in distribution system. Section 3 describes the wind turbine penetration level. Section 4 describes the voltage stability index used to study voltage stability margin of the system. Section 5 describes power loss reduction. Section 6 describes optimization technique used to choose optimal size and location of WT. Results obtained considering voltage stability margin, voltage profile, and power loss are presented in Sections 7 and 8. Finally, Section 9 summarizes the main conclusions.
2. Distribution System and Load Flow
Power flow problems can be solved by several methods and are classified as either direct or iterative. The direct methods employ the direct solutions related to the linear system, but actually all methods are iterative because the basic problem involves the solution of a system of non-linear equations.
Although conventional power flow methods such as Newton-Raphson (NR) technique, Fast decoupled (FD) load flow and Gauss-Seidel are well developed in dealing with the transmission and sub-transmission sections of the power system networks, they are considered to be inefficient in handling distribution networks. This is because the (DS) is different in several ways from its transmission counterpart.
§ DS has a strictly radial topology nature or weakly meshed networks in contrast with transmission systems which are “tightly” meshed networks.
§ DS is a low voltage system having low X/R ratio sections and a wide range of reactance and resistance values. The practical DSs low X/R ratio sections may cause both the NR and FD conventional methods to diverge. The line impedance angles are small enough to deteriorate the dominance of the NR Jacobian main diagonal, making it prone to singularity. Such a low X/R value would also prevent the Jacobian matrix from being decoupled and simplified.
§ DS may consist of a tremendously large number of sections and buses spread throughout the network.
All of these characteristics strongly suggest that DS is to be classified as an ill-conditioned power system.
In this paper a popular solution method used to treat ill conditioning is the forward-backward sweep (FB) method. The forward/backward method has been developed by D. Thukaram  ; this method employs Ohm’s law, Kirchhoff Voltage Law and Kirchhoff Current Law to calculate nodal voltages and branch currents.
The backward stage determines the branch currents from the end node towards the root node with a constant nodal voltage, while the forward stage is used to determine the nodal voltages by assuming the source nodal voltage is constant. The forward-backward substitution was replaced by the forward-backward sweep technique. The technique used in this paper is proposed in appendix A for more details refer to ref  .
3. Wind Turbine Penetration Level
Wind turbine Penetration Level (PL) used to study its impact on the power system operation. This PL gives percentage WT generated power contribution in load supply is given by:
PWT: wind turbine active power.
Pl: total loadsactive power.
The paper studied the impact of WT at different PL and different power factor (p.f) i.e. (lead and lag) on distribution system operation.
4. Voltage Stability
The voltage instability problem in a distribution network may spread to the transmission network causing a blackout as happened in Brazilian system in June 1997  , US/Canada blackout, the Scandinavian blackout and the Italian blackout, which all happened in 2003  . It is attracting more attention from researchers around the world.
There are several definitions of voltage stability. One definition developed by IEEE and CIGRE that Voltage stability refers to the ability of a power system to maintain steady voltages at all buses in the system after being subjected to a disturbance from a given initial operating condition  .
The term voltage collapse is also often used. It is the process by which the sequence of events accompanying voltage instability leads to a blackout or abnormally low voltages in a significant part of the power system.
Voltage Stability Indices
The main objective of Voltage Stability Indices (VSIs) is to estimate the distance from the current operating point to the system voltage marginally stable point. Numerical indices help operators to monitor how close the system is to collapse or to initiate automatic remedial action schemes to prevent voltage collapse.
Most of the VSIs that have been proposed are based on steady state power flow.
An index, which can be evaluated at all buses in radial distribution systems, was presented by M. Charkravorty and D. Das in  .
The distribution system can be represented by the equivalent of two bus system as showing in Figure 1; the Equation represents the voltage stability index SI is given by:
: Voltage stability index of bus i + 1 ().
n: the total number of buses.
i: the branch number.
r (i): resistance of branch i.
x (i): reactance of branch i.
V (i): voltage of bus i.
V (i + 1): voltage of bus i + 1.
P (i + 1): total real power load fed through bus i + 1.
Q (i + 1): total reactive power load fed through bus i + 1.
When SI approaches 0.0, that means the system is unstable.
Figure 1. Simple two bus system.
5. Power Losses
Before the WT connected to the system showing in Figure 1, the total active power losses are:
Once the WT is connected at the system, three different situations may take place:
・ situation 1:
・ situation 2:
・ situation 3:
For situations 1 and 2, the injection of PWT reduces the Pi drawn from bus i. In other words, the WT reduces the load’s value at bus i + 1. However, the WT in situation 3 acts as in situations 1 and 2, plus it will reverse a portion of the injected PWT back to bus i. Reducing the power loss and maintaining the voltage within acceptable limits are vital in WT applications. Although the first and second situations’ results in some cases are significant, situation 3 leads to better results if the optimal WT size and location are employed.
To evaluate the impact of installing WT at distribution feeder on power loss reduction (PLR) the following relation used:
PLOSS: System losses without WT
: System loss after WT installation
6. Optimal Size and Location of WT
Optimal location and size of WT generation obtained using Grid Search Algorithm  - . The Grid Search Algorithm is applied by adding WT to each bus, changing the size of WT from 0% of total load power to 100% of total load power with the step size of 0.05% of total load power. In this algorithm, the objective is to minimize total power losses of the system (PlossT) by injected active power of WT (PWT) for WT placement The main constraints are to restrict the voltages along the radial system within 1 ± 0.05 pu, as in Equation (5-a) and Equation (5-b).
To minimize subject to
The flowchart of the Grid Search algorithm to determine the optimum size and location of WT is proposed in the appendix B Figure B-1.
7. First Case Study
The first test case is a 9-bus   , single feeder, radial distribution system shown in Figure 2. This system has no laterals. The rated line voltage of the system is 23 kV. The details of the feeder and the load characteristics are given in Table C-1 in Appendix C.
Using FB sweep analysis SI index and PLR are evaluated for different PL of WT.
Wind turbine installed with different penetration levels at bus 9 (weakest branch of distribution feeder based on the value of the VSI). It’s found that as the penetration level increase system voltage stability enhanced and power loss decreased as shown in the following Figures 3-9.
Figure 2. Single line diagram of 9-bus feeder.
Figure 3. Bus voltages at different PLs (0.9 lead p.f.).
Figure 4. λ-V curves at bus-9 with different PLs (lead p.f.).
Figure 5. Bus voltages at different PLs (0.9 lag p.f.).
Figure 6. SI index at different PLs (0.9 lead p.f.).
Figure 7. SI index at different PLs (0.9 lag p.f.).
Figure 8. λ-V curves at bus-9 with different PLs (lag p.f.).
Results show that the impact of WT also depends on its power factor i.e. (in case of lead p.f enhance voltage stability and PLR more than lag p.f case).
WT installed at bus 9 with different PLs listed in following Table 1 for both 0.9 lead and lag p.f.
The Figures 3-10 show the impact of connection WT with different sizes and locations on voltage profile, SI index, λ-V curves, active losses, and reactive losses; while Figure 11 shows the relation between penetration level and power loss of the system.
Tables 2-8 give the numerical values of the impact of connection WT with different sizes and locations on the system. While Table 9 and Table 10 gives the optimal size and location to install WT in the system.
Optimal penetration level at bus 9 to minimize the total system losses is shown in Figure 11.
8. Second Case Study
Another 33-bus radial distribution test system has been used   . This test system has a main feeder and four laterals (sub-feeders). The data of the feeder is presented in Table C-2 in appendix C, the single line diagram shown in Figure 12. The system rated voltage is 11 kV.
Study the voltage stability enhancement and PLR by installing WT at bus 5 of distribution feeder with different PL% listed in Table 11.
The Figures 13-20 show the impact of connection WT with different sizes and locations on voltage profile, SI index, λ-V curves, active losses, and reactive losses; while Figure 21 shows the relation between penetration level and power loss of the system.
Figure 9. Power losses of 9-bus feeder at different PLs (lead p.f.).
Table 1. Different PLs of WT.
Figure 10. Power losses of 9-bus feeder at different PLs (lag p.f.).
Figure 11. Total system real power loss at different PLs.
Table 2. Bus Voltages at different PLs (0.9 lead p.f.).
Table 3. Bus voltages at different PLs (0.9 lag p.f.).
Table 4. SI index at different PLs (0.9 lead p.f.).
Table 5. SI index at different PLs (0.9 lag p.f.).
Table 6. Power losses of 9-bus feeder at different PLs (lead p.f.).
Table 7. Power losses of 9-bus feeder at different PLs (lag p.f.).
Table 8. Voltage and voltage index of bus-9 in all cases.
Table 9. Optimal size of WT at various buses.
Table 10. Optimal location for installing WT on system.
Table 11. Different PLs of WT.
Figure 12. Single line diagram of feeder 33-bus laterals and sub-lateral.
Figure 13. Bus voltages at different PLs (0.9 lead p.f.).
Figure 14. Bus voltages at different PLs (0.9 lag p.f.).
Figure 15. SI index at different PLs (0.9 lead p.f.).
Figure 16. SI index at different PLs (0.9 lag p.f.).
Figure 17. λ-V curves at bus-5 with different PLs (lead p.f.).
Figure 18. λ-V curves at bus-5 with different PLs (lag p.f.).
Figure 19. Power loss of 33-bus at different PLs (lead p.f.).
Figure 20. Power loss of 33-bus at different PLs (lag p.f.).
Figure 21. Total system real power losses at different PLs.
Tables 12-18 give the numerical values of the impact of connection WT with different sizes and locations on the system, while Table 19 and Table 20 give the optimal size and location to install WT in the system.
Optimal penetration level at bus 5 to minimize the total system losses is 110% as shown in Figure 21.
Table 12. Bus voltages at different PLs (0.9 lead p.f.).
Table 13. Bus voltages at different PLs (0.9 lag p.f.).
Table 14. SI index at different PLs (0.9 lead p.f.).
Table 15. SI index at different PLs (0.9 lag p.f.).
Table 16. Power loss of 33-bus at different PLs (lead p.f.).
Table 17. Power loss of 33-bus at different PLs (lag p.f.).
Table 18. Voltage and voltage index of bus-5 in all cases.
Table 19. Optimal location for installing WT on the system.
Table 20. Optimal size of WT at various buses.
The paper studied the impact of installing wind turbine in radial distribution system with different PLs on voltage stability and power loss reduction, the optimum size and location of wind turbine are also determined in this paper. The study is applied on the 9-bus and 33-bus test systems. The results show that as penetration level increases voltage stability and power loss reduction are enhanced especially with lead p.f. operation. Voltage stability index reflects a guide to choose the suitable bus for installing wind turbines DG, and these results are supported by PLR and load margin of the system.
Power flow solution
A). Power Flow Solution for Single Line Feeders
Using the system in Figure A1, the solution steps are summarized as follows:
1). Find the sum of active and reactive power loads for all buses as well as the sum of all resistances and inductive reactances of each section connecting two buses.
2). Assume that sending end real power, reactive power and voltage to be approximated by:
V0 = l + j0 p.u (3-A)
where the incremental increases Pfactor and Qfactor are given by:
Equation (4-A) and (5-A) are due to the fact that in reality, there are no lossless systems and there is rarely a feeder with only two buses. These approximation factors are added to reduce the required number of iterations required for a solution i.e., get closer to the exact loss values, then, use these values P0&Q0 in the initial iteration.
3). Apply the following power flow equations to the feeder.
Pfactor: Real power loss approximation.
Qfactor: Reactive power loss approximation.
μp: Active power multiplier, set to zero when there is no active power source or set to 1 when there is active power source.
Figure A1. Radial distribution feeder model including wind turbine.
μ0: Reactive power multiplier, set to zero when there is no reactive power source or set to ± 1 when there is a reactive power source.
APi + 1:Active Power magnitude injected at bus i + 1
RPi + 1:Reactive Power magnitude injected at bus i + 1
4). If the absolute values of Pn & Qn of the last bus are zero or within an acceptable tolerance, the power flow solution is acceptable. Otherwise, go to the next step.
5). For the first bus in the main feeder set:
Then, use these new initial values and repeat step 3.
A). Power Flow Solution for Lateral and Sub-lateral Feeders:
1). Find the sum of real and reactive power loads on the sub-laterals (if available) and represent them as loads on their laterals. Do the same for laterals and represent them as loads on the main feeder.
2). Apply the recursive power flow solution algorithm as in (single line feeders) for the main feeder.
3). Set the voltage of the far bus that represents a lateral to be equal to the (V0) of this lateral. Solve the lateral individually as described in the case of the single line feeder.
4). For sub-laterals (if available), set the voltage of the far bus that represents a sub-lateral to be equal to the (V0) of this sub-lateral. Solve the sub-lateral individually as described in the case of the single line feeder.
5). Get the total real and reactive powers injected into the sub-lateral and represent them again as a single load on its lateral.
6). Apply the power flow solution to the lateral individually again. If there is another sub-lateral go back to step 4 and solve the second far bus on the lateral that represents the sub-lateral. Otherwise, go to next step.
7). Find the total real and reactive powers injected into the lateral (use data from last power flow run) and represent them again as a load on the main feeder.
8). Apply the power flow recursion to the main feeder. If there are more laterals, go back to step 3 and solve for the second far bus that represents the lateral on the main feeder. If not, go to the next step.
9). With the bus voltages found using the previous step, solve the laterals and sub-laterals (if available) individually. Then the power flow solution is reached.
Figure B-1.Flowchart of the grid search algorithm.
Table C-1.9-Bus feeder systems data.
Table C-2. 33-Bus feeder systems data.