电力系统全过程动态仿真的组合数值积分算法研究.pdf

《电力系统全过程动态仿真的组合数值积分算法研究.pdf》由会员分享，可在线阅读，更多相关《电力系统全过程动态仿真的组合数值积分算法研究.pdf（7页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、第 29 卷 第 28 期 中 国 电 机 工 程 学 报 Vol.29 No.28 Oct.5,2009 2009 年 10 月 5 日 Proceedings of the CSEE 2009 Chin.Soc.for Elec.Eng.23 文章编号：0258-8013(2009)28-0023-07 中图分类号：TM 74 文献标志码：A 学科分类号：47040 电力系统全过程动态仿真的组合数值积分算法研究 宋新立，汤涌，刘文焯，仲悟之，吴国旸，刘涛(中国电力科学研究院，北京市 海淀区 100192)Mixed Numerical Integral Algorithm for Full

2、 Dynamic Simulation of the Power System SONG Xin-li,TANG Yong,LIU Wen-zhuo,ZHONG Wu-zhi,WU Guo-yang,LIU Tao(China Electric Power Research Institute,Haidian District,Beijing 100192,China)ABSTRACT:Full dynamic simulation of the power system is such a method in which the electro-mechanical transient,an

3、d medium-term and long-term dynamic phenomena are unified.For the large difference between the time constants of dynamic simulation models,it is a typical nonlinear stiff system simulation,and the numerical integral algorithm plays a key role in simulation.Based on the analysis of the requirements f

4、or solving nonlinear stiff systems and the main problems in the existing numerical integral algorithm for dynamic simulation of the power system,a new mixed integral algorithm is proposed.Integrating the advantages of both implicit trapezoidal integration with fixed time step and Gear variable time

5、step algorithm,the new algorithm can greatly improve the speed of transient stability in full dynamic simulation and easily deal with the discontinuous equations in control systems.The results from two simulation cases show that the proposed algorithm is effective and feasible.KEY WORDS:full dynamic

6、 simulation;electro-mechanical transient;medium-term and long-term dynamic;stiff system;numerical integral algorithm;discontinuity process 摘要：全过程动态仿真是将电力系统机电暂态和中长期动态过程有机地统一起来的一种数字仿真。仿真模型中的时间常数差别很大，是典型的刚性非线性系统。数值积分算法是全过程动态仿真的核心。分析非线性刚性系统求解的特点和电力系统动态仿真中现有数值积分算法存在的主要问题，提出一种新的适合全过程动态仿真的组合数值积分算法。该算法有机地结合

7、了固定步长隐式梯形积分法和变步长吉尔(Gear)法的优点，克服现有变步长吉尔法在机电暂态过程中计算效率低下和间断处理复杂的问题。算例仿真与分析的结果 基金项目：“十一五”国家科技支撑计划重大项目(2008BA A13B03)；国家重点基础研究发展计划项目(973 项目)(2004CB217901)。Key Project of the National Eleventh-Five Year Research Program of China(2008BAA13B03);The National Basic Research Program of China(973 Program)(2004C

8、B217901).表明了新算法的有效性和可行性。关键词：全过程动态仿真；机电暂态；中长期动态；刚性系统；数值积分算法；间断处理 0 引言 电力系统全过程动态仿真是将电力系统机电暂态、中期和长期过程有机地统一起来进行数字仿真，能够描述电力系统受到扰动之后整个连续的动态过程，是电力系统规划设计、调度运行和科学研究的重要仿真计算工具之一1-3。随着中国特高压骨干电网的快速建设和智能电网的加快实施，规模巨大的全国性交直流互联电力系统即将形成，这使得电力系统特性发生变化，系统的静态和动态行为变得更加复杂。因此，为建设中国统一坚强的智能电网，迫切需要研究和分析这种非线性超大规模电力系统动态特性机制、严重事

