潮流计算课程设计

1、课程设计任务书学生姓名：专业班级：指导教师：工作单位：题目：电力系统潮流分析初始条件：系统如图所示元件导纳参数为：y12=0.5-j3,y23=0.8-j4,y13=0.75-j2.5(1)、根据给定的运行条件，确定上图所示电力系统潮流计算式各节点的类型和待求量；(2)求节点导纳矩阵 Y;(3)给出潮流方程或功率方程的表达式；(4)当用牛顿-拉夫逊法计算潮流时，给出修正方程和迭代收敛条件。时间安排：6 6 月 7 7 日：熟悉设计任务；6 6 月 8 8 日：收集相关资料；6 6 月 9 9 日设定设计原理；6 6 月 1212 日到 1515 日：计算分析及结果分析；6 6 月 1616 日

2、：写设计报告。指导教师签名：年月日系主任(或责任教师)签名:目录摘要 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摘要电力系统的出现使高效、无污染、使用方便易于调控的电能得到广泛应用,推动了社会生产各个领域的变化，开创了电力时代，开启了第二次科技革命。电力系统的规模和发展水平成为一个国家经济发展水平的标志之一。至今人类文明

3、的主流发展方向依然与电力有着不可分割的联系。潮流计算是电力网络设计及运行中最基本的计算， 对电力网络的各种设计方案及各种运行方式进行潮流计算，可以得到各种电网各节点的电压，并求得网络的潮流及网络中各元件的电力损耗，进而求得电能损耗。在数学上是多元非线性方程组的求解问题，求解的方法有很多种。牛顿一拉夫逊法是数学上解非线性方程式的有效方法，有较好的收敛性。将牛顿法用于潮流计算是以导纳矩阵为基础的，由于利用了导纳矩阵的对称性、稀疏性及节点编号顺序优化等技巧，使牛顿法在收敛性、占用内存、计算速度等方面都达到了一定的要求。关键词：电力系统潮流计算牛顿一拉夫逊算法1.设计意义潮流计算是电力系统分析中的一种

4、最基本的计算，也是最重要的计算。他的任务是对给定运行条件确定系统运行状态，如各母线上的电压(幅值及相角)、网络中的功率分布及功率损耗等。对于正在运行的电力系统，通过潮流计算可以判断电网电压母线、支路电流和功率是否越限，如果有越限，就应采取措施，调整运行方式。对于电力系统，进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。潮流计算还可以为继电保护和自动装置定整计算、电力系统故障计算和稳定计算等提供原始数据。潮流计算的结果是电力系统稳定计算和故障分析的基础。具体表现在以下方面：(1)在电网规划阶段，通过潮流计算，合理规划电源容量及接入点，合理规划网架，选择无功补偿方案，满足规划水平

5、的大、小方式下潮流交换控制、调峰、调相、调压的要求。(2)在编制年运行方式时，在预计负荷增长及新设备投运基础上，选择典型方式进行潮流计算，发现电网中薄弱环节，供调度员日常调度控制参考，并对规划、基建部门提出改进网架结构，加快基建进度的建议。(3)正常检修及特殊运行方式下的潮流计算，用于日运行方式的编制，指导发电厂开机方式，有功、无功调整方案及负荷调整方案，满足线路、变压器热稳定要求及电压质量要求。(4)预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案。以上这些均是电力系统在实际工作中进行潮流计算的意义。同时潮流计算对于电力系统的检测也有积极的影响。综合各方面的用途，我们不

6、难发现：潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。有着重要价值和广泛前景。2设计要求如图所示的电网：元件导纳参数为：yi2=0.5j3,y23=0.8j4,yi3=0.75j2.51）根据给定的运行条件，确定图中电力系统潮流计算时各节点的类型、待求量；2）求节点导纳矩阵；3）给出潮流方程或功率方程的表达式；4）当用牛顿一拉夫逊法计算潮流时，给出修正方程和迭代收敛条件；系统如图 1 所示图 1 简单电力系统3.设计环节1设计思路此题目为负载电力系统潮流计算模型。首先写出节点导纳矩阵，并分析各节点的类型，找出待求量。然后，确定潮流方程。最后进行潮流计算，而最后一步，可利用牛顿一

