潮流计算课程设计

《电力系统分析》课程设计说明书课程设计任务书学生姓名：专业班级：指导教师：工作单位：题目:电力系统潮流分析初始条件：系统如图所示要求完成的主要任务:元件导纳参数为：5.275.0,48.0,35.0y132312jyjyj（1）、根据给定的运行条件，确定上图所示电力系统潮流计算式各节点的类型和待求量;（2）求节点导纳矩阵Y;（3）给出潮流方程或功率方程的表达式；（4）当用牛顿-拉夫逊法计算潮流时，给出修正方程和迭代收敛条件。时间安排：6月7日：熟悉设计任务；6月8日：收集相关资料；6月9日：设定设计原理；6月12日到15日：计算分析及结果分析；6月16日：写设计报告。指导教师签名：年月日系主任（或责任教师）签名：年月日321010.51.02+j123y13y12y《电力系统分析》课程设计说明书目录摘要………………………………………………………………………………11设计意义………………………………………………………………………22设计要求…………………………………………………………………………33设计环节…………………………………………………………………………43.1设计思路……………………………………………………………………43.2节点类型…………………………………………………………………43.3待求量……………………………………………………………………43.4导纳矩阵…………………………………………………………………43.5潮流方程…………………………………………………………………53.6牛顿—拉夫逊算法………………………………………………………63.6.1牛顿拉夫逊算法原理……………………………………………63.6.2修正方程…………………………………………………………73.6.3收敛条件…………………………………………………………93.6.4牛顿—拉夫逊法潮流计算程序框图……………………………94结果分析………………………………………………………………………115小结……………………………………………………………………………12参考文献…………………………………………………………………………13《电力系统分析》课程设计说明书摘要电力系统的出现使高效、无污染、使用方便易于调控的电能得到广泛应用，推动了社会生产各个领域的变化，开创了电力时代，开启了第二次科技革命。电力系统的规模和发展水平成为一个国家经济发展水平的标志之一。至今人类文明的主流发展方向依然与电力有着不可分割的联系。潮流计算是电力网络设计及运行中最基本的计算，对电力网络的各种设计方案及各种运行方式进行潮流计算，可以得到各种电网各节点的电压，并求得网络的潮流及网络中各元件的电力损耗，进而求得电能损耗。在数学上是多元非线性方程组的求解问题，求解的方法有很多种。牛顿—拉夫逊法是数学上解非线性方程式的有效方法，有较好的收敛性。将牛顿法用于潮流计算是以导纳矩阵为基础的，由于利用了导纳矩阵的对称性、稀疏性及节点编号顺序优化等技巧，使牛顿法在收敛性、占用内存、计算速度等方面都达到了一定的要求。关键词：电力系统潮流计算牛顿—拉夫逊算法《电力系统分析》课程设计说明书1.设计意义潮流计算是电力系统分析中的一种最基本的计算，也是最重要的计算。他的任务是对给定运行条件确定系统运行状态，如各母线上的电压（幅值及相角）、网络中的功率分布及功率损耗等。对于正在运行的电力系统，通过潮流计算可以判断电网电压母线、支路电流和功率是否越限，如果有越限，就应采取措施，调整运行方式。对于电力系统，进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。潮流计算还可以为继电保护和自动装置定整计算、电力系统故障计算和稳定计算等提供原始数据。潮流计算的结果是电力系统稳定计算和故障分析的基础。具体表现在以下方面：(1)在电网规划阶段,通过潮流计算,合理规划电源容量及接入点,合理规划网架,选择无功补偿方案,满足规划水平的大、小方式下潮流交换控制、调峰、调相、调压的要求。(2)在编制年运行方式时,在预计负荷增长及新设备投运基础上,选择典型方式进行潮流计算,发现电网中薄弱环节,供调度员日常调度控制参考,并对规划、基建部门提出改进网架结构,加快基建进度的建议。(3)正常检修及特殊运行方式下的潮流计算,用于日运行方式的编制,指导发电厂开机方式,有功、无功调整方案及负荷调整方案,满足线路、变压器热稳定要求及电压质量要求。(4)预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案。以上这些均是电力系统在实际工作中进行潮流计算的意义。同时潮流计算对于电力系统的检测也有积极的影响。综合各方面的用途，我们不难发现：潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。有着重要价值和广泛前景。《电力系统分析》课程设计说明书2设计要求如图所示的电网：元件导纳参数为：5.275.0,48.0,35.0y132312jyjyj1）根据给定的运行条件，确定图中电力系统潮流计算时各节点的类型、待求量；2）求节点导纳矩阵；3）给出潮流方程或功率方程的表达式；4）当用牛顿—拉夫逊法计算潮流时，给出修正方程和迭代收敛条件；系统如图1所示图1简单电力系统321010.51.02+j123y13y12y《电力系统分析》课程设计说明书3.设计环节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所确定的矩阵。支路导纳矩阵：表示一个电路中各支路导纳参数的矩阵。其行数和列数均为电路的支路总数。二端口导纳矩阵：对应于二端口网络方程，由二端口参数组成节点导纳矩阵：以导纳的形式描述电力网络节点注入电流和节点电压关系的《电力系统分析》课程设计说明书矩阵。它给出了电力网络连接关系和元件特性的全部信息，是潮流计算的基础方程式。本例应用结点导纳矩阵具体计算时，根据如下公式：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.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°1.25-j5.50.5-j30.75-j2.50.5-j31.3-j70.8-j40.75-j2.50.8-j41.55-j6.5YnjjijiiiVYVjQP1.*1112Uj2215.0jq《电力系统分析》课程设计说明书=(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同真解仍然有误差。为了进一步逼近真解，可以反复进行迭代计算，迭代计算通式是迭代过程的收敛判据为21kkxxf或式中，1和2为预先给定的小正数。牛顿-拉夫逊法实质上就是切线法，是一种逐步线性化的方法，此法不仅用于求单变量方程，也适用于多变量非线性代数方程的有效方法。0133qp0`00xfxfx0`00001xfxfxxxxkkkxfxfxx`1k《电力系统分析》课程设计说明书牛顿法至少是二阶收敛的，即牛顿法在单根附近至少是二阶收敛的，在重根附近是线性收敛的。牛顿法收敛很快，而且可求复根，缺点是对重根收敛较慢，要求函数的一阶3.6.2修正方程计算1、2节点的不平衡量iiiPVQ和、0VVV20222S022节点3是平衡节点，其电压iijfeiV是给定的，故不参加迭代。根据给定的容许误差510，按收敛判据kikikiVQP2,,max进行校验，以上节点1、2的不平衡量都未满足收敛条件，于是继续以下计算。修正方程式为VJW(n=3)T22211]VPQP[WTfefe2211V222222122122222212122121111121211111fVeVfVeVfPePfPePfQeQfQeQfPePfPePJ以上雅可比矩阵J中的各元素值是通过求偏导数获得的，对PQ节点来说，isisQP和是给定的，因而可以写出202GGG21312113131010101010101101101jjjjjjjjjjsseBfGffBeGePPPP101BBB1QQQQ1312113131010101010101101101jjjjjjjjjjsseBfGefBeGf《电力系统分析》课程设计说明书()()0()()0iijijiijjijjisjjjjijiijijijjjijjiisijjjijipfffeGeGePBBQQfffGeeGeBB对PV节点来说，给定量是isisVP和，因此可以列出2222()()0()0iisijijiijjijjjijjijiiisiifffeGeGePPBBfVVe当ji时,雅可比矩阵中非对角元素为22()0iiijiijijjiiijiijijjjjPQGeBfefPQBeGffeUUef