网络潮流计算课程设计

1广东工业大学华立学院课程设计（论文）课程名称题目名称__________________________________学生学部（系）专业班级学号学生姓名指导教师2011年6月9日目录摘要........................................................................................................................................21.设计目的与要求.................................................................................................................321.1设计目的.................................................................................................................31.2设计要求.................................................................................................................31.3设计题目.................................................................................................................42.题目解析.............................................................................................................................42.1设计思路.................................................................................................................52.2基本公式和变量分类.............................................................................................52.3电力系统潮流计算计算结果..................................................................................72.3.1节点导纳矩阵...................................................................................................72.3.2简化雅可比矩阵...............................................................................................82.3.3修正、迭代.........................................................................................................93.程序设计.............................................................................................................................93.1MATLAB编程说明................................................................................................93.2潮流计算流程图....................................................................................................103.3潮流计算源程序....................................................................................................104.小结...................................................................................................................................14参考文献..............................................................................................................................14摘要关键词：电力系统潮流计算牛顿—拉夫逊算法在如今的社会，电力已经成为人们必不可少的需求，而建立结构合理的大型3电力系统不仅便于电能生产与消费的集中管理、统一调度和分配，减少总装机容量，节省动力设施投资，且有利于地区能源资源的合理开发利用，更大限度地满足地区国民经济日益增长的用电需要。电力系统建设往往是国家及地区国民经济发展规划的重要组成部分。电力系统的出现，使高效、无污染、使用方便、易于调控的电能得到广泛应用，推动了社会生产各个领域的变化，开创了电力时代，发生了第二次技术革命。电力系统的规模和技术水准已成为一个国家经济发展水平的标志之一。电力系统稳态分析包括潮流计算（或潮流分析)和静态安全分析。潮流计算针对电力革统各正常运行方式，而静态安全分析则要研究各种运行方式下个别系统元件退出运行后系统的状况。其目的是校验系统是否能安全运行，即是否有过负荷的元件或电压过低的母线等。原则上讲，静态安全分析也可U用潮流计算来代替。但是一般静态安全分析需要校验的状态数非常多，用严格的潮流计算来分析这些状态往往计算量过大，因此不得不寻求一些特殊的算法以满足要求。潮流计算在数学上是多元非线性方程组的求解问题，求解的方法有很多种。牛顿—拉夫逊法是数学上解非线性方程式的有效方法，有较好的收敛性。将牛顿法用于潮流计算是以导纳矩阵为基础的，由于利用了导纳矩阵的对称性、稀疏性及节点编号顺序优化等技巧，使牛顿法在收敛性、占用内存、计算速度等方面都达到了一定的要求。41.设计目的与要求1.1设计目的1.掌握电力系统潮流计算的基本原理；2.潮流计算是电力系统分析中的一种最基本的计算，他的任务是对给定运行条件确定系统运行状态，如各母线上的电压（幅值及相角）、网络中的功率分布及功率损耗等。具体表现在以下方面：(1)在电网规划阶段,通过潮流计算,合理规划电源容量及接入点,合理规划网架,选择无功补偿方案,满足规划水平的大、小方式下潮流交换控制、调峰、调相、调压的要求。(2)在编制年运行方式时,在预计负荷增长及新设备投运基础上,选择典型方式进行潮流计算,发现电网中薄弱环节,供调度员日常调度控制参考,并对规划、基建部门提出改进网架结构,加快基建进度的建议。(3)正常检修及特殊运行方式下的潮流计算,用于日运行方式的编制,指导发电厂开机方式,有功、无功调整方案及负荷调整方案,满足线路、变压器热稳定要求及电压质量要求。(4)预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案。3.掌握并能熟练运用一门计算机语言（MATLAB语言或FORTRAN或C语言或C++语言）；4.采用计算机语言对潮流计算进行计算机编程计算。1.2设计要求1）根据给定的运行条件，确定图中电力系统潮流计算时各节点的类型、待求量；2）求节点导纳矩阵；3)赋予各节点电压变量初值后，求解不平衡量；4）形成雅克比矩阵；5）给出潮流方程或功率方程的表达式；56）当用牛顿—拉夫逊法计算潮流时，给出修正方程和迭代收敛条件；7）上机编程调试；8）计算分析给定系统潮流分析并与手工计算结果作比较分析。1.3设计题目2.题目解析62.1设计思路首先写出节点导纳矩阵，并分析各节点的类型，找出待求量。然后，确定潮流方程。最后进行潮流计算，而最后一步，可利用牛顿—拉夫逊法潮流分析。其次此电力系统是一个5节点，4支路的电力网络。其中包含3个PQ节点，一个PV节点，和一个平衡节点。由于牛顿拉夫逊法所需解的方程组最少，最后确定采用牛顿拉夫逊法。2.2基本公式和变量分类（1）.节点电压U和节点导纳矩阵Y。（2）.变量分类。在潮流问题中，任何复杂的电力网和电力系统都可以归结为以下元件（参数）组成。1）.发电机（注入电流或功率）。2）.负载（负的注入电流或功率）。3）.输电线支路（电抗、电阻）。4）.变压器支路（电阻、电抗、变化）。5）.变压器对地支路（导纳和感纳，本例中忽略）。6）.母线上的对地支路（阻抗或导纳，本例中忽略）。7）.线路上的对地支路（一般为线路电容导纳）。（3）.节点功率方程。有节点电压方程展开节点功率方程，此为电力系统的潮流方程的一般形式：1(1,2,)niijiiiiijjSPjQUIUYUi…,n1()(123niiiijjjiPjQIYUiU、、、...、n)（1-1）令iiiijijijUejfYGjB，将此式带入（1-1），将实部、虚部分开，得到直角坐标形式的节点功率方程：11()()0()()0jniiiijjijjiijjijjjniiiijjijjiijjijjjPPeGeBffGfBeQQfGeBfeGfBe（1-2）潮流方程具有的特点是：①他能表征电力系统稳态运行特性；②其为一组非线性方程，只能用迭代方法求其数值解；③方程中的电压U和导纳Y即可表示为直角坐标，又可表示为极坐标。因而潮流方程有多种表达方式——极坐标形式、直角坐标形式和混合坐标形式。7潮流计算的约束条件，即电压、相角和功率的约束条件。同时iU，ijY为电压和支路导纳jjjijijijUejfYGjB。但（1-2）是节点电压的线性方程组，iP，iQ分别为注入有功功率和无功功率。所以只要找出节点电压ie，jf这两个物理量，就能算出功率误差iP，iQ。对于PV节点，电压iU已知，则：222222()isiisiiUUUUef（1-3）所以对于n个节点的网络，其中（m-1）个PQ节点，一个不平衡点，其不平衡方程数目：0iP，i=1，2，…，n，则有（n-1）个方程。0iQ，i=1，2，…，n，则有（m-1）个方程。20iU，i=m+1，m+2，…，n，则有（n-m）个方程。(4)修正方程：形成雅可比矩阵。N-R法的思想是`()()0FxFxx；本例()PjQFx；对F(x)求偏导的式（1-2）、式（1-3），即式（1-2）、式（1-3）中的0iP、0iQ、U是多维变量的函数，对多维变量求偏导（iiPe、ijPe、iiPf、ijPf、iiQe、ijQe、iiPe、…），并以矩阵的形式表达称为雅可比矩阵。当j=i时，对角元素为8111122()()()()22niijjijjiiiiiiiijiniijjijjiiiiiiiijiniijjijjiiiiiiiijiniijjijjiiiiiiiiijiiiiiiiPGeBfGeBfNePGfBeBeGfHfQGfBeBeGfLeQGeBfGeBfJfUeeUf