当前位置:首页 >> 电力/水利 >>

电力系统运行方式分析和计算


电力系统运行方式分析和计算 设计报告

专 班 学 姓

业: 级: 号: 名:

电气工程及其自动化 11 级电气 1 班 201130213166 杨玉豪 201130212176 潘鸣

华南理工大学电力学院 2015-01-05

0、课程设计题目 A3:电力系统运行方式

分析和计算
姓名: 一、一个 220kV 分网结构和参数如下:
#1

指导教师:

500kV变电站 220kV变电站

30km 火电厂

#2 11km #4 11km #5 25km 16km 9km #6 20km #3

G
500kV 站(#1)的 220kV 母线视为无穷大母线,电压恒定在 230kV。 图中,各变电站参数如下表: 编号 类型 #1 500kV 站 #2 220kV 站 #3 220kV 站 #4 220kV 站 #5 220kV 站 #6 220kV 站 各变电站负荷曲线基本一致。日负荷曲线主要参数为: ? 日负荷率:0.85,日最小负荷系数:0.64

220kV 最大负荷,MVA 平衡节点 230+j40 210+j25 300+j85 410+j110 220+j30

各线路长度如图所示。所有线路型号均为 LGJ-2*300,基本电气参数为: 正序参数: r = 0.054Ω /km, x = 0.308Ω /km, C = 0.0116 ?F/km; 零序参数: r0 = 0.204Ω /km, x0 = 0.968Ω /km, C0 = 0.0078 ?F/km; 40? C 长期运行允许的最大电流:1190A。 燃煤发电厂 G 有三台机组,均采用单元接线。电厂 220kV 侧采用双母接线。发电机组 主要参数如下表(在 PowerWorld 中选择 GENTRA 模型) :
升 压 变 机 组 台 数 1 1 1 单台 容 量 (M W) 300 300 250 额定电 压 ( EV ) 10.5 10.5 10.5 功 率 因 数 0.85 0.85 0.85 容 量 MV A 350 350 300 1.8 1.8 2.1 0.18 0.18 0.2 1.2 1.2 1.5 8 8 7 7 7 6 0.00004 0.00003 0.00003 0.298 0.305 0.321 10.22 10.32 9.38 300 300 250 120 120 100 Xd Xd’ Xq Td0’ TJ= 2H Pmax Pmin (MW)

t/(MW2? h)

ai,2 ai,1 ai,0
t/(MW ?h) t/h

(MW)

升压变参数均为 Vs%=10.5%,变比 10.5kV/242kV。不计内阻和空载损耗。 稳定仿真中发电机采用无阻尼绕组的凸极机模型。不考虑调速器和原动机模型。不考虑 电力系统稳定器模型。励磁系统模型为:

该模型在 PowerWorld 中为 BPA_EG 模型,主要参数如下: KA=40 TA=0.1 TA1=0.1 KF=0.05 TF=0.7 VRmax=3.7 VRmin=0.0 发电厂按 PV 方式运行,高压母线电压定值为 1.05VN。考虑两种有功出力安排方式: ? 满发方式: 开机三台,所有发电机保留 10%的功率裕度; ? 轻载方式: 仅开 250MW 机组,且保留 10%的功率裕度; ? 发电厂厂用电均按出力的 7%考虑。 二、 设计的主要内容:

1、根据负荷变化和机组出力变化,拟定至少两种典型运行方式; (完成) 2、进行参数计算和标幺化,形成两种典型运行方式的潮流计算参数; (完成) 3、用 Matlab 编制潮流计算程序,可任选一种潮流计算方法; (完成) 4、用所编制的潮流程序完成典型运行方式的潮流计算,进行电压和网损分析; (完成) 5、 用 PowerWorld 软件进行潮流计算并与自己编制的软件计算结果进行校核和分析; (完 成) 6、用所编制的潮流程序完成大方式的“N-1”潮流校核,进行线路载流能力和电压水平 分析; (完成) 7、用 Matlab 编制三相短路的短路容量计算程序; (完成) 8、对主要 220kV 母线进行三相短路容量测算,并与 PowerWorld 的计算结果进行校核; (完成) 9、自行选择 2-3 种故障方案,用 PowerWorld 进行稳定计算,给出摇摆曲线,并计算故 障的极限切除时间。 10、 假定电网公司下发给燃煤发电厂 G 的日发电计划曲线如下图, 按照等微增率准则对 三台机组进行经济负荷分配,同时采用 matlab 中的 quadprog 函数对三台机组进行负荷 优化分配,并对两种分配结果进行分析比较。要求给出三台机组的日发电计划曲线。 11、编制课程设计报告

P/ M W
795 600 400

0
三、 设计要求和设计成果:

6

9

21 24

t/h

1、每两位同学为一组,自行分工,但任务不能重复; 2、每位同学对自己的设计任务编写课程设计说明书一份; 3、一组同学共同完成一份完整的设计报告; 4、设计说明和报告应包含: ? 以上设计任务每一部分的计算过程和结果分析; ? 所编制的潮流、短路和机组经济负荷分配源程序(主要语句应加注释) ; ? 潮流计算结果(潮流图) ? 稳定计算的功角曲线等;

1、电力系统参数计算及其标幺化
1.1 电力系统等效电路图

π 型等效电路

1.2 运行方式拟定
1.满载发电负荷最大运行方式:

发电厂: 满发。 取发电机容量的 10%为裕量, 再按已知保留出力的 7%作为厂用电, 即发出功率为总容量的 83.7%。 负荷:采用最大负荷计算。 2.满载发电负荷最小运行方式: 发电厂: 满发。 取发电机容量的 10%为裕量, 再按已知保留出力的 7%作为厂用电, 即发出功率为总容量的 83.7%。 负荷:将最大负荷与日最小负荷系数相乘,得负荷最小值。

1.3 线路参数计算
线路参数给定如下: 正序参数: r = 0.054Ω /km x = 0.308Ω /km C = 0.0116 ?F/km; 零序参数: r0 = 0.204Ω /km x0 = 0.968Ω /km C0 = 0.0078 ?F/km; 线路长度:L12:30km; L23:20km; L24:11km; L36:9km; L45:11km;L6G:16km; L5G:25km; 1.线路参数有名值计算: 按照双回路线路参数考虑,应用如下公式进行有名值计算: R=0.5rl X=0.5xl 0.5B=ωcl 正序参数: 线 L12 L23 L24 L36 L45 L6G L5G 路 R/ 0.81 0.54 0.297 0.243 0.297 0.432 0.675 Ω X/ 4.62 3.08 1.694 1.386 1.694 2.464 3.85 Ω 0.5 0.0001 0.00007 0.00004 0.00003 0.00004 0.00005 0.00009 B/ 093 289 009 280 009 831 111 Ω 零序参数: 线 L12 L23 L24 L36 L45 L6G L5G 路 R/ 3.06 2.04 1.122 0.918 1.122 1.632 2.55 Ω X/ 14.52 9.68 5.324 4.356 5.324 7.744 12.1 Ω 0.5 0.00007 0.00004 0.00002 0.00002 0.00002 0.00003 0.00006 B/ 351 901 695 205 695 921 126 Ω 2.线路参数标幺值计算: 所选基准电压:230KV;基准容量:100MW。即: = 100, = 230 , = 529 应用如下公式进行标幺值计算: ? = 2

2 ? = B 正序参数: 线 L12 路 R 0.00153 1 X 0.00873 3 0.5 0.05781 B 97 零序参数: 线 L12 路 R 0.00578 4 X 0.02744 8 0.5 0.03888 B 69 L23 0.00102 1 0.00582 2 0.03855 88 L23 0.00385 6 0.01829 9 0.02592 63 L24 0.00056 1 0.00320 2 0.02120 76 L24 0.00212 1 0.01006 4 0.01425 66 L36 0.00045 9 0.00262 0 0.01735 12 L36 0.00173 5 0.00823 4 0.01166 45 L45 0.00056 1 0.00320 2 0.02120 76 L45 0.00212 1 0.01006 4 0.01425 66 L6G 0.00081 7 0.00465 8 0.03084 60 L6G 0.00308 5 0.01463 9 0.02074 21 L5G 0.00127 6 0.00727 8 0.04819 72 L5G 0.00482 0 0.02287 3 0.03240 65

1.4 发电机参数计算
采用 10.5KV 作为发电机端的基准电压,230kV 为系统侧的基准电压。 将三台机组分别赋予编号,两个 300MW 机组为 1、2 号,剩余一台 250MW 机组为 3 号。 1.功率输出 PG1=PG2=300×83.7%=251.1MW,PG3=250×83.7%=209.25MW QG1=QG2= ( 0.85 )2 ? 251.12 =155.617803 MVar QG3= (
209.25 2 ) 0.85 251.1

? 209.252 =129.681503MVar

相应的标幺值: PG1*=PG2*=2.511000, PG3*=2.092500 QG1*=QG2*=1.556178, QG3*=1.296815 2.机端电压 VG=10.5kV, 取基准电压 10.5kV,VG*=1 3.相关电抗值归算 应用如下公式进行归算: X * ? X ? VSNN ? VSB2
B 2

机组 1 2

发电容量 300MW 300MW

Xd* 0.5100 0.5100

Xd’* 0.051000 0.051000

Xq* 0.408000 0.408000

3

250MW

0.7140

0.068000

0.510000

1.5 变压器参数计算
分别将与 1、2、3 号机组相连的变压器分别编号为 1、2、3 号。 根 据 题 目 计 算 XT , 相 关 计 算 公 式 为 : k*=(242*10.5)/(230*10.5)=1.05 变压器 容量/MVA 变比 k* XT/Ω 1 350 0.950413 17.569200 2 350 0.950413 17.569200 3 300 0.950413 20.497400
S% X T * ? V100 ? VSNN ? VSB2 B 2



XT* 0.033212 0.033212 0.038747

1.6 负荷参数计算
已知的日负荷率 0.85 和日最小负荷系数 0.64,可得到以下数据: 负荷有名值: 节点 类型 Pmax/MW Qmax/MW Pmin/MW #1 平衡节点 #2 PQ 230 40 147.2 #3 PQ 210 25 134.4 #4 PQ 300 85 192 #5 PQ 410 110 262.4 #6 PQ 220 30 140.8 #7 PQ 0 0 0 #8 PV 300 185.9233 251.1 #9 PV 300 185.9233 251.1 #10 PV 250 154.9361 209.25 负荷标幺值: 节点 类型 电压 Pmax/MW Qmax/MW Pmin/MW #1 平衡节 1 点 #2 PQ 1 2.3 0.4 1.472 #3 PQ 1 2.1 0.25 1.344 #4 PQ 1 3 0.85 1.92 #5 PQ 1 4.1 1.1 2.624 #6 PQ 1 2.2 0.3 1.408 #7 PQ 1 0 0 0 #8 PV 1.05 3 1.859233 2.511 #9 PV 1.05 3 1.859233 2.511 #10 PV 1.05 2.5 1.549361 2.0925

Qmin/MW 25.6 16 54.4 70.4 19.2 0 155.6178 155.6178 129.6815 Qmin/MW 0.256 0.16 0.544 0.704 0.192 0 1.556178 1.556178 1.296815

2、潮流计算的 Matlab 编程及 PowerWorld 仿真
2.1 牛顿拉夫逊法 Matlab 计算程序
牛顿拉夫逊法是计算潮流时常见的方法,该方法具有广泛的应用。该方法是 牛顿在 17 世纪提出的一种在实数域和复数域上近似求解方程的方法。多数方程 不存在求根公式,因此求精确根非常困难,甚至不可能,从而寻找方程的近似根 就显得特别重要。 方法使用函数 f(x)的泰勒级数的前面几项来寻找方程 f(x) = 0 的根。牛顿迭代法是求方程根的重要方法之一,其最大优点是在方程 f(x) = 0 的单根附近具有平方收敛,而且该法还可以用来求方程的重根、复根。 相应程序详见附件 1. 2.1.1 满发满载运行方式下计算 在这种运行方式下,三台机组全部开启,各自保留 10%作为功率裕度,重用 负荷占发电量 7%,负荷为最大负荷。此方式下,节点 1 为平衡节点,节点 2、3、 4、5、6、7 为 PQ 节点,节点 8、9、10 为 PV 节点。 经过 Matlab 的运算,我们得到如下结果:

2.1.2 满发轻载运行方式下计算 经 Matlab 运算,我们得到:

2.2 PowerWorld 仿真电路图
2.2.1 满发满载方式运行下电力系统的潮流计算 满发满载方式运行下,发电厂三台机组均满载,取发电机容量 10%为裕量, 并保留出力的 7%作厂用电, 即发出功率为总容量的 83.7%,负荷采用最大负荷计 算。满发满载电路图:

编辑模式下,#1 为平衡节点,电压恒定为 1.05VN。为了反映这一特性,#1 应接一台容量无穷大的发电机。具体参数设置如图 3 所示(机组出力暂不确定, 可填为 100MW。 )

接在平衡节点处的发电机参数设置 #7 为 PV 节点,P=711.45MW,V=1.05VN,正常运行时各机组发出的有功功率 P1=P2 =251.1MW,P3=209.25MW。具体参数设置如图 4、图 5 所示。

接在#7 处的发电机#8、#9 参数设置

接在#7 处的发电机#10 参数设置 其余节点为 PQ 节点,将基准电压设置为 230kV,带上相应大小的负荷。最 后,用输电线路将各节点连接起来。 潮流仿真结果如下: 各节点电压 节点编号 基准电压 kV 标幺电压 实际电压(kV) 相角 (度) 1 230 1.00000 230.0000 0 2 230 0.9992 229.816 -1.70 3 230 1.0021 230.483 -1.88 4 230 0.9985 229.655 -2.02 5 230 0.9998 229.954 -2.08 6 230 1.0041 230.943 -1.80 7 230 1.0090 232.07 -1.39 各发电机运行状态 节点编号 发电机编号 有功出力(MW) 无功出力(MVar) 电压标幺值 1 1 664.18 -100.66 1.0000 8 8 251.10 139.61 1.0500 9 9 251.10 139.61 1.0500 10 10 209.25 119.20 1.0500 各支路状态 首端节 末端节 首端有功 首端无功功 首端视在功 有功损 无功损耗 点编号 点编号 功率(MW) 率(Mvar) 率(MVA) 耗(MW) (Mvar) 1 2 332.1 -50.3 335.9 1.72 -1.75 2 3 43.3 -60.8 74.6 0.05 -7.42 2 4 172.1 -7.8 172.3 0.17 -3.28 3 6 -61.8 -65.9 90.3 0.04 -3.28 4 5 21.9 -47 51.8 0.01 -4.15

5 7 -183.1 6 7 -171.8 电网线损率计算: 线损率按以下公式计算: 线损率= 支路 #1#2 #2#3 线 损 0.518% 0.115% 率 电网网损率的计算:

-97.8 -77.6

207.6 188.5

0.54 0.28

-6.65 -4.63

线路首端有功功率-线路末端有功功率 ?100% 线路首端有功功率

#2#4 0.099%

#3#6 0.065%

#4#5 0.046%

#5#7 #6#7 0.295% 0.163%

网损率按以下公式计算: 网损率=

发电机输出功率PD -负荷PL ?100% 发电机输出功率PD

其中, 发电机输出功率包括发电厂和平衡节点的功率,发电机输出功率和负荷均 只取有功分量。则发电机输出功率 PD=711.45MW+664MW=1375.45MW 总负荷 PL =230MW+210MW+300MW+410MW+220MW =1370MW 网损率=(PD-PL)/PD×100%=0.396% 2.2.2 满发轻载方式运行下电力系统的潮流计算 满发轻载方式运行下,发电厂三台机组均满载,取发电机容量 10%为裕量, 并保留出力的 7%作为厂用电, 即发出功率为总容量的 83.7%,负荷采用最小负荷 计算。

满发轻载电路图 潮流仿真结果如下: 各节点电压

节点编号 基准电压 kV 标幺电压 实际电压(kV) 相角 (度) 1 230 1.00000 230.0000 0 2 230 1.0056 231.2880 -0.49 3 230 1.0093 232.1390 -0.39 4 230 1.0061 231.4030 -0.58 5 230 1.0079 231.8170 -0.51 6 230 1.0113 232.5990 -0.26 7 230 1.0157 233.6110 0.17 各发电机运行状态 节点编号 发电机编号 有功出力(MW) 无功出力(MVar) 电压标幺值 1 1 167.57 -169.22 1.0000 8 8 251.10 118.21 1.0500 9 9 251.10 118.21 1.0500 10 10 209.25 100.86 1.0500 各支路状态 首端节 末端节 首端有功 首端无功功 首端视在功 有功损 无功损耗 点编号 点编号 功率(MW) 率(Mvar) 率(MVA) 耗(MW) (Mvar) 1 2 83.8 -84.6 119.1 0.2 -10.47 2 3 -37.6 -60.8 71.3 0.05 -7.56 2 4 47.6 -26.4 54.4 0.02 -4.20 3 6 -104.9 -61.0 121.3 0.07 -3.17 4 5 -48.4 -49.4 69.2 0.03 -4.16 5 7 -179.6 -80.4 196.8 0.48 -7.15 6 7 -175.3 -67.4 187.8 0.28 -4.75 电网线损率计算: 线损率按以下公式计算: 线损率= 支路 #1#2 #2#3 线 损 0.239% 0.133% 率 电网网损率的计算
线路首端有功功率-线路末端有功功率 ?100% 线路首端有功功率

#2#4 0.042%

#3#6 0.067%

#4#5 0.062%

#5#7 #6#7 0.267% 0.160%

网损率按以下公式计算: 网损率=

发电机输出功率PD -负荷PL ?100% 发电机输出功率PD

其中, 发电机输出功率包括发电厂和平衡节点的功率,发电机输出功率和负荷均 只取有功分量。则发电机输出功率 PD=711.45MW+167.571MW=879.021MW 总负荷 PL =147.2MW+134.4MW+192MW+262.4MW+140.8MW =876.8MW 网损率=(PD-PL)/PD×100%=0.253%

2.3 Matlab 和 PowerWorld 计算的结果比较
以大方式的结果为例:电压标幺值很接近;电角度有一定差别,但 一般小数点第一位都相同;线损有一定差别,但也都不大。

造成结果差别的原因, 首先, 两种程序输入参数上存在误差, Matlab 的参数会保留一定位小数,不能达到最精确值,而 PowerWorld 会将线 路参数强制以美制单位保存, 每次查看和修改参数自动转换成公制单位 造成误差, 而且这种误差还会随着查看与修改的次数累积; 使用 Matlab 的算法较为理论和简单, 系统中各个元件和部分的反馈和相互作用没什 么体现,而 PowerWorld 仿真的算法更为精确和复杂,如发电机的模型 更加详细。

3、大方式下“N-1”潮流校核
使用满发大负荷运行方式,即:三台机组保留 10%容量备用,发出功率的 7% 为厂用电,出力为机组总容量的 83.7%;负荷使用提供的 220kV 最大负荷。 N-1 情况下, 将 3 发电机组及其变压器统一以#7 节点代表进行计算。依次分 别断开各条支路中双回线的一条线,计算节点电压和各个功率及损耗;根据设计 要求:40°C 长期运行允许的最大电流:1190A。

3.1 Matlab“N-1”潮流校核
N-1 故障可理解为双回线中的一回被切除,只剩一回在维持运行。这种情况 会造成电力系统潮流状态的变更,需要进行潮流校核来检验系统的稳定性。 此方式中, 我们依次分别将其中一条双回线参数改成相应的单回线参数带入 Matlab 计算。 经过 Matlab 运算,得到以下结果: N-1 故障下各节点电压 故障线 #1 #2 #3 #4 #5 #6 #7 路 L12 1.000 1.011 1.020 1.011 1.016 1.026 1.039 L23 1.000 1.007 1.022 1.008 1.014 1.027 1.039 L24 1.000 1.008 1.018 1.009 1.014 1.024 1.037 L36 1.000 1.007 1.016 1.008 1.014 1.026 1.039 L45 1.000 1.008 1.018 1.008 1.016 1.024 1.039 L57 1.000 1.005 1.018 1.003 1.005 1.025 1.041 L67 1.000 1.006 1.013 1.008 1.014 1.018 1.040 N-1 故障各线路电流 故障线 L12 L23 L24 L36 L45 L57 L67 路 L12 2.84 0.96 1.37 0.67 1.16 2.09 1.7 L23 3.65 0.57 1.59 0.84 1.03 2.61 1.39 L24 3.69 0.75 1.66 0.7 1.02 2.73 1.22 L36 3.65 1.03 1.35 0.9 1.08 2.65 1.32 L45 3.65 0.84 1.35 0.63 1.25 2.47 1.49 L57 3.65 0.78 1.97 0.95 1.21 2.21 1.97

L67 3.65 0.98 1.04 1.06 1.08 2.95 1.47 根据基准值的选取,可以得到题目已知的 40℃长期允许最大电流 1190A 的 标幺值为 2.737。

3.2 PowerWorld“N-1”潮流校核

断线电路图(以 L12 为例) L12 断开一路: 节点 电压 基准 值 (kV) 标幺 值 有名 值 (kV) 功角 (°) 线路 电流值 (kA) 1 23 0 1 23 0 0 12 2 230 0.9978 5 229.50 6 -3.4 23 节点电压和功角 3 4 230 1.0008 7 230.20 0 -3.6 230 0.9971 5 229.34 4 -3.7 5 230 0.9985 4 229.66 3 -3.8 45 6 230 1.0028 9 230.66 5 -3.5 57 1.0078 8 231.81 3 -3.1 7 230

一回线相电流 24 36

1.6874 0.1923 0.4337 0.2307 0.1354 0.5244

67 0.4743

分析:由于题目限定 40°C 长期运行允许最大电流是 1190A,而表中 12 回线路的断线电流为 1687.4A,已经超过极限,说明该线路载流能力较 差, 支路 12 不能承受断一回线路的电流,而对于各母线的电压水平均无 大变化。 L23 断开一路: 节点 电压 基准 值 (kV) 标幺 值 有名 值 (kV) 功角 (°) 线路 电流值 (kA) L24 断开一路: 节点 电压 基准 值 (kV) 标幺 值 有名 值 (kV) 功角 (°) 线路 电流值 (kA) L36 断开一路: 1 23 0 1 23 0 0 12 0.843 2 230 0.9991 2 229.79 8 -1.7 23 0.2199 节点电压和功角 3 4 230 1.0018 9 230.43 4 -1.94 230 0.9977 1 229.47 3 -2.27 5 230 0.9991 3 229.80 1 -2.29 45 6 230 1.0038 0 230.87 3 -1.89 57 1.0085 8 231.97 4 -1.53 7 230 1 23 0 1 23 0 0 12 2 230 0.9987 3 229.70 9 -1.7 23 节点电压和功角 3 4 230 1.0032 3 230.74 3 -1.99 230 0.9981 9 229.58 3 -2.03 5 230 0.9997 3 229.93 8 -2.11 45 6 230 1.0049 5 231.14 0 -1.90 57 1.0094 2 232.16 7 -1.45 7 230

一回线相电流 24 36

0.8411 0.2945 0.4539 0.2228 0.1581 0.5148

67 0.4791

一回线相电流 24 36 0.773

0.2003 0.1188 0.5629

67 0.4312

节点 电压 基准 值 (kV) 标幺 值 有名 值 (kV) 功角 (°) 线路 电流值 (kA) L45 断开一路: 节点 电压 基准 值 (kV) 标幺 值 有名 值 (kV) 功角 (°) 线路 电流值 (kA) L57 断开一路: 节点 电压 基准 值 (kV)

1 23 0 1 23 0 0 12

2 230 0.9989 4 229.75 6 -1.7 23

节点电压和功角 3 4 230 1.0013 1 230.30 2 -1.89 230 0.9983 4 229.61 8 -2.01

5 230 0.9998 2 229.95 8 -2.06 45

6 230 1.0048 3 231.11 0 -1.75 57

7 230

1.0093 8 232.15 7 -1.35

一回线相电流 24 36

0.8419 0.181

0.4198 0.4055 0.1352 0.5383

67 0.4524

1 23 0 1 23 0 0 12 0.842

2 230 0.9989 9 229.76 9 -1.7 23 0.199

节点电压和功角 3 4 230 1.0020 8 230.47 9 -1.89 230 0.9980 2 229.54 4 -2.01

5 230 1.0003 1 230.07 1 -2.11 45

6 230 1.0041 3 230.95 1 -1.82 57

7 230

1.0091 7 232.11 0 -1.41

一回线相电流 24 36

0.4265 0.2296 0.2246 0.5208

67 0.4700

1 23 0

2 230

节点电压和功角 3 4 230 230

5 230

6 230

7 230

标幺 值 有名 值 (kV) 功角 (°)

1 23 0 0

0.9980 8 229.55 7 -1.69

1.0021 4 230.49 2 -1.75

0.9962 3 229.13 3 -2.08

0.9964 9 229.19 2 -2.20

1.0046 3 231.06 5 -1.62

1.0104 7 232.40 8 -1.11

线路 电流值 (kA) L67 断开一路: 节点 电压 基准 值 (kV) 标幺 值 有名 值 (kV) 功角 (°)

12

23

一回线相电流 24 36

45

57

0.8397 0.1870 0.5350 0.3212 0.1628 0.8095

67 0.5740

1 23 0 1 23 0 0

2 230 0.9985 4 229.66 4 -1.70

节点电压和功角 3 4 230 1.0002 8 230.06 4 -1.95 230 0.9981 4 229.57 2 -1.97

5 230 0.9998 2 229.95 9 -1.98

6 230 1.0017 3 230.39 8 -1.91

7 230

1.0098 5 232.26 6 -1.19

线路 电流值 (kA)

12

23

一回线相电流 24 36

45

57

0.8409 0.2074 0.371 0.1565 0.1359 0.5855

67 0.7971

3.3 Matlab 与 PowerWorld“N-1”潮流校核对比:
通过就“N-1”潮流校核分别进行 Matlab 和 PowerWorld 计算、仿真,我们 发现结果差距较大,其中 PowerWorld 仿真下只有 L12 电流值超过给定值,而 Matlab 计算下过流线路较多。经过反复计算、仿真,我们大胆认为中间存在的 较大差异并非本身操作失误, 而是软件在运算中机理不同而生:节点导纳矩阵形 成过程不同。

4、220kV 母线的三相短路容量测算
选用满发重荷的运行方式,平衡机基准容量为 100MVA,发电机基准容量为 各自的额定容量, 变压器改为原来的标幺变比和标幺阻抗,对非平衡母线的所有

220kv 的母线进行三相短路测算。

4.1 短路容量的 Matlab 计算
相应 Matlab 计算程序详见附件 2. 4.1.1 满发满载运行方式下短路容量计算 经过 Matlab 运算,分别得到短路电流和短路容量:
节点 电流 容量 #1 13.787 1390 #2 14.995 1500 #3 14.998 1050 #4 15.093 1510 #5 15.139 1510 #6 15.069 1510 #7 15.291 1530

4.1.2 满发轻载运行方式下短路容量计算 经过 Matlab 运算,分别得到短路电流和短路容量:
节点 电流 容量 #1 9.693 970 #2 10.343 1030 #3 10.370 1040 #4 10.387 1040 #5 10.432 1040 #6 10.427 1040 #7 10.591 1060

4.2 短路容量的 PowerWorld 仿真

短路容量仿真操作 4.2.1 满发满载运行方式下短路容量计算 经过 PowerWorld 仿真,分别得到短路电流和短路容量:
节点 容量 #1 #2 #3 #4 #5 #6 #7

1433.7

1475.5

1475.5

1476.4

1478.8

14.795

1492.6

4.2.2 满发轻载运行方式下短路容量计算 经过 PowerWorld 仿真,分别得到短路电流和短路容量:
节点 容量 #1 #2 #3 #4 #5 #6 #7

961.5

986.4

989.5

986.5

990.4

992.8

1001.6

4.3 两种计算方法结果比较

使用 Matlab 和 PowerWorld 进行计算和仿真的结果有一定相关性, 但是存在一定差别。 我们认为除了两者输入各种参数的误差之外,还可 能是导纳矩阵算法之间的差距。PowerWorld 形成的节点导纳矩阵与我 们处理的不同,短路电流的计算准确。Matlab 的算法是我们根据书本 概念按自己的理解写成的程序代码,核心算法也遵照书本所述,可能是 我们对负荷模型处理不够准确造成了这种差别。

5、PowerWorld 的网络稳定扫描
如果短路点越靠近发电机节点,则短路故障对电力系统扰动越大。因此在以 下稳定扫描仿真中,均把短路点尽可能靠近母线 7 上。若在此情况下,发电机摇 摆曲线仍能保持暂态稳定,则电力系统稳定性符合要求。

5.1 暂态参数设置
先按照要求添加发电机的暂态模型,设置好发电机的暂态参数。线路的零 序阻抗参数。具体参数设置如下图:

发电机#8、#9 GENTRA 模型参数

发电机#10 GENTRA 模型参数

发电机#8、#9、#10 BPA_EG 模型参数

发电机升压变压器零序参数设定

线路零序参数设定(以 L67 为例)

5.2 故障方案
故障均设置在 t=1.0s 发生,仿真时间为 20s。 故障方案一:母线 7 三相短路,1.10s 切除故障线路。

故障方案一设计

故障方案一运行

从故障切除时间为 0.1s 的各发电机转子角摇摆曲线,我们可以看出个发电 机的转子角由于故障而产生扰动,而在切除故障之后逐渐趋于稳定,各转子角的 差(功角)趋于稳定,所以系统暂态稳定。 故障极限切除时间求解: 用二分法通过不断地设置切除时间, 知道系统在故障切除之后不能保持暂态 稳定性,到达临界值,从而确定大致的系统故障的极限切除时间。 故障切除时:1.4218s

2.4218s 功角图

2.4219s 功角图 故障方案二:线路 6-7 三相短路,1.20s 切除故障线路。

从故障切除时间为 0.2s 的各发电机转子角摇摆曲线,我们可以看出个发电 机的转子角由于故障而产生扰动,而在切除故障之后逐渐趋于稳定,各转子角的 差(功角)趋于稳定,所以系统暂态稳定。 故障极限切除时间求解: 用二分法通过不断地设置切除时间, 知道系统在故障切除之后不能保持暂态 稳定性,到达临界值,从而确定大致的系统故障的极限切除时间。 故障切除时:1.4217s

2.4217s 功角图

2.4218s 功角图

5.3 方案总结
由此可以得出故障的极限切除时间的大致区间为[1.4217, 1.4218],我们可 以取故障的极限切除时间为:1.4217s。

6、发电机组的负荷经济分配
6.1 遵守等微增率原则的经济分配方案
F(P)=ai2p2+ai1p+ai0 F1(P1)=0.00004P2+0.298P+10.22 F2(P2)=0.00034P2+0.305P+10.32 F3(P3)=0.00004P2+0.321P+9.38 150≤P1≤300 150≤P2≤300 100≤P3≤250

由题目可知,三个发电厂的微增率表达式如下所示: Λ 1=0.00008P1+0.298 Λ 2=0.00006P1+0.305 Λ 3=0.00006P1+0.321 令λ 1=λ 2=λ 3, ①P1+P2+P3=400, P1=150; P2=150; P3=100. ②P1+P2+P3=600, P1=300; P2=2000; P3=100.

③P1+P2+P3=795, P1=300; P2=300; P3=195.

6.2 使用 Matlab 函数进行优化的分配方案
使用 Matlab 优化的过程中,我们按要求使用 qua 函数进行优化分配, 这个函数专门适用于进行二次优化,使用广泛,准确度高。 详情请见附件 3. 得出以下结果: ①P1=1.7999×102 ②P1=300 P2=200 P2=1.2×102 P3=100 P2=2.99999×102 P3=1.95×102 P3=1.00×102

③P1=2.999999×102

7、课程设计总结
关于电力系统的计算,很重要的一个步骤就是形成节点导纳矩阵。 该矩阵的正确与否能很大程度上决定运算结果的正确性。 在计算和仿真 的过程中,PowerWorld 的导纳矩阵形成的算法更精确,跟 Matlab 上常 用的代码算出来的结果有所不同, 很可能这些不同造成了两种方法的结 果差异,比如对“N-1”潮流计算,结果的大相径庭使我们对此产生了 深刻的认识。 关于 PowerWorld,参数的正确设置非常重要,标幺值和有名值的 区分很关键,只要这几点做好,得出来的结果比 Matlab 运算的更加可 靠。 经过这次课程设计,我们学会了常用的电力系统仿真软件的使用。 提前熟悉了很多行业知识, 为将来就业后可能遇到的跟专业知识相关的 简单情况模拟和处理,提高了我们处理问题的能力。 经过这次课程设计, 我们初步了解了进行电力系统分析的大致步骤 和方法, 将电力系统课程所学知识联系起来解决实际问题。虽然在完成 设计任务的过程中所涉及到的工具在之前的学习中没有学到, 但应用到 的理论、 原理却是跟我们书本, 甚至是跟实际的电力系统研究息息相关。 这加深了我们对以前所学知识的理解。总的来说,这次课程设计的经历 是是一笔宝贵的财富,让我们获益匪浅。

8、附件
附件一:牛顿拉夫逊法 Matlab 计算程序
% 牛拉法计算潮流程序 %---------------------------------------------------------------------% B1 矩阵:1、支路首端号;2、末端号;3、支路阻抗;4、支路对地电纳 % 5、支路的变比;6、支路首端处于 K 侧为 1,1 侧为 0 % B2 矩阵:1、该节点发电机功率;2、该节点负荷功率;3、节点电压初始值 % 4、PV 节点电压 V 的给定值;5、节点所接的无功补偿设备的容量 % 6、节点分类标号:1 为平衡节点(应为 1 号节点) ;2 为 PQ 节点;3 为 PV 节点; %----------------------------------------------------------------------clear all; format long; n=input('请输入节点数:nodes='); nl=input('请输入支路数:lines='); isb=input('请输入平衡母线节点号:balance='); pr=input('请输入误差精度:precision='); B1=input('请输入由各支路参数形成的矩阵:B1='); B2=input('请输入各节点参数形成的矩阵:B2='); Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1= zeros(nl); %-----------------------------------------------------------------for i=1:nl %支路数 if B1(i,6)==0 %左节点处于 1 侧 p=B1(i,1);q=B1(i,2); else %左节点处于 K 侧 p=B1(i,2);q=B1(i,1); end Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5)); %非对角元 Y(q,p)=Y(p,q); %非对角元 Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4); %对角元 K 侧 Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4); %对角元 1 侧 end %求导纳矩阵 disp('导纳矩阵 Y='); disp(Y) %------------------------------------------------------------------G=real(Y);B=imag(Y); %分解出导纳阵的实部和虚部 for i=1:n %给定各节点初始电压的实部和虚部 e(i)=real(B2(i,3)); f(i)=imag(B2(i,3));