9、故特征及其智能稳定控制策略的全过程动态稳定仿真技术。全过程动态仿真主要涉及发电机及其励磁系统和动力系统、动态负荷、电力电子元件等众多动态元件和输电网络等组成的非线性动态系统4-5。描述这一非线性动态系统的是一组高阶的微分方程组和代数方程组。全过程动态仿真采用时域仿真法计算电力系统动态元件及其控制系统的动态过程。仿真中模型的时间常数差异很大，混合着快速和慢速动态过程，是典型的刚性非线性动态系统。时域仿真计算就是采用适当的数值积分方法求解微分代数方程组的初值问题。因此，数值积分方法是整个全过程动态仿真技术的基础和核心。电力系统全过程动态仿真用于研究电力系统在受到干扰之后系统较长时间的机电过渡过程。

10、仿真时间从几秒到数十分钟甚至若干个小时，时间跨度大。因此，需要采用自动变步长积分方法，在系24 中 国 电 机 工 程 学 报 第 29 卷 统的快变阶段(机电暂态)使用小步长计算，而在慢变阶段(中长期动态)使用大步长。现有全过程动态仿真程序的数值积分方法大多采用吉尔类变步长方法，例如瑞典的SIMPOW程序6、法国和比利时的EUROSTAG程序7-8和中国电力科学研究院开发的全过程动态仿真程序9。这种方法的优点是暂态过程及中长期动态过程可以采用统一的模型和 数值积分方法，在中长期动态过程中可以大步长进行仿真。但在应用实践中发现这种方法的主要缺点和问题为：1）机电暂态过程中计算步长过小，导致仿真

