电科院PSD-BPA稳定程序简介

1、中国电力科学研究院系统分院/系统所/软件室2008年3月杭州,psd-bpa暂态稳定程序,简介,规模：10000节点算法：微分方程采用隐式梯形积分法；代数方法采用三角分解法；微分方程和代数方程交替求解；,程序流程示意图,第一部分,主要模型,发电机模型1/7示意图,x、y为阻尼绕组f为励磁绕组g为由于q轴铁芯内部构成的短路阻尼线圈,发电机模型2/7方程,发电机模型3/7类型,考虑阻尼绕组，tq00，xq不等于xq（隐极机详细模型，6阶）考虑阻尼绕组，tq0=0，xq等于xq（凸极机详细模型，5阶）没有阻尼绕组，tq00，xq不等于xq（隐极机详细模型，4阶）；没有阻尼绕组，tq0=0，xq等于x

2、q（凸极机详细模型，3阶）没有阻尼绕组，td0=999，eq恒定模型经典模型，xq不等于xq（隐极机详细模型，4阶）；没有阻尼绕组，td0=999，eq恒定模型经典模型，xq等于xq（凸极机详细模型，3阶）e恒定的经典模型；无穷大母线，emws=999999；,发电机模型4/7类型列表,发电机模型5/7数据卡,每个发电机必须填写一张mf卡（经典模型填写mc卡）；考虑阻尼绕组，每台发电机填写m卡；在case卡中填写通用的阻尼绕组参数；不填写发电机数据卡，必须填写ln卡，作为负荷处理。,发电机模型6/7mf/mc卡,识别码id和功率比例（多台连同一母线用）发电机动能emws参数参数的基准容量mva

3、base(缺省系统基准容量)ra、xd、xq暂态电抗xd、xq暂态时间常数tq0、td0饱和相关参数xl、sg1.0、sg1.2阻尼转矩d,发电机模型7/7m/case卡,m卡填写x“d、x”q、t“d0、t”q0，基准容量为对应mf卡中的基准容量；case卡填写通用阻尼绕组参数次暂态、暂态电抗比值xfact汽轮机t“d0、t”q0水轮机t“d0、t”q0根据tq0判断使用的时间常数没有m卡的发电机采用此参数，f1卡26可控制所有都采用此参数,励磁模型（31个）,1968年ieee励磁模型ea/eb/ec/ed/ee/ef/eg/ej/ek1981年ieee励磁模型fa/fb/fc/fd/fe

4、/ff/fg/fh/fj/fk/fl电机工程学会励磁模型fm/fn/fo/fp/fq/fr/fs/ft/fu/fvfx,pss模型,sf/sp/ss/sgsh、sh/sh+si/si+st,原动机、调速器模型,调速器和原动机结合在一起的模型gg/gh/gc单独的液压调速器模型gs/gl/gw单独的电调模型ga+gigi+gjgk、gd/gz单独的原动机模型ta/tb/tc/td/te/tf锅炉的主汽压力变化模型gx,风电系统模型,固定转速风电系统模型变速恒频风电系统模型,直流系统模型,两端直流详细模型，该模型可以采用定熄弧角控制和余裕开关，以及各种调制方式和直流故障；两端直流系统简化模型，假定

5、直流端是一个理想系统。两端直流换相失败直流调制，包括大方式、小方式、双侧频率调制、逆变侧熄弧角调制新增cigre的直流系统模型和南方电网的直流模型；,负荷模型-1,静态负荷模型，有恒功率、恒电流、恒阻抗和与频率变化相关的部分组成动态负荷模型，即感应马达考虑配电网的综合负荷模型,负荷模型-2,p1p2p3p41，q1q2q3q41可以按照区域、分区或节点填写，优先级顺序：节点、分区、区域若综合负荷中有一定比值的异步马达者，剩余比值的负荷必须填负荷卡恒功率和恒电流负荷在低电压下将转为恒定阻抗负荷，计算控制继续卡f1卡的第2831列可以指定该电压值，缺省值为0.5pu,负荷模型-3,新静态负荷模型（