V(i)=B2(i,4); %PV 节点电压给定模值 end for i=1:n %给定各节点注入功率 S(i)=B2(i,1)-B2(i,2); %i 节点注入功率 SG-SL B(i,i)=B(i,i)+B2(i,5); %i 节点无功补偿量 end %-------------------------------------------------------------------P=real(S);Q=imag(S); %分解出各节点注入的有功和无功功率 ICT1=0;IT2=1;N0=2*n;N=N0+1;a=0; %迭代次数 ICT1、a;不满足收敛要求的节 点数 IT2 while IT2~=0 % N0=2*n 雅可比矩阵的阶数;N=N0+1 扩展 列 IT2=0;a=a+1; for i=1:n if i~=isb %非平衡节点 C(i)=0;D(i)=0; for j1=1:n C(i)=C(i)+G(i,j1)*e(j1)-B(i,j1)*f(j1);%Σ (Gij*ej-Bij*fj) D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1);%Σ (Gij*fj+Bij*ej) end P1=C(i)*e(i)+f(i)*D(i);% 节 点 功 率 P 计 算 ei (Gij*ej-Bij*fj)+fiΣ (Gij*fj+Bij*ej) Q1=C(i)*f(i)-e(i)*D(i);% 节 点 功 率 Q 计 算 fi (Gij*ej-Bij*fj)-eiΣ (Gij*fj+Bij*ej) %求 i 节点有功和无功功率 P',Q'的计算值 V2=e(i)^2+f(i)^2; %电压模平方 % 以 下 针 对 非 PV 节 点 来 求 取 功 率 差 及 Jacobi 矩 阵 元 ----------------------------if B2(i,6)~=3 %非 PV 节点 DP=P(i)-P1; %节点有功功率差 DQ=Q(i)-Q1; %节点无功功率差 % 以 上 为 除 平 衡 节 点 外 其 它 节 点 的 功 率 计 --------------------------------------

