电力系统潮流计算课程设计.docx
《电力系统潮流计算课程设计.docx》由会员分享，可在线阅读，更多相关《电力系统潮流计算课程设计.docx（58页珍藏版）》请在冰豆网上搜索。
电力系统潮流计算课程设计
一、问题重述
1.1课程设计要求
1、在读懂程序的基础上画出潮流计算基本流程图
2、通过输入数据，进行潮流计算输出结果
3、对不同的负荷变化，分析潮流分布，写出分析说明。
4、对不同的负荷变化，进行潮流的调节控制，并说明调节控制的方法，并
列表表示调节控制的参数变化。
5、打印利用DDRTS进行潮流分析绘制的系统图，以及潮流分布图。
1.2课程设计题目
1、系统图：
两个发电厂分别通过变压器和输电线路与四个变电所相连。
变电所1
变电所2
变电所4
变电所3
10kV母线
35kV母线
10kV母线
35kV母线
一次侧电压220kV
一次侧电压220kV
线路长为60km
线路长为80km
线路长为100km
线路长为80km
线路长为100km
线路长为80km
母线1
母线2
。
。
。
。
。
。
。
。
。
。
。
。
。
。
。
。
。
。
。
。
。
。
。
。
。
。
母线3
2*QFQ-50-2
2*QFS-50-2
2*TQN-100-2
2*TQN-100-2
电厂一
电厂二
2、发电厂资料：
母线1和2为发电厂高压母线，发电厂一总装机容量为（
400MW），母线3为
机压母线，机压母线上装机容量为（
100MW），最大负荷和最小负荷分别为50MW
和30MW；发电厂二总装机容量为（200MW）。
3、变电所资料：
（1）变电所1、2、3、4低压母线的电压等级分别为：
10KV35KV10KV35KV
（2）变电所的负荷分别为：
50MW40MW50MW60MW
1/33
（3）每个变电所的功率因数均为cosφ=0.85；
（4）变电所2和变电所4分别配有两台容量为75MVA的变压器，短路损耗414KW，短路电压（%）=16.7；变电所1和变电所3分别配有两台容量为63MVA的变压器，短路损耗为245KW，短路电压（%）=10.5；4、输电线路资料：
发电厂和变电所之间的输电线路的电压等级及长度标于图中，单位长度的电
阻为0.17，单位长度的电抗为0.402，单位长度的电纳为2.78*10-6S。
1.3课程设计基本内容
1.对给定的网络查找潮流计算所需的各元件等值参数，画出等值电路图。
2.输入各支路数据，各节点数据利用给定的程序进行在变电所在某一负荷情况下的潮流计算，并对计算结果进行分析。
3.跟随变电所负荷按一定比例发生变化，进行潮流计算分析。
1）4个变电所的负荷同时以2%的比例增大；
2）4个变电所的负荷同时以2%的比例下降
3）1和4号变电所的负荷同时以2%的比例下降，而2和3号变电所的负荷同时以2%的比例上升；
4.在不同的负荷情况下，分析潮流计算的结果，如果各母线电压不满足要
求，进行电压的调整。
（变电所低压母线电压10KV要求调整范围在
9.5-10.5之间；电压35KV要求调整范围在35-36之间）
5.轮流断开环网一回线，分析潮流的分布。
6.利用DDRTS软件，进行绘制系统图进行上述各种情况潮流的分析，并进行结果的比较。
7.最终形成课程设计成品说明书。
二、问题分析
2.1节点设置及分类
根据系统图可知此系统为两端供电网路，将母线1，2设为节点1，10，将变电所1、2、3、4的高低压侧分别设为节点2、3、4、5、6、7、8、9。
并且，将节点1设为平衡节点，将节点10设为PV节点，其余节点设为PQ节点。
2.2参数求取
设定基准值SB
100MVA，UB
SB
484根据题目
220KV，所以ZB
2
UB
原始资料，计算发电厂、变压器及线路的参数。
（1）运用下列公式计算变压器参数：
2
2
PkUN
RT
UK%UN
2
XT
1000SN
100S
N
2/33
ZTRT
XT
ZT
ZT
（2）计算线路参数
ZB
ZR
jX
（r
jx）L
（3）变电所负荷分别为：
变电所1
SL=50+j30.987
变电所2
SL=40+j27.79
变电所3
SL=50+j30.987
变电所4
SL=60+j37.18
将参数整理，见下表：
首端号
末端号
阻抗有名值
阻抗标幺值
电纳有名值
电纳标幺值
1
2
8.5+j20.1
0.0176+j0.0415
j0.000556
j0.2691
1
4
13.6+j32.16
0.0218+j0.0664
j0.0002224
j0.1076
1
6
13.6+j32.16
0.0218+j0.0664
j0.0002224
j0.1076
2
3
1.495+j40.335
0.0031+j0.0833
0
0
4
5
1.78+j53.885
0.037+j0.1113
0
0
4
6
10.2+j24.12
0.0211+j0.0498
j0.0001668
j0.0807
6
7
1.495+j40.335
0.0031+j0.0833
0
0
6
8
6.8+j16.08
0.0140+j0.0332
j0.0004448
j0.2153
8
9
1.78+j53.885
0.0037+j0.1113
0
0
8
10
8.5+j20.1
0.0176+j0.0415
j0.000556
j0.2691
（4）计算变压器分接头变比
变压器有5个抽头，电压调节范围为UN
2*2.5%，UN对应的分接头开始
时设变压器高压侧接主接头,降压变压器5个分接头时的非标准变比k*以备调压
时选用
k*1
220
（1
5%）
10
0.955
11
220
k*2
220
（1
2.5%）
10
0.932
11
220
k*3
220
10
0.909
11
220
k*4
220
（1
2.5%）
10
0.886
11
220
k*5
220
（1
5%）
10
0.864
11
220
对变电所低压母线为35KV时，非标准变比与10KV时相同。
2.3计算方法
利用牛顿拉夫逊法进行求解，用MATLAB软件编程，可以求解系统潮流分布根据题目的不同要求对参数进行调整，通过调节变压器变比和发电厂的电压，
3/33
求解出合理的潮流分布，最后用PSAT进行潮流分析，将两者进行比较。
2.4牛顿—拉夫逊法
1、牛顿—拉夫逊法概要
首先对一般的牛顿—拉夫逊法作一简单的说明。
已知一个变量X函数为：
f（X）0
（0）
到此方程时，由适当的近似值X出发，根据：
（n1）
f（X
X
X（n）
f（X
（n）
（n）
）（n1,2,......）
）
（n）
反复进行计算，当X满足适当的收敛条件就是上面方程的根。
这样的方法就是所谓的牛顿—拉夫逊法。
（n）
这一方法还可以做下面的解释，设第n次迭代得到的解语真值之差，即X的误差为时，则：
f（X（n）
）0
把f（X（n）
）在X
（n）
附近对
用泰勒级数展开
2
f（X（n）
）
f（X（n））
f（X（n））
f（X（n））......0
2!
上式省略去2以后部分
f（X（n））f（X（n））0
（n）
X的误差可以近似由上式计算出来。
f（Xf（X
（n）
（n）
）
）
（n）
比较两式，可以看出牛顿—拉夫逊法的休整量和X的误差的一次项相等。
用同样的方法考虑，给出n个变量的n个方程：
f1（X1,X2,,Xn）0
f2（X1,X2,,Xn）0
fn（X1,X2,,Xn）0
对其近似解X1得修正量X1可以通过解下边的方程来确定：
4/33
f1（X1,X2,,Xn）
f1
f1
f1
x1
x2
xn
X1
f2（X1,X2,
Xn）
f2
f2
f2
X2
x1
x2
xn
fn（X1,X2,,Xn）
fn
fn
fn
Xn
x1
x2
xn
式中等号右边的矩阵
fn都是对于X1,X2,
Xn的值。
这一矩阵称为雅可比
xn
（JACOBI）矩阵。
按上述得到的修正向量
X1,
X2,,Xn后，得到如下关系
XnXn
Xn
这比X1,X2,,Xn更接近真实值。
这一步在收敛到希望的值以前重复进行，
一般要反复计算满足
maxX1
n1
X1
n1,X2
n1
X2
n1,,Xnn1
Xnn1
为预先规定的小正数，
Xnn
1是第n次迭代Xn的近似值。
2、牛顿法的框图及求解过程
（1）用牛顿法计算潮流时，有以下的步骤：
①给这各节点电压初始值e（0）,f（0）；
②将以上电压初始值代入公式，求修正方程的常数项向量
P（0）,Q（0）,（V2）（0）；
③将电压初始值在带入上述公式，求出修正方程中系数矩阵的各元素。
④解修正方程式e（0）,f
（0）；
⑤修正各节点电压e
（1）
e（0）
e（0），f
（1）
f
（0）
f（0）；
⑥将e
（1），f
（1）在带入方程式，求出P
（1）,
Q
（1）,（V2）
（1）；
（k）（k）
⑦检验是否收敛，即maxPi,Qi
如果收敛，迭代到此结束，进一步计算各线路潮流和平衡节点功率，并打印输出结果。
如果不收敛，转回②进行下次迭代计算，直到收敛为止。
（2）程序框图如下：
5/33
启动
输入原始数据
形成节点导纳矩阵
分解各节点初始电压的实部和虚部
迭代次数K=0
（k）
（k）
（k）
（k）
求PQ节点的Pi,Qi
，求PV节点的Pi
，
Ui
置节点号i=0
是
雅克比矩阵是否形成，
i>n
否
求
得
雅
克
比
矩
阵
各
元
素
（k）
（k）
Jij
（k）
（k）
（k）
Sij
（k）
Hij
Nij
Lij
Rij
增大节点号
i=i+1
把雅克比矩阵单位化
求解修正方程，得ei（k），fi（k）
回带各电压新值，K=K+1
求解最大修正量|ei（k）|max，|fi（k）|max
否
是否收敛
是
计算输出电压大小及相角，节电功率及支路损耗
停止
三、问题求解
3.1等值电路的计算
电压是衡量电力系统电能质量的标准之一。
电压过高或过低，都将对人身及其用电设备产生重大的影响。
保证用户的电压接近额定值是电力系统调度的基本
6/33
任务之一。
当系统的电压偏离允许值时，电力系统必须应用电压调节技术调节系统电压的大小，使其维持在允许值范围内。
本文经过手算形成了等值电路图，并编写好了程序得出节点电压标幺值，使其满足所要求的调整范围。
我们首先对给定的程序输入部分作了简要的分析，程序开始需要我们确定输入节点数、支路数、平衡母线号、支路参数矩阵、节点参数矩阵。
（1）为了保证整个系统潮流计算的完整性，我们把凡具有母线及发电机处均选作节点，这样，可以确定有10个节点，节点号见等值电路图。
（2）确定完节点及编号后，各条支路也相应确定了，网络中总计有13条支路，我们对各支路参数进行了计算。
根据所给实际电路图和题中的已知条件，有以下公式计算各输电线路的阻抗和对地支路电容的标幺值和变压器的阻抗标幺
值。
选择电压基准值为UB
220KV和功率基准值SB100MVA
所以ZB
UB
484
SB
3.2依据题目要求及原始资料画出系统等值电路图：
3
S=50+j30.987
5
S=40+j27.79
7
S=50+j30.987
9
S=40+j24.8
1.495+j40.335
1.78+j53.885
1.495+j40.335
1.78+j53.885
1:
k
1:
k
1:
k
1:
k
1
2
3
4
2
4
6
8
10.2+j24.12
6.8+j16.08
j0.0000834
j0.0004448
8.5+j20.1
13.6+j32.16
8.5+j20.1
13.6+j32.16
j0.000556
j0.0001112
j0.000556
j0.0001112
1
10
U=231V
发电厂一
δ=0
P=200MW
发电厂二
U=231V
图1系统等值电路
利用psat软件可得计算所需系统图如图2所示。
7/33
图2
3.3潮流计算
3.3.1B1、B2矩阵的形成：
根据所求参数，以及
B1矩阵的含义，列写B1矩阵如下:
B1=[1
2
8.5+20.1i
0.000556i
1
0;
1
4
13.6+32.16i
0.0002224i
1
0;
1
6
13.6+32.16i
0.0002224i
1
0;
2
3
1.495+40.335i
0
1
1;
4
5
1.78+53.885i
0
1
1;
4
6
10.2+24.12i
0.0001668i
1
0;
6
7
1.495+40.335i
0
1
1;
6
8
6.8+16.08i
0.0004448i
1
0;
8
9
1.78+53.885i
0
1
1;
8
10
8.5+20.1i
0.000556i
1
0];
由各个变电所负荷功率可以计算出总功率为200MW，而发电厂一、二的总装机容量分别为400MW和200MW。
令发电厂二的功率为200MW，为了减小线路上的损耗，令发电机的电压为额定电压的1.05倍。
根据前面叙述的节点分类，形成B2矩阵如下：
B2=[0
0
231
231
0
1;
0
0
220
0
0
2;
0
50+30.987i
220
0
0
2;
8/33
0
0
220
0
0
2;
0
40+27.79i
220
0
0
2;
0
0
220
0
0
2;
0
50+30.987i
220
0
0
2;
0
0
220
0
0
2;
0
60+37.18i
220
0
0
2;
200
0
231
231
0
3];
3.3.2编写程序并运行
从结果中可以观察到各个节点电压标幺值分别为:
节点
1
2
3
4
5
6
7
8
9
10
电压
1.0500
1.0335
1.0055
1.0324
1.0190
1.0287
1.0006
1.0283
0.9817
1.0500
折算到有名值以后，可以观察到节点3、5、7、9的电压都与题目给定的范围相比偏高，因此调节变压器分接头和发电厂电压，通过对系统进行多次调整，最终得到合理结果，调整结果如下:
调节方法：
调节方法电厂1电压
电厂2电压
分接头1
分接头2分接头3
分接头4
未调整
1.05
1.05
1
1
1
1
调整后
1.0409
1.0409
1.05
1.025
1.025
1.025
调节结果：
节点号
1
2
3
4
5
6
7
8
9
10
未调整
1.0500
1.0335
1.0055
1.0324
1.0190
1.0287
1.0006
1.0283
0.9817
1.0500
调整后
1.0409
1.0240
0.9454
1.0159
0.9562
1.0161
0.9620
1.0170
0.9435
1.0409
由上表观察到，进行最后一次调节时节点
3，5，7，9的电压均在题目允许的范
围内，对线路损耗进行分析，统计每次调整后各个支路的有功损耗，记录于下表：
功率损耗
未调整
调整后
（1,2）支路
0.475067-28.0836i0.488386-27.5332i
（1,4）支路
0.217892-11.1554i0.425216-10.3807i
（1,6）支路
0.95644-9.3679i
1.1139-8.7545i
（2,3）支路
0.10571+2.8521i
0.11959+3.2266i
（4,5）支路
0.060308+1.8257i
0.095422+2.8887i
（4,6）支路
0.51037-7.3677i
0.45013-7.2695i
（6,7）支路
0.10676+2.8803i
0.11549+3.116i
（6,8）支路
2.74909-16.2733i
2.75176-15.7397i
（8,9）支路
0.19012+5.7553i
0.20583+6.2309i
（8,10）支路6.44415-13.8237i6.52583-13.0626i
总损耗
11.8159-72.7582i
12.2915-67.478i
由此可以看出，随着变压器分接头的调高有功损耗逐渐增大，符合实际情
况，而最终调节使得电压在规定范围内时，有功损耗增加了0.4756MW，相对与实际损耗12.29MW小很多，可以认为是合理的。
具体的潮流分布见下表
9/33
各条支路的首端功率Si各条支路的末端功率Sj
S（1,2）=50.608+6.68043iS（2,1）=-50.1196-34.2136i
S（1,4）=-2.65728+34.5734iS（4,1）=3.0825-44.9541i
S（1,6）=-35.6592+49.1553iS（6,1）=36.7731-57.9099i
S（2,3）=50.1196+34.2136iS（3,2）=-50-30.987i
S（4,5）=40.0954+30.6787iS（5,4）=-40-27.79i
S（4,6）=-43.1779+14.2754iS（6,4）=43.628-21.5448i
S（6,7）=50.1155+34.103iS（7,6）=-50-30.987i
S（6,8）=-130.5166+45.35175iS（8,6）=133.2683-61.09145i
S（8,9）=60.2058+43.4109iS（9,8）=-60-37.18i
S（8,10）=-193.4742+17.6806i
S（10,8）=200-30.74322i
各节点的功率S为（节点号从小到大排列）：
0.1229+0.9041i
0.0000-0.0000i-0.5000-0.3099i
0.0000-0.0000i
-0.4000-0.2779i
-0.0000-0.0000i-0.5000-0.3099i
-0.0000+0.0000i
-0.6000-0.3718i
2.0000-0.3074i
得到的各节点的电压标幺值，电压角度，节点注入无功，支路
升级会员 