P-Q分解法潮流计算实验

xxxx实验报告学生姓名：学号：专业班级：实验类型：□验证□综合■设计□创新实验日期：2010.10.16实验成绩：一、实验项目名称P-Q分解法潮流计算实验二、实验目的与要求：目的：电力系统分析的潮流计算是电力系统分析的一个重要的部分。通过对电力系统潮流分布的分析和计算，可进一步对系统运行的安全性，经济性进行分析、评估，提出改进措施。电力系统潮流的计算和分析是电力系统运行和规划工作的基础。潮流计算是指对电力系统正常运行状况的分析和计算。通常需要已知系统参数和条件，给定一些初始条件，从而计算出系统运行的电压和功率等；潮流计算方法很多：高斯-塞德尔法、牛顿-拉夫逊法、P-Q分解法、直流潮流法，以及由高斯-塞德尔法、牛顿-拉夫逊法演变的各种潮流计算方法。本实验采用P-Q分解法进行电力系统分析的潮流计算程序的编制与调试，获得电力系统中各节点电压，为进一步进行电力系统分析作准备。通过实验教学加深学生对电力系统潮流计算原理的理解和计算，初步学会运用计算机知识解决电力系统的问题，掌握潮流计算的过程及其特点。熟悉各种常用应用软件，熟悉硬件设备的使用方法，加强编制调试计算机程序的能力，提高工程计算的能力，学习如何将理论知识和实际工程问题结合起来。要求：编制调试电力系统潮流计算的计算机程序。程序要求根据已知的电力网的数学模型（节点导纳矩阵）及各节点参数，完成该电力系统的潮流计算，要求计算出节点电压、功率等参数。三、主要仪器设备及耗材每组计算机1台、相关计算软件1套四、实验内容：1.理论分析：P-Q分解法是从改进和简化牛顿法潮流程序的基础上提出来的，它的基本思想是：把节点功率表示为电压向量的极坐标方程式，抓住主要矛盾，以有功功率误差作为修正电压向量角度的依据，以无功功率误差作为修正电压幅值的依据，把有功功率和无功功率迭代分开来进行。牛顿法潮流程序的核心是求解修正方程式，当节点功率方程式采取极坐标系统时，修正方程式为：VVLNJHQP/或展开为：VVLJQVVNHP//(4)以上方程式是从数学上推倒出来的，并没有考虑电力系统这个具体对象的特点。电力系统中有功功率主要与各节点电压向量的角度有关，无功功率则主要受各节点电压幅值的影响。大量运算经验也告诉我们，矩阵N及J中各元素的数值相对是很小的，因此对牛顿法的第一步简化就是把有功功率和无功功率分开来进行迭代，即将式(4)化简为：VVLQHP/(5)这样，由于我们把2n阶的线性方程组变成了二个n阶的线性方程组，因而计算量和内存方面都有改善。但是，H，L在迭代过程中仍然不断变化，而且又都是不对称矩阵。对牛顿法的第二个化简，也是比较关键的一个化简，即把式(5)中的系数矩阵简化为在迭代过程中不变的对称矩阵。众所周知，一般线路两端电压的相角差是不大的(通常不超过10~20度)，因此可以认为：cos1sinijijijijGB(6)此外，与系统各节点无功功率相应的导纳LiB必定远远小于该节点自导纳的虚部，即：2iLiiiiQBBV因此，2iiiiQVB(7)考虑到以上关系后，式(5)中系数矩阵中的元素表达式可以化简为：22iiiiiijijijiiiiiijijijHVBHVVBLVBLVVB(8)这样，式(5)中系数矩阵可以表示为：2111121211222212122221122nnnnnnnnnnnVBVVBVVBVVBVVBVBHLVVBVVBVB(9)进一步可以把它们表示为以下矩阵的乘积：111211122122120000nnnnnnnnBBBVVBBBHLVVBBB(10)将它代入(5)中，并利用乘法结合率，我们可以把修正方程式变为：1111111212222221221200nnnnnnnnnnVPVBBBVPBVBBVPVBBB(11)及11111121222221221200nnnnnnnnnVQVBBBVQBVBBVQVBBB(12)将以上两式的左右两侧用以下矩阵左乘11200nVVV=121/1/1/00nVVV就可得到11221111121222212212nnPVnPnVPnnnnnnVVBBBBVBBVBBB(13)及112211112122212212nnQVnQnVQnnnnnVVBBBBVBBVBBB(14)以上两式就是P-Q分解法达到修正方程式，其中系数矩阵只不过是系统导纳矩阵的虚部，因而是对称矩阵，而且在迭代过程中维持不变。它们与功率误差方程式)3,2,1(sincos1niBGVVPPijijijijnjjjiisi（15）)3,2,1(cossin1niBGVVQQijijijijnjjjiisi（16）构成了P-Q分解法迭代过程中基本计算公式，其迭代步骤大致是：(1)根据求得的Y矩阵形成有功迭代和无功迭代的简化雅可比矩阵```,BB。(2)给定各节点电压相角初值和各节点电压初值(0)(0),;iiV(2)根据(15)计算各节点有功功率误差iP，并求出/;iiPV(3)解修正方程式(13)，并进而计算各节点电压向量角度的修正量i(4)修正各节点电压向量角度i；)1()1()(kikiki(17)(5)根据式(16)计算各节点无功功率误差iQ，计算时电压相角用最新的修正值，并求出/;iiQV(6)解修正方程式(14)，求出各节点电压幅值的修正量iV(7)修正各节点电压幅值iV)1()1()(kikikiVVV(18)(8)返回(2)进行迭代，直到各节点功率误差及电压误差都满足收敛条件。2.理论数据：在上图所示的简单电力系统中，网络各元件参数的标么值如下：40.010.012jz,01528.0210120jyy,1.1,3.013kjz01920.0410140,50.012.014jyyjz01413.0,40.008.042024024jyyjz系统中节点1、2为PQ节点，节点3为PV节点，节点4为平衡节点，已给定P1s+jQ1s=-0.30-j0.18P2s+jQ2s=-0.55-j0.13P3s=0.5V3s=1.10V4s=1.05∠0°容许误差ε=10-5节点导纳矩阵：各节点电压：节点efvζ1.0.984637-0.0085960.984675-0.5001722.0.958690-0.1083870.964798-6.4503063.1.0924150.1289551.1000006.7323474.1.0500000.0000001.0500000.000000各节点功率：节点PQ1-0.300000-0.1800002–0.550000-0.13000030.500000-0.55130540.3678830.26469853.实验程序：n=input('pleaseentertheshortvaluen:');k=zeros(n,n);z=zeros(n,n);Y=zeros(n,n);yd=zeros(n,n);y=zeros(n,n);%其中yd(i,j)表明i,j结点间与i结点接地电阻，y(i,j)表明i,j结点间正常联接电阻。%输入数据z(1,2)=0.10+0.4*i;z(1,3)=0.3*i;z(1,4)=0.12+0.5*i;z(2,4)=0.08+0.4i;yd(1,2)=0.01528*i;yd(2,1)=0.01528*i;yd(1,4)=0.01920*i;yd(4,1)=0.01920*i;yd(2,4)=0.01413*i;yd(4,2)=0.01413*i;k(1,3)=1.1;%数据处理%form=1:nforj=1:nifz(m,j)~=0y(m,j)=1/z(m,j);y(j,m)=y(m,j);endendendform=1:nforj=1:nifk(m,j)~=0y(m,j)=k(m,j)/z(m,j);y(j,m)=y(m,j);yd(m,j)=(k(m,j)-1)*k(m,j)/z(m,j);yd(j,m)=(1-k(m,j))/z(m,j);endendendform=1:nforj=1:nifm==jY(m,j)=sum(y(m,:))+sum(yd(m,:));elseY(m,j)=-y(m,j);Y(j,m)=Y(m,j);endendendYA=[-0.3,-0.55,0.5,0;-0.18,-0.13,0,0;1,1,1.1,1.05;0,0,0,0];G=real(Y);B=imag(Y);B1=B([1,2,3],[1,2,3]);B2=B([1,2,],[1,2,]);fork1=0:100form=1:(n-1)sum=0;forj=1:nh=A(3,m)*A(3,j)*(G(m,j)*cos(2*pi/360*(A(4,m)-A(4,j)))+B(m,j)*sin(2*pi/360*(A(4,m)-A(4,j))));sum=sum+h;endop(1,m)=A(1,m)-sum;endV1=A([3],[1,2,3]);a=op./V1;a=a*inv(-B1)*180/pi;os=V1.\a;A([4],[1,2,3])=A([4],[1,2,3])+os;form=1:2sum=0;forj=1:nw=A(3,m)*A(3,j)*(G(m,j)*sin(2*pi/360*(A(4,m)-A(4,j)))-B(m,j)*cos(2*pi/360*(A(4,m)-A(4,j))));sum=sum+w;endoq(1,m)=A(2,m)-sum;endV2=A([3],[1,2]);b=oq./V2;b=b*inv(-B2);V2=V2+b;A([3],[1,2])=A([3],[1,2])+b;ifmax(max(abs(op)),max(abs(oq)))0.00001break;endendsum=0;sum1=0;sum2=0;forj=1:nx=A(3,4)*A(3,j)*(G(4,j)*cos(2*pi/360*(A(4,4)-A(4,j)))+B(4,j)*sin(2*pi/360*(A(4,4)-A(4,j))));sum=sum+x;c=A(3,4)*A(3,j)*(G(4,j)*sin(2*pi/360*(A(4,4)-A(4,j)))-B(4,j)*cos(2*pi/360*(A(4,4)-A(4,j))));sum1=sum1+c;d=A(3,3)*A(3,j)*(G(3,j)*sin(2*pi/360*(A(4,3)-A(4,j)))-B(3,j)*cos(2*pi/360*(A(4,3)-A(4,j))));sum2=sum2+d;endA(1,4)=sum;A(2,4)=sum1;A(2,3)=sum2;disp('PQVS');disp(A');4..组织调试结果：disp(Y):disp(A'):五、思考讨论题或体会或对改进实验的建议1．潮流计算有几种方法？简述各种算法的优缺点。答:潮流计算目前比较主要的方法有三种:高斯迭代法（高斯塞德尔法），牛顿拉夫逊法以及P-Q分解法。高斯迭代法是直接迭代，对初值要求比较低，程序简单，内