ATPEMTP的应用讲解PPT

1、武汉大学电气工程学院武汉大学电气工程学院武汉大学电气工程学院武汉大学电气工程学院1. 1. 电力系统暂态计算的特点电力系统暂态计算的特点武汉大学电气工程学院武汉大学电气工程学院3电力系统的暂态计算的特点电力系统的暂态计算的特点 电力系统从一个稳定状态到另一个稳定状态的过渡过电力系统从一个稳定状态到另一个稳定状态的过渡过程称作电力系统的暂态。通常负荷的波动、发电机出力程称作电力系统的暂态。通常负荷的波动、发电机出力的波动之类的系统的微小变动不作为暂态考虑。的波动之类的系统的微小变动不作为暂态考虑。 引起暂态的原因主要有：引起暂态的原因主要有：（1）雷击电力设备。）雷击电力设备。（2）短路等电气事

2、故。）短路等电气事故。（3）断线等机械事故。）断线等机械事故。（4）断路器的投切操作，如线路充电、变压器充电、）断路器的投切操作，如线路充电、变压器充电、 重合闸、甩负荷等。重合闸、甩负荷等。（5）隔离开关的投切操作，如母线投切等。）隔离开关的投切操作，如母线投切等。暂态的定义暂态的定义武汉大学电气工程学院武汉大学电气工程学院4电力系统数字仿真的分类电力系统数字仿真的分类电力系统数字仿真电力系统数字仿真稳态计算稳态计算暂态计算暂态计算潮流、故障及潮流、故障及系统系统稳定和稳定和电压稳定的静稳定计算电压稳定的静稳定计算等等有效值有效值范畴范畴系统稳定和电压稳系统稳定和电压稳定的动稳定计算定的动稳

3、定计算有效值有效值范畴范畴瞬时瞬时值值范畴范畴过电压过电压、ssr、trv、ssr、trv、变压器涌流、高次谐变压器涌流、高次谐波波及保护控制及保护控制计算计算等等电力系统的暂态计算的特点电力系统的暂态计算的特点武汉大学电气工程学院武汉大学电气工程学院5电力系统各种现象的变化速度和计算范围电力系统各种现象的变化速度和计算范围电力系统的暂态计算的特点电力系统的暂态计算的特点武汉大学电气工程学院武汉大学电气工程学院61. 制定试验电压标准。制定试验电压标准。2. 决定设备的绝缘水平。决定设备的绝缘水平。3. 调查事故原因，寻找对策。调查事故原因，寻找对策。4. 计算过电压发生概率，预测事故率。计算

4、过电压发生概率，预测事故率。5. 检查设备的动作责能，如断路器的暂检查设备的动作责能，如断路器的暂 态恢复电压和零点偏移。态恢复电压和零点偏移。6. 继电保护和控制设备的整定等等。继电保护和控制设备的整定等等。 暂态仿真的目的暂态仿真的目的电力系统的暂态计算的特点电力系统的暂态计算的特点武汉大学电气工程学院武汉大学电气工程学院71. 频率范围广，从直流到数十兆赫芝。频率范围广，从直流到数十兆赫芝。 下图是下图是cigre对各种过电压的频率分类。对各种过电压的频率分类。 暂态仿真的特点暂态仿真的特点电力系统的暂态计算的特点电力系统的暂态计算的特点武汉大学电气工程学院武汉大学电气工程学院82. 元

5、件模型因计算目的而异元件模型因计算目的而异以单相变压器为例以单相变压器为例 暂态仿真的特点暂态仿真的特点(a)工频过电压计算用工频过电压计算用电力系统的暂态计算的特点电力系统的暂态计算的特点武汉大学电气工程学院武汉大学电气工程学院92. 元件模型因计算目的而异元件模型因计算目的而异(以单相变压器为例以单相变压器为例)暂态仿真的特点暂态仿真的特点(b)操作过电压计算用操作过电压计算用(c)雷过电压计算用雷过电压计算用电力系统的暂态计算的特点电力系统的暂态计算的特点武汉大学电气工程学院武汉大学电气工程学院103. 行波现象和分布参数行波现象和分布参数 设幅值为设幅值为i的正弦波电流以波速的正弦波电

