Abstract
Power system state estimation is the important module in monitoring the power system network. The state estimator provides states of the power system by processing the measurements placed optimally across the power system network. However, temporary un-observability may occur due to unexpected loss of measurement devices or communication links. In the present paper, the reliability of the static state estimation has been tested under three different scenarios, viz. over-determined, critically determined, and under-determined systems. The state estimation problem has been solved using JADE-adaptive differential evolution algorithm utilizing both the conventional and the synchronized phasor measurements. From the test results it has been determined that the proposed technique determines the solution even when the power system is partially observable and also detects the un-observable buses present in the network. The results thus obtained using JADE have been compared with conventional weighted least square and improved self-adaptive particle swarm optimization-based SE techniques. On the bases of various performance indices, the results show the effectiveness, reliability, and accuracy of the proposed algorithm when compared to other state estimation techniques.
Similar content being viewed by others
Explore related subjects
Discover the latest articles, news and stories from top researchers in related subjects.References
Abdallah EN, Ghazala AA, Hanafy N (2005) Power system state estimation using genetic algorithms. In: Proceedings of 2005 middle east power systems conference, pp. 669–676
Abur A, Expósito AG (2004) Power system state estimation: theory and implementation, 1st edn. CRC Press, New York
Aravindhababu P, Neela R (2008) A reliable and fast-decoupled weighted least square state estimation for power systems. Electr Power Compon Syst 36(12):1200–1207
Basetti V, Chandel AK (2015) Hybrid power system state estimation using Taguchi differential evolution algorithm. IET Sci Meas Technol 9(4):449–466
Bi TS, Qin XH, Yang QX (2008) A novel hybrid state estimator for including synchronized phasor measurements. Electr Power Syst Res 78(8): 1343–1352. http://www.sciencedirect.com/science/article/pii/S0378779607002374
Brest J et al (2006) Self-adapting control parameters in differential evolution: a comparative study on numerical benchmark problems. IEEE Trans Evol Comput 10(6):646–657
Brown M et al (2014) Enhanced bad data processing by phasor-aided state estimation. IEEE Trans Power Syst 29(5):2200–2209
Chakhchoukh Y, Ishii H (2016) Enhancing robustness to cyber-attacks in power systems through multiple least trimmed squares state estimations. IEEE Trans Power Syst 31(6):4395–4405
Chakrabarti S et al (2010) Inclusion of PMU current phasor measurements in a power system state estimator. IET Gener Transm Distrib 4(10):1104–1115
Christensen GS, Soliman SA (1989) A new technique for linear static state estimation based on weighted least absolute value approximations. J Optim Theory Appl 61(1):123–136
Contreras-hernandez EJ, Cedefio-maldonado JR (2006) A self-adaptive evolutionary programming approach for power system state estimation. In: Proceedings of 2006 IEEE international midwest symposium on circuits and systems conference, San Juan, pp. 571–575
Das S, Suganthan PN (2011) Differential evolution: a survey of the state-of-the-art. IEEE Trans Evol Comput 15(1):4–31
De Almeida MC, Garcia AV, Asada EN (2012) Regularized least squares power system state estimation. IEEE Trans Power Syst 27(1):290–297
Deb K (2004) Revolutionary optimization by evolutionary principles
del Valle Y et al (2008) Particle swarm optimization: basic concepts, variants and applications in power systems. IEEE Trans Evol Comput 12(2):171–195
Deng R, Xiao G, Lu R (2017) Defending against false data injection attacks on power system state estimation. IEEE Trans Industr Inf 13(1):198–207
Hossam-Eldin AA, Abdallah EN, El-nozahy MS (2009) A modified genetic based technique for solving the power system state estimation problem. Int J World Acad Sci Eng Technol 3(7):307–316
Kavasseri R, Srinivasan SK (2011a) Joint placement of phasor and conventional power flow measurements for fault observability of power systems. IET Gener Transm Distrib 5(10):1019–1024
Kavasseri R, Srinivasan SK (2011b) Joint placement of phasor and power flow measurements for observability of power systems. IET Gener Transm Distrib 26(4):1929–1936
Korres GN, Manousakis NM (2011) State estimation and bad data processing for systems including PMU and SCADA measurements. Electr Power Syst Res 81(7):1514–1524. doi:10.1016/j.epsr.2011.03.013
Korres GN, Manousakis NM (2012) A state estimator including conventional and synchronized phasor measurements. Comput Electr Eng 38(2):294–305
Li, W et al (2011) Review and research trends on state estimation of electrical power systems. In: 2011 Asia-Pacific power and energy engineering conference, pp. 1–4
Liang J, Sankar L, Kosut O (2016) Vulnerability analysis and consequences of false data injection attack on power. IEEE Trans Power Syst 31(5):3864–3872
Liu S, Wang J (2009) An improved self-adaptive particle swarm optimization approach for short-term scheduling of hydro system. In: 2009 International asia conference on informatics in control, automation and robotics, CAR 2009, pp 334–338
Logic N, Heydt GT (2005) An approach to network parameter estimation in power system state estimation. Electr Power Compon Syst 33(11):37–41
Logic N, Kyriakides E, Heydt GT (2007) Lp state estimators for power systems. Electr Power Compon Syst 33(7):37–41
Lu Z et al (2014) Distributed agent-based state estimation considering controlled coordination layer. Int J Electr Power Energy Syst 54:569–575
Madtharad C, Premrudeepreechacharn S, Watson NR (2003) Power system state estimation using singular value decomposition. Electr Power Syst Res 67(2):99–107
Mallick S et al (2013) Optimal static state estimation using improved particle swarm optimization and gravitational search algorithm. Int J Electr Power Energy Syst 52:254–265
Mili L, Phaniraj V, Rousseeuw PJ (1991) Least median of squares estimation in power systems. IEEE Trans Power Syst 6(2):511–523. http://ieeexplore.ieee.org/lpdocs/epic03/wrapper.htm?arnumber=76693
Minot A, Lu YM, Li N (2016) A distributed gauss-newton method for power system state estimation. IEEE Trans Power Syst 31(5):3804–3815
Monticelli A (1999) State estimation in electric power systems—a generalized approach, 1st edn. Kluwer Academic Publishers, Boston
Rakpenthai C et al (2010) An interior point method for WLAV state estimation of power system with UPFCs. Int J Electr Power Energy Syst 32(6):671–677
Safdarian A, Fotuhi-Firuzabad M, Aminifar F (2012) A non-iterative approach for AC state estimation using line flow based model. Int J Electr Power Energy Syst 43(1):1413–1420
Sodhi R, Srivastava SC, Singh SN (2010) Phasor-assisted hybrid state estimator. Electr Power Compon Syst 38(5):533–544
Storn R, Price K (1996) Minimizing the real functions of the ICEC’96 contest by differential evolution. In: Proceedings of 1996 IEEE evolutionary computation conference, pp 842–844
Wang Y, Cai Z, Zhang Q (2011a) Differential evolution with composite trial vector generation strategies and control parameters. IEEE Trans Evol Comput 15(1):55–66
Wang Y, Cai Z, Zhang Q (2011b) Differential evolution with composite trial vector generation strategies and control parameters. IEEE Trans Evol Comput 15(1):55–66
Zhang J, Sanderson AC (2009) JADE: Adaptive differential evolution with optional external archive. IEEE Trans Evol Comput 13(5):945–958
Zhao J et al (2016) Power system real-time monitoring by using PMU-based robust state estimation method. IEEE Trans Smart Grid 17(1):300–309
Zimmerman RD, Murillo Sánchez CE, Thomas RJ (2011) MATPOWER: Steady-state operations, planning, and analysis tools for power systems research and education. IEEE Trans Power Syst 26(1):12–19
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflicts of interest
The authors declare that there is no conflict of interest.
Human participants or Animals performed
This article does not contain any studies with human participants or animals performed by any of the authors.
Additional information
Communicated by V. Loia.
Appendix A
Appendix A
Differential evolution is a stochastic optimization technique, and the solution obtained is dependent on the set of variables generated randomly by searching over different promising areas (Das and Suganthan 2011; Brest et al. 2006; Wang et al. 2011a). The process of estimating the states of un-observable buses has been explained with an example of IEEE 14-bus system as given below:
JADE uses a population P of size \(N_P \) that consists of floating point encrypted chromosomes. Each chromosome \(p_j \) is a vector that contains D decision variables. Similar to all other evolutionary algorithms (EAs), JADE begins with randomly generating chromosomes of size \(N_P \) with D decision parameters within the promising area constrained by the specified minimum and maximum parameters limits. For example, let us consider the initial value of the jth parameter in the ith chromosome is generated according to (11). Here the size of decision parameters for IEEE 14-bus system is 27 (14 voltage magnitudes and 13 voltage angles) with boundary constraints given in (4) (Das and Suganthan 2011; Wang et al. 2011a). Once the initial population is generated, the fitness of each individual chromosome is calculated according to (2). This population evolves over G generations to attain an optimal solution. Further, a parent population from the current generation is termed as a target vector.
where \(p_j^{\min } \hbox { and }p_j^{\max }\) are minimum and maximum values of state variable \(p_j \), and \(\alpha _j \) is uniformly distributed random number in the range of 0–1.
After initialization, JADE is carried out with three simple cycles of stages, viz. differential vector-based mutation, crossover, and selection.
1.1 Difference vector mutation
Another population known as donor vector \(v_{i,G} \) population is obtained by performing difference vector mutation on each individual (i) of the current population \(\{x_{i,G} |i=1,2,\ldots ,N_P \}\). In order to perform difference vector mutation, for each chromosome in the current population three other chromosomes are selected. These three chromosomes are selected as follows. First to improve the convergence characteristic, best chromosome \(x_\mathrm{best,G}^p \) is selected from current population, second chromosome \(\tilde{x}_{r2,G} \) is selected from historical data archive A which consists of recently explored inferior solutions to provide the information of progress direction and for improving the diversity of the population, and the third chromosome \(x_{r1,G} \) is selected randomly from the current population (Zhang and Sanderson 2009). The mutation vector strategy with external archive is achieved in the following manner:
where the indices \(r_1 \hbox { and }r_2 \) are two randomly selected distinct integers and \(F_i \) is the scaling factor.
Now, let us consider a chromosome \(x_{32,300} \), which represents the 32nd chromosome of 300th generation in target population of size, say, 100. For performing difference vector mutation, corresponding to chromosome \(x_{32,300} \), three other chromosomes were selected, viz. best chromosome \(x_\mathrm{best,300}^p\) from current population; second chromosome \(\tilde{x}_{116,300} \) which denotes 116th chromosome of 300th generation, randomly selected from the union of current population and historical data archive \((P\cup A)\); and third chromosome \(x_{66,300} \) which is 66th chromosome of 300th generation, randomly selected from current population. These three chromosomes appear as shown below:
Now difference vector mutation with a scaling factor F of 0.604909433 is performed according to (12). The resultant 32nd chromosome \(v_{32,300} \) of donor population for 300th generation is as given below.
1.2 Binomial crossover
The binomial crossover operation is performed on donor population to enhance the exploration capability of the present technique. This is accomplished as follows. For each decision variable in the population, a random number is generated between 0 and 1. If the generated random number is less than or equal to crossover rate (CR) then the decision variable is inherited from the donor vector \(v_{32,300}\), else decision variable is inherited from the target vector \(x_{32,300} \) and is defined as (Das and Suganthan 2011; Liu and Wang 2009; Brest et al. 2006; Wang et al. 2011b):
where \(i=1,\ldots ,N_{p}, j=1,\ldots ,D\), and q is a randomly chosen index \(\in \left\{ {1,\ldots ,N_p } \right\} \). The obtained trial chromosome \(u_{32,300} \) is given as follows:
Now if both target chromosome \(x_{32,300} \) and the trial chromosome \(u_{32,300} \) obtained after binary crossover are compared with the true values then the error between \(x_{32,300} \) and true values comes out to be 0.5402. This is high compared to error obtained using trial chromosome \(u_{32,300} \) which is 0.5319.
1.3 Selection
Finally, after performing the binomial crossover operation the selection process is performed between the target vector and the trial vector chromosomes depending upon the fitness of the each individual in order to select better individual to the next generation. Selection operation is described as (Das and Suganthan 2011; Liu and Wang 2009; Brest et al. 2006; Wang et al. 2011b):
By using selection operation, trial chromosome \(u_{32,300} \) has been found to be fittest chromosome and is then selected for the next generation.
At each generation, the control parameters CR and F for each target vector are updated by utilizing th e normal and the Cauchy distributions, respectively (Zhang and Sanderson 2009). The crossover rate \(\text {CR}_i \) of each chromosome \(x_i \) is achieved according to (15), and next, the value of \(\text {CR}_i \) is trimmed to \(\left[ {0,1} \right] \).
The mean \(\mu _\mathrm{CR} \) of normal distribution is initially initialized to 0.5, and this value is then updated for each iteration according to (16).
where \(S_\mathrm{CR} \) represents the set of all effective crossover rates \(\text {CR}_i\)’s at each generation G, c indicates a positive constant value between 0 and 1, and \(\hbox {mean}_A (\cdot )\) shows the arithmetic mean (Zhang and Sanderson 2009).
Similarly, the scaling factor \(F_i \) of each chromosome \(x_i \) is generated independently using (17) at every generation G and next the value \(F_i \) is trimmed to 1 if \(F_i \ge 1\) or restored if \(F_i \le 0\).
The location parameter \(\mu _F \) is initially set to 0.5, and at the end of each generation the parameter is updated as:
where \(\hbox {mean}_L (\cdot )\) is the Lehmer mean and \(S_F \) denotes the set of all successful scaling factors in generation G (Zhang and Sanderson 2009).
Rights and permissions
About this article
Cite this article
Basetti, V., Chandel, A.K. & Subramanyam, K.B.V.S.R. Power system static state estimation using JADE-adaptive differential evolution technique. Soft Comput 22, 7157–7176 (2018). https://doi.org/10.1007/s00500-017-2715-3
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00500-017-2715-3