Σ Σ





% 求 取 Jacobi 矩 阵 ---------------------------------------------------------for j1=1:n if j1~=isb&j1~=i %非平衡节点&非对角元 X1=-G(i,j1)*e(i)-B(i,j1)*f(i); % dP/de=-dQ/df X2=B(i,j1)*e(i)-G(i,j1)*f(i); % dP/df=dQ/de

X3=X2; % X2=dp/df X3=dQ/de X4=-X1; % X1=dP/de X4=dQ/df p=2*i-1;q=2*j1-1; J(p,q)=X3;J(p,N)=DQ;m=p+1; % X3=dQ/de J(p,N)=DQ 节点无功功率差 J(m,q)=X1;J(m,N)=DP;q=q+1; % X1=dP/de J(m,N)=DP 节点有功功率差 J(p,q)=X4;J(m,q)=X2; % X4=dQ/df X2=dp/df elseif j1==i&j1~=isb %非平衡节点&对角元 X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/de X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dP/df X3=D(i)+B(i,i)*e(i)-G(i,i)*f(i); % dQ/de X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);% dQ/df p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;% 扩展列△ Q m=p+1; J(m,q)=X1;q=q+1;J(p,q)=X4;J(m,N)=DP;%扩展列△P J(m,q)=X2; end end else % 下 面 是 针 对 PV 节 点 来 求 取 Jacobi 矩 阵 的 元 素 ----------------------------------------DP=P(i)-P1; % PV 节点有功误差 DV=V(i)^2-V2; % PV 节点电压误差 for j1=1:n if j1~=isb&j1~=i %非平衡节点&非对角元 X1=-G(i,j1)*e(i)-B(i,j1)*f(i); % dP/de X2=B(i,j1)*e(i)-G(i,j1)*f(i); % dP/df X5=0;X6=0; p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV; % PV 节 点电压误差 m=p+1; J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; % PV 节 点有功误差 J(m,q)=X2; elseif j1==i&j1~=isb %非平衡节点&对角元 X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);% dP/de X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);% dP/df X5=-2*e(i); X6=-2*f(i); p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV; % PV 节 点电压误差 m=p+1;

