Poonam Singhal et al Int. Journal of Engineering Research and Applications ISSN : 2248-9622, Vol. 4, Issue 7( Version 3), July 2014, pp.121-133 RESEARCH ARTICLE www.ijera.com OPEN ACCESS Transient Stability Enhancement of a Multi-Machine System using Particle Swarm Optimization based Unified Power Flow Controller Poonam Singhal1, S. K. Agarwal2, Narender Kumar3 Abstract In this paper an attempt has beenmade to investigate the transient stability enhancement of both SMIB and Multi-machine system using UPFC controller tuned by Particle Swarm Optimization. Power injection modelfor a series voltage source of UPFC has been implemented to replace UPFC by equivalent admittance. The admittance matrix of the power system is then modified according to the power injection model of UPFC. To mitigate the power oscillations in the system, the required amount of series voltage injected by UPFC controller has been computed in order to damp inter area & local mode of oscillations in multi-machine system. Keywords:UPFC, Transient Stability, Converters, SMIB, Multi-machine modeling, Capacitor Dynamics,UPFC(Unified Power Flow Controller), FACTS (Flexible AC Transmission System), Power injection model. Particle Swarm Optimization. I. INTRODUCTION Modern interconnected power system comprise a number of generators, a large electrical power network & a variety of loads& such a complex system is vulnerable to sudden faults, load changes & uncertain nature resulting in instability & so the designing of such a large interconnected power systems to ensure secure & stable operation is a complicated problem. The power transmission network is an important element of the power system which is responsible for causing instability, sudden voltage collapse during transient conditions. In recent years, technologically developments in Flexible AC Transmission Systems or FACTS devices [1] provided better performance in enhancing power system. The FACTS envisage the use of a solid state power converter technology to achieve fast & reliable control of power flow in transmission line. The greatest advantage of using a FACTS device in a power transmission network is to enhance the transient stability performance by controlling the real & reactive power flow during fault conditions. The FACTS devices are classified into two categories; the shunt type comprises of SVC [21]&STATCOM while the series comprises of TCSC & TCPST. UPFC is the most versatile belong to FACTS family which combines both shunt & series features. It is capable of instantaneous control of transmission line parameters [5].The UPFC can provide simultaneous control of all or selectively basic parameters of power system [6, 7, 20] (transmission voltage, line impedance and phase angle) and dynamic compensation of AC power system. In this paper UPFC is used to inject thevariable voltage in series with the line, whose magnitude can varies from 0 to www.ijera.com Maximumvalue (depends on rating) and phase angle can differ from 0 to 2π which modulates the line reactance to control power flow in transmission line. By regulating the power flow in a transmission line, the speed of oscillations of the generator can be damped effectively. II. UPFC DESCRIPTION The UPFC constitute of two voltage source converters linked with common DC link. One of the converters is connected in series (known as series converter or SSSC) with the line via injection transformer and provides the series voltage and other converter (called shunt converter or STATCOM) is connected in parallel via shunt transformer to inject the shunt current and also supply or absorb the real power demand by series converter at common dc link. The reactive power exchanged at the ac terminal of series transformer is generated internally by voltage source converter connected in series. The voltage source converter connected in shunt can also generate or absorb controllable reactive power, thus providing shunt compensation for the line independent of the reactive power exchange by the converter connected in series. The UPFC, thus, can be summarize as an ideal ac to ac power converter in which real power can flow freely in either direction between the ac terminal of the two converters and reactive power can be generated or absorbed by the two converters independently at their own ac terminals. Figure 1 shows the equivalent circuit of the UPFC in which the UPFC can be represented as a two port device with controllable voltage source Vse in series with line and controllable shunt current source Ish. The voltage across the dc capacitor is 121|P a g e Poonam Singhal et al Int. Journal of Engineering Research and Applications ISSN : 2248-9622, Vol. 4, Issue 7( Version 1), July 2014, pp. maintained constant because the UPFC as a whole does not generate or absorb any real power. In this paper, the UPFC is used to damp the power oscillations and thus improves the transient stability of both SMIB and multi-machine power system.UPFC parameters have been optimized using Bacterial Foraging in paper [16]. Particle Swarm Optimization technique, another heuristic technique based on bird and fish flock movement behavior is implemented for optimizing the parameters of UPFC efficiently and effectively in real time. In this paper, section II describes about the UPFC. Section III explains the brief description of the SMIB test systemtaken whereas section IV explains the mathematical modeling of the SMIB system equipped with UPFC Controller. Section V briefly describes the multi-machine power system under study. Section VI illustrates the modeling ofmulti-machine power system. Modeling of unified power flow controller for load flow is illustrated in Section VII. Section VIII illustrates the simulation results ofboth SMIB system and multi-machine system without and with UPFC controller based on PSO and section IX concludes and expresses the future scope of work. I1 I2 www.ijera.com Gen Bus1 Bus2 UPFC Fig 2. IV. MATHEMATICAL MODELING a. Synchronous Machine Model Mathematical models of synchronous machine vary from elementary classical model to more detailed one. Here, the synchronous generator is represented by third order machine model [2,4]. Different equations for stator, rotor, excitation system etc. are as follows: Stator: Vq Eq' rs iq xd' id Vd Ed' rs id xq' iq V1 vd jv q I id ji q S V1 I * P v d id v q iq Q v q i d v d iq Where, V2 V Ish 1 d dt Pm Pe dt M Vse Let D and K=0 rs is rotor winding resistance Fig. 1.Equivalent Circuit of UPFC III. SMIB SYSTEM UNDER STUDY For analysis of the UPFC for damping the power swings, parameters of the system like generator rotor angle, theta, terminal voltage, active power and reactive power are considered. A 200MVA, 13.8KV, 50Hz generator supplying power to an infinite bus through two transmission circuits as shown in fig(2) is considered. The network reactance shown in fig. is in p.u. at 100 MVA base. Resistances are assumed to be negligible. The system is analyzed with different initial operating condition, with quantities expressed in p.u. on 100MVA, 13.8 KV base. P=0.4p.u, 0.8p.u, 1.2p.u and 1.4p.u are different operating conditions considered. The other generator parameters in p.u. are given in appendix. The three phase fault at infinite bus bar is created for 100msec duration & simulation is carried out for 10sec. to examine the transient stability of the study system. www.ijera.com x d' is d-axis transient reactance xq' is q-axis transient reactance E d' is d-axis transient voltage Eq' is q-axis transient voltage Rotor: ' Tdo dEq' dEq' dt dt E q' E f xd xd' id E f Eq' xd xd' id ' Tdo Where, Tdo' is d-axis open-circuit transient time constant Tqo' is q-axis open-circuit transient time constant E f is field voltage Torque Equation Te Eq' iq Ed' id xq' xd' id iq Excitation System d ( E fd ) K e (Vref Vt ) E fd dt Te Te 122|P a g e Poonam Singhal et al Int. Journal of Engineering Research and Applications ISSN : 2248-9622, Vol. 4, Issue 7( Version 1), July 2014, pp. -0.6≤Efd≤ 0.6 Where Gen b. UPFC Controller Structure in PQ Mode P-Q demand on load side is met by controlling series voltage injection. In order to achieve the system stability, it is required to control the in phase & quadrature component of the series injected voltage after fault near infinite bus-bar. 1. Series Converter Let the voltage injected by the series converter is V L. In d-q frame of reference, VL can be written as VLd VL sin , VLq VL cos VL Series Transformer Shunt Transformer Shunt Inverter ∆QLine u 1, u 2 DC Link Series Inverter Control q ∑ P (+) ref u1=∆Vq, u2=∆Vp Fig.4 PI controller of the study system in PQ Mode Ө Fig. 3.d-q representation of series converter In phase and quadrature components of V L are responsible forreactive and active power flow in line. VLp VLpo V p VLq VLqo Vq X X i Vbd VLd X d' X 1 iq X 1isq ' d Vbd2 Vbq2 Pref Vd VLd Qref Vq Pline(-) d Vb Qref(+) ∆PLine Vt Vbq VLq e ∑ 2 Where, tan 1 vd vq ' q QLine(-) Series Inverter V Lq VLd 2 www.ijera.com id V Lq id 1 d X 1isd (6) i sd Vq VLq iq i sq i sd Vd V Ld iq i sq Where, i isd VLpo VLd sin tan 1 d iq isq i isd VLqo VLq sin tan 1 d iq isq T Vq k p1P ki1 P 0 P Pref P VLq cos tan 1 id iq VLd cos tan 1 id iq isd isq isd isq 2. Shunt Converter Let ish be the current injected by shunt voltage source converter which is in same phase as that of generator terminal voltage, hence it will not supply or absorb reactive power & its aim is to provide the real power demand of series power voltage source converter. Let lossless UPFC device is considered, then * R V1 I1* V2 I 2 0 And with losses, to maintain the voltage across the capacitor, shunt power should be equal to sum of series power and capacitor power. dVdc Vt ish VLd id isd VLq iq isq Vdc C dt From the above equation, ish can be obtained. Where, ishq ish cos ishd ish sin , q Vt Ө ish T V p k p 2 Q ki 2 Q 0 Q Qref Q d Fig.5. d-q representation of shunt converter c. Transmission Line currents The transmission line current it is split into d-q components represented as itd hand itq. itd id ishd www.ijera.com , itq iq ishq 123|P a g e Poonam Singhal et al Int. Journal of Engineering Research and Applications ISSN : 2248-9622, Vol. 4, Issue 7( Version 1), July 2014, pp. d. Capacitor Dynamics The difference between the shunt power and series power is the capacitor power. Mathematically, Capacitor power=Shunt power– Series power dVdc C Vdc Psh Pse dt Psh Vt ish Pse VLd itd VLq itq dVdc C dt Bus3 Bus7 Ln Bus5 Bus1 Bus6 ln ln Ln UPFC -2 3 G 3 ln ln YL3 G1 YL1 G2 Bus2 ln2 Vdc Vt ish VLp itd VLq ttq YL2 YL4 DESCRIPTION OF MULTIMACHINE SYSTEM A multi-machine power system (8 bus system) with twoUPFC‟s connected are considered as shown in fig 2. Three generators are connected at buses 1, 2 and 3.Generator 1 is hydro & G2 & G3 are Thermalgenerators.The parameters of all the generators are given in appendix. Bus number 1 is taken as slack bus. The transmission line parameters are also given in appendix.UPFC1 is connected between bus number 4 and 5, while UPFC2 is connected between bus number 7 and 8. Loads are connected at bus number 1, 2 and 3. They are represented in terms of admittances YL1, YL2, YL3 and YL4& are computed from load bus data. Using power Injection model[19] of UPFC1 and UPFC2, they are represented by equivalent admittances and the admittance matrix of the power network is then modified. Simple AVR are connected to each generator. The operating conditions taken are: p1=4.5, q1=1.5, p2=1.3, q2=0.6,p3= 1.0,q3=0.5. Fault is created at the Centre of line3-4at 0.5 sec & cleared at 0.6 sec which results in power deviation in the transmission line. The real power deviation in the transmission line generates the component VSEQ of series voltage injection in quadrature with line current and the reactive power deviation in line generates the in phase component VSEP. kp and ki relates these components to power deviations. These PI controller parameters are then tuned by PSO technique to damp the inter area and local area oscillations effectively. VI. Bus4 Bus8 UPFC -1 3 ish then can be computed from the above equation V. www.ijera.com MATHEMATICAL MODELING A. Synchronous Machine Model Mathematical models of synchronous machine vary from elementary classical model to more detailed one. Fig. 6 Here, the synchronous generator is represented by third order machine model. The differential equations governing the dynamics of each machine (1,2 and3) are same as for SMIB system modeling given above. . B. Power Injection Model of UPFC Two voltage source model of UPFC shown in fig. (1) is converted into two power injections in rectangular form for power flow studies [3, 27]. The advantage is that it maintains the symmetric characteristics of the admittance matrix. Vi Pij j i Zi0 Ii0 Pdc Qij c Fig. V SH Fig. 7 1. Shunt Converter: Let the current injected by the shunt converter is I i 0 & voltage source is VSH Shown in fig.(7). The shunt side of UPFC is converted into power injection at bus bar i only. * V V Si 0 Pio jQio Vi i SH Zio Where www.ijera.com Vse Z io Rio jX io 124|P a g e Poonam Singhal et al Int. Journal of Engineering Research and Applications ISSN : 2248-9622, Vol. 4, Issue 7( Version 1), July 2014, pp. www.ijera.com Pi inj Bij Vi sin and Qi inj Bij Vi cos Similarly, 2 2 Vse Vi Zij Pij Vj q i j Pi0 Qi0 S j V j I SE V j VSE Bij * * S j V j jBij Vi e j Bij Vi V j sin ij jBij Vi V j cos ij * Pj inj Bij Vi V j sin ij Qij Q j inj Bij Vi V j cosij Fig.8 Where Since the shunt reactive compensation capability of UPFC is not utilized that is the UPFC shunt converter is assumed to be operating at unity power factor [20].Its main function is to transfer the real power demand of series converter through the dc link, so Pio Vi I SH and Qio 0 ij i j Based on the explanation above, the injection model of a series connected voltage source can be represented by two independent loads [3] as shown in fig (27). Vi 1. Series Converter: Let the ideal voltage injected by the series converter is Vse and reactance Xij be present between two buses (i,j) in the power system shown below[26]. The series side of UPFC is then converted into two power injections at buses i& j. Vi Vj Xij Vse Vj Xij Pj Bij Vi V j sin ij Pi Bij Vi sin 2 Q j Bij Vi V j cos ij Qi Bij Vi cos 2 Fig. 11 Further the real power associated with converter-1 can be written as Pio Vi I SH Where I SH is in phase current with the bus voltage Fig. 9 Vi The Norton equivalent of the above circuit is shown as: Xij Vi Vj Pdc is the power transfer from shunt side to series side. Vi V SE V j Pdc Re V SE X ij * V V se V j Re V e j i i X ij On solving, ISE Pdc Bij Vi V j sin ij Bij Vi sin 2 Fig.10 I SE When power loss inside the UPFC is neglected, than Pio=Pdcand modified injection model is formulated as shown in fig.(11) VSE in parallel with the line Z ij Where, Z ij 0 jX ij , V S i Pi jQi Vi SE X ij Since I VSEQ tan 1 ijd I V SEP ijq tan 1 * www.ijera.com tan 1 Vid V iq j Pj Bij Vi sin jBij VS cos 2 Vj i * Vi VSE Bij and VSE Vi e j Si Vi jBij Vi e j Xij Vi Bij 1 / X ij 2 Pi0 Qi0 Qj Pi Qi 125|P a g e * Poonam Singhal et al Int. Journal of Engineering Research and Applications ISSN : 2248-9622, Vol. 4, Issue 7( Version 1), July 2014, pp. Fig.12 & the reference frame of the ith machine are related Si inj Si Sio Pi Pio jQi by the transformation [4]. S j inj S j Pj jQ j Vi Xij Xse VDi cos i sin i Vdi V Qi sin i cos i Vqi Vj Pj Bij Vi V j sin ij Pi Bij Vi sin Vi I SH 2 Q j Bij Vi V j cos ij Qi Bij Vi cos 2 Fig. 13 I ijd and I ijq are d-q axis transmission line currents. Vid The component in phase and quadrature component of Vseis responsible for real and reactive power flow in line & hence finally mitigate the power oscillations. P1 Pref 1 P1 I Di cos i sin i I di I Qi sin i cos i I qi Where δi is the angle between the machine reference axes (di, qi) and the common reference axes (D,Q) as shown: and Viq are d-q axis voltage of ith bus of UPFC. di D VD i Q1 Qref 1 Q1 , P2 Pref 2 P2 Vi Q2 Qref 2 Q2 , For Machine 1, VSEP1 kp1Q1 ki1 Q VSEQ1 kp2P1 ki2 P1 1, For Machine 2, VSEP2 kpp1Q2 kii1 Q2 , Vqi qi Vd i δi 2. Dynamics of Capacitor Fig. 14 The dynamics of the D.C voltage neglecting losses can be represented by dVdc Pio Pdc dt dVdc 1 2 Vi I SH Bij Vi V j sin ij Bij Vi sin dt C Vdc Above equation will be written for both UPFC‟s. VII. PROCEDURE FOR MULTIMACHINE POWER SYSTEM SIMULATION The digital simulation reads the initial nodes specifications and generates the steady state load flow solution assuming the reference bus voltage in p.u as 1∟0 , corresponding to a common reference 0 rotating representing Further the voltage behind the transient reactance (Eq‟), transient reactance (xdi‟), quadrature axis reactance (xqi) of the ith machine are related to its By putting the values of Pio and Pdc, in above equation,the dynamics of the D.C link is represented as: frame Q VQi VSEQ2 kpp2P2 kii2 P2 CVdc www.ijera.com at synchronous speed thus terminal voltage components in the common reference frame as 0 0 xqi cos i sin i I Di cos i sin i VDi ' ' Eqi xdi 0 sin i cos i I Qi sin i cos i VQi Hence machine angle δi can be computed as x qi I Qi V Di tan i x qi I Di VQi This voltage is then utilized in the 3rd order machine model and the differential equations are solved. Using this machine angle, the voltage behind the transient reactance of machine are solved. Moreover for an „n‟ number of machines & „m‟ number of load buses in a power network, the algebraic equations are written in compact form: 1 1 E m x mTmn I Tmn V I YmV 1 1 E m x mTmn YmV Tmn V Q axis. The common reference frame www.ijera.com 126|P a g e Poonam Singhal et al Int. Journal of Engineering Research and Applications ISSN : 2248-9622, Vol. 4, Issue 7( Version 1), July 2014, pp. Hence the generator terminal voltage components in the common reference frame and the machine internal voltages are related by 1 1 V U Tmn xmTmn Ym Tmn Em www.ijera.com A. Algorithm Step 1: Initialize the number of particles,minimum and maximum limits for the d-dimension search space mm ax , mm in anditmax. Step2: Set initial position of particles x ik Where, xik m m in m m ax m m in * rand Where U is the unity matrix. Both UPFC‟s are represented by an equivalent Step 3: Set initial velocity vi of particles admittances derived from their power injection Where, models : vm in 0.1 * mm ax mm in YLi YLj Pi jQ j Vi 2 Pj jQ j V j2 Where Pi Bij Vi sin Vi I SH 2 Qi Bij Vi cos 2 Pj Bij Vi V j sin ij Q j Bij Vi V j cos ij When fault is created in any line, the admittance matrix is then modified according to position of fault& taking into account of UPFC admittances. All system equations are converted into common frame of reference which is rotating at synchronous speed. VIII. PARTICLE SWARM OPTIMIZATION BASICS PSO was introduced by Eberhart and kennedy in1995 [22]. It is a heuristic & stochastic based optimization technique. PSO can be used on optimization problem that are partially irregular, noisy, change over time, etc. It is developed from swarm intelligence and is based on the research of bird and fish flock movement behavior .The particle swarm optimization consists of swarm of particles which are initialized with a population of random candidate solution in the multidimensional search space. During their flying movement they follow the trajectory according to their own best flying experience (pbest) & best flying experience of the group (gbest). During this process, each particle modify its position & velocity according to shared information to follow the best trajectory leads to optimum solution & this technique is simple & very few parameters need to be determined. The choice of PSO parameters can have a large impact on optimization performance. Selecting PSO parameters that yield good performance has therefore been the subject of much research. www.ijera.com vidk vm in vm ax vm in * rand vm ax 0.1* mm ax mm in Step 4: Select initial pibestk, gbestk Step5: calculate cost function f (xi) taking initial position into consideration. Step 6: If f (xi) > f(pibestk) Then Update the value of xik by pibestk Else retain the pibestk as pibestk end if; Step 7: If f(pibestk) > f(gbestk) Then Update value of pibestk→gbestk Else go to step 8. step 8: If maximum number of iterations has been done then storegbestk as best position value. Else Increase order of iteration by one. Update particle velocity &position: vik 1 vik c1* rand1 pbestik xik c 2 * rand 2gbest k xidk xik 1 xik vik 1 Where, c1 and c 2 represents the acceleration factors, rand1 and rand 2 represents distributed random numbers between (0,1).First part of equation (1) depicts the previous velocity of the particle, the second part is a positive cognitive component & third part is a positive social component as described in [23]. Repeat Step 4 to step 7 Else if, End. IX. FORMULATION OF AN OBJECTIVE FUNCTION Objective function formulated is based on the optimization parameters. It is worth mentioning that the PIparameters of UPFC are tuned using PSO to minimize the power system oscillations after a disturbance so as to improve the transient stability. These oscillations are reflected in 127|P a g e Poonam Singhal et al Int. Journal of Engineering Research and Applications ISSN : 2248-9622, Vol. 4, Issue 7( Version 1), July 2014, pp. the deviations in the generator rotor speed ∆w1,∆w2 and∆w3 in the present study. The objective function J for SMIB system is formulated as the minimization of J t t , x dt where t 2 t, x denotes the 0 rotor speed deviation. And for multi-machine system, t J t 1 t , x 2 t , x 3 t , x dt 2 2 2 0 In the above equations, ∆w1 (t,x), ∆w2 (t,x) and ∆w3 (t,x) denotes the rotor speed deviations of generator 1,2& 3 for a set of controller parameters x and (note that, here x represents the parameters to be optimized i.ekp&ki, the parameters of PI controller) and t is the time range of the simulation. For objective function calculation, the time-domain simulation of the power system model is carried out for the simulation period. It is aimed to minimize this objective function in order to improve the system response in terms of the settling time and overshoots. Flow chart of PSO X. www.ijera.com Simulation The SMIB system with UPFC is simulated using PSO technique in MATLABenvironment. The number of parameters to be optimized are four (kp1, kp2 , ki1& ki2). The analysis is carried out with 3 phase fault at Infinite bus for 0.1 sec& system is run for 10 sec. The result have been shown with controller parameters tuned by PSO (each particle is assigned with a set of 4 variables to be optimized and assigned with random value within the universe of disclosure). The following are the simulation results of the system without controller and with PSO based UPFC controller for different operating conditions. Case 1: P=0.4p.u, Q=0.6p.u START Initialize the number of Particle n, Particle position xi (mmax,mmin) and Particle velocity Vi (Vimax, Vimin) and max.iter Set Initial Velocity vik, pbestk,gbestk→ Calculate function f(xi) If f(xi)>f(pibestk) NO YES Pibestk = pibestk Update xik→pibestk If f(xi)>f(gbestk) Update Particle Position and velocity YES Updatexik→gbestk NO Count = count+1 Is Count= itmax YES gbestk is the best position value STOP www.ijera.com 128|P a g e Poonam Singhal et al Int. Journal of Engineering Research and Applications ISSN : 2248-9622, Vol. 4, Issue 7( Version 1), July 2014, pp. www.ijera.com 2.6 2.4 2.2 delta(rad) 2 1.8 1.6 1.4 1.2 1 0.8 0 1 2 3 4 5 6 7 8 9 10 time(sec) Fig.18 (i) 6 Speed deviation(rad/sec) 4 Case 2: P=1.2p.u, Q=0.6p.u 2 0 -2 -4 -6 0 1 2 3 4 5 6 7 8 9 10 time(sec) Fig. 18(ii) In case of multi-machine power system simulation, bus 1 is taken as a slack bus.A fault is created at the centre of the transmission line connecting bus numbers 3 and 7 & simulation is carried out without and with PSO based UPFC Controller foroperating conditions P1=4.5p.u, Q1=1.5p.u, P2=1.3p.u, Q2=0.6p.u, P3=1.0p.u, Q3=0.5p.u. Both inter area & local mode of oscillations are shown below. Fig. 16(ii) Case 3: P=1.4p.u, Q=0.6p.u, PSO based UPFC controller provides the power system stability as shown in fig.18 while without UPFC, system becomes unstable as shown in fig.17 300 delta(rad/sec) 250 200 150 100 50 0 0 1 2 3 4 5 6 7 8 9 10 time(sec) Fig.17 www.ijera.com 129|P a g e Poonam Singhal et al Int. Journal of Engineering Research and Applications ISSN : 2248-9622, Vol. 4, Issue 7( Version 1), July 2014, pp. www.ijera.com 0.3 delta3-delta2(rad) 0.2 0.1 0 -0.1 -0.2 -0.3 -0.4 0 5 10 15 20 25 30 35 40 45 time(sec) Fig. 19 (vi) z 0.05 delta2-delta1(rad) 0 -0.05 -0.1 -0.15 -0.2 -0.25 -0.3 -0.35 -0.4 0 5 10 15 20 25 30 35 40 45 50 30 35 40 45 50 time(sec) Fig. 19 (iv) 0 delta3-delta1(rad) -0.05 -0.1 -0.15 -0.2 -0.25 -0.3 -0.35 -0.4 -0.45 0 5 10 15 20 25 time(sec) Fig. 19 (v) www.ijera.com 130|P a g e 50 Poonam Singhal et al Int. Journal of Engineering Research and Applications ISSN : 2248-9622, Vol. 4, Issue 7( Version 1), July 2014, pp. www.ijera.com speed deviation, machine angle, generator terminal voltage, active power and reactive power are settled down quickly after the disturbance in an SMIB system and angular speed deviation, rotor angle, & voltage oscillations for inter area (m/c 2 vs. m/c 1 and m/c 3 vs. m/c1) and local area (m/c 3 vs. m/c 2) are damped out effectively during fault in line 3-7 in multimachine power system. This paper has a vital scope for future studies. Artificial Intelligence Methods like genetic algorithm, Fuzzy Logic technique, ant colony optimization technique and Hybrid BFO-Particle Swarm optimization technique can be used to tune the UPFC controller parameters for better transient stability enhancement. Moreover additional auxiliary signals can be added as a supplementary signals for better transient stability enhancement. APPENDIX Generator Data(All reactance are on 100 MVA base) for SMIB system: Generator Data ' xd' 0.17p.u, xq' 0.3p.u, x do 0.15 p.u, xd 1.9 p.u, xq 1.6,Tdo‟ = 6.314 p.u, H= 5s M/c1 0.05 VSEP1(p.u0 0.04 0.03 0.02 0.01 0 -0.01 0 1 2 3 4 5 6 7 8 9 10 time(sec) M/c1 0 VSEP2(p.u) -0.01 -0.02 -0.03 -0.04 -0.05 -0.06 0 1 2 3 4 5 6 7 8 9 10 time(sec) XI. CONCLUSION The UPFC controller parameters (kp1, ki1, kp2, ki2) of UPFC Controller in an SMIB system and (kp1, kp2, kpp1, kpp2 ki1, ki2, kii1, kii2) in multi-machine system are tuned effectively for transient stability enhancement. The response curve shows that the www.ijera.com Exciter data: KE=50, TE= 0.1s Generator Data (All reactance are on 100 MVA base) for Multi-machine system: Generator 1 xd1 = 0.1468p.u , xd1‟ = 0.0608p.u , xq1= 0.0969p.u ,xq‟=0.0969,Tdo1‟ = 8.96, H1=5 Generator 2 xd1 = 0.8958p.u , xd1‟ = 0.1198p.u , xq1= 0.8645p.u,xq‟=0.1969, Tdo1‟ = 0.1969, H2=4 Generator 3 xd1 =1.3125p.u, xd1‟ = 0.1813p.u, xq1=1.2578 p.u,xq‟=0.25, Tdo1‟ = 0.25, H3=4 Exciter Data: KE1= 25,TE1= 0.05sec,KE2=25, TE2=0.02sec, KE3=25, TE3= 0.06sec Transmission line Data: z24=j0.068, z37=j0.068, z84=jz37, z78=0.1*z37, z45=j0.07, z56=0.1*z45, z61=j0.07 Equivalent Load Reactance: YL5=0, YL6=0, YL7=0, YL8=0, YL1=6.261-j1.044, YL2=0.0877-j0.029, YL3=0.0877-j0.0292, YL4=0.9690-j0.3391 UPFC Data: C = 3500 µF, BSE=142.85siemens, V dc base=31.5 kV PSO based Controller Data for SMIB system: kp1=0.3145, ki1=2.7387, kp2=0.5297, ki2=1.9115 PSO based Controller Data for Multi-machine system: kp1 = 0.045, ki1 =0.3245, kpp1 =0.0511, kii1 =0.2211 kp2=0.015,ki2=0.1021, kpp2=0.0115, kii2=0.3101 PSO Parameters: 131|P a g e Poonam Singhal et al Int. Journal of Engineering Research and Applications ISSN : 2248-9622, Vol. 4, Issue 7( Version 1), July 2014, pp. c1=0.2,c2=0.4, Number of Particles=30,wmax=0.5, wmin=-0.5 [14] References [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] N. G. Hingorani, L. Gyugyi, “Understanding FACTS”, IEEE Press, 2001. P. Kundur, “Power System Stability and Control”, Tata McGraw-Hill, 2006. Y. H. Song, A. T. Johns, “Flexible AC Transmission Systems (FACTS)”, IET, 2009. Power System control & Stability by P.M Anderson & A.A. Foaud. Ali Ajami, S H Hosseini, G B Gharehpetian, “Modelling and Controlling of UPFC for power system transient studies”, ECTI Transations on Electrical Eng., Electronics, And Communications Vol.5, No.2, August 2007. L.Gyugi, C.D. Schauder, “Unified Power Flow Controller: A New Approach To Power Transmission Control:, IEEE Transactions on Power Delivery, Vol.10, No.2, pp 1085-1093, April 1995. H.Fujita, Y.Watanabe, H.Akagi,”Control and Analysis of Unified Power Flow Controller”, IEEE Transaction on Power Electronics, Vol.14, No.6, pp. 1021-1027, November 1999. Pranesh Rao. M. L. Crow, Zhiping Yang, “STATCOM Control for Power System Voltage Control Applications”, IEEE Transaction on Power Delivery, Vol 15, No.4, pp 1311-1317, October 2000 S.H. Hosseini, A. Ajami, “Transient StabilityEnhancement of AC Trasmission System UsingSTATCOM," TENCON'02, October 2002BeijingChina L. Gyugyi, C.D. Schauder, K.K. Sen, “Static Synchronous Series Compensator: A SolidState Approach to the Series Compensation of Trasmission Line”, IEEE Transactions on on Power Delivery, Vol. 12, No. 1, pp. 406413,1997 S.H. Hosseini, A. Ajami, “Dynamic Compensation of Power System Using a Mulilevel Sataic Synchronous Series Compensator”, 8th International Iranian Conference on Electrical Engineering May 2000 A. Nabavi-Niaki, M. R. Iravani, “Steady State and Dynamic Models of Unified Power Flow Controller for Power System Studies,” IEEE Transactions on Power System, Vol. 11, No. 4, November 1996, pp. 1937-1943 C.Shauder and H.Mehta, “Vector analysis and control of advanced static VAR compensators”, IEE Proc, 140, No. 4, July 1993 www.ijera.com [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] www.ijera.com R. Mihalic et. al., “Improvement of Transient Stability Using Unified power Flow Controller,” IEEE Transactions on Power Delivery, Vol. 11, No. 1, January 1996, pp. 485-492 S. Limyingcharoen, U.D. Annakkage and N.C.Pahalawaththa, “Effects of Unified Power FlowController on Transient Stability,” IEE Proc.Gener. Transm. Distrib., Vol. 145, No. 2, March 1998, pp. 182-188 Poonam Singhal, S. K. Agarwaland Narender Kumar,”Optimization of UPFCController ParametersUsing BacterialForaging Technique forEnhancing Power System Stability”,Vol.5.No.2 (March2014),IJOAT. Mahmoud H. M, M. A. Mehanna and S. K. Elsayed, “Optimal Location of Facts Devices to Enhance the Voltage Stability and Power Transfer Capability”, Journal of American Science, 2012. Theja,B.S. ; Rajasekhar, A. ; Kothari, D.P., “An intelligent coordinated design of UPFC based power system stabilizer for dynamic stabilityenhancement of SMIB power system” IEEE International Conference on Power Electronics, Drives and Energy Systems (PEDES), 2012 Flexible AC Transmission Systems(FACTS), by Yong Hua Song & Allan T.Johns. Poonam Singhal,S.K.Agarwal & Narender Kumar,”Transient Analysis Enhancement of a Power System using Unified Power flow controller,IREMOS, Volume 6,no5,Dec2013. Poonam Singhal,S.K.Agarwal & NarenderKumar,”StabilityEnhancement of a Multi-Machine Power System using Static VAR Compensator.IJAIR, Vol2,Issue2,2013. J. Kennedy and R. Eberhart, “Particle swarm optimization,” in Proc.IEEE Int. Conf. Neural Networks, 1995, pp. 1942– 1948., Asanga Ratanweera, Saman k. Halgamuge, Harry C. Watson,”Self Organising Hierarchical Particle Swarm Optimizater with time-varying Acceleration Coefficient Sedraoui, K.; Fnaiech, F.; Al-Haddad, K., “Solutions for Power Quality Using the Advanced UPFC - Control Strategy and Case Study”IREE Volume: 3 , Issue: 5, pp. 811-819, Sept-Oct, 2008 Iraj Rahimi Pordanjani, Hooman Efranian Mazin, Gevorg B. Gharehpetian, “Effects of STATCOM, TCSC, SSSC and UPFC on Static Voltage Stability”IREE Volume: 4 , Issue: 4, part B, pp. 1376-1384, Nov-Dec, 200f9 132|P a g e Poonam Singhal et al Int. Journal of Engineering Research and Applications ISSN : 2248-9622, Vol. 4, Issue 7( Version 1), July 2014, pp. [26] [27] www.ijera.com A.nabavi Niaki and M.R.Iravani “SteadyState and Dynamic Models of UPFC for Power System Studies”IEEE Transactions on power systems,vol 11,.no.4,Nov1996 N.Dizdarevic and G.Andersson“Power Flow regulation by use of UPFC’s Injection model” Authors’ information 1 Associate Professor, Deptt. of Electrical Engineering, YMCAUST, Faridabad, India. 2 Professor and Head, Deptt.of Electronics Engineering, YMCAUST, Faridabad, India. 3 Professor, Deptt. of Electrical Engineering, DTU, Delhi, India. 4 Assistant Professor, Deptt. of EEE, MAIT, Delhi, India. Poonam Singhalwas born in 1964.She is currently Associate Professor in YMCA UniversityofScience& Technology, Faridabad, Haryana (India).She did her B.ScEngg from REC Rourkela,Odisha (India)&M.Tech (Power System & Drives) from YMCA University of Science &Technology, Faridabad, Haryana (India). Her main area of interest are FACTS, Power System Operation control & Simulation. S. K Agarwalwas born in 1961. He is currently Professor in Department of Electronics Engineering, YMCA University of Science &Technology, Faridabad and Haryana (India) .He did his B.Tech from REC Calicut,Kerala(India) & M.E. (Controls & Instrumentation) from Delhi Technological University Delhi(India) &Ph.D from JamiaMiliaIslamia University Delhi(India). Prof. S. K. Agarwal has many publications in National/International Journals & presented papers in International conferences. Prof.Agarwal‟s main area of interest are Power System & Controls. Narender Kumarwas born in Aligarh(India). He received his B.Sc.Engg.& M.Sc. Engg in Electrical Engineering from A.M.U. Aligarh in 1984 & 1986 respectively. Presently, he is Professor of Electrical Engg. in Delhi Technological University, Delhi. Prof.Narender Kumar has published various papers in national and international conferences and journals. Prof.NarenderKumar‟s main area of interest are Power System Operation Control &Simulation,SSR,and AGC etc. www.ijera.com 133|P a g e

© Copyright 2017 ExploreDoc