Real-Time State Estimation in the Cameroon Power System

An important tool for the energy management system (EMS) is state estimation. Based on measurements taken throughout the network, state estimation gives an estimation of the state variables of the power system while checking that these estimates are consistent with the measurements. Currently, in the Cameroon power system, state estimates have been provided by ad hoc supervisory control and data acquisition (SCADA) systems. A disadvantage is that the measurements are not synchronised, which means that state estimation is not very precise during dynamic phenomena in the network. In this paper, real-time phasor measurement units (PMUs) that provide synchronised phasor measurements are proposed for integration into the power system. This approach addresses two important issues associated with the power system state estimation, namely, that of measurement accuracy and that of optimization of the number of measurement sites, their location, and the importance given to their measurements on the dynamic state estimation.


Introduction
The electricity sector in Cameroon has been undergoing gradual deregulation in the last few years.AES-SONEL, which took over the generation, transmission, and distribution of electricity in the country in 2001, disposes of a total installed capacity of 933 MW (721 MW from hydroelectricity sources and 212 from thermal sources including 24 MW of isolated capacity) [1].The transportation network has 480 km of 225 kV lines, 337 km of 110 kV lines, and 1064 km of 90 kV lines, with 24 sub-stations.
The need for better utilization of power system facilities has increased significantly in the last decade and will continue, as the electricity market is becoming more complex and competitive.More real-time data from substations and generating plants is required to support SCADA (supervisory control and data acquisition) and EMS (energy management system) in utility control centres.There has been a significant increase in interest in planning, operating, equipment, and maintenance personnel for real-time events and archived data.This wide interest demonstrates the need for the use of real-time state estimators.
The system control centre in Mangombe is responsible for the load dispatch and reliability of the main 225 kV grid.Tanyi and Mbinkar have proposed that four area control centres can be located throughout the regions covering the southern interconnected grid (SIG) to be responsible for the operation and maintenance procedures on the transmission, subtransmission, and distribution networks in their respective areas [2].There are presently a few telemetry systems, and all generating plants and major substations have separate remote terminal units (RTUs) and/or remote supervisory and telemetry systems reporting to both the system and area control centres.
Almost all the RTUs and older remote supervisory and telemetry systems connect to the control centres over dedicated telecommunication circuits on a one-on-one basis usually over telecommunication media owned by Cameroon Telecommunications (CAMTEL).CAMTEL has an extensive microwave telecommunication network for protective relaying, remote supervisory control, and operational voice requirements.There are also many power line carrier telecommunication links.Both the microwave and the power line carrier equipment are moving towards digital technology with the deregulation process.Some stations are not accessible by the telephone common carrier, and most of these use point-to-multipoint cellular radio communications.
State estimation forms an important part of the modern day control centre and substation automation functions.State estimation calculates the voltage magnitudes and angles at all the buses of the power system.The accuracy of the state estimation output is very crucial, as it forms the basis for the EMS functions like security analysis, voltage stability analysis, optimal load dispatch, and so forth.The power system is not a static system, but it changes with time, though very slowly.This is essentially due to the continuous variation of loads.Once the loads change, the generation has to keep up to those changes, and this triggers a change in the generations.This in turn leads to a change in flows and injections at all the buses of a power system.These dynamic changes in the power system are difficult to be captured by the traditional static state estimation (SSE) in place.
The accuracy of the state estimation depends upon the quality of the communication links and the measurements used for the state estimation.Traditionally, line flows, injections, and voltage measurements were the most commonly used measurements for state estimation.In the 1980s, a new measuring instrument called the phasor measurement unit (PMU) was developed [3].The specificity of the device is that it can directly measure the voltage phasors and current phasors at a given bus.In other words, it can directly measure voltage magnitude and its angle, which forms the state vector.Hence, if all the buses were to have PMUs installed and their measurements were transmitted to the control centre, we can have a complete real-time measurement of the state vector of the power system, instead of estimation.But due to cost considerations, PMUs cannot be installed at all the locations of a power system.Hence, the state estimators need to process the normal SCADA data along with the PMU data to obtain a good estimate of the states of the power system at a given instance.
In this paper, a measurement placement algorithm for the PMUs is proposed.The proposed algorithm is applied to find the positions of measurement placement in the southern interconnected grid (SIG) of the Cameroon power system.The number of PMU measurements is also very important for the accuracy of the state estimation.The results obtained are based on simulations carried out using an evaluation version of PowerWorld software.

