电力系统课程设计附分段法程序模板.doc

电力系统课程设计附分段法程序

47

资料内容仅供参考，如有不当或者侵权，请联系本人改正或者删除。

电力系统分析
课程设计

题 目: 潮流计算与故障分析
院系名称: 电气工程学院 专业班级:
学生姓名:
学生姓名:
指导教师:
指导教师意见:

成绩:
签名: 年 月 日
年 1 月 15 日

目录
1绪论 1
1.1背景 1
1.2国内外发展现状 1
1.3 课题目标 1
1.4设计要求 1
2 方案论证 2
2.1电力系统潮流计算的原理 2
2.2 电力系统故障分析的原理 3
3潮流计算与分析 4
3.1设计指标 4
3.2潮流计算 4
3.3计算结果分析 8
3.4暂态稳定定性分析 9
3.5暂态稳定定量分析 12
3.6运行结果与分析 17
4心得总结 29
5主要参考文献 30
附录 31

1绪论
1.1背景
电力系统潮流计算是研究电力系统稳态运行情况的一种基本电气计算。它的任务是根据给定的运行条件和网络结构确定整个系统的运行状态, 如各母线上的电压(幅值及相角)、 网络中的功率分布以及功率损耗等。电力系统潮流计算的结果是电力系统稳定计算和故障分析的基础。
1.2国内外发展现状
近20多年来, 潮流算法的研究依然非常活跃, 可是大多数研究都是围绕改进牛顿法和P-Q分解法进行的。另外, 随着人工智能理论的发展, 遗传算法、 人工神经网络、 模糊算法也逐渐被引入潮流计算。可是, 到当前为止这些新的模型和算法还不能取代牛顿法和P-Q分解法的地位。由于电力系统规模的不断扩大, 对计算速度的要求不断提高, 计算机的并行计算技术也将在潮流计算中得到广泛的应用, 成为重要的研究领域。
1.3 课题目标
学会使用电力系统分析软件。经过电力系统分析软件对电力系统的运行进行实例分析, 加深和巩固课堂教学内容。
根据所给的电力系统, 绘制短路电流计算程序, 经过计算机进行调试, 最后成一个切实可行的电力系统计算应用程序, 经过自己设计电力系统计算程序不但能够加深学生对短路计算的理解, 还能够锻炼学生的计算机实际应用能力。
熟悉电力系统分析综合这门课程, 复习电力系统潮流计算和故障分析的方法。了解Simulink在进行潮流、 故障分析时电力系统各元件所用的不同的数学模型并在进行不同的计算时加以正确选用。学会用Simulink, 经过图形编辑建模, 并对特定网络进行计算分析。
1.4设计要求
系统的暂态稳定性是系统受到大干扰后如短路等, 系统能否恢复到同步运行状态。图1为一单机无穷大系统, 分析在f点发生短路故障, 经过线路两侧开关同时断开切除线路后, 分析系统的暂态稳定性。若切除及时, 则发电机的功角保持稳定, 转速也将趋于稳定。若故障切除晚, 则转速曲线发散。

