^{1}

^{*}

^{1}

^{1}

Transient analysis of 33 KV power transmission line stability of Egi communi-ty is considered in this research work with the aim of reducing the frequency of fault occurrence and voltage collapse in the network. The supply is taken from Egi generating station located at Total Nigeria Limited Gas Plant Obite at voltage level of 33 KV to Egi communities. This work focuses on the transient nature of network stability since transient fault is the most dangerous in elec-trical systems. The swinging of the generator rotor in the event of transient three-phase short circuit fault can be monitored by the circuit breakers and the protective relays which causes mal-functioning of the circuit breakers and pro-tective relays leading to abnormal behavior of the network. Therefore, data obtained from the power station were used as a case study of Independent Power Producer (IPP) in Nigeria. For investigation of the power angle, angular velocity, rotor angle differential changes, and angular velocity differential changes, an electrical transient analyzer tool was employed (ETap version 16.00) for circuit breaker and protective relay time setting of (0.00, 0.05, 0.10, 0.15, 0.20, 0.25, 0.30, 0.35, 0.40, 0.45, 0.50, 0.55, 0.60). The work used the Trapezoidal numerical technique for data analysis. The graphs were plotted using Matlab R2015a and the results obtained showed that when a symmetrical three-phase short circuit fault occur at one or any of the feeders, the fault must be cleared as quick as possible through the coordination of the circuit breakers and protective relays. For this research work, 17 cycles corresponding to relay time setting of t = 0.34 s were recommended and at each cycle, changes in time with respect to changes in rotor angle, angular velocity, rotor differential and angular velocity differential were calculated on the power network simultaneously. The results demonstrated that the Trapezoidal method is numerically stable, accurate and has faster respond time when compared to Modified Euler and swing equation techniques in event of fault occurrence in network.

Electrical energy demand in Egi community is continuously increasing each year in both the domestic and industrial sectors as most villages are hooked up to Egi network and industries and businesses are established to meet the demand of the teeming population. The establishment of Egi Gas Turbines by Total Exploration and Production Nigeria Limited (TEPNG) to supply electrical power to the host communities has improved the economic activities of the people; hence electrical energy will continue to be the single most important energy in these communities as it plays a crucial role in the development of any society. To distribute electrical power from Egi Gas Turbine to the neighboring communities, transmission and distribution lines were built and they are far away from the power station and from each other. The power station output is connected and transmitted through the Egi Transmission and Distribution Network (ETDN) to the respective communities [

The interconnected electric power network (ETDN) present various challenges to Total Exploration and Production Nigeria Limited, the operator of the power station, these problems are due to inherent nature of power system through planning, construction, operation and control. One of these challenges is transient stability or voltage stability. The Transient instability in Egi transmission network is caused by a number of factors which ranges from atmospheric phenomenon like lightening, solar fares, geomagnetic interference, sudden addition of inductive load, switching load ON and OFF, switching of power lines especially during repairs, interruption of fault current by breakers etc. The instability in the power network had led to the collapse of the feeders of Egi Gas Turbine on several occasions [

Steady and efficient electricity supply from TEPNG resulted in massive industrialization within the Egi communities, as more welding shops, salon shops, grinding shops and businesses are evolving everyday which translates to more inductive load on the network, the resultant effect is that by 2007, the network could no longer be synchronized. This led to voltage instability in the system which results to constant load shedding and delivering of poor power supply to the communities [

1) Complete collapse of feeder 1 and 2 of Egi Gas Turbines as a result of instability in the network.

2) Continuous addition of inductive load on the system by establishment of welding and fabrication workshops, wood sawing workshops, grinding shops and palm/kernel mills resulting to instability in the network.

3) Increasing demand of power for domestic use due to influx of population and expansion of the communities by building new houses which add more load to the network causing system instability.

4) Establishment of high local micro businesses like salon and barbing shops, frozen fish/meat shops, cooled drink and restaurants etc. are also on the increase which add load that affect system stability.

5) The current load demand exceeds the power generation capacity hence leading to total system collapse.

