International Journal Of Computational Engineering Research (ijceronline.com) Vo2. 2 Issue.

Load Flow Analysis of Distribution System Including Wind Turbine Generating System Models
1,2

P.Srihari1, G.Srinivasa Rao2

Equipments are underpinning this strategy. Wind energy is converted into mechanical power by wind turbines. This mechanical power is fed to the generator directly or through gear system. This cumulative arrangement of turbine and Generator is called a wind turbine generating unit (WTGU). A large number of wind turbine generating systems (WTGSs) are already in operation and many new systems are being planned. The integration of WTGS with the power systems is taking place at both the transmission and the distribution voltage levels. This growth in wind generation has spurred investigations, to understand the behavior of the wind turbine generating Systems (WTGSs) as well as their impact on the power grid. As distributed generation penetrates into the distribution network its impact on the voltage profile and total power losses needs to be analyzed by appropriate load flow methods. Due to high R/X ratio of distribution lines, the performance of load flow methods used for the transmission network are inadequate for distribution systems. These methods do not make use of the structural advantage (radial / weakly mesh) of distribution system. This led to the development of separate load flow methods for distribution system. These methods are broadly classified into three categories: direct methods, backward / forward sweep methods and Newton-Raphson based methods [1]. The backward/forward sweep methods have an edge over the others as they are simple, flexible and the computational time requirement is low [1]. The basic backward/forward sweep method described in [2] is applicable for radial networks only. The practical distribution systems are weakly meshed. The compensation based distribution load flow algorithm proposed in [3] address this issue. It is improved to adapt the PV bus modeling in [4]. The method was further extended to the distribution networks containing DGs, where these are modeled as constant PQ or PV nodes [5].

,Post Graduate Student , Department of EEE, V R Siddhartha Engg. College (Autonomous), Vijayawada-7.

Abstract The power world is sauntering towards ecofriendly distributed generation (DG), their integration into the distribution network postures challenges to existing load flow techniques. This paper is concerned with developing models of various types of wind turbine generating systems (WTGSs) used as distributed generation (DG) sources and demonstrating their application for steady state analysis. As wind energy is used as a DG power source, its impact on the distribution system needs to be analyzed. This requires a load flow algorithm capable of simulating both radial and weakly mesh systems along with Wind Turbine Generating Systems (WTGSs). The compensation based distribution load flow algorithm is extended including WTGS modeling in the present work. The application of the proposed models for the load flow analysis of radial systems having WTGS has been demonstrated. Based on these studies we evaluate the impact of wind based DG on the voltage profile and losses of radial distribution networks. Simulation studies have been conceded out on radial distribution systems having WTGS as DG sources to illustrate the application of the proposed models.

Keywords Distributed Generation, Distribution load flow, Wind Turbine Generating Systems, Forward and Backward sweep, Voltage Profile, Losses. I.INTRODUCTION
The supply of electric power to homes, offices, schools, factories, stores, and nearly every other place in the modern world is now taken for granted. Electric power has become a fundamental part of the infrastructure of contemporary society, with most of todays daily activity based on the assumption that the desired electric power is readily available. The power systems which provide this electricity are some of the largest and most complex systems in the world. They consist of three primary components: the generation system, the transmission system, and the distribution system. Each component is essential to the process of delivering power from the site where it is produced to the customer who uses it. Wind energy continues to be one of the fastest growing energy technologies and it looks set to become a major generator of electricity throughout the world. The role of wind energy in electrical power generation is increasing day by day and new advances in power electronic

Ii. Wtgs Modeling

The mechanical power converted from wind energy by a wind turbine is described as (1a)

Issn 2250-3005(online)

September| 2012

Page 1521

International Journal Of Computational Engineering Research (ijceronline.com) Vo2. 2 Issue. 5

(1b) Where, Pm is the mechanical power developed in watts, is density of air (kg/ ), V is wind speed(m/s), is power coefficient, is pitch angle, is tip speed ratio, is gear ratio and A is the area swept by the rotor is the angular velocity(rpm),R is the turbine rotor radius(m) and ,x are constants. The variation of power coefficient with a variation of , is non linear in nature. It can be estimated by the relation. (1c) Where, (2)

From the equivalent circuit, given, in Fig.1,the reactive power balance can be written as follows (4) Substituting the voltage equation, Eq.(3) into Eq(4), we get (5) Similarly, from the equivalent circuit, the active power balance can be written as (6) And, substituting voltage equation, Eq.(3), into Eq.(6), we get the active power demand of the machine as follows (7) Equ.s (5) &(7) have a straight solution and depend on the terminal voltage magnitude, mechanical power input and the equivalent circuit [6-8].