图1 单机无穷大系统
发电机的参数: SGN=352.5MWA,PGN=300MW,UGN=10.5Kv,,,,,,,,,Rs=0.0028,H(s)=4s;TJN=8s,负序电抗: 。
变压器T-1的参数:STN1=360MVA,UST1%=14%,KT1=10.5/242;
变压器T-2的参数:STN2=360MVA,UST2%=14%,KT2=220/121;
线路的参数: l=279km,UN=220K,XL=0.41欧/km, 欧/km, 线路的零序电抗为正序电抗的5倍。运行条件如下: =115KV, =250W, 功率因数为0.95。 当0.029s发生单相故障时, 断路器0.29s动作切除故障, 判定发电机转速是否达到稳定? 系统是否稳定?
2 方案论证
进行电力系统的潮流分析有两种方案进行选择: 牛顿—拉夫逊法; P—Q法。我们采用最基本的牛顿—拉夫逊法, 这种方法简单快捷, 容易掌握。
2.1电力系统潮流计算的原理
电力系统潮流计算是电力系统分析中的一种最基本的计算, 是对复杂电力系统正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运行状态的计算, 即节点电压和功率分布, 用以检查系统各元件是否过负荷。各点电压是否满足要求, 功率的分布和分配是否合理以及功率损耗等。对现有的电力系统的运行和扩建, 对新的电力系统进行规划设计以及对电力系统进行静态和稳态分析都是以潮流计算为基础。潮流计算结果可用如电力系统稳态研究, 安全估计或最优潮流等对潮流计算的模型和方法有直接影响。实际电力系统的潮流技术那主要采用牛顿—拉夫逊法。
潮流计算是电力系统分析最基本的计算。除它自身的重要作用之外, 潮流计算还是网损计算、 静态安全分析、 暂态稳定计算、 小干扰静态稳定计算、 短路计算、 静态和动态等值计算的基础。电力系统潮流计算也分为离线计算和在线计算两种, 前者主要用于系统规划设计和安排系统的运行方式, 后者则用于正在运行系统的经常监视及实时控制。
2.2 电力系统故障分析的原理
电力系统稳定性是指当系统在某一正常运行状态下受到某种扰动后, 能否经过一定的时间后恢复到原来的运行状态或者过渡到一个新的稳定运行状态。如果能够回到原来的运行状态或者建立一个新的稳定运行状态, 则认为系统在该运行状态下是稳定的。反之, 若系统不能够回到原来运行状态或者不能建立一个新的稳定运行状态, 则说明系统的状态变量没有一个稳态值, 而是随时间不断增大或振荡, 系统是不稳定的。电力系统稳定性被破坏后, 将造成大量用户供电中断, 甚至导致整个系统的瓦解, 后果极为严重。因此, 保持电力系统的稳定性, 对于电力系统安全可靠运行, 具有非常重要的意义。
在电力系统运行过程中, 时常会发生故障, 其中大多数是短路故障( 简称短路) 。所谓短路, 是指电力系统正常运行情况以外的相与相之间或相与地( 或中性线) 之间的连接。发生短路时, 由于系统中总阻抗大大减少, 因而短路电流可能达到很大数值( 几万安至十几万安) 。这样大的电流所产生的热效应和机械效应会使电气设备受到破坏; 同时短路点的电压降到零, 短路点附近的电压也相应地显著降低, 使此处的供电系统受到严重影响或被迫中断; 若在发电厂附近发生短路, 还可能使整个电力系统运行解列, 引起严重后果。

3潮流计算与分析
3.1设计指标
3.1.1网络参数及运行参数计算
取, 。为使变压器不出现非标准变比, 各段基准电压为,
3.1.2各元件参数归算后的标么值如下:
, , , , ,

3.1.3 运算参数的计算结果如下:
; ;

3.2潮流计算
上述单机无穷大系统在发电机保持不变的条件下的潮流计算如下:

①当保持时

取正号得

②当保持时

③当保持时

④当保持时

3.3计算结果分析
由上述计算结果可作出图2功率特性图和图3电势变化特性图。
计算结果说明和的差别是很小的。同时, 和在各种计算条件下的差别也是很小的。从的变化能够看到, 保持常数的电势不同, 调节器的放大系数也不同。

图2 功率特性计算结果

( a) ( b) (c)
图3 各电势变化特性

3.4暂态稳定定性分析
如图4所示为一正常运行是的简单电力系统及其等值电路, 发电机经过变压器和双回路向无限大系统送电。发电机在正常运行、 故障以及故障切除后3种状态下的功角特性曲线如图5所示

图4 简单电力系统及其等值电路
a)正常运行方式及其等值电 b) 故障情况及其等值电路 c)故障切除后及其等值电路

