华中科技大学现代电力系统分析潮流计算作业文档格式.docx
《华中科技大学现代电力系统分析潮流计算作业文档格式.docx》由会员分享，可在线阅读，更多相关《华中科技大学现代电力系统分析潮流计算作业文档格式.docx（20页珍藏版）》请在冰豆网上搜索。
潮流计算分为离线计算和在线计算两大类。
离线计算可以用于电力系统调度，并确定系统的运行方式；
离线计算的结果还可以用于电力系统规划方案的分析以及优化系统的运行状态；
此外离线潮流计算可以作为初值，用于配合系统的故障分析以及稳定性分析。
在线计算可以用于能量管理系统EMS中，主要用作静态分析；
在线计算结果还可以帮助调度员合理安排系统潮流以及进行系统的安全性分析；
此外在线潮流计算还可以帮助优化系统的潮流及运行状态。
3.电力系统潮流计算常规方法
牛顿-拉夫逊法
电力系统潮流计算本质上是多元非线性方程组的求解问题，而求解的方法有很多种。
最常见的方法则是牛顿—拉夫逊法、高斯-塞德尔法、快速解耦法及直流潮流法。
牛顿—拉夫逊法能有效地解决非线性方程组，有较好的收敛性；
它利用了导纳矩阵的对称性、稀疏性以及节点标号顺序优化的技巧，已经成为电力系统潮流计算的最常用方法。
因此，本次潮流计算将重点了解牛顿—拉夫逊法。
牛顿—拉夫逊法实质上就是切线法，是一种逐步线性化的方法。
设有非线性方程组
其近似解为
，
，...，
。
设近似解与精确解分别相差
，…，
，则可得到
将上式中的多元函数在初值附近分别展开成泰勒级数。
则可得到
忽略上式中
的二次及以上阶次的各项，整理便可以得到
上式是对于修正量
的线性方程组，称为牛顿法的修正方程式。
利用高斯消去法或三角分解法可以求出修正量
然后对初始解近似解进行修正
经过k次迭代即可将上式写成如下形式
如此反复的迭代，直至满足收敛判据
或
将牛顿-拉夫逊法用于潮流计算，要求潮流方程写成形如方程式的形式，由于节点电压可以采用不同的坐标表示，牛顿—拉夫逊潮流计算也将相应地采用不同的计算公式。
节点电压用直角坐标表示时的牛顿—拉夫逊潮流计算
采用直角坐标系时，节点电压可表示为
导纳矩阵元素表示为
则可得到节点相应的有功功率和无功功率为
假定系统中的第1，2，…，m号节点为PQ节点；
假定系统中的第m+1，m+2，…，n-1号节点为PV节点，第i个节点的给定功率设为
和
，对节点可列写如下方程
由于第n号节点为平衡节点，其电压
是给定的，故不参加迭代。
式总共包含了2（n-1）个变量，而由可知，待求的变量也是2（n-1）个。
还注意到已经具备了方程组的形式。
因此即可得到如下的修正方程式
式中
而J则为雅克比矩阵。
有了这些，即可按照上面的牛顿—拉夫逊法计算电力系统的潮流。
由于极坐标下的牛顿—拉夫逊潮流计算的推导过程与上类似，这里就不在详尽介绍。
牛顿—拉夫逊法潮流计算程序框图
牛顿—拉夫逊法潮流计算程序框图如下图所示。
图3-1牛顿—拉夫逊法潮流计算程序框图
4.选用的潮流计算的系统
本次选用的潮流计算系统为美国西部电网WSCC三机九节点系统。
该系统比较经典，经常被业内人士用来作为分析的典型例子。
其接线图如下图所示
图4-1三机九节点系统连接图
统的基准值和频率取为：
系统中设备和线路的参数如表1~4所示：
表4-1发电机参数标幺值
参数名称
VN/kV
xd
x'
d
xq
q
xl
T'
d0
q0
H
D
G1
0s
G2
18
6s
G3
表2变压器参数的标幺值
变压器名称
T1
T2
T3
电抗X（pu）
变比K
1
表3线路参数的标幺值
线路名称
line1
line2
line3
line4
line5
line6
阻抗
Z（pu）
R
L
导纳B（pu）
表4负荷参数
负荷名称
LumpA
LumpB
LumpC
有功/MW
125
90
100
无功/Mvar
50
30
35
5.利用软件matpower计算潮流
matpower简介
Matpower更准确地讲，应该是matlab中的一个工具箱，而不是一个独立的软件。
它是由卡奈尔大学电气学院电力系统工程研究中心的RAY，CARLOS和甘德强在ROBERTTHOMAS的指导下开发出来的。
最初的基于MATLAB的电力潮流和最优潮流代码是为POWERWEB项目的需要而编写的。
Matpower是一个基于matlabm文件的组建包，用来解决电力潮流和优化潮流的问题。
它致力于为研究人员和教育从业者提供一种易于使用和可更新的仿真工具。
Matpower的设计理念是用尽可能简单、易懂，可更新的代码来实现最优秀的功能。
基态潮流计算
基态潮流计算条件说明及数据输入
本次潮流计算时选取Bus1为平衡节点，Bus2，Bus3为PV节点，其余的节点为PQ节点；
取发电机G2、G3的输出有功功率为其额定功率，则依据matpower用户手册即可输入以下数据（数据说明见注释）：
function[baseMVA,bus,gen,branch,areas,gencost]=case9hbyq
%定义case9hbyq为一个含变压器的三机九节点电力系统的潮流算例
version='
2'
;
%定义系统的基值为100MVA
baseMVA=100;
%bus_itypePdQdGsBsareaVmVabaseKVzoneVmaxVmin
%%母线参数，定义母线的序号、节点类型（PV、PQ、Swing节点）、注入的有功/无功、分区、
%%母线电压标幺值、基值、相角、电压标幺值的上下限等
bus=[
1300002101;
220000110181;
3200001101;
4100002102301;
519030002102301;
6100001102301;
7110035001102301;
8100001102301;
9112550002102301;
];
%busPgQgQmaxQminVgmBasestatusPmaxPminPc1Pc2Qc1min
%Qc1maxQc2minQc2maxramp_agcramp_10ramp_30ramp_qapf
%%电机参数，定义电机的有功无功出力，无功输出上下限、端口电压标幺值、电机基值等
gen=[
100300300110012501000000000000;
20300300110013001000000000000;
30300300110012701000000000000;
%%线路参数，定义两母线之间线路的阻抗和导纳值等
%branchdata
%fbustbusrxbrateArateBrateCratioanglestatusangminangmax
branch=[
1400250250250101-360360;
45250250250001-360360;
56150150150001-360360;
3600300300300101-360360;
67150150150001-360360;
78250250250001-360360;
8200250250250101-360360;
89250250250001-360360;
94250250250001-360360;
%%-----OPFData-----%%
%%areadata
%arearefbus
areas=[
15;
%%最优潮流计算所需的条件
%%generatorcostdata
%1startupshutdownnx1y1...xnyn
%2startupshutdownnc（n-1）...c0
gencost=[
21500035150;
2200003600;
23000031335;
潮流计算结果
将上述数据输入到文件中保存并运行后，在命令行窗口中输入runpf（'
case9hbyq'
）指令即可得到基态潮流计算的结果。
实际仿真结果见附录一。
表5-1母线的潮流计算结果
Bus
Voltage
Generation
Load
Mag（pu）
Ang（deg）
P（MW）
Q（MVAr）
*
-
2
3
4
5
6
7
8
9
Total：
表5-2线路的潮流计算结果
Branch
From
To
FromBusInjection
ToBusInjection
Loss（I^2*Z）
Total：
基态潮流结果计算分析
由基态潮流计算结果可以得到，该三机九节点系统的潮流计算采用的是牛顿-拉夫逊算法，总共迭代了4次、共计秒就收敛了，且各母线电压水平及有功无功功率在约束条件内。
其次，可知发电机总共发出的功率为
，而负荷总共消耗了
显然
造成有功功率不等的主要原因是输电线路存在电阻，部分有功消耗在了输电线路上，由潮流计算结果可得知
；
三台发电机的有功输出分别为、和。
三台发电机之所以有这样的有功输出，是因为的前提假设是与发电机1相连的母线为平衡节点，与发电机2、3相连的节点为PV节点，这意味着发电机节点的有功出力可变，而发电机2和3的有功出力不变。
在给定初始数据时，已经假定发电机2和3按照额定的有功功率出力和，由有功功率平衡可知，此时的发电机1的有功出力只能为。
造成无功功率不等的主要原因是变压器和输电线路会产生或消耗部分无功功率。
由基态潮流计算结果很容易验证发电机发出的无功功率加上输电线电纳产生的无功功率之和与线路上消耗的无功功率加上负荷消耗的无功功率之和。
当然，改变初始条件，如发电机出力及部分母线电压时，所算出的潮流会和此种情况有所不同，但是它们都有相同规律，那就是系统的有功功率和无功功率要平衡，且满足极限运行条件。
最优潮流计算
在运用matpower计算本次潮流计算之前，首先要清楚什么是最优潮流。
这里所谓的最优潮流的指标就是电力系统的经济运行，即在满足极限运行条件下，电力系统用最少的成本，发出能满足负荷需求的功率。
假定有n个火电厂，其燃料消耗特性分别为
系统的总负荷为
，暂不考虑网络中的功率损耗，假定各个发电厂的输出功率不受限制，则系统负荷在n个发电厂间的经济分配问题可以表示为：
在满足下式
的条件下，使目标函数
为最小。
这是多元函数求条件极值的问题。
由拉格朗日乘数法可知，当其满足以下条件时，在发出满足负荷功率需求的条件下，所消耗的燃料最少，即经济性最好。
当然，任一发电厂的有功功率和无功功率都不应该超出它的上、下限，各节点电压也必须维持在其允许的变化范围内，及
最优潮流计算条件说明及数据输入
和基态潮流计算一样，最优潮流计算时依旧选取Bus1为平衡节点，Bus2，Bus3为PV节点，其余的节点为PQ节点；
取发电机G2、G3的输出有功功率为其额定功率。
需要在此说明的是generatorcostdata是最优潮流计算的必须条件。
最优潮流的理论结果
由中输入的数据可得到三机九节点系统中三台发电机
的燃料特性分别为
其中有
结合、与即可求得
最优潮流仿真计算结果
在Matpower中，有专门的函数计算系统的最优潮流。
在命令行窗口中输入runopf（'
实际仿真结果见附录二。
表5-1优化后母线的潮流计算结果
表5-2优化后线路的潮流计算结果
最优潮流结果计算分析
除了可以得到与基态潮流计算结果类似的分析外，还可以得到以下分析结果。
线路损耗明显比基态潮流计算结果的小。
由最优潮流计算结果可以得到，三台发电机的有功功率输出分别为、和，与理论计算结果比较接近。
之所以会有些差别是因为输电线路存在电阻，而电阻上会消耗有功功率，因此实际功率平衡条件并没有得到充分的满足。
因此，理论计算结果和实际计算结果接近，但并不完全一致。
但这并不意味着理论计算结果毫无意义，当用软件计算电力系最优统潮流时，完全可以利用理论上的结果来判断潮流计算的结果是否合理或正确。
6.感想与小结
本次潮流计算让我收获颇多，我不仅学会了利用matpower计算电力系统潮流和最优潮流的方法，还懂得了与人多交流的重要性，更领悟到了一种学术研究的方式。
通过本次潮流计算，我了解到了电力系统潮流计算及其意义，在此基础上，了解电力系统潮流计算的模型以及常规的潮流计算的方法，掌握并熟练使用电力系统潮流计算软件Matpower。
达到了本次潮流计算的最初目的。
除此之外，我还懂得了与人交流的重要性。
在本次潮流计算之前，我完全不知道该如何使用matpower，因此不得不慢慢地详读matpower手册。
然而手册上的某些参数的说明并不完善，且无例子可参靠，因此学习进程相当缓慢。
后来在一次和师兄师姐的聊天中，知道了他们用过matpower，于是便向他们请教，最终解决了疑问，快速地完成了潮流计算。
从这次潮流计算中，我学会了与人交流，与人分享，交流和分享将会推动集体的进步，会让自己受益颇多。
最后，我还领悟到了一种学术研究的方式。
在开始做潮流计算之前，我便有个疑问，我们为什么要做潮流计算？
潮流计算该如何算？
因此，围绕这个问题，我又开始了自己的思考，并小有所获。
这种收获在我的报告中有所体现，我的报告的书写逻辑很明确，首先介绍什么是潮流计算，潮流计算有何意义，常用的计算方法有哪些，最后才过渡到利用matpower计算电力系统的潮流这一块。
这意味着在干一件事之前，我们首先要搞清楚做这件事有何意义，因为做一件没有意义的事情只会浪费我们的时间。
对于一个课题，我们认识到它是一个问题甚至比我们该如何去做更重要。
总而言之，我将好好利用这一次的收获，帮助我在未来科研的道路上越走越远，越走越顺。
参考文献
[1]何仰赞,温增银.电力系统分析:
上册[M].华中科技大学出版社,2002.
[2]丁晓莺,王锡凡.最优潮流在电力市场环境下最新发展[J].电力系统自动化,2002,26（13）:
1-7.
[3]徐恒娇,王洪诚,胡江航,等.基于Matpower的潮流计算方法[J].物联网技术,2013,3
（1）:
43-45.
[4]肖宇,董珉.应用Matlab计算电力系统的最优潮流[J].机械工程与自动化,2006（4）:
60-62.
附录一、基态潮流计算文件
>
runpf（'
）
MATPOWERVersion,14-Dec-2011--ACPowerFlow（Newton）
Newton'
smethodpowerflowconvergedin4iterations.
Convergedinseconds
===============================================================================
|SystemSummary|
Howmany?
Howmuch?
P（MW）Q（MVAr）
----------------------------------------------------------------------
Buses9TotalGenCapacityto
Generators3On-lineCapacityto
CommittedGens3Generation（actual）
Loads3Load
Fixed3Fixed
Dispatchable0Dispatchableof
Shunts0Shunt（inj）
Branches9Losses（I^2*Z）
Transformers3BranchCharging（inj）-
Inter-ties2TotalInter-tieFlow
Areas2
MinimumMaximum
---------------------------------------------------------
VoltageMagnitude.@bus9.@bus6
VoltageAngledeg@bus9deg@bus2
PLosses（I^2*R）-MW@line8-9
QLosses（I^2*X）-MVAr@line8-2
|BusData|
BusVoltageGenerationLoad
#Mag（pu）Ang（deg）P（MW）Q（MVAr）P（MW）Q（MVAr）
--------------------------------------------
升级会员 