15号-潮流计算课程设计--

武汉理工大学《电力系统分析》课程设计说明书11.设计意义潮流计算是电力系统分析中的一种最基本的计算，也是最重要的计算。他的任务是对给定运行条件确定系统运行状态，如各母线上的电压（幅值及相角）、网络中的功率分布及功率损耗等。对于正在运行的电力系统，通过潮流计算可以判断电网电压母线、支路电流和功率是否越限，如果有越限，就应采取措施，调整运行方式。对于电力系统，进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。潮流计算还可以为继电保护和自动装置定整计算、电力系统故障计算和稳定计算等提供原始数据。潮流计算的结果是电力系统稳定计算和故障分析的基础。具体表现在以下方面：(1)在电网规划阶段,通过潮流计算,合理规划电源容量及接入点,合理规划网架,选择无功补偿方案,满足规划水平的大、小方式下潮流交换控制、调峰、调相、调压的要求。(2)在编制年运行方式时,在预计负荷增长及新设备投运基础上,选择典型方式进行潮流计算,发现电网中薄弱环节,供调度员日常调度控制参考,并对规划、基建部门提出改进网架结构,加快基建进度的建议。(3)正常检修及特殊运行方式下的潮流计算,用于日运行方式的编制,指导发电厂开机方式,有功、无功调整方案及负荷调整方案,满足线路、变压器热稳定要求及电压质量要求。(4)预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案。以上这些均是电力系统在实际工作中进行潮流计算的意义。同时潮流计算对于电力系统的检测也有积极的影响。综合各方面的用途，我们不难发现：潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。有着重要价值和广泛前景。武汉理工大学《电力系统分析》课程设计说明书22设计要求如图所示的电网：元件导纳参数为：5.275.0,48.0,35.0y132312jyjyj1）根据给定的运行条件，确定图中电力系统潮流计算时各节点的类型、待求量；2）求节点导纳矩阵；3）给出潮流方程或功率方程的表达式；4）当用牛顿—拉夫逊法计算潮流时，给出修正方程和迭代收敛条件；系统如图1所示图1简单电力系统321010.51.02+j123y13y12y武汉理工大学《电力系统分析》课程设计说明书33.设计环节3.1设计思路此题目为负载电力系统潮流计算模型。首先写出节点导纳矩阵，并分析各节点的类型，找出待求量。然后，确定潮流方程。最后进行潮流计算，而最后一步，可利用牛顿—拉夫逊法潮流分析。3.2节点类型电力系统潮流计算中，节点一般分为如下几种类型：PQ节点：节点注入的有功功率无功功率是已知的PV节点：节点注入的有功功率已知，节点电压幅值恒定，一般由无功储备比较充足的电厂和电站充当；平衡节点：节点的电压为1*exp(0°)，其注入的有功无功功率可以任意调节，一般由具有调频发电厂充当。对于本题目，节点分析如下：节点1给出有功功率为2.，无功功率为1，PQ节点。节点2给出有功功率为0.5，电压幅值为1.0，PV节点。节点3电压相位是0，电压幅值为1，平衡节点。3.3待求量节点1待求量是P，Q；节点2待求量是Q，；节点3待求量是U，。3.4导纳矩阵导纳矩阵分为节点导纳矩阵、结点导纳矩阵、支路导纳矩阵、二端口导纳矩阵。结点导纳矩阵：对于一个给定的电路(网络)，由其关联矩阵A与支路导纳矩阵Y所确定的矩阵。武汉理工大学《电力系统分析》课程设计说明书4支路导纳矩阵：表示一个电路中各支路导纳参数的矩阵。其行数和列数均为电路的支路总数。二端口导纳矩阵：对应于二端口网络方程，由二端口参数组成节点导纳矩阵：以导纳的形式描述电力网络节点注入电流和节点电压关系的矩阵。它给出了电力网络连接关系和元件特性的全部信息，是潮流计算的基础方程式。本例应用结点导纳矩阵具体计算时，根据如下公式：Yii=yi0+∑yijYik=-yik由题给出的导纳可求的节点导纳矩阵如下：5.525.1131211jyyy12210.53yyj13310.752.5yyj73.1232122jyyy23320.84yyj5.655.1323133jyyy进而节点导纳矩阵为：3.5潮流方程对n个节点的网络，电力系统的潮流方程一般形式是（i=1,2，…，n）其中Pi=PGi-PLdi,Qi=QGi-QLdi,即PQ分别为节点的有功功率无功功率。代入得潮流方程：1.25-j5.50.5-j30.75-j2.50.5-j31.3-j70.8-j40.75-j2.50.8-j41.55-j6.5YnjjijiiiVYVjQP1.*武汉理工大学《电力系统分析》课程设计说明书5=(1.25-j5.5)·11U+(0.5-j3)·21+(0.75-j2.5)01°=(0.5-j3)·11U+(1.3-j7)·21+(0.8-j4)·01°=(0.75-j2.5)·11U+(0.8-j4)·21+(1.55-j6.5)·01°3.6牛顿—拉夫逊算法3.6.1牛顿拉夫逊算法原理：牛顿法（NewtonMethodf(x)=0(Newton)法，就是将非线性方程线性化的一种方法。它是解代数方程和超越方程的有效方法之一。设有单变量非线性方程0xf，给出解的近似值0x，它与真解的误差为0x，则00xxx将满足0xf，即000xxf将上式左边的函数在0x附近展成泰勒级数，如果差值0x很小，0x二次及以上阶次的各项均可略去得：000`000xxfxfxxf这是对于变量的修正量0x的线性方程式，成为修正方程，解此方程可得修正量用所求得的0x去修正近似解，便得修正后的近似解1x同真解仍然有误差。为了进一步逼近真解，可以反复进行迭代计算，迭代计算通式是1112Uj2215.0jq0133qp0`00xfxfx0`00001xfxfxxxxkkkxfxfxx`1k武汉理工大学《电力系统分析》课程设计说明书6迭代过程的收敛判据为21kkxxf或式中，1和2为预先给定的小正数。牛顿-拉夫逊法实质上就是切线法，是一种逐步线性化的方法，此法不仅用于求单变量方程，也适用于多变量非线性代数方程的有效方法。牛顿法至少是二阶收敛的，即牛顿法在单根附近至少是二阶收敛的，在重根附近是线性收敛的。牛顿法收敛很快，而且可求复根，缺点是对重根收敛较慢，3.6.2修正方程计算1、2节点的不平衡量iiiPVQ和、0VVV20222S022节点3是平衡节点，其电压iijfeiV是给定的，故不参加迭代。根据给定的容许误差510，按收敛判据kikikiVQP2,,max进行校验，以上节点1、2的不平衡量都未满足收敛条件，于是继续以下计算。202GGG21312113131010101010101101101jjjjjjjjjjsseBfGffBeGePPPP101BBB1QQQQ1312113131010101010101101101jjjjjjjjjjsseBfGefBeGf武汉理工大学《电力系统分析》课程设计说明书7修正方程式为VJW(n=3)T22211]VPQP[WTfefe2211V222222122122222212122121111121211111fVeVfVeVfPePfPePfQeQfQeQfPePfPePJ以上雅可比矩阵J中的各元素值是通过求偏导数获得的，对PQ节点来说，isisQP和是给定的，因而可以写出()()0()()0iijijiijjijjisjjjjijiijijijjjijjiisijjjijipfffeGeGePBBQQfffGeeGeBB对PV节点来说，给定量是isisVP和，因此可以列出2222()()0()0iisijijiijjijjjijjijiiisiifffeGeGePPBBfVVe当ji时,雅可比矩阵中非对角元素为22()0iiijiijijjiiijiijijjjjPQGeBfefPQBeGffeUUef武汉理工大学《电力系统分析》课程设计说明书8当ji时,雅可比矩阵中对角元素为:111122()()()()22niijjijjiiiiiijiniijjijjiiiiiijjniijjijjiiiiiijiniijjijjiiiiiijjiijiiiPGeBfGeBfePGfBeGfBefQGfBeGfBeeQGeBfGeBffUeeUff代入数值后的修正方程为05.012020073.135.05.0325.15.535.05.525.12211fefe求解修正方程得1015.003611.02547.02211fefe计算节点电压的第一次近似值3.6.3收敛条件1015.01015.001013611.03611.007453.02547.01020212020212010111010111fffeeefffeee一轮迭代结束，根据收敛条件收敛判据kikikiVQP2,,max，若等式成立，结果收敛，迭代结束，计算平衡节点的功率和线路潮流计算，否则继续计算雅可比矩阵，解修正方程，直到满足收敛判据。3.6.4牛顿—拉夫逊法潮流计算程序框图具体框图如武汉理工大学《电力系统分析》课程设计说明书9图2牛顿—拉夫逊法潮流计算程序框图开始形成节点导纳矩阵输入原始数据设节点电压(0)(0)iief，i=1,2…,n,is置迭代次数0k置节点号i=1按式（3-3），（3-4）计算雅克比矩阵元素按式（3-2）计算PQ节点的()kiP，()kiQ，PV节点的()kiP，()2kiU求解修正方程式，得()kie，()kif雅克比矩阵是否已全部形成？计算平衡节点及PV节点功率求()max||ke，()max||kf迭代次数k=k+1i=i+1()()maxmax||,||kkef？潮流计算完成计算各节点电压的新值：(1)()()kkkieee(1)()()kkkifff武汉理工大学《电力系统分析》课程设计说明书104.结果分析给定节点电压初值0,0.1030201030201fffeee,四次笔算迭代过程后，得到节点电压和不平衡功率的变化情况分别于表1、2所示（510）：迭代计数k节点电压111jfeV222jfeV333jfeV10.7453-j0.36111-j0.1015120.4131-j0.35100.9901-j0.1479131.2973-j0.37971.008