图5 简单系统正常运行、 故障及故障切除后的功率特性
①正常运行
正常运行时发电机的功率特性曲线为, 此时向无穷大系统输送的功率与原动机输出的机械功率相等( 假设扰动后保持不变) , 图5中的a点即正常运行发电机的运行点, 此时功角为。
②故障阶段
发生短路后功率特性立即降为, 但由于转子的惯性, 转子角度不会立即变化, 发电机的运行点由a点突变至b点, 输出功率显著减少, 而原动机机械功率不变, 故产生较大的过剩功率。故障情况愈严重, 功率曲线幅值愈低( 三相短路时为零) , 则过剩功率愈大。在过剩转矩的作用下发电机转子将加速, 其相对速度( 相对于同步转速) 和相对角度逐渐增大, 使运行点由b点向c点移动。如果故障一直存在, 则始终存在过剩转矩, 发电机将不断加速, 最终与无限大系统失去同步。
③故障及时切除
实际上, 短路故障后继电保护装置将迅速动作切除故障路线。假设在c点时将故障切除, 则发电机的功率特性变为, 发电机的运行点从c点突然变至e点( 同样由于不能突变) 。这时, 发电机的输出功率比原动机的机械功率大, 使转子受到制动, 转子速度逐渐减慢。但由于此时的速度已经大于同步转速, 因此相对角度还要继续增大。假设制动过程延续到f点时转子转速才回到同步转速, 则角不再增大。可是, 在f点是不能持续运行的, 因为这时机械功率和电磁功率仍不平衡, 前者小于后者。转子将继续减速, 开始减小, 运行点沿功率特性由f点向e、 k点转移。在达到k点以前转子一直减速, 转子速度低于同步转速。在k点虽然机械功率与电磁功率平衡, 但由于这时转子速度低于同步转速, 继续减小。但越过k点以后机械功率开始大于电磁功率, 转子又加速, 因而一直减小到转速恢复同步转速后又开始增大。此后运行点沿着开始第二次振荡。如果振荡过程汇总没有任何能量损耗, 则第二次又将增大至f点的对应角度, 以后就一直沿着往复不已的振荡。实际上, 振荡过程中总有能量损耗, 或者说总存在着阻尼作用, 因而振荡逐渐衰弱, 发电机最后停留在一个新的运行点k上持续运行。K点即故障切除后功率特性与的交点。图6画出了上述振荡过程中负的过剩功率, 转子角速度和相对角度随时间变化的情形。图中考虑了阻尼作用。

图6 振荡过程
④故障切除过晚
如果故障线路切除得过晚, 如图7所示。这时在故障线路切除前转子加速已比较严重, 因此当故障线路切除后, 在到达与图5中相应的f点时转子转速仍大于同步转速。甚至在到达h点时转速还未降至同步转速, 因此就将越过h点对应的角度。而当运行点越过h点后, 转子又立即承受加速转矩, 转速又开始升高, 而且加速度越来越大, 将不断增大, 发电机和无限大系统之间最终失去同步。失步过程如图8所示。

图7 故障切除过晚的情形 图8 失步过程
由上可见, 快速切除故障是保证暂态稳定的有效措施。
由以上分析可知, 系统暂态稳定与否是和正常运行的情况( 决定和大小) 以及扰动情况( 发生什么故障、 何时切除) 紧密相关。为了准确判断系统在某个运行方式下受到某种扰动后能否保持暂态稳定, 必须经过定量的分析计算。
3.5暂态稳定定量分析
⑴系统转移阻抗和功率特性计算, 在图1单机无穷大系统中, 当f点发生两相短路时的负序和零序等值网络如图9 a) b)所示

图9 序网及短路时的等值电路图
a)负序网络 b)零序网络 c)短路是的等值电路

两相接地时的短路附加电抗为

短路时的等值电路如图c)所示, 系统的转移电抗和功率特性分别为

故障切除后系统的转移电抗和功率特性分别为

