电力系统基础仿真算法对比分析研究

1、摘要微分代数方程的解算方法是电力系统仿真技术的基础，不同的解算方法得到仿真模型的解算结果之间存在差异，这对模型的计算分析及评估有重要影响。基于当前较为常用的商用软件平台中使用的三种分析方法一EMTP类软件的节点分析法、用于 Matlab/Simulink 的状态空间法以及用于 RT-LAB平台的状态空间节点法，结合各算法的特性，基于典型 RLC 电路模型进行仿真，通过仿真结果说明各算法之间的差异，并理论分析了结果差异产生的 原因。关键词：微分代数方程；隐式梯形法；状态空间法；状态空间节点法0引言当前电力电子技术的迅速发展及电网中电力电子器件的大量引入 使得电力系统中许多计算和控制问题日益复杂。

2、在一定程度上来说，电 力试验几乎是不可能直接进行的，因此适时地、合理地借助电力系统的 仿真技术是非常必要的。而电力系统仿真技术对计算机运行速率以及数 值计算方法具有较强的依赖性。随着计算机技术的不断发展，电力系统仿真技术成为电力系统规划、保护、调度甚至故障研究的重要工具，同 时仿真结果的准确性对指导电力系统安全、可靠、经济运行具有重要的 影响，而仿真结果的准确性主要依赖于基础算法的研究工作。目前电力系统仿真技术基础算法有大量的研究成果，当前较为成熟的电力系统仿真平台中使用的建模分析方法大体可分为：节点分析法(nodal analysis method )、状态空间法 (state-space

3、techniques)和状态空间节点法(state space nodal, SSN )。这三种分析方法使用 不同的数值计算方法对电路进行求解。本文就这三类方法针对具体的 RLC谐振电路仿真结果显现的差异性进行对比分析研究，展现各种算法之间的差异性。1三种方法基本原理简析数值积分方法是电力系统暂态稳定性分析计算的基本方法。数值积分法主要分为显式数值积分法和隐式数值积分法。当用于求解微分方程时，显式算法求解是对时间进行差分，不存在迭代和收敛的问题，因此 求解效率高，但计算精度不高，并需设定合适的时间步长求解以保证系 统的稳定，较大或较小的时间步长都会导致求解时间很长；隐式算法求解与时间无关，采用

4、的是牛顿迭代法，因此存在迭代收敛问题，这类算 法的计算精度较高。1.1 节点分析法利用节点分析法进行电路求解的仿真软件统称为EMTP (Electromagnetic Transient Program )类软件。代表性的有离线仿真软件PSCAD/EMTDC、EMTPRV ,实时数字仿真系统 RTDS等。EMTP 类软件的实质是利用隐式梯形法解算电路暂态参数的过程。 对隐式梯形 法的理解：对于一常微分方程，即，=/(工)(I)曲在t- At到t积分步长内的隐式梯形积分公式为x(7) = x(/-A/) + /x(f) + /x(f-Af)j(2)隐式梯形法的核心思想是把网络中的所有元件模型都用

5、不变的等值电阻和反映历史记录的等值电流源来代替。然后根据这个等值电路建立节点方程式，对其形成的导纳矩阵求逆，即可求出所需的解1。隐式梯形积分法是电力系统仿真分析中行之有效的求解算法，主要 由于它数值稳定性好，可以采用大步长，可通过联立求解同步发电机的差分方程和网络的代数方程来消除接口误差。但是，当电力网络中存在断路器操作时，采用这种方法将使非状态变量产生不正常的摆动，造成数值震荡。针对这种现象，研究人员将CDA ( critical dampingadjustment )引入到隐式梯形积分方法中，并取得较好的效果。与显式数值积分法相比其计算量要大。经分析可知，隐式梯形积分 法的误差小，局部截断