6、ieee模型及变种）p1p2p3p4+p51，q1q2q3q4+q51需要两张数据卡，即la或lb卡及其继续卡l+卡，l+卡必须紧接于la或lb卡之后,负荷模型-4,感应马达模型mi卡机械转矩系数，a、b用卡填写，c由初始化得出初始化之后s0在0.0010.2之间，功率因素在0.50.98之间在稳定数据文件中，马达卡应填写在发电机相关数据卡后面马达滑差达到1时的处理方式0稳定程序停止计算1堵转（缺省）2切除马达低压释放的电压值vi和ti：当马达机端电压降到vi以下，经ti延时后切除马达,负荷模型-5,新感应马达模型（ml、mj、mk）填写tj（惯性时间常数）、pper（马达功率占母线功率的百分

7、比）、以及kl（初始滑差s0下的负载率）或初始滑差s0分区和区域负荷填写pmin形成马达的最低负荷功率值（不适用于ml卡）如果原马达卡mi卡、分区马达卡mj卡、区域马达卡mk卡、mb卡和单个节点形式ml卡同时存在，优先级为mi卡、ml卡、mb卡、mj卡、mk卡。不含有马达的节点数据卡mb卡对mj和mk指定分区或者区域的指定节点不形成马达。也可以指定初始滑差：第79列填写字符“s”，26-29列填写初始滑差值，格式为f4.4。填写的初始滑差s0必须在0.001和0.2之间马达动能计算方法,负荷模型-6,负荷模型-7,考虑配电网支路的综合负荷模型（le）可以按照区域、分区或节点填写rd配电网支路电

8、阻，xd配电网支路电抗pper：本支路有功功率占该节点总负荷有功功率的比例，缺省为1qper：本支路无功功率占该节点总负荷无功功率的比例，缺省同pperpg、qg：与本支路相连的所有发电机的有功和无功总出力pfac，配电网支路静态负荷的功率因数，剩下的是无功补偿容量同一个负荷节点可以有多个配电网支路综合负荷模型数据卡，各个配电网支路之间采用le卡中17列的支路顺序号区分与配电网支路相连的静态负荷模型、马达卡，必须在lalbmimlmjmk中第3列填写对应的配电网支路顺序号，即对应的le卡中17列的顺序号le卡填写实例：,静止无功补偿,静止无功补偿器svc控制系统模型静止无功发生器statcom

9、模型静止无功补偿装置辅助信号系统模型,串补模型,固定串补模型可控串补模型,各种稳定措施模型,自动减负荷（低压、低频减载），最多可以设置10级电气制动低压自动切除电抗器电容器、电抗器自动投切控制切机快关阀门,继电保护模型,低频负荷、发电机、线路保护串联电容器间隙保护阻抗继电器省缺的距离继电器功率摇摆继电器远方继电器,第二部分,故障模型,基本的故障形式ls卡,mde1，线路出口三相短路mde2，母线三相短路mde3，线路中间三相短路mde4，切机、切负荷mde5，直流操作及两端直流功率修改*mde6，多端直流操作mde7，发电机机械功率阻尼mde8，快速关闭汽门mde9，复故障*及触发暂态稳定器，

10、冻结静补输出mde10，失磁mde11，励磁参考电压修改mde12，线路参数修改*mde13，投入线路*,复故障包含的故障子类型,nde1，单相短路nde2，两相短路nde3，三相短路nde4，两相对地短路nde5，串联电容单相击穿nde6，三相断线nde7，单相断线nde8，两相断线nde9，串联电容两相击穿nde11，异名相相间短路*nde12，异名相短路接地*,简化的故障形式flt卡,typ1，三相短路故障跳三相typ2，单相瞬时短路故障typ3，单相永久性短路故障typ4，三相短路单相拒动故障typ5，母线三相短路故障跳所有出线typ11，发电机快关typ21、22、23，切负荷故障t