11、效率很低，例如，对于中国东北华北华中联网系统的仿真计算，全过程动态仿真程序在快变的机电暂态阶段，其计算速度要比现有暂态稳定程序慢 56 倍；2）算法难以处理模型中的间断环节。本文分析刚性系统求解的特点和现有电力系统动态仿真中数值积分方法的主要问题，提出一种新的适于电力系统机电暂态及中长期动态过程的组合数值积分算法，并通过算例说明新算法的有效性和可行性。1 刚性系统求解的特点 1.1 求解时的数值稳定性要求 对于非线性系统：d()/d(,(),0,tttttT=yfy (1)式中：为待求的 m 维函 T11()(),(),.,()mty ty tyt=y数向量；t为时间；T为仿真总时间。如果是刚

12、性系统，则雅可比(雅可比)矩阵fi/yi的特征值i(i=1,2,.,m)的实部小于 0，且实部绝对值的最大值和最小值差别很大.10。对于刚性系统的数值积分求解，为保证方法的数值稳定性，选取步长 h 很重要，需要使得每个步 长值ih都落在所采用的数值积分方法的绝对稳定区域内，即满足：,1,2,.,iihhim=。当方法的绝对稳定区域有限时，例如，常见的2 阶龙格库塔法等显式方法，步长要限制在 1/max Reii的数量级，即系统最小时间常数的数 量级。如果仿真时间T(即数值积分区间0,T)与最大时间常数同级，则积分步数N=T/h就与刚性比同级，是一个很大的数。因此，要克服这种步长限制，要求数值方

13、法的绝对稳定区域包含左半复平面，即要求方法具有A稳定性。而具有A稳定性的方法必是隐式的，且其最大阶数是 2 阶。梯形积分法是具有最小误差常数的A稳定方法10。1.2 求解时的收敛性要求 由于具有A稳定的数值积分法是隐式的，因此，每积分 1 步必须伴随 1 个隐式方程的求解过程。非线性方程的求解通常有简单迭代法和牛顿迭代法 2种。现以梯形积分法为例说明刚性系统中求解的特点11。式(1)的梯形求解公式为,1,11(,)(,)22i ni ninninnhhyyf tf t+=+yy (2)式中n为计算步数。当使用简单迭代求解时，方 程(2)的收敛性条件为 12h (3)式中为相应的雅可比矩阵按模最

14、大的特征值(电力系统全过程动态仿真中最大模值的数量级为102)。对于刚性系统来说，为满足收敛条件，h被限制在很小的范围内(电力系统全过程动态仿真中h的数量级为102)，实际上这与显式方法的稳定域对h的限制是一样的，所以，简单迭代法不能满足数值积分中使用较大步长的要求，必须使用牛顿法迭代求解。牛顿法求解时，公式(2)改写为 1,1,11()(,)(,)022ini ni ninninnhhFyyf tf t+=yyy(4)则牛顿迭代可表示为(1)()()1(),1,111()()lllli ni nininyyFF+=yy+)(5)式中：l为是迭代次数；为F(),1(lii nF y+i关于变量

15、y的雅可比矩阵，并在()1ln+y取值。实际积分计算中 雅可比矩阵的元素值经常是慢变化的，对于1步或步长不变的若干步，矩阵不必重新计算，仅在不收敛或误差较大时重新计算。该迭代方法称为拟牛顿法。迭代初值通过显式积分方法容易预测得到。由此可知，对于电力系统全过程动态仿真来说，若要使用较大的仿真步长，则不能使用简单迭代求解，而采用牛顿法联立求解是有效方法之一。2 现有数值积分算法的分析 2.1 电力系统全过程动态仿真中微分代数方程的一般形式 电力系统全过程动态仿真中需要求解的微分代数方程组可表示为 11212d(,)d(,)ttt=yf yyg yy0 (6)第 28 期 宋新立等：电力系统全过程动

16、态仿真的组合数值积分算法研究 25 式中：微分方程表示电力系统元件的动态特性，是系统的状态方程；代数方程表示电力系统元件的静态特性，主要是系统的网络方程；y1为n1个状态向量(微分变量)；y2为n2个代数向量(代数变量)。微分代数方程组的阶数为n1和n2之和。电力系统全过程动态仿真对象通常是大规模交直流电力系统，其规模可达成千上万个母线和支路、数千台发电机及其控制系统、数十条直流输电线路。因此，式(6)中的方程阶数常达数万阶以上。由于微分代数方程组(6)的代数约束方程的雅可比矩阵是非奇异的，所以电力系统全过程动态中的微分代数方程组的秩为1，是一个隐式的非线性微分方程组，能够使用刚性微分方程理论

17、和算法进行数值积分求解12。2.2 基于吉尔法的变步长积分方法 吉尔法是公认的求解刚性问题的有效方法之一，特别是对于计算精度要求不高的系统。对于变步长的吉尔积分法，其计算步骤主要包括预测、校正迭代、截断误差计算和自动变阶变步长控制4 步13。吉尔法能够自启动，起步时使用1阶。由于只有2阶及以下的吉尔法的稳定域能够覆盖复平面的左半平面，所以使用的最大阶数为2阶。校正迭代采用拟牛顿法，因而吉尔法满足刚性系统求解的数值稳定性和收敛性要求。电力系统微分代数方程组的特征值多为复数，时域解表现为振荡曲线。经典吉尔法应用到电力系统微分代数方程时，则会遇到绝对稳定域覆盖复平面的右半平面，导致在有的情况下出现计

18、算危险性问题8，即仿真出的电力系统电气量或控制量的阻尼比实际结果要强，从而使仿真出现较大偏差，得出错误的结论。虽然选择高精度或小步长可避免这种问题，但代价是对于本来不存在计算危险性的系统，仿真速度也很慢。因此，电力系统全过程动态仿真中需要对其微分代数方程组的微分和代数变量分别处理，对微分变量使用隐式亚当姆斯(Adams)法(2阶隐式亚当姆斯法就是梯形积分法)；对代数变量使用吉尔法。这种方法的实现步骤与经典的吉尔法完全相同，所不同的是校正和截断误差计算中的系数需要根据变量的性质取不同的值8。2.3 吉尔法使用中的步长控制问题 通过大量的电力系统动态仿真实践发现，对于吉尔法中迭代误差和变步长控制策

19、略，很难找到一种合理的参数，使得在不同的电力系统运行方式和不同的故障形式下，对机电暂态和中长期动态过程都具有较满意的仿真步长。如果误差控制严格，机电暂态过程步长容易过小，例如，小于1 ms。虽然对于中长期动态过程，这种严格的误差控制能够使仿真采用合理的大步长，但在暂态阶段仿真速度很慢；反之，如果误差控制条件放宽，则虽然机电暂态的初期步长控制的比较合适，但在随后的机电暂态及慢变(中长期动态)过程中，累积误差会很大，校正迭代中往往不收敛，需要强制减小步长，计算量反而增加很多；此外，从编程实现方面考虑，仿真中自动更改控制参数实现较困难。所以，现有的吉尔算法很难找到一个误差和步长控制的计算方法，以使不

20、同的仿真算例在快变和慢变过程都能使步长得到合理的控制。实践中，为保证计算可靠，常常采用误差控制严格的方法，以使得仿真精度有保证，但大大限制了仿真速度。电力系统全过程动态仿真中往往有频繁的自动装置动作，产生很多的机电暂态过程，导致在较长的仿真时间内积分计算只能使用很小的步长。2.4 吉尔法使用中的间断问题分析 采用吉尔法时，当某个变量积分区间包含间断点时，在预测阶段不能预知右端函数会发生间断，即由一个函数变换为另一个函数，导致在校正阶段或截断误差计算中该变量的误差很大，因此，算法需要将步长缩小，同时雅可比矩阵也要重新计算，直到积分区间不包含间断为止。然后因误差减小，步长放大，又将间断点包含在积分

21、步中间，因此又将缩小步长。这样反复多次，最后以小步长“爬”过间断点。对于小规模刚性系统的仿真，例如航天领域等，间断问题已有不少解决方法13。对于电力系统全过程动态仿真来说，控制系统中存在大量的限幅、死区等非线性环节14-15，例如，一台发电机的励磁系统和调速系统模型中，间断环节数量常超过10个，而对于一个具有千台以上发电机模型的大规模电力系统来说，这种环节的数量非常大。小步长引起的计算步数增加及反复的雅可比矩阵修改和重新分解都使计算量增加，计算效率低下。目前电力系统动态仿真中有2种解决这一问题的措施：1）先假定函数是连续的，正常积分，该步完成后，如果发现跨过间断点，则需要计算出断点出现的准确时

22、间，并抛弃刚刚完成的积分步，然后把步长调整到断点处，以顺利通过断点，之后换为另一个函数。该方法在SIMPOW程序中被采用，其优点是间断计算准确，但缺点是对大规模电力系统26 中 国 电 机 工 程 学 报 第 29 卷 来说，控制系统间断点数量非常多，反复回退到断点处，使得积分效率很低，而且编程实现时判断处理很复杂。2）迭代中遇到间断时，把发生间断的控制系统模块独立出来，单独实行小步长积分计算完成该积分区间步，同时修改该模块求解变量的截断误差计算公式，以使该模块变量不影响整个系统的误差判断。这种方法的主要问题在于单独模块的小步长积分方法及间断点过后的处理。如果采用不变阶的固定步长计算，则虽然编

23、程实现简单，但是算法中变量的2阶导数不易计算出；如果采用变阶变步长算法，很难保证每次小步长后计算方法都是2阶。变量的2阶导数不能计算出来将导致该模块的后续积分计算误差增大，使得整个系统积分步长减小，因而存在与不处理间断导致步长小的类似问题。另外，该方法需要独立的模块积分计算，也存在编程实现复杂的问题。间断问题大多数出现在系统快变阶段，此时，变量变化剧烈，含限幅等非线性环节的微分方程的右端函数经常切换。在电力系统控制系统模型中，因含间断的环节数量众多，采用上述2种措施进行编程开发的难度很大，程序的可靠性难以保证。2.5 固定步长的隐式梯形积分法 现有的电力系统暂态稳定仿真程序中固定步长的隐式梯形

24、积分法已经比较成熟，得到了广泛的应用，例如国内电力系统常用的PSD-BPA电力系统分析程序和PSASP电力系统综合分析程序16-17。这种算法的优点是：1）具有A稳定性，能够自启动计算；2）机电暂态过程中积分步长通常为较小(0.01 0.02 s)，能够满足电力系统机电暂态过程仿真要求(最大仿真时间一般在2030 s以内)；3）采用简单迭代法求解，控制系统的间断环节易处理。所以，该方法具有编程简单可靠、便于扩展的优点。算法的具体描述见文献18。但是这种适于机电暂态过程的积分方法不能直接推广应用到中长期动态仿真程序。原因是该方法采用简单迭代求解，由1.2节的迭代收敛性要求知，仿真步长较小，而中长

25、期动态的慢速过程时间很长，长达数分钟甚至数十分钟，如果仿真步长一直很小，则仿真步数达数万步以上，仿真效率很低。3 组合数值积分算法的构造 3.1 基本思路 固定步长的梯形积分法和变步长的数值积分法都具有A稳定性，且都是自启动的，计算中2种方法可以相互切换。由此，构造新算法的基本思路在于在仿真中可使2种积分方法得到有机结合，扬长避短，根据电力系统动态过程的特点自动选择合适的积分方法：1）在电力系统全过程仿真的机电暂态过程中采用固定步长的梯形积分法，动态元件的微分方程和电力网络的代数方程进行简单迭代求解；2）在中长期动态过程中采用变步长的吉尔法，微分方程和代数方程联立求解；3）固定步长和变步长2种

26、方法在仿真中依据一定的策略自动切换，从而在保证数值稳定性和仿真精度的前提下，大大缩短仿真时间，提高了程序的计算效率。3.2 组合算法的切换策略 积分方法切换策略是算法的关键点，其依据机电暂态过程和中长期动态过程的如下不同特点：1）机电暂态过程属于系统的快变阶段，电网母线电压较低(一般在0.8 pu以下)，求解变量值变化剧烈。此时，如果采用吉尔法，则步长会较小，一般小于0.01 s，甚至常常低于0.001 s，求解缓慢；如果采用固定步长的梯形积分法，步长可固定在0.010.02 s，且采用简单迭代求解，不需要联立求解大型方程组，所以既能保证精度又有相对较快的速度。2）中长期动态过程属于系统的慢变

27、阶段，正常情况下电网最低母线电压较高(一般在0.8 pu以上)，求解变量值变化缓慢。此时，如果采用吉尔法，则步长会较大，一般大于0.05 s，甚至达数秒以上；如果采用固定步长的梯形积分法，则迭代次数较少，一般为23次(固定步长梯形积分的迭代次数少也是系统处于慢速动态过程的一个明显特征19)。由以上2特点知，根据母线电压、步长和固定步长积分法的迭代次数可以判断系统处于机电暂态还是中长期动态过程。此外，为保证能可靠地判断系统的动态过程,一种算法要持续计算一定步数。由上分析可构造组合数值积分算法的切换策略：1）当满足如下条件之一时，切换到固定步长积分法。系统发生故障或操作之后导致网络结构改变或进行了

28、切机切负荷等动态元件的投入或切除，例如输电线路发生短路故障而跳开、发电机切机等；在采用变步长的数值积分法时，步长小于0.01 s，且持续30步。2）当满足如下条件之一时，切换到变步长积第 28 期 宋新立等：电力系统全过程动态仿真的组合数值积分算法研究 27 分法。采用固定步长的积分法的迭代次数只有2次，且持续40步；电网最低母线电压大于0.85 pu，迭代次数只有2次，且持续20步。3.3 组合算法的程序步骤 基于新的组合数值积分算法的电力系统全过程动态仿真程序框图如图1所示。本时步计算成功？有故障或操作？求取 t+时刻积分初值 积分方法选择 结束 时间到？t=t+h,n=n+1结束 截断误

29、 差小于 容许值？变阶 变步长 判断 收敛？校正计算 预测，k2=0 截断误 差计算 k2=k2+1最大迭代次数？判断是否已到最小步长？变阶 变步长N Y YYNN N N YYYN k1=0 最大迭代 次数？判断 收敛？解网络 方程 求注入 电流 计算微 分方程 k1=k1+1 N Y Y N 积分方法？开始积分，t=0，n=0 初始化 固定步长的 梯形积分法 变步长吉尔法 开始 N 图 1 基于新算法的全过程动态仿真程序流程图 Fig.1 Flow chart of the full dynamic simulation program of the power system based

30、on the new method 4 算例验证 为验证本文所提出的新数值积分算法在机电暂态及中长期动态过程中的有效性，文中列举了2个大型电力系统仿真算例，使用吉尔法和组合数值积分算法，比较了2种算法的仿真过程：1）机电暂态稳定的仿真步长与仿真效率。该算例取自中国南方交直流互联电网2006年夏大运行方式，含2 843个交流母线、1 489条交流输电线路、3回直流输电线路、343台发电机和676个负荷点。求解的微分代数方程组阶数为15 224。仿真的故障形式为受端广东电网内的博罗横沥500 kV交流线路三相永久性故障，这是该方式下电网中比较严重的三相永久性故障之一。仿真时间为3 s，0.2 s时

31、该线路发生故障，0.3 s线路两侧开关跳开。该算例主要验证组合数值积分算法仿真机电暂态过程的有效性。组合数值积分算法和吉尔法的计算结果完全相同，此处仅以图2所示的故障处节点电压U的变化曲线为例进行比较，但2种算法的仿真步长h的变化完全不同，如图3所示。1.00.50.0012 3t/s U/pu 新数值积分算法吉尔法 图 2 母线电压的仿真曲线 Fig.2 Curves of the bus voltage simulated 0.080.040.00012 3t/s h/s 新数值积分算法 吉尔法 图 3 仿真步长的比较 Fig.3 Comparison of the time step i

32、n simulation 采用吉尔法的仿真过程为：故障期间及故障切除后的0.4 s内，吉尔算法为保证迭代收敛性采用较小的步长计算了209步，平均仿真步长在0.001 5 s左右。0.6 s之后的2.4 s内吉尔法增加步长，仿真了191步，平均步长为0.012 s。增加步长过程中，由于控制系统间断环节的影响，导致仿真步长反复减小。例如，受端系统2台发电机的励磁电压UF变化曲线如图4所示，励磁电压多次达到了限幅值，即对应的微分方程右端项出现了多次间断，相应的仿真步长出现大幅下降。整个仿真研究的3 s内，仿真了424步，耗费机时16.2 s。012 3t/s UF/pu 发电机(NPS)发电机(CP

33、K)6402 图 4 发电机励磁电压的变化 Fig.4 Changes of the excitation voltage in generator models 采用组合数值积分算法的仿真过程为：相应的故障期间及故障切除后的0.4 s内，采用固定步长的梯形积分法步长为0.01 s，只计算了40步。整个28 中 国 电 机 工 程 学 报 第 29 卷 仿真研究的3 s内，仿真了291步，耗费机时5.2 s。可见，对故障严重情况下的机电暂态过程，采用组合数值积分算法的仿真效率要明显高于吉尔积分算法。2）机电暂态及中长期动态的仿真。该算例取自中国西北电网2006年夏大运行方式，含889个交流母线

34、、387条交流输电线路、179台发电机和255个负荷点。仿真中陕西电网火电机组考虑了锅炉、汽轮机以及协调控制系统模型。研究时间为150 s，采取的故障形式为：1.0 s时陕西电网韩城二厂1号机组退出运行，系统失去400 MW的发电容量；31.0 s时西电东送断面上的330 kV桃西线发生单相接地故障；31.1 s时该单相故障转为两相对地永久性短路；31.2 s时线路两侧开关跳开故障相；31.3 s时线路两侧开关重合闸；31.4 s时线路重合闸不成功跳开三相；31.5 s时安稳装置跳开断面上的其他3回330 kV输电线路，使得陕西电网与青海、甘肃、宁夏电网解列。解列后的陕西电网功率缺额近1 20

35、0 MW。低周减载切除的总负荷为924 MW，系统最低频率约为48.85 Hz。仿真步长的变化如图5所示；秦岭火电机组的机械功率P变化如图6所示，可看出受锅炉主蒸汽压力变化的影响，中长期过程中机组出力在100 s左右存在一个低谷。0 50 100 150t/s h/s 新数值积分算法 Gear 法 0.4 0.0 0.2 图 5 仿真步长的变化 Fig.5 Changes of the time step in simulation 0 50 100 150t/s P/MW 新数值积分算法 Gear 法 120 100 110 图 6 发电机机械功率的变化 Fig.6 Mechanical p

36、ower of generators 仿真中系统在1 s时发生切机，由于故障较轻，2种算法的仿真步长都增加很快，最大值在0.1 s左右。3134.5 s以机电暂态过程为主，经历了线路复故障、跳开线路和低周减载等引起的机电暂态过程，采用吉尔法积分的仿真步长在此期间较小，有173个积分步长小于0.001 s；采用新数值积分算法的仿真在此阶段采用固定步长的梯形积分法，步长为0.01 s。35 s后系统进入中长期动态过程，2种算法的仿真都能自动采用较大步长(最大为0.1 s左右)。该算例说明组合数值积分算法对机电暂态及中长期动态仿真具有较强的适应性，在中长期动态仿真中能自动采用较大步长仿真，有效地提高

37、了仿真效率。5 结论 本文提出了一种能适于电力系统机电暂态及中长期动态(全过程)仿真的组合数值积分算法。该算法有机地结合了固定步长的隐式梯形积分法和变步长的吉尔法的优点，并通过一定的切换策略，使二者在仿真中自动切换：在电力系统的机电暂态过程中自动采用隐式梯形积分算法；而在中长期动态过程中自动采用变步长吉尔算法。采用组合数值积分算法的电力系统全过程动态仿真算例说明，组合算法能有效地解决现有的变步长吉尔法在电力系统机电暂态阶段存在的计算速度过慢和间断环节处理复杂的问题，从而大大提高电力系统全过程动态仿真程序的仿真效率和实用性。参考文献 1 宋新立，汤涌，卜广全，等面向大电网安全分析的电力系统全过程

38、动态仿真技术J电网技术，2008，32(22)：23-28 Song Xinli，Tang Yong，Bu Guangquan，et alFull dynamic simulation for the stability analysis of large power system J Power System Technology，2008，32(22)：23-28(in Chinese)2 舒印彪，张文亮，周孝信，等特高压同步电网安全性评估J中国电机工程学报，2007，27(34)：1-6 Shu Yinbiao，Zhang Wenliang，Zhou Xiaoxin，et alSecur

39、ity evaluation of UHV synchronized power gridJProceedings of the CSEE，2007，27(34)：1-6(in Chinese)3 赵良，郭强，覃琴，等特高压同步电网稳定特性分析J中国电机工程学报，2008，28(34)：47-51 Zhao Liang，Guo Qiang，Qin Qin，et alAnalysis on stability characteristics of UHV synchronized power gridJProceedings of the CSEE，2008，28(34)：47-51(in Ch

40、inese)4 杨志新，李乃湖，陈珩华东电力系统调度员培训仿具器的动态全过程准实时仿真J中国电机工程学报，1996，18(3)：205-210 Yang Zhixin，Li Naihu，Chen HengQuasi real-time dynamic silmulation for DTS of east China power systemJProceedings of the CSEE，1996，18(3)：205-210(in Chinese)5 汤涌，宋新立，刘文焯，等电力系统全过程动态仿真中的长过程动态模型：电力系统全过程动态仿真软件开发之三J电网技术，2002，26(11)：20-

41、25 Tang Yong，Song Xinli，Liu Wenzhuo，et alPower system full dynamic simulation Part III：long term dynamic modelsJ Power System Technology，2002，26(11)：20-25(in Chinese)6 ABB CorpSIMPOW dynmic simulation user manualRVsters,Sweden：ABB Corp.，2003 第 28 期 宋新立等：电力系统全过程动态仿真的组合数值积分算法研究 29 7 Stubbe M，Bihain A，

42、Deuse JSTAG：a new unified software program for the study of the dynamic behaviour of electrical power systemsJIEEE Transactions on Power Systems，1989，4(1)：129-138 8 Astic J Y，Bihain A，Jerosolimski M The mixed Adams-BDF variable step size algorithm to simulate transient and long term phenomena in pow

43、er systemsJ IEEE Transactions on Power Systems，1994，9(2)：929-935 9 汤涌，宋新立，刘文焯，等电力系统全过程动态仿真的数值方法：电力系统全过程动态仿真软件开发之一J电网技术，2002，26(9)：7-12 Tang Yong，Song Xinli，Liu Wenzhuo，et alPower system full dynamic simulation Part I：numerical methodJPower System Technology，2002，26(9)：7-12(in Chinese)10 徐绪海，朱方生刚性微分方

44、程的数值方法M武汉：武汉大学出版社，1997：36-47 11 刘德贵，费景高，韩天敏刚性大系统数字仿真方法M郑州：河南科学技术出版社，1996：410-411 12 宋晓秋微分代数方程的数值仿真算法J计算机工程与设计，2000，21(5)：58-60 Song Xiaoqiu Numerical algorithm for diferential-algebraic equations JComputer Engineering and Design，2000，21(5)：58-60(in Chinese)13 刘德贵，费景高动力学系统数字仿真算法M北京：科学出版社，2000：246-253

45、，335-384 14 邵伟，徐政 励磁系统参数整定研究J 中国电机工程学报，2002，22(10)：12-17 Shao Wei，Xu Zheng Study on parameter setting of excitation systems JProceedings of the CSEE，2002，22(10)：12-17(in Chinese)15 王官宏，陶向宇，李文锋，等原动机调节系统对电力系统动态稳定的影响J中国电机工程学报，2008，28(34)：80-86 Wang Guanhong，Tao Xiangyu，Li Wenfeng，et al Influence of tur

46、bine governor on power system dynamic stabilityJProceedings of the CSEE，2008，28(34)：80-86(in Chinese)16 汤涌，卜广全，印永华，等PSD-BPA 暂态稳定程序用户手册R北京：中国电力科学研究院，2007 17 中国电力科学研究院综合程序课题组电力系统分析综合程序用户手册R6.2 版北京：中国电力科学研究院，2007 18 汤涌电力系统稳定计算隐式积分交替求解J电网技术，1997，21(2)：1-3 Tang YongAn implicit integration alternating sol

47、ution method for power system dynamic stability simulationJ Power System Technology，1997，21(2)：1-3(in Chinese)19 Arrillaga J，Waston N RComputer modeling of electrical power systemsM2nd edNew York：John Wiley&Sons Ltd.，2001：256-257 收稿日期：2009-06-25。作者简介：宋新立(1971)，男，硕士，高级工程师，主要从事电力系统仿真软件开发和电力系统分析技术研究，；汤涌(1959)，男，博士，教授级高级工程师，博士生导师，主要从事电力系统仿真与分析研究；刘文焯(1972)，男，硕士，主要从事电力系统分析软件开发和电力系统分析技术研究；仲悟之(1979)，男，博士，工程师，主要从事电力系统分析软件开发和分析技术研究；吴国旸(1974)，男，硕士，高级工程师，主要从事电力系统分析软件开发和分析技术研究；刘涛(1984)，男，硕士研究生，研究方向为电力系统仿真软件开发。宋新立 (责任编辑 谷子)