时域仿真法暂态稳定分析知识分享

8时域仿真法暂态稳定分析8.1弓I言电力系统暂态稳定分析的主要目的是检查系统在大扰动下(如故障、切机、切负荷、重合闸操作等情况)，各发电机组间能否保持同步运行，如果能保持同步运行，并具有可以接受的电压和频率水平，则称此电力系统在这一大扰动下是暂态稳定的。在电力系统规划、设计、运行等工作中都要进行大量的暂态稳定分析，因为系统一旦失去暂态稳定就可能造成大面积停电，给国民经济带来巨大损失。通过暂态稳定分析还可以研究和考察各种稳定措施的效果以及稳定控制的性能，因此有很大的意义。当电力系统受到大扰动时，发电机的输入机械功率和输出电磁功率失去平衡，引起转子的速度及角度的变化，各机组间发生相对摇摆，其结果可能有两种不同情况。一种情况是这种摇摆最后平息下来，系统中各发电机仍能保持同步运行，过渡到气个新的运行状态，则认为系统在此扰动下是暂态稳定的。另一种情况是这种摇摆最终使一些发电机之间的相对角度不断增大，也就是说发电机之间失去了同步，此时系统的功率及电压发生强烈的振荡，对于这种情况，我们称系统失去了暂态稳定。这时，应将失步的发电机切除并采取其他紧急措施。除此以外，系统在大扰动下还可能出现电压急剧降低而无法恢复的情况，这是另一类失去暂态稳定的形式，也应采取紧急措施恢复电压，恢复系统正常运行。这两大类暂态稳定问题分别称为功角型和电压型暂态稳定问题，并且常互相影响，互相关联。为了防止在大扰动下系统失去暂态稳定，在电力系统中需要根据预想的典型大扰动，分析系统在这些典型扰动下的暂态稳定性，这就是电力系统暂态稳定分析的基本任务，其中最大量的分析是功角稳定问题。现代电力系统一方面采用了先进技术和装置来改善系统的暂态稳定性，如快速高顶值倍数的励磁系统、快关汽门、制动电阻、静止无功补偿装置、高压直流输电技术等等；但另-方面又出现了一些对暂态稳定不利的因素，例如：大型机组参数恶化，其相应的暂态电抗xd增大和惯性时间常数Tj相对减少；超高压长距离重负荷输电线路的投入；同杆并架线路的增加等等。此外，有些措施对第一摇摆稳定有利，但对系统后续摇摆中的阻尼性能及相应的系统稳定性带来不利影响，因此要注意稳定措施的全局规划及协调。电力系统暂态稳定分析目前主要有两种方法，即时域仿真(timesimulation)法，又称逐步积分(stepbystep)法，以及直接法(directmethod)，又称暂态能量函数法(transientenergyfunctionmethod)o时域仿真法将电力系统各元件模型根据元件间拓扑关系形成全系统模型，这是一组联立的微分方程组和代数方程组，然后以稳态工况或潮流解为初值，求扰动下的数值解，即逐步求得系统状态量和代数量随时间的变化曲线，并根据发电机转子摇摆曲线来判别系统在大扰动下能否保持同步运行，即暂态稳定性。时域仿真法由于直观，可适应有几百台机、几千条线路、几千条母线的大系统，可适应各种不同的元件模型和系统故障及操作，因而得到广泛应用。本章介绍时域仿真法暂态稳定分析，而直接法暂态稳定分析在下一章中介绍。8・2简化模型时域仿真法暂态稳定分析本节采用简化的元件模型来介绍时域仿真法暂态稳定分析的基本原理、步骤，并提出采用复杂元件模型时会出现的问题。

