直角坐标系下牛顿法潮流计算.doc
- 文档编号:115868
- 上传时间:2022-10-03
- 格式:DOC
- 页数:27
- 大小:844.58KB
直角坐标系下牛顿法潮流计算.doc
《直角坐标系下牛顿法潮流计算.doc》由会员分享,可在线阅读,更多相关《直角坐标系下牛顿法潮流计算.doc(27页珍藏版)》请在冰豆网上搜索。
武汉理工大学《电力系统分析》课程设计说明书
直角坐标系下牛顿法潮流计算
1电力系统潮流计算
潮流计算是电力系统分析中的一种最基本的计算,它的任务是对给定的运行条件确定系统的运行状态,如母线上的电压(幅值及相角)、网络中的功率分布及功率损耗等。
在电力系统规划设计和现有电力系统运行方式的研究中,都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性.可靠性和经济性。
此外,电力系统潮流计算也是计算系统动态稳定和静态稳定的基础。
2节点导纳矩阵的形成
在图1(a)的简单电力系统中,若略去变压器的励磁功率和线路电容,负荷用阻抗表示,便可以得到一个有5个节点(包括零电位点)和7条支路的等值网络,如图1(b)所示。
将接于节点1和4的电势源和阻抗的串联组合变换成等值的电流源和导纳的并联组合,变得到图1(c)的等值网络,其中和分别称为节点1和4的注入电流源。
图1电力系统及其网络
以零电位点作为计算节点电压的参考点,根据基尔霍夫定律,可以写出4个独立节点的电流平衡方程如下:
(2-1)
上述方程组经过整理可以写成
(2-2)
式中,;;;;;;;。
一般的,对于有个独立节点的网络,可以列写个节点方程
(2-3)
也可以用矩阵写成
(2-4)
或缩写为
(2-5)
矩阵称为节点导纳矩阵。
它的对角线元素称为节点的自导纳,其值等于接于节点的所有支路导纳之和。
非对角线元素称为节点、间的互导纳,它等于直接接于节点、间的支路导纳的负值。
若节点、间不存在直接支路,则有。
由此可知节点导纳矩阵是一个稀疏的对称矩阵。
3牛顿-拉夫逊法潮流计算
3.1牛顿-拉夫逊法的基本原理
牛顿—拉夫逊法(Newton—Raphson法)是求解非线性方程代数方程组的有效迭代计算方法。
在牛顿—拉夫逊法的每一次迭代过程中,对非线性方程通过线性化处理逐步近似。
下面以单变量加以说明。
设有单变量非线性方程
(3-1)
求解此方程时。
先给出解的近似值它与真解的误差为,则将满足方程,即
(3-2)
将(3-8)式左边的函数在附近展成泰勒级数,于是便得
(3-3)
式中,……分别为函数在处的一阶导数,….,n阶导数。
如果差值很小,3-9式右端的二次及以上阶次的各项均可略去。
于是,3-9便简化为
=0(3-4)
这是对于变量的修正量的现行方程式,亦称修正方程式。
解此方程可得修正量
(3-5)
用所求的去修正近似解,变得
(3-6)
由于3-10是略去高次项的简化式,因此所解出的修正量也只是近似值。
修正后的近似解同真解仍然有误差。
但是,这样的迭代计算可以反复进行下去,迭代计算的通式是
(3-7)
迭代过程的收敛判据为
(3-8)
或(3-9)
式中,为预先给定的小正数。
这种解法的几何意义可以从图3-1得到说明。
函数y=f(x)为图中的曲线。
f(x)=0的解相当于曲线与x轴的交点。
如果第k次迭代中得到,则过点作一切线,此切线同x轴的交点便确定了下一个近似值。
由此可见,牛顿-拉夫逊法实质上就是切线法,是一种逐步线性化的方法。
应用牛顿法求解多变量非线性方程组3-1时,假定已给出各变量的初值,….,令,,…..分别为各变量的修正量,使其满足方程3-2即
(3-10)
将上式中的n个多元函数在初始值附近分别展成泰勒级数,并略去含有,,……,二次及以上阶次的各项,便得
.
(3-11)
方程式3-17也可以写成矩阵形式
(3-12)
方程式3-18是对于修正量,,……,的线性方程组,称为牛顿法的修正方程式.利用高斯消去法或三角分解法可以解出修正量,,……,。
然后对初始近似值进行修正
(i=1,2,….,n)(3-13)
如此反复迭代,在进行k+1次迭代时,从求解修正方程式
(3-14)
得到修正量,,,并对各变量进行修正
(i=1,2,…,n)(3-15)
式3-20和3-21也可以缩写为
(3-16)
和(3-17)
式中的X和分别是由n个变量和修正量组成的n维列向量;F(X)是由n个多元函数组成的n维列项量;J是n阶方阵,称为雅可比矩阵,它的第i、j个元素是第n个函数对第j个变量的偏导数;上角标(k)表示阵的每一个元素都在点处取值。
迭代过程一直到满足收敛判据
(3-18)或(3-19)
为止。
和为预先给定的小正数。
将牛顿-拉夫逊法用于潮流计算,要求将潮流方程写成形如方程式3-1的形式。
由于节点电压可以采用不同的坐标系表示,牛顿-拉夫逊法潮流计算也将相应的采用不同的计算公式。
3.2节点电压用直角坐标表示是的牛顿-拉夫逊法潮流计算
采用直角坐标时,节点电压可表示为
导纳矩阵元素则表示为
将上述表示式代入的右端,展开并分出实部和虚部,便得
(11-45)
假定系统中的第1,2,3···,m号节点为PQ节点,第i个节点的给定功率设为和,对对该节点可列写方程
(i=1,2,···,m)(11-46)
假定系统中的第m+1,m+2,···,n-1号节点为PV节点,则对其中每一个节点可以列写方程
(i=m+1,m+2,···,n-1)(11-47)
第n号节点为平衡点,其电压是给定的,故不参加迭代。
式(11-46)和式(11-47)总共包含了2(n-1)个方程,待求的变量有也是2(n-1)个。
我们还可看到,方程(11-46)和式(11-47)已经具备了方程组(11-34)的形式。
因此,不难写出如下的修正方程式
(11-48)
式中
上述方程中雅克比矩阵的各元素,可以对式(11-46)和式(11-47)求偏导数获得。
当时
(11-49)
当时
(11-50)
修正方程式(11-48)还可以写成分块矩阵的形式
(11-51)
式中,和都是二维列向量;是介方阵。
对于PQ节点
(11-52)
对于PV节点
(11-53)
从表达式(1-7)~(1-11)可以看到,雅克比矩阵有以下特点:
(1)雅克比矩阵各元素都是节点电压的函数,它们的数值将在迭代过程中不断的改变。
(2)雅克比矩阵的子块中的元素的表达式只用到导纳矩阵中的对应元素。
若,则必有。
因此,式(1-9)式中分块形式的雅克比矩阵同节点导纳矩阵一样稀疏,修正方程的求解同样可以用稀疏矩阵的求解技巧。
(3)无论在式(1-6)或式(1-9)中雅克比矩阵的元素或子块都不具有对称性。
用牛顿-拉夫逊法计算潮流的流程框图示于图11-27。
首先要输入网络的原始数据以及各节点的给定值并形成节点导纳矩阵。
输入节点电压初值和,置迭代计数k=0。
然后开始进入牛顿法的迭代过程。
在进行第k+1次迭代时,其计算步骤如下:
(1)按上一次迭代计算出的节点电压值和,利用式(11-46)和式(11-47)计算各类节点的不平衡量、和。
(2)按条件(11-43)校验收敛,即
<
如果收敛,迭代到此结束,转入计算各线路潮流和平衡节点的功率,并打印输出计算结果。
不收敛则继续计算。
(3)利用式(11-49)和式(11-50)计算雅克比矩阵的各元素。
(4)解修正方程式(11-48)求节点电压的修正量和。
(5)修正各节点的电压
(6)迭代计数加1,返回第一步继续迭代过程。
迭代结束后,还要计算出平衡节点的功率和网络中的功率分布。
输电线路功率的计算公式如下:
3.3牛顿-拉夫逊法潮流计算程序框图
否
是
输出
输入原始数据
计算平衡节点功率及全部线路功率
k=k+1
解修正方程式,求
计算雅可比矩阵各元素
是否收敛
计算及
k=0
给定节点电压初值
形成节点导纳矩阵
图2牛顿-拉夫逊法潮流计算程序框图
4C语言实现直角坐标表示的牛顿拉夫逊的计算
4.1计算节点导纳矩阵
该电力系统的等值电路如图所示
图3电力系统等值电路
根据节点方程和等值电路图,计算节点导纳矩阵的C语言子程序如下。
G[1][1]=z12r/(z12r*z12r+z12m*z12m)+k*k*z13r/(z13r*z13r+z13m*z13m)+z14r/(z14r*z14r+z14m*z14m);
B[1][1]=-z12m/(z12r*z12r+z12m*z12m)-k*k*z13m/(z13r*z13r+z13m*z13m)-z14m/(z14r*z14r+z14m*z14m)+y140+y120;
G[2][2]=z12r/(z12r*z12r+z12m*z12m)+z24r/(z24r*z24r+z24m*z24m);
B[2][2]=-z12m/(z12r*z12r+z12m*z12m
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 直角 坐标系 牛顿 潮流 计算