Steady State Simulation of 33 kV Power Grid

An example is presented in this paper relating to power problems in Sandakan power network. Sandakan is a suburb in east coast of Sabah state of Malaysia. The problems were reported with power flow and N-1 contingency in terms of blackout after main grid supply outages with overload and high fault current on distribution system. This paper focuses on analysis of steady state stability of 33 kV power grid using load flow, contingency analysis and voltage stability (P-V Curve). The analysis is done by using industrial grade power software called Power System Simulation for Engineers (PSS/E). The power flow result showed that there are three generators generating out of limits. Contingency results showed that three transformer branches and three distribution branches are affected. There are nine weakest buses that violate contingency voltage deviation criterion and cannot withstand more load power when N-1-1 case in the branch of bus 8 to bus 28 is under outage.


Introduction
In Sabah grid, the power demand is increasing annually but the generated capacities is less than the power demand, especially in the east coast [1].The current equipments in transmission system such as electricity cable are getting old and are found operated closer to their limits of stability and cannot withstand with increased power supply [2].Due to these, the grid system is exposed to disturbances or contingencies which can cause system collapse and blackout.Sabah faced two major serious blackouts [1].The most severe blackout happened in 2014 for 10 hours of state-wide blackout.The collapse is triggered by flashover which is from conductor of 132 KV transmission line.Another one is blackout in the whole east coast due to outage of 275 KV transmission line.
As stated in Sabah grid Code [3], most of transmission lines, distribution and transformers fulfil N-1 contingency requirement, but in fact studies showed that there is no N-1 contingency in aging existing transformer and existing line configuration [2].Therefore, when N-1 contingency happens, overload conditions occur on those transformers, distribution and transmission lines.
Contingency analysis has been developed by [4] using sensitivity factors to approximate power flow on branches whereas voltage performance index is used by [5] to approximate the contingency voltage on a certain bus in a power system.For voltage stability part, [6] reviewed four commonly used voltage stability analysis tools which are PV/QV curve analysis, L index, Modal analysis and V/Vo index.Authors have done comparison of accuracy results on IEEE bus power system.This paper focuses on steady state stability for distribution level of power grid.Thus, 33 KV Sandakan network of Sabah Grid System is chosen to simulate steady state stability which consists of load flow simulation, contingency simulation and P-V curve.Contingency scenarios are created to test its overall steady state stability of the grid in terms of contingency voltage deviation violation and percentage overload.Moreover, P-V analysis is performed to determine the weakest buses in the network.The process of analysing the stability can be daunting and challenging if the power network is highly complex, large size and non-linear.The process takes a lot of time in the calculations to access all the power variables and contingencies [7].Therefore, a Power System Simulation for Engineers (PSS/E) software is utilized to perform all power flow computations in this steady state analysis.

Existing N-1 Network
The existing model of power system network shown in Figure 1 is modelled by using (PSS/E).This network is disconnected from the main grid power supply, especially the supply from 275 KV transmission line.Therefore, in that case, the network itself is assumed as N-1 under outages of main grid and is considered as external N-1 condition.Date of bus names, Transformers, generators, loads and branches regarding powers and impedances are given in the tables at the Appendix.
The grid system has the following major components:

Load Flow Analysis
Load flow analysis is used to calculate voltage, current flows, active and reactive There are four power elements that load flow used to calculate.These are: 1) Generators (swing buses): generated active and reactive power flows 2) Distribution and transformer branches: current flows, active and reactive power flows, percentage voltage drop, power factor, and power losses.
3) Buses (load buses and generator buses): bus voltage, active and reactive power flows, and current flows.
4) Loads: active and reactive power flows, current flows, and percentage loading.
Table 1 shows the power flow results of generated active power and reactive power for each in-service generating unit which act as swing bus.P max and P min are the maximum and minimum generated capacity for each generator.Same goes to Q min and Q max which represent reactive power capacity.With respect to PGen, all generators generated within their active power and reactive power limits except for generators named KB_6.6,SB_6.6, and BN_11.
Table 2 shows the total generated power by all generators before and after the Table 3 shows the load flow results for all in-service branches (transformer and distribution).Note: DB is distribution branch, TB is transformer branch.
Transformer branches are found higher active and reactive power losses, higher percent of voltage drop, higher current flows and higher active and reactive power flows than that in distribution branches.Only transformer branches have higher loadings problem compared to distribution branches.For distribution branches, DB10_1 is recorded with highest power losses.For transformer branches, TB3 is recorded with highest reactive power loss and TB5 and TB6 share the highest active power losses.TB2 is recorded with highest percent of voltage drop among all branches, followed by DB2_1 and DB2_2.DB14_1, DB14_2, DB17_1 and DB17_2 are recorded with zero among all variables.The voltage range is from 0.9 pu to 1.1pu.That means all bus voltages are safe and secure under N-1 case which is disconnected from main grid supply.
This load flow analysis shows that all bus voltages are slightly higher than 100% and 2 load buses are overloaded as well as 3 generators generating out of limits.

Contingency Analysis
To compute the branch power flows after certain level of outage, contingency sensitivity factors are used to approximate the change in line flows and the changes in generation in a power system.It is one of fastest way to calculate possible overloads in a power system network.The main two sensitivity factors are Generation Shift Factors (GSF) and Line Outage Distribution Factors (LODF).
For GSF part, the generation factor is defined as changes in power flow in particular line when a change in power generation at reference bus occurs.
Generation shift factors,   In PSS/E software, the process of performing contingency analysis is done automatically and comprehensively without manually tripping each line.Before carrying out N-1 contingency analysis, three types of files are created.They are *mon.file,*sub.file and *con.file.Each file is described in Table 6. Figure 2 shows the process of creating these files.

*con.file
It is used to trip line or power elements to create contingency scenarios.Three types of contingencies: N-0 (system intact), N-1 (single power element outage) and N-2 (two power elements outage).In this case, N-1 is chosen.Before doing any N-1 contingency analysis, it is compulsory to comply the contingency voltage range is set within ±10% of nominal voltage at steady state level for 33 kV distribution network.Since the nominal voltage is 1.0 pu, therefore the range is from 0.9 pu to 1.1pu.Besides, there are two types of contingency case analyzed as follows: 1) External contingency (N-1): is the power network cut off supply from main grid 2) Internal contingency (N-1-1): is another contingency happens after external contingency occurred Branch flows and overload condition are determined for affected branches in every contingency case.Also, determines the buses that violate the contingency voltage deviation criterion which is the changes in voltage cannot rise more than 0.06 pu and drop more than 0.03 pu.
From Table 7, transformer branches which are marked in italic are recorded higher than that in distribution branches.The most severe overload cases happen on the transformer branches near to SB_6.6 and KB_6.6 generators.This is because the lost generation from incomer main grid is supplied by these two supplies, causing them to generate more power demand.
In the contingency analysis, all branches except six branches are considered safe and remained unaffected throughout all contingency cases (N-1-1 and N-1).
Table 8 shows that the list of 31 unaffected branches out of 37 branches involved in contingency analysis.The branches highlighted in italic are transformer branches and the rest are distribution branches.
There is no contingency voltage deviation violation report in all contingency cases except for case when outage of transformer branch bus 8 to bus 28.The bus that cannot withstand N-1 or N-1-1 is if they violate the voltage deviation criterion.The violation case can be shown visually in Figure 3.The affected buses highlighted in red are the buses that cannot withstand N-1.The transformer branch outage is highlighted in black.Table 9 shows the affected buses during the contingency case named SINGLE 8-28 with their respective contingency voltage, initial voltage, and the voltage deviation limits.Contingency results showed 3 transformer branches and 2 distribution branches are affected.