电力系统基本上是由发电机、励磁系统、原动机及调速器以及网络和负荷组成的。其相互联系示于图8-1。其中发电机分为二部分，即转子运动方程部分和电磁回路方程部分。转子运动方程反映了当发电机输入机械功率Pm和输出电功率p不平衡时引起发电机转速«和转子角^的变化。发电机转速信号送入调速系统和参考速度比较，其偏差作为调速器的控制输入量，以控制原动机的输出机械功率Pm。发电机转子角^则用于进行发电机dq坐标下电量和网络同步坐标下电量间的接口。发电机的电磁回路方程即发电机定子、转子绕组在dq坐标下的电压方程，它以励磁系统输出励磁电压E§（文献中常用E舟）为输入量，发电机端电压和电流经坐标变换，可跟同步坐标下网络方程接口，并联立求解。所解得的机端电压U,反馈回励磁系统，励磁系统将机端电压和参考电压U浏比较，以控制发电机励磁电压Ef。而发电机的输出电磁功率［将影响转子运动的功率平衡及转子速度和角度的变化。网络一般表示为节点导纳阵形式，网络除和发电机相连外，还和负荷相连。图8-1中只画出了实际网络和一台发电机、一个负荷之间的联系。实际的电网有许多发电机和负荷，通过网络互相联系和互相影响，造成了电力系统暂态稳定分析的复杂性。原所叽及段电肌转子话功有寿原所叽及段电肌转子话功有寿图8-1电力系统基本组成部分及相互联系示意图暂态稳定分析由于主要研究发电机转子摇摆特性，主要和网络中的工频分量有关，故发电机可忽略定子暂态而采用实用模型，而网络采用准稳态模型，负荷则采用第4章所介绍的静态模型或（和）计及机械暂态或机电暂态的动态模型。为了突出电力系统暂态稳定分析的基本原理和步骤，本节对发电机采用经典二阶模型，忽略凸极效应，并设暂态电成'后的暂d态电动势E'幅值恒定，从而忽略励磁系统的动态，以简化分析。应当指出，E'恒定已计及了励磁系统的一定作用，即认为励磁系统足够强，从而能保证X'后的暂态电动势E'恒定。d另外，本节中忽略调速器和原动机动态作用，即认为机械功率Pm为定常值。在上述模型及相应假定下，发电机忽略定子绕组内阻时的定子电压标幺值方程为U=戒-jX'杈 （8-1）G dG式中，U，I为发电机端电压及流出的电流，均为同步坐标下的复数量；救=E/为GG暂态电动势，E'=const.。式（8-1）是同步坐标下的复数线性代数方程。发电机的转子运动方程为（标幺值，下同）：

式中E血。前式中E血。前="d8——=①一1、dt(8-2)P=Re(U|IG)P=const.当将式(8-1)、式(8-2)和网络方程联立求解时，可解出沸，I,叽&GG对于负荷，设采用最简单的线性负荷模型，从而对于三相对称负荷有苗=Z件 或1=y/l (8-3)式中，zlYl分别为负荷等值阻抗和导纳；乌，I分别为负荷电压及其吸收的电流。若设网络节点导纳阵方程为(8-4)YlI=1(8-4)式中，矽和闻别为节点电压和各节点注入网络的电流。对于发电机节点，I相应元为I；G对于负荷节点，I相应元为-1；对于网络节点，I相应元为零。L式(8-1)〜式(8-4)构成了系统的基本方程。这是一组联立的微分方程组和代数方程组。暂态稳定时域仿真分析的核心是当tn时刻的变量值已知时，如何求出tn+1时刻的变量值，以便由10时的变量初值(一般是潮流计算得的稳态工况下变量值)，逐步计算出t1，12,...时刻的变量值，并在系统有操作或发生故障时作适当处理。下面先介绍上述简化模型下，。〜tg时段的计算方法。对于式(8-1)〜式(8-4)的简化模型电力系统，可将负荷阻抗并入导纳阵，这只要修正负荷接点对应的导纳阵对角元，从而负荷接点转化为网络节点，式(8-4)中相应节点的注入电流化为零。同时将各发电机方程(8-1)改写为导纳方程形式，即(8-5)1=芝一生=Y建一YI1塑1一YU(8-5)G jX'jX' GGGGGGd d式中，Y=-L.，为发电机暂态导纳，式(8-5)的等值电路如图8-2所示。显然可把Y并GjXr Gd入网络导纳阵，即修正发电机节点相对应的导纳阵对角元，则联立求解发电机和网络方程的问题就转化为在发电机节点有注入电流I=YgE时，网络方程(已将Yg和yl并入导纳阵)的求解问题。网络方程的求解本质上是求解一组复数线性代数方程，可用高斯消去法。由于系统无操作时，导纳阵不变，故可预先对导纳阵作三角分解，存储因子表，然后每一时步根据各节点注入的电流求解各节点电压。在计算每一时步各发电机的等值注入电流戊时，由G