11、yp31、32、33，计算三相短路电流和单相短路电流typ35、36、37，计算潜供电流typ41、42、43，计算工频过电压,复故障的修改(mde=9),不再需要增加小开关支路；线路断线故障时，原来程序故障侧无效，新程序改为有效（1-前侧；2-后侧；3-两侧）；增加了在线路任何一点发生短路故障的功能，可以在故障卡58-63列填写故障位置百分比；增加了异名相短路接地和异名相短路的故障形式（1.1.14）。,复故障实例（一）,0周波，线路n11-n12的n11侧a相发生单相短路故障lsn11230.n12230.90.01115周波，线路n11-n12的两侧a相断开，用单相断线故障模拟lsn11

12、230.n12230.95.071350周波，线路n11-n12的两侧a相重合闸重合，用单相断线故障消失模拟lsn11230.n12230.-950.071355周波，故障没有消失，线路n11-n12的两侧三相断开，用三相断线故障模拟lsn11230.n12230.955.061355周波，线路n11-n12的两侧三相断开，短路故障消失lsn11230.n12230.-955.0111,复故障实例（二）,0周波，线路n11-n12距n11侧40%处a相发生单相短路故障lsn11230.n12230.90.040.1115周波，线路n11-n12两侧跳a相lsn11230.n12230.95.0

13、71350周波，单相重合闸lsn11230.n12230.-950.071355周波，重合闸不成功，跳两侧三相lsn11230.n12230.955.0613lsn11230.n12230.-955.040.111,异名相短路故障(ls,mde=9_11/12),异名相短路故障实例,例如：线路”n11(230kv)-s11(230kv)”的a相与线路”n12(230kv)-s12(230kv)”的b相短路故障发生卡lsn11230.s11230.90.01111ls+n11230.s11230.9n12230.s12230.2故障消失lsn11230.s11230.-95.1111ls+n11

14、230.s11230.-9n12230.s12230.2,填写异名相短路故障注意的问题,必须同时填写ls卡和ls+卡；ls卡填写主要的故障信息，如故障时间、故障位置、故障形式等，同时必须在68列填写5-32列线路的故障相；ls+卡用于填写短路对应的另外一条线路的信息，在40-66列填写，以及填写40-66列所指线路的故障相，在68列填写；ls卡和ls+卡中的36-37列必须都同时填写9或-9，9表示故障发生，-9表示故障消失；对于故障消失，也必须同时填写ls卡和ls+卡，除了36-39列为-9外，其它位置的数据应与故障发生对应的ls和ls+卡相同。,修改线路参数卡(ls，mde=12),功能：

15、修改某一线路的参数填写ls卡（mde=12），其中主要填写修改的时间、新的线路参数值等；可以同时对正序、零序进行修改；注意：对应线路不能有其他故障和操作，可以进行多次修改。,修改线路参数卡实例,lslishan50525.fengy50525.125.0.00037.00379.8586.85861lslishan50525.fengy50525.125.0.00089.00654.6536.65360,投入线路数据卡(ls，mde=13),功能：投入某一线路填写ls卡（mde=13），其中主要填写投入的时间、投入线路的参数等；一般同时指定正序、零序参数，也可以只指定正序或零序；注意：对应线路

16、不能有其他故障和操作，只能投入一次。,投入线路数据卡实例,lslishan50525.fengy50525.135.0.00037.00379.8586.85861lslishan50525.fengy50525.135.0.00089.00654.6536.65360,简化的三相短路故障卡(flt,1),用于简化三相短路故障的填写；只需要填写故障位置、故障发生时间、断开时间即可；断开时间填写相对于发生时间的时间差；如果故障线路高抗采用新的填写方法，程序自动处理高抗；否则，线路三相断开时需要用ls卡切除。,简化的三相短路故障卡图例,简化的三相短路故障卡实例,采用ls卡填写的三相短路如下：ls母

17、线b230.母线1230.15.0ls-母线b230.-母线1230.-110.0采用flt卡填写的对应的形式如下：flt母线b230.母线1230.115.5.5.采用ls卡填写的线路中间三相短路为：ls母线b230.母线1230.30.040ls-母线b230.-母线1230.-35.040如果采用flt卡填写，对应的形式如下：flt母线b230.母线1230.110.5.5.40,简化的单瞬故障卡(flt,2),用于填写线路单相短路、单相跳开、重合闸重合成功的过程；需要填写故障相、故障发生时间、断开时间、重合时间；断开、重合时间填写相对于发生时间的时间差；如果故障线路高抗采用新的填写方法