Mathematical Formulation
where  is a ( * 1) measurement vector obtained from PMUs,  is a ( * 1) true state vector, ℎ(⋅) is a ( * 1) vector of nonlinear functions,  is a ( * 1) measurement error vector,  is the number of measurements, and  is the number of state variables.The usual choice for state variables is the bus voltage phase angles and magnitudes, while the measurements are active and reactive power flows, node injections and voltage magnitudes.
As in load flow calculations, it has been found that state estimation algorithms based on decoupled versions of (1) behave adequately for the usual power networks [4].Consequently, the following decoupled model is usually adopted: where  (  * 1) and V ( V * 1) are the vectors of true voltage magnitudes and phase angles and  and  are subscripts indicating partitions of vectors and matrices corresponding to active and reactive measurements, respectively;   = −1,  V = , and  is the number of network nodes.

Problem Formulation.
The state estimation problem of a power system can be considered as the following weighted least squares (WLS) method optimization problem [5]: where [] is called a measurement covariance matrix.Usually, the distribution of each measurement error is assumed to be Gaussian noise with zero mean and independent; thus, [] becomes a diagonal matrix [6].These weights are chosen as proportional to the accuracy of measurement.Normally, the covariance of the measurements is unknown and often assumed to be an identity matrix since the same instrumentation is used to obtain them.The solution of (3), in WLS sense, is obtained by In this paper, the voltage phasors at all busbars are chosen as state variables, since they allow the branch currents, shunt currents, and currents injected into the busbar to be determined [5].The measured values are the busbar voltage phasors, the injection current phasors, and the line current phasors.Based on these measured values, all possible relations of the estimated state variable vector V and the measurement vector  in the presence of measurement error w can be expressed as where  Equation ( 5) indicates all possible relationships between the measured values and the estimated state variables [5].In practice, the measuring devices will not be placed at all possible positions for economic reasons.Thus, it is necessary to find an appropriate measurement placement position set such that all state variables can still be estimated with a minimum number of measuring devices.

Operating Principles of a PMU
As shown in Figure 1, the analogue power signal is converted into digital data by the analog to digital converter.For example, if we want to measure the voltage of the signal, then samples are taken for each cycle of the waveform, and then the fundamental frequency component is calculated as follows using discrete Fourier transform (DFT) [7]: where  is the total number of samples in one period,  is the phasor, and   are the waveform samples.
From the above expression, we can see that the phasor information is dependent on the number of samples of the waveform and will give a correct representation of the fundamental frequency component even in the presence of transient signals.The samples are taken over one cycle or over multiple cycles of the power signal.
Data windowing is used to smoothen the frequency spectrum and minimise overshoot of the signal caused due to discontinuities of the signal in the time domain [7].
The optical isolator shields the PMU system from the very high bus voltages.The unit is composed of a number of phasors that capture measurements of analogue voltage, current waveform, and the line frequency.After that, the phasor measurements are digitized by an analogue to digital converter and stamped with the creation time provided by a clock at the network control centre (NCC).The clock pulses are typically supplied by GPS, but optical fibre cables have proven to be better [5].These pulses are used for synchronization of multiple PMUs with a precision of maximum 1 microsecond difference.Afterwards, the data is transferred to a phasor data concentrator, where data with the same creation time are encapsulated in one packet and then transmitted as a single stream to the phasor data server which can be a PSGuard [2] or traditional SCADA.PMU provides a dynamic system observation of the network, because the measurements are taken with a high sampling rate from geographically distant locations, and they are then grouped together according to the time stamp provided.Phasor measurement unit transmits samples in different sizes.The sample size depends on the number of phasors in a unit.Samples and packets have the same meaning when talking about PMU transfer rate.The required transfer rate differs from a 50 Hz system to a 60 Hz system.For example, a 60 Hz system has a rate of up to 60 samples per second, while the 50 Hz one has up to 50 samples per second [8].
The PMUs measure voltage, current, and frequency phasors using the discrete Fourier transform (DFT) and can detect transients or surges within milliseconds of their occurrence [9].
PMUs use the IEEE 1344 data format for communication with the central monitoring station.In addition to the propagation delay of the particular link, the message format of the PMU and data rate of the link determine the communication delay in the system.Furthermore, there is also a processing delay due primarily to the window size of the DFT.As shown in Figure 2, PSGuard system utilizes phasors which are measured by PMUs resulting in an accuracy of <1 s [9].A communication system links PMUs, which are located in substations to the PSGuard system for data transmission.Data exchange can be established between PSGuard and other control and protection systems to allow for optimum data sharing and control actions.This calls for an efficient communication network for the SIG.There are already several projects in progress for the installation of fibre optic cables [10].