系统极限切除角计算
应用等面积定则, 可求的极限切除角为

式中, 临界角
以上的角度都是用弧度表示的, 换算成度数为

根据分段计算法计算极限切除时间
发电机惯性时间常数

取为0.05s,

第一个时间段

第二个时间段

第三个时间段开始瞬间, 故障被切除, 故

在上述简单系统中, 短路故障期间发电机摇摆曲线即转子的运动方程为
⑴
这是两个一阶的非线性常微分方程, 它们的起始条件是已知的, 即
; ;
当计算出故障期间的曲线后, 就可由曲线找到与极限切除角相应的极限切除时间。
如果问题是已知切除时间, 而需要求出曲线来判断系统的稳定性, 则当曲线计算到故障切除时, 出于系统参数改变, 以致发电机功率特性发生变化, 必须开始求解另一组微分方程, 即
⑵
式⑵的起始条件为 ; ;
其中, 为给定的切除时间, 、 为与时刻对应的和, 可由故障期间的曲线和曲线求得( 和都是不能突变的) 。这样, 由式⑵可继续求得和随时间变化的曲线。一般讲, 在计算几秒钟内的变化过程时, 如果始终不超过, 而且振荡幅值越来越小, 则系统是暂态稳定的。
求解式⑴的微分方程组, 得到系统故障期间的曲线如图10所示。从图中可查得对应的极限切除角的极限切除时间为0.2416s。

图10 系统故障期间的曲线
求解式⑵的微分方程组, 得到故障切除后系统的曲线如图11所示。从图中能够看出, 到0.52s时即开始减少, 最大角度为, 系统是稳定的。

图11 故障切除后系统的曲线
⑵在图1单机无穷大系统中, 当f点在0.014s发生单相短路时

单相短路时的短路附加电抗为

短路时系统的转移电抗和功率特性分别为

故障切除后系统的转移电抗和功率特性分别为

由于, 故无法确定极限切除角。而在故障不切除时, 其临界角为

对应的
不切除故障时的面积

即可能的减速面积大于可能的加速面积, 系统不切除故障也能保持暂态稳定, 故无法也没有必要确定极限切除角。因此断路器在0.14s时故障切除与否, 发电机转速都能达到稳定, 系统都能保持暂态稳定。
3.6运行结果与分析
>> Untitled6
程序计算结果如下:
输出不平等功率:DP=
Columns 1 through 10
-0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000

Columns 11 through 20
-0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000

Columns 21 through 30
-0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 0.9746 0.9746

Columns 31 through 40
0.9746 0.9745 0.9745 0.9745 0.9744 0.9743 0.9743 0.9742 0.9741 0.9740

Columns 41 through 50
0.9739 0.9737 0.9736 0.9735 0.9733 0.9732 0.9730 0.9728 0.9726 0.9724

Columns 51 through 60
0.9722 0.9720 0.9717 0.9715 0.9712 0.9710 0.9707 0.9704 0.9701 0.9698

Columns 61 through 70
0.9695 0.9692 0.9689 0.9686 0.9682 0.9678 0.9675 0.9671 0.9667 0.9663

Columns 71 through 80
0.9659 0.9655 0.9651 0.9646 0.9642 0.9637 0.9633 0.9628 0.9623 0.9618

Columns 81 through 90
0.9613 0.9608 0.9603 0.9597 0.9592 0.9586 0.9581 0.9575 0.9569 0.9563
Columns 91 through 100
0.9557 0.9551 0.9545 0.9539 0.9533 0.9526 0.9519 0.9513 0.9506 0.9499

Columns 101 through 110
0.9492 0.9485 0.9478 0.9471 0.9464 0.9456 0.9449 0.9441 0.9434 0.9426

Columns 111 through 120

0.9418 0.9410 0.9402 0.9394 0.9386 0.9377 0.9369 0.9360 0.9352 0.9343

Columns 121 through 130
0.9334 0.9325 0.9317 0.9307 0.9298 0.9289 0.9280 0.9270 0.9261 0.9251