The induction generators are generally used in wind turbine for electricity generation and modeled with standard equivalent circuit from using the concept of a rotating transformer [69]. In the study of [6], Feijoo and Cidras proposed two models based on the steady- state equivalent circuit of the induction generator. In the study, wind farms with asynchronous generators are modeled as PQ or RX buses in the load ow analysis. When the WTGS node is modeled as a PQ-bus, produced active and consumed reactive power. The other model, in that study, is called RX bus model in which active and reactive powers are calculated by using equivalent circuit parameters of the induction machine. Both are solved iteratively using a Newton type algorithm. Implementation of RX model into distribution system load ow analysis is found to be problematic, therefore, excluded from this study and the PQ model is used for the analyses. Model-I for WTGS: Let us consider an induction motor equivalent circuit referred to the stator side as given in Fig. 1. In the gure, Vs and Vr stand for the magnitude of the stator voltage and referred rotor voltage, respectively. Rsc and Xsc show short-circuit equivalent resistance and reactance, respectively. Rm resistance represents core losses and Xm shows magnetizing reactance of the motor. Ps and Q s stand for the active and reactive power demand of the motor, respectively. From gure, recognizing the fact that Qr = 0, the referred rotor voltage (Vr) can be obtained by solving well-known biquadratic equation given in [10] by neglecting the short-circuit losses as follows. (3)

Fig.1. Induction motor equivalent circuit. Model-II for WTGS: Similarly, the referred rotor voltage (Vr) of the induction motor, given in Fig. 1, can also be obtained by using modied bi-quadratic equation given in [12] as follows (8a) (8b)

(8c) From Fig.1,recognizing the fact that =0,the referred rotor voltage can be written as follows (9a) (9b)
September| 2012 Page 1522

Issn 2250-3005(online)

International Journal Of Computational Engineering Research (ijceronline.com) Vo2. 2 Issue. 5

(9c) Substituting the voltage equation, Eq.(9c) into Eq.(4), we get

power output of the WTGU for a given terminal voltage and wind speed. Pitch regulated fixed speed WTGU[7] In this class of WTGU, the pitch angle controller regulates the wind turbine blade angle ( ) according to the wind speed variations. This designed power output Pe of the WTGU with wind speed is provided by the manufacturer in the form of a power curve. Hence, for a given wind speed Pe can be obtained from the power curve of the WTGU, but Qe needs to be computed. With the slip known, the reactive power output Qe is calculated from the induction generator equivalent circuit. Pitch regulated fixed speed WTGU power output (Pe and Qe) (12)

(10a) (10b)

Similarly, Substituting the voltage equation, Eq.(9c) into active power balance, Eq.(6),we get

(11a) (13) (11b) Likewise, Equations (10) and (11) have also a straightforward solution and depend on terminal voltage magnitude, mechanical power input and the equivalent circuit parameters of the induction machine except the rotor slip. It is noted that for the induction generators the mechanical power is negative due to the change on the ow of power which is transferred from rotor (bus-r) to the stator (bus-s). Hence, the mechanical power input, Pm, must be taken in negative sign for both developed models. Since there are two solutions of real and reactive power outputs of the developed Model-II[11], the required root of the equation, given in Equations (10) and (11), must be determined. For this reason, variation of the active and reactive power output of the generator with different mechanical power input (Pm) are calculated using both expressions, for a hypothetical machine. On the other hand, the solution of Equations (10b) and (11b) remain on the feasible region. Therefore, the solution of Equations (10b) and (11b) is unique for the reactive and active power outputs of the induction generator, respectively. They can facilitate the computation of real and reactive powers of the induction generator for a specified mechanical power input and terminal voltage, in a simple way. Presently various types of WTGU have been installed and they can be broadly classified into three categories, namely fixed, semi-variable and variable speed types[7]. The models developed here for the WTGU are intended to obtain the
Issn 2250-3005(online)

Semi Variable Speed WTGU [7] This class of WTGU consists of a pitch controlled wind turbine and a wound rotor induction generator. The rotor circuit of the generator is connected to an external variable resistance. There are two controllers, a pitch controller and rotor resistance controller. These two controllers are designed to operate in a coordinated manner. This design guarantees that the active power output is equal to the maximum power at wind speeds below nominal and equal to rated power above nominal wind speeds. For this class of WTGU also, the manufacturer provides the designed real power out-put versus wind speed characteristics [7]. We note that the expression for Pe and Qe can be recast as a quadratic function of a new variableR2/s. Hence, even when R2 and s are unknown the quantityR2 /s (say Req) can be computed by solving the quadratic equation in Req involving Pe. To compute Qe, this value of Req is used in the modified expression for Qe. The quadratic equation for Req is as follows: (14) Where (15)