J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6; % PV 节 点有功误差 J(m,q)=X2; end end end end end %以上为求雅可比矩阵的各个元素及扩展列的功率差或电压差 --------------------------------------for k=3:N0 % N0=2*n (从第三行开始,第一、二 行是平衡节点) k1=k+1;N1=N; % N=N0+1 即 N=2*n+1 扩展列△P、△Q 或△U for k2=k1:N1 % 从 k+1 列的 Jacobi 元素到扩展列的△ P、△Q 或△U J(k,k2)=J(k,k2)./J(k,k);% 用 K 行 K 列对角元素去除 K 行 K 列后的非对角元素进行规格化 end J(k,k)=1; % 对角元规格化 K 行 K 列对角元素赋 1 % 回 代 运 算 ------------------------------------------------------------------if k~=3 % 不是第三行 k > 3 k4=k-1; for k3=3:k4 % 用 k3 行从第三行开始到当前行的前一行 k4 行消去 for k2=k1:N1 % k3 行后各 行上三角元素 J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去运算(当前行 k 列元素消为 0) end %用当前行 K2 列元素减去当前行 k 列元素乘以第 k 行 K2 列元素 J(k3,k)=0; %当前行第 k 列元素已消为 0 end if k==N0 %若已到最后一行 break; end % 前 代 运 算 -----------------------------------------------------------for k3=k1:N0 % 从 k+1 行到 2*n 最后一行 for k2=k1:N1 % 从 k+1 列到扩展列消去 k+1 行后各行下 三角元素 J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去运算 end %用当前行 K2 列元素减去当前行 k 列元素乘以第 k 行