6、流以波速300m/s传播，传播，ab两点间距离两点间距离为为300m。 两点电流的最大差值：两点电流的最大差值： 50hz ：3.1410-4i 100khz ：0.628i 在高频领域，变电站的母线也在高频领域，变电站的母线也要考虑波的传播过程，并当成分要考虑波的传播过程，并当成分布参数线路处理。布参数线路处理。 暂态仿真的特点暂态仿真的特点电力系统的暂态计算的特点电力系统的暂态计算的特点武汉大学电气工程学院武汉大学电气工程学院114. 非线性元件和开关非线性元件和开关 在暂态计算中，需要考虑在暂态计算中，需要考虑元件的非线性和开关的动元件的非线性和开关的动作。最常见的非线性特作。最常见的非

7、线性特性，有变压器的励磁特性、性，有变压器的励磁特性、避雷器的电压避雷器的电压-电流特性等。电流特性等。开关也可以看成一种非线开关也可以看成一种非线性元件。性元件。 右图是避雷器的电压右图是避雷器的电压-电电流流特性。特性。暂态仿真的特点暂态仿真的特点10ka10ka20ka20ka1000kv避雷器v-i特性1000kv避雷器v-i特性mama1620kv1620kv电力系统的暂态计算的特点电力系统的暂态计算的特点武汉大学电气工程学院武汉大学电气工程学院125. 元件参数的频率元件参数的频率特性特性 在暂态计算中，需在暂态计算中，需要考虑元件参数的要考虑元件参数的频率特性。最重要频率特性。最

8、重要的是线路参数的频的是线路参数的频率特性。另外，变率特性。另外，变压器参数也有频率压器参数也有频率特性。特性。暂态仿真的特点暂态仿真的特点电力系统的暂态计算的特点电力系统的暂态计算的特点武汉大学电气工程学院武汉大学电气工程学院136. 时间跨度的要求时间跨度的要求 稳态计算的对象是一个时间断面，而暂态计算要稳态计算的对象是一个时间断面，而暂态计算要模拟一个时间过程。数字计算机不可能连续地模拟模拟一个时间过程。数字计算机不可能连续地模拟暂态现象，只能在离散的时间点暂态现象，只能在离散的时间点(步长步长t)求解，这求解，这将会导致累积误差。避免这类误差的累积是暂态计将会导致累积误差。避免这类误差

9、的累积是暂态计算程序的重要课题之一。算程序的重要课题之一。 鉴于暂态模拟的上述特点，暂态模拟比稳态计算鉴于暂态模拟的上述特点，暂态模拟比稳态计算不论是程序编制还是应用，难度都要大得多。不论是程序编制还是应用，难度都要大得多。暂态仿真的特点暂态仿真的特点电力系统的暂态计算的特点电力系统的暂态计算的特点武汉大学电气工程学院武汉大学电气工程学院14（1）明确计算目的，确定主频率。）明确计算目的，确定主频率。（2）确定作为计算对象的系统范围。）确定作为计算对象的系统范围。（3）所使用的模型是否合适（模型的前提条件和）所使用的模型是否合适（模型的前提条件和相关理论）？相关理论）？（4）选择的参数是否合适

10、？分析各种参数值对计）选择的参数是否合适？分析各种参数值对计算结果的影响程度。算结果的影响程度。（5）计算结果是否符合定性的（从物理意义判断）计算结果是否符合定性的（从物理意义判断的）或者近似计算的预测结果。如果是事故分的）或者近似计算的预测结果。如果是事故分析，是否符合事故记录。析，是否符合事故记录。暂态仿真的注意事项暂态仿真的注意事项电力系统的暂态计算的特点电力系统的暂态计算的特点武汉大学电气工程学院武汉大学电气工程学院2. 2. emtp的介绍的介绍 武汉大学电气工程学院武汉大学电气工程学院16l全称：全称： electro-magnetic transients programl世界上

