范文一:电力系统电磁暂态分析
电力系统电磁暂态分析
绪 论
一、 电力系统暂态过程概述
电网参数
?元件参数: 发电机、 变压器、 线路等各个元件的属性参数( R 、 L 、 C 、 K )
运行参数: 反映元件当前运行状态的电气或机械? 参数( U 、 I 、 φ 、 ω )。
电网的运行状态
?从控制管理的角度分类: 正常、 警戒、 紧急、 崩溃、恢复。
?按运行参数的变化来分类:稳态(正常)、暂态(故障)。 暂态的起因:突然的扰动(短路、设备故障、雷电、操作等)。
暂态过程包括二个过程:
?电磁暂态: 线路变压器电压电流等电气运行参数
量的快速变化;(故障分析、 元件保护快速切除);
?机电暂态: 电机角位移、 角速度等机械运行
参数量的慢速过程。(稳定性分析,电网安稳装
置)。
电力系统分析:
?电力系统稳态分析
?电力系统暂态分析(电磁暂态分析、 机电暂态分
析)
a) 电力系统电磁暂态分析(电力系统故障分析):
研究交流电力系统发生短路后电磁暂态过程。
b) 电力系统机电暂态分析(电力系统稳定性):
电力系统受到各种扰动后的机电暂态过程。
研究假设: 快速与慢速过程的解耦,突出关键和便于分析
?电磁暂态分析中一般忽略机电参数变化; ?机电暂态分析则对电磁参数作近似假设。 二、 电力系统暂态过程的研究方法
?现场试验: 风险大, 费时费事, 非必要时不采用;
?模拟试验: 动模试验(物理模拟)、 数学模拟(直流计算台、交流计算台)、暂态网络分析仪
( TNA ,
电力系统物理模型与计算机技术结合产物)。 ?数字仿真: 数学模型、 计算方法、 程序编制、 问题求解。 EMTP ( Electro-magnetic Transient
Program )( UBC 、 BPA 、 EPRI/DCG ); 实时数 字仿真 RTDS ( Real Time Digital Simulation )。
范文二:电力系统电磁暂态分析总复习
《电力系统电磁暂态分析》总复习
第一章 故障分析的基本知识概念
1、 扰动(故障)引起电磁暂态、机电暂态。电磁暂态
是电压电流等电气运行参数的快速变化过程。机电
暂态是角速度等机械运行参数的慢速变化。电力系
统电磁暂态分析是研究交流电力系统发生短路(断
线)后的电磁暂态过程。
2、 故障类型:4种短路类型、 2种断线类型。 对称故障
(三相短路) 、 不对称故障 (不对称短路、 断线故障) 。
不对称短路(横向故障) 、断线故障(纵向故障、非
全相运行) 。简单故障、复杂故障。
3、 B S 相同,但 B B B U I Z 、 、 不一样
① 精确计算法:选定某电压等级的基准电压, 然后
按变压器实际变比归算到各个电压等级,再在各个
电压等级把有名值归算为标幺值。
② 实用计算方法 (近似计算法) :取变压器变比为 平均电压之比。
B B B B B B B B B 2B B B S() U () () () S , US S U B
B B B MVA KV I KA Z I I U I Z S Ω===av 、 、 、 满足关系:
选定:、 =U则:
③ 元件参数的归算
不同基准容量和电压的元件标幺值归算(发电机、 变压器、电抗器等)
4、 无限大电源供电的三相突然短路分析。
1 短路电流含有二种分量:稳态基频分量、直流自
由分量。
2 稳态分量与短路前相比增大,直流自由分量逐渐
衰减为 0(出现原因:电感电流不能突变) 。
3
4 电流瞬时值。在短路后 0.01秒出现。
(3) 1.81.9f
imp imp imp i I K =冲击系数 =~ 5 有效值:为时刻 t 为中心一个周期内瞬时电流的
均方根。最大短路电流的有效值为:
(3)I f imp =
6 短路功率:即是某支路的短路电流与额定电压构
成的三相功率。 *
(3)(3)* ==f f D av D B B S I S S I S =2*(B)*() *() 2*(B)*() *() X XN N B B N N N B N B B B B N av N N N
N U U S I X X S U I U S I U U U X X S I ????????= ? ????? ? ?????????
????=== ? ?????=当 时(实用计算中):=
第二章 同步发电突然三相短路分析
1、 Park 变换:
① 一 种 旋 转 坐 标 变 换 :abc dqo ,
abc dqo 基频交流 直流
② Park 变换引入的意义:由于定子、转子之间存 在相对运动, 定子各个绕组的磁路会发生周期性
的变化,故其电感系数(自感和互感)或为 1
倍或为 2倍转子角的周期函数 (转子角本身是时
间的三角周期函数) ,故电压 (磁链 ) 方程是一组 变系数 的微分方程。通过 Park 变换,定子的静
止三相绕组被三个随转子旋转的等效绕组 dd 绕
组、 qq 绕组、 00绕组所代替, 由于这些绕组 (包
括原来转子绕组) 之间均不存在相对运动, 故在
Park 坐标系下的各个绕组的自感和互感均为常
数,故电压 (磁链 ) 方程是一组 常系数 的微分方
程。
③ 各个绕组自电抗(0d q f D Q x x x x x x 、 、 、 、 、 ) 和互电抗 ad aq x x 、 (电枢反应电抗) 。
2、 同步发电机的稳态计算
基本方程 (来自于 Park 形式的磁链电压方程的简 化) :
() 0() 0 0q q d q d d q d q q d q d q d q q ad f
E U r I jx I U r I jx I r E U x I U x I E x I ?=++???=++?
=?=+???=-?
=
当忽略电阻 时可写成标量形式:
注意:发电机空载电势
实用计算中(已知 q U I E → 、 ) : ① 隐级机:d q x x =,则复数(向量)直接运算:
() q d d E U r I jx I U jx I =++≈+
② 凸 级 机 :
d q x x ≠则 :计 算 虚 构 电 势 :() Q q q E U r I jx I U jx I =++≈+ 可 得 到
q U δ
轴与 间的夹角 。
() () sin() q Q d Q d q d q E E x x I E x x I ?δ=+-=+-+
cos sin() q d d q q q d d d E U rI x I U x I U x I δ?δ=++≈+=++
3、
① 无阻尼绕组:
由于暂态电势 ' q E 正比于励磁绕组磁链,故在短路 前后保持不变,机端短路时 0d q U U U ===,故:暂态电流(电流基频分量的起始值)
' I' ' q d E x =
② 有阻尼绕组:
由于次暂态电势 q
d E E E ''''''、 、 在短路前后均保持不 变,机端短路时 0d q U U U ===,故:次暂态电流 (电流基频分量的起始值) :
d q
d d q E E x x === 0q d q d q d q E U I x U j I x ?''=+???=-
?
q d q d d q d q E U I x E U j I x ?''''=+???''''=-
?
③ 稳态电流
如果励磁调节器没有动作,则 q E 与短路前相等。
4、 同步发电机的三相短路暂态过程
量,电磁感应产生的为自由交流分量
① 定子绕组:
? 基频分量。 强制分量、暂态分量(d T ') 、 d 轴
次暂态分量(d T '') 、 q 轴次暂态分量(q T '') 。
? 直流分量(a T )
? 倍频分量(a T )
② 励磁绕组:
? 直流分量。 强制分量、自由直流分量(d T ') 、
自由直流分量(d T '') 。
? 基频分量(a T ) 。
③ D 绕组:
? 直流分量 (d T '') 。
? 基频分量(a T ) 。
④ Q 绕组:
? 直流分量 (q T '') 。
? 基频分量(a T ) 。
Iq
d
E x ∞=定子稳态电流:
5、 定子基频分量
6、 定子全电流
I I I e I I e I I I e
I I --∞∞-?=-+-+???=?=≈
第三章 三相短路的实用计算
1、 起始次暂态电流 I ''计算 (短路电流基频交流分量的
初始值 ) 。 (冲击电流、短路电流最大有效值、短路 容量) :等值法 和 叠加法。
注意:确定系统各个元件的次暂态参数(发电机、 电动机、 线路、 变压器、 负荷) 选取 B B av S U U =、 计 算等值网络中各个元件的参数 (正序) 。 另外注意网 络变换及其化简(串联、并联、星三角形变换、多
个有源支路的并联等)
2、 任一时刻短路电流计算 (交流分量的有效值) :运算 曲线法 。计算步骤:
① 制定等值电路(发电机电抗用 d x ,略去电阻、 并联支路以及负荷) ;
② 网络变换化简求取各个等值发电机对短路点的 转移电抗; (可合并远离短路点的同类型电源) ③ 将转移电抗转化为计算电抗。 (无限大电源不需 要) 。
④ 查运算曲线求得各个电源的短路电流的标幺值。 (无限大电源的短路电流为转移电抗的倒数) ⑤ 将查得的短路电流标幺值按电源基准容量以及 平均电压转换为有名值后相加。
3、 计算机求解法。叠加法
① 先作潮流计算, 求节点电压和支路电流的正常分 量;
② 用节点导纳矩阵方程,在 f 节点注入单位电流, 求解线性方程,得到 f 列的阻抗矩阵元素。 ③ 用 f 点的自阻抗元素求故障点次暂态电流。 ④ 用节点 i 与 f 点之间的互阻抗元素,求节点 i 的 节点电压的故障分量;
⑤ 合并节点电压的正常分量和故障分量得到实际
节点电压;
⑥ 求得各个支路的次暂态电流。
第四章 对称分量法及其各个元件的序参数和序网
1、 对称分量法:
2、 各个电网元件的正、负、零序阻抗
① 旋 转 元 件 (发 电 机 、 电 动 机 、 调 相 机 ) :
(0)(2)(1)x
x x ≠≠
② 静 止 磁 耦 合 元 件 (输 电 线 、 变 压 器 ) : (0)
(2)(1)x x x ≠=
③ 静 止 无 磁 耦 合 元 件 (电 抗 器 、 电 容 器 ) : (0)(2)(1)x
x x ==
3、 变压器的零序电抗和等值电路
① 0(/) YN d Y ?、 :零序电抗为 ) 2() 1(I (0)x x x x x II ==≈ +, ?侧相当于接地。
② 0(/) YN Y Y 、 y 零序电抗为
∞≈=) 0(I (0)x x m x +, y 侧开 路。
③ 00(/)
YN n Y Y 、 y :零 序 电 抗 为 ) 2() 1(I (0) x x x x x II ==≈ +。 Yn 侧无通路则开路,有
通路则接入外接电抗接地。
④ YN 、 d 、 y (Y Y //0
?) :类似于二绕组分析 ⑤ YN 、 d 、 yn (00//Y Y ?) :类似于二绕组分析 ⑥ YN 、 d 、 d (??//0
Y ) :类似于二绕组分析 4、 输电线零序阻抗的大小取决于输电线路的结构和参
数。与正负序阻抗相比的大小变化。单回架空输电 线、双回架空输电线、架空地线的零序阻抗的大小 变化比较。
5、 零序网的绘制
① 首先在故障点画出零序电源; 查明零序电流可以 流过的通路作为画出零序网的依据 (注意:中性 点是否接地、绕组接线方式) 。
② 从短路点开始, 将零序通路上的各个元件 (变压 器、 发电机、 线路等) 用其零序等效电路代替 (不 通零序电流的不画出) (注意中性点的接地阻抗 画出时要取实际阻抗的 3倍) 。
第五章 不对称分量法故障的计算
1、 不对称故障短路的计算分析
① 三类不对称故障的复合序网(含经阻抗短路) 。 其中包含了正序 /负序 /零序网络、等值序阻抗、 开路电压、 三个序电压与序电流关系方程、 边界 条件等方面的内容。 基于该复合序网求取序电流 和序电压以及相量图。
② 掌握三种不对称短路的短路电流大小 (与三相短 路比较) ,故障和非故障相电压大小比较(与正 常电压比较) 。
③ 掌握正序等效定则以及应用计算曲线求取故障 处正序电流(以及短路电流) 。
④ 掌握利用节点导纳矩阵求取节点阻抗矩阵 f 列 元素的方法。 并用自阻抗 (等值序阻抗) 求故障 处序电流、 序电压, 用互阻抗求任意节点序电压 以及支路序电流。 合成相电压和相电流要注意可 能的相位变化。
2、 非全相运行的计算分析
① 处理思想:类似于不对称短路的分析方法, 把断 口 (不对称处) 的不对称电压、 电流分解为三序 分量, 用对称分量法形成针对断线断口的三个序 网等值电路,再结合边界条件联立求解即可。
② 基于断口开路电压的求解方法
③ 基于潮流结果的求解方法(叠加法)
④ 基于节点导纳矩阵 (节点阻抗矩阵) 求取等值序 网参数的方法, 求取任意节点对断线端口的转移 阻抗, 求取任意节点序电压和任意支路序电流的 方法。合成相电压和相电流的方法。
范文三:电力系统电磁暂态分析0-绪论
电气工程学院
学习目标:掌握对称分量法。
内容提要:
电力元件和电网运行状态,短路后的电流和电压。
一次元件:
1、发电设备,指同步(隐极和凸极)、异步、双馈发电机和光伏 发电等电源
2、变电设备,指各种类型的变压器
小型:0.63MVA及以下;中型:0.8~6.3MVA;大型:8~63MVA;特大型:90MVA及以上(绕组的电阻、电抗比值)
自耦和普通变压器(等效电路)
三绕组和两绕组(铁芯)
三相和由单相构成的三相(铁芯)
升压和降压(绕组排列)
3、输配电设备,包括输电线路、母线、断路器等
电缆和架空线路
多股导线、扩径导线、分裂导线(载流面积、电晕损耗、线路电抗) 钢筋混凝土杆塔和铁塔
换位(三相平衡)
柔性交流输电(大容量、非线性的电力电子元件)
直流输电系统
4、用电设备,指各种负荷
综合负荷P+jQ
动态模型 ,是t、U、f、dU/t、 df/t的函数,适用于U、f急剧变化时的情况。 静态模型 ,是U、f的函数,适用于f变化缓慢时的情况。
二次元件:
对一次系统起调节、控制功能的元件,如励磁、调速、继电保护等。
以 运行参数 表示:
1、功率S、P、Q,电压U,电流I,相位θ,频率f 。
2、断路器、隔离刀闸等的状态。
3、包含运行参数的函数或集合。
受 元件参数 的影响,包括电阻、电抗、电导、输入阻抗、转移阻抗、 变压器变比、时间常数、放大倍数等。
稳态 指系统运行参数不变化。
实际当系统运行参数在很小范围变化时,通常认为电力系统处于稳态。 暂态 指系统运行参数连续变化。
电力系统从稳态进入暂态的原因是遭受 扰动 的影响。 扰动包括故障和正常操作。
机电暂态
基于输入机械功率和输出电磁功率的平衡情况,分析发电机转子的运动规律,考 察机组和外部系统之间维持同步、联合运行的能力,又称为 稳定性问题 。
一般会将系统电路部分进行等值简化处理,重点研究电磁参数和转子运动参数的 关系。
电磁暂态
研究一次系统中电磁部分的运行参数变化规律,主要求解暂态过程中的电流、电 压相量。
由于电磁参数的变化速度大大快于机械运动参数的变化速度,一般忽略电磁暂态 与机电暂态之间的相互影响,假设全系统的频率不变,各发电机电动势相量也不 变。
电力系统动态模拟
暂态电流、电压的变化过程
1、物理模拟,指电力系统动态模拟系统
2、数字模拟,基于计算机的模拟系统,有实时和非实时的区分
电力系统故障分析
求解某时刻的暂态电流、电压相量
包括人工计算方法和计算机软件计算方法
动态非线性模型
用于电力系统动态模拟,模型复杂,计算难度相对较大
静态线性模型
用于电力系统故障分析,模型简单,计算难度相对较小
1、严格遵循实际电网的特性建模、求解会使计算十分复杂,效率低下。
2、电力系统中存在不确定性因素,合理简化和近似计算不可避免。
3、实际工程中,允许范围内的误差不会影响后续研究的正确性 。
范文四:电力系统电磁暂态概述
第 1章 电力系统电磁暂态概述
1.1 电力系统电磁暂态现象 .................................................................................... 1.2 电力系统电磁暂态分析的目的 .......................................................................... 1.3 电力系统电磁暂态研究的方法 ........................................................................ 1.4电力系统电磁暂态仿真的特点 ........................................................................... 1.5电力系统数字仿真 ............................................................................................... 思考与练习题
1.1 电力系统电磁暂态现象 . ......................................................................................... 2 1.2 电力系统电磁暂态分析的目的 . ............................................................................. 4 1.3 电力系统电磁暂态研究的方法 . ............................................................................. 5 1.4 电力系统电磁暂态的特点 . ..................................................................................... 7 1.4.1 频率范围广 . ............................................................................................................ 7 1.4.2 元件模型因计算目的而异 . .................................................................................... 8 1.4.3 行波现象和分布参数 . .......................................................................................... 10 1.4.4 非线性元件和开关操作 . ...................................................................................... 16 1.4.5 元件参数的频率特性 . .......................................................................................... 17 1.4.6 时间跨度的要求 . .................................................................................................. 18 1.5 电力系统数字仿真 . ............................................................................................... 18 1.5.1 电力系统数字仿真的分类 . ................................................................................ 18 1.5.2 电力系统数字仿真的优点 . ................................................................................ 20 1.5.3 电力系统数字仿真软件 . .................................................................................... 21
第 1章 电力系统电磁暂态概述
1.1电力系统电磁暂态现象
电力系统稳态运行时, 发电厂发出的功率与用户所需要的功率及电网中损耗 的功率相平衡, 系统的电压和频率都是稳定的。 但电力系统在运行过程中常常会 发生故障或需要进行操作,常见的电力系统故障有:雷击电力设备等雷害故障, 短路、接地故障和谐振等电气故障,断线等机械故障。常见的电力系统操作有: (1)断路器的投切操作,如合空载线路、合空载变压器、切空载线路、重 合闸、甩负荷等。
(2)隔离开关的投切操作,如母线投切等。
电力系统发生故障或进行操作时, 系统的运行参数发生急剧变化, 系统的运 行状态有可能急促地从一种运行状态过渡到另一种运行状态, 也有可能使正常运 行的电力系统局部甚至全部遭到破坏, 其运行参数大大偏离正常值, 如不采取特 别措施, 系统很难恢复正常运行, 这将给国民经济生产和人民生活带来严重的后 果。
电力系统运行状态的改变, 不是瞬时完成的, 而要经历一个过渡状态, 这种 过渡状态称为暂态过程。 电力系统的暂态过程通常可以分为电磁暂态过程和机电 暂态过程。 电磁暂态过程指电力系统各元件中电场和磁场以及相应的电压和电流 的变化过程, 机电暂态过程指由于发动机和电动机电磁转矩的变化所引起的电机 转子机械运动的变化过程。
虽然电磁暂态过程和机电暂态过程同时发生并且相互影响, 但由于现代电力 系统规模的不断扩大, 结构愈益复杂, 需要考虑的因素繁多, 再加上这两个暂态 过程的变化速度相差很大, 要对它们统一分析是十分复杂的工作, 因此在工程上 通常近似地对它们分别进行分析。 例如, 在电磁暂态过程分析中, 由于在刚开始 的一段时间内, 系统中的发电机和电动机等转动机械的转速由于惯性作用还来不 及变化, 暂态过程主要决定于系统各元件的电磁参数, 故常不计发动机和电动机 的转速变化, 即忽略机电暂态过程。 而在静态稳定性和暂态稳定性等机电暂态过 程分析中, 转动机械的转速已有了变化, 暂态过程不仅与电磁参数有关, 而且还 与转动机械的机械参数(转速、角位移)有关,分析时往往近似考虑或甚至忽略 电磁暂态过程。 只在分析由发动机轴系引起的次同步谐振现象、 计算大扰动后轴
系的暂态扭矩等问题中,才不得不同时考虑电磁暂态过程和机电暂态过程。
下面以一个简单开关接通 RL 电路的例子, 以便获得对在电力系统暂态时起 关键作用的物理过程的充分了解。
一个正弦波电压接通到一个电感与电阻串联的电路上,如图 1-1所示。这实 际上是一个高压断路器闭合到短路的输电线路或短路的电缆的最简单单相表示 法。电压源 ) sin() (θω+=t E t e m 代表连接的同步发电机的电动势。电感 L 包括发 电机的同步电感、电力变压器的漏电感与母线、电缆与输电线的电感,电阻 R 表示供电电路的电阻损耗。
图 1-1 正弦波电压源接通到 RL 串联电路
假设 s 0=t 时合闸,应用基尔霍夫电压定律,得到电路方程
) sin(θω+=+t E dt
di L Ri m (1-1) 该方程为一阶常系数、 线性、 非齐次常微分方程, 其解就是合闸电路的全电 流,它由两部分组成:稳态分量和暂态分量,即
aP P i i i +=
其中稳态分量为
) tan sin() (122R L t L R E i m
P ωθωω--++=(1-2)
暂态分量,也就是合闸电流的自由分量 , 记为
a T t pt aP Ce
Ce i -==(1-3) 式中, L R p -=,为特征方程 0=+pL R 的根; R
L p T a =-=1, 为暂态分量电 流衰减的时间常数; C 为由初始条件决定的积分常数。
在开关闭合之前, 电感 L 中的磁通为 0, 根据磁通守恒定律,
在闭合的瞬间,
即
0) tan sin() (122=-++
-R L L R E C m ωθω
由此得到 ) tan sin() (122R L L R E C m ωθω--+-=
(1-4)
从而得到合闸的全电流表达式为 t
L R m
m
e R L L R E R L t L R E t i ) /(122122) tan sin() () tan sin() () (----+-+-++=ωθωωθωω(1-5)
式(1-5)中的暂态分量含有 t L R e ) /(-衰减项,也称为直流分量,其系数为常 数, 数值大小取决于电流合闸瞬间, 在 ) tan (1R L
ωθ--为 πk (其中 2, 1, 0=k ) 时,
直流分量为 0,电流立即进入稳态,换言之,不存在暂态振荡过程。但当开关闭 合电路不在 ) tan (1
R L ωθ--为 πk (其中 2, 1, 0=k )时,合闸过程将引起电磁暂态 过程,在 ) tan (1
R L ωθ--为 ?90时,暂态过程将达到最大电流,如图 1-2所示。
图 1-2 开关合闸的暂态过程电流波形
1.2电力系统电磁暂态分析的目的
电磁暂态过程分析的主要目的在于分析和计算故障或进行操作后可能出现 的暂态过电压和过电流, 以便对电力设备进行合理设计。 通常情况下,
电力系统
电磁暂态产生的过电压在确定设备绝缘水平中起决定作用, 据此制定高电压试验 电压标准,确定已有设备能否安全运行,并研究相应的限制和保护措施。此外, 对于研究电力系统新型快速保护装置的动作原理及其工况分析, 故障测距原理与 定点方法以及电磁干扰等问题, 也常需要进行电磁暂态过程分析。 另外, 调查事 故原因,寻找对策;计算电力系统过电压发生概率,预测事故率;检查电气设备 的动作责能, 如断路器的暂态恢复电压和零点偏移; 检查继电保护和安全自动装 置的响应等,也离不开电磁暂态过程的计算和模拟。
电磁暂态过程变化很快,一般需要分析和计算持续时间在 ms 级,甚至是 μs 级以内的电压、电流瞬时值变化情况,因此,在分析中需要考虑元件的电磁耦 合, 计及输电线路分布参数所引起的波过程, 有时甚至要考虑三相结构的不对称、 线路参数的频率特性以及电晕等因素的影响。
1.3电力系统电磁暂态研究的方法
为了保证电力系统运行的可靠性、 安全性和经济性, 在电力系统设计、 运行、 分析和研究中必须全面地了解实际系统的电磁暂态特性。 目前, 研究电力系统电 磁暂态过程的手段有 3种:
(1)系统的现场实测方法。
(2)应用暂态网络分析仪(Transient Network Analyzer ,简称 TNA )的物 理模拟方法。
(3)计算机的数字仿真(或称数值计算)方法。
系统的现场实测方法是在实际的电力系统上直接进行试验和研究, 六十年代 之前经常要在实际电力系统进行短路、 操作等试验, 这种试验对电力系统的考验 是真实和严格的, 以确保电力系统运行的可靠性、 安全性和经济性, 但是系统的 现场实测方法会对电力系统的正常运行和电气设备带来很大危害, 短路点的电弧 有可能烧坏电气设备, 很大的短路电流通过设备会使发热增加, 当持续一定时间 后, 可能使设备过热而损伤; 很大的短路电流引起的电动力有可能使设备变形或 遭到不同程度的破坏; 操作试验过程中产生的过电压可能引起电气设备载流部分 的绝缘损坏, 加剧绝缘材料的老化。 即便如此, 实测对于研究电力系统电磁暂态 过程仍是非常重要的,它一方面验证 TNA 及数字仿真的准确性,为系统安全运 行提供依据;另一方面可以全面研究系统各类元件的参数特性,为 TNA 及数字
仿真提供更精确的原始数据。
系统的现场实测常常会遇到困难, 有些困难甚至是不可能解决的, 利用模型 系统进行试验和分析就成为一种非常有效的途径。 暂态网络分析仪就是一种用于 研究电力系统动态特性的物理模型系统。 TNA 方法多用于模拟操作过电压和交 流过电压的暂态现象,同时通过改变元件特性, TNA 也可用来模拟更高频率下 的暂态现象。 它是在相似理论的指导下, 把实际电力系统的各个部分, 如同步发 电机、变压器、输配电线路、电力负荷等按照相似条件设计、建造并组成一个电 力系统模型,这样将一个高电压、大电流、体积庞大的电力系统,按照一定的比 例转化为一个低电压、 小电流、 体积较小的模拟试验台, 在模拟台中出现的电磁 暂态现象, 电压和电流的波形与它模拟的电力系统是一样的, 用这种模型代替实 际电力系统进行各种正常运行与故障状态的试验和分析。 与系统的现场实测相比, TNA 方法对电力系统的正常运行和电气设备不产生影响,为了缩小模拟装置的 尺寸, 节省电感元件和电容元件, 减少模拟设备的昂贵费用, 并考虑到现有的技 术条件、模拟精度要求等,选择恰当的比例尺是非常重要的。 TNA 具有物理意 义清晰,易于理解和使用的优点,可以多次重复试验现象,便于观察和研究,北 美不少大的电力公司都将 TNA 作为培训新员工的一种工具。
随着现代电力系统的发展, 电力系统的规模和复杂程度发生很大变化, 采取 物理模型的动态模拟方法受到很大限制。 与此同时, 数字计算机和数值计算技术 飞速发展, 数字计算机的性能价格比不断提高, 出现了用数学模型代替物理模型 的新型模型系统。电力系统数字仿真(Digital Simulation of Power System)就是 将电力系统的电源、 网络和负荷元件建立其数学模型, 用数学模型在数字计算机 上进行实验和分析的过程。 电力系统数字仿真的主要步骤为建立各元件数学模型、 建立数字仿真模型和进行仿真试验。 建立数学模型是处理物理原型与数学模型之 间的关系, 建立数字仿真模型是处理数学模型与计算之间的关系。 有些数学模型 是利用数字计算机和模拟计算机的混合数学模型系统。 电力系统数字仿真是一门 新兴的技术科学, 它的产生和发展是同现代科学技术发展分不开的, 数字仿真与 实际系统试验和动态物理模拟相比, 不仅节省了大量的人力、 物力和财力, 而且 不受外部条件的限制, 几乎不受系统规模和时间跨度的约束, 甚至不受各种暂态 现象频率范围的限制 (理论上它可以对各类暂态过程进行计算, 但是, 它需要相
关设备真实的频率特性,有时候,这种频率特性是很难得到的) 。具有无可比拟 的灵活性,能达到试验不可达到的广度和深度。譬如我国南北联网这样的课题, 地理上相距数千公里, 跨越了几个大电网, 没有办法用试验来分析联网可能出现 的问题, 但通过数字仿真发现南北联网可能会出现低频震荡问题。 今天实际系统 的现场实测方法主要是为了建立数学模型,取得数学模型的参数。
1.4电力系统电磁暂态的特点
1.4.1频率范围广
电力系统中暂态现象的研究所涉及的频率范围广,从直流到大约 50MHz 的 范围。 高于系统频率的暂态现象通常涉及到电磁暂态, 而低于系统频率的暂态现 象主要涉及到机电暂态过程。 表 1-1给出了多种暂态现象的起因以及它们通常的 频率范围。
表 1-1电力系统暂态的起因及频率范围
网络中每个元件的模拟都要与所研究的特定暂态现象的频率范围相符合。 当所研究现象的频率大于 1 MHz 时,如 GIS 中由于隔离开关操作所引起的快速 暂态现象, 则不仅在母线上产生波的传播, 而且施加在变压器、 支柱绝缘子以及 在某些情况下管形母线上的弯管处, 它们非常小的电容和电感对模拟结果都将产 生非常重要的影响。
表 1-1中所列电磁暂态现象的频率范围可以分成 4组,对应于各暂态现象 的频率范围之间存在着重叠,图 1-3是国际大电网会议(CIGRE )对各种过电压 的频率分类; 各类的频率范围是与其所表示的过电压波形的实际陡度相关的。 研 究者必须清楚自己的研究对象所在的频率范围, 确定被模拟设备的频率特性, 只 有这样,才能得到满意的电磁暂态分析计算结果。
图 1-3各种过电压的频率范围
通常,频率越高,所考虑的现象 (如过电压 ) 在时间上空间上的衰减越快,因 此所考虑的物理范围 (模拟范围 ) 越小,模拟时间越短。相对地,在工频或与此接 近的频率领域,为了掌握现象的性质,需要大范围长时间的模拟。图 1-4表示电 力系统数字仿真中各种计算所考虑现象的时间幅度和计算涉及的系统规模。
图 1-4电力系统各种现象的变化速度和计算范围
1.4.2元件模型因计算目的而异
电力系统由各种不同的元件所组成,
元件的动态特性对于系统的暂态过程有
0.11101001k 10k 100k 1M 10M
100M
0.1Hz-3kHz
Ⅰ暂时过电压
Ⅱ缓波前过电压
Ⅲ快波前过电压
Ⅳ特快波前过电压 频率 (Hz)
1
10100
1微秒 1分 1毫秒 1秒 1小时 系统规模 (发电机数 )
直接的影响。为此,首先需要研究各元件的动态特性,建立它们的数学模型。在 此基础上, 根据系统的具体结构, 即各元件之间的相互关系, 组成全系统的数学 模型, 然后采用适当的数学方法进行求解, 这便是电力系统暂态分析的一般方法。 然而, 由于各元件的动态响应有所不同, 系统各种暂态过程的性质也不相同。 因此, 在不同目的暂态过程分析中, 所考虑的元件种类和对它们数学模型的要求 并不相同。 例如, 在电磁暂态过程分析中, 所研究的暂态过程持续时间通常较短, 在此情况下, 一些动态响应比较缓慢的元件, 如原动机及调速系统等的影响往往 可以忽略不计, 而发电机定子回路和电力网中的电磁暂态过程则需加以考虑。 相 反, 在电力系统稳定性分析中, 则通常忽略发电机定子回路和电力网中的电磁暂 态过程, 而将线路和变压器等元件用它们的等值阻抗来描述。 另外, 就同一种系 统暂态过程来说, 对于不同的分析精度和速度要求, 元件所用数学模型的精确程 度也不相同。 一般地说, 在进行系统规划和设计时, 暂态分析的精度要求可以适 当降低,这时各元件可以采用较粗略的数学模型,以便提高分析速度。因此,在 建立元件数学模型时, 不但需要研究它们的精确模型. 而且需要考虑各种简化模 型,以适应不同的需要。
在建立元件模型时还必须注意研究对象所处的频率范围。例如,在计算交 流过电压时, 变压器采用通常的以互感及绕组漏感和电阻表示的模型, 如图 1-5(a)所示; 但在计算操作过电压时, 除了上述要素外, 还需要考虑绕组的对地电容和 端子间电容及绕组间电容,如图 1-5(b)所示;而在计算雷过电压时,变压器模型 通常用冲击电容表示,无需考虑电感和电阻要素,如图 1-5(c)所示。当研究现象 的频率很高时, 变压器和互感器的杂散电容、 引线的微小电感对计算结果都有举 足轻重的影响。
C 12/2
C C 12
(a) 交流过电压计算用 (c) 雷过电压计算用
(b) 操作过电压计算用
图 1-5变压器的模型
1.4.3行波现象和分布参数
电力系统采用长线路将能源中心发出的电能输送给各电力用户, 长线路的具 体形式有架空输电线路和电缆线路两种,每微段的线路都呈现自感和对地电容, 即线路是具有分布参数特性的电路元件。 当电力系统中某一点突然发生雷电过电 压或操作过电压时, 这一变化并不能立即在系统其他各点出现, 而是以电磁波的 形式按一定的速度从电压或电流突变点向系统其他部位传播。 例如, 当架空输电 线路遭受雷击时, 雷击点导线将产生雷电过电压, 该电压将沿着导线向两侧传播。 这个沿线路传播的电压以及与其相伴而行的电流波称为行波。 当行波到达变电站 或其他节点时, 由于电路参数的改变, 将引起波的折射和反射。 这种在分布参数 电路中产生的暂态过程本质上是电磁波的传播过程,简称波过程。
实际电力系统采用三相交流或双极直流输电, 属于多导线线路, 而且沿线路 的电场、 磁场和损耗情况也不尽相同, 因此所谓的均匀无损单导线线路实际上是 不存在的。 但为了揭示线路波过程的物理本质和基本规律, 可暂时不考虑线路的 电阻和电导损耗, 并假定沿线线路参数处处相同, 即首先研究均匀无损单导线中 的波过程。
1. 波传播的物理概念
假设有一无限长的均匀无损单导线,见图 1-6(a), t =0时刻合闸直流电源, 形成无限长直角波,单位长度线路的电容、电感分别为 C 0、 L 0,线路参数看成
是由无数很小的长度单元 x ?构成,如图 1-6(b)所示。
图 1-6均匀无损的单导线
(a )单根无损线首端合闸于 E ; (b )等效电路
合闸后,电源向线路电容充电,在导线周围空间建立起电场,形成电压。靠 近电源的电容立即充电, 并向相邻的电容放电, 由于线路电感的作用, 较远处的 电容要间隔一段时间才能充上一定数量的电荷, 并向更远处的电容放电。 这样电 容依次充电, 沿线路逐渐建立起电场, 将电场能储存于线路对地电容中, 也就是 说电压波以一定的速度沿线路 x 方向传播。 随着线路的充放电将有电流流过导线 的电感, 即在导线周围空间建立起磁场, 因此和电压波相对应, 还有电流波以同 样的速度沿 x 方向流动。 综上所述, 电压波和电流波沿线路的传播过程实质上就 是电磁波沿线路传播的过程, 电压波和电流波是在线路中传播的伴随而行的统一 体。
2. 波动方程及解
为了求出无损单导线线路行波的表达式, 令 x 为线路首端到线路上任意一点 的距离。 线路每一单元长度 d x 具有电感 L 0d x 和电容 C 0d x , 如图 1-7所示, 线路 上的电压和电流 i 都是距离和时间的函数。
图 1-7均匀无损单导线的单元等值电路
根据节点电流方程 ∑=0i 可知,
x x
i i t u x C i d d 0??++??= 根据回路电压方程 ∑=0u 可知,
x x
u u t i x L u d d 0??++??= 整理得
00=??+??t
u C x i (1-6) 00=??+??t i L x u (1-7) 由方程 (1-6)对 x 再求导数,由方程 (1-7)对 t 再求导数,然后消去 i ,并用类似 的方法消去 u 得。
220022t
u C L x u ??=??(1-8) 220022t
i C L x i ??=??(1-9) 其中 00C L 、 为单位长度电感和电容。
通过拉普拉斯变换将 u (x , t ) 变换成 U (x , S ) , i (x , t ) 变换成 I (x , S ) ,并假 定线路电压和电流初始条件为零, 利用拉氏变换的时域导数性质, 将式 (1-8)、 式 (1-9)变换成
0) , () () , (222=-??S x U S R x
S x U (1-10) 0) , () () , (222=-??S x I S R x
S x I (1-11) 其中 v
S S R ±=) (。 根据 2阶齐次线性微分方程性质,令 0
01C L v =
,则式(1-10) 、 (1-11)解 为 x v S b x v
S f e S U e S U S x U ) () () , (+=-(1-12) x v S b x v S f e S I e
S I S x I ) () () , (+=-(1-13) 将以上频域形式解变换到时域形式为
) () () , (v
x t i v x t i t x i b f ++-=(1-14)
) () () , (v
x t u v x t u t x u b f ++-=(1-15) 式 (1-14)、式 (1-15)就是均匀无损单导线波动方程的解。
3. 波速和波阻抗
在波动方程中定义
01C L v = v 为波传播的速度。对于架空线路
s m v /1031
800?==εμ
即沿架空线传播的电磁波波速等于空气中的光速度。 而一般对于电缆, 波速 s m v /105. 18?≈, 其传播速度低于架空线, 因此减小电缆介质的介电常数可提高 电磁波在电缆中传播速度。
定义波阻抗
Ω=-==00C L i u i u Z b b f f
一般对单导线架空线而言, Z 为 500Ω左右, 考虑电晕影响时取 400Ω左右。 由于分裂导线和电缆的 0L 较小而 C 0较大,故分裂导线架空线路和电缆的波阻抗 都较小,电缆的波阻抗约为十几欧姆至几十欧姆不等。
波阻抗 Z 表示了线路中同方向传播的电流波与电压波的数值关系,但不同 极性的行波向不同的方向传播, 需要规定一个正方向。 电压波的符号只取决于导 线对地电容上相应电荷的符号, 和运动方向无关。 而电流波的符号不但与相应的 电荷符号有关, 而且与电荷运动方向有关, 根据习惯规定:沿 x 正方向运动的正 电荷相应的电流波为正方向。 在规定行波电流正方向的前提下, 前行波与反行波 总是同号,而反行电压波与电流波总是异号,即
Z i u Z i u b
b f f
-== 必须指出, 分布参数线路的波阻抗与集中参数电路的电阻虽然有相同的量纲,
但物理意义上有着本质的不同:
(1)波阻抗表示向同一方向传播的电压波和电流波之间比值的大小;电磁 波通过波阻抗为 Z 的无损线路时,其能量以电磁能的形式储存于周围介质中, 而不像通过电阻那样被消耗掉。
(2)为了区别不同方向的行波, Z 的前面应有正负号。
(3)如果导线上有前行波,又有反行波,两波相遇时,总电压和总电流的 比值不再等于波阻抗,即
Z u u u u Z i i u u i u b
f b f b f b f ≠-+=++= (4)波阻抗的数值 Z 只与导线单位长度的电感 L 0和电容 C 0有关,而与线 路长度无关。
4. 前行波和反行波
下面用行波的概念来分析波动方程解的物理意义。
首先讨论式 (1-10),电压 u 的第一个分量 ) (v
x t u f -。设任意波形的电压波 ) (v
x t u f -沿着线路 x 传播, 如图 1-8所示, 假定当 t =t 1时刻线路上任意位置 1x 点 的电压值为 u a ,当时间 t =t 2时刻时 (t2>1t ) ,电压值为 u a 的点到达 2x ,则应满足
v
x t v x t 2211-=-即 ) (1212t t v x x -=-
图 1-8行波运动
由于 v 恒大于 0, 且由于 (t 2>1t ) , 则 0) (12>-x x , 由此可见 ) (v
x t u f -
表示前
行波; 同样的方法可以证明 ) (v
x t u b +表示沿 x 反方向行进的电压波, 称为反行电 压波。 ) (v x t i f -, ) (v
x t i b +的证明过程类似。 为方便将式 (1-14)和式 (1-15)可写成
b f i i i +=(1-16)
b f u u u +=(1-17)
由式 (1-16)和式 (1-17)可知, 线路中传播的任意波形的电压和电流传播的前行 波和反方向传播的反行波, 两个方向传播的波在线路中相遇时电压波与电流波的 值符合算术叠加定理, 且前行电压波与前行电流波的符号相同, 反行电压波与反 行电流波的符号相反。
5. 行波的折射和反射
当波沿传输线传播,遇到线路参数发生突变,即波阻抗发生突变的节点时, 都会在波阻抗发生突变的节点上产生折射和反射。
如图 1-9所示,当无穷长直角波 E u if =沿线路 1达到 A 点时后,在线路 1上除 f u 、 f i 外又会产生新的行波 b u 、 b i ,因此线路上总的电压和电流为
??
???+=+=b f b f i i i u u u 111111(1-18)
图 1-9波通过节点的折反射 设线路 2为无限长, 或在线路 2上未产生反射波前, 线路 2上只有前行波没 有反行波,则线路 2上的电压和电流为
??
???==f f i i u u 2222(1-19) 然而节点 A 只能有一个电压电流,因此其左右两边的电压电流相等,即 21u u =, 21i i =,因此有
??
???+=+=b f f b f f i i i u u u 112112(1-20) 将 11Z i u f if
=, 222Z i u f f =, 11Z i u b
ib -=, E u if =代入上式得, ???????=+-==+=
E E Z Z Z Z u E E Z Z Z u b f βα2112121222(1-21) 其中 α, β分别为折射与反射系数。 α, β计算如式(1-22)所示。
?????????
+=+-=+=
βαβα122112212Z Z Z Z Z Z Z (1-22) 另外, 线路是用集中参数还是分布参数等值, 除跟线路长度有关, 还与暂态 过程的频率有关。设线路长 300m(约一档距 ) ,幅值为 I 的正弦波电流以光速 (300m/μs) 传播,如图 1-10所示, AB 两点间的传播时间 Δt 为 1μs 。如果是工 频 50Hz ,两点的电流的差值最大为 3.14×10-4I ,这样可以看成同一值。但如果 是 100kHz ,其差值最大可达到 0.628I 。因此在高频领域,即使距离很短,例如 变电站的母线, 也要考虑波的传播过程, 即当成分布参数线路处理。 变压器有时 也要当作分布参数线路处理。
图 1-10线路上的波过程
1.4.4非线性元件和开关操作
电力系统的暂态过程往往是因状态的变化而造成的。 这种变化可以是断路器 正常或故障操作而引起触头的闭合或开断; 可以是雷电入侵波或操作过电压引起 有间隙避雷器间隙击穿或电流过零时电弧的熄灭; 也可以是系统发生故障造成相
对地或相间突然短接等。 在暂态计算中把电路中节点之间的闭合和开断用广义的 开关操作来表示。 因此, 开关的计算模型以及正确处理开关操作所引起系统状态 变化的程序方法,是电力系统电磁暂态计算的重要组成部分。
电力系统中大部分元件属于线性元件, 或可以近似地认为是线性元件, 但也 有一些元件具有明显的非线性特性, 这些特性对暂态过程产生明显的影响。 典型 的非线性元件有避雷器的非线性电阻,如图 1-11所示;变压器或电抗器等铁磁 元件因铁心饱和而形成的非线性电感; 以及断路器、 保护间隙的电弧电阻等。 因 此, 在暂态计算程序中应包括计及这些非线性元件特性的数学模型, 并且含有一 定的求解非线性电路的数值分析方法。 对工程计算来说, 还希望计算模型和分析 方法尽可能实用,以便在尽可能短的计算时间里,得到具有一定准确度的结果。
10kA 20kA 1000kV避雷器V-I特性
mA 1620kV
图 1-11避雷器的电压 -电流特性
在实际计算中, 经常采用被称为分段线性化的方法来处理非线性元件, 即把 非线性元件的特性用几段具有不同斜率的直线线段来表示, 把非线性元件局部等 值为线性元件。
1.4.5元件参数的频率特性
在电力系统电磁暂态分析过程中,一个元件的特性模拟,不只是要作出正 确的等值电路, 还要模拟它的频率特性, 因为这些频率特性有时对暂态现象有着 决定性的影响。
在暂态计算中, 通常需要考虑频率特性的元件是架空线路和电缆。 架空线路 的正序电感 L 1实际上是常数,在导线的趋肤效应不显著时,正序电阻 1R 基本上 也是常数。零序电感 0L 和零序电阻 0R 则因大地回路的趋肤效应而与频率密切相 关。 图 1-12所示为架空线路电阻和电感的频率特性。 变压器参数也有频率特性, 但通常没有考虑。
图 1-12架空线路电阻和电感的频率特性
1.4.6时间跨度的要求
稳态计算的对象是一个时间断面, 而暂态计算要模拟一个时间过程。 数字计 算机不可能连续地模拟暂态现象, 只能在离散的时间点 (步长 Δt) 求解, 这将会导 致累积误差。如何减少这类误差的积累是暂态仿真程序的重要课题。
鉴于暂态计算的上述特点, 暂态计算比稳态计算不论是程序编制还是应用难 度都要大得多。
1.5电力系统数字仿真
1.5.1电力系统数字仿真的分类
根据原型系统、 数学模型和数字计算机三者的特征可以把电力系统数字仿真 分成各种不同的类型。
按照原型系统状态变化的时间过程, 可分为连续系统仿真和离散事件系统仿 真。连续系统仿真的系统状态量随时间连续变化,它的数学模型是一组方程式, 包括连续模型 (用微分方程描述 ) 、离散时间模型 (用差分方程描述 ) 和连续与离散 混合模型。 离散事件系统仿真的系统状态量只在一些时间点上由某随机事件的驱 动而发生变化, 这类系统在两个事件之间其状态量保持不变, 它的数学模型一般 只用流程图或网络图描述。
按照仿真目的, 可分为以分析研究为目的的研究用系统仿真和以培训运行人 员为目的的培训用系统仿真。 研究用电力系统数字仿真, 如电力系统电磁暂态计 算程序(EMTP ) ,它可用于研究由开关操作、故障和雷击等引起的电磁暂态、 电磁谐振和机电振荡, 也可用于研究交直流换流器、 控制系统和继电保护装置等 的特性。 除此以外, 还有大量适合于专门功能的电力系统数字仿真程序, 如电力 1
系统综合程序(BPA )等。培训用电力系统数字仿真,如电力系统调度员培训仿 真系统(DTS ) 、变电站培训仿真系统等,利用计算机及相关设备,将电力系统 完整的模拟出来, 并可以在上面进行正常操作训练及故障排除训练。 培训用的仿 真是为了训练系统运行、 调度人员对系统环境的反应和判断能力, 因此要求仿真 的环境尽可能逼真,而对于仿真精度,只由培训的要求决定。
研究用电力系统数字仿真又可分为系统稳态计算和暂态计算两大类。 当需要 研究电力系统处于相对平衡状态的运动特性时, 采用系统稳态仿真; 当研究系统 处于受扰动状态的运动特性时, 则采用系统暂态仿真。 两者的数学模型不同, 仿 真方法也不相同。 稳态仿真中有潮流计算、 故障计算, 以及稳定计算和电压稳定 计算中的静稳定计算。暂态仿真中有过电压计算、次同步振荡(SSR )计算、暂 态恢复电压(TRV )计算、高次谐波计算,及稳定计算和电压稳定计算中的暂态 稳定性计算。电力系统数字仿真按照研究电力系统运行状态的分类如图 1-13所 示。
图 1-13电力系统数字仿真按照研究电力系统运行的状态的分类示意图
按照计算方法, 电力系统数字仿真又可以分成有效值计算和瞬时值计算两大 类。 有效值计算用于大规模系统的比较长时间的状态模拟, 而瞬时值计算用于局 部系统的短时间的暂态模拟。 有效值计算通常只使用正序回路, 而瞬时值计算需 要使用三相回路。但是,在包含电力电子设备的情况,如高压直流系统,由于晶 闸管等开关的频繁动作,即使是长时间模拟也需要采用瞬时值计算。另一方面, 由于硬件设备和计算技术的进步,大规模系统的瞬时值在线 (实时 ) 计算也成为可 能了。 稳态计算和暂态计算中的动稳定计算属于有效值计算范畴, 其它的暂态计 算都属于瞬时值范畴。
1.5.2电力系统数字仿真的优点
数字仿真在电力系统研究中得到快速的发展, 除了计算机技术和软件技术的 飞速发展外, 电力系统数字仿真的独特优点是促使其快速发展的重要因素, 这些 优点有:
(1)数字仿真不受被研究电力系统的规模和复杂性的限制。世界各国都在 不断扩大电力系统的规模, 大多数工业发达国家都建立了自己的全国统一电力系 统, 有些相邻国家间还建立了跨国联合电力系统。 我国已实现跨区域电力系统的 互联, 依托三峡工程, 实现以长江三峡为中心的全国统一电力系统。 随着规模的 增大, 电力系统的结构也变得更加复杂。 这些规模庞大和结构复杂的大系统, 试 验和研究的现场实测方法已很难进行。在电力系统动态模拟上做几十台发电机、 几十条输电线路的电力系统暂态过程试验也是相当困难, 而采用电力系统数字仿 真就不存在这些困难, 可以进行数百台发电机和上千条输电线路的大型电力系统 数字仿真。
(2)保证被研究系统的安全性。电力系统的故障试验、稳定性破坏试验、 核电站控制试验等, 直接在原型系统上做实验有很大的危险性, 甚至是不允许的, 这时,用电力系统数字仿真的方法是唯一可行的途径。
(3)提高系统试验的经济性。在实际电力系统上做试验要暂停部分用户供 电,需要配备各种测量设备、测量通道、通信工具,要求很多运行、调度人员和 测试人员密切配合,花费大量人力、物力和财力,因此这种试验很难实现。如果 用数字仿真做试验, 所需费用要少得多。 而且, 数字仿真试验的设备一般都可重 复使用,只需少数计算人员参加,试验时间很短。
(4)增强对电力系统发展的预测性。需要对未来电力系统的特性做预测性 的分析和研究, 这些工作在实际电力系统中难以实现, 而系统数字仿真可以对设 计方案进行大量试验和计算, 进行经济技术比较和优化, 还可以对未来系统的假 设条件的合理性进行验证。 电力系统规划的方案是靠仿真得到的; 新元件的接入、 运行方式的确定是用仿真结果作为依据的; 新方法研究、 新装置设计、 参数确定 是用仿真来确认的。 电力系统仿真软件试验已经成为电力系统设计、 规划和运行 阶段不可或缺的部分。
1.5.3电力系统数字仿真软件
世界各国的电力系统数字仿真软件众多,目前国内外获得广泛应用的电力 系统仿真软件主要有 3大类:
基于瞬时值计算的离线仿真软件,如 EMTP 、 PSCAD/EMTDC、 PSAPAC 等;
基于有效值计算的离线仿真软件,如 BPA (包括中国版 BPA ) 、 PSASP 、 PSS/E、 NETOMAC 等;
基于瞬时值计算的实时仿真软件,如 RTDS 、 HYPERSIM 。
1. 电磁暂态分析程序(EMTP )
电磁暂态分析程序 EMTP (Electro-Magnetic Transients Program) 是国际公认 的电力系统电磁暂态分析的标准程序,其创始人是加拿大 UBC 大学任教的 Dommel 教授, 目前 EMTP 有三个版本, 即 BPA-EMTP 、 ATP-EMTP 和 DCG-EMTP 。 BPA-EMTP 是最早由美国邦纳维尔电力局 (Bonneville Power Adminstration, BPA ) 无偿提供的 EMTP 版本, 其用户现在大多已转用 ATP-EMTP 。 ATP-EMTP 是 BPA 的 Scott-Meyer 以自己的业余时间和资金开发的 BPA-EMTP 的替代程序, ATP-EMTP 坚持无偿提供的原则, 在全世界拥有最多的用户, 是目前国际上主流 版本的 EMTP 程序。 DCG-EMTP 是 1981年成立的 DCG (EMTP 合作开发组织) 开发的 EMTP 版本,需有偿使用。
EMTP 具有分析功能多和元件模型全等优点, 可以用于电力系统的稳态和暂 态仿真分析,系统可由集中参数、分布参数元件、线性与非线性元件、具有频率 相关参数的线路、 各种类型开关、 电力电子元件、 变压器及电机、 多种类型电源、 控制电路的任意组合构成,只要是电路计算的范畴,对研究对象几乎没有限制。 EMTP 的计算精度经过了 IEEE 和 CIGRE 等国际权威组织的认定,因此计算结 果的可信性很高。
实际上, EMTP 是开发其他的电磁暂态程序, 如 EMTDC 、 RTDS 、 HYPERSIM 的基础。
2.直流电磁暂态程序(PSCAD/EMTDC)
Dennis Woodford 博士于 1976年在加拿大曼尼托巴直流研究中心开发完成 了 EMTDC (Electro-Magnetic Transients in DC System)的初版,它既可以研究交 直流电力系统问题, 又能完成电力电子仿真及其非线性控制, 是一个离线仿真的
电磁暂态计算程序, 它有精确的直流元件模型、 方便的数据输入方式以及强大的 数据分析功能,是进行直流系统分析和工程研究的有力工具。事实上, EMTP 程 序所采用的电力系统模型和技术都可以应用于 EMTDC 中。 PSCAD (Power System Computer Aided Design)是其图形用户界面, PSCAD 的开发成功,使得 用户能更方便地使用 EMTDC 进行电力系统分析, 使电力系统复杂部分可视化成 为可能, 而且软件可以作为实时数字仿真器的前置端。 可模拟任意大小的交直流 系统,在对直流系统电磁暂态仿真方面有绝对优势。操作环境为:UNIX OS, Windows95, 98,NT; Fortran 编辑器;浏览器和 TCP/IP协议。主要功能如下:
(1)可以分析系统中断路器操作、故障及雷击时出现的过电压;
(2)可对包含复杂非线性元件(如直流输电设备)的大型电力系统进行全 三相的精确模拟,其输入、输出界面非常直观、方便;
(3)进行电力系统时域或频域计算仿真;
(4)电力系统谐波分析及电力电子领域的仿真计算;
(5)实现高压直流输电、 FACTS 控制器的设计。
3. PSAPAC 程序
PSAPAC 由美国 EPRI 开发,是一个全面分析电力系统静态和动态性能的软 件工具。功能:
DYNRED (Dynamic Reduction Program) :网络化简与系统的动态等值,保 留需要的节点。
LOADSYN (Load Synthesis Program) :模拟静态负荷模型和动态负荷模型。 IPFLOW (Interactive Power Flow Program) :采用快速分解法和牛顿-拉夫 逊法相结合的潮流分析方法, 由电压稳态分析工具和不同负荷、 事故及发电调度 的潮流条件构成。
TLIM (Transfer Limit Program) :快速计算电力潮流和各种负荷、事故及发 电调度的输电线的传输极限。
DIRECT :直接法稳定分析软件弥补了传统时域仿真工作量大、 费时的缺陷, 并且提供了计算稳定裕度的方法,增强了时域仿真的能力。
LTSP (Long Term Stability Program) :LTSP 是时域仿真程序, 用来模拟大型 电力系统受到扰动后的长期动态过程。 为了保证仿真的精确性, 提供了详细的模
型和方法。
VSTAB (V oltage Stability Program) :该程序用来评价大型复杂电力系统的电 压稳定性, 给出接近于电压不稳定的信息和不稳定机理。 为了估计电压不稳定状 态,使用了一种增强的潮流程序,提供了一种接近不稳定的模式分析方法。 ETMSP (Extended Transient Midterm Stability Program) :EPRI 为分析大型 电力系统暂态和中期稳定性而开发的一种时域仿真程序。 为了满足大型电力系统 的仿真, 程序采用了稀疏技术, 解网络方程时为得到最合适的排序采用了网络拓 扑关系并采用了显式积分和隐式积分等数值积分法。
SSSP(Small-signal Stability Program):该程序有助于局部电厂模式振荡和站 间模式振荡的分析,由多区域小信号稳定程序 (MASS)及大型系统特征值分析程 序(PEALS )两个子程序组成。 MASS 程序采用了 QR 变换法计算矩阵的所有特 征值, 由于系统的所有模式都计算, 它对控制的设计和协调是理想的工具; PEALS 使用了两种技术:AESOPS 算法和改进 Arnoldi 方法,这两种算法高效、可靠, 而且在满足大型复杂电力系统的小信号稳定性分析的要求上互为补充。
4. BPA 程序
BPA 程序是美国联邦政府能源部下属邦纳维尔电力局计算方法开发组自 20世纪 60年代初期开发的大型电力系统离线分析程序。该程序采用稀疏矩阵技巧 的牛顿—拉夫逊法, 并将梯形积分法运用于暂态稳定的计算, 形成较为稳定的数 值解。 BPA 程序的结构分为潮流程序和稳定程序两部分。
BPA 潮流程序主要用来计算电力系统潮流。该程序中的负荷模型包含恒定 功率负荷、恒定电流负荷和恒定阻抗负荷模型。可以根据某节点上 P 和 Q 的扰 动量, 计算系统中各节点灵敏度、 线路灵敏度和网损灵敏度值。 程序的输出具有 内容详细和格式灵活的特点, 既可以有选择地列表输出原始数据、 计算结果和潮 流分析报告, 也可以应用单线图格式潮流图形程序及地理接线图格式潮流图形程 序输出。
BPA 稳定程序含有 9种传统励磁模型和 1981年 IEEE 提出的 11种新励磁模 型, 可模拟多种类型的直流型励磁机、 交流型励磁机及静态型励磁机, 可以进行 多端直流的模拟。 程序可以在屏幕上输出最大摇摆角, 还可以给出对应的两台发 电机名。
BPA 程序的主要功能是进行大型交直流混合电力系统潮流和暂态稳定计算, 同时还能进行短路电流计算和电网静态等值分析等。 BPA 配备有较完善的辅助工 具, 包括单线图和地理接线图格式潮流图程序、 稳定曲线作图工具。 具有计算规 模大、计算速度快、数据稳定性好、 功能强等特点。 BPA 稳定程序中包括详细的 发电机和各种励磁系统模型, 但没有提供用户自定义功能, 程序中 HVDC 、 FACTS 元件及其控制模型不够完善。 另外, BPA 的数据格式要求比较严格, 与一些国际 上通用的机电暂态仿真程序之间的数据互换比较困难。
中国版 BPA 程序是由中国电力科学院引进、 消化、 吸收美国 BPA 程序 1983年 9月版本的基础上开发而成。从 1984年开始在我国推广应用以来,已在我国 电力系统规划设计、 调度运行和试验研究等各部门得到了广泛的应用, 成为我国 电力系统分析计算的重要工具之一。
5.电力系统分析综合程序(PSASP )
电力系统分析综合程序 (Power System Analysis Software Package, PSASP) 是 中国电力科学研究院开发的一套具有高度集成性、 开放性的大型软件包。 PSASP 与 Execl 、 AutoCAD 、 MATLAB 等通用的软件包分析工具有方便的接口, 可充分 利用其它软件包的资源。 该软件在我国高校研究人员和电力系统现场都有广泛应 用。
PSASP 结构分为三层,第一层是公用数据和模型的资源库,其中包括 :电网 基础库、 固定模型库、 用户自定义模型库和用户程序库等。 第二层是基于资源库 的应用程序包,包括稳态分析、故障分析、机电暂态分析和暂态稳定计算。第三 层是计算结果库和分析工具, 软件进行各种分析计算后, 生成的结果数据以多种 形式输出或转换为 Excel 、 AutoCAD 、 MATLAB 等其他数据格式。
PSASP 的功能主要有稳态分析、故障分析和机电暂态分析。稳态分析包括 潮流分析、网损分析、最优潮流和无功优化、静态安全分析、谐波分析和静态等 值等。 故障分析包括短路计算、 复杂故障计算及继电保护整定计算。 机电暂态分 析包括暂态稳定计算、电压稳定计算、控制参数优化等。
6. PSS/E软件
电力系统仿真软件 PSS/E(Power System Simulator for Engineering)是由美 国电力技术公司(PTI )开发的商业软件,主要用于电力系统仿真和计算。 PSS/E
程序国际影响较大, 在国内 PSS/E以其独有的魅力正受到各高校、 研究机构的青 睐。
PSS/E包含了电力系统机电暂态分析计算的常见模块,其主要特点有:仿真 规模非常大, 有利用超大规模系统的计算; 灵活的模型自定义功能; 强劲的交互 式计算过程控制等。 另外, PSS/E数据格式标准能与多种仿真程序进行数据转换。 在版本更新方面, PSS/E做的更好,随着版本的不断更新, PSS/E模型库中的各 种元件模型更加完善。
7. NETOMAC程序
NETOMAC (Network Torsion Machine Control)是德国西门子公司研制的电 力系统仿真软件。 NETOMAC 仿真分析软件功能强大,较好的结合了机电暂态 和电磁暂态两大应用领域, 也具有网络等值、 频域分析、 参数辨识和优化等功能。 但 NETOMAC 程序的使用存在掌握较为困难,潮流计算收敛性较差,无法方便 绘制地理接线图格式的潮流图,输入数据与 BPA 转换较为困难,用户数据输入 /输出界面不够友好等不足。在世界各国尤其是欧洲拥有众多用户。
NETOMAC 软件程序结构和特点:
(1)程序元件模型健全 , 该软件模型库几乎包含了当前电力系统所有的元 件。
(2)仿真频带宽,该软件既可模拟雷电波的过程,又能进行电磁暂态、机 电暂态、稳态等电力系统过程的仿真计算。
(3)功能多且强,该软件可进行潮流、短路、稳定、动态等值、电动机启 动、参数辨识、机组轴系扭振、优化潮流等多种计算。
(4) NETOMAC 的微机版本可以计算包含数千条线路和数百台发电机的大 型网络,能较好地满足大规模实际电力系统仿真的需要。
8. 实时数字仿真器 (RTDS )
RTDS (Real Time Digital Simulators)是加拿大曼尼托巴 RTDS 公司开发的一种实 时全数字电磁暂态电力系统模拟装置, 它是电力系统规模的不断扩大和计算机并 行处理技术及数字仿真技术飞速发展的产物。其算法与 EMTP 仿真程序相同, RTDS 的图形用户界面是 RSCAD 。 RTDS 提供了几乎所有的传统电力系统元件模 型,并且可实现物理和数字混合仿真,对控制设备和保护装置进行测试。因此,
RTDS 仿真系统可应用于电力系统的规划及优化、可行性研究、保护系统的设计 与测试、系统稳态及电磁暂态分析及其教育培训等诸多方面。 RTDS 的数据文件 比较复杂, 目前 RSCAD 自带了 BPA-RSCAD 和 PSS/E-RSCAD的数据转换程序, 但转换程序均不完善。
9.全数字电力系统实时仿真系统(HYPERSIM )
HYPERSIM 是加拿大魁北克 TET 公司开发的一种基于并行计算技术、采用 模块化设计、面向对象编程的电力系统全数字仿真系统。 HYPERSIM 目前具有 Unix , Linux 、 Windows 等 3种版本, 即可在 Sun Unix工作站或 Linux/Windows PC机上进行离线仿真,也可在 SGI 超级计算机或者 Linux PC-Cluster上与实际电力 系统器件连接而进行实时在线仿真。 HYPERSIM 提供了电磁仿真的灵活性,比 传统的模拟仿真器更加灵活,简单,廉价,可用于电力系统、电力电子及电机拖 动分析,用于 HVDC 及 FACTS 设备动态性能测试,以及控制系统性能测试、继 电保护和重合闸装置闭环测试。
HYPERSIM 硬件采用基于共享存储器的多 CPU 超级并行处理计算机如 SGI3000(SGI350)或多 CPU 的并行计算用的 Alpha 工作站。主要用于电力系 统电磁暂态仿真,仿真的规模可以相当大,也可以用于装置试验。其中基于 SGI3200服务器的 HYPERSIM 也可用于直流系统动态特性仿真, HYPERSIM 软 件的核心是 EMTP 程序。
习题
1-1 什么是电力系统电磁暂态现象?
1-2 分析电力系统电磁暂态现象的目的是什么?
1-3 电力系统电磁暂态研究的方法有哪些?各有什么特点?
1-4 电力系统电磁暂态现象的特点有哪些?
1-5 电力系统数值计算中著名的仿真程序有哪些?
范文五:电力系统暂态分析讲义
第八章 电力系统暂态稳定
第一节 暂态稳定概述
暂态稳定 荷突增) 、甩负荷(负荷突减)等。
能保持暂态稳定:扰动后,系统能达到稳态运行。 分析暂态稳定的时间段:
起始:0~1s ,保护、自动装置动作,但调节系统作用不明显,发电机采用
q
E '、 P T 恒定模型; 中间:1~5s, AVR 、 P T 的变化明显,须计及励磁、调速系统各环节; 后期:5s~mins ,各种设备的影响显著,描述系统的方程多。 本书中重点讨论 起始 阶段。
基本假定:⑴ 网络中, ω=ω0 (网络等值电路同稳态分析) ⑵ 只计及正序基波分量, 短路故障用正序增广网络表示
一.物理过程分析
发电机采用 E ’ 模型。 故障前:
2
21T L
T d
I x x x x x +++'= 电源 电势节点 到系统的 直接 电抗 δsin I
I x U
E P '=
故障中,
?
++'++++'=x x
x x x x x x x T L
T d
T L
T d
II ) )(() 2
() (2122
δsin II
II x U
E P '=
21T L T d
III x x x x x +++'= δsin III
III x U
E P '=
功角特性曲线为:
● 故障发生后的过程为:
运行点变化 原因 结果
a → b 短路发生 PT >PE , 加速, ω上升, δ增大
b→ c ω上升, δ增大 ω>ω0 ,动能增加
c→ e 故障切除 PT e → f 动能释放 减速,当 ωf =ω0,动能释放完毕, δm 角达最大 f → k PT ① 若最大摇摆角 h m δδ<,系统可经衰减的振荡后停止于稳定平衡点 k="" ,系统="">,系统可经衰减的振荡后停止于稳定平衡点> ② 暂态稳定分析与初始运行方式、故障点条件、故障切除时间、故障后状态 有关。 ③ 电力系统暂态稳定分析是计算电力系统故障及恢复期间内各发电机组的 功率角 i δ的变化情况(即 δ–t曲线) ,然后根据 i δ角有无趋向恒定(稳定) 数值,来判断系统能否保持稳定,求解方法是非线性微分方程的数值求解。 二.等面积定则 ●故障中,机组输入的机械功率 T P >发 电机输出的电磁功率 E P ,发电机加速, II T J P P dt d T -=2 20δω δδδδδδδδδd d d d dt d dt d dt d dt d dt d ====) (22 δδδωd P P d T E T J ) (0 -=∴ 积分得: ?-==-c d P P T T II T c J c J δδδδωδδω0 ) (21) (212 2020 左侧=转子在相对运动中动能的增量; 右侧=过剩转矩对相对位移所做的功―― T P 线下方的阴影面积――称为加速 面积; ●故障切除后 III T J P P dt d T -=2 2 0δω () ?-=-m d P P T III T m J δδδδδω0 ) (212 020 ∵ m δδ=时, 0ωω=f , ∴ 0=m δ ?-=m c d P P T T III c J δδδδω) (210 右侧=制动转矩对相对角位移所做的功 =T P 线上方的阴影面积(称为减速面积) ● 因减速过程中,转速恢复同步转速(即加速过程中的动能释放完毕) 时 δ角达最大,所以加速面积=减速面积――等面积定则。 等面积定则的应用 ① 在三状态暂态稳定分析中确定极限切除角 δ角摇摆越过 δh ,则 P T >PE ,将使 δ角 继续扩大,系统失去暂态稳定。 ∴ δh 是最大允许摇摆角。 当加速面积与允许的减速面积相等时, ??-=-cm h cm d P P d P P T III II T δδδδδδ0 ) () ( ??-=-cm h cm d P P d P P T IIIM IIM T δδδδδδδδ0 ) sin () sin ( ∴ IIM IIIM IIM h IIIM h T cm P P P P P --+-= 9cos cos ) (cos δδδδδ 暂态稳定判据 1:cm c δδ≤,系统能保持暂态稳定,否则不能保持暂态稳定 。 ②加速面积、减速面积的应用例――单相重合闸的作用分析 重合闸动作分析: 暂态稳定判据 2:实际加速面积 <允许的减速面积>允许的减速面积> 三.微分方程的数值解法(摇摆曲线法) 故障中,转子运动方程为: ) sin (1) 1(0δωωωδ II T J x U E P T dt d dt d '-=-= 初始条件:IM T P P t 10sin , 1 , 0-====δδω 故障切除后, ) sin (1) 1(0δωωωδ III T J x U E P T dt d dt d '-=-= 初始条件:c c c t t δδωω=== , , 是已知初始条件,可进行数值求解的非线性 简单系统,当 δ达 δ m 后开始下降,说明功角特性曲线上运行点开始往平衡点 k 移动,所以 暂态稳定判据 3:简单系统中,当 δ达 δ m 后开始减小,则系统能保持暂态稳 定; δ>180°,系统不能保持暂态稳定。 δ-t 曲线的计算方法(微分方程数值解法) (二)改进欧拉法 四 提高暂态稳定的措施 ● 缩短电气距离(即提高静态稳定的措施) ● 减少功率差额 : 1.保护装置 ① 快速切除故障,可减小加速面积; ② 采用自动重合闸,重合成功可增大减速面积; 2.提高发电机输出的电磁功率 ① 强励,提高端电压,增大有功输出; ② 电气制动,消耗发电机有功功率; 3. 快速关闭汽门 4.系统失去稳定后的措施 ⑴ 适当设置解列点; ⑵ 异步运行再同步。 考核要求: 1、电力系统暂态稳定的物理过程。 2、等面积定则分析电力系统暂态稳定。 3、提高暂态稳定的措施。 例题: 习题 一、简答题 1. 写出五条提高电力系统暂态稳定运行的主要措施。 1) 采用自动重合闸装置。 2) 发电机装设强行励磁装置。 3) 电气制动。 4) 变压器中性点经小电阻接地。 5)快速切除故障。 12. 什么是电力系统暂态稳定 ? 电力系统的瞬态稳定性, 是指电力系统在某一运行状态下受到较大的干扰后, 能够过渡到一个新的运行 状态或恢复到原来的运行状态的能力。 13. 试述等面积定则的基本含义。 在加速期间积蓄的动能增量全部耗尽,即加速面积和减速面积大小相等,这就是等面积定则。 14. 解释采用减少原动机输出功率的方法能够改善系统的暂态稳定性。 减少了加速面积,增大了减速面积。 15.试用等面积定则分析当自动重合闸成功时为什么可以提高电力系统暂态稳定性? 增大了减速面积。 二、计算题 1. 设已知系统短路前、短路时、短路切除后三种情况的以标幺值表示的功角特性曲线:P 1=2sin δ、 P 2=0.5sin δ、 P 3=1.5sin δ及输入发电机的机械功率 P T =1。求极限切除角。 (6分 ) 设已知系统短路前、 短路时、 短路切除后三种情况的以标幺值表示的功角特性曲线: δsin 21=P 、 δsin 5. 02=P 、 δsin 5. 13=P 及输入发电机的机械功率 1=T P 。求极限切除角。 由 T P P =1得:rad 5236 . 0300=?=δ 由 T P P =3得:rad h 4112. 219. 138810. 41180=?=?-?=δ 由 3366. 0cos cos ) (cos 230 230=--+-= m m m h m h T cm P P P P P δδδδδ 所以, ?=33. 70cm δ 3. 简单系统及参数标么值如图,若在 K 点发生三相短路,故障后经一段时间切除故障线路,求使得系 统保持暂态稳定的极限切除角。 ?0X d ′ =j0.2 转载请注明出处范文大全网 » 电力系统电磁暂态分析