Power system analysis and integration of the proposed Nigerian 750-kV power line to the grid reliability

The present situation of power generation in Nigeria obviously represents a challenge to our ability for rethinking the delivery of energy at maximum efficiency. Previous research on the existing Nigerian 330-kV network grid, recommended that the network be transformed from radial to ring because of high losses inherent in it and the voltage insecurity. In this study, the existing 330-kV network was reconfigured based on the identified regions mapped out for upgrade to form a ringed 750-kV super grid. The bus voltages of some of the buses in the existing 330-kV were upgraded to 750-kV and new transmission lines added to create an integrated super grid with a ring structure as compared to the radial nature of the existing 330-kV grid. These proposed buses have been selected for upgrade based on the fact that they are positioned in critical areas within the topology of the grid that transforms the existing radial structure to a ring one. The method is also cheaper than making the entire network a 750-kV system. Load-flow analysis was carried out on the existing 330-kV Nigerian Grid and the proposed Nigerian 750-kV integrated into the existing grid using Newton–Raphson algorithm. The results analysis of the new network revealed a significant reduction of 30.2% power loss. This was validated using the code-based MATLAB and Power World Simulation model-based software. Contingency analysis was also carried out on both grids using the Power World Simulator. The study revealed that the 750-kV super grid was able to mitigate the losses experienced on the existing grid significantly with better voltage profiles in all the buses. It also revealed that the new network (330-kV and 750-kV integrated) performed better to the single line contingency analysis with less violations occurring and no unsolvable cases.


Introduction
The operation of transmission lines at high voltage level, above 500-kV, has been embraced by the developed countries as a standard [1][2][3][4][5]. A reliable power supply is an essential part of any nation thriving toward economic growth and national development. Nigeria, being the giant of Africa and home to the largest economy in Africa as of 2014 [6], cannot boast of a steady power supply due to the state of the power sector operations. It is no longer news that the national power system has been identified with frequent outages. This advertently caused the cost of living to increase as citizens now have to consider the cost of running their private generators. Furthermore the negative role it plays in industrial and economic development has to be considered. Industries and companies refrain from production activities in the nation because of the outrageous cost of running their power sources. This has also caused some of the existing industries to collapse due to the unavailability of electricity thereby resulting in the loss of jobs in those industries and further increase in the hardship faced by the populace of this country [7].
Over the years, there have been reforms in the power sector to curb the many deficiencies which the power sector experiences. Notable progress has been made in the generating sector after the deregulation of the generating wing of PHCN. Regardless of that, the problem of power outages and unavailability of consistent power still plague innocent citizens. The reason is that the grid's weak infrastructure is characterized by high line losses and unbalance of reactive power between generation and load center due to forced outages and system disturbances [8][9][10][11], and cannot support the effective transmission of power to the areas where it is needed for utilization. The Nigerian power transmission system currently utilizes a 330-kV grid for power transmission across the country. The network is more of radial nature, and its infrastructure comprises effective 31 buses, 6,000 MW total installed power generation capacity, and 4,889.2 km length of transmission lines [12,13]. Several losses have already occurred in the existing grid inferring the inadequacy of the current transmission to deal effectively with the generation. The nation currently suffers the problem of under capacity and poor transmission due to the huge losses experienced on our transmission lines [14][15][16].
In light of this observed challenges, the major objective of the ongoing reform is to expand the existing grid network to a proposed mesh structure and increase the number of buses to a total of 49 buses; 32 buses from existing 330-kV plus 17 buses of the proposed 750-kV. However, this expansion remains a proposed one because no technical analysis shows its viability The rapid increase in the demand for electricity as a result of population growth, industrial development and rise in consumer electrical appliances has necessitated the stepping up of generation and demands similar of the transmission capabilities of the grid network to deliver quality power supply to the consumers. The present Nigerian transmission grid network is prone to evacuating and dispatching reliable and quality electricity supply and maintaining an operational standard of security to prevent any collapses [17]. Moreover, there are ever-increasing power demands and restrictions on land and space to build more transmission corridors globally [18]. All the above-aforementioned reasons set the platform for the modification of power system analysis of a 750-kV transmission line and integration to the grid to support the increase in generation and demand by citizens.
The power system analysis of the 750-kV transmission grid integrated into existing 330-kV involves the load-flow analysis and the contingency analysis for system security and reliability. The load-flow analysis of this system serves as an important tool in determining the results that are necessary for planning, design and operating of power system under steady state The power-flow solution is therefore used to evaluate the bus voltage, branch current, real power flow, reactive power flow for specified generation and load conditions. The idea behind this research work is that at 750-kV, the current through the transmission lines would be less thus a decline in the losses experienced across the transmission lines [19]. This research will show the comparison between the load-flow analysis of the 330-kV grid and the 750-kV integrated into the grid. The resulting data obtained from the outcome will be used to infer the appropriate solution for the advancement and modification of the existing transmission grid for maximum efficiency and continuous sustainable energy supply.