于E的相角部随时间而变，需由转子运动方程计算确定，故实用中可根据o时刻的先用某种微分方程的数值求解法来估算t时刻的①和&，如设n+1 n+1 n+1由式(8-2)取q①+c&q①+c&h=3+h(P-P)/T①n+1 n8»8+n+1 nnn①+①"2"+1(8-6)式(8-6)又称作是微分方程(8-2)在t〜t时段上的差分代数方程，从而可得E=EM，nn+1 n+1 n+1则各发电机t时刻等值电流源I 可求，可进而求解网络方程得洪，然后可根据式(8-5)n+1 G,n+1 n+1计算发电机端电流I，并计算发电机的电磁功率P广Re吨切。这样计算得的tn+1时刻的变量精度可能较差，必要时可进行校正和迭代计算，以改善精度。图8-2经典模型发电机等值电路图简化模型的电力系统暂态稳定分析的步骤和流程框图见图8-3。下面对其作简要说明。暂态稳定分析首先输入原始数据，这包括系统元件的模型、参数、网络拓扑信息、扰动过程信息、稳定分析要求(如计算步长、仿真总时间、失稳判据等)、打印输出要求，另外还应输入暂态分析的初始稳态工况，一般为潮流计算结果。此即流程框图中框①。然后根据潮流及原始数据计算各代数变量和状态变量的初值，及E'和Pm的稳态值，采用简化模型时E'和P在暂态过程中保持不变。此即流程框图中框②。对于负荷节点，m潮流中已计算得负荷有功功率PL0、无功功率Ql0、及负荷母线电压乌。，则由PL0+问L0=匕肾°F (8-7)可计算负荷等值导纳yl。对于发电机节点，潮流中已计算得发电机发出的有功P°、无功功率QG°及端电压UG°，则由PG0-jQG0=吃1G° (8-8)计算I，再由式(8-1)计算E0=EN8°=[I°+jX'd吃，得E'及8°，电磁功率P°=时吃IG0)=Pm°，E'和Pm°在暂态中保持不变。此外3°=1(p.u)，至此系统暂态分析的初值计算毕。根据网络元件参数及网络拓扑关系形成网络稳态工况下节点导纳阵，也可直接从潮流输出中读入。将负荷等值导纳七及发电机内部暂态导纳K=,并入导纳阵，对导纳L GJX‘d阵作因子表计算。此即流程框图框③。将时钟指针tn置零，根据扰动过程参数，判别当前有无扰动发生。若有扰动则需要根据扰动参数修改导纳阵及微分方程，并设tn时刻状态量不突变，据扰动后系统代数方程计算t时刻的代数量，作为t：~、时步的初值，此即流程框图中框④和⑤；若tn时刻无扰动则直接转入框⑥。此流程框图另作文件处理图8-3简化模型暂态稳定分析流程框图作t：~tn+1时段计算，求取tn+1时刻的状态量和代数量，前面对此已作介绍，不予重复。此即框⑥的工作。若本时步末要求打印输出结果，则转框⑦作相应处理，否则判别是否要停机：包括由于仿真总时间到而要求停机及据失步判据已判明系统失步不必继续计算而停机。若要停机则作相应结尾处理而停机，否则表明系统还应继续仿真下去，则更新时标，转去下一步计算。此即流程框图中框⑧一⑩的工作。下面对实际的暂态稳定分析中的主要问题作一初步讨论，以便在后续章节中加以解决。发电机凸极效应和采用高阶模型时的问题当计及发电机凸极效应时，X'。X'，dq因此定子电压方程不能表示为与经典模型相似的同步坐标下的复数方秋8-1)，而需分别建立定子d绕组、q绕组的电压方程，并在联网时作特殊处理，这包括凸极效应处理和dq-xy坐标变换。此外发电机采用三阶及三阶以上实用模型时，要计及励磁系统动态，需将发电机和励磁系统微分方程联立作数值计算。当进一步计及原动机和调速器动态时，还要加入其相应的微分方程一起作数值计算。负荷采用非线性静态模型或动态模型时的问题当负荷采用非线性静态模型时，在联网计算中需要求解非线性代数方程组，从而增加了分析计算的复杂性。实用计算时，要对负荷和网络的接口作特殊处理，以便计算各时段的网络潮流。当负荷采用动态模型时，联网计算需将其微分方程差分化，化为相应计算时步的差分代数方程，再和网络方程联立求解，动态负荷和网络接口时也要作适当处理。微分方程求数值解的数值稳定性问题暂态稳定分析要求求解联立的微分方程组和代数方程组，对于tn-tn+i时步计算通常将微分方程根据某种数值积分准则或根据泰勒级数化为差分代数方程，从而由。及过去时刻的系统变量求取t川时刻的状态量和代数量。若采用不同的数值积分方法(如改进尤拉法、龙格-库塔法、隐式梯形积分法等等)，则数值积分误差的传递规律，或者说数值稳定性将不同，有些方法在一定条件下会使分析结果严重畸变。此外，采用不同的数值积分方法还会影响计算的处理过程以及计算的精度和时间。微分方程和代数方程交替求解时的“交接误差”问题在求解系统的微分方程组和代数方程组时，有些算法对微分方程和代数方程交替求解，即对于系统方程组(8-9)式中，y为状态矢量；z为代数矢量；/、g为适当维数的函数。若设tn时刻的yn和气已解出，并据式(8-9)的第二式，用某种数值积分法估计状态量丁在t时刻的值y(0)，再将y(0)TOC\o"1-5"\h\zn+1 n+1 n+1代入式(8-9)的第一式通过求解代数方程计算z(0)，这样求得的y(0)和z(0)一般不能严格满足n+1 n+1 n+1式(8-9)的第二式。为改善精度，可进一步根据yn，气,yn%，zn+1和式(8-9)第二式，作yn+1的校正计算，得校正后的y，然后再代入式(8-9)第一式计算与校正后的y相对应的n+1 n+1z，如此迭代直到计算收敛。显然这种计算方法对代数方程和微分方程交替求解，计算结n+1果不能同时满足式(8-9)中的两组方程，从而造成所谓的“交接误差”，若多次迭代又会增加机时。为了消除“交接误差”，可把式(8-9)中的代数方程和差分化的系统微分方程联立求解，但求解过程较复杂，因为一般要求解一组非线性代数方程组。故障及操作的处理问题当系统发生故障或操作(切机、切负荷、切除线路等等)时，系统节点导纳阵和微分方程组要作相应的修正。由于系统状态量在过程中不发生突变，而代数量则在操作瞬间要发生突变，故还要根据操作后的系统代数方程求解突变后的代数量。还应指出，在发生不对称操作和故障时，还要根据序网理论和故障分析理论作相应处理，在复杂故障时，处理更为复杂。下面几节将分别介绍上述问题的实际处理方法，然后再在此基础上介绍采用复杂元件模型时系统暂态稳定分析的典型计算方法和步骤。8.3发电机节点的处理和机网接口计算发电机节点的处理和机网接口计算与发电机采用的模型有关，也和联网计算采用的方法有关。目前的发电机节点处理方法大体上可分为以下4类。发电机采用经典模型时的处理方法。这一类处理方法已在上一节简化模型暂态稳定分析中作了介绍，即化为图8-2所示发电机等值导纳K二二和发电机等值电流源GJX‘d%=yE相并联的形式。将Yg并入导纳阵，无操作时导纳阵不变，而每时步据发电机转子角^更新发电机注入网络的等值电流源戊，即可求解网络方程，计算节点电压。G考虑凸极效应的直接解法。其实质是将网络复数线性代数方程实、虚部分开，增阶化为打同步坐标下的实数线性代数方程，并将发电机方程由dq坐标化为打坐标，再和网络方程联立求解，最终是在实数域内求解线性代数方程。这种解法对负荷非线性适应能力差，且发电机方程由dq坐标据转子角^转化为坐标，引起导纳阵中发电机节点相应的对角(2X2)子块由于^变化而为非定常元素，每一时步要重新计算因子表，机时多且内存要增加一倍，目前在实用的暂态稳定分析程序中已不再采用此法，但这种方法物理概念清楚，不需迭代，求解网络方程为求解实线性代数方程组，也有一定优点。后面将对之作简单介绍。考虑凸极效应的迭代解法。该方法特点是力求在复数域中求解线性代数方程来实现网络方程求解，并要求导纳阵元素在无操作时保持定常，而不随发电机转子角^而变化，从而克服了直接解法的缺点。但发电机的凸极效应及转子角变化对机网接口计算的影响，都要通过修正发电机注入网络的电流源来计及，而电流源计算还同t时刻的节点电压值有关。n+1

