Professor Subhas Chandra Mukhopadhyay
Exeley Inc. (New York)
Subject: Computational Science & Engineering, Engineering, Electrical & Electronic
eISSN: 1178-5608
SEARCH WITHIN CONTENT
Husham Idan Hussein ^{*} / Ghassan Abdullah Salman / Ahmed Majeed Ghadban
Keywords : RDPS, RDPS, PSO, Neural network, Hybrid algorithm, State estimation, Monitoring system, Optimization methods, PMU
Citation Information : International Journal on Smart Sensing and Intelligent Systems. Volume 12, Issue 1, Pages 1-10, DOI: https://doi.org/10.21307/ijssis-2019-005
License : (BY-NC-ND-4.0)
Received Date : 17-January-2019 / Published Online: 05-September-2019
This study presented a new hybrid algorithm to improve the state estimation (SE) of radial distribution power systems (PSs). The proposed particle swarm optimization–neural network (PSO–NN) algorithm constructed an independent and fast monitoring system with high accuracy that can detect abnormal conditions or failures in a PS. In this study, PSO was adopted to determine the appropriate weights of the NN model. The speed and accuracy of PSO with the NN model were evaluated in the SE of power system networks (PSNs). The information obtained through SE was used to enhance the operations and customer service delivery in terms of energy consumption and power quality in PSNs. Capacitor banks were installed to reduce the losses and improve the voltage profiles. The PSO–NN algorithm was assessed on IEEE (9, 33, and 69) bus standards. Simulation results proved that the new technique can be tested on any distribution network because of its accurate and efficient SE. Results indicated that the PSO–NN algorithm had better performance than the phasor measurement units.
SE
state estimation
PS
power system
PSN
power system networks
PSO–NN
particle swarm optimization–neural network
DPN
distribution power network
RDPSNs
radial distribution power system networks
DMS
distribution management system
GTP
graph theoretic procedure algorithm
DSSE
distribution system state estimation model
WLS
weighted least square
WLAV
weighted least absolute value
PMUs
phasor measurement units
ANN
artificial neural network
EMS
energy management system
Energy management systems (EMSs) have become obsolete as control frameworks due to the expansion of interconnected electric power frameworks and the multifaceted nature of intensity framework structures (Lu et al., 1993; Gastoni et al., 2004; Pajic and Clements, 2005). State estimation (SE) performs an essential job in EMSs and provides a dependable and predictable framework. This condition requires the improvement of power system (PS) operations. Energy transmission from a small number of nodes (substations) to a multitude of consumers is performed by the electric distribution system. Future system networks of electric power distribution will be ‘smart,’ more efficient, better secured, and capable of distributing more reliable and affordable electric power than today’s system network (Smart Grid, 2008; Husham and Ahmed, 2018). The magnitudes and angles of branch current or system node voltages often define the SE of a distribution PS. A system of nonlinear equations is formed with the available information about the reactive and real power and their states. PSSE is defined as ‘a data processing algorithm that converts redundant meter readings and other available information into an SE of an electric PS.’ As shown in Figure 1, ample measurement redundancies in the transmission network has allowed system observation and processing of bad data in SE after more than four decades of development.
A review of important algorithms and DSSE requirements in distribution management systems is provided in the study of Jako (2009). DSSE’s future as part of the vision for a smart grid and as an emerging technique is identified. Two SE techniques that use compressed power measurements are presented in the study of Anggoro and Chan (2017). The primary computational complexities of the two methods are discussed. Direct and indirect SE approaches can accurately estimate voltage states by randomly projecting only half of the available power information, and they have similar performance with the same complexity order. WLS and WLAV methods of power SE are described in the study of Shafiul et al. (2014). The proposed estimation method can still estimate the PS state because it has fewer variables than WLAV and a wider performance range than WLS although the system is unobservable. A complete formula is provided for single and three-phase DSSE optimization problems on the basis of the system voltage drop formula and system quasi-symmetric matrix, as presented in the study of Syed et al. (2017). DSSE is suitable for large-scale radial-exploited distribution feeders by removing interconnecting nodes with null-injected reactive and active powers, as proposed in the study of Hornik et al. (1989).
The artificial neural network (ANN) is fed with adopted phase angles and voltage magnitudes as inputs. This process greatly improves the precision of load active power margin estimation for the New England 39-bus system. Phasor measurement units (PMUs) proposed in the study of Seyed and Mohammad (2014) can provide phase angles and voltage magnitudes for real-time applications. In addition to normal steady-state conditions, the ANN-based system can estimate voltage stability margins under different possibilities. To achieve a good ANN performance, the ANN algorithm is used in this work as a mapping tool for estimating the available voltages and margin angles of the system. Feed forward networks can approximate any measurable function to the desired accuracy level (Chakrabarti and Jeyasurya, 2004; Zhou et al., 2010; Paulo et al., 2013; Paulo et al., 2015; Ghassan et al., 2018). This advantage is tackled in this study. The proposed ANN approach is compact and efficient because it is fed with the input features of the post-power flow. Installing PMUs at every node and connecting them to fast and reliable telecommunication channels to gather the required inputs are inefficient. A limited number of PMUs should be placed at important locations. Thus, a suboptimal search method is proposed for the placement of PMUs. The placement of PMUs for detecting bad data is proposed in the studies of Popovic et al. (1998) and Husham (2016). This study indicates that system observability is improved by critically eradicating measurements when the number of strategically placed PMUs is minimum. This process enhances bad data detection. The concept of robust distances was presented by Chen and Abur (2008) to detect the difference between bad and good leverage measurements. They posited that a robust estimator performs better than a WLS estimator with the presence of bad leverage data point and bad data.
This study presented a new hybrid algorithm that uses PMUs to evaluate the SE of radial distribution in PSs. The proposed hybrid algorithm is a combination of particle swarm optimization (PSO) and NN model in which PSO was used to determine the appropriate weights of the NN model. The speed and accuracy of PSO and ANN were evaluated in the SE of power system networks (PSNs). The PSO–NN algorithm was tested on three radial distribution PSNs (RDPSNs) of IEEE (9, 33, and 69) buses. The networks were designed using MATLAB/PSAT package. The hybrid algorithm code was developed using a MATLAB program. Capacitor banks were installed as VAR compensators to reduce the power losses and improve the voltage profiles of RDPSNs.
Conventional SE is the same as the objective function of distribution SE, which is expressed as follows:
Where:
The difference between the measured and calculated state variables is minimized by using this function. One of the state variables is used as a load value at each section rather than a voltage utilized by conventional state estimators. Therefore, the reactive and active state power values are used as state variables. The input values V and δ are also used as state variables. In the study of Narvaez and Grijalva (2008), the state variables are obtained among the boundaries.
The constraints of the objective function are defined as follows:
The ratio of the target section to the total load of the target network and the total power inputted into the target network are used to calculate the center value of the bound at each load. The average input is the center value of the bound of each variable. Fast power flow distribution is utilized to calculate the current and voltage (Shigenori et al., 2001).
A PMU can measure the current phasors of some or all of the lines connected to the installed bus and the voltage phasors of the installed bus. Efficient placement strategies (benefits derived from PMUs) should be applied to justify the costs of bringing the measurement system up to date. PMU placement methods are facing a huge challenge from the increasing dependence on renewable energy sources that are constantly added to modern distribution grids (Fukuyama et al., 1996).
A distribution grid is converted to an active system and its power path flow is changed by distribution generators. Literature review of papers (Manousakis et al., 2011; Ren and Jordan, 2018) published in the past addressed PMU placements in their classification of PMU placement problems on the basis of the algorithms and methods utilized to determine the difficulty of selecting the location and ideal number of installed PMUs. Graph theoretic procedure (GTP), integer programming, PSO, and generic algorithm methods have different constraints and different requirements for PMU installation objectives, which can be used to achieve different goals in different cases.
Rules that can be applied to PMU systems to make them observable are as follows:
Rule 1: the current phasor and bus voltage of all incident branches of installed PMU at a bus are all known.
Rule 2: the voltage phasor at the other end of one end of a branch can be obtained using Ohm’s law when the current and voltage phasor at the other end of the branch are known.
Rule 3: the current phasor through a branch can be calculated when the voltage phasors at two ends of the branch are known, as shown in Figure 2. In this study, a GTP algorithm was utilized to determine the location of PMUs. This process was conducted to show the observability on the basis of the topology and to ensure full observability on the optimal placement of PMUs.
A simple approach for determining the locations and sizes of capacitor banks in RD feeders was presented in this study. In this approach, the total reactive loads and total reactive losses in all sections of the feeders were the same as the amount of required compensation for any feeder (Manousakis et al., 2012).
We installed capacitor banks on a bus through trial and error method to minimize power losses and improve voltage profiles. In the distribution network, the compensation of reactive power is a vital issue. Transformers, converters, and induction motors are combined to consume reactive power (More and Jadhav, 2013). A reasonable amount of reactive power is not transmitted through the network because the voltage drop and energy losses are increased. Thus, reactive power should be generated close to its consumers (Huang et al., 2014; Mohamed et al.). Figure 3 illustrates the effect of reactive power compensation.
Compensating devise Q _{ c } is connected to the consumer terminals, and the power delivered through a distribution line is as follows:
where S _{ L } is the power delivered; P the real power; Q and Q _{ c } the reactive power and reactive power compensator.Each NN in PSO–NN consists of its velocity and network position. The weight of the NN (W) is related to the network position. The network’s velocity (∆ W) is equivalent to updating the NN’s weight. PSO’s main role in the NN is to attain the best set of weights (particle positions) where numerous particles are aiming for the best solution.
A PS is composed of numerous particles. Each particle monitors its own position, best position, velocity, best fitness, current fitness, and the particles nearby. To operate the NN successfully, the position vector should match the weight vector of the network. The fitness value should match the forward transmission through the network. The particle’s new solutions are guided by using its best neighbor and global best particle. The particle’s position serves as the answer at the end (Hamada et al., 2008). The PSO’s main task in the NN is to secure the best set of weights in a situation where numerous NNs compete for the best solutions. The velocity vector is given by the following equation:
In each iteration the position of each particle is updated by using the velocity vector of the following equation:
where ${X}_{i}^{k}$ is the particle position at (k) iteration; ${X}_{i}^{k+1}$ the particle position at (k + 1) iteration; ${V}_{i}^{k}$ the particle velocity of at (k) iteration; ${V}_{i}^{k+1}$ the particle velocity at (k + 1) iteration; w the weight parameter; r _{1} and r _{2} the random number; c _{1} and c _{2} the learning factors [0, 1].The development of the PSO–NN model with its angle and voltage magnitude is expressed as follows:
A population of NNs with various initial weights was assembled for the PSO–NN.
For every NN (particle), an iteration of the training data set was conducted, and the sum of squared errors was calculated.
All the errors of NNs were compared to find the best NN in the neighborhood.
The weights of the network were recorded, and Step 6 was performed when one of the networks exceeded the required minimum errors. Otherwise, Step 5 was performed.
The PSO algorithm was applied for each NN to update its position, that is, its velocity vector (Δ W) and weight (w).
The set error goal was obtained from the trained PSO–NN.
The testing pattern was applied to evaluate the performance of the trained PSO–NN.
The PSO–ANN model aimed to estimate bus voltage and bus voltage angles.
A distribution system is characterized with short lines that have high R/X ratios and systems that are unbalanced. This system experiences more losses than a transmission network. An RDS normally spreads over a large area (Husham et al., 2018). This condition makes many customers’ domestic loads difficult to monitor and measure. The network’s observability and state estimator’s capability to provide realistic estimates are also hampered. This serious challenge can be tackled by using PMUs. This study compared the results of the proposed hybrid algorithm (PSO–NN) with the results of PMUs.
In this study, a new algorithm was tested on PSNs as an observational system to determine the increase and decrease of bus voltage and their angles. Some differences were observed between the values because the networks were monitored rather than performing aggregate analysis, as shown in Figures 7, 8, and 10.
The proposed PSO–NN algorithm was designed on MATLAB. RPDSN IEEE (9, 33, and 69) bus standards were developed and simulated using a MATLAB/PSAT toolbox. PMUs were applied on the distribution networks by using the toolbox. Then, a GTP algorithm was adopted to find the PMU’s location. The capacitor bank’s size and location were determined through trial and error method. For each test, the number of capacitors and PMUs varied. As shown in Table 1, the measurements had low MSE.
The proposed PSO–NN algorithm and its installed PMU (Manousakis et al., 2012) were evaluated on three tests.
A PSAT package is used to mount a PMU on the IEEE 9-bus radial distribution network that has been designed. Two capacitor banks with sizes of 0.045 and 0.05 p.u are installed on buses 4 and 7, respectively, to improve bus voltage. Then, four PMUs are fixed on buses 2, 4, 6 and 8.
As shown in Figures 4 and 5, the results of the PSO–NN algorithm reveal that the state estimator has a good validation. As shown in Table 1, the input value is consistent with the output (target) with only a small error margin. R = 0.99996 and R = 0.99991 are the regressions for bus voltage and δ bus voltage, respectively, as illustrated in Figure 6.
The IEEE 33-bus network is designed using the PSAT package. Four capacitor banks, each with a size of 1 p.u, are mounted on buses 18, 20, 24, and 29. In total, 14 PMUs are fitted on buses 1, 11, 13, 15, 17, 19, 21, 24, 28, 3, 30, 32, 6, and 9. The PSO–NN algorithm used obtains a good correspondence between the input and output (target) with only a minor error, as shown in Table 1. Figures 7 and 8 show a full comparison between the actual and estimated values. R = 0.98354 and R = 0.98606 are the regressions for bus voltage and δ bus voltage, respectively, as shown in Figure 9.
Four capacitor banks are installed on buses (24, 45, 49, and 63), each with a size of 0.55 p.u to improve the voltage profiles. In total, 27 PMUs are installed on buses 12, 15, 17, 19, 21, 23, 26, 3, 32, 34, 38, 40, 42, 45, 49, 5, 51, 54, 58, 60, 62, 64, 66, 69, 7, 69, and 9 to ensure network observability. As shown in Figures 10 and 11, similar results were obtained between the input and target when the PSO–NN algorithm was tested in the network. As presented in Table 1, errors appeared in the results. As shown in Figure 12, R = 0.99934 and R = 0.99583 were the regressions for the bus voltage and δ bus voltage, respectively.
As previously mentioned, the simulation results proved that the hybrid PSO–NN algorithm can be tested on any distribution network because of its accurate and efficient SE of the input values and target.
A new algorithm was developed in this study for the SE of angle bus voltage and voltage magnitude of RDPSNs. Voltage was a good predictor of the proximity to voltage collapse, and phase angle was a good predictor of power flow. Capacitor banks were installed to enhance the voltage profiles. PMUs were installed to make the networks observable. Results indicated that redundancy occurred and the accuracy of the estimated states improved when synchronized PMUs were added into the measurement sets. Three distribution networks were used in the SE of IEEE (9, 33, and 69) buses. The simulation results showed that the proposed algorithm was accurate and fast and only required few variables. A small MSE difference was observed between the target output and the output obtained by utilizing PMUs and PSO–NN, as shown in the regression plots for each network. The target output was represented by a dotted line. MATLAB software and PSAT package were used to obtain the results. The SE problem was transformed to a constrained problem and solved using the PSO–NN algorithm. The effectiveness of the proposed method was confirmed by the simulation results.