Methodology adopted for the work
The test system for this research work is a proposed 750-kV transmission grid. The grid is an integrated system comprising the existing 330-kV grid and the selected identified regions to upgrade to the 750-kV transmission voltage. Thirteen buses have been identified for the transmission voltage level upgrade to reduce the cost. The various locations of the upgraded buses are shown on the Nigerian site plan map of Fig. 1. These buses have been selected for an upgrade because they are positioned in critical areas within the grid's topology, transforming the existing radial structure to a ring one. A step-wise approach has been applied to upgrade this system due to the cost of implementing upgrading all the buses to 750 kV at once. In an attempt to ring the entire grid with the proposed thirteen buses mapped out from the super grid, two more lines were added: Mambilla-Jalingo and Jalingo-Gombe. This is in anticipation of power that will soon be evacuated from the Mambilla Hydropower plant, which is expected to be completed by 2018. Thus, the name and distance of every line as mapped out from the supper grid map are presented in Table 1. The transmission line data for the 16 lines of the upgraded 750-kV shown in Table 1 are calculated from the first principle, while the bus and generator data of the existing 330-kV grid were used as obtained from PHCN.
The single-line diagram of the upgraded 750-kV lines, as presented in Table 1, is shown in Fig. 1. The integration of this upgraded 750-kV grid with the existing 330-kV grid is illustrated in Fig. 2.

Load-flow analysis of the integrated Nigerian 750-kV transmission system
To carry out any comparative analysis between the integrated 750-kV system that formed a new network with the existing 330-kV system and the existing 330-kV grid system, it is necessary to perform the load-flow analysis of the system first. It is essential because the load-flow results are required for loss of generating unit, loss of transmission line, or failure of any single component [21]. The load-flow analysis is used to evaluate bus voltages, branch currents, and real and reactive power flow under steady-state balanced conditions. For this research work, load-flow analysis based on the Newton-Raphson technique was performed using Power World Simulator on run mode. The simulation software was used to determine the bus voltages, real and reactive power flow losses on the various lines under normal conditions for the 330-kV and the 750-kV integrated grids. The real and reactive power losses obtained were also verified using the MATLAB software to show the authenticity of the results obtained from the simulation of the network.

Newton-Raphson method
This method is a successive approximation procedure which is based on an initial estimate of the unknown and the use of Taylor's series expansion.
For a typical bus of power system, the net current entering bus I is given by Eq. (1).  Complex power is obtained from Eqs. (4) and (5) The real and imaginary powers are respectively given by Eqs. (6) and (7). (1) The matrix form of Newton-Raphson method is given as shown in Eq. (8) where; ΔP and ΔQ are respectively the real and reactive power mismatch, Δ|V | and Δ are respectively bus voltage magnitude and angle, J1 through J4 are called Jacobian sub-matrices.

