Research Article  Open Access
Improvement of Power Flow Calculation with Optimization Factor Based on Current Injection Method
Abstract
This paper presents an improvement in power flow calculation based on current injection method by introducing optimization factor. In the method proposed by this paper, the PQ buses are represented by current mismatches while the PV buses are represented by power mismatches. It is different from the representations in conventional current injection power flow equations. By using the combined power and current injection mismatches method, the number of the equations required can be decreased to only one for each PV bus. The optimization factor is used to improve the iteration process and to ensure the effectiveness of the improved method proposed when the system is illconditioned. To verify the effectiveness of the method, the IEEE test systems are tested by conventional current injection method and the improved method proposed separately. Then the results are compared. The comparisons show that the optimization factor improves the convergence character effectively, especially that when the system is at high loading level and R/X ratio, the iteration number is one or two times less than the conventional current injection method. When the overloading condition of the system is serious, the iteration number in this paper appears 4 times less than the conventional current injection method.
1. Introduction
Power flow studies are necessary for planning, operating, economic scheduling, and other analysis, such as transient stability, voltage stability, and contingency studies. The tasks of power flow calculation are to solve the steadystate operating conditions of power systems based on the operation modes and the wiring of the systems. Different power flow methods have emerged and continue to be developed with the changing needs.
The conventional power flow solution comprises power equations expressed in terms of rectangular or polar coordinates. Many important contributions have been reported in this area [1–3]. To improve the convergence characteristic of GaussSeidel method, the NewtonRaphson method was introduced. The NR method was once considered the state of the art powerflow technique and widely accepted in industry applications. However, the main disadvantage of the NewtonRaphson method is the necessity for factorizing and updating the Jacobian matrix during the iterative solution process [4]. To solve this problem, Fast Decoupled power flow method was proposed to speed up the iteration process of the NR method and decrease the required minimum memory storage. Nevertheless, the convergence rate of the decoupled NR is influenced by the illconditioned case, such as when a system has high ratio of line R/X [5–8]. Some other methods have been presented in other forms, such as the use of sequence component frame [9–11] and the method based on the loop frame of reference [12].
In [13], a new power flow method based on current injection was presented, in which the current injection equations are written in rectangular coordinates. The Jacobian matrix is composed of block matrices and has the same structure as the nodal admittance matrix. By using this method, the Jacobian matrix can be updated faster than using conventional NR power flow method in the case of PQ buses. However, in the case of the existence of PV buses, the current injection mismatch power flow method was presented to improve the convergence character of current injection method [14, 15]. The improvement in [16] increased the number of required equations to three for each PV bus. The paper [17] presented a new representation of PV buses in revised current injection power flow method to reduce the number of required equations and to improve the convergence character of well and illconditioned systems, while the elements related to PV buses have to be changed during the iteration process. In [18] a development of power flow calculation for illconditioned system was presented. The fact that “the Taylor series expansion of the load flow equations is expressed up to the third term completely and the final term has the same form but different variable as the first term” was used in [18]. Following the main ideas presented in [18], which is called Iwamoto’s method in [15], a new second order power flow method was proposed. This method is useful for heavily loaded and overloaded systems as well as illconditioned distribution systems.
This paper improved the current injection method presented in [17] following the main idea in [18]. The optimization presented in [18] is briefly described here firstly to make the paper clear. Then the improved method, results, and conclusions are presented.
2. Current Injection Power Flow
The basic current injection method proposed in [17] is used in this paper. The combined power and current injection mismatches power flow formulation can be calculated from the following:
where : real parts of current mismatch at bus , : imaginary parts of current mismatch at bus , : active power mismatch at bus , : real voltage component correction at bus , : imaginary voltage component correction at bus , and : voltage angle correction at bus .
2.1. Equations for PQ Buses
The current mismatch for a given bus is where : complex conjugated voltage phasor at bus , : specified active and reactive net powers at bus , : bus admittance matrix element, and : voltage phasor at bus . Consider where , : active and reactive powers of generators for bus and , : active and reactive powers of loads for bus .
Equation (2) can be expanded into its real and imaginary components: where : real voltage component at bus and : imaginary voltage component at bus .
Equations (4) and (5) are written in compact forms: where , : specified real and imaginary parts of current at bus , , : calculated real and imaginary parts of current at bus .
Using the NewtonRaphson solution algorithm, the elements in (1) for all buses as being of the PQ type are given by The diagonal elements are given by The elements , , , and are presented in Appendix of [13].
The following equation is expression of current mismatches in (1):
2.2. Representation of PV Buses
The calculated active power at assumed bus can be calculated by where , : voltage phase angle at bus .
In the Jacobian matrix of (1), the elements related to buses of PV type are given as follows.
Diagonal elements
Off diagonal elements
3. Description of Iwamoto’s Method
The method presented in [18] is described briefly here to make the improved method proposed in next part clearly. The conventional power flow equations in the rectangular coordinates are The Taylor series expansion of (13) turns out to be where : estimate of , : Jacobian matrix, and : error (correction vector).
Moving all the righthand side of (14) to the lefthand side
To obtain the solution of this equation, it denotes by , the step size optimization factor. Applying this factor to this equation, then it follows that Or It can be rewritten as where The following cost function is considered to determine the value of the in a least squared sense: The optimization solution is given by Namely, where The state vector for the next iteration is given by
4. Improvement in Current Injection Method
4.1. Parameters for PQ Buses
The improved method is based on the current injection method by defining a step size optimization factor. The equation (5) is the expression of imaginary current injection related to bus . Neglecting all terms of order higher than two, we then get Hence, the following equations can be got based on the theory in part 3: In a similar way, the following equations can be obtained Thus, the parameters can be obtained as follows
4.2. Parameters for PV Buses
From (1) for a PV bus one has Using the same method as parameters for PQ buses, the following equations can be obtained
5. Results
This paper used MATLAB to implement the program of power flow method proposed. IEEE 30bus and 118bus test system are tested here. Firstly the validity of the improved method proposed is proven. Here the systems are tested as symmetrical systems. Tables 1 and 2 show a part of the results. Then the results obtained by conventional current injection and the improved method proposed, respectively, are compared. Figures 1 and 2 show the iteration characteristics of IEEE 30bus test system using the two different methods and Figures 3 and 4 show the results of IEEE 118bus system in the same way.
 