由于t时刻的节点电压正待计算，而不预知，因此t时刻相应的电流源也不能预先准确n+1 n+1计算，故要通过迭代，逐步逼近准确值，这就是迭代解法的本质。迭代解法相对于直接解法有节省内存、因子表定常、计算速度快、便于适应非线性负荷模型等特点，但计算中每时步计算需要迭代，并有迭代误差。迭代解法在一些实用暂态稳定分析程序中仍在使用，并常和改进欧拉法求解微分方程相结合。后面将对此方法作进一步介绍。考虑凸极效应的牛顿法。牛顿法是求解非线性代数方程组的优良方法，有良好的收敛性能，已广泛用于电力系统潮流计算。当发电机计及凸极效应，负荷计及非线性，系统中元件微分方程化为差分代数方程后，全网的代数方程联立，实质上是要求解一组非线性代数方程，故也可采用牛顿法求解。相对于直接解法和迭代解法，用牛顿法进行机网接口计算编程复杂，因为要计算雅可比矩阵元素，而雅可比矩阵元素随时间而变化，故计算机时也较多。但其最大优点是对非线性元件模型的适应性好，可将微分方程的差分代数方程和系统代数方程联立求解，无“交接误差”，故计算精度高、累计误差小，因而在暂态稳定分析中广泛应用。它常和隐式梯形积分法求解微分方程相结合，后面将对隐式梯形积分法作进一步介绍。8.3.1考虑凸极效应的直接解法当发电机计及暂态凸极效应，即X'。X'时，发电机定子方程就不能用式(8-1)的简单dq复数关系来表示，必须对d轴、q轴等值绕组分别列方程。当发电机采用四阶(或三阶)实用模型时，定子电压方程为(三阶模型时，下列方程中E,=0,xq为Xq)（8-10）（8-11）\U=E'+X'I-rI（8-10）（8-11）d d qq ad\U=E'-X'I-rIq q dd aq当发电机采用五阶、六阶实用模型时，定子电压方程为U=E"+X〃I-rId d qq adU=E"-X"I-rIq q dd aq由于式(8-10)和式(8-11)有相同形式，故下面将以发电机三阶、四阶实用模型为例，即用式(8-10)讨论机网接口计算问题。将式(8-10)写成导纳阵形式IdId=ra-X广rqa-1E'd-Ud-1raX广rqaE'dIq」X，dE，q-Udr2+X，X，-X'dE'q（8-12）为了和网络方程接口，需将dq坐标化为xy同步坐标，对式(8-12)二边右乘坐标变换阵sin6cos6『：T=-cos6sin6，则fy=Tfq，从而式（8-12）化为rI]1rI]1rrX，]「E'-U]def「GB]「Ee-U]x=T aqT-1xxxxIfxIr2+X'X'-X'rE，-UBG1EE-Uya dqdaLyyyyyy式中(8-13)E'=E'sin8+E'cos8xd qE'=-E'cos8+E'sin8yd1qG= [r-(X'-X')sin28/2]X r2+X'X'adqa dqj-1 1(8-14)JB= (X'+X')+(X'-X')cos28(8-14)x r2+X'X'\_2dqd/ _G= 1 [r+(X'-X')sin28/2]yr2+X'X'adqa1dq「1B= -—(X'+X')+(X'-X')cos28/2yr2+X，X' 2dqdd，显然式(8-13)中导纳阵GXByB显然式(8-13)中导纳阵GXByG是转子角8的函数，且Gx^G,Bx」-By，不具备yG-BD厂形式，无法直接将式(8-13)化为复数方程，然后与网络方程联立，在复数域中求BG解。为了便于机网接口，直接解法中先把n个节点的网络复数线性代数方程仲=酉曾阶化为2n维的实线性代数方程。「I_Ix1y」1「「G11LB11-B-11GMAM「G_1iLB1i-B1iGM1’AM「G1nLB1n-B-1nG1nM1_U一x1U,My11「i〕「G-B一「G-B一「G-B一_「U]xi——i1AiiiiAininxiIi1LBGLBGLBG」__UL对i1Mi1MiiMiiMinMin_yiM「I-「G-B-「G-B-「G-B-_「U-xnn1AniniAnnnnxnIn1LBGLBGLBG__Uyn_1—n1n1nininnnn_Lyn」式中，G+jB=Y为Y阵中i行j列元素，/+jI=我和U+jU=或分别为院ijijij xiyiixi yii中第i个元素。为不失一般性，设式(8-13)所描写的发电机接于网络第i个节点，则式(8-13)中I和y「U]「I]「U]「I]Ux即为式(8-15)中和U，将式(8-13)代入式(8-15)中第i个节点方程，消去yLyi」Lyi-1■yi并将i节点方程整理为I'.defGBEyG-BUG+G-B+BUxi—xxx—ijijxj+iixiixxiI'BGE'BGUB+BG+GU.yiJyyyj=1ijijyjiiyiiyyij归(8-16)显然，当根据系统微分方程预估本计算时步末的E',E',8等状态量的值后，则可根据dq

