时域仿真法暂态稳定分析

电力系统暂态稳定分析的主要目的是检查系统在大扰动下(如故障、切机、切负荷、重合闸操作等情况)，各发电机组间能否保持同步运行，如果能保持同步运行，并具有可以接受的电压和频率水平，则称此电力系统在这一大扰动下是暂态稳定的。

在电力系统规划、设计、运行等工作中都要进行大量的暂态稳定分析，因为系统一旦失去暂态稳定就可能造成大面积停电，给国民经济带来巨大损失。

通过暂态稳定分析还可以研究和考察各种稳定措施的效果以及稳定控制的性能，因此有很大的意义。

当电力系统受到大扰动时，发电机的输入机械功率和输出电磁功率失去平衡，引起转子的速度及角度的变化，各机组间发生相对摇摆，其结果可能有两种不同情况。

一种情况是这种摇摆最后平息下来，系统中各发电机仍能保持同步运行，过渡到新的运行状态，则认为系统在此扰动下是暂态稳定的。

另一种情况是这种摇摆最终使一些发电机之间的相对角度不断增大，也就是说发电机之间失去了同步，此时系统的功率及电压发生强烈的振荡，对于这种情况，我们称系统失去了暂态稳定。

这时，应将失步的发电机切除并采取其他紧急措施。

除此以外，系统在大扰动下还可能出现电压急剧降低而无法恢复的情况，这是另一类失去暂态稳定的形式，也应采取紧急措施恢复电压，恢复系统正常运行。

这两大类暂态稳定问题分别称为功角型和电压型暂态稳定问题，并且常互相影响，互相关联。

为了防止在大扰动下系统失去暂态稳定，在电力系统中需要根据预想的典型大扰动，分析系统在这些典型扰动下的暂态稳定性，这就是电力系统暂态稳定分析的基本任务，其中最大量的分析是功角稳定问题。

现代电力系统一方面采用了先进技术和装置来改善系统的暂态稳定性，如快速高顶值倍数的励磁系统（供给同步发电机励磁电流的电源及其附属设备统称为励磁系统。

它一般由励磁功率单元和励磁调节器两个主要部分组成。

励磁功率单元向同步发电机转子提供励磁电流；而励磁调节器则根据输入信号和给定的调节准则控制励磁功率单元的输出）、快关汽门(实现两种减功率方式:短暂减功率和持续减功率1、短暂减功率用于系统故障初始的暂态过程,减少扰动引起的发电机转子过剩动能以防止系统暂态稳定破坏。

2、持续减功率用于防止系统静稳定破坏、消除失步状态、限制设备过负荷和限制频率升高。

)、制动电阻（1、保护变频器电机在快速停止过程中，由于惯性作用，会产生大量的再生电能，如果不及时消耗掉这部分再生电能，就会直接作用于变频器专用型制动电阻，变频器的直流电路部分，轻者，变频器会报故障，重者，则会损害变频器； 2、保证电源电网络，制动电阻将电机快速制动过程中的再生电能直接转化为热能，这样再生电能就不会反馈到电源电网络中，不会造成电网电压波动，从而起到了保证电源网络的平稳运行的作用。

）、静止无功补偿装置、高压直流输电技术等等；但另一方面又出现了一些对暂态稳定不利的因素，例如：大型机组参数恶化，其相应的暂态电抗dX '增大和惯性时间常数J T 相对减少；超高压长距离重负荷输电线路的投入；同杆并架线路的增加等等。

此外，有些措施对第一摇摆稳定有利，但对系统后续摇摆中的阻尼性能及相应的系统稳定性带来不利影响，因此要注意稳定措施的全局规划及协调。

电力系统暂态稳定分析目前主要有两种方法，即时域仿真(time simulation)法，又称逐步积分(step by step)法，以及直接法(direct method)，又称暂态能量函数法(transient energy function method)。

时域仿真法将电力系统各元件模型根据元件间拓扑关系形成全系统模型，这是一组联立的微分方程组和代数方程组，然后以稳态工况或潮流解为初值，求扰动下的数值解，即逐步求得系统状态量和代数量随时间的变化曲线，并根据发电机转子摇摆曲线来判别系统在大扰动下能否保持同步运行，即暂态稳定性。

时域仿真法由于直观，可适应有几百台机、几千条线路、几千条母线的大系统，可适应各种不同的元件模型和系统故障及操作，因而得到广泛应用。

电力系统基本上是由发电机、励磁系统、原动机及调速器以及网络和负荷组成的。

其相互联系示于图8-1。

其中发电机分为二部分，即转子运动方程部分和电磁回路方程部分。

转子运动方程反映了当发电机输入机械功率m P 和输出电功率e P 不平衡时引起发电机转速ω和转子角δ的变化。

发电机转速信号送入调速系统和参考速度比较，其偏差作为调速器的控制输入量，以控制原动机的输出机械功率m P 。

发电机转子角δ则用于进行发电机dq 坐标下电量和网络xy 同步坐标下电量间的接口。

发电机的电磁回路方程即发电机定子、转子绕组在dq 坐标下的电压方程，它以励磁系统输出励磁电压f E 为输入量，发电机端电压和电流经坐标变换，可跟同步坐标下网络方程接口，并联立求解。

所解得的机端电压t U 反馈回励磁系统，励磁系统将机端电压和参考电压ref U 比较，以控制发电机励磁电压f E 。

而发电机的输出电磁功率e P 将影响转子运动的功率平衡及转子速度和角度的变化。

网络一般表示为节点导纳阵形式，网络除和发电机相连外，还和负荷相连。

图8-1中只画出了实际网络和一台发电机、一个负荷之间的联系。