0—slack bus, 1—PQ bus, 2—PV bus. 
 
0—slack bus, 1—PQ bus, 2—PV bus. 
Tables 3 and 4 compare the conventional current injection method and the improved method proposed here when the system has high R/X ratios and compare the two methods when the system is operating in overloaded condition.


(a) Validity of the Method. Tables 1 and 2 show the results of the iteration and present the validity of the improved method. The voltage amplitudes and angles are reasonable and acceptable. At the same time the feasibility of the method is proved.
(b) Convergence Characteristic. The iteration convergence criterion for the program of the improved method proposed is that the maximum active power mismatch and the maximum reactive power mismatch are both less than the tolerance. The results presented in Figures 1–4 show that the convergence characteristic of the improved method proposed is better than the conventional current injection method.
(c) High R/X Ratio. Many power flow methods always have bad convergence characteristic if the system has high R/X ratio. In Table 3 the R/X ratios are divided into five levels. When the system has high X/R ratio, both the two methods have good convergence characteristic, but with the increase of R/X ratio the advantage of the improved method proposed becomes obvious; the iteration number is one or two times less than the conventional current injection method.
(d) Overloading. The comparison of the iteration number for different loading cases is shown in Table 4. In this case the system is in overloading condition. In Table 4, five loading cases are presented. There has been little difference between the two methods when the overloading of the system is not serious. However, when the condition becomes more serious the iteration number of the improved method proposed appears 4 times less than the conventional current injection method. Even more, when the convergence characteristic of the improved method proposed is still acceptable when overloading is serious, the conventional current injection method shows no convergence.
6. Conclusions
This paper has presented improvement in power flow calculation based on current injection method by introducing an optimization factor. Unlike other current injection methods, the PV buses are represented by power mismatches here. The results have demonstrated the good performance of the improved method proposed in this paper. By using optimization factor the iteration character got improved. Even when the system has high R/X ratio or operates in overloaded conditions, good results could be got by using this method.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was supported in part by the National Natural Science Foundation of China (no. 51205046) and the Fundamental Research Funds for the Central Universities (nos. CDJZR170008, 106112014CDJZR175501). The constructive comments provided by the anonymous reviewers and the editors are also greatly appreciated.
References
 H. YingYi and W. FuMing, “Development of threephase Newton optimal power flow for studying imbalance/security in transmission systems,” Electric Power Systems Research, vol. 55, no. 1, pp. 39–48, 2000. View at: Publisher Site  Google Scholar
 B. Stott, “Review of loadflow calculation methods,” Proceedings of the IEEE, vol. 62, no. 7, pp. 916–929, 1974. View at: Publisher Site  Google Scholar
 A. G. Expósito and E. R. Ramos, “Augmented rectangular load flow model,” IEEE Transactions on Power Systems, vol. 17, no. 2, pp. 271–276, 2002. View at: Publisher Site  Google Scholar
 C. Xing, SteadyState Analysis of Power Systems, vol. 262, Chia Electric Power Press, 3rd edition, 2007.
 W.M. Lin and J.H. Teng, “Threephase distribution network fastdecoupled power flow solutions,” International Journal of Electrical Power and Energy System, vol. 22, no. 5, pp. 375–380, 2000. View at: Publisher Site  Google Scholar
 V. M. da Costa, M. L. de Oliveira, and M. R. Guedes, “Developments in the analysis of unbalanced threephase power flow solutions,” International Journal of Electrical Power and Energy Systems, vol. 29, no. 2, pp. 175–182, 2007. View at: Publisher Site  Google Scholar
 Y. Song, X. Li, and W. Cai, “Adaptive and faulttolerant reactive power compensation in power systems via multilevel STATCOMs,” International Journal of Innovative Computing, Information and Control, vol. 9, no. 8, pp. 3403–3413, 2013. View at: Google Scholar
 D. Nlu and Y. Wei, “A novel socialenvironmentaleconomic dispatch model for thermal/wind power generation and application,” International Journal of Innovative Computing, Information and Control, vol. 9, no. 7, pp. 3005–3014, 2013. View at: Google Scholar
 K. L. Lo and C. Zhang, “Decomposed threephase power flow solution using the sequence component frame,” IEE Proceedings C: Generation Transmission and Distribution, vol. 140, no. 3, pp. 181–188, 1993. View at: Publisher Site  Google Scholar
 M. AbdelAkher, K. M. Nor, and A. H. A. Rashid, “Improved threephase powerflow methods using sequence components,” IEEE Transactions on Power Systems, vol. 20, no. 3, pp. 1389–1397, 2005. View at: Publisher Site  Google Scholar
 I. Ngamroo, “Dynamic events analysis of Thailand and Malaysia power systems by discrete wavelet decomposition and short term fourier transform based on GPS synchronized phasor data,” International Journal of Innovative Computing, Information and Control, vol. 9, no. 5, pp. 2203–2228, 2013. View at: Google Scholar
 T. Chen and N. Yang, “Loop frame of reference based threephase power flow for unbalanced radial distribution systems,” Electric Power Systems Research, vol. 80, no. 7, pp. 799–806, 2010. View at: Publisher Site  Google Scholar
 P. A. N. Garcia, J. L. R. Pereira, S. Carneiro, and V. M. Da Costa, “Threephase power flow calculations using the current injection method,” IEEE Transactions on Power Systems, vol. 15, no. 2, pp. 508–514, 2000. View at: Publisher Site  Google Scholar
 P. A. N. Garcia, J. L. R. Pereira, S. Carneiro Jr., M. P. Vinagre, and F. V. Gomes, “Improvements in the representation of PV buses on threephase distribution power flow,” IEEE Transactions on Power Delivery, vol. 19, no. 2, pp. 894–896, 2004. View at: Publisher Site  Google Scholar
 C. A. Ferreira and V. M. Da Costa, “A second order power flow based on current injection equations,” International Journal of Electrical Power and Energy Systems, vol. 27, no. 4, pp. 254–263, 2005. View at: Publisher Site  Google Scholar
 V. M. da Costa, N. Martins, and J. L. R. Pereira, “Developments in the newton raphson power flow formulation based on current injections,” IEEE Transactions on Power Systems, vol. 14, no. 4, pp. 1320–1326, 1999. View at: Publisher Site  Google Scholar
 S. Kamel, M. AbdelAkher, and F. Jurado, “Improved NR current injection load flow using power mismatch representation of PV bus,” International Journal of Electrical Power and Energy Systems, vol. 53, no. 1, pp. 64–68, 2013. View at: Publisher Site  Google Scholar
 S. Iwamoto and Y. Tamura, “A load flow calculation method for illconditioned power systems,” IEEE Transactions on Power Apparatus and Systems, vol. 100, no. 4, pp. 1736–1743, 1979. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Lei Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.