没有合适的资源?快使用搜索试试~ 我知道了~
电力系统稳态分析-牛顿拉夫逊法.pdf
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
0 下载量 112 浏览量
2022-10-27
23:59:25
上传
评论
收藏 665KB PDF 举报
温馨提示
试读
24页
电力系统稳态分析-牛顿拉夫逊法.pdf电力系统稳态分析-牛顿拉夫逊法.pdf
资源推荐
资源详情
资源评论
0 引言
潮流是配电网络分析的基础,用于电网调度、运行分析、操作模拟和设计
规划,同时也是电压优化和网络接线变化所要参考的内容.潮流计算通过数值仿
真的方法把电力系统的详细运行情况呈现给工作人员,从而便于研究系统在给
定条件下的稳态运行特点。随着市场经济的发展,经济利益是企业十分看重
的,而线损却是现阶段阻碍企业提高效益的一大因素。及时、准确的潮流计算
结果,可以给出配电网的潮流分布、理论线损及其在网络中的分布,从而为配
电网的安全经济运行提供参考。从数学的角度来看,牛顿—拉夫逊法能有效进行
非线性代数方程组的计算且具有二次收敛的特点,具有收敛快、精度高的特
点,在输电网中得到广泛应用.随着现代计算机技术的发展,利用编程和相关软
件,可以更好、更快地实现配电网功能,本文就是结合牛顿—拉夫逊法的基本原
理,利用 C++程序进行潮流计算,计算结果表明该方法具有良好的收敛性、可靠
性及正确性。
1 牛顿-拉夫逊法基本介绍
1.1 潮流方程
对于 N 个节点的电力网络(地作为参考节点不包括在内),如果网络结构
和元件参数已知,则网络方程可表示为:
(1-1)
式中,为 N*N 阶节点导纳矩阵;为 N*1 维节点电压列向量;为 N*1 维节点
注入电流列向量。如果不计网络元件的非线性,也不考虑移相变压器,则为对
称矩阵。
电力系统计算中,给定的运行变量是节点注入功率 ,而不是节点注入电流,
这两者之间有如下关系:
(1—2)
式中,为节点的注入复功率,是 N*1 维列矢量;为的共轭;是由节点电压的共
轭组成的 N*N 阶对角线矩阵.
由(1-1)和(1—2),可得:
1
上式就是潮流方程的复数形式,是 N 维的非线性复数代数方程组。将其展
开,有:
j=1,2,…。,N (1-3)
式中, 表示所有和相连的节点,包括。
将节点电压用极坐标表示,即令,代入式(1—3)中则有:
故有:
i=1,2,…,N (1-4)
式(1—4)是用极坐标表示的潮流方程。
而节点功率误差:
(1-5)
(1-6)
式中:,为节点给定的有功功率及无功功率。
1.2 牛顿-拉夫逊法基本原理
1。2。1 牛拉法的一般描述
牛拉法是把非线性方程式的求解过程变成反复对相应的线性方程式的求解
过程,即非线性问题通过线性化逐步近似,这就是牛拉法的核心。下面以非线
性方程式的求解过程来进行说明.
设电力网络的节点功率方程一般形式如下:
(1-7)
式中,为节点注入功率给定值;为对应的物理量和节点电压之间的函数表达
式;为节点电压。写成功率偏差的形式:
(1-8)
应用牛拉法求解如下。在给定的初值处将式(1—8)作一阶泰勒展开:
定义为潮流方程的雅克比矩阵,为在处的值,则有:
用修正就得到的新值。如果用 k 表示迭代次数,写成一般的表达式,有:
(1—9)
对于潮流收敛的情况,应比更接近于解点.收敛条件为:
由简单迭代法收敛性分析的结论知,越接近解点,牛顿-拉夫逊法收敛越快,
它具有二阶收敛速度。由图 1。1 可以直观地了解牛拉法的步骤:
图 1.1 牛顿-拉夫逊法的几何解释
2
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 表示出来.
图 2。1 牛顿法计算潮流的程序框图 图 2.2 电压幅值和角度求解步骤框图
当不符合收敛的条件“amontk〉1"时,即认为计算不收敛。具体程序见附录。
2。2 计算步骤
下面讨论的是极坐标形式的牛顿法求解过程,大致分为以下几个步骤:
① 形成节点导纳矩阵;
② 给各节点电压设初值();
③根据式(1-12)、(1-13)生成雅克比矩阵(H、N、M、L);
④ 将节点电压初值代入式( 1-5)、式(1-6),求出修正方程式的常数项向
量;
3
⑤ 求解修正方程,得到电压幅值和角度;
⑥判断是否收敛,若收敛,计算平衡节点和线路功率;
⑦输出结果,并结束。
3 算例
3.1 系统模型
本文以图 3.1 所示电力网络为例,调用基于牛顿—拉夫逊法的 C++程序。
图 3.1 系统模型
其中节点 4 设为平衡节点,电压标幺值为 1。05,计算误差为 0。
000001。
3。2 输入与输出
将图 3.1 所示模型的相关数据放在 data。dat 文件中
图 3.2 输入节点和支路数据
对各个数字含义的解释如下:
网络模型有四个节点,四条支路,编号见图 3。1.第一个零下面三行数为
支路参数,分别表示三条支路的起始和终止节点编号,后面的为电阻、电抗和电
纳,电导均为 0,例如:1 2 0。1 0.4 0。01528。第二个零下面的为变压器支
路,各数字意义同支路参数.接下去三行均为节点参数,分别表示注入有功功率
和无功功率。
调用 text。cpp 文件,得到运行结果,见图 3。3 和图 3。4。
图 3.3 运行结果 1
图 3。4 运行结果 2
3。3 结果分析
将上述仿真结果整理为表格 3.1、3。2,其中“+”表示节点 i 输出功率给
节点 j,“—”表示节点 j 输出功率给 i(纵向为 i,横向为 j)。
表 3.1 节点有功功率输入与输出
节点号
1
2
3
1
0
-0.24316
0。5
2
+0.245981
0
0
4
3
—0。5
0
0
4
-0.046563
-0.312949
0
4
节点号
1
2
3
4
根据表格计算:
0。0482143
1
0
0。0110505
0.097016
0.10464
0.319671
2
-0.014708
0
0
0.160255
0
3
-0。029001
0
0
0
0
4
—0。136187
—0.14036
0
0
表 3.2 节点无功功率输入与输出
节点 1 有功功率:0+0。245981—0。5—0.046563=—0。300582
无功功率:0-0.014708-0.029001-0.136187=-0.179896
节点 2 有功功率:—0。24316+0+0-0.312949=—0。556109
无功功率:0。0110505+0+0-0.14036=-0。1293095
节点 3 有功功率:0.5+0+0+0=0.5
无功功率:0.097016+0+0+0=0。097016
节点 4 有功功率:0。0482143+0。319671+0+0=0。3678853
无功功率:0。10464+0。160255+0+0=0。264895
根据已知条件,两个 PQ 节点的注入有功、无功分别为:
P1=0。3,Q1=0.18;
P2=0。55,Q2=0.13
潮流计算误差:
可见,误差均在允许范围内.
线路损耗:
3.4 结论
通过上面的分析与计算,验证了程序的正确性。由于编写过程的不足,线
路损耗没能直接计算出来,而是需要手算,比较遗憾。程序在运行过程中 ,需要
区分三种不同的节点,这由子程序保证实现 .相比于快速分解法,牛拉法程序较
为复杂,但更精确一点,潮流误差较小。
4 总结
本文基于牛顿-拉夫逊潮流算法的基本原理,利用 C++编程计算了一个 4 节点
简单电力网络的潮流,并验证了运行成果,误差在允许范围之内。因为牛拉法计
算过程中要不断生成新的雅各布矩阵 ,所以相对来说占用内存较多 ,但收敛速度
5
剩余23页未读,继续阅读
资源评论
G11176593
- 粉丝: 6642
- 资源: 3万+
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功