(16)

(17) Given wind speed u w and the terminal voltage V

September| 2012

Page 1523

International Journal Of Computational Engineering Research (ijceronline.com) Vo2. 2 Issue. 5

1. For the given uw obtain Pe from the power curve of the WTGU (provided by the manufacturer). 2. Compute 3. Knowing Req , compute Reactive power Qe as (19) (18)

III. LOAD FLOW ANALYSIS USING Wtgss

As the wind generator is of induction type the bus connected to wind generator can be modeled as PQ bus in load flow algorithm. The load flow algorithm including wind generator is explained through flow chart.
Initialize the bus voltages of the system, and mechanical power of the WTGSs for a specified wind speed by using Equation (1)

The Pitch regulated fixed speed WTGS model and Semi variable speed WTGS model considered here are the Vestas unit of 1MW rating .The Induction generator circuit parameters are given in Table I. TABLE I Induction generator circuit parameters Type of Pitch Regulated Semi Variable WTGS Fixed Speed Speed Model R1(pu) Xl1(pu) R2(pu) Xl2(pu) Xm (pu) Xc (pu) 0.005986 0.08212 0.01690 0.107225 2.5561 2.5561 0.005671 0.15250 0.00462 0.096618 2.8985 2.8985

Set the iteration counter (it=1) Calculate active and reactive power outputs of WTGSs by using: Eqs: (5) and (7) for Model-I Eqs: (10b) and (11b) for Model-II (12) and (13) for PRFS Calculate loads and WTGSs current as

Calculate active and reactive power transferred from each branch including losses and each line currents by using loads and WTGSs powers and currents, respectively. Calculate each line receiving end bus voltages using forward voltage formulation of particular load flow algorithm.

TABLE II Induction generator circuit parameters (Model I & II) 1. 2. 3. 4. Generator short circuit resistance(ohm) Short circuit reactance(ohm) Open circuit resistance (ohm) Open circuit reactance(ohm) 7.81 18.76 1230 500.6

The system reaches convergence?

NO
it = it+1

YES
Print Results.

5.

System rms base voltage(kV)

12.66

Stop

Fig. 2. Flow chart for distribution load flow with WTGS models IV. Results And Discussion A 33-bus distribution system is used to evaluate the effect of WTGS on the power losses, voltage prole of radial distribution system, and on the performance of load ow algorithms. Sweep-based algorithms are used and they all have been coded in

Issn 2250-3005(online)

September| 2012

Page 1524

International Journal Of Computational Engineering Research (ijceronline.com) Vo2. 2 Issue. 5

Matlab. Transformer parameters are taken from [7]. Tables III & IV represents the load ow results of the systems with and without WTGS modeled by using both developed models, pitch regulated and semi variable speed W T G S m o d e l s are also depicted. Similarly, in load ow analyses, initial voltage magnitude at all buses is considered to be the same as the source bus (1.0 pu) for test system and a tolerance of 10-4 on convergence criteria is used. M echanical power input is selected as Pm =1.0 MW. From the tables, it is clearly seen that the inclusion of a WTGS signicantly, alters the power intake from the main substation, and total power losses of the test systems are also decreased. Moreover, the voltage magnitude of the bus at which the WTGS i s connected is increased due to injected active power. The Induction machine circuit parameters and the data is taken from [11] for Model I & Model II, whereas for the remaining models i.e., Pitch regulated fixed speed and Semi variable speed is taken from [7].It is clearly observed that from the Tables III & IV developed models are more advantageous than the Pitch regulated fixed speed WTGS and Semi variable speed WTGS models. is the WTGS bus voltage. is the mechanical power input(1MW). is the WTGS active power generation. is the WTGS reactive power demand. is the active power loss. is the reactive power loss. Semi Variable Speed(SVS) WTGS. Pitch Regulated Fixed Speed(PRFS) WTGS. ABLE III Load flow result of the 33-bus system for Pitch Regulated and Semi Variable Speed WTGS models

Bus No. 6 18 25 32 6 18 25 32 6 18 25 32

WTGS Model _ _ _ _ SVS SVS SVS SVS PRFS PRFS PRFS PRFS

Ps (MW) _ _ _ _ 1 1 1 1 1 1 1 1

Qs (MVAr) _ _ _ _ 0.624 0.640 0.616 0.650 0.5665 0.5640

Pl (MW)

Ql (MVAr) 0.9331 0.8959

0.369

0.247

0.9648 0.8788 0.9614 0.9243

0.316

0.219

0.9737 0.9246 0.9593 0.9218