Measurement Placement Algorithm
The proposed measurement placement algorithm for the PMUs was described using three IEEE test systems [5] and is now being applied to the SIG.It is assumed that the power network is decomposed into  areas [4].The areas are connected through boundary buses which belong simultaneously to both adjacent areas.

Decomposition Technique.
It is assumed that the entire network under consideration is decomposed into  nonoverlapping observable areas.Each subnetwork consists of internal buses, which are adjacent to branches belonging only to this subnetwork, and boundary buses.These buses are connected to other areas by tie lines.For example, decomposition of a network into four nonoverlapping areas is shown in Figure 3.To obtain decoupled state estimation equations of these areas, constraints of measurement placement are as follows.
(i) There is no injection current measurement at the boundary bus.
(ii) There is no line current measurement on the tie line.
(iii) Each subnetwork has at least one busbar voltage measuring device.
Based on the previous conditions, parameters of the tielines do not appear in the estimation equations, and the measurement model of the power system state estimation can be written as  where   , V  , and ℎ  denote the measurement vector, the state variable vector, and the measurement matrix of the th subnetwork, respectively.  represents the measurement error vector related to   .These   , V  , ℎ  , and   are formulated according to (5) for the th subnetwork.
The state variables of each sub-network can be estimated independently.Since the size of the sub-network is smaller, determining an optimal measurement placement can be performed with less computation effort using either the minimum condition number criteria [11] or other approaches.The optimal measurement placement of each sub-network is selected independently to ensure that the number of measurements is equal to the number of state variables.Thus, the measurement placement obtained from the minimum condition number criteria can be rearranged to minimize the number of placement sites by the following heuristic algorithm [11].

Reduction of Placement Sites
Step 1.Based on the placement position list obtained from the measurement placement algorithm, determine the bus where either an injection current or a busbar voltage measuring Determine the minor bus with the maximum connection number of those branches which have line current measurements.
Determine if there is a line current measuring device on the branch connected to the major buses.
Install injection current or busbar voltage measurement devices while observing the constraints (i), (ii), and (iii) above.
Move line current measuring device on the connected branches close to the selected bus.device is installed.These buses are called major buses.Other buses are called minor buses.

Maximum connection number of line current measuring branches equals to one?
Step 2. If there is a line current measuring device on the branch connected to the major buses, the device is moved close to the major buses.
Step 3. From the branches with line current measuring devices which are not connected to major buses, determine the minor bus with the maximum connection number of those branches which have line current measurements.Then, the line current measuring device on the connected branches is moved close to the selected bus.Note that this minor bus will not be considered again in the next iteration.
Step 4. Repeat Step 3 until the maximum connection number of line current measuring branches equals to one.
Figure 4 illustrates the algorithm for the reduction of placements sites with a flowchart.