K2 列元素 J(k3,k)=0; %当前行第 k 列元素已消为 0 end %是第三行 k=3 % 第 三 行 k=3 的 前 代 运 算 ---------------------------------------------------for k3=k1:N0 %从第四行到 2n 行(最后一行) for k2=k1:N1 %从第四列到 2n+1 列(即扩展列) J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);%消去运算 (当前 行 3 列元素消为 0) end %用当前行 K2 列元素减去当前行 3 列元素乘以第三行 K2 列元素 J(k3,k)=0; %当前行第 3 列元素已消为 0 end end end %上面是用线性变换方式高斯消去法将 Jacobi 矩阵化成单位矩阵 %---------------------------------------------------------------------------------for k=3:2:N0-1 L=(k+1)./2; e(L)=e(L)-J(k,N); %修改节点电压实部 k1=k+1; f(L)=f(L)-J(k1,N); %修改节点电压虚部 end %修改节点电压--------------------------for k=3:N0 DET=abs(J(k,N)); if DET>=pr %电压偏差量是否满足要求 IT2=IT2+1; %不满足要求的节点数加 1 end end ICT2(a)=IT2; %不满足要求的节点数 ICT1=ICT1+1; %迭代次数 end %用高斯消去法解"w=-J*V" disp('迭代次数:'); disp(ICT1); disp('没有达到精度要求的个数:'); disp(ICT2); for k=1:n V(k)=sqrt(e(k)^2+f(k)^2); %计算各节点电压的模值 sida(k)=atan(f(k)./e(k))*180./pi; %计算各节点电压的角度 E(k)=e(k)+f(k)*1i; %将各节点电压用复数表示 else