0.285 0.572 0.564

0.197 0.9731 0.9202

TABLE IV Load flow result of the 33-bus system with and without WTGS models.(Model I & Model II) Bus No. 6 18 25 32 6 18 25 32 WTGS Model _ _ _ _ II II I I Ps (MW) _ _ _ _ 0.919 0.919 0.992 0.992 Qs (MVAr) _ _ _ _ 0.923 0.923 0.017 0.017 0.181 0.129 0.369 0.247 Pl (MW) Ql (MVAr) 0.9331 0.8959 0.9648 0.8788 0.9680 0.9295 0.9750 0.9449

Issn 2250-3005(online)

September| 2012

Page 1525

International Journal Of Computational Engineering Research (ijceronline.com) Vo2. 2 Issue. 5

33-BUS RADIAL DISTRIBUTION SYSTEM

Voltage Magnitude

1 0.95 0.9 Without WTGS 0.85 0 5 10 With WTGS 25 30 35 15 20 Bus Number

the voltage profile is analyzed. The effectiveness of linking two buses can be studied in terms of decrement in losses and improvement in voltage profile. And the same load flow analysis can be done for other models i.e., Variable speed WTGS and Stall regulated fixed speed WTGS are also useful for the future planning of the system.

References
[1] A. Losi, M. Russo, "Object Oriented Load Flow for Radial and Weakly Meshed Distribution Networks," IEEE Trans. Power Syst., vol. 18, no.4, pp. 1265-1274, Nov. 2003. W.H. Kersting, Distribution System Modeling and Analysis, CRC Press ,New York,2002. D.Shirmohammadi, H.W.Hong, A.Semlyen, and G.X.Luo,A compensation based power flow for weakly meshed distribution and transmission networks, IEEE Trans. Power Syst., vol.3, pp. 753-762, may 1998. C. S. Cheng, D. Shirmohammadi, A three-phase power flow method for real-time distribution system analysis, IEEE Trans. Power Syst., vol10, pp. 671679, May 1995. Y. Zhu, K. Tomsovie, Adaptive Power Flow Method for Distribution Systems with Dispersed Generation, IEEE Trans. Power Del, vol. 17, pp. 822-827, July 2002. Feijoo A, Cidras J. Modeling of wind farms in the load ow analysis, IEEE Transaction on Power Systems 2000;15(1):1106. Divya KC, Rao PSN. Models for wind turbine generating systems and their application in load ow studies, Electric Power Systems Researc2006;76:844-56. Cutsem TV, Vournas C. Voltage stability of electric power systems, Boston: Kluwer Academic Publishers Springer; 1998. Hatziargyriou ND, Karakatsanis TS, Papadopoulos M. Probabilistic load ow in distribution systems containing dispersed wind power generation, IEEE Transaction on Power Systems 1993;8(1):15965. Haque MH. Load ow solution of distribution systems with voltage dependent load models, Electric Power Systems Research 1996;36(3):1516. Ulas Eminoglu Modeling and application of wind turbine generating system (WTGS) to distribution systems, Renewable Energy, vol.34,no.11, pp.24742483, Nov.2009. Satyanarayana S, Ramana T, Sivanagaraju S, Rao GK. An efficient load flow solution for radial distribution network including voltage dependent load models. Electric Power Components and Systems 2007;35(5):539-51.

Fig. 3 Voltage profile variation of the 33-bus system for PRFS WTGS model
33-BUS RADIAL DISTRIBUTION SYSTEM

[2] [3]

Voltage Magnitude

1 0.95 0.9 With WTGS 0.85 0 5 10 Without WTGS 15 20 Bus Number 25 30 35

[4]

[5] [6] [7]

Fig. 4 Voltage profile variation of the 33-bus system for SVS WTGS model

Voltage Magnitude

33-BUS RADIAL DISTRIBUTION SYSTEM 1 0.95 0.9

[8] [9]

Without WTGS
0.85 0 5 10 15

With WTGS
20 25 30 35

Bus Number
Fig. 5 Voltage profile variation of the 33-bus system (with and without WTGS models) V. CONCLUSION In this paper, two new models[11] for wind turbine generating systems are provided. The main advantage of the developed analytical models for WTGS is that they facilitate the computation of real and reactive power outputs for a specied mechanical power input and terminal voltages, in a simple way such that it does not require computation of any system parameters (i.e., rotor slip).Whereas the other proposed models in [7] rotor slip and power curves plays a crucial role. These developed models [11] can overcome the computational complexity and the computational time problem. The WTGS impact on the system losses and
[10] [11]

[12]

Issn 2250-3005(online)

September| 2012

Page 1526