Columns 131 through 140
0.9242 0.9232 0.9222 0.9212 0.9202 0.9192 0.9182 0.9172 0.9162 0.9151

Columns 141 through 150
0.9141 0.9130 0.9120 0.9109 0.9098 0.9087 0.9076 0.9065 0.9054 0.9043

Columns 151 through 160
0.9031 0.9020 0.9009 0.8997 0.8986 0.8974 0.8962 08950 0.8938 0.8926

Columns 161 through 170
0.8914 0.8902 0.8890 0.8878 0.8865 0.8853 0.8841 0.8828 0.8815 0.8803

Columns 171 through 180
0.8790 0.8777 0.8764 0.8751 0.8738 0.8725 0.8712 0.8699 0.8686 0.8672

Columns 181 through 190
0.8659 0.8645 0.8632 0.8618 0.8605 0.8591 0.8577 0.8563 0.8550 0.8536

Columns 191 through 200
0.8522 0.8508 0.8494 0.8480 0.8465 0.8451 0.8437 0.8423 0.8408 0.8394

Columns 201 through 210
0.8379 0.8365 0.8350 0.8336 0.8321 0.8306 0.8292 0.8277 0.8262 0.8247

Columns 211 through 220
0.8233 0.8218 0.8203 0.8188 0.8173 0.8158 0.8143 0.8128 0.8113 0.8097

Columns 221 through 230
0.8082 0.8067 0.8052 0.8037 0.8021 0.8006 0.7991 0.7976 0.7960 0.7945

Columns 231 through 240
0.7930 0.7914 0.7899 0.7884 0.7868 0.7853 0.7838 0.7822 0.7807 0.7791

Columns 241 through 250
0.7776 0.7761 0.7745 0.7730 0.7714 0.7699 0.7684 0.7668 0.7653 0.7638

Columns 251 through 260
0.7622 0.7607 0.7592 0.7576 0.7561 0.7546 0.7531 0.7516 0.7500 0.7485

Columns 261 through 270
0.7470 0.7455 0.7440 0.7425 0.7410 0.7395 0.7380 0.7365 0.7351 0.7336

Columns 271 through 280
0.7321 0.7306 0.7292 0.7277 0.7263 0.7248 0.7234 0.7219 0.7205 0.7191

Columns 281 through 290
0.7177 0.7162 0.7148 0.7134 0.7120 0.7106 0.7093 0.7079 0.7065 0.7052

Columns 291 through 300
0.2135 0.2099 0.2064 0.2028 01992 0.1956 0.1921 0.1886 0.1851 0.1816

输出功角增量:DQ=
Columns 1 through 10
0 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000
Columns 11 through 20
-0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000

Columns 21 through 30
-0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 0.0008

Columns 31 through 40
0.0023 0.0039 0.0054 0.0070 0.0086 0.0101 0.0117 0.0132 0.0148 0.0163

Columns 41 through 50
0.0179 0.0194 0.0210 0.0225 0.0241 0.0256 0.0272 0.0288 0.0303 0.0319

Columns 51 through 60
0.0334 0.0350 0.0365 0.0381 0.0396 0.0412 0.0427 0.0443 0.0458 0.0474

Columns 61 through 70
0.0489 0.0505 0.0520 0.0535 0.0551 0.0566 0.0582 0.0597 0.0613 0.0628

Columns 71 through 80
0.0644 0.0659 0.0674 0.0690 0.0705 0.0721 0.0736 0.0751 0.0767 0.0782

Columns 81 through 90
0.0797 0.0813 0.0828 0.0843 0.0859 0.0874 0.0889 00905 0.0920 0.0935

Columns 91 through 100
0.0950 0.0966 0.0981 0.0996 0.1011 0.1027 0.1042 0.1057 0.1072 0.1087

Columns 101 through 110
0.1102 0.1118 0.1133 0.1148 0.1163 0.1178 0.1193 0.1208 0.1223 0.1238