Calculations of the transmission line data for 750-kV lines
The determination of the line parameters is carried out using the per-unit system. Typical base value of 1,000 MVA was adopted for the 750-kV transmission grid system using Eq. (11).
Thus, the impedance and admittance base values are obtained respectively as shown in Eqs. (12) and (13), respectively.
(15) X p.u = X actualvalue Z base ; The calculated per-unit values of the resistance, reactance, and susceptance were used to multiply the individual line lengths to generate the line parameters for the proposed 750-kV, as shown in Table 2

Power world simulator
Power world simulator (PWS) is the software used in this research work for simulation purposes. Both load-flow analysis and contingency analysis were carried out using this software. PWS is very user-friendly, and it is interactive simulation software. The simulator allows full-colored animated one-line diagrams with zooming and panning capability, which helps to increase the user's understanding. Load-flow simulation was carried out on the existing 330-kV and the integrated 750-kV super grid. The results obtained from both simulations served as the basis for comparison in terms of MW and MVar losses, and the number of voltages violated buses in the grid.

Load-flow simulation of the existing 330-kV transmission grid
To carry out load flow using the PWS, the single-line diagram of the existing 330-kV grid was drawn in the edit mode of the simulator. The main simulation was done in the run mode and the run mode executed in the simulation using N-R iterative method in order to obtain bus voltages, phase angles, real and reactive power flows, and line losses. The results were obtained after inputting line data, load data and generator data obtained from PHCN into the dialogue box of the PWS in edit mode. For this simulation, Egbin was set as the slack bus and it was used to simulate the load flow of the existing grid. Figure 3 shows the existing grid in the run mode of the power world program. The simulation of the network obtained the unknown bus angles and voltage levels at each bus.

Load-flow simulation of the 330-kV and integrated 750-kV transmission grid
In executing the load-flow simulation of the 750-kV transmission grid, a similar method was used as that of the 330-kV grid. The 750-kV grid integrated to the existing 330-kV was redrawn in the edit mode of the power world simulator as shown in Fig. 4. The line, bus and generator data were entered in the edit mode of the simulator. The main simulation was done in the run mode and the run mode executed for the simulation using N-R iterative method in order to obtain bus voltages, phase angles, real and reactive power flows, and line losses. The edit mode of the load-flow simulation is as presented in Fig. 5. The simulation of the network obtained the unknown bus angles and voltage levels at each bus. MATLAB was used to verify the losses recorded from the load-flow simulation of both the existing Nigerian 330-kV and the integrated 750-kV grids using the power world simulator. The simulation was done by imputing the same load and transmission line data for the 330-kV and 750-kV grids used on power world simulator into a MATLAB script developed to perform load-flow analysis and print total losses.

Power system contingency analysis
Contingencies are characterized as possibly harmful disruptions that happen during the steady-state operation of   Research Article SN Applied Sciences (2021) 3:864 | https://doi.org/10.1007/s42452-021-04847-3 a power system. A power system operating in its normal mode is prone to experiencing contingencies such as sudden loss of line or generator, sudden increase or decrease in power demand. These contingencies cause transmission line overloading or bus voltage violations. It represents an important tool to study the effect of components outages in power system security during operation and planning [22]. The contingency that analyses the model failure events of all credible outages, one after the other in sequence have been studied and executed on the power systems. It is also important for this research because it shows how both the 330-kV and integrated systems will respond to the given outage conditions, thus, providing validation for the integrated 750-kV and a foundation for the implementation of the system in the Nigerian national grid. This research used a simple contingency format to analyze the new test system and the existing 330-kV system. The contingency analysis implored was a simple single-line contingency analysis. The simulation was carried out using the PWS. The contingency analysis sequence of execution can be described in a simple way by the following flowchart.