end %计算各输出量-----------------------------------------------------disp('各节点的实际电压标幺值 E 为:'); disp(E); %显示各节点的实际电压标幺值 E 用复数表示 disp('-----------------------------------------------------'); disp('各节点的电压大小 V 为:'); disp(V); %显示各节点的电压大小 V 的模值 disp('-----------------------------------------------------'); disp('各节点的电压相角 deg 为:'); disp(sida); %显示各节点的电压相角 for p=1:n C(p)=0; for q=1:n C(p)=C(p)+conj(Y(p,q))*conj(E(q)); %计算各节点的注入电流的共 轭值 end S(p)=E(p)*C(p); %计算各节点的功率 S = 电压 X 注入电流的共轭值 end disp('各节点的功率 S 为:'); disp(S); %显示各节点的注入功率 disp('-----------------------------------------------------'); disp('各条支路的首端功率 Si 为:'); for i=1:nl p=B1(i,1);q=B1(i,2); if B1(i,6)==0 Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5))... -conj(E(q)))*conj(1./(B1(i,3)*B1(i,5)))); Siz(i)=Si(p,q); else Si(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)./B1(i,5))... -conj(E(q)))*conj(1./(B1(i,3)*B1(i,5)))); Siz(i)=Si(p,q); end disp(Si(p,q)); SSi(p,q)=Si(p,q); ZF=['S(',num2str(p),',',num2str(q),')=',num2str(SSi(p,q))]; disp(ZF); disp('-----------------------------------------------------'); end disp('各条支路的末端功率 Sj 为:'); for i=1:nl p=B1(i,1);q=B1(i,2);