7、拉夫逊法潮流分析。1节点类型电力系统潮流计算中，节点一般分为如下几种类型：PQ 节点：节点注入的有功功率无功功率是已知的PV 节点：节点注入的有功功率已知，节点电压幅值恒定，一般由无功储备比较充足的电厂和电站充当；平衡节点：节点的电压为 1*exp(00),其注入的有功无功功率可以任意调节，一般由具有调频发电厂充当。对于本题目，节点分析如下：节点 1 给出有功功率为 2.,无功功率为 1,PQ 节点。节点 2 给出有功功率为 0.5,电压幅值为 1.0,PV 节点。节点 3 电压相位是 0,电压幅值为 1,平衡节点。1待求量节点 1 待求量是 P,Q;节点 2 待求量是 Q,6;节点 3 待求

8、量是 U,每。1导纳矩阵导纳矩阵分为节点导纳矩阵、结点导纳矩阵、支路导纳矩阵、二端口导纳矩阵。结点导纳矩阵：对于一个给定的电路(网络)，由其关联矩阵 A 与支路导纳矩阵 Y 所确定的矩阵。支路导纳矩阵：表示一个电路中各支路导纳参数的矩阵。其行数和列数均为电路的支路总数。二端口导纳矩阵：对应于二端口网络方程，由二端口参数组成节点导纳矩阵：以导纳的形式描述电力网络节点注入电流和节点电压关系的矩阵。它给出了电力网络连接关系和元件特性的全部信息，是潮流计算的基础方程式。本例应用结点导纳矩阵具体计算时，根据如下公式：Yi=yi0+Z2yijYik=-yik由题给出的导纳可求的节点导纳矩阵如下：yii=y

9、i2yi3=1.25-j5.5y12=y2i=0.5-j3yi3=y3i=0.75-j2.5y22-y2iy23=i.3J7y23-y32-0.8j4y33=y3iy32-i.55-j6.5进而节点导纳矩阵为：j5.50.5-j30.75-j2.5Y=0.5-j3i.3-j70.8-j4_0.75-j2.50.8-j4i.55-j6.51潮流方程对 n 个节点的网络，电力系统的潮流方程一般形式是n.八YjVjji其中 P=PGi-PLdi,Qi=QGi-QLdi,即 PQ 分别为节点的有功功率无功功率代入得潮流方程：2-jij-=(1.25-j5.5)U/6i+(0.5-j3)+(0.75-j

10、2.5)亿0Ui、105一io。-jq2=(0.5-j3)-U1cl+(1.3-j7)-12+(0.8-j4)-10Pi-jQi*Vi(i=i,2,，n)1,2p3-q3=(0.75-j2.5)-UiZ6i+(0.8-j4)1/&2+(1.55-j6.5)1/01.03.6牛顿一拉夫逊算法牛顿拉夫逊算法原理：牛顿法(NewtonMethod)：解非线性方程f(x)=0的牛顿(Newton)法,就是将非线性方程线性化的一种方法。它是解代数方程和超越方程的有效方法之设有单变量非线性方程f(x)=0,给出解的近似值 x(它与真解的误差为xf),则*=*严)+8胞满足f(x)=0,即f仅？)+

11、Axf)=0将上式左边的函数在 xf 谕近展成泰勒级数，如果差值 AxC 依小，AxC 匕次及以上阶次的各项均可略去得：f仪()+Ax?)=f(x(加f(x()3xC)=0这是对于变量的修正量Axf0粕线性方程式，成为修正方程，解此方程可得用所求得的AxC去修正近似解，使得x1=x。x0=x。一fx0修正后的近似解 x(H 真解仍然有误差。为了进一步逼近真解，可以反复进行迭代计算，迭代计算通式是fxk7-7迭代过程的收敛判据为|f(x(k)卜或|Ax(k卜42式中，工1和工2为预先给定的小正数。牛顿-拉夫逊法实质上就是切线法，是一种逐步线性化的方法，此法不仅用于求单变量方程，也适用于多变量非线

12、性代数方程的有效方法。牛顿法至少是二阶收敛的，即牛顿法在单根附近至少是二阶收敛的，在重根附近是线性收敛的。修正量xf)=-fx0fx0牛顿法收敛很快，而且可求复根，缺点是对重根收敛较慢，要求函数的一阶导数存在。修正方程计算 1、2 节点的不平衡量AR、AQi和 AVj33物(0)=以-睢)=%即(Gije&-Bijfj(0+fi(0f，)+Beg.jgj-.=-2-G11G12Gl3=一2一0二一2Qi0)=Q1s-Q10)=Q1s-I|f1(0fe6jej0)-BfS)-e:6jfj0)+3的()ILj1j1二1一B11B12B13)=1一0二-1V式)=|V2S-|V20=0节点

