Abstract It is very important to obtain the active power output of each generator by solving the optimal power flow (OPF) for the economic operation of power system. The traditional methods mainly include simplified model, convex relaxation and intelligent algorithm, but these methods also have some shortcomings such as imprecision and dependence on the initial value. Holomorphic embedding method (HEM) is a new method used in power flow calculation and it ensures convergence through proper holomorphic embedding forms. HEM has the advantages of non-iteration, recursion, analyzability, no dependence on the initial value, and the obtained solution is a runnable point. Based on these above advantages, HEM is used to solve the OPF flow, and a non-iterative OPF calculation method based on HEM is proposed. Firstly, taking a 3-bus power system as an example, the objective function is the lowest generating cost, and the active power output of the generator is optimized. Secondly, constraints and KKT conditions of the optimization model are listed. Based on the basic construction rules of HEM, the holomorphic embedding forms of constraints and KKT conditions are constructed. Thirdly, according to the basic properties of holomorphic functions, the recursive equations of power series coefficients of analytic expressions of bus voltage and generator active power output are derived. Finally, coefficients can be obtained by the cross recursion of the above two sets of equations, and the generating cost is obtained. Simulation results on IEEE 3-bus, 4-bus, 30-bus and 39-bus New England power systems show that the HEM is convergent when solving the OPF. The calculation result of the active power output of generators and the bus voltage are close to the reference value, and the error is less than 1% in the first two power systems. In the four systems, the convergence accuracy is achieved at the 6th, 8th, 10th and 12th order respectively. For the result of the objective function, the error of the results calculated by HEM reaches to 10-6 in the first two systems. The source of error is the order of recursion, and the process of each recursion is to make the result gradually approximate to the exact solution. However, the error increases with the increase of the scale of the system. In New England 39-bus system, the error is larger because inequality constraints are not considered in the system, which makes the active power output more ideal, so that the generating cost and active power loss are smaller. In the case of a small system scale, the holomorphic embedding method proposed in this paper has a shorter computing time than the interior point method (IPM). However, the computing time in IEEE 30-bus and New England 39-bus systems is longer than the IPM. Finally, the feasibility of the HEM in IEEE 118 bus power system is verified. The following conclusions can be drawn through the simulation analysis: (1) HEM can be used to obtain the analytical expressions of the bus voltage and the active power output of the generator. The coefficients of the power series are obtained through the cross recursion of the KKT condition and constraints in holomorphic embedding form. (2) Due to the constraint of power flow included in the model, the optimal result obtained is operable in the real system, which has the advantage of being independent of the initial value. (3) The calculation results of HEM meet the accuracy requirements, and it is proved that HEM is feasible in the calculation of OPF.
Liu Chengxi,Xu Shenkai,Lai Qiupin. Non-Iterative Optimal Power Flow Calculation Based on Holomorphic Embedding Method[J]. Transactions of China Electrotechnical Society, 2023, 38(11): 2870-2882.
[1] 国家发展改革委, 国家能源局. 电力发展“十三五”规划[R]. 2016. [2] 林哲, 胡泽春, 宋永华. 最优潮流问题的凸松弛技术综述[J]. 中国电机工程学报, 2019, 39(13): 3717-3728. Lin Zhe, Hu Zechun, Song Yonghua.Convex relaxation for optimal power flow problem: a recent review[J]. Proceedings of the CSEE, 2019, 39(13): 3717-3728. [3] Carpentier J.Contribution to the economic dispatch problem[J]. Bulletin de la Societe Francoise des Electriciens, 1962, 3(8): 431-447. [4] 刘科研, 盛万兴, 李运华. 互联电网的直流最优潮流分解算法研究[J]. 中国电机工程学报, 2006, 26(12): 21-25. Liu Keyan, Sheng Wanxing, Li Yunhua.Research on decomposition algorithm of DC optimal power flow in large scale interconnection power grids[J]. Proceedings of the CSEE, 2006, 26(12): 21-25. [5] 王峥, 丛培杰. 基于改进前推回代法的辐射状配电网潮流计算[J]. 东北电力技术, 2008, 29(2): 7-10. Wang Zheng, Cong Peijie.Power flow calculation for radial distribution network based on forward and backward substitution method[J]. Northeast Electric Power Technology, 2008, 29(2): 7-10. [6] 麻秀范. 含分布式电源的配电网规划与优化运行研究[D]. 北京: 华北电力大学, 2013. [7] Zhang Ning, Cheng Jiangnan, Wang Yi.Probabilistic optimal energy flow of district multienergy systems: an MPLP-based online dictionary-learning approach[J]. IEEE Transactions on Industrial Informatics, 2019, 15(9): 4867-4877. [8] Pan Shixian, Li Zhigang, Zheng J H, et al.A non-iterative solution method for DC optimal power flow based on holomorphic embedding[C]//2021 IEEE 5th Conference on Energy Internet and Energy System Integration (EI2), Taiyuan, China, 2022: 1865-1870. [9] 朱泽安, 王承民, 马明, 等. 基于拉格朗日函数鞍距分配的广义内点法[J]. 电力系统自动化, 2018, 42(5): 18-24. Zhu Zean, Wang Chengmin, Ma Ming, et al.Generalized interior-point method based on Lagrange function saddle space distribution[J]. Automation of Electric Power Systems, 2018, 42(5): 18-24. [10] Gan Lingwen, Low S H.An online gradient algorithm for optimal power flow on radial networks[J]. IEEE Journal on Selected Areas in Communications, 2016, 34(3): 625-638. [11] 唐坤杰, 董树锋, 朱炳铨, 等. 大规模输配一体化系统牛顿法潮流计算性能分析及改进方法[J]. 电力系统自动化, 2019, 43(6): 92-99. Tang Kunjie, Dong Shufeng, Zhu Bingquan, et al.Performance analysis and improvement of Newton method for power flow calculation of large-scale integrated transmission and distribution network[J]. Automation of Electric Power Systems, 2019, 43(6): 92-99. [12] Luo Zhiquan, Ma W K, So A M C, et al. Semidefinite relaxation of quadratic optimization problems[J]. IEEE Signal Processing Magazine, 2010, 27(3): 20-34. [13] 巨云涛, 黄炎, 张若思. 基于二阶锥规划凸松弛的三相交直流混合主动配电网最优潮流[J]. 电工技术学报, 2021, 36(9): 1866-1875. Ju Yuntao, Huang Yan, Zhang Ruosi.Optimal power flow of three-phase hybrid AC-DC in active distribution network based on second order cone programming[J]. Transactions of China Electrotechnical Society, 2021, 36(9): 1866-1875. [14] 李洪美, 崔翰韬, 万秋兰. 考虑电动汽车充电策略的配网重构二阶锥规划模型[J]. 中国电机工程学报, 2015, 35(18): 4674-4681. Li Hongmei, Cui Hantao, Wan Qiulan.Distribution network reconfiguration based on second-order conic programming considering EV charging strategy[J]. Proceedings of the CSEE, 2015, 35(18): 4674-4681. [15] 郭清元, 吴杰康, 莫超, 等. 基于混合整数二阶锥规划的新能源配电网电压无功协同优化模型[J]. 中国电机工程学报, 2018, 38(5): 1385-1396. Guo Qingyuan, Wu Jiekang, Mo Chao, et al.A model for multi-objective coordination optimization of voltage and reactive power in distribution networks based on mixed integer second-order cone programming[J]. Proceedings of the CSEE, 2018, 38(5): 1385-1396. [16] Sennewald T, Linke F, Reck J, et al.Parameter optimization of differential evolution and particle swarm optimization in the context of optimal power flow[C]//2020 IEEE PES Innovative Smart Grid Technologies Europe (ISGT-Europe), Hague, Netherlands, 2020: 1045-1049. [17] 刘明凯, 王占山, 邢彦丽. 基于强化多目标差分进化算法的电-气互联系统最优潮流计算[J]. 电工技术学报, 2021, 36(11): 2220-2232. Liu Mingkai, Wang Zhanshan, Xing Yanli.Enhanced multi-objective differential evolutionary algorithm based optimal power flow calculation for integrated electricity and gas systems[J]. Transactions of China Electrotechnical Society, 2021, 36(11): 2220-2232. [18] 王浩翔, 赵冬梅, 陶然, 等. 基于分解的多目标进化算法的含MMC-HVDC交直流混合系统最优潮流研究[J]. 电工技术学报, 2020, 35(17): 3691-3702. Wang Haoxiang, Zhao Dongmei, Tao Ran, et al.Study on optimal power flow for AC/DC hybrid system incorporating MMC-HVDC based on MOEA/D[J]. Transactions of China Electrotechnical Society, 2020, 35(17): 3691-3702. [19] Moeini-Aghtaie M, Abbaspour A, Fotuhi-Firuzabad M, et al.A decomposed solution to multiple-energy carriers optimal power flow[J]. IEEE Transactions on Power Systems, 2014, 29(2): 707-716. [20] Koh C S, Ryu J S, Fujiwara K.Convergence acceleration of the Newton-Raphson method using successive quadratic function approximation of residual[J]. IEEE Transactions on Magnetics, 2006, 42(4): 611-614. [21] Li Mingyang, Zhao Qianchuan, Luh P B.Decoupled load flow and its feasibility in systems with dynamic topology[C]//2009 IEEE Power & Energy Society General Meeting, Calgary, AB, Canada, 2009: 1-8. [22] 姜涛, 张勇, 李雪, 等. 电力系统交直流潮流的全纯嵌入计算[J]. 电工技术学报, 2021, 36(21): 4429-4443, 4481. Jiang Tao, Zhang Yong, Li Xue, et al.A holomorphic embedded method for solving power flow in hybrid AC-DC power system[J]. Transactions of China Electrotechnical Society, 2021, 36(21): 4429-4443, 4481. [23] Liu Chengxi, Wang Bin, Hu Fengkai, et al.Online voltage stability assessment for load areas based on the holomorphic embedding method[C]//2018 IEEE Power & Energy Society General Meeting (PESGM), Portland, 2018: 1. [24] Liu Chengxi, Sun Kai, Wang Bin, et al.Probabilistic power flow analysis using multi-dimensional holomorphic embedding and generalized cumulants[C]// 2019 IEEE Power & Energy Society General Meeting (PESGM), Atlanta, GA, USA, 2019: 1. [25] 刘承锡, 赖秋频, 姚良忠, 等. 基于全纯嵌入法的电力系统解耦模型:薄弱节点辨识[J]. 中国电机工程学报, 2022, 42(5): 1736-1748. Liu Chengxi, Lai Qiupin, Yao Liangzhong, et al.Power system decoupling model based on holomorphic embedding method: weak bus identification[J]. Proceedings of the CSEE, 2022, 42(5): 1736-1748. [26] Yao Rui, Liu Yang, Sun Kai, et al.Efficient and robust dynamic simulation of power systems with holomorphic embedding[J]. IEEE Transactions on Power Systems, 2020, 35(2): 938-949. [27] Kundur P, Balu N J, Lauby M G.Power system stability and control[M]. New York: McGraw-Hill, 1994 [28] Trias A, Marín J L.A Padé-Weierstrass technique for the rigorous enforcement of control limits in power flow studies[J]. International Journal of Electrical Power & Energy Systems, 2018, 99: 404-418. [29] 陈珩. 电力系统稳态分析[M]. 3版. 北京: 中国电力出版社, 2007.
2023, Vol. 38 