11、应用最为广泛的电磁暂态计算标准程序nemtp的介绍的介绍武汉大学电气工程学院武汉大学电气工程学院17（1）1962年年emtp的创始人的创始人dommel在他的博士论文在他的博士论文中首次将水力学中流体计算的中首次将水力学中流体计算的schnyder-bergeron法用于电压和电流的行波计算。法用于电压和电流的行波计算。（2）1966年年dommel就职美国邦纳维尔电力局就职美国邦纳维尔电力局,1968年完成年完成emtp的原型的原型(约约4000行行)。（3）1973年年scott meyer接替接替emtp的开发和管理工的开发和管理工作。作。（4）现有）现有bpa-emtpbpa-emt

12、p、atp-emtpatp-emtp和和dcg-emtpdcg-emtp三个版本。三个版本。latp-emtpatp-emtp：无偿使用：无偿使用, ,在全世界拥有最多的用户。最在全世界拥有最多的用户。最新版本为第新版本为第2.12.1版，源程序约版，源程序约21.821.8万行。图形接口为万行。图形接口为atpdrawatpdraw，最新版本为第，最新版本为第5.45.4版。版。ldcg-emtpdcg-emtp：有偿使用，最新版本是：有偿使用，最新版本是20052005年发行的年发行的rvrv版。图形接口为版。图形接口为emtpworksemtpworks。emtp的开发经过的开发经过ne

13、mtp的介绍的介绍武汉大学电气工程学院武汉大学电气工程学院18（1）将潮流计算使用的节点导纳矩阵和稀疏矩）将潮流计算使用的节点导纳矩阵和稀疏矩阵处理方法用于暂态计算程序，通过求解逆矩阵处理方法用于暂态计算程序，通过求解逆矩阵来求各节点电压。阵来求各节点电压。（2）在稳态计算时，）在稳态计算时，emtp把非线形元件当作把非线形元件当作线形元件，把具有分布参数的线路当作具有集线形元件，把具有分布参数的线路当作具有集中参数的中参数的型电路，通过求解下式所示的复数型电路，通过求解下式所示的复数节点导纳方程式来求得各节点的复数电压。节点导纳方程式来求得各节点的复数电压。emtp的解法的解法 ivynem

14、tp的介绍的介绍武汉大学电气工程学院武汉大学电气工程学院19（3）在暂态计算时，）在暂态计算时，emtp采用梯形积分将微采用梯形积分将微分方程式离散化，将所有的元件用等值电导和分方程式离散化，将所有的元件用等值电导和等值电流源表示。通过求解下式所示的节点方等值电流源表示。通过求解下式所示的节点方程式来求得各节点电压的瞬时值。程式来求得各节点电压的瞬时值。emtp的解法的解法ivg（4）将）将schnyder-bergeron法用于电压和电流法用于电压和电流的行波计算，成功地解决了分布参数线路的电的行波计算，成功地解决了分布参数线路的电路计算问题。路计算问题。 nemtp的介绍的介绍武汉大学电气

15、工程学院武汉大学电气工程学院20emtp的解法的解法用下图说明建立节点电导矩阵的方法用下图说明建立节点电导矩阵的方法nemtp的介绍的介绍武汉大学电气工程学院武汉大学电气工程学院21l电阻电阻emtp的解法的解法)()(1)(2112tvtvrtil电感电感dttdiltv)()(1313)()()(2)(133113tthisttvtvltti)()(2)()(311313ttvttvltttitthistnemtp的介绍的介绍武汉大学电气工程学院武汉大学电气工程学院22l电容电容emtp的解法的解法dttdvcti)()(1414)()()(2)(144114tthisttvtvtcti)