P-V (Load Power-Load Voltage) Curve Analysis
P-V curve is commonly used as voltage stability analysis tool to analyze maximum additional load power that a bus can sustain before its voltage collapses, after the load power exceeds its power limits.Figure 4 shows a simple P-V curve.When load power exceeds the power limit of a bus, the bus voltage will start to drop until reaching the critical point, which is voltage collapse point when the bus reaches the maximum additional power.PSS/E is able to simulate different P-V curve for different level of contingency.So, in this paper, P-V curve is simulated for external contingency and internal contingency cases for all 26 in-service buses.The maximum additional power transfer and voltage collapse point are recorded for each buses in every contingency cases.
All buses are simulated with P-V curves for every contingency case.Thus, 988 (38 contingency cases X 26 buses) P-V curves are simulated.Table 10 shows the list of number of contingency cases with the overall maximum additional power transfer among all buses.Among 38 contingency cases, 15 of them are recorded with the highest maximum additional power transfer with 312.50 MW, followed Among all P-V curves, almost all buses except generator buses have their voltage collapse point at below 0.9 pu.That is not reasonable that the steady state operating voltage is within 10% of the nominal voltage, from 0.9 pu to 1.1 pu.
Since the weakest case is SINGLE 8-28 (1), in order to determine weakest buses, the maximum incremental power transfer at 0.9 pu of the nominal voltage is determined for each bus.Any voltage which is outside of the tolerance is not acceptable.Figure 5 shows the finding of maximum additional power transfer at Figure 5. Maximum additional power transfer at 0.9 pu in PV curve.
0.9 pu.Table 11 shows the maximum additional power transfer at 0.9 pu for each bus with respective voltage collapse point at contingency case named SINGLE 8-28.
From Table 11, the generator buses (KB_6.6,SB_6.6,BN_11, GN_11, LD1_11, LD2_11, LD3_11 and LD4_11) have maintained constant bus voltage along all incremental power transfer.This is because no losses involved when the generators generate the output voltage to the nearest buses.The loss will be very small.Any bus that withstand below 50% of the 143.75 MW is considered as weakest buses.Therefore, the weakest buses are SB_33, KB_33, SD_33, TS_33, Table 11.List of all buses with their respective maximum additional power transfer at 0.9 puin weakest contingency case (Single 8-28).

Figure 1 .
Figure 1.Existing power network drawn in PSS/E.
generation which takes place on i th bus For LODF, the line outage distribution factor is defined as the change in power flow on i th line during pre-contingency line flow on  th line.Line outage distribution factors, power flow on  th line o i P = Power line flow on i th line before contingency of i th line

Figure 2 .
Figure 2. Process of creating these files.

Figure 3 .
Figure 3. N-1-1 case with affected buses when outage of transformer branch from bus 8 to bus 28.

Table 2 .
Total generated power in sandakan before and after N-1 outage.

Table 4
and Table5show the load flow results for all in-service buses.Load buses are represented by bus code 1 while generator buses are represented by bus code 2. UB_33 and BM_33 are recoded zero in active and reactive power flows and current flows due to no loads connected to them.Generator buses are recorded higher current flows compared to load buses.All bus voltages are within contingency voltage range stated by Distribution Code of Energy Commission.

Table 3 .
Distribution and transformer branch results.

Table 4 .
Load flows results for all in-service buses.

Table 5 .
Load flow results for all in-service loads.

Table 6 .
Three types of PSS/E files with descriptions.

Table 7 .
List of affected branches with respective total number of involved cases and overload percent.

Table 9 .
Contingency voltage deviation violation reports. .H. Teck, N. Barsoum by 6 cases with 306.25 MW and 5 cases with 287.50 MW.The only one case with the lowest maximum additional power transfer (143.75MW) happened when outage of transformer branch from bus 8 to bus 28.That means all buses can only withstand incremental power transfer up to 143.75 MW when this outage happens.That is why this case has only voltage violation cases, referring to Figure 4. PV curve.K

Table 8 .
They cannot withstand higher power flows compared to other cases.

Table 10 .
List of number of contingency cases and the maximum additional power transfer.

Table A2 .
Transformer branch data used in PSS/E.

Table A3 .
Machines data used in PSS/E.

Table A4 .
Load data used in PSS/E.

Table A5 .
Branch/distribution line data used in PSS/E.