Illustration and Results
The flowchart presented in Figure 4 is illustrated using the six-bus power system as shown in Figure 5, in which there are five placement sites and two generators [5], and then applied to the high voltage network (225 kV/90 kV) of the SIG.

Six-Bus Power
System.The six-bus power system with five PMU placement sites and two generators is shown in Figure 5.
Step 1. Determine the major buses and minor buses (see Table 1).
Step 2. Move devices close to major buses (see Table 2).
Step 3. We can see that bus number 2 has the maximum connection number of three, since it is connected to buses 3, 4, and 6.Thus, the bus no. 2 will be chosen in this iteration.
Then, the device is moved close to bus no. 2.
Step 4. After moving the device, buses 5 and 6 have the connection number of one.Thus, the iteration procedure is stopped, and the placement positions of the six-bus power system become as shown in Figure 6.Notice that the number of placement sites reduces to three (see Table 3).

Application to the SIG.
The SIG has two main generating stations located at Edea and Songloulou which are interconnected to the high voltage network using 225 kV and 90 kV transmission lines.There are five main substations-Bekoko, Logbaba, Mangombe, Oyomabang, and Nkongsamba [2].
Figure 7 shows the interconnection of the busbars of these substations with tie lines.The generator at Bekoko represents  the Limbe thermal plant.The entire network is divided into three areas.
The load dispatch centre is located at Mangombe and is responsible for monitoring and controlling the entire electrical network through regional control centres located at the other substations.The communication network between the different substations consists of optical fibre cables.The algorithm previously described is applied to the SIG to determine the optimal number and placement sites for the PMUs.Step 1. Determine the major buses and minor buses (see Table 4).
Step 2. Move devices close to major buses (see Table 5).
Step 3. It is seen that bus no. 1 has the maximum connection number of two, since it is connected to buses 3 and 4.But the devices are already moved to bus no. 1.So, the algorithm ends.The number of placement sites is consequently reduced from seven to four.These are located at Songloulou, Edea, Bekoko, and Nkongsamba.

Real-time Simulation of the SIG.
A load flow analysis of the reduced network with an evaluation version of Power-World is seen in picture form in Figure 9.What stands out is that the line from Bekoko to Nkongsamba is heavily loaded beyond its dynamic and thermal limits.This makes the SIG unstable and unreliable, and the transmission losses must be high.There are parallel lines between the generating points and the major load centres in Douala and Yaounde.However, these parallel lines are of two different transmission voltages (90 and 225 kV) and capacities.The 90 kV line is unable to carry the whole load especially to Yaounde should the 225 kV line fail.This is what makes the SIG very unreliable.The situation is better (but perhaps only marginally so) for Douala because although there are lines to Logbaba (a 225 kV line from Songloulou, a 225 kV line from Mangombe, and a 90 kV line from Edea), the components are old, interconnection transformers are old and/or overloaded, and one of the two   225 kV lines is overloaded.The result of this simulation stresses the need for transmission line upgrade and renewal within the SIG [12] in order to enhance its reliability.

Discussion
Conventionally, we want the PMU placement to be optimal in the sense that it makes all buses in the system observable with a minimum number of PMUs.This is especially true for largescale complex systems [13].The proposed algorithm is applied to the six-bus test systems under balanced conditions in which single-phase representation is adequate.This method is referred to as Madtharad's algorithm [5].The true state of the network is computed by using power-flow calculation.
When this algorithm is applied to the SIG, starting from seven initial installation sites, it is seen that the number of placement sites has been reduced to four.The basic rule of PMU placement is that, when a PMU is placed at a bus, it can measure the voltage phasor at that bus, as well as at the buses at the other end of all the incident lines, using the measured current phasor and the known line parameters [14].It is assumed in this study that the PMU has a sufficient number of channels to measure the current phasors through all branches incident to the bus at which it is placed.Figure 8 illustrates the optimal placement of the PMUs within the SIG.
Other researchers, notably [3,5,15], have shown that location of the PMU is important to obtain a better state estimate for a power system.Though a PMU at any location provides more accurate estimates than the case when there is no PMU, there are certain PMU locations in the network which seem to provide much higher reduction in error than at any other bus in the network.This PMU placement algorithm can therefore be followed to quickly find out the optimal location of PMUs, especially in the case of large networks.