实际的电网有许多发电机和负荷，通过网络互相联系和互相影响，造成了电力系统暂态稳定分析的复杂性。

图8-1 电力系统基本组成部分及相互联系示意图暂态稳定分析由于主要研究发电机转子摇摆特性，主要和网络中的工频分量有关，故发电机可忽略定子暂态而采用实用模型，而网络采用准稳态模型，负荷则采用静态模型或机械暂态或机电暂态的动态模型。

为了突出电力系统暂态稳定分析的基本原理和步骤，对发电机采用经典二阶模型，忽略凸极效应，并设暂态电抗d X '后的暂态电动势E 'E 幅值恒定，从而忽略励磁系统的动态，以简化分析。

应当指出，E '恒定已计及了励磁系统的一定作用，即认为励磁系统足够强，从而能保证d X '后的暂态电动势E '恒定。

另外，忽略调速器和原动机动态作用，即认为机械功率m P 为定常值。

在上述模型及相应假定下，发电机忽略定子绕组内阻时的定子电压标幺值方程为G d G I X j E U '-'= （8-1）式中，G U ，G I 为发电机端电压及流出的电流，均为同步坐标下的复数量；δ∠'='E E 为暂态电动势，E '=const.。

式(8-1)是同步坐标下的复数线性代数方程。

发电机的转子运动方程为(标幺值)：⎪⎩⎪⎨⎧-=-= 1ωδωdt d P P dt d T em J （8-2）式中) Re(*G Ge I U P = m P =const.当将式(8-1)、式(8-2)和网络方程联立求解时，可解出G U ,G I ,ω,δ。

对于负荷，设采用最简单的线性负荷模型，从而对于三相对称负荷有 L L L I Z U = 或 L L L U Y I =（8-3）式中，L L Y Z ，分别为负荷等值阻抗和导纳；L L I U ，分别为负荷电压及其吸收的电流。

若设网络节点导纳阵方程为I U Y =（8-4)式中，I U 和分别为节点电压和各节点注入网络的电流。

对于发电机节点，I相应元为G I ；对于负荷节点，I相应元为L I -；对于网络节点，I 相应元为零。

式(8-1)~式(8-4)构成了系统的基本方程。

这是一组联立的微分方程组和代数方程组。

暂态稳定时域仿真分析的核心是当n t 时刻的变量值已知时，如何求出1+n t 时刻的变量值，以便由0t 时的变量初值(一般是潮流计算得的稳态工况下变量值)，逐步计算出,,21t t …时刻的变量值，并在系统有操作或发生故障时作适当处理。

下面先介绍上述简化模型下，n t ~1+n t 时段的计算方法。

对于式(8-1)~式(8-4)的简化模型电力系统，可将负荷阻抗并入导纳阵，这只要修正负荷接点对应的导纳阵对角元，从而负荷接点转化为网络节点，式(8-4)中相应节点的注入电流化为零。

同时将各发电机方程(8-1)改写为导纳方程形式，即G G G G G G d G d G U Y I U Y E Y X j U X j E I -'-'='-''==def (8-5)式中，d G X jY '=1，为发电机暂态导纳，式（8-5）的等值电路如图8-2所示。

显然可把G Y 并入网络导纳阵，即修正发电机节点相对应的导纳阵对角元，则联立求解发电机和网络方程的问题就转化为在发电机节点有注入电流E Y I G G '=' 时，网络方程(已将G Y 和L Y 并入导纳阵)的求解问题。

网络方程的求解本质上是求解一组复数线性代数方程，可用高斯消去法。

由于系统无操作时，导纳阵不变，故可预先对导纳阵作三角分解，存储因子表，然后每一时步根据各节点注入的电流求解各节点电压。

在计算每一时步各发电机的等值注入电流G I ' 时，由于E' 的相角δ随时间而变，需由转子运动方程计算确定，故实用中可根据n t 时刻的n n n e n m P P δω,,,,,，先用某种微分方程的数值求解法来估算1+n t 时刻的1+n ω和1+n δ，如设nn t t h -=+1，由式(8-2)取 ⎪⎩⎪⎨⎧⎥⎦⎤⎢⎣⎡-++≈-+=+≈+++ 12/)(11,,1h T P P h h n n n n J n e n m n n n n ωωδδωωωω （8-6）式(8-6)又称作是微分方程(8-2)在n t ~1+n t 时段上的差分代数方程，从而可得11++∠'='n n E E δ ，则各发电机1+n t 时刻等值电流源1,+'n G I 可求，可进而求解网络方程得1+n U ，然后可根据式(8-5)计算发电机端电流G I ，并计算发电机的电磁功率) Re(G *I U P G e =。

这样计算得的1+n t 时刻的变量精度可能较差，必要时可进行校正和迭代计算，以改善精度。

图8-2 经典模型发电机等值电路图简化模型的电力系统暂态稳定分析的步骤和流程框图见图8-3。

下面对其作简要说明。

(1) 暂态稳定分析首先输入原始数据，这包括系统元件的模型、参数、网络拓扑信息、扰动过程信息、稳定分析要求(如计算步长、仿真总时间、失稳判据等)、打印输出要求，另外还应输入暂态分析的初始稳态工况，一般为潮流计算结果。

此即流程框图中框①。

(2) 然后根据潮流及原始数据计算各代数变量和状态变量的初值，及E 'E 和m P 的稳态值，采用简化模型时E '和m P 在暂态过程中保持不变。

此即流程框图中框②。

对于负荷节点，潮流中已计算得负荷有功功率0L P 、无功功率0L Q 、及负荷母线电压0L U ，则由2000L L L L U Y jQ P =+ (8-7)可计算负荷等值导纳L Y 。