18、，程序自动处理高抗；否则，线路三相断开时需要用ls卡切除。,简化的单瞬故障卡图例,简化的单瞬故障卡实例,采用ls卡填写的形式如下：lsn11230.n12230.90.0111lsn11230.n12230.95.0713lsn11230.n12230.-95.0111lsn11230.n12230.-950.0713如果采用flt卡填写，对应的形式如下：fltn11230.n12230.2110.5.5.50.50.,简化的单永故障(flt,3),模拟单相永久性短路故障，即线路某一相单相短路、故障相两侧断开、故障相两侧重合闸重合、重合闸重合失败两侧三相断开；需要填写故障相、故障位置、故障发生

19、时间、断开时间、重合时间、三相断开时间；断开时间填写相对于发生时间的时间差；如果故障线路高抗采用新的填写方法，程序自动处理高抗；否则，线路三相断开时需要用ls卡切除。,简化的单永故障图例,简化的单永故障实例,用ls卡填写lsn11230.n12230.90.0111lsn11230.n12230.95.0713lsn11230.n12230.-950.0713lsn11230.n12230.-955.0111ls-n11230.-n12230.-155.0用flt卡填写fltn11230.n12230.3110.5.5.50.50.55.55.,简化的三相短路单相拒动故障(flt,4),模拟线

20、路三相短路开关单相拒动的严重故障，即线路三相短路、然后两侧开关断开，但其中一相开关拒动、经过一定时间后后备保护动作切除线路，程序在模拟时单相拒动侧为近距离故障侧；需要填写故障相、故障位置、故障发生时间、断开时间、三相断开时间；断开时间填写相对于发生时间的时间差；如果故障线路高抗采用新的填写方法，程序自动处理高抗；否则，线路三相断开时需要用ls卡切除。,简化的三相短路单相拒动故障图例,简化的三相短路单相拒动故障实例,用ls卡填写0周波存在三相短路lsn11230.n12230.90.03115周波存在单相短路、两相断线、三相断线，没有三相短路lsn11230.n12230.-95.0311lsn

21、11230.n12230.95.0841lsn11230.n12230.95.0612lsn11230.n12230.95.011150周波存在三相断线，没有单相短路、两相断线、三相断线lsn11230.n12230.-950.0841lsn11230.n12230.-950.0612lsn11230.n12230.-950.0111ls-n11230.-n12230.-150.0用flt卡填写fltn11230.n12230.410.5.5.50.,简化的母线故障(flt,5),模拟某一母线发生短路故障、跳开与该母线相连的所有出线，短路故障可以是单相、两相、三相短路；需要填写短路故障形式、故

22、障发生时间、断开时间、三相断开时间；断开时间填写相对于发生时间的时间差；如果故障线路高抗采用新的填写方法，程序自动处理高抗；否则，线路三相断开时需要用ls卡切除。,简化的母线故障图例,简化的母线故障实例,采用ls卡填写lsn21230n19230.20.0ls-n21230-n19230.-22.0ls-n21230-n16230.-12.0ls-n21230-n22230.-12.0采用flt卡填写fltn21230230.2.,新增的切负荷故障(flt,21/22/23),按照给定的百分比或总切负荷大小切除静态负荷和马达负荷（不能用来切发电机）；可以按照节点、分区或区域的形式切负荷，在32

23、-33列分别填写21、22、23；可用于切除并联无功补偿；可以切除含有配电网的综合负荷模型。,新增的切负荷故障实例,切节点负荷（50mw，25mvar）fltbus1230.21w5.50.25.切节点负荷（50%，50%）fltbus1230.215.50.50.切分区负荷（30%，20%）fltz1225.30.20.切区域负荷（40%，30%）fltarea235.40.30.切节点马达负荷（60%）fltbus1230.1215.60.切节点无功补偿（20mvar）fltbus1230.21w5.20.,新增的切负荷故障注意事项,只能切除负荷，不能增加，切除的负荷大于实际值，全部切除，