6、误差为 o(h(3),且数值稳定性好。1.2 状态空间法状态空间法是通过状态变量描述来建立系统内部状态变量与外部 输入变量和输出变量之间的关系。状态变量法属于一般性建模方法，列出的电路状态方程和输出方程组成为式中，A为系统矩阵，B为输入矩阵，C为输出矩阵，D为直接传 递矩阵，它们是由系统的结构和参数所定出的常数矩阵。典型的应用状态空间法求解电路的仿真软件有 Simulink、PLECS等。以Simulink为 例，其数值分析方法具体来分包含两大类：显式 Runge-Kutta (RK) 方法和隐式刚性解法屿。Runge-Kutta(RK)方法可以构造高阶精度的方法来求解初值，同时可以避免函数f

7、的偏导数计算。低阶显式RK法通常通过半离散化的方法来求解偏微分方程，但目 前对于系统中含有初值问题的常微分方程通常倾向于采用中高阶来解 决。针对一个系统来考虑时，稳定性要比精确度更重要。而显式 RK方 法中，高阶RK法的稳定性要比低阶稳定性差。所以对于低精度要求的 系统采用低阶方法更有效。Simulink中汇集了各种求解常微分方程数值的方法，这些方法分为 两大类：可变步长类算法和固定步长类算法。在相同的误差允许条件下，变步长计算可以降低计算次数、缩短仿真时间。另外变步长仿真计算可 以提供误差控制和过零检测，应用较为广泛。Simulink提供了 7种连续的变步长常微分方程求解器：Ode45、Od

8、e23、Ode113、Ode15s、Ode23s、Ode23t、Ode23tb。另外，Simulink中包含2种离散的常微分方程求解器：Backward Euler和Tustin。这两种算法也属于隐式一龙格库塔法。Backward Euler (后向欧拉法)是一种具有较高稳定性的算法( L- 稳定，属于A-稳定范畴)，并且由于其差分方程中不出现 t-At时刻的 非状态变量，在此过程中不会因为非状态变量的突变产生等值注入源，因此该算法不会产生数值震荡。但是，相比于隐式梯形积分法，后向欧 拉法的精度要差。Tustin算法又称梯形法或双线性变换（A-稳定）。从误差或者精度的角度考虑，Tustin算法

9、要优于欧拉法。表1各变步长求解器的适用条件4,710求解赤券型辆型翼粮别特点稳定集理中等电方卓鼾舒.大¥散情况F形可使HI一ode23低等染式，不芋求解希.近次能卡限用在在杯点刚性IJ图W件看轲岛学息式.窘!kit耕fh翻熊奥或校产的话.比3壮有败巫阶方法必旧低等制中黑宅/求髀器，Mft或用奥WI5Kx或成不耶右敬时1晚阶有也L-ptnhle.修m低¥单年求笔=1亶小二彩祕1耻加1蒯不是很有敢的刖性网典.速度快W支 防方拈Lbsiabk岫2小刚H他等何题是坐嗖刷慢.希基俊符数政康院时,前11M求解代数慵分月理 115ALi2 Hfi/yfAA-Hiihk低等比值丈杵附.求明

10、刚作战尔心程RI DRE：却办法LiuNcTable 1 Applicable conditions of solvers with variable step size4,7 101.3 状态空间节点法实时数字仿真软件 RT-LAB是在Matlab/Simulink的基础上发展而来的，通过对电路进行并行处理达到实时化的目的。对于一个没有长 输电线、节点密集、电力电子设备种类和数量众多的电力系统，利用传 统的电力系统模型解耦方法很难保证在确定的步长和精度的前提下完 成模型的解算。RT-LAB在状态空间模型的基础上，结合节点分析法的 优点，提出专门针对以状态空间方程建模的大系统解耦方法一状态空间

11、 节点法11-13。SSN建模方法的基本原理是将电路模型分割成若干任意大小的子 网络，并使用状态空间法对各子网络进行建模和计算，求解出各个子网络的等效电路。在各个子网络之间使用经典的节点法建立即时求解的节 点网络方程，并在同一个仿真步长内，对各个子网络所解算出的等效电 路进行联解得到整个电路的全局结果。这样就避免了解耦所带来的人为延时，并提高了仿真的精度和数值稳定性。另外，适当的划分状态空间群组，不但可以有效减少系统的等效电 气节点数，还能减小每个群组的开关数量及预计算量，节省存储空间， 加速仿真。由于子网络的规模较小，因此子网络内部可以采用高阶的算 法进行状态空间方程的求解，例如RT-LAB

12、中常用的ARTEMIS “art5 ” 算法。考虑到稳定性，art5法具备L-稳定性(属于A稳定范畴)，使 得art5算法能够抵抗类似梯形算法的不稳定性的产生，在消除数值振 荡方面有一定的优势。2案例分析本文中采用了带有阶跃信号控制理想开关的RLC谐振电路来验证三种平台(PSCAD/EMTDC,Simulink , RT-LAB)中不同数值算法之间 的差异，如图1所示。RLC串联谐振电路及参数设置为：串联电阻R=1 X10-4Q,串联电感L=1 X10-5 H,串联电容C=6 X10-6F,阶跃信号发生 器的设置为 Step Time 为 0.1s。EMTDC、Simulink 以及 RT-L

13、AB 中 的离散算法设置的步长均为 t=1 X10-5S, Simulink中的连续模型设置Ideal ituh理想开共-i-udc，i；波霹酒时元阱图1 RLC串联谐振电路及参数设置* LX图2 RT-LAB仿真测算模型Fig. 2 RLC model simulated by RT-LAB根据已知的电路，当理想开关闭合后，可以得到该电路的微分方程:将式（4）转化为状态空间方程可得:&SN Intel tdue BlockS1W愉岷位5电压Uc。图2为采用RT-LAB中的SSN建模方法给出的RLC模型。一 MJemyjy 二 c. Tr1 OmiM-kku reincnl的最大步长限

14、制为 t=1 x10-4s。其中示波器输出的是电容器C两端的Fig. 1 RLC series-resonant circuit and parameter setting状态空间方程离散化后的精确解为：Xi =/" / + J 泮'f&zndE为了获得离散的表达式，需要对矩阵指数e招进行估计。分别利用隐式梯形法、后退欧拉法和 art5算法对矩阵指数进行近似：(1)隐式梯形积分法十-A /& = 2(7)1-A 2(2)后退欧拉法也之一!(S)I-MA(3) art5 算法卷“，5 + (必)2。'=/ - 3AM / 5 +/ 20- (/ 60通过

15、对RLC谐振电路进行离散化后可知，基于不同的仿真算法，系统最后得到的递归方程不同，从而拥有不同的收敛速度和暂态过程。下面对不同仿真算法的仿真结果进行对比分析。3结果分析对上述电路以及对应参数列出的微分方程进行求解，可得其解的形式如下：0而卜二j+C1e-l2vmi)+100 以0”带入所测得的两个点，可求解常数C、C2的值。针对三种仿真平台中的不同算法进行设置，图 3图6分别给出了4种不同离散算法：EMTDC中的隐式积分算法、Simulink-Discrete 中 的两种状态空间法 (Tustin、Backward Euler )以及RT-LAB中的art5 算法的解算结果，其中图5插图给出了