式(8-14)计算G,G,B式(8-14)计算G,G,B,B及E',E'；然后根据式(8-16)可计算发电机注入网络的等值电流7一GB_源IX；再用发电机2X2等值导纳阵xB_xGyiyyxyxyxy根据式(8-16)修正网络导纳阵i节点相应的2X2子块GiiBii-BiiGii，并将式(8-16)代替式(8-15)中的第i号节点方程，对各个发电机节点均作相似处理后，便可求解网络方程。这种机网接口求解方法称为直接解法，其主要优点是物理概念清楚，不需迭代。网络方程求解为求解一组实线性代数方程。其主要问题是增阶处理使内存要增加一倍；G,B,G,B非定常，从而网络方程系数矩阵经修正后也非定常，因此每一时步要作三角分解，计算量较大；另外，这种方法对非线性负荷适应性略差。因此，这种解法目前逐步为下面介绍的方法所取代。8.3.2考虑凸极效应的迭代解法设发电机采用3〜6阶实用模型，与直接解法相似以式(8-10)为基础进行讨论，将之化为式(8-13)和式(8-14)表达的导纳参数形式。观察式(8-14)可知，Gx,B,G,By中的定常部分r -(X'+X')具有形式，其中G= a,B=一d q。故可用复数形式G+jB具有r2+XdX "XdX；)下面作具体推导。表示；而Gx,Bx,G,By中随8变化部分，可同电动势一起用非定常的电流源表示，从而引出相应的联网迭代解法。下面作具体推导。将式(8-13)改写为r11x=1ydefgr11x=1ydefg-BirE’x

BG-IEry-Ux-Uy1 . .c、a+2 d qr2+X'X'dq-(