24、不过切；节点形式，可以按照百分比或者总量（有名值）切除；分区、区域形式只能按照百分比切除；切除并联无功负荷时，必须按照有名值填写，可以增加可以切除，因此有正、负号；不能用于切发电机,计算短路电流功能(flt,31/32/33),目的：简化用稳定程序近似计算短路电流功能：计算单相或三相接地的短路电流；可以按照节点、分区、区域形式填写，分别填写31/32/33；可以指定分区、区域中计算的电压等级范围；注意：稳定程序可以近似计算短路电流，与专用的短路电流程序考虑的问题有一定的差别，结果会有所不同。,计算短路电流功能实例,fltbus1230.313fltbus1230.311fltbus1230.3

25、130.01fltz1323fltarea333fltarea333230.fltz1321230.fltz13230.01230.525.,计算短路电流功能输出实例,flt母线a230.311flt母线a230.313*共有1个节点计算三相短路电流，计算结果如下：节点名称短路电流(pu)短路电流(a)短路容量(mva)母线a230.0024.4195166129.8332441.952*共有1个节点计算单相短路电流，计算结果如下：节点名称短路电流(pu)短路电路(a)母线a230.001.60978404.089,计算短路电流功能输出实例,fltgz323220.*分区gz中有38个节点需要

26、计算三相短路电流，具体结果如下：节点名称短路电流(pu)短路电流(a)短路容量(mva)zhugang220.00129.48667933981.44112948.668zhannech220.0089.37519123454.9048937.519xintang220.00112.72731829583.25011272.731wuxianb220.0078.78638520676.0647878.639wuxiana220.00118.21281431022.81811821.281shijing220.00105.09585627580.51010509.586runzhou220.00

27、80.29019221070.7118029.019lumingb220.0094.22231324726.9459422.231luminga220.0076.93368520189.8577693.369lianzhon220.00112.07042729410.86111207.043jiahe220.00110.66001129040.72111066.001hengyund220.00124.69883032724.95712469.884hengyun220.00118.12229930999.06411812.229chashan220.0071.55772418779.0337

28、155.772bishan220.00104.86088627518.84610486.089,计算潜供电流功能(flt,35/36/37),计算线路a相接地、a相两侧断路器跳开时的短路电流，即潜供电流；可以按照线路、分区或区域的形式填写，可以指定分区、区域中的电压等级范围；程序计算时对一条线路的两侧都进行计算，并输出两侧对应的潜供电流；注意：线路高抗对潜供电流有明显的影响，如果有高抗，必须采用新的填写方法，否则，计算结果不正确。,计算潜供电流功能图例,计算潜供电流功能实例,fltn11230.s11230.35fltz136fltarea37fltz136230.525.fltarea372

29、30.525.,计算潜供电流功能输出实例,flt0136*分区01中有6个线路需要计算潜供电流，具体结果如下：支路名称前侧潜供电流(a)后侧潜供电流(a)母线1230.0母线a230.01.9251.996母线1230.0母线b230.01.8171.730母线2230.0母线a230.03.5813.080母线2230.0母线c230.01.6751.629母线3230.0母线b230.04.0663.242母线3230.0母线c230.02.2952.491,计算工频过电压功能(flt,41/42/43),计算线路两种形式的工频过电压，即：线路一侧a相接地短路、故障侧开关三相断开，计算b、

