蛙跳算法在暂态稳定性计算中应用.docVIP

蛙跳算法在暂态稳定性计算中应用 　　摘 要：电力系统暂态稳定性分析是解决电力系统稳定问题的基础。数值积分方法是进行电力系统暂态稳定性分析的一类基本和比较可靠的方法。如何实现暂态稳定性的实时和超实时分析计算，是现代大规模电网实时分析与控制研究的重要课题。文章提出的蛙跳算法易于执行，与同阶的传统RK法相比，其计算过程更为简捷。 关键词：暂态稳定性分析；数值积分方法；蛙跳算法 1 蛙跳算法简介[1][2] 如果Hamilton系统可以写成如下的形式： （1） 那么称该系统是一个可分的Hamilton系统，微分方程组（1）写成如下形式： （2） 即 （3） 对于形如式（3）的Hamilton系统，Neri已经提出了一类比较好的构造辛差分格式的方法。该方法非常简单易懂，但实际推导非常复杂，在构造6阶辛差分格式时，就几乎不可能实现了，至少难度很大。针对以上问题，Haruo Yoshida提出了一种构造高阶显式辛差分格式的方法。 令 Hamilton系统方程如下： （4） 定义微分算子DG如下： 则Hamilton系统方程（4）可以写成： （5） 在t=？子时，式（5）的精确解为 （6） 又因为 （7） 我们可以将式（6）写成 （8） 如果存在实数ci，di，i=1，2，…，k使得 （9） 令近似解为 （10） 那么 （11） 那么近似解（10）具有 阶精度，且式（10）中的exp（ci？子DT）、exp（di？子DV）均为辛变换，这样就得到一个n阶辛差分格式（10） 可分Hamilton系统n阶显式辛格式的数学表达式如下： （12） 要想构造可分Hamilton系统的显式辛格式，关键在于确定系数ci，di，i=1，2，…，n的值，Haruo Yoshida已经给出了确定系数的方法。下面给出1～4阶显式辛格式的系数。 令 当n=1时，有 c=1，d=1 当n=2时，有 或 当n=3时，有 或 当n=4时，有 或 其中 上述显式辛算法易于执行，与同阶的传统RK法相比，其计算过程更为简捷，但此方法仅适合于可分系统。 2 蛙跳算法的暂态稳定计算公式研究[3][4][5][6] 若发电机采用经典模型，负荷采用恒阻抗表示，采取简化处理后，便可作出全系统的等值电路，这将是一个多电势源的线性网络.此线性网络的导纳型节点方程为 （13） 式中， 是各发电机输出电流的列向量； 是各发电机电势的列向量。 YG是仅保留发电机电势源节点和参考节点，而其他节点经过网络变换全部消去后的等值网络的节点导纳矩阵，YG可由潮流计算用的节点导纳矩阵修改后得到。 在潮流计算中是以发电机端点作为发电机节点。现在应在每一台发电机节点i后面，通过发电机内阻抗ZGi支路，增加一个电势源节点i’，其注入电流IGi等于原来发电机节点的注入电流，而原发电机节点i的注入电流则等于零。接入ZGi和增加节点i’后，应对原潮流计算用的节点导纳矩阵进行修改，发电机有几台，修改后的导纳矩阵将增加几阶。 对于负荷节点，当负荷用恒定阻抗表示时，可在负荷节点并联接入负荷的等值阻抗，并令原负荷节点注入电流■k=0。由于负荷阻抗是并联接在负荷节点与参考点之间的，所以网络的节点数不增加，但导纳矩阵的负荷节点的自导纳应变为Y’kk=Ykk+YLk。 如果原网络有N个节点，其中发电机节点有n个，则作了上述修改后的导纳矩阵将有N+n阶。若对节点重新编号，把发电机电势源节点编为1，2，…，n号，则其余n+1，n+2，…，n+N号节点都成了无注入电流的节点。将修改后的导纳矩阵分块，节点方程可写成： （14） 根据方程式消元的方法，将上式展开，消去VN，即可求得（13），其中 （15） 将式（13）代入到发电机功率算式■Gi=PGi+JQGi=■GiI*Gi中，经整理后可得 （16） 式中 （17） 系统中每一台发电机的摇摆方程可写成如下形式： （18） 式中Pei为每台发电机的电磁功率，即式（16）中的PGi，将式（16）和式（17）代入上式，可将发电机摇摆方程写成如下的可分形式： （19） 式中 其中Mi为第i台发电机的惯性时间常数；Pmi为第i台发电机的机械功率；Gij为收缩导纳矩阵的电导；Bij为收缩导纳矩阵的电纳；Ei为第i台发电机的暂态电动势。 若系统中有m台发电机，记 其中 则发电机摇摆方程可写成 （20） 由式（3）和（1