16、()(2)()(411414ttvttvtcttitthistnemtp的介绍的介绍武汉大学电气工程学院武汉大学电气工程学院23l无损分布参数线路无损分布参数线路emtp的解法的解法)()()()(ctxzfctxzfvctxfctxfi)(2ctxzfziv)()()()(515151tzitvtzitv)()(1)(15115thisttvzti)()(1)(51515titvztthistnemtp的介绍的介绍武汉大学电气工程学院武汉大学电气工程学院24对节点对节点1emtp的解法的解法)()()()()(115141312tititititi)()()()()(2)(2)(1)(122

17、115141314321thisttthisttthisttitvtctvlttvrtvztcltr对对n个节点个节点 histtitvg)()(nemtp的介绍的介绍武汉大学电气工程学院武汉大学电气工程学院25(1)单相电阻、电感和电容串联电路，单相或多单相电阻、电感和电容串联电路，单相或多相相型电路型电路 这些是这些是emtp最基本的元件，可以构成发电机、最基本的元件，可以构成发电机、变压器、线路、负荷或电力系统其他元件，也可变压器、线路、负荷或电力系统其他元件，也可以单独使用，如用作断路器的合闸电阻、接地电以单独使用，如用作断路器的合闸电阻、接地电阻、并联电抗器、变压器的杂散电容等等。阻

18、、并联电抗器、变压器的杂散电容等等。 emtp的元件的元件nemtp的介绍的介绍武汉大学电气工程学院武汉大学电气工程学院26(2)架空线路和电缆架空线路和电缆 在在emtp中可以用常规中可以用常规型电路、无损分布参型电路、无损分布参数电路、无畸变分布参数电路、带集中电阻的分数电路、无畸变分布参数电路、带集中电阻的分布参数电路或者具有频率相关参数的分布参数电布参数电路或者具有频率相关参数的分布参数电路来模拟架空线路和电缆。使用最多的是常规路来模拟架空线路和电缆。使用最多的是常规型电路、无损分布参数电路、带集中电阻的分布型电路、无损分布参数电路、带集中电阻的分布参数电路和考虑频率相关的参数电路和考

19、虑频率相关的j.marti模型。模型。emtp的元件的元件nemtp的介绍的介绍武汉大学电气工程学院武汉大学电气工程学院27(3)变压器变压器 在在emtp中可以用单相饱和变压器模型、三中可以用单相饱和变压器模型、三相相饱和变压器模型或者多相耦合饱和变压器模型或者多相耦合r-l电路来模拟各电路来模拟各种变压器。种变压器。 通常用通常用3台单相饱和变压器模型来模拟三相变台单相饱和变压器模型来模拟三相变压器。如要模拟磁滞特性，需要在模型外并联压器。如要模拟磁滞特性，需要在模型外并联磁滞特性支路。在高频领域要考虑杂散电容的磁滞特性支路。在高频领域要考虑杂散电容的影响。影响。 emtp的元件的元件ne

20、mtp的介绍的介绍武汉大学电气工程学院武汉大学电气工程学院28(4)非线形元件非线形元件 emtp准备了各种非线形元件模型，用来模拟准备了各种非线形元件模型，用来模拟变压器和电抗器的非线形电感、避雷器的非线形变压器和电抗器的非线形电感、避雷器的非线形电阻等。电阻等。emtp中的非线形元件模型分为二大类，中的非线形元件模型分为二大类，一类是准非线形一类是准非线形(pseudo-nonlinear)元件，一类元件，一类是正规的非线形元件。对于准非线形元件，用分是正规的非线形元件。对于准非线形元件，用分段线形模型来表示非线形特性。对于正规的非线段线形模型来表示非线形特性。对于正规的非线形元件，采用补

21、偿法求解。形元件，采用补偿法求解。 emtp的元件的元件nemtp的介绍的介绍武汉大学电气工程学院武汉大学电气工程学院29(5)(5)简单电压源和电流源简单电压源和电流源 emtp将经常遇到的电源函数作为内部函数，将经常遇到的电源函数作为内部函数，其中包括阶跃函数其中包括阶跃函数(type-11)、斜角函数、斜角函数(type-12,13)、正弦函数、正弦函数(type-14)、冲击函数、冲击函数(type-15)、简化换流阀、简化换流阀(type-16)、tacs控制电源控制电源(type-17)和和tacs电源电源(type-60)。还有可接于。还有可接于两个节点间的电源两个节点间的电源(

22、type-18) 。 emtp的元件的元件nemtp的介绍的介绍武汉大学电气工程学院武汉大学电气工程学院30(6)同步电机同步电机 在在atp-emtp中有两种同步电机模型，即中有两种同步电机模型，即基于基于dq0坐标的坐标的type-59模型和基于相坐标的模型和基于相坐标的type-58模型。後者避免了模型。後者避免了dq0坐标转换成相坐标转换成相坐标引起的数值不稳定，但由于每个步长都需坐标引起的数值不稳定，但由于每个步长都需要更新矩阵并重新三角形化，因此计算时间要更新矩阵并重新三角形化，因此计算时间比前者长。两者的轴系模拟相同。比前者长。两者的轴系模拟相同。emtp的元件的元件nemtp的

23、介绍的介绍武汉大学电气工程学院武汉大学电气工程学院31(7)通用电机通用电机 emtp的通用电机模型可以模拟各种同步电的通用电机模型可以模拟各种同步电机、各种异步电机和各种直流电机，功能十分机、各种异步电机和各种直流电机，功能十分广泛，但使用不太方便。另外，通用电机模型广泛，但使用不太方便。另外，通用电机模型也是使用也是使用dq0坐标的，存在坐标转换引起数值坐标的，存在坐标转换引起数值不稳定的可能性。不稳定的可能性。 笔者已开发了基于相坐标的异步电机模型。笔者已开发了基于相坐标的异步电机模型。 emtp的元件的元件nemtp的介绍的介绍武汉大学电气工程学院武汉大学电气工程学院32(8)(8)开

24、关开关 emtpemtp中有时控开关、间隙开关、中有时控开关、间隙开关、tacstacs开关开关（二极管、闸晶管、（二极管、闸晶管、tacstacs信号控制开关）和信号控制开关）和测量开关。时控开关主要用于模拟断路器、隔测量开关。时控开关主要用于模拟断路器、隔离开关和短路。测量开关用于输出支路的电流、离开关和短路。测量开关用于输出支路的电流、功率或者向功率或者向tacstacs传递变量。为了计算操作过电传递变量。为了计算操作过电压，还准备了统计开关和规律化开关。为了详压，还准备了统计开关和规律化开关。为了详细模拟开关的物理特性，需要在开关上增加其细模拟开关的物理特性，需要在开关上增加其他的元件

25、，如断口间的杂散电容、电弧电阻。他的元件，如断口间的杂散电容、电弧电阻。 emtp的元件的元件nemtp的介绍的介绍武汉大学电气工程学院武汉大学电气工程学院33(9)控制系统（控制系统（tacs或或models) tacs用功能块用功能块(传递函数、特殊装置和信号源传递函数、特殊装置和信号源)表表示控制系统。由于示控制系统。由于tacs的解要在下一时步传递给的解要在下一时步传递给网络，以及在网络，以及在tacs内部的时延，有时会引起数值内部的时延，有时会引起数值不稳定和解的不正确。另外，需用户输入不稳定和解的不正确。另外，需用户输入tacs变变量初始值。笔者的量初始值。笔者的tacs改良版已解

