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 . Generally, transient instability is a major challenge in power system operation and control. At present, feeders 1 and 2 have just been reinstalled into operation which collapsed as a result of instability in the network. In order to reduce the frequency and duration of load shedding, a total of 3 MW was formally transferred to the Rivers State Generation Plant at Omoku. It is on this basis that the following problems were identified:
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 . For stability of power system, the fault clearing time (FCT) and the critical clearing time (CCT) has to be determined and specified .
2.1. Materials and Method
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. Figure 1 represent the equivalent machine model of a synchronous generator, the excitation voltage E' is behind the transient reactance . This research work will adopt the application of implicit method using Trapezoidal Rule in resolving transient stability of 33 KV transmission network of Egi community, Nigeria and compare it with Modified Euler and the nominal Swing equation techniques. Figure 2 is the power transfer curve. Figure 3 is the single line diagram of the power supply from Egi gas turbine to the communities. Figure 4 is the simulated network from Egi gas turbine to 33 Kv transmissions and distribution network. Figure 5 is the graph of incremental changes of machine rotor angle against time, while Figure 6 is comparison of changes in machine rotor of Swing, Modified Euler and Trapezoidal methods against changes in time. Figure 7 is the simulated swing curve of the rotor angle of the case study.
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  and .
In normal operating condition, the relative position of the rotor axis and the stator magnetic field are fixed . The angle between the rotor axis and the stator field is known as the torque angle, if the loading of the generator is large, the torque angle δ will invariably be large . If a load torque TL 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 Figure 2.
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 . The oscillation of the torque angle δg last for several cycles and it can be assumed that Figure 1 shows satisfactory representation of the electrical system.
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  and .
The excitation voltage behind the transient reactant of a synchronous generator is expressed as:
Figure 1. Equivalent machine model .
Figure 2. Power transfer curve .
Figure 3. Single-line diagram of power supply from Egi station to the communities.
where = excitation voltage;
V = bus voltage;
= transient reactance  and .
The total reactance during the period of transient stability is the addition of transformer reactance , the line reactance and the transient reactance ,
where = Line reactance;
= transformer reactance;
= transient reactance  and .
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,
= Accelerating power;
= Mechanical power;
= 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.
= generator voltage;
= Electromagnetic voltage developed by generator;
= Generator reactance;
= 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
= Generator voltage;
= Load voltage;
= Line reactance;
= Load torque .
Appling the principle of rotation, the swing equation that relates the inertia of prime mover and generator is:
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,
This is known as the swing equation and the trapezoidal rule is used to evaluate the load angle ( ).
2.2. Bus Data
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.
Table 1. Network parameters.
Table 2. Generation system parameters (Egi Gas Turbine).
Table 3. Bus loading for the transmission network.
Table 4. Study case designation, (Egi Gas Turbine).
2.3. Calculation of Line Network Parameters
Using MVA base of 20 MVA
= generator reactance + line reactance + transformer reactance  and .
Real Power (S) :
Excitation Voltage  and :
Therefore, we have:
Therefore, we have the initial operating condition as follows :
2.4. Application of Swing Equation Technique to Transient Stability
Applying the values into the swing equation technique :
where  and .
Substituting into Equation (8) we have:
Integrating both sides:
Integrating the second time:
The iteration is repeated for thirtieth cycles.
2.5. Application of Trapezoidal Rule to Transient Stability Analysis
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 :
= machine rotor angle evaluated at n + 1;
= Latest rotor angle;
= differential rotor angle evaluated at n;
= differential rotor angle evaluated at n + 1;
= angular speed evaluated at n + 1;
= Latest angular speed;
= differential angular speed evaluated at n;
= differential angular speed evaluated at n + 1;
, interval between cycles and is constant throughout the period of iteration.
 and . (12)
= Synchronous speed
where H = inertia constant;
= mechanical power;
= electrical power.
During fault,   and .
Iteration at t = 0
Substituting in Equation (12)
Again, from Equation (13)
Applying Trapezoidal rule at 0
The iteration is repeated for thirtieth cycles.
Calculation of Critical Clearing Time (CCT)
Critical clearing time, :
= critical clearing time ;
= critical clearing angle = 89.63˚ = 1.5637 rad;
= initial angle = 23.2˚ = 0.4047 rad;
= synchronous frequency = 314.1593 rad/sec;
= mechanical power or shaft power = 31.6 MW;
H = inertia constant of generators = 3 MJ/MVA.
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.
3. Results and Discussion
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 Figure 5 shows the characteristics behavior of the machine rotor angle (δ) of the Trapezoidal method against time (t). The machine rotor angle in degree gradually increased from when the fault was initiated at t = 0 and rotor angle of 23.2˚ up to maximum allowable stability angle of 90˚ and continued up to 161.44˚ at t = 0.6 s. The network is very unstable and cannot be synchronized because the rotor angle keeps increasing even when the speed has reduced. The increase in machine rotor angle beyond 90˚ is an indication of instability in the network which will trigger loss of synchronization if corrective action is not taken. It is shown that the fault was cleared at the seventeenth cycles at t = 0.34 s (FCT) and rotor angle of 85.05˚, but the machine rotor angle continues to increase up to 161.44˚, even when the angular speed has decreased suggesting that the machine and the network is unstable and the generator is running under heavy loads which will lead to system collapse. Table 6 shows the incremental time against the rotor angle for swing, modified Euler and trapezoidal methods.
Figure 4. Simulated Network from Egi Gas Turbine to 33 KV Transmission network.
Table 5. Incremental time (t) with changes in rotor angle (δ) for trapezoidal method.
Figure 5. Graph of machine rotor angle (δ) against time for trapezoidal method.
3.1. Analysis of the Simulated Case
Figure 6 is the presentation of the machine rotor angle analysis in the case study at different operating time. Rotor angle (rd1) represent the case of fault which lasted for several time resulting to instability in the system, (rd2) is the rotor angle in the case of fault but was cleared very quick after 1.5 cycles which enable
Table 6. Incremental time (t) against rotor angle (δ) for swing, modified euler & trapezoidal methods.
Figure 6. Presentation of comparison of swing, modified euler and trapezoidal curves.
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
Table 7. Changes in time (t) against changes in rotor angle (δ1, δ2, δ3, δ4) obtained from simulation.
Figure 7. Simulated swing curve of rotor angle against time (t).
through the coordination action of relays and circuit breakers, the network was immediately brought to synchronism at 5.5 cycles after clearing the fault.
4. Conclusion and Recommendations
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.
 Omorogiuwa, E. and Elechi, P. (2015) Economic Effects of Technical and Non- Technical Losses in Nigeria Power Transmission System. IOSR Journal of Electrical and Electronics Engineering, 10, 61-69.