Columns 111 through 120
0.1253 0.1268 0.1283 0.1298 0.1313 0.1328 0.1343 0.1358 0.1373 0.1388

Columns 121 through 130
0.1403 0.1418 0.1433 0.1448 0.1463 0.1477 0.1492 0.1507 0.1522 0.1537

Columns 131 through 140
0.1551 0.1566 01581 0.1596 0.1610 0.1625 0.1640 0.1654 0.1669 0.1684

Columns 141 through 150
0.1698 0.1713 0.1727 0.1742 0.1756 0.1771 0.1785 0.1800 0.1814 0.1829

Columns 151 through 160
0.1843 0.1858 0.1872 0.1886 0.1901 0.1915 0.1929 0.1944 0.1958 0.1972

Columns 161 through 170
0.1986 0. 0. 0.2029 0.2043 0.2057 0.2072 0.2086 0.2100 0.2114

Columns 171 through 180
0.2128 0.2142 0.2156 0.2170 0.2184 0.2198 0.2212 0.2226 0.2239 0.2253

Columns 181 through 190
0.2267 0.2281 0.2295 0.2309 0.2322 0.2336 0.2350 0.2363 0.2377 0.2391

Columns 191 through 200
0.2404 0.2418 0.2432 0.2445 0.2459 0.2472 0.2486 0.2499 0.2513 0.2526

Columns 201 through 210
0.2539 0.2553 0.2566 0.2579 0.2593 0.2606 0.2619 0.2632 0.2646 0.2659

Columns 211 through 220
0.2672 0.2685 0.2698 0.2711 0.2724 0.2737 0.2750 0.2763 0.2776 0.2789

Columns 221 through 230
0.2802 0.2815 0.2828 0.2841 0.2854 0.2867 0.2879 0.2892 0.2905 0.2918

Columns 231 through 240
0.2930 0.2943 0.2955 0.2968 0.2981 0.2993 0.3006 0.3018 0.3031 0.3043

Columns 241 through 250
0.3056 0.3068 0.3080 0.3093 0.3105 0.3117 0.3130 0.3142 0.3154 0.3166

Columns 251 through 260
0.3179 0.3191 0.3203 0.3215 0.3227 0.3239 0.3251 0.3263 0.3275 0.3287

Columns 261 through 270
0.3299 0.3311 0.3323 0.3335 0.3347 0.3359 0.3370 0.3382 0.3394 0.3406

Columns 271 through 280
0.3417 0.3429 0.3441 0.3452 0.3464 0.3475 0.3487 0.3499 0.3510 0.3522

Columns 281 through 290
0.3533 0.3545 0.3556 0.3567 0.3579 0.3590 0.3601 0.3613 0.3624 0.3635

Columns 291 through 300
0.3648 0.3678 0.3709 0.3740 0.3769 0.3793 0.3812 0.3825 0.3832 0.3836

Column 301
0.3837

输出功角:Q=
Columns 1 through 10
4.0069 4.0069 4.0069 4.0069 4.0069 4.0069 4.0069 4.0069 4.0069 4.0069

Columns 11 through 20
4.0069 4.0069 4.0069 4.0069 4.0069 4.0069 4.0069 4.0069 4.0069 4.0069

Columns 21 through 30
4.0069 4.0069 4.0069 4.0069 4.0069 4.0069 4.0069 4.0069 4.0069 4.0077

Columns 31 through 40
4.0100 4.0139 4.0193 4.0263 4.0349 4.0450 4.0567 4.0699 4.0846 4.1010

Columns 41 through 50
4.1188 4.1383 4.1593 4.1818 4.2059 4.2316 4.2588 4.2875 4.3178 4.3497

Columns 51 through 60
4.3831 4.4180 4.4546 4.4926 4.5322 4.5734 4.6161 4.6604 4.7062 4.7535

Columns 61 through 70