30、c相电压；线路一侧a相接地短路、故障侧a、c相开关断开，计算c相电压；可以按照线路、分区或区域的形式填写，可以指定分区、区域中的电压等级范围；程序计算时对一条线路的两侧都进行计算，并输出两侧对应的两种情况的工频过电压；注意：线路高抗对工频过电压有明显的影响，如果有高抗，必须采用新的填写方法，否则，计算结果不正确。,计算工频过电压功能图例,计算工频过电压功能实例,fltn11230.s11230.41fltz142fltarea43fltz142230.525.fltarea43230.525.,计算工频过电压功能输出实例,flt0142*分区01中有6个线路需要计算工频过电压，具体结果如下(p

31、u)：支路名称/-a相短路、abc三相断开-/a相短路、ac两相断开/前侧b相前侧c相后侧b相后侧c相前侧c相后侧c相母线1230.0母线a230.02.010192.182552.169742.304312.130402.20604母线1230.0母线b230.02.079832.506602.110122.259052.490052.15967母线2230.0母线a230.02.200082.503872.306542.580032.437792.47686母线2230.0母线c230.02.027552.169462.158132.289852.120422.18649母线3230.0母

32、线b230.02.272642.726752.384842.708132.695492.62766母线3230.0母线c230.02.074622.233092.346202.552412.161742.43962,使用故障卡时注意的问题,原理的故障卡如果填写正确，新版本程序不需要修改仍然能够继续使用；基本故障卡ls卡和简化故障卡flt卡可以联合使用；flt卡仅用于填写特定形式故障，目的是简化常用故障的填写方法，如果填写的故障不能采用flt卡填写，可以采用ls卡填写或flt、ls卡的组合方式填写；同一条支路在同一个时刻不能同时存在多个故障形式，例如，一条支路同时存在单相断线，又同时存在线路三相

33、断线；对于串补支路，如果在稳定数据中填写串补的数据卡rz卡，则不能填写任何形式的故障操作卡；部分功能需要考虑线路高抗的影响，如果采用新形式填写，程序能够自动处理高抗；否则，计算结果可能不正确。,第三部分,计算控制说明,case数据卡,指定潮流方式名，必须与对应的潮流数据一致指定存储计算数据相隔的步数，数据点比较多时用制定采用新、旧减负荷的标志填写通用阻尼绕组参数负荷负序电纳比例值,计算控制卡ff卡,计算步长dt和计算时间endt改变计算步长dtc、istp指定迭代误差toli、delang指定网络方程和微分方程的迭代次数交、直流积分步长的比例指定低电压提出运行，0.4pu30周波指定最大角度限

34、制noanglim指定马达滑差达到1的处理方式,计算控制卡ff卡,infbus，所有发电机为无穷大nopp，没有次暂态nodq，发电机全部为经典模型nosat，所有发电机无饱和nogov，无调速器ieqpc，所有发电机eq恒定noex，无励磁nosc，无pssnoload，负荷都为恒阻抗mftomg、mgtomf，发电机模型转换,计算控制卡f0卡,控制计算过程监视曲线：发电机最大功角差曲线*指定发电机功角差曲线最低电压曲线*、最高电压曲线指定节点电压曲线最低频率曲线*、最高频率曲线指定节点频率曲线,计算控制卡f1卡,计算过程控制，为ff卡的继续指定缺省发电机阻尼因子强制发电机采用通用阻尼参数指

35、定负荷低压转换为恒阻抗的电压值强制sp和sg卡之间的转换指定调速器gh/gs死区,第四部分,输出部分,输出可以输出的变量,母线：电压、频率、有功无功负荷等发电机：角度、速度偏差、励磁电压、机械功率、电磁功率、无功功率、励磁电流、pss信号等线路：有功、无功、电流等直流：点燃角、熄弧角、电压、电流、功率等串补：电抗、触发角、电流等前20个低电压、过电压、低频、高频节点系统电压最低的节点、功角差最大的发电机,输出可以输出的文件类型,文本输出文件.out辅助文件.swx作图文件.cur输出代码1-7表示输出到不同文件1/3/5/7，输出到out文件2/3/6/7，输出到cur文件3，输出到swx文件

