电力系统课程设计用C语言进行复杂电力系统的潮流计算.doc

电力系统课程设计报告
班级:
学号:
姓名:
2013年3月 7日
用C语言进行复杂电力系统的潮流计算
电力系统的潮流计算是对电力系统分析的最基本步骤也是最重要的步骤,是指在一定的系统结构和运行条件下,确定系统运行状态的计算,也即是对各母线(节点)电压,各元件(支路)传输电线或功率的计算。通过计算出的节点电压和功率分布用以检查系统各元件是否过负荷,各点电压是否合理,以及功率损耗等。
即使对于一个简单的电力系统,潮流计算也不是一件简单就可以完成的事,其运算量很大,因此如果对于一个大的、复杂的电网来说的话,由于其节点多,分支杂,其计算量可想而知,人工对其计算也更是难上加难了。特别是在现实生活中,遇到一个电力系统不会像我们期望的那样可以知道它的首端电压和首端功率或者是末端电压和末端功率,而是只知道它的首端电压和末端功率,更是使计算变的头疼万分。为了使计算变的简单,我们就可以利用计算机,用C语言编程来实现牛顿-拉夫逊(Newton-Raphson)迭代法,最终实现对电力系统潮流的计算。

用牛顿-拉夫逊迭代法进行电力系统潮流计算的相关概念

如图所示的电力网络,将节点i和j的电压用i和j表示,它们之间的支路导纳表示为yij,那么有基尔霍夫电流定律可知注入接点I的电流i(设流入节点的电流为正)等于离开节点I的电流之和,因此有

j


I

(1-1)
(1-2)
如令则可将(1-2)改写为:
I=1,2,…,n. (1-3)
上式也可以写为: I=YU (1-4)
其中Y为节点导纳矩阵,也称为稀疏的对称矩阵,它是n×n阶方阵。对角元Yii称为自导纳,它等于与该节点I直接相连的所有支路导纳总和;非对角元Yij(i≠j)称为互导纳或转移导纳,它等于连结节点I,j支路导纳的负数,且有Yij=Yji,当节点I,j之间没有支路直接相连时,Yij=Yji=0。
电力系统的分析计算中,往往要作不同的运行方式下的潮流计算,如果系统发生变化,如投切一条线路或一台变压器,由于改变了一条支路的状态或参数只影响该支路两端节点的自导纳和他们之间的互导纳,因而对每一种运行方式不必重新形成导纳矩阵,只需对原有导纳矩阵作相应的修改即可。

在实际的电力系统中,已知的条件往往不是节点的注入电流而是负荷和发电机功率,且这些功率一般不随节点的电压变化而变化,而节点的电流则是随电压的变化而变化的,因此在已知节点导纳矩阵的情况下,必须用已知的节点功率来替代未知的节点注入电流,才能求出节点电压,每一个节点的注入功率方程式为:
G
I

(1-5)
(1-6)
(1-7)
节点注入电流用功率和电压表示为:
(1-8)
功率方程可以表示为:
(1-9)

对于有n个节点的电力网络,可以列出n个功率方程,由图可知一个节点有四个变量:注入有功功率Pi,注入无功功率Qi,节点电压幅值Ui和相角。n个节点有4n个变量,但只有2n个关系式,所以为了使潮流有确定解,必须给定其中2n个变量。根据给定节点变量的不同,可以有以下三种类型的节点:
PQ节点:给定注入功率Pi,Qi,即已知PGi,PLi,QGi,QLi,待求Ui,δi。例如:降压变电所母线(负荷节点),固定出力的发电厂母线。
PV节点:给定了注入有功功率Pi(PGi,PLi),Ui和QLi,待求QGi(Qi),δi。例如:有一定无功电源的降压变电所母线,有一定储备的发电厂母线。
平衡节点:给定了Ui,δi和PLi ,QLi,待求PGi ,QGi,即Pi,Qi,用来平衡全电网的功率,通常在一个独立的电力系统中只设一个平衡节点。
-拉夫逊迭代法
牛顿-拉夫逊迭代法将解非线性方程组的过程转化为反复求与之相对应的线性方程的求解过程。
对于一个n维非线性方程组: n=1,2,3,…,n
假定其初值为x1(0),x2(0),…,xn(0),也即其近似解,它与真值之间的误差为也即各变量与真解之间的修正量。
将这n个方程式都在初值的附近展开成Taylor级数且忽略二次项及高次项,则可得修正方程
, I=1,2,…n. (1-10)
将修正方程写成矩阵形式:
(1-11)
其中令J= ,称之为雅可比(Jacobi)方阵。
它的第I行,第j列交点的元素为第I个函数对第j个变量xj的偏导数在点(x1(0),x2(0),…,xn(0))的值,所以方程组是线性方程,可用于求出,从而得到新的近似解,
(1-12)
于是得到一般迭代式:
(1-13)
于是得到近似解:
(1-14)
迭代一直进行到Max{|yi-fi(x1(0),x2(0),…