A Multi-Rate Parallel Simulation Method for Active Distribution Network Based on Characteristic Matrix Partition Equivalence and Adaptive Interpolation Switching
Lou Guannan1, Jiang Xiaoyu1, Yang Zhichun2, Gu Wei1, Yang Fan2
1. School of Electrical Engineering Southeast University Nanjing 211096 China;
2. Electric Power Research Institute State Grid Hubei Electric Power Co. Wuhan 430074 China
With the rapid development of power electronics technology and increasing penetration of distributed generations, modern active distribution networks exhibit complex dynamic characteristics due to the multiple time scales of system operation, which poses challenges to the simultaneous demand of computational efficiency and accuracy in the conventional simulation method with a single simulation-rate. This paper proposes a multi-rate parallel simulation method based on the characteristic matrix partition equivalence and adaptive interpolation switching. Firstly, a state space matrix is constructed for the whole active distributed system, and the system is divided into multiple subsystems with different dynamic characteristics based on the magnitude of the eigenvalues of the state matrix and their corresponding natural frequencies. The iterative computation of the high-speed subsystem is equivalently calculated within a certain period of time through binomial expansion and multidimensional matrix summation formulas, and then the block equivalence of the system's state matrix is achieved. This method transforms the original full-region state quantity iterative computation into the independent simulation solutions for subregions at small step sizes and equivalent compensation at large step sizes, simplifying the original high-order matrix iterative computation to the low-order matrix iterative computation, thereby significantly reducing the simulation complexity while ensuring simulation accuracy. For the problem of asynchronous data interaction between regions, the error generation principle and main influencing factors of the traditional linear interpolation method are analyzed through the magnification step length vector method. Based on the second-order Taylor expansion, the interpolation algorithm is improved by introducing undetermined coefficients. A multi-step adaptive interpolation simulation algorithm is proposed, which can adjust relevant interpolation parameters according to the system frequency and simulation step length requirements, significantly reducing interpolation errors. Finally, using the improved IEEE123 node distribution network example, the feasibility and effectiveness of the proposed multi-rate parallel simulation method are verified through the comparison of partition equivalence/interpolation calculations and simulation results. The specific innovations of this article are as follows: Based on the discrete state space matrix, the system is considered as a multiple-rate counterpart due to different dynamic characteristics, and then the influence of coupling between fast and slow system regions on the accuracy of multi-rate simulation is evaluated through diagonalizing the matrix. The equivalent decoupling of the partitioned subsystem models is implemented to achieve independent resolution of subregions in small steps and equivalent compensation in large steps, significantly reducing the computational load while ensuring simulation accuracy. For the problem of asynchronous data interaction in multi-rate regions, the error caused by the conventional linear interpolation methods are analyzed by using the magnification step length vector method, and the relationships between the error and the system frequency and step length are also derived. Based on the second-order Taylor expansion, a multi-rate multi-step adaptive interpolation algorithm is proposed, which adjusts the interpolation algorithm parameters according to the simulation requirements for the first time to meet the minimum error condition, and performs optimal switching among multiple interpolation algorithms to improve the inter-region interaction coordination and simulation accuracy. Through simulation comparison, the following conclusions can be drawn:
(1) This method evaluates the impact of fast and slow partition coupling on the accuracy of multi rate simulation numerical values, and decouples the subsystem model equivalently to achieve independent solution of the subregion in small steps and equivalent compensation in large steps, reducing the order of the state iteration matrix during simulation and significantly improving efficiency while ensuring computational accuracy.
(2) On the basis of analyzing the principle and influencing factors of linear interpolation errors, an adaptive interpolation switching simulation algorithm is proposed. Relevant interpolation parameters are adjusted according to the system frequency and simulation step size requirements, and interpolation algorithm optimization is carried out to improve the accuracy and stability of parallel simulation.
楼冠男, 蒋啸宇, 杨志淳, 顾伟, 杨帆. 基于特征矩阵分区等值和自适应插值切换的有源配电网多速率并行仿真方法[J]. 电工技术学报, 0, (): 9032-32.
Lou Guannan, Jiang Xiaoyu, Yang Zhichun, Gu Wei, Yang Fan. A Multi-Rate Parallel Simulation Method for Active Distribution Network Based on Characteristic Matrix Partition Equivalence and Adaptive Interpolation Switching. Transactions of China Electrotechnical Society, 0, (): 9032-32.
[1] 蔡瑶, 卢志刚, 潘尧, 等. 计及多重差异的交直流混合多能微网多时间尺度优化调度[J/OL]. 电工技术学报: 1-19[2023-11-11]. https://doi.org/10.19595/j.cnki.1000-6753.tces.230645.
Cai Yao, Lu Zhigang, Pan Yao, et al. Multi-time-scale Optimal Scheduling of AC-DC Hybrid Multi-energy Microgrid Considering Multiple Differences[J]. Transactions of China Electrotechnical Society: 1-19[2023-11-11].https://doi.org/10.19595/j.cnki.1000-6753.tces.230645.
[2] 马钊, 安婷, 尚宇炜. 国内外配电前沿技术动态及发展[J]. 中国电机工程学报, 2016, 36(6): 1552-1567.
Ma Zhao, An Ting, Shang Yuwei.State of the art and development trends of power distribution technologies[J]. Proceedings of the CSEE, 2016, 36(6): 1552-1567.
[3] 张玉莹, 曾博, 周吟雨, 等. 碳减排驱动下的数据中心与配电网交互式集成规划研究[J/OL]. 电工技术学报: 1-17[2023-11-11]. https://doi.org/10.19595/j.cnki.1000-6753.tces.221646.
Zhang Yuying, Zeng Bo, Zhou Yinyu, et al. Research on Interactive Integration Planning of Data Centers and Distribution Network Driven by Carbon Emission Reduction[J/OL]. Transactions of China Electrotechnical Society: 1-17[2023-11-11]. https://doi.org/10.19595/j.cnki.1000-6753.tces.221646.
[4] Hadjsa N.有源智能配电网[M]. 陶顺,肖湘宁,彭骋, 译. 北京: 中国电力出版社, 2012.
[5] 王成山, 王丹, 周越. 智能配电系统架构分析及技术挑战[J]. 电力系统自动化, 2015, 39(9): 2-9.
Wang Chengshan, Wang Dan, Zhou Yue.Framework analysis and technical challenges to smart distribution system[J]. Automation of Electric Power Systems, 2015, 39(9): 2-9.
[6] 江岳文, 罗泽宇, 程诺. 基于线性化方法的交直流混合配电系统网架规划[J/OL]. 电工技术学报: 1-16[2023-11-11]. https://doi.org/10.19595/j.cnki.1000-6753.tces.222363.
Jiang Yuewen, Luo Zeyu, Cheng Nuo. Network Planning of AC/DC Hybrid Distribution System Based on Linearization Method [J]. Transactions of China Electrotechnical Society: 1-16. [2023-11-11].https://doi.org/10.19595/j.cnki.1000-6753.tces.222363.
[7] 何绍民, 张喆, 卢倚平, 等. 基于计算前沿面的实时仿真数值积分并行构造及其数值模型解耦加速方法[J]. 电工技术学报, 2023, 38(16): 4246-4262.
He Shaomin, Zhang Zhe, Lu Yiping, et al.Numerical model decoupling acceleration method with numerical integration parallelism construction based on computation front in real-time simulation[J]. Transactions of China Electrotechnical Society, 2023, 38(16): 4246-4262.
[8] Hariri A, Faruque M O.A hybrid simulation tool for the study of PV integration impacts on distribution networks[J]. IEEE Transactions on Sustainable Energy, 2017, 8(2): 648-657.
[9] Armstrong M, Marti J R, Linares L R, et al.Multilevel MATE for efficient simultaneous solution of control systems and nonlinearities in the OVNI simulator[J]. IEEE Transactions on Power Systems, 2006, 21(3): 1250-1259.
[10] Tomim M A, Martí J R, De Rybel T, et al.MATE network tearing techniques for multiprocessor solution of large power system networks[C]//IEEE PES General Meeting, Minneapolis, MN, USA, 2010: 1-6.
[11] 邓潘, 盛万兴, 刘科研, 等. 有源配电网典型暂态模型及区域分割并行仿真方法研究[J]. 电网技术, 2020, 44(4): 1211-1219.
Deng Pan, Sheng Wanxing, Liu Keyan, et al.Methods of transient model and parallel simulation of region segmentation in active distribution network[J]. Power System Technology, 2020, 44(4): 1211-1219.
[12] 王照琪, 唐巍, 张博, 等. 基于优化分网策略的有源配电网多速率并行暂态仿真分析[J]. 电网技术, 2020, 44(2): 673-682.
Wang Zhaoqi, Tang Wei, Zhang Bo, et al.Multi-rate parallel transient simulation of active distribution network based on optimization decomposition strategy[J]. Power System Technology, 2020, 44(2): 673-682.
[13] 刘科研, 叶华, 裴玮, 等. 基于特征值分析的配电网多速率暂态仿真性能评估方法研究[J]. 电网技术, 2020, 44(11): 4088-4096.
Liu Keyan, Ye Hua, Pei Wei, et al.Research on computational performance of multi-rate transient simulation in distribution network based on eigen-value analysis[J]. Power System Technology, 2020, 44(11): 4088-4096.
[14] 魏松韬. 复杂有源配电网电磁暂态多速率并行仿真理论与方法[D]. 南京: 东南大学, 2022.
[15] Heffernan M D, Turner K S, Arrillaga J, et al. Computation of AC-DC system disturbances - part I. interactive coordination of generator and convertor transient models[J]. IEEE Power Engineering Review, 1981, PER-1(11): 15-16.
[16] 王耀坤, 赵煜, 梁英, 等. 基于信息物理系统的配电网故障处理仿真研究[J]. 智慧电力, 2019, 47(12): 103-109.
Wang Yaokun, Zhao Yu, Liang Ying, et al.Simulation on distribution network fault handling based on CyberPhysical system[J]. Smart Power, 2019, 47(12): 103-109.
[17] 李鹏, 曾凡鹏, 王智颖, 等. 基于FPGA的有源配电网实时仿真器I/O接口设计[J]. 中国电机工程学报, 2017, 37(12): 3409-3417, 3668.
Li Peng, Zeng Fanpeng, Wang Zhiying, et al.I/O interface design for FPGA based real-time simulator of active distribution networks[J]. Proceedings of the CSEE, 2017, 37(12): 3409-3417, 3668.
[18] 陈颖, 宋炎侃, 黄少伟, 等. 基于GPU的大规模配电网电磁暂态并行仿真技术[J]. 电力系统自动化, 2017, 41(19): 82-88.
Chen Ying, Song Yankan, Huang Shaowei, et al.GPU-based techniques of parallel electromagnetic transient simulation for large-scale distribution network[J]. Automation of Electric Power Systems, 2017, 41(19): 82-88.
[19] 周斌, 汪光森, 李卫超, 等. 基于FPGA的电力电子系统电磁暂态实时仿真通用解算器[J]. 电工技术学报, 2023, 38(14): 3862-3874.
Zhou Bin, Wang Guangsen, Li Weichao, et al.An FPGA-based general solver for electromagnetic transient real-time simulation of power electronic systems[J]. Transactions of China Electrotechnical Society, 2023, 38(14): 3862-3874.
[20] 占鹏. 基于多核CPU的电力系统多速率电磁暂态仿真[D]. 北京: 华北电力大学, 2018.
[21] Benigni A, Monti A, Dougal R A.Latency-based approach to the simulation of large power electronics systems[J]. IEEE Transactions on Power Electronics, 2014, 29(6): 3201-3213.
[22] Li Yupeng, Shu Dewu, Hu Jingwei, et al.A multi-area Thevenin equivalent based multi-rate co-simulation for control design of practical LCC HVDC system[J]. International Journal of Electrical Power & Energy Systems, 2020, 115: 105479.
[23] 韩佶, 董毅峰, 苗世洪, 等. 基于MATE的电力系统分网多速率电磁暂态并行仿真方法[J]. 高电压技术, 2019, 45(6): 1857-1865.
Han Ji, Dong Yifeng, Miao Shihong, et al.Multi-rate electromagnetic transient parallel simulation of power system based on MATE[J]. High Voltage Engineering, 2019, 45(6): 1857-1865.
[24] Bartel A, Günther M.A multirate W-method for electrical networks in state-space formulation[J]. Journal of Computational and Applied Mathematics, 2002, 147(2): 411-425.
[25] 李珂, 顾伟, 柳伟, 等. 基于FPGA的变流器并行多速率电磁暂态实时仿真方法[J]. 电力系统自动化, 2022, 46(13): 151-158.
Li Ke, Gu Wei, Liu Wei, et al.Real-time parallel multi-rate electromagnetic transient simulation method for converters based on field programmable gate array[J]. Automation of Electric Power Systems, 2022, 46(13): 151-158.
[26] do Couto Boaventura W, Semlyen A, Iravani M R, et al. Robust sparse network equivalent for large systems: part I-methodology[J]. IEEE Transactions on Power Systems, 2004, 19(1): 157-163.
[27] 郑雷, 王斌, 喻敏. 基于分形插值的电力系统故障录波数据压缩与还原算法研究[J]. 计算机测量与控制, 2014, 22(4): 1185-1188.
Zheng Lei, Wang Bin, Yu Min.Research about compression and reproduction algorithm of digital recorded data from a faulted power system based on fractal interpolation[J]. Computer Measurement & Control, 2014, 22(4): 1185-1188.
[28] Xu Junjun, Wu Zaijun, Yu Xinghuo, et al.An interval arithmetic-based state estimation framework for power distribution networks[J]. IEEE Transactions on Industrial Electronics, 2019, 66(11): 8509-8520.