36、,输出相关数据卡,90卡，表示开始*mh卡，输出控制卡*bh/b卡，母线输出卡bv卡，暂态最低电压数据卡gh/ghc/g卡，发电机输出卡lh/l/lc卡，线路输出卡dh/d卡，直流输出卡rh/r卡，串补输出卡py卡，振荡频率、阻尼比输出卡99卡，输出结束卡*ogm卡，最大发电机功角差输出卡subegin-suend输出总结卡,输出控制卡mh卡,控制输出的开始tstrt和结束时间maxt指定画图的时间刻度edt指定画图点相隔的时间步数nskp输出各变量最大最小值的起始结束时间指定输出每一步电压越限的电压值指定输出不良视在阻抗的角度aref控制是否显示计算曲线ifplot控制是否删除sol文件和o

37、pt文件,输出输出数据卡位置,90必须是第一张，表示开始mh卡在90卡后面mh卡后面是母线、发电机、线路、直流、串补相关的数据卡对于任何一个类型相关的数据卡，首先必须是控制卡，然后是输出卡；例如发电机，必须首先是gh卡，然后是g卡，所有的g卡放在一起；其它类似99放在后面，表示结束,输出寻找最低电压数据卡,自动寻找指定分区内指定电压等级节点中暂态电压最低和低电压时间最长的节点填写bvbv+卡与节点输出卡b卡填写在一起，放在bh卡后面输出两个数据表：暂态电压最低的节点和低电压时间最长的节点,输出寻找功角差最大发电机,自动寻找指定分区内，各计算时间功角差最大的发电机填写ogm卡填写在99卡的后面,

38、计算阻尼比、振荡频率功能,方法一：填写py卡，可以计算发电机、线路、节点等部分相关变量曲线的振荡频率和阻尼比；方法二：在发电机输出卡g卡、线路输出卡l卡的第80列填写一个非零值，程序将自动计算对应的发电机功角、线路有功功率的衰减系数、振荡频率和阻尼比。程序在稳定结果文件out文件末尾集中输出所有计算的衰减因子和振荡频率。挑选出2个主要的振荡模式输出，通常情况下，第一个振荡模式为主要的振荡模式，有时第二个是主要的振荡模式可以采用mhc卡选择计算的时间段,计算阻尼比、振荡频率功能,线路幅值1衰减系数1振荡频率1阻尼比1母线b230.0母线1230.0101.339-0.2671.0490.041母

39、线a230.0母线1230.035.851-0.2811.0530.042,输出结果总结卡,subegin开始，suend结束，放在ogm卡后面，输出部分计算数据可以输出的数据：线路、变压器功率（初值）断面功率（初值）母线电压（初值和暂态最低值）发电机暂态过程最大摆角bv卡指定的暂态电压最低的节点名、低电压时间ogm卡功角最大的发电机根据电压和阻尼比自动判断系统是否稳定优点：关心的结果集中输出，简单明了,输出结果总结卡实例,subeginsu断面潮流lbus1525.bus2525.12840.01407.12628.07080lbus1525.bus2525.22840.01407.1262

40、8.07080su暂态最低电压bbus1525.3su发电机的最大摆角ggen115.83blgmsuend,输出结果总结卡实例,断面潮流1000mw，暂态最低电压0.76pu，发电机的最大摆角156度，最低电压值是0.68，对应节点”bus2525.”，低于0.75持续时间30周波，功角差最大的发电机是”gen115.8”和”gen210.5”，最大功角差为88度。,第五部分,常见问题说明,输入数据格式说明,采用固定格式输入，填写数据必须在要求的列字符串：标识为“ai”，其中的i为一整数，表示字符串的长度整数：标识为“ii”，其中i为整数的位数，可以填写位数小于或等于i的整数。浮点数：标识为

41、“fm.n”，其中的m为输入数据占有的最大列数，n为小数点后面的位数。,输入数据格式说明浮点数,加入了小数点，则程序读入后的数值就是所填写的数，即使输入浮点数的格式不与要求的相同如果输入的数据没有小数点，则按照缺省处理，读入后自动截取n位小数，如果输入数据达不到n为，则前面补0填写时应该：应按照要求格式在指定列内规范化填写；最好填写小数点，这样即可以减少出错又可以增加可读性；如果不填写小数点，程序将缺省处理，最好满格填写。,发电机部分参数基准容量问题,发电机卡mf卡29-32列的容量为发电机部分参数的基准容量，可能需要折算的参数：发电机m卡、mf中的电阻和电抗pss模型sp卡和sg卡中的kqs