16、局部的放大效果。测试结果表明，不同的解算方法得到的解算结果差异性较大。160 340 -下RO L0.00.20.40.6C.S1.0图3 EMTDC的解算结果Fig. 3 Simulation results of EMTDCiso1(K)0 o.d0.20.4口与。禺1 舟/三图4 Tustin的解算结果Fig. 4 Simulation results of Tustin图5 Backward Euler的解算结果Fig. 5 Simulation results of Backward Euler图6 Artemis 中的SSN求解结果Fig. 6 Simulation results

17、 with SSN solvers通过与理论推导结果进行对比，EMTDC方法（二阶隐式梯形积分法）的仿真结果和基于Simulink平台的Tustin算法（二阶隐式梯形积分法）与真解相近，且二者给出的仿真结果的形貌差别不大（电路电气 振荡频率、相位和衰减及振荡衰减后的稳态值），但振荡初始幅值存在 差异。art5算法在振荡的幅值上差异更为明显， 同时与真解相比，art5算 法能够较为迅速地衰减电路中产生的数值振荡（100 ms内原型震荡衰 减完毕），art5为一种L稳定算法，其特性会衰减较高频率的原型振荡， 此外，可能与SSN建模方法使得在同一个仿真步长内将计算分为两个 阶段计算而又未采用迭代计算

18、有一定的关系。如上所述，Backward Euler在此过程中不会因为非状态量的突变 而产生等值注入源，能够对数值振荡起到很好的抑制效果。在Simulink 仿真平台中，Ode15s、Ode23t、Ode23tb 三种解算 方法的参数设置中存在 Solver reset method 选项（Fast和Robust ）, 根据不同的设置方式，得出的结果差异性较大。其中Ode15s和Ode23tb两种算法设为Fast条件下输出的结果出现了非常大的冲击值（此处结果并未列出），与正常解算值之间差别甚大，显著异于真解， 因而报告只给 Robust选项下的结果。图7给出了 Ode15s、Ode23s 和O

19、de23tb的解算结果。明显的，三种算法均未能给出准确的解，原 型振荡明显不是按照指数规律衰减，与实际物理过程不符，算法加入了人为的衰减因子。另外，Ode23t算法与真解相近，但原型震荡初始幅值偏大。2W1500 0.0070.4060.Kt.O施图7 Simulink中三种连续模式求解器的解算结果Fig. 7 Simulation results of three kinds of continuous solvers in Simulink此外，Ode23t算法在设置Solver reset method这一参数和限制最大步长的不同设置均存在较大差异，如图8所示，其中插图给出了仿真时间为1

20、.02.0s时间段内的仿真结果输出波形。其中，图 8 (a)和 图8 (b)分别为 Max step size设置为auto时的Fast和Robust仿 真结果，图8 (c)和图8 (d) 分别为Max step size 设置为1 x10-4s 时的Fast和Robust仿真结果，1.21.3s后原型振荡不再衰减，甚至 稍有增幅并维持恒定幅值持续振荡，这是不合理的，主要是由于步长在 较大的情况下截断和累积误差较大造成。t/s(a)0.00.20.40.60£LOt/s(b)0.0ON0.40.60.81.0/s (c)00020.40.60.81.0图8不同条件设置下的 Ode23

21、t算法的仿真结果Fig. 8 Simulation results of Ode23t that the parameters was set different(d)图9低频RLC模型仿真结果Fig. 9 Simulation results of low frequency RLC circuit图9给出了低频RLC模型仿真结果。其中RLC谐振电路模型参数 为：串联电阻R=1 X10-4Q,串联电感L=1 X10-3H,串联电容C=6 X10- 4F,阶跃信号发生器的设置为 Step Time为0.1s, Max step size 设置 为 1 X10-4s。结合前面的结果分析可以看出，Tustin算法针对高低频谐振电路的 仿真结果并无明显差异，在处理低频电路时，其振荡衰减速度略低、震 荡幅值略大；Backward Euler算法、SSN算法以及Ode23tb算法在 处理高低频谐振电路的仿真时，其结果存在明显的差异：B