Investigating Possible Induction Generator Effects Due to Sub-Snchronous Resonances APPLICATION NOTES This application note deals with an investigation of possible induction generator effect triggered b sub-snchronous resonance frequencies during transient events. The following will be addressed in this application note: Harmonic impedance profiles of a power sstem network. Series compensation of lines and the resulting changes in the harmonic impedance profile. Induction generator effect phenomena. Machine model to stud induction generator effect phenomena. Voltage amplification problems due to induction generator effect. Sensitivit of sstem and machine parameters. The single line view of the PSCAD model of this network is shown below. A simple single machine, infinite bus tpe sstem is selected for simplicit, but the techniques discussed are tpical and applicable in a tpical investigation. Field Voltage held constant S2M S / H out in hold vfld cfld Initial Angle of source = -14.9 deg. Initial Ramp up time = 0.2 sec. Machine enabled at = 1.1 sec. - 100 [Hz] Z(f) Ef0 Ef If 376.99 telec Tmstd Te Tm S w Tm A V VMac #2 #1 VH IH 6.0832 [ohm].4323 [H] 21.977 [uf] VC VI 5178 002 TIME S2M Fault applied at 1.5sec. Duration 75 msec..0345 [H] 483 [ohm] Initial Angle of source = -3.813 deg. Initial Ramp up time = 0.2 sec. Source Magnitude = 477.8 kv. Timed Fault Logic ISW ABC->G Figure 1 Single line representation of the model used for the stud 500 kv sstem bus 002 5178 Initial Angle of source = -3.813 deg. Initial Ramp up time = 0.2 sec. Source Magnitude = 477.8 kv. Figure 2 Power sstem representation. The main power network is modeled as an impedance behind a voltage source as shown in Figure 2. The source impedance can be determined from the short circuit level at this bus. If the stud required a more accurate representation of the network frequenc response, more buses behind the 500 kv sstem bus must be added to the model. This is addressed in a separate application note. The transmission line from the 500 lv bus to the generating station is represented b simple R,L elements representative of the fundamental frequenc data. Once again, a detailed, fre-
6.0832 [ohm].4323 [H] 21.977 [uf] VC Figure 3 Simple representation of a series compensated line. quenc dependent line model ma be used in a practical stud but, for the purpose of this application note, this simple representation easil ields to verifing the sensitivit of line loss etc. to sub-snchronous effects. The series capacitor represents a series compensation. The value used here is representative of approximatel 75% compensation. The snchronous machine is modeled with all of its detailed parameters in place. The rotor circuit plas an important role during sub-snchronous events and thus, it is important to investigate the sensitivit of such parameters during the stud. Figure 4 Options available in the machine model for initialization. The simulation should be first initialized to a specific load flow situation. This can be achieved b entering the bus voltages and magnitudes at appropriate locations. In this example, the bus voltage information at the 500 kv source and at the machine terminal ma be used to initialize the simulation. During the initialization process, the machine can be made to act as a voltage source, operating at the specified magnitude and phase angle. This is achieved using the source to machine conversion feature of the machine model as shown in Figure 3. Initiall, the signal S2M from the timer is zero and the machine model will act as a voltage source during this period. When this value is changed to 1 at a specific time (when all initial transients have settled), the model will act as a machine, governed b the equations relating terminal voltages to the winding currents. A constant field voltage input is assumed for this stud. The constant field voltage is based on the initialized value (Ef0) computed b the machine model during the initialization period. This fixed value will ensure the same stead state operation after the machine model is switched from a source to a machine. The machine mechanical dnamics are not modeled for this stud. This can be realized in a number of was. Figure 5 Defining the machine speed through an external signal. 1. Run the machine in the Lock rotor mode b making this input entro equal to zero during the simulation. This is shown in Figure 4 and the machine speed with be equal to 1 pu (ie. snchronous speed). 2. An arbitrar speed can be specified b enabling the multimass option as shown in Figure 5. Figure 6 shows the parameter box in the machine model to enter the bus voltage angle and its magnitude for the specific power flow condition. Once the sstem is set up and initialized as outlined below, it will run to the specified stead state as can be verified b the results shown in Figure 7. Sstem harmonic impedance profile A tpical harmonic impedance profile of a high voltage network is shown in Figure 8. Figure 6 Initial terminal conditions at the machine bus. Tpicall, the resonance points are at super-snchronous frequencies (ie. higher than 60 Hz in a 60 Hz sstem). During
Figure 7 Stead state operation of the sstem. Figure 8 Profile at a 230 kv bus on the Northern California sstem. transient events, such as breaker operations, faults and fault clearance, the voltage and current waveforms would displa such frequencies. Since there are no driving forces (generators) to sustain such frequencies, the will be eventuall damped out at a rate determined b the sstem losses and loads. Figure 9 shows high frequenc transients during a breaker operation for the simple sstem shown in the same figure. It also shows the effect of losses on transient damping. 0.80 0.60 0.40 0.20 0-0.20-0.40-0.60-0.80 0.150 0.100 50 00-50 -0.100-0.150-0.200 Ea1 Ia1 Main : Graphs Ia2 0.180 0.200 0.220 0.240 0.260 0.280 0.300 0.80 0.60 0.40 0.20 0-0.20-0.40-0.60-0.80 0.150 0.100 50 00-50 -0.100-0.150-0.200 Figure 9 Transient response of a simple RLC circuit. Ea1 Ia1 Main : Graphs Ia2 0.493 0.214-0.279 Min -0.305 Max 0.496 67 01-66 Min -22 Max 67 0.1960 0.2000 0.2040 0.2080 0.2120 0 2030 0 2050 501 08 R=0 BRK Timed Breaker Logic Open@t0 503 Hz Ia2 0.1 Ea2 5 BRK Losses increase damping This sstem had a harmonic resonance at around 503 Hz and can be measured using the PSCAD model shown in Figure 10. This can be used to plot impedance profiles at different locations as shown in Figure 7. Series compensation of lines and the resulting changes in the harmonic impedance profile The addition of new equipment to the existing sstem will naturall effect the harmonic impedance profile. Sub-snchronous effects are of concern if the harmonic resonance points gets shifted to frequencies lower than the rated sstem frequenc. It is well known that the addition of series capacitors to compensate the transmission line reactance can give rise to this situation. - 100 [Hz] Z(f) VI 500 kv sstem bus 5178 002 Figure 10 Measuring the harmonic impedance in a network model in PSCAD.
Z - uncompensated Impedance profile Z-compensated Figure 11 Impedance profile with and without series compensation. Figure11 shows the impedance profile at the 500 kv bus for the sstem shown in Figure 1. The series compensation has resulted in a resonance point at around 40 Hz. Following a disturbance, the currents and voltages around this point will show slow transients at around 40 Hz. Figure12 shows such a sub-harmonic response in a sstem which had a resonance point at around 9 Hz. Subsnchronous currents and voltages in the network following a disturbance. 9.67 Hz Figure 12 Waveforms showing sub harmonic currents and voltages in a sstem that had a resonance point around 10 Hz. Induction generator effect phenomena The interaction of the sub harmonic currents and the voltage with the machine can result in the Induction generator Effect. The sub-harmonic currents will produce a rotating mmf which will assume a frequenc corresponding to the same frequenc. The rotor circuit, which is rotating at or near the rated snchronous speed, responds to the sub-harmonic mmf in a manner similar to an induction machine. Since the machine speed is greater than the sub-harmonic mmf rotation, the effect is similar to an induction machine in a generating mode where the slip is negative. This can be understood Articles and submissions addressing the use of PSCAD in the real world are alwas welcome. 2006 Manitoba HVDC Research Centre Inc
b examining the basic induction machine theor and the resulting stead state equivalent circuits as shown in Figure 13. Is 0.1 0.1 0.1 (1-s)/sRr Ws Wr S : = Ws R r s Effective rotor resistance Negative for generator action Figure 13 Equivalent circuit of a tpical induction machine. If the combination of the stator resistance and the resistance presented b the network, as seen at the machine terminals, is low enough, the effective resistance can be negative. This is the condition for positive feedback effect, refered to as the induction generator effect. Figure 14 Data entered in the PSCAD machine model. Machine model to stud induction generator effect phenomena The detailed machine model in PSCAD is suitable to stud possible induction generator effects. The data required for the machine model is shown in Figure14. VI 5178 002 The mechanical dnamics are not included in the example and the field winding voltage is assumed to be constant over the stud interval. However, these details can be easil included in the model. Voltage amplification problems due to induction generator effect A three phase fault is applied at the 500 kv bus at 1.5 s. The fault duration is 75 ms. The fault component and the timed fault logic component in the Master Librar are used to simulate the fault. This model is shown in Figure 15. Figure 16 shows different current and voltage waveforms upon the clearance of the fault. The transients die out and the sstem reaches a stead state. Results when the line resistance is lowered are shown in Figures 17 and 18. With the line resistance set to 2.0832 Ohms, the transients are sustained for a longer period. With the line resistance set to 832 Ohms, the terminal voltage grows to ver high values and this is a result of the induction generator effect. Anasis of the current waveform using the FFT component of PSCAD (see the PSCAD case) shows a prominent 40 Hz subharmonic component. Timed Fault Logic ISW.0345 [H] 3 2 2 1 1 3.0 2.0 - -2.0-3.0-483 [ohm] ABC->G Figure 15 PSCAD models used to simulate a fault. 2.0 3.0 Figure 16 Simulation results following a three phase fault clearance at 1.5 s. The fault duration is 75 s.
40 Hz 2 Vmac 3 2 2 1 1 [8] 5.40995773019-0 -0 2.0 3.0 Figure 17 Results when the line resistance is 2.0832 Ohms. 4 4 3 3 2 2 1 1 Figure 19 Harmonics in the current waveform. The rotor circuit interaction with the sub-harmonic stator current results in the induction generator effect. Thus, it ma be advisable to investigate the sensitivit of certain machine parameters to the simulation results. Figure 20 shows the results where the line resistance is as in the results in Figure 17 (i.e. Line R=2.0832 Ohms). The field time constant Tdo was lowered from 4.3 to 1.3 (just to see the effect). This results in a larger rotor circuit resistance and correspondingl, a negative equivalent resistance of a larger magnitude. Figure 20 shows that the sstem displaes a negativel damped induction generator effect. 3 3 2 2 1 1 8.0 6.0 2.0-2.0 - -6.0-8.0-1 2.0 3.0 Figure 18 Results when the line resistance is 832 Ohms. 3.0 2.0 - -2.0-3.0-2.0 3.0 Figure 20 Simulation results with the filed time constant set to a lower value. Prepared b Dharshana Muthumuni Manitoba HVDC Research Centre Inc.