C语言进行潮流计算.docxVIP

电力系统课程设计 C语言潮流计算 学院：电气工程 班级：电 092 班 学号： 0912002020 学生姓名： 闵 凯 2013.3.7 电力系统的潮流计算是对电力系统分析的最基本步骤也是最重要的步骤，是指在一 定的系统结构和运行条件下，确定系统运行状态的计算，也即是对各母线（节点）电压，各 元件（支路）传输电线或功率的计算。 通过计算出的节点电压和功率分布用以检查系统各元 件是否过负荷，各点电压是否合理，以及功率损耗等。 即使对于一个简单的电力系统，潮流计算也不是一件简单就可以完成的事， 其运算量很 大，因此如果对于一个大的、复杂的电网来说的话，由于其节点多，分支杂，其计算量可想 而知，人工对其计算也更是难上加难了。 特别是在现实生活中， 遇到一个电力系统不会像我 们期望的那样可以知道它的首端电压和首端功率或者是末端电压和末端功率， 而是只知道它 的首端电压和末端功率， 更是使计算变的头疼万分。 为了使计算变的简单， 我们就可以利用 计算机，用C语言编程来实现牛顿-拉夫逊（Newton-Raphson）迭代法，最终实现对电力系 统潮流的计算。 用牛顿-拉夫逊迭代法进行电力系统潮流计算的相关概念 1.节点导纳矩阵 如图所示的电力网络，将节点i和j的电压用U i和U j表示,它们之间的支路导纳表示为 yj,那么有基尔霍夫电流定律可知注入接点 I的电流I i（设流入节点的电流为正）等于离开节点 I的电流之和，因此有 yu上式也可以写为 yiJ(Ui-uJ)I i 八 I ij yu 上式也可以写为 yiJ(Ui-uJ) I i 八 I ij n -z j =0$ yj u =YH Ii :I = YU -yj =Yj则可将（1-2）改写为： YjUj I=1,2,…,n. (1-1) (1-2) (1-3) (1-4) 其中Y为节点导纳矩阵，也称为稀疏的对称矩阵，它是 nM阶方阵。对角元 Yu称为自 导纳，它等于与该节点 I直接相连的所有支路导纳总和；非对角元 Yj( i并)称为互导纳或 转移导纳，它等于连结节点 I, j支路导纳的负数，且有 Yj=Yji，当节点I, j之间没有支路 直接相连时，Yj=Yji=O。 电力系统的分析计算中，往往要作不同的运行方式下的潮流计算，如果系统发生变化， 如投切一条线路或一台变压器，由于改变了一条支路的状态或参数只影响该支路两端节点的 自导纳和他们之间的互导纳， 因而对每一种运行方式不必重新形成导纳矩阵， 只需对原有导 纳矩阵作相应的修改即可。 ?潮流计算的功率方程 在实际的电力系统中， 已知的条件往往不是节点的注入电流而是负荷和发电机功率， 且 这些功率一般不随节点的电压变化而变化， 而节点的电流则是随电压的变化而变化的， 因此 在已知节点导纳矩阵的情况下， 必须用已知的节点功率来替代未知的节点注入电流， 才能求 I i -為 Yij U I i -為 Yij U j j吕 (1-6) (1-7) 节点注入电流用功率和电压表示为： 「_ Si _(PGi -Pu)-j(QGi -Qu) _ P - jQi (18) 一氏一 Ui 一 l?i -功率方程可以表示为： P - iQ? n P 丛八 YjUj (1-9) 对于有n个节点的电力网络，可以列出 对于有n个节点的电力网络，可以列出 n个功率方程，由图可知一个节点有四个变量: 注入有功功率Pi，注入无功功率Qi，节点电压幅值 U i和相角：打。n个节点有4n 注入有功功率Pi，注入无功功率 只有2n个关系式，所以为了使潮流有确定解，必须给定其中2n个变量。根据给定节点变量的 不同，可以有以下三种类型的节点： PQ节点：给定注入功率 Pi, Qi,即已知PGi, PLi, QGi, Qu，待求Ui, $。例 女口：降压变电所母线(负荷节点)，固定出力的发电厂母线。 PV节点：给定了注入有功功率 Pi ( PGi , PLi), Ui和Qu，待求QGi (Qi), §。 例如：有一定无功电源的降压变电所母线，有一定储备的发电厂母线。 平衡节点：给定了 Ui，§和PLi，QLi，待求PGi，QGi，即Pi，Qi，用来平衡 全电网的功率，通常在一个独立的电力系统中只设一个平衡节点。 ?牛顿-拉夫逊迭代法 牛顿-拉夫逊迭代法将解非线性方程组的过程转化为反复求与之相对应的线性方程的 求解过程。 对于一个n维非线性方程组：fi (x1, x2,..., xn^ yi n=1,2,3,…,n 假定其初值为Xi(0),X2(0),…,Xn(0),也即其近似解，它与真值之间的误差为 .-■：x1(0)^■■:x2(0),...^■■:xn(0)也即各变量与真解之间的修正量。 将这n个方程式都在初值的附近展开成 Taylor级数且忽略二次项及高次项