The existing total load consumption would continue to increase, hence, the need to carry out forecast study for effective operation of the plant [

This work investigated the behaviour of transient stability of 33 KV transmission network of Egi community with emphasis on the behaviour of 9.455 MVA synchronous generators after a large disturbance on the power network. The generation, transmission and distribution data were collected and analyzed. The data of the power network collected and the results generated were analyzed. An Etap analyzer was used and results were generated. The changes in time (t) against changes in power angle and changes in time (t) against changes in angular velocity were investigated using analytical tool. The swing equation model and Trapezoidal rule were used to determine the accelerating power of the synchronous generators when the load was suddenly removed. The coordination of the network relays and the circuit breakers is crucial to the stability of power network; the relays are programmed to sense the fault and transmitted to the circuit breakers to isolate any faulty part of the system to avoid system collapse. The faster these sensors act in micro-second, the more stable the network. Electrical fault in power network can be categories into three (3) distinct parts namely, the pre-fault condition, fault condition and post fault condition [

Electrical faults in power systems are non-linear in nature therefore requires non-linear solutions.

In determining the transient stability of power system, the active power transferred by generator plays important role since the stability of the system is influenced by the angle stability of the generators which are connected to a complex network of inter-connected power system [

In normal operating condition, the relative position of the rotor axis and the stator magnetic field are fixed [_{L} is applied suddenly, the generator may lose synchronism (fall out of step) even if the load torque is less than maximum torque δ_{max} the generator can develop. The maximum value of the load torque that can be applied suddenly to the generator without the generator losing synchronism can be obtain from the dynamic behavior of the synchronous generator.

Let consider an unloaded synchronous generator connected to a power system, if losses are neglected, the torque angle δ_{g} is zero. If a load torque δ_{L} is suddenly applied to the shaft of the synchronous generator, the generator speed slows down and the torque angle δ_{g} increases. As the torque angle δ_{g} increases, the generator develops more torque δ_{g} to meet the load torque δ_{L}. The torque angle δ_{g} reaches the value of the load torque δ_{L} at which the load torque δ_{L} equals the torque angle δ_{g} of the generator as shown in

However due to inertia the torque angle δ_{g} increases beyond the load torque δ_{L} and the generator develops more torque than required for the load. The deceleration decreases and the torque angle δ_{g} reaches a maximum value δ_{max} and then swings back. The torque angle δ_{g} oscillates around the load angle δ_{L}, due to damping in the system, the torque angle will settle down to the required value of load angle δ ′ L . The oscillation of the torque angle δ_{g} last for several cycles and it can be assumed that

The oscillation in the torque angle δ_{g} as a function of time can be obtained in solving the differential equation that describes the behavior of the system [

The excitation voltage behind the transient reactant of a synchronous generator is expressed as:

E ′ g = V + J X ′ d I (1)

where E ′ g = excitation voltage;

V = bus voltage;

X ′ d = transient reactance [

The total reactance during the period of transient stability is the addition of transformer reactance X t , the line reactance X l and the transient reactance X ′ d ,

X e q = X l + X t + X ′ d (2)

where X l = Line reactance;

X t = transformer reactance;

X ′ d = transient reactance [

The equation which describe the relative change of position of the rotor (the torque angle) with respect to the stator as a function of time is the swing equation which can be expressed in terms of toque or power of the machine, hence,

P a = P s + P e (3)

P a = Accelerating power;

P s = Mechanical power;

P e = Eletromagnetic power [

The electromagnetic or electrical power is equivalent to the product of the generator voltage, the electromagnetic voltage developed by the generator, the sine of the torque angle and the reciprocal of generator reactance.

P e = V g V m X g sin δ g (4)

V g = generator voltage;

V m = Electromagnetic voltage developed by generator;

X g = Generator reactance;

δ g = generator torque angle [

The shaft power (mechanical power) is the power transferred or developed by shaft of the generator due to load and is given by

P s = E g E l X d sin δ l (5)

E g = Generator voltage;

E l = Load voltage;

X d = Line reactance;

δ l = Load torque [

Appling the principle of rotation, the swing equation that relates the inertia of prime mover and generator is:

K j d 2 δ d t 2

where j = moment of inertia of prime mover and generator;

Kj = inertia constant of machine [

The power balance equation of a synchronous rotating generator can be expressed as,

E g E l X e q sin δ l = K j d 2 δ d t 2 + V g V m X e q . sin δ g

P max sin δ l = K d d δ d t + K j d δ 2 d t 2 (6)

This is known as the swing equation and the trapezoidal rule is used to evaluate the load angle ( δ l ).

The network bus loading in the case study is the maximum and minimum load range from Egi Gas turbine data base. The generator and network data used for this work are presented in Tables 1-4.

Bus | From | To | Length (Km) | Resistance R (pu) | Inductance X (pu) |
---|---|---|---|---|---|

1 | Plant | Itu-Ogba | 21 | 0.000259 | 0.03498 |

3 | Plant | Ohali-elu | 35 | 0.000145 | 0.050553 |

4 | Plant | Akabuka | 10 | 0.0000399 | 0.01210 |

6 | Plant | Ama | 10 | 0.000132 | 0.03882 |

S/N | Parameters | Value Assign |
---|---|---|

1 | H Constant | 5.673 p. u |

2 | MW | 30 MW |

3 | MVAR | 22.4 MVAR |

4 | Frequency (f) | 50 Hz |

5 | Ammature Resistance (R_{a}) | 0.0252 p.u |

6 | Alternator internal Reactance (X_{d}) | 0.76 p. u |

7 | Alternator Inductance (L) | 0.000242 p. u |

8 | Generator Internal Voltage (E) | 5.0000 |

9 | Infinite Bus Voltage (V_{inf}) | 1 |

10 | Turbine Power | 31.6 MW |

Bus | MW (max) | MW (Min) | MVAR (Max) | MVAR (Min) |
---|---|---|---|---|

1 | 1.8 | 0.0007 | 0.00144 | 0.00056 |

3 | 1.5 | 0.0004 | 0.0012 | 0.00032 |

4 | 2.5 | 0.0009 | 0.002 | 0.000672 |

6 | 1.2 | 0.0004 | 0.00096 | 0.00032 |

S/N | Bus | Designation | Description |
---|---|---|---|

1 | 1 | Feeder 1 | Uzo-uzimini/Ahiahu communities |

2 | 3 | Feeder 3 | Umur-obor communities. |

3 | 4 | Feeder 4 | EgiEtiti (excluding Obite) |

4 | 6 | Feeder 6 | Obite/Ede/Ama |

Using MVA base of 20 MVA

P e = 28.4 20 = 1.42 p . u

P m = 31.6 20 = 1.58 p . u

Reactance:

X e q = generator reactance + line reactance + transformer reactance [

X e q = 0.19 p . u + 0.013 p . u + 0.03 p . u

X e q = 0.35 p . u

Real Power (S) [

S = P e P f ∠ cos − 1 ( ∅ ) (7)

where P f = ∅ = 0.8 [

S = 1.42 0.8 ∠ cos − 1 ( 0.8 ) = 1.775 ∠ 36.87

Current:

I = S * V * = 1.775 ∠ − 36.87 1.0 ∠ 0

I = 1.775 ∠ − 36.87

Excitation Voltage [

E ′ g = V + J X ′ d I

Therefore, we have:

E ′ g = 1 ∠ 0 + ( j 0.35 ) ( 1.775 ∠ − 36.87 )

E ′ g = 1 + j 0 + ( 0.35 ∠ 90 ) ( 1.775 ∠ − 36.87 )

E ′ g = 1 + j 0 + 0.3728 + j 0.4970

E ′ g = 1.3728 + j 0.4970

E ′ g = 1.3728 p . u

Therefore, we have the initial operating condition as follows [

H = 5.673 p .u

ω 0 = ω L = 2 π f o = 2 × π × 50 = 314.1593 rad / sec

δ 0 = δ L = 23.2 ∘

1 rad = 57.32 ∘

Therefore, 23.2 ∘ in rad = 23.2 ∘ 57.32 ∘ = 0.4047 rad [

Applying the values into the swing equation technique [

H d 2 δ π f d t 2 = P m − P e (8)

where H d 2 δ π f d t 2 = P a = accelerating power [

Substituting into Equation (8) we have:

5.673 π × 50 d 2 δ d t 2 = ( 1.58 − 0 ) (during fault)

5.673 π × 50 d 2 δ d t 2 = 1.58

d 2 δ d t 2 = 1.58 × π × 50 5.673

Integrating both sides:

d δ d t = ∫ π × 50 5.673 d t × 1.58

d δ d t = ∫ π × 50 5.673 d t + 0

Integrating the second time:

d δ ( t ) = π × 50 11.346 t 2 + δ 0

d δ ( t ) = π × 50 11.346 t 2 + 0.4047

δ ( t ) = π × 50 11.346 t 2 + 0.4047

The iteration is repeated for thirtieth cycles.

Applying the Trapezoidal rule and the two mechanical equations that describe the rotor dynamics, the trapezoidal rule can be investigated for thirtieth cycles as follow [

δ n + 1 = δ L + Δ t 2 ( d δ n d t + d δ n + 1 d t ) (10)

ω n + 1 = ω L + Δ t 2 ( d ω n d t + d ω n + 1 d t ) (11)

where,

δ n + 1 = machine rotor angle evaluated at n + 1;

δ L = Latest rotor angle;

d δ n d t = differential rotor angle evaluated at n;

d δ n + 1 d t = differential rotor angle evaluated at n + 1;

ω n + 1 = angular speed evaluated at n + 1;

ω L = Latest angular speed;

d ω n d t = differential angular speed evaluated at n;

d ω n + 1 d t = differential angular speed evaluated at n + 1;

Δ t = 0.02 s , interval between cycles and is constant throughout the period of iteration.

d δ n d t = ( ω L − 2 π f o ) [

2 π f o = Synchronous speed

d ω n d t = 1 M ( P m − P e ) [

where M = H π f o

d ω n d t = 1 H π f o ( P m − P e )

d ω n d t = 50 π H ( P m − P e ) (13)

where H = inertia constant;

P m = mechanical power;

P e = electrical power.

During fault, P e = 0 [

Iteration at t = 0

Evaluation of ( d δ d t , d ω d t , δ , ω )

Substituting in Equation (12)

d δ 0 d t = ( ω L − 2 π f o )

d δ 0 d t = ( ω L − 314.1593 )

Here ω L = 314.1593 ,

d δ 0 d t = 314.1593 − 314.1593

d δ 0 d t = 0

Again, from Equation (13)

d ω 0 d t = 50 π H ( P m − P e )

d ω 0 d t = 50 × 22 7 11.346 ( 1.58 − 0 )

d ω 0 d t = 13.8501 ( 1.58 − 0 )

d ω 0 d t = 21.8831

Applying Trapezoidal rule at 0

δ ( 0 ) = δ L + Δ t 2 ( d δ 0 d t )

δ ( 0 ) = 0.4047 + 0 2 (0)

δ ( 0 ) = 0.4047 rad = 23.2 ∘

ω ( 0 ) = ω L + Δ t 2 ( d ω 0 d t )

ω ( 0 ) = 314.1593 + 0 2 ( 21.8831 )

ω ( 0 ) = 314.1593 rad / sec

[ δ ( 0 ) = 0.4047 rad = 23.2 ∘ ω ( 0 ) = 314.1593 rad / sec ]

δ L = 0.4047 rad

ω L = 314.1593 rad / sec

The iteration is repeated for thirtieth cycles.

Calculation of Critical Clearing Time (CCT)

Critical clearing time, t c r [

t c r = ( δ c r − δ 0 ) ⋅ 4 H ω s P m [

where,

t c r = critical clearing time [

δ c r = critical clearing angle = 89.63˚ = 1.5637 rad;

δ 0 = initial angle = 23.2˚ = 0.4047 rad;

ω s = synchronous frequency = 314.1593 rad/sec;

P m = mechanical power or shaft power = 31.6 MW;

H = inertia constant of generators = 3 MJ/MVA.

Substituting,

t c r = ( 1.5637 − 0.4047 ) × 4 × 3 314.1593 × 31.6

t c r = 0.037 s ≈ 0.04 s

It implies that for effective fault clearing operation of the relay coordination and circuit breakers, the relay setting time should be 0.04 s. This will clear the three-phase fault as fast as possible for the network to be in synchronism with the machine to avert system collapse.

The results of the rotor angle and the angular velocity are represented for detail analysis of the power network. The results show the inherent nature of the network variables when subjected to external disturbance and occurrence of fault on any part of the network. The results are therefore presented as follows (Figures 4-7 and Tables 5-7).

The graph of

Time | δ (rad) | δ (deg) |
---|---|---|

0 | 0.4047 | 23.2 |

0.04 | 0.4091 | 23.44 |

0.08 | 0.4376 | 25.08 |

0.12 | 0.5076 | 29.10 |

0.16 | 0.6126 | 35.11 |

0.20 | 0.7527 | 43.15 |

0.24 | 0.9278 | 53.18 |

0.28 | 1.1379 | 65.23 |

0.32 | 1.3830 | 79.28 |

0.36 | 1.6631 | 95.33 |

0.40 | 1.9512 | 111.84 |

0.44 | 2.1939 | 125.76 |

0.48 | 2.3915 | 137.08 |

0.52 | 2.5537 | 146.38 |

0.56 | 2.6917 | 154.29 |

0.60 | 2.8165 | 161.44 |

Time (s) | δ (Swing, deg.) | δ_{2} (Mod. Euler, deg.) | δ (Trapezoidal, deg) |
---|---|---|---|

0 | 23.25 | 23.25 | 23.25 |

0.02 | 23.67 | 23.45 | 23.25 |

0.04 | 24.47 | 25.46 | 23.49 |

0.06 | 26.05 | 31.71 | 23.82 |

0.08 | 28.28 | 47.77 | 25.08 |

0.10 | 31.14 | 79.13 | 26.84 |

0.12 | 34.63 | 111.30 | 29.10 |

0.14 | 38.76 | 121.61 | 31.81 |

0.16 | 43.52 | 106.13 | 35.11 |

0.18 | 48.92 | 41.49 | 38.88 |

the rotor angle to returned to stability, (rd3) is the rotor angle in the case of fault up to 3.5 cycles which was cleared for the network to returned to stability while (rd4) is the rotor angle at fault up to 5.5 cycles which lasted for much longer time to be cleared before the network settle to stability

Appendix B1 shows when a three-phase fault occurred at feeder 1 and the result of the simulation present the swinging of generator rotor angle which indicate unstable state. The swing curves show that the fault region is near to generator. Appendix B2 shows generator rotor angle when the three-phase fault was cleared at 1.5 cycles and the network returned to stability region.

Appendix B3 shows generator rotor angle with fault cleared at 5.5 cycles, initially the network was unstable at the occurrence of the fault at t = 0, but

t | δ_{1} | δ_{2} | δ_{3} | δ_{4} |
---|---|---|---|---|

0 | 31.69 | 31.69 | 31.69 | 31.69 |

0.05 | 31.90 | 32.33 | 33.10 | 32.99 |

0.10 | 32.99 | 33.77 | 35.78 | 37.30 |

0.15 | 35.11 | 35.75 | 38.89 | 41.62 |

0.20 | 38.15 | 38.04 | 42.05 | 45.78 |

0.25 | 42.21 | 40.18 | 45.49 | 49.25 |

0.30 | 47.19 | 42.33 | 46.92 | 51.57 |

0.35 | 53.16 | 43.98 | 47.97 | 52.76 |

0.40 | 60.02 | 44.45 | 47.85 | 51.62 |

0.45 | 67.85 | 44.22 | 46.56 | 49.33 |

0.50 | 76.55 | 42.80 | 44.26 | 45.82 |

through the coordination action of relays and circuit breakers, the network was immediately brought to synchronism at 5.5 cycles after clearing the fault.

This research work described the behavior of transient stability of 33 Kv transmission network of Egi community with emphasis on the behavior of 9.455 MVA synchronous generators after a large disturbance on the power network. The generation, transmission and distribution data were collected and analyzed. The data of the power network collected and the results generated were analyzed. An Etap analyzer was used and results were generated. The changes in time (t) against changes in power angle (δ) and changes in time (t) against changes in angular velocity (ω) were investigated using analytical tool. The swing equation model and Trapezoidal rule were used to determine the accelerating power of the synchronous generators when the load was suddenly removed. Stability could not be obtained even when network fault was cleared at t = 17 cycles = 0.34 s (FCT) and power angle of 87.05˚, above 17 cycles the swinging caused by network disturbance causes network instability. The results of the swinging curves and load flow were also presented at the appendices.

Successful operation of power system plants depends on the ability of the power engineers to adopt and apply the recommendations resulting from transient and load flow stability studies to prevent network collapse. If these studies are taken seriously, it will prevent total network break down and collapses.

The following points should be noted from the investigation and study of the transient stability analysis of 33KVEgi transmission network for considerations:

1) Continuous evaluation of reactive power and the power angle (machine rotor angle).

2) Fast and high-speed operation of protective relays and circuit breakers with fault clearing time set at t = 0.04s.

3) Expansion of the generating units to meet the load demand.

4) Installation of time-synchronizing phasor measurement units to monitor network coordination.

5) Minimization of transfer reactance to increase power.

The authors declare no conflicts of interest regarding the publication of this paper.

Elechi, P., Ahiakwo, C.O. and Okwu, U.N. (2020) Transient Stability Analysis of 33 KV Transmission Network of Egi Community, Nigeria. Journal of Power and Energy Engineering, 8, 18-35. https://doi.org/10.4236/jpee.2020.83002