42、，sh卡中的kp，sh+卡中的kpm和kpe，si卡中的kr（pss卡中不指定情况下）励磁模型ed卡中的ki；所有f开头的励磁模型卡中的rc和xc；励磁模型fd卡中的ki；励磁模型fj卡中的kc；励磁模型fl卡中的xl和ki；双轴汽轮发电机调速器和原动机模型gc卡中的dh对于实际的发电机模型，电阻电抗参数的基本容量一般都是以发电机额定容量为基准的，如果mf卡中29-32列填写的容量值不等于发电机额定容量，必须根据基准容量对参数进行折算；如果mf卡中的29-32列的基准容量填写发电机的基准容量，则可以避免折算问题，减少出错。,小开关电抗问题,小开关的电抗通常取0.00010.0004pu，如果太

43、小可能会引起计算出错；,计算步长的选取,不同的步长对计算结果稍有影响，特别是极限情况下；发电机采用详细模型、有马达负荷等条件下，由于存在小的时间常数，如部分励磁时间常数、发电机次暂态时间常数等，计算步长应取较小；一般建议采用0.5周波或1.0周波，缺省为1周波可以先采用小步长，计算一段时间后采用较大步长计算,采用马达负荷模型注意的问题,建议采用新形式的马达模型mlmjmk填写马达数据；对于容量很小的马达，例如1mw以下的马达，对系统稳定性通常影响较小，但是容易引起数值计算问题，因此建议不采用马达负荷；步长应采用较小步长，一般采用0.5周波或1.0周波；,计算不平衡信息,加强了稳定计算前初始状态

44、的检查，可能会出现原来能够计算但是现在计算中止的现象；错误信息：初始化不平衡、计算中止；电压差别最大的节点为“节点名”，潮流电压为“电压值”，稳定计算前电压为“电压值”；该错误信息中的“节点名”或附近节点数据有错误，应检查是否有如下情况：(1)潮流中有无功出力、稳定中没有发电机卡(2)潮流中存在电抗很小的支路（0.0001）(3)潮流收敛性较差,稳定曲线输出修改,原来程序在进行稳定计算完成之后，需要按任意键直到所有输出曲线画面显示完毕后，才能形成稳定曲线文件（*.cur）。现在程序对此进行改进，不再需要按任意键，在计算完成后并显示第一个输出曲线画面时，曲线文件已经形成，如果不需要看所有的曲线，

45、则直接关闭计算窗口即可,监视曲线最大发电机功角差修改,稳定计算过程中，通常会动态显示三条监视曲线，分别为最大发电机功角差曲线、最低母线电压曲线、最低频率曲线，其中最大发电机功角曲线为全系统所有发电机的最大功角差。如果整个系统是由多个同步网组成的（它们之间可能没有联系或采用直流相连），此时显示功角差曲线时对应最大功角差的两台发电机可能不在一个同步网内，这样显示的功角差实际上没有意义，同时可能会导致误解。因此，本版本对此进行了修改，最大功角差曲线显示的是同步网中功角差最大的发电机。如果有故障操作，则仅显示故障操作所在同步网中的最大功角差；如果没有故障操作，则所有同步网中寻找各自的最大功角差，显示其中的最大值。,稳定发电机数据处理方式的变化,如果潮流中发电机的有功无功出力为0，则自动删除稳定数据中该发电机相关的数据卡，给出警告信息（原程序不删除，参与计算）；如果稳定数据中发电机及其相关数据卡的发电机在潮流数据中不存在，则自动删除该发电机相关数据卡，并给出警告信息；（原程序给出错误信息，中止计算）；由于新程序和旧程序对此处理方式不同，注意对计算结果的影响。,调速器模型gs限幅的修改,第六部分,风电系统模型,固定转速风电模型,固定转速风电模型,风速作用在风力机叶片上产生转矩，该转矩驱动轮毂转动。然而，由于风能密度低，此时风力机转动的