if B1(i,6)==0 Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5))... -conj(E(p)))*conj(1./(B1(i,3)*B1(i,5)))); Sjy(i)=Sj(q,p); else Sj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)*B1(i,5))... -conj(E(p)))*conj(1./(B1(i,3)*B1(i,5)))); Sjy(i)=Sj(q,p); end disp(Sj(q,p)); SSj(q,p)=Sj(q,p); ZF=['S(',num2str(q),',',num2str(p),')=',num2str(SSj(q,p))]; disp(ZF); disp('-----------------------------------------------------'); end disp('各条支路的功率损耗 DS 为:'); for i=1:nl p=B1(i,1);q=B1(i,2); DS(i)=Si(p,q)+Sj(q,p); disp(DS(i)); DDS(i)=DS(i); ZF=['DS(',num2str(p),',',num2str(q),')=',num2str(DDS(i))]; disp(ZF); disp('-----------------------------------------------------'); end

附件二:三相短路容量计算的 Matlab 程序代码
%短路容量计算程序 %-------------------------------------------------------------------% B1 矩阵:1、支路首端号;2、末端号;3、支路阻抗;4、支路对地电纳 % 5、支路的变比;6、支路首端处于 K 侧为 1,1 侧为 0

% B2 矩阵:1、该节点发电机功率;2、该节点负荷功率;3、节点电压初始值 % % 4、PV 节点电压 V 的给定值;5、节点所接的无功补偿设备的容量 6、节点分类标号:1 为平衡节点(应为 1 号节点);2 为 PQ 节点;