Conclusion
In this paper, the application algorithm for PMU measurement placement for power system state estimation is presented.The proposed algorithm finds an optimal measurement placement based on minimum condition number criterion of each subnetwork.Since the problem size is smaller, this decomposition improves the computational time required to obtain types and positions of the available measurement compared to the previously applied methods based on SCADA which consider the whole network.From the results obtained, the decomposition technique provides a promising strategy to solve the measurement placement problem of a large power system state estimation.In all cases considered [16], PMUs have a beneficial effect on the accuracy of the state estimation.Further studies will consider the dynamic model of the network.
The main contribution of this work lies in investigating the feasibility of using Madtharad's algorithm for the PMU placement problem in the Cameroon power system.Future work will include additional constraints into the PMU placement problem, such as the existence of conventional measurements, user-defined measurement redundancy at the buses, and the consideration of measurement uncertainty.These constraints are difficult to handle by conventional optimization methods.The promising results presented in this paper will encourage the researchers and utilities in using Madtharad's algorithm for the larger power systems.

Glossary
AES SONEL: The company that manages the transmission and distribution of electricity in Cameroon Bekoko: Area within the city of Douala with the second largest electricity substation Busbar: An electrical conductor onto which electric power is injected into the network or removed for distribution Cameroon: Country in west central Africa CAMTEL: Cameroon Telecommunications PLC Douala: Cameroon's economic capital Edea: City in the littoral region of Cameroon where about 40% of the country's electricity is generated EMS: Energy management system Limbe: City in the south-west region of Cameroon where the country's largest thermal (heavy fuel oil) electricity plant is located Logbaba: Area within the city of Douala with the largest electricity substation Mangombe: City in the littoral region of Cameroon where the electricity network load dispatch centre is located NCC: Network control centre Nkongsamba: City in the littoral region of Cameroon with a substation connecting the longest transmission lines Oyomabang: Area within the city of Yaounde with the largest electricity substation, and some thermal plants PMU: Phasor measurement unit PowerWorld: Software used in the modeling and analysis of electric power systems PSGuard: Power system wide area monitoring system based on ABB's process control system and Windows Server technology Real-time: Being able to process data as it comes in SCADA: Supervisory control and data acquisition SIG: Southern interconnected grid of the Cameroon power system Songloulou: Area in the littoral region of Cameroon where about 45% of the country's electricity is generated Substation: A branch of a main electrical power station where electric power is converted, redistributed, or modified in voltage level Tie-line: Electricity transmission line connecting two adjacent busbars (buses) or substations Yaounde: Cameroon's political capital.

Figure 1 :
Figure 1: Basic block diagram of a PMU.

Figure 3 :
Figure 3: Decomposition of a network into four nonoverlapping sub-networks.
measuring device closer to the major bus Divide the entire network into r nonoverlapping areas.

Figure 4 :
Figure 4: Flowchart for reduction of PMU placement sites.

Figure 5 :
Figure 5: Six-bus power system with five PMU placement sites.

Figure 6 :
Figure 6: Optimal placement sites for PMUs in the six-bus power system.

Figure 7 :
Figure 7: Installation of measuring devices within the SIG.

Figure 8 :
Figure 8: Optimal placement sites for PMUs within the SIG.

Figure 9 :
Figure 9: Simplified real-time model of the SIG.
[]is an identity matrix,   ,   , and   are the measured values of the busbar voltage, the injection current, and the line current subvectors, respectively,   and   are the bus admittance submatrices of the measured and the calculated busbar voltage related to , respectively,   and   are the bus admittance submatrices of the measured and the calculated busbar voltage related to   , respectively, V  and V  are the measured and the calculated busbar voltage subvectors, respectively, and   ,   , and   are the corresponding measurement error or noise subvectors.