课设论文-(修复的)

电力系统潮流计算课程设计任务书（1）设计名称：电力系统潮流计算课程设计设计性质：理论计算，计算机仿真与验证计划学时：两周一、设计目的1.培养学生独立分析问题、解决问题的能力；2.培养学生的工程意识，灵活运用所学知识分析工程问题的能力3.编制程序或利用电力系统分析计算软件进行电力系统潮流分析。二、原始资料1、系统图：IEEE30节点。2、原始资料：见IEEE30节点标准数据库三、课程设计基本内容：1.采用PSAT仿真工具中的潮流计算软件计算系统潮流；1)熟悉PSAT仿真工具的功能；2)掌握IEEE标准数据格式内容；3)将IEEE标准数据转化为PSAT计算数据；2.分别采用NR法和PQ分解法计算潮流，观察NR法计算潮流中雅可比矩阵的变化情况，分析两种方法计算潮流的优缺点；3.分析系统潮流情况，包括电压幅值、相角，线路过载情况以及全网有功损耗情况。4.选择以下内容之一进行分析：1)找出系统中有功损耗最大的一条线路，给出减小该线路损耗的措施，比较各种措施的特点，并仿真验证；2)找出系统中电压最低的节点，给出调压措施，比较各种措施的特点，并仿真验证；3)找出系统中流过有功功率最大的一条线路，给出减小该线路有功功率的措施，比较各种措施的特点，并仿真验证；5.任选以下内容之一作为深入研究：（不做要求）1)找出系统中有功功率损耗最大的一条线路，改变发电机有功出力，分析对该线路有功功率损耗灵敏度最大的发电机有功功率，并进行有效调整，减小该线路的损耗；2)找出系统中有功功率损耗最大的一条线路，进行无功功率补偿，分析对该线路有功功率损耗灵敏度最大的负荷无功功率，并进行有效调整，减小该线路的损耗；3)找出系统中电压最低的节点，分析对该节点电压幅值灵敏度最大的发电机端电压，并有效调整发电机端电压，提高该节点电压水平；四、课程设计成品基本要求：1.绘制系统潮流图，潮流图应包括：1)系统网络参数2)节点电压幅值及相角3)线路和变压器的首末端有功功率和无功功率2.撰写设计报告，报告内容应包括以下几点：1)本次设计的目的和设计的任务；2)电力系统潮流计算的计算机方法原理，分析NR法和PQ分解法计算潮流的特点；3)对潮流计算结果进行分析，评价该潮流断面的运行方式安全性和经济性；4)找出系统中运行的薄弱环节，如电压较低点或负载较大线路，给出调整措施；5)分析各种调整措施的特点并比较它们之间的差异；6)结论部分以及设计心得；五、考核形式1.平时表现：3人一组，考查本组学生平时纪律、态度等；2.报告质量：设计成品的完成质量、撰写水平等；3.答辩考核：参照设计成品，对计算机方法进行电力系统潮流计算的相关问题等进行答辩；4.考核最终成绩由平时表现、报告成绩、答辩成绩组成。电力系统潮流计算摘要：电力系统潮流计算是电力系统分析的基础内容。通过对电力网络的潮流计算我们可以对系统的安全性，经济型，稳定性进行分析判断。并且能够准确的，快速的对电力系统中出现的各种问题进行调整。潮流计算是电力系统最基本，最常用的计算。Newton-Raphson法是数学上解非线性方程式的有效方法，有较好的收敛性。本次课程设计的主要内容就是基于PSAT的潮流计算和分析。针对问题一：通过题目中给的ieee格式的数据进行分析和理解，将其转换为PSAT标准格式的数据。然后使用PSAT利用N-R法对题目中的6机30节点的电力网络进行潮流计算。针对问题二：通过对PSAT中的N-R法的中断控制可以得到N-R法每次迭代的雅克比矩阵。通过对这些雅克比矩阵的分析得到雅克比矩阵的变化规律，即每次迭代出来的雅克比矩阵的个元素都在减小并且变化率也越来越小。利用P-Q法再次对题目中的数据进行潮流计算，得到的潮流的结果与N-R法求得的结果一样，计算的时间短于N-R法。然而，利用N-R法计算潮流其收敛性要好于P-Q分解法。针对问题三：利用PSAT中的simulink功能画出此6机30节点的系统。将目中的数据以PSAT标准格式填入此网络中。并对此网络进行仿真运行得到voltagemangnitudes图，voltageangles图，lineflows图。通过对这三个图的分析得到此网络的电压水平和潮流情况。针对问题四：通过对voltagemangntudes图的观察可以看出30节点的电压较低，需要提高30节点的电压。因此我采用了4种方法来提高30节点的电压。(1)提高8节点发电机的机端电压，提高30节点的电压。(2)提高27和28节点之间的变压器变比，提高30节点的电压。(3)在27节点和30节点之间串联电容，提高30节点的电压。(4)在母线10或24处并联电容器，提高30节点的电压。关键词：PSAT潮流计算牛顿-拉夫逊法P-Q分解法二．课程设计准备工作2.1牛顿-拉夫逊法原理及应用牛顿-拉夫逊法按照电压的表示方法不同，又分为直角坐标形式和极坐标形式，牛顿-拉夫逊法潮流计算具有二阶收敛特性，计算中收敛速度较快，但是当导纳矩阵阶数较高时，初值敏感性问题突出，如果初值选取的合适，则计算会非常快速准确。若初值选取不合适则计算有可能会不收敛。牛顿-拉夫逊法设计原理：设有非线性方程组1211)...(yxxxfn，，，；2212)...(yxxxfn，，，；……nnnyxxxf)...(21，，，；其近似解00201...nxxx，，，。近似解与精确解分别相差，，，nxxx...21则有如下：102021011...yxxxxxxfnn，，，；202021012...yxxxxxxfnn，，，；……nnnnyxxxxxxf0202101...，，，；对上述每一项式按泰勒级数展开忽略ix的高次方，可得如下修正方程组：nnnnnnnnnnnnxxxxfxfxfxfxfxfxfxfxfxxxfyxxxfyxxxfy...........................210020102022012010210110020100201220020111，，，，，，，，，；简写为：xJf；其中：J称为函数if的雅可比矩阵。求解功率方程步骤如下：iinjjjijijQPUYU1；代入ijijijjBGY及)sin(cosiiiijUU得：jQPeUjBGeUijinjjijijjiii)(1；将上式的实部与虚部分开列写如下：)sincos(1njjijijijijjiiBGUUP；)cossin(1njjijijijijjiiBGUUQ；iU及i分别为i节点的电压幅值及相角，iP及iQ分别为i节点注入的有功功率、无功功率。修正方程式如下：nPnppnppnnnnppppnpnpnpnpppUUUUHnnHHHNHNHNHNHJJHHJJHHLJLJNHNHLJLJNHNHPPQPQP2221112211221132221111222221212222212112121111121211112211有功功率不平衡量、无功功率不平衡量iP、iQ分别由下式计算：njjijijijijjiiiBGUUPP1sincosnjjijijijijjiiiBGUUQQ1cossin1、当ij时雅可比矩阵各个元素分别为：iiiiPH；iiiiQJ；iiiiiUUPN；iiiiiUUQL；2、当ij时雅可比矩阵各个元素为：jiijPH；jiijQJ；jjiijUUPN；jjiijUUQL；牛顿-拉夫逊法潮流计算步骤有以下几步：(1)形成节点导纳阵BY。(2)设各节点的电压的初值(0)i)0()0()0(、或、iiiUfe。(3)利用各节点的电压初值计算修正方程的不平衡量2)0()0()0(iiiUQP以及、。(4)利用各节点电压初值求出雅克比矩阵的各个元素。、以及、、、)0()0()0()0()0()0(ijijijijijijSRLJNH(5)解修正方程，求出各节点电压的变化量，即修正量(0)i)0()0()0(、或、iiiUfe。如果符合要求则跳出，如果不符合要求则继续下一步。(6)计算各节点电压的新值，即修正后的值。(7)御用各节点电压的新值自第三步开始下一次迭代。(8)计算平衡节点的功率和线路功率。牛顿-拉夫逊法计算框图：图一：牛顿-拉夫逊法潮流计算流程图输入原始数据形成节点导纳阵启动设节点电压初值，置迭代次数对PQ节点：计算kiP、kiQ对PV节点：计算kiP、2kiU置节点号i=1雅克比矩是否已全部形成计算雅克比矩阵元素增大节点号i=i+1否解修正方程式，由kiPkiQ2kiU和J求各节点电压变量kie、kif；i=1,2,…,n,i≠s是求出maxke、maxkf迭代是否收敛，maxke、maxkf？计算平衡节点功率和线路功率停止是计算各节点电压新值kikikieee1kikikifff1i=1,2,…,n;i≠s增大迭代次数，k=k+1否2.2P-Q分解法原理P-Q分解法是极坐标牛顿-拉夫逊法的一种简化算法快速分解法，有两个主要特点:(1)降阶在潮流计算的修正方程中利用了有功功率主要与节点电压相位有关,无功功率主要与节点电压幅值有关的特点,实现P-Q分解,使系数矩阵由原来的2N×2N阶降为N×N阶,N为系统的节点数(不包括缓冲节点)。(2)因子表固定化利用了线路两端电压相位差不大的假定,使修正方程系数矩阵元素变为常数,并且就是节点导纳的虚部。由于以上两个特点,使快速分解法每一次迭代的计算量比牛顿法大大减少。P-Q分解法只具有一次收敛性,因此要求的迭代次数比牛顿法多,但总体上快速分解法的计算速度仍比牛顿法快。快速分解法只适用于高压网的潮流计算,对中、低压网,因线路电阻与电抗的比值大,线路两端电压相位差不大的假定已不成立,用快速分解法计算,会出现不收敛问题。P-Q分解法的计算步骤与牛顿-拉夫逊法的计算步骤略有不同，在开始迭代之前就形成了系数矩阵'B和''B并得到了他们的逆阵。P-Q分解法计算步骤如下：(1)形成节点导纳矩阵BY。(2)形成系数矩阵'B,'B。(3)给定电压初值0iU，0i。(4)求功率不平衡量iiUP0。(5)接修正方程式得0。(6)求功率不平衡量iiUQ0。(7)解修正方程式得0jU。(8)修正电压初值00,iiU。(9)判断若满足条件则计算功率,若不满足条件则返回第3步重新求解。P-Q分解法的计算框图如下：启动输入原始数据形成节点导纳阵形成'B阵，并求其逆阵形成''B阵，并求其逆阵设节点电压：)0(iU和)0(i（i=1,2,3,…n）设置迭代次k=0置pk=0，qk=0计算)()(kikiUP，i=1,2,…,nis迭代是否收敛，max)()(kkUP？是置pk=0qk=0？是计算平衡节点功率sS~和线路功率ijS~计算)(ki，(i=1,2,…n,is)计算新的)(ki值，)1()()1(kikiki增大迭代次数kk+1置1qk置1pk增大迭代次数kk+1计算)()(kkUQ,(i=1,2,…m,is)迭代是否收敛，max)()(kkUQ？计算)(kiU，(