1、电力系统分析（上）实验讲义第 21页 共 21 页实验一：节电导纳矩阵的形成一实验目的掌握节点导纳矩阵形成的方法二实验学时：2学时三实验原理与方法n 个独立节点的网络，n 个节点方程 。式中的即为节点导纳矩阵。 1自导纳 具体说，就等于与节点i相连的所有支路导纳的和。2互导纳 即给节点i加单位电压，其余节点全部接地，由节点j注入网络的电流。节点导纳矩阵的特点：（1） 直观易得阶数：等于除参考节点外的节点数n；对角元：等于该节点所连导纳的总和；非对角元Yij：等于连接节点i、j支路导纳的负值。（2） 稀疏矩阵，非对角元素中有大量的零元素。（3） 对称矩阵。3非标准变比变压器在包括变压器的输电线路
2、中，变压器线圈匝数比为标准变比时，变压器的高、低压两侧的电压和电流值用线圈匝数比来换算是不成问题的。但是变压器线圈匝数比为不等于标准变比时需要加以注意。ZTk:11212图中是按标准变比换算出来的变压器高、低压侧的电压和电流，理想变压器的线圈匝数比k：1表示变压器线圈匝数比对标准变比的比值。由图可得：上面的电压电流关系用形等值网络表示有两种： (a) (a) (b)对于用导纳表示的形等值网络，从1-1端口看进去的节点自导纳为：，和k等于1时相同。从2-2端口看进去的节点自导纳为：，是标准变比时导纳的 k2倍。互导纳, 是标准变比时导纳的 k倍。由以上可见，当有非标准变比变压器时，可按如下次序形
3、成节点导纳矩阵。（1） 先不考虑非标准变比（认为k=1），求导纳矩阵。（2） 再把接入非标准变比变压器的节点的自导纳加上，其中YT是从变压器相连接的另一端节点来看变压器的漏抗的倒数。（3） 由接入非标准变比变压器的对端节点来看自导纳不变。（4） 变压器两节点间的互导纳加上。4系统变更时的修正（1）从原有的节点上引出新的支路（输电线路或变压器），在这一支路另一端设新的节点。（2）在原有的支路上并联新的支路。（3）在没有支路直接相连的两个原有节点间附加新的支路。（4）原有变压器的变比或者分接头位置发生变化时。下面分别讨论这几种变更情况。（1） 从原节点i增加新的节点j和新的阻抗为z的支路时，节电导
4、纳矩阵的阶次增加一阶。自导纳和互导纳变化如下：（2） 在原有节点i和j间增加阻抗为z的新支路时，节点导纳矩阵的阶次不变，自导纳和互导纳分别变化为：（3）在没有支路直接相连的两个原有节点间附加新的支路。（4）变压器变比由k变成k时，用前面介绍的非标准变比变压器的处理方法，进行如下变化：5节点导纳矩阵的计算方法根据以上的讨论，导纳矩阵的计算归结如下： （1） 导纳矩阵的阶数等于电力系统网络的节点数。（2） 导纳矩阵各行非对角元素中非零元素的个数等于对应节点所连的不接地支路数。（3） 导纳矩阵的对角元素即自导纳等于相应节点所连支路的导纳之和。（4） 导纳矩阵非对角元素等于节点i,j之间导纳的负数（即
5、i,j之间阻抗倒数的负数）。当i,j之间有多条并联支路时，求时应求所有并联支路导纳的代数和的负数。四实验内容及步骤：内容：（1）形成图(a)所示网络的导纳矩阵。 （2）将图（a）中变压器的变比改为1.03，再重新求导纳矩阵。 （3）形成图(b)所示网络的导纳矩阵。 (a) (b)步骤：节点导纳矩阵形成的步骤：形成节点导纳矩阵的原始数据为了形成节点导纳矩阵，必须知道电力系统的接线图。网络接线由节点及连接两个节点的支路确定。实际上，只要输入了各支路两端的节点号，就相当于输入了系统的接线图。除了系统的接线图外，还要知道系统中各支路的阻抗值，如果线路对地有电容则需输入电纳。此外，对变压器支路还要知道他
6、的变比及变比在哪一侧。这样，一条支路一般需要输入6个数据，即i,j,z,bc,t,it,其中i,j是支路两端节点号，z为支路阻抗，bc为线路电纳，t为变压器支路的变比。it为折算到哪一侧的标志（如果支路的首端i处于高压侧则请输入“1”，否则请输入“0”）,在程序中用矩阵B来进行输入（其中矩阵的行数为支路数，列数为上述6个数据）。当支路为变压器支路时，t填实际变比，当支路为线路时，t为1，当支路为接地支路时，t为0。程序根据t是否为零作为区分接地支路与不接地支路的标志，或者把接地支路作为节点注入电流源的已知量来输入。矩阵X是由各节点的节点号与该节点的接地阻抗构成。形成节点导纳矩阵的程序框图及清单
7、。图(a)程序示例：程序运行结果：五实验报告要求（1）实验报告要求手写。（2）写清楚实验名称、目的。（3）简述计算机形成节点导纳矩阵基本原理。（4）按前述的实验内容要求，编程形成节点导纳矩阵，打印程序清单和计算结果。实验二：牛顿-拉夫逊法潮流计算一实验目的掌握的牛顿-拉夫逊法潮流计算的计算机方法。二实验学时：2学时三实验原理与方法1.牛顿拉夫逊迭代法非线性方程组： 其近似解为，与精确解相差有将上式按泰勒级数展开由此可得线性方程或修正方程为：其矩阵形式为：非线性代数方程的牛顿法迭代格式为：2直角坐标系下N-R法潮流计算：节点电压方程带入得：修正方程：其中：（i=1,2n,is）（i=1,2,m,
8、is） （i=m+1,m+2n,）雅可比矩阵各元素：非对角元素(ij)对角元素(i=j)四实验内容及步骤：内容：（1）求如1所示系统的潮流分布。 （2）求图2所示系统的潮流分布。（图中的元件参数以阻抗形式给出）。 图1图2步骤：N-R法潮流计算的步骤：需输入的原始数据（1）节点数n,支路数nl,平衡母线节点号isb（一般固定为1）,误差精度pr.(2) 由支路参数形成的矩阵B1矩阵B1的每行是由下列参数构成的：某支路的首端号P;某支路的末端号Q,且PQ;支路的阻抗（R+jX）;支路的对地容抗；支路的变比K；折算到哪一侧的标志（如果支路的首端P处于高压侧则请输入“1”，否则请输入“0”）。（3）
9、各节点参数形成的矩阵B2矩阵B2的每行是由下列参数构成的：节点所接发电机的功率SG；节点负荷的功率SL；节点电压的初始值；PV节点电压U的给定值；节点所接的无功补偿设备的电纳；节点分类标号igl。(igl=1为平衡节点；igl=2为PQ节点；igl=3为PV节点) （4）输入由节点号及其对地阻抗形成的矩阵X。程序框图、步骤及清单（1）程序流程图如下：（2）程序的步骤：输入原始数据和信息：y、Pis、Qis、Uis、约束条件；形成节点导纳矩阵YB；设置各节点电压初值ei(0), fi(0) ；将初始值代入，求不平衡量Pi(0), Qi(0), Ui2(0)；计算雅可比矩阵各元素（Hij（0）、L
10、ij（0）、Nij（0）、Jij（0）、Rij（0）、Sij（0）；解修正方程，求 ei（0）， fi（0）；求节点电压新值ei(1) =ei(0)ei(0), fi(1) = fi(0) fi(0)；判断是否收敛：Max| Pi(k) |, Max| Qi(k) |；若不收敛，重复迭代第4、5、6、7步，直到满足第8步的条件；求平衡节点的功率和PV节点的Qi及各支路的功率。（3）图1程序示例： 程序运行结果如下：五实验报告要求（1）试验报告要求手写。（2）写清楚实验名称、目的。（3）简述N-R法潮流计算的计算机算法的基本原理。（4）以前述给出的系统为对象，编程进行潮流计算，打印程序清单和计算结果。
