文档介绍：该【电力系统稳态分析-牛顿拉夫逊法 】是由【木木在江边】上传分享，文档一共【31】页，该文档可以免费在线阅读，需要了解更多关于【电力系统稳态分析-牛顿拉夫逊法 】的内容，可以使用淘豆网的站内搜索功能，选择自己适合的文档，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此文档到您的设备，方便您编辑和打印。精品文档
精品文档
精品文档
0引言
潮流是配电网络分析的基础,用于电网调度、运行分析、操作模拟和设计规划,同时也是电压优化和网络接线变化所要参考的内容。潮流计算通过数值仿真的方法把电力系统的详细运行情况呈现给工作人员,从而便于研究系统在给定条件下的稳态运行特点。随着市场经济的发展,经济利益是企业十分看重的,而线损却是现阶段阻碍企业提高效益的一大因素。及时、准确的潮流计算结果,可以给出配电网的潮流分布、理论线损及其在网络中的分布,从而为配电网的安全经济运行提供参考。从数学的角度来看,牛顿—拉夫逊法能有效进行非线性代数方程组的计算且具有二次收敛的特点,具有收敛快、精度高的特点,在输电网中得到广泛应用。随着现代计算机技术的发展,利用编程和相关软件,可以更好、更快地实现配电网功能,本文就是结合牛顿-拉夫逊法的基本原理,利用C++程序进行潮流计算,计算结果表明该方法具有良好的收敛性、可靠性及正确性.
1牛顿-拉夫逊法基本介绍
对于N个节点的电力网络(地作为参考节点不包括在内),如果网络结构和元件参数已知,则网络方程可表示为:
(1—1)
式中,为N*N阶节点导纳矩阵;为N*1维节点电压列向量;为N*1维节点注入电流列向量。如果不计网络元件的非线性,也不考虑移相变压器,则为对称矩阵。
电力系统计算中,给定的运行变量是节点注入功率,而不是节点注入电流,这两者之间有如下关系:
(1—2)
精品文档
精品文档
精品文档
式中,为节点的注入复功率,是N*1维列矢量;为的共轭;是由节点电压的共轭组成的N*N阶对角线矩阵.
由(1-1)和(1—2),可得:
上式就是潮流方程的复数形式,是N维的非线性复数代数方程组。将其展开,有:
j=1,2,….,N(1—3)
式中,表示所有和相连的节点,包括。
将节点电压用极坐标表示,即令,代入式(1-3)中则有:
故有:
i=1,2,…,N(1—4)
式(1—4)是用极坐标表示的潮流方程.
而节点功率误差:
(1—5)
(1—6)
式中:,为节点给定的有功功率及无功功率。
1。2牛顿—拉夫逊法基本原理
精品文档
精品文档
精品文档
牛拉法是把非线性方程式的求解过程变成反复对相应的线性方程式的求解过程,即非线性问题通过线性化逐步近似,这就是牛拉法的核心。下面以非线性方程式的求解过程来进行说明。
设电力网络的节点功率方程一般形式如下:
(1-7)
式中,为节点注入功率给定值;为对应的物理量和节点电压之间的函数表达式;为节点电压。写成功率偏差的形式:
(1—8)
(1-8)作一阶泰勒展开:
定义为潮流方程的雅克比矩阵,为在处的值,则有:
,写成一般的表达式,有:
(1—9)
对于潮流收敛的情况,:
由简单迭代法收敛性分析的结论知,越接近解点,牛顿-拉夫逊法收敛越快,它具有二阶收敛速度。由图1。1可以直观地了解牛拉法的步骤:
精品文档
精品文档
精品文档
图1。1牛顿-拉夫逊法的几何解释
1。2。2极坐标的牛顿-拉夫逊法
在极坐标中,有如下的形式:
(1-10)
共2n-r个方程,状态变量为:
共2n-r个待求量。r个PV节点的电压幅值给定,不需求解。潮流雅克比矩阵的维数是(2n—r)*(2n—r),结构如下:
上式右侧的对电压幅值的偏导数项中的电压幅值的阶数减少了1,为使雅克比矩阵的各部分子矩阵具有一致的形式,在实际计算中,常将该项乘以电压幅值,并选取作为待求的修正量,则雅克比矩阵可写成:
精品文档
精品文档
精品文档
(1-11)
将式(1-10)和(1—11)代入式(1—9)的修正方程即可求得x的修正量,用它修正x直到为止。
将式(1-11)用下式表示:
其中每个字块的计算公式如下:
对角元素:(1—12)
非对角元素:(1—13)
2牛顿法潮流计算步骤
2。1程序流程图
在了解了牛拉法的原理之后,明确程序编写思路,如图2。1、2。2所示。其中图2。1中的“计算电压幅值和角度”步骤较多,单独用图2。2表示出来.
精品文档
精品文档
精品文档
当不符合收敛的条件“amontk>1”时,。
2。2计算步骤
下面讨论的是极坐标形式的牛顿法求解过程,大致分为以下几个步骤:
①形成节点导纳矩阵;
②给各节点电压设初值();
③根据式(1-12)、(1-13)生成雅克比矩阵(H、N、M、L);
④将节点电压初值代入式(1—5)、式(1-6),求出修正方程式的常数项向量;
⑤求解修正方程,得到电压幅值和角度;
⑥判断是否收敛,若收敛,计算平衡节点和线路功率;
精品文档
精品文档
精品文档
⑦输出结果,并结束.
3算例
,调用基于牛顿—拉夫逊法的C++程序。
其中节点4设为平衡节点,电压标幺值为1。05,计算误差为0。000001。
精品文档
精品文档
精品文档
对各个数字含义的解释如下:
网络模型有四个节点,四条支路,编号见图3。1。第一个零下面三行数为支路参数,分别表示三条支路的起始和终止节点编号,后面的为电阻、电抗和电纳,电导均为0,例如:120。10。。第二个零下面的为变压器支路,各数字意义同支路参数。接下去三行均为节点参数,分别表示注入有功功率和无功功率。
,得到运行结果,见图3。.
图3。3运行结果1
精品文档
精品文档
精品文档
将上述仿真结果整理为表格3。1、,其中“+”表示节点i输出功率给节点j,“-”表示节点j输出功率给i(纵向为i,横向为j).
节点号
1
2
3
4
1
0
+
—0。5
—
2
-0。24316
0
0
—0。312949
3
0
0
0
4
0。319671
0
0
表3。2节点无功功率输入与输出
节点号
1
2
3
4
1
0
—0。014708
-
-
2
0。0110505
0
0
—
3
0。097016
0
0
0
4
0。160255
0
0
根据表格计算:
节点1有功功率:0+--=-0。300582
无功功率:0-0。014708—-0。136187=-
精品文档
精品文档
精品文档
节点2有功功率:—+0+0—0。312949=—
无功功率:+0+0—0。14036=—0。1293095
节点3有功功率:0。5+0+0+0=0。5
无功功率:0。097016+0+0+0=
节点4有功功率:0。0482143++0+0=0。3678853
无功功率:++0+0=0。264895
根据已知条件,两个PQ节点的注入有功、无功分别为:
P1=0。3,Q1=0。18;
P2=0。55,Q2=0。13
潮流计算误差:
可见,误差均在允许范围内。
线路损耗:
3。4结论
通过上面的分析与计算,,线路损耗没能直接计算出来,而是需要手算,比较遗憾。程序在运行过程中,需要区分三种不同的节点,这由子程序保证实现。相比于快速分解法,牛拉法程序较为复杂,但更精确一点,潮流误差较小。