%Yd 为修改后的节点导纳矩阵 %----------------------------------------------------------------------

clear all; format long;

g1=input('300MW 发电机数:g1='); g2=input('250MW 发电机数:g2='); n=input('请输入节点数:n='); nl=input('请输入支路数:nl='); B1=input('请输入由各支路参数形成的矩阵:B1='); B2=input('请输入各节点参数形成的矩阵:B2=');

Y=zeros(n);

% Y 为修改前节点导纳矩阵 for i=1:nl if B1(i,6)==0 p=B1(i,1);q=B1(i,2); else p=B1(i,2);q=B1(i,1); end %左节点处于 K 侧 %支路数 %左节点处于 1 侧

Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5)); Y(q,p)=Y(p,q);

%非对角元 %非对角元 %对角元 K 侧 %对角元

Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4); Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4); 1侧 end

%---------------------------------------------------------%Y2-Y6 为各 PQ 节点负荷的导纳

Y1=0;Y2=conj(B2(2,2));Y3=conj(B2(3,2));Y4=conj(B2(4,2));Y5=conj(B2(5, 2)); Y6=conj(B2(6,2)); Xd300=0.51j/0.950413^2; XT300=0.033212j/0.950413^2; Xd250=0.714j/0.950413^2; XT250=0.038747j/0.950413^2; Y7=g1/(XT300+Xd300)+g2/(XT250+Xd250); %处理相应的负荷及机组部分的导纳 %------------------------------------------------------------

C=[Y1,Y2,Y3,Y4,Y5,Y6,Y7]; Yd=Y;

for i=1:n Yd(i,i)=Yd(i,i)+C(i); end disp(Yd); Z = inv(Yd); %求节点阻抗矩阵 %修改各节点自导纳

for j=1:n I(j) = 1/Z(j,j); S(j)=abs(I(j)); Sn(j)=S(j)*100; %短路电流有名值 end %计算完毕---------------------------------------------------%电压故障前电压标幺值为 1

disp('各节点短路时的短路电流幅值标幺值') disp(abs(I)) disp('短路容量有名值 Sn='); disp(Sn);

附件三:使用 Matlab 中 quadprog 函数进行优化的代码
function qua( %主程序 clear all; format long; )

H=input('请输入由发电机参数 a2 除以 2 形成的对角矩阵:a2='); %输入 a2 f=input('请输入由发电机参数 a1 形成的行矩阵:a1='); %输入 a1 UB=input('请输入由发电机最大输出功率形成的列矩阵:UB='); %输入 UB LB=input('请输入由发电机最小输出功率形成的列矩阵:LB='); %输入 LB Beq=input('请输入负荷功率之和:Pd='); %输入 Beq Aeq=[1,1,1] %线性方程的系数约束矩阵 [P,F]=quadprog(H,f,[],[],Aeq,Beq,LB,UB,[],[]); disp('各发电机的负荷功率 P 为:'); disp(P); disp('3 台发电机的耗量总和 F 为:'); disp(F+10.22+10.32+9.38); %加上 a0 系数


相关文章:
电力系统运行方式分析和计算
电力系统运行方式分析和计算 设计报告 专班学姓 业: 级: 号: 名: 电气工程及其自动化 11 级电气 1 班 201130213166 杨玉豪 201130212176 潘鸣 华南理工大学电力...
北京交通大学电力系统运行方式及潮流分析实验报告
实验目的 1、掌握电力系统主接线电路的建立方法 2、掌握辐射形网络的潮流计算方法; 3、比较计算机潮流计算与手算潮流的差异; 4、掌握不同运行方式下潮流分布的特点...
电力系统运行方式及潮流分析实验
掌握电力系统主接线电路的建立方法 2、掌握辐射形网络的潮流计算方法 3、比较计算机潮流计算与手算潮流的差异 4、掌握不同运行方式下潮流分布的特点 二、 实验内容...
电力系统的最大运行方式和最小运行方式
在设计变、配电站选择开关电器和确定继电保护装置整定值时,往往 需要根据电力系统不同运行方式下的短路电流值来计算和校验所选用电器的稳定度和继电保护装置的灵敏...
《电力系统分析》_总结(概念、公式、计算)
电力系统分析》_总结(概念、公式、计算)_理学_高等教育_教育专区。自考 复习...10.在我国 110kv 以上电力系统的中性点运行方式一般为中性点直接接地方式。 11...
电力系统运行方式及潮流分析实验报告
及潮流分析实验 1、掌握电力系统主接线电路的建立方法 2、掌握辐射形网络的潮流计算方法; 3、比较计算机潮流计算与手算潮流的差异; 4、掌握不同运行方式下潮流分布...
电力系统 运行方式及潮流分析实验结题报告(实验一)
实验目的 1、掌握电力系统主接线电路的建立方法; 2、掌握辐射形网络的潮流计算方法; 3、比较计算机潮流计算与手算潮流的差异; 4、掌握不同运行方式下潮流分布的...
实验1 电力系统运行方式及潮流分析实验
电力系统运行方式及潮流分析实验一、 实验目的 电力系统规划、设计、运行分析离不...二、 实验内容 辐射形网络的潮流计算,分析负荷变化时线路两端电压的变化情况,...
电力系统分析上作业答案
电力系统经消弧线圈接地时,应采用的补偿方式为(C) ...电力系统潮流计算采用的数学模型是(A) 。 A、节点...电力系统静态稳定性是指(A ) A、正常运行的电力...
电力系统分析
第 11 章 电力系统简单不对称故障的分析和计算 1.电力系统简单不对称故障包括...出现了不正常断开或三相阻抗不相等的情况,这种不对称运行方式称为非全 相运行。...
更多相关标签:
电力系统运行方式 | 电力系统的运行方式 | 电力系统最大运行方式 | 电力系统运行方式 pdf | 电力系统运行方式实验 | 电力系统正常运行方式 | 电力系统最小运行方式 | 电力系统潮流计算 |