26、决这些问题。改良版已解决这些问题。 models用类似用类似fortran的高级语言表示控制系统的高级语言表示控制系统,更灵活，更直观。但也存在时延和初始值的问题，更灵活，更直观。但也存在时延和初始值的问题，且所需的计算时间多于且所需的计算时间多于tacs。emtp的元件的元件nemtp的介绍的介绍武汉大学电气工程学院武汉大学电气工程学院34（）从概念上讲，（）从概念上讲，emtp可应用于任何电路的可应用于任何电路的电磁暂态现象计算。因此不仅电力系统，其他电磁暂态现象计算。因此不仅电力系统，其他的领域也可以应用，如用于磁浮铁路导轨线卷、的领域也可以应用，如用于磁浮铁路导轨线卷、核反应离子束发生装置的暂态现象仿真。核反应离子束发生装置的暂态现象仿真。（）（）emtp也可应用于电力系统的稳态计算。也可应用于电力系统的稳态计算。（）发电机轴扭振现象（）发电机轴扭振现象(ssr)的模拟是的模拟是emtp应用的扩张。应用的扩张。（）从概念上讲，（）从概念上讲，emtp不适合电磁场计算。不适合电磁场计算。 emtp的应用范围的应用范围nemtp的介绍的介绍武汉大学电气工程学院武汉大学电气工程学院35