13、3 是平衡节点，其电压 Vi=e+jfi是给定的，故不参加迭代。根据给定的容许误差名=1。1 按收敛判据 max,Pi(k)AQi(k)qVi2(k)E 进行校验，以上节点 1、2 的不平衡量都未满足收敛条件，于是继续以下计算。修正方程式为W=-JV(n=3)_2_TW=RQ1卬2V2V-卜自f1:e2f2T毡 P扭R处PR1的5flce2cf213-Q1凶一Q1困*无242咨胆蟠.g%无2济2势V223；3V-V22困外出2讲2以上雅可比矩阵 J 中的各元素值是通过求偏导数获得的， 对 PQ 节点来说,pis和Qis是给定的，因而可以写出Pi=Pis-。工(Ge-Bufrf声(Gufj+Bu

14、e0jiWjWAQi=Qis-f2(Ge一Bijfj)十e,z(Gufj十Be)=0jiWjW对 PV 节点来说，给定量是pjDVis，因此可以列出&Pi=Pis-e(Ge-Buf,)f7(Guf上+8)=jiWjjiWj2222ri=Ve-(i+fJ=0,当 j#i 时，雅可比矩阵中非对角元素为fj於PcAQi=-9西Bjfi)wfjwe2“：U-Bijei-Gijfi当 j=i 时，雅可比矩阵中对角元素为：MPreif：P7rQi:e：Qin=一(Gijej-Bijfj)-Giiei-Biifij1n=，(GijfjBM-GifiBiiej=inL=(Gijfj，Bijej)-Gi

15、ifi，BiiSjin=-(Gij6-Bijfj)G，BiifijiUi2迫Ui2-fT=-2fi代入数值后的修正方程为-1.25-5.50.531AB1-21-5.51.253-0.53-10.53-1.3-70.5一00-20f2_110_求解修正方程得|-0.2547-0.3611I00.1015-计算节点电压的第一次近似值收敛条件裾=se=1-0.2547=0.7453fj=f10.-：f10=0-0.3611=-0.3611=02.62=1-0=1f21）=f20;f20）=0-0.1015=0.1015一轮迭代结束，根据收敛条件收敛判据max,Pi（k）AQi（k）AVi2（k）8

16、,若等式成立，结果收敛，迭代结束，计算平衡节点的功率和线路潮流计算，否则继续计算雅可比矩阵，解修正方程，直到满足收敛判据牛顿一拉夫逊法潮流计算程序框图具体框图如图 2 所示:.沿W2.计算平衡节点及 PVPV 节点功率图 2 2 牛顿一拉夫逊法潮流计算程序框图.结果分析给定节点电压初值 e/)=ef)=e3=1.0,f*)=f2f)=f)=0,经过四次笔算迭代过程后，得到节点电压和不平衡功率的变化情况分别于表 1 和表 2 所示(取w=10)：迭代计数 k节点电压V1+jf1V2=e2+jf2V3=e3+jf310.7453-j0.36111-j0.1015120.4131-j0.35100.

17、9901-j0.1479131.2973-j0.37971.0083-j0.0185140.8217-j0.36440.9986-j0.08801表 3.1 迭代过程中节点电压变化情况迭代计数 k节点/、平衡量PAQIP2V220-2-10.501-0.1482-0.9769-0.0726-0.01032-0.0902-0.6071-0.0480-0.00223-0.6272-4.3251-0.3610-0.01714-0.1816-1.2510-0.1042-0.0049表 3.2 迭代过程中节点不平衡量变化情况结果值与我的小组同学基本一样， 也在我的预期之内。 考虑到只用一种算法来算不全面

18、，我又尝试了 matlab 编程和 C 语言编程来计算待求量，得到了基本一致的结果。由此确定牛顿法具有很好的二次收敛性，是求解多元非线性方程的正确算法.小结这次的电力系统分析课程设计让我对平时所学的专业知识有了更深刻更具体的了解， 明白了理论知识必须与实践相结合才能更好的发挥作用。在不停的翻书上网查资料的过程中，我积累了大量的潮流计算以及电力系统的知识，全面透彻的了解了相关知识。同时在安装使用 matlab 和 VC 的过程中熟练掌握了相关的用法，为将来操作此类软件打下了坚实的基础。与人交流合作的能力也得到了提升。课设结束的时候， 我最大的感慨是学习这些知识的时候没有投入全部的精力， 以致课设刚开始时花费了大量时间补充基础知识。计算机的操作还需提升，英语也要加强。从这些事我下决心以后上课时必须认真听讲，不能虚度光阴。我坚信这次的课设为我的毕业论文和未来打下了很好的基础， 使我能更从容的面对以后生活和学习上的种种挑战。同学与老师的帮助也是必不可少的，在我遇到困难的时候，得到了来自同学的协助与老师的指