Result analysis and discussion
In this section, the load-flows and contingency analysis results that were carried out on both the existing 330-kV and the configured new network (i.e., with the 750-kV integrated power line grid) are compared and analyzed. Thus, the load-flow results obtained for unknown bus angle, voltage magnitude, and real and reactive power losses at each bus using MATLAB and PWS in run-mode were analyzed and discussed for performance evaluation. Also, the simple single-line contingency analysis results employed in the study were analyzed and discussed to provide validation for the integrated 750-kV and a foundation for the implementation of the system in the Nigerian national grid.  The analysis of Fig. 6 shows that 11 buses had voltage violations under steady-state operation. These buses include Birinin-Kebbi; Akangba; Katampe; Onitsha; New Heaven; Alaoji; Shiroro; Kano; Kaduna; Jos; and Gombe, with the bus at Gombe recording the lowest voltage of 0.36pu.
Also, Fig. 7 showed that 14 power lines recorded losses higher than 4 MW with critical value along line 20-19, representing okapi to Onitsha that recorded

Load-flow analysis of integrated 750-kV super grid (new model)
The load-flow analysis of the new model using the Power Word Simulator produced the values of the bus voltage magnitudes, voltage angles, and real and reactive power flow losses on the various lines under normal conditions. Figure 9 shows the voltage profile of the N-R in PWS load-flow with bus voltage violations at Birnin-Kebbi, New Heaven, Alaoji, Kano, Gombe, Jalingo and Mambilla buses. The real, and reactive power losses are graphically represented in Figs. 10 and 11, respectively. The analysis of the results showed that under steady-state operation, 7 buses were slightly below the specified voltage tolerance, with the bus at New Heaven recorded the lowest voltage of 0.883pu. Analysis of Fig. 10 showed that 8 lines recorded losses higher than 4 MW and these lines include 1-2; 3-5; 8-10; Also, the total real and reactive power losses recorded were 173.21 and − 846.39 Mvar respectively. These show that the combined new model of 330-kV and 750-kV super grid was able to reduce the real power losses by 30.2% as shown in Fig. 12 and the reduction in the reactive power losses shows shortages in reactive power. While Fig. 13 shows the combined real power losses of 330-kV and 750-kV transmission line. Table 3 shows a side-by-side comparison of the total losses obtained when the load-flow analysis was simulated in both the PWS and MATLAB software. The result in Table 3 validates the authenticity of the losses recorded using the Power World Simulator. The Real Power Loss differences between 330-kV and 750-kV grids were 78.44 MW and 69.74 MW in PWS and MATLAB, respectively. This shows that they are within the same range. The percentage loss differences between both grids were also similar with 30.2% and 37.87% in PWS and MATLAB, respectively. This validates that loss reduction will occur in the integrated 750-kV Super Grid. Tables 4 and 5 show the results of single-line contingency analysis carried out on the existing 330-kV grid and the combined 330-kV grid integrated with 750-kV super grid respectively. The contingency results were evaluated by subjecting the various lines to a single contingency analysis to examine the effect of loss of any single line on the network. The analysis of the results in Tables 4 and 5 show the following: When the existing 330-kV was subjected to contingency analysis, a total of 201 voltage violations with 2 unsolved cases were recorded as shown in Table 4 The integrated 750-kV super grid responded better to the single-line contingency analysis by recording only 70 voltage violations with 0 unsolved case as presented in Table 5.

Conclusion and recommendation
So far, it has been established that the present Nigerian 330-kV power grid is characterized by lower power availability to consumers. This is a consequence of some identified poor power qualities like power losses, voltage violations, and network topology. These reasons set the platform for the motivation to develop a new network system combining the existing 330-kV grid and some buses selected for upgrade to 750-kV super grid integrated to the existing 330-kV grid to transform the existing radial structure to a ring topology at a reasonable cost. The performance outcome of the power system analysis revealed a significant improvement in the voltage profile of voltage violated buses in the existing 330-kV grid. There was also a considerable reduction in transmission line losses by 30.2% compared to the storage, and impacts of the dynamic thermal rating on wind-integrated power networks to increase the reliability enhancement of the modified model.

Acknowledgements
The authors wish to acknowledge the management of Covenant University for her part sponsorship and support toward the success of this research work. " Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.