matlab电力系统潮流计算

华中科技大学信息工程学院 课程设计报告书题目: 电力系统潮流计算 专 业： 电气工程及其自动化 班 级： 学 号： 学生姓名： 指引教师： 11 月 10 日 学 号 学生姓名 专业（班级）电气工程及其自动化（）设计题目电力系统潮流计算设计技术参数1. 变压器的阻抗，励磁损耗。2. 线路的阻抗，无功损耗。3. 变压器及线路的功率损耗4. 变压器及线路的电压降落5. 各节点的功率和电压设计要求1.计算过程涉及理论计算和MATLAB仿真两大部分。2.报告里面的理论分析（手工计算）过程必须要有具体推导过程，再把仿真过程放到相应的部分中去。工作量共24页；4000字工作计划第15周：对系统模型分析和资料收集，找到计算的思路，总结计算措施；第16周：对整个模型进行具体计算，并完毕word文档。参考资料1何仰赞, 温增银电力系统分析（第三版）M. 华中科技大学，2 刘同娟. MATLAB在电路分析中的应用.电气电子教学学报.3王守相，刘玉田 电力系统潮流计算研究现状-山东电力技术1996年05期指引教师签字 11月12日 学生姓名： 学号： 专业（班级）： 班 课程设计题目： 电力系统潮流计算 指引教师评语： 成绩： 指引教师： 年 月 日信息工程学院课程设计成绩评估表摘 要电力系统稳态分析涉及潮流计算和静态安全分析。本文重要运用的事潮流计算，潮流计算是电力网络设计与运营中最基本的运算，对电力网络的多种设计方案及多种运营方式进行潮流计算，可以得到多种电网各节点的电压，并求得网络的潮流及网络中的各元件的电力损耗，进而求得电能损耗。本位就是运用潮流计算具体分析，并有MATLAB仿真。核心词： 电力系统 潮流计算 MATLAB仿真AbstractElectric power system steady flow calculation and analysis of the static safety analysis. This paper, by means of the calculation, flow calculation is the trend of the power network design and operation of the most basic operations of electric power network, various design scheme and the operation ways to tide computation, can get all kinds of each node of the power grid voltage and seek the trend of the network and the network of the components of the power loss, and getting electric power. The standard is to use the power flow calculation and analysis, the specific have MATLAB simulation.Key words: Power system; Flow calculation; MATLAB simulation目 录1 任务提出与方案论证22 总体设计32.1潮流计算等值电路32.2建立电力系统模型32.3模型的调试与运营33 具体设计43.1 计算前提43.2手工计算74设计图及源程序114.1MATLAB仿真114.2潮流计算源程序115 总结19参照文献20 1 任务提出与方案论证潮流计算是在给定电力系统网络构造、参数和决定系统运营状态的边界条件的状况下拟定系统稳态运营状态的一种基本措施,是电力系统规划和运营中不可缺少的一种重要构成部分。可以说,它是电力系统分析中最基本、 最重要的计算,是系统安全、 经济分析和实时控制与调度的基本。常规潮流计算的任务是根据给定的运营条件和网路构造拟定整个系统的运营状态，如各母线上的电压（幅值及相角）、网络中的功率分布以及功率损耗等。潮流计算的成果是电力系统稳定计算和故障分析的基本。在电力系统运营方式和规划方案的研究中，都需要进行潮流计算以比较运营方式或规划供电方案的可行性、可靠性和经济性。同步，为了实时监控电力系统的运营状态，也需要进行大量而迅速的潮流计算。因此，潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。在系统规划设计和安排系统的运营方式时，采用离线潮流计算；在电力系统运营状态的实时监控中，则采用在线潮流计算。是电力系统研究人员长期研究的一种课题。它既是对电力系统规划设计和运营方式的合理性、 可靠性及经济性进行定量分析的根据，又是电力系统静态和暂态稳定计算的基本。潮流计算经历了一种由手工到应用数字电子计算机的发展过程，目前的潮流算法都以计算机的应用为前提用计算机进行潮流计算重要环节在于编制计算机程序，这是一项非常复杂的工作。对系统进行潮流分析，本文运用 MATLAB中的Simpower Systems工具箱设计电力系统，在simulink 环境下，不仅可以仿真系统的动态过程，还可以对系统进行稳态潮流分析。2 总体设计SimpowerSystems使用Simulink环境，可以将该系统中的发电机、变压器，线路等模型联结起来，形成电力系统仿真模拟图。在加人测量模块，并对各元件的参数进行设立后，用measurement和sink中的仪器可以观测各元件的电压、电流、功率的大小。2.1潮流计算等值电路 2.2建立电力系统模型在Simulink中按照电力系统原型选择元件进行建模。所建立的模型和建立的措施在具体设计中详述。在电力系统模型的建立工程中重要波及到的是：元器件的选择及其参数的设立；发电机选型；变压器选择；线路的选择；负荷模型的选择；母线选择。 2.3模型的调试与运营建立系统模型，并设立好参数后来，就可以在Simulink环境下进行仿真运营。运营的具体成果和分析也在具体设计中详述。3 具体设计3.1 计算前提一方面是发电机的参数计算，先对5个发电厂简化为5台发电机来计算。发电机G1：发电机G2：发电机G3：发电机G4：发电机G5：另一方面是变电站的参数计算，我们还是对7个变电站简化为7台变压器来计算。变压器T1：变压器T2：（双并联）变压器T3：（四并联）变压器T4：（双并联）变压器T5：变压器T6：（两个三绕组变压器并联）变压器T7：（双并联）再次是传播线参数计算，5条传播线的具体计算如下。根据教材查得线路L1：线路L2：线路L3：（双回路）线路L4：线路L5：（双回路）3.2手工计算FLR1：FLR2：FLR3：FLR4：FLR5:计算每一种FLR的功率分布和电压分布计算如下：FLR1：FLR2：功率分布： 电压分布：FLR3：功率分布： 电压分布：FLR4：功率分布： 电压分布：FLR5:这里我们先将f点和发电机G5当做电源，通过和构成两端供电网络以g点作为运算负荷进行计算。电压分布：4设计图及源程序4.1MATLAB仿真有关的原始数据输入格式如下：1、B1是支路参数矩阵，第一列和第二列是节点编号。节点编号由小到大编写。2、对于具有变压器的支路，第一列为低压侧节点编号，第二列为高压侧节点编号，将变压器的串联阻抗置于低压侧解决，第三列为支路的串列阻抗参数，第四列为支路的对地导纳参数，第五烈为含变压器支路的变压器的变比，第六列为变压器与否与否具有变压器的参数，其中“1”为具有变压器，“0”为不具有变压器。3、B2为节点参数矩阵，其中第一列为节点注入发电功率参数；第二列为节点负荷功率参数；第三列为节点电压参数；第六列为节点类型参数，其中“1”为平衡节点，“2”为PQ节点，“3”为PV节点参数。4、X为节点号和对地参数矩阵。其中第一列为节点编号，第二列为节点对地参数。4.2潮流计算源程序close allclear alln=input(请输入节点数:n=);n1=input(请输入支路数:n1=);isb=input(请输入平衡节点号:isb=);pr=input(请输入误差精度:pr=);B1=input(请输入支路参数:B1=);B2=input(请输入节点参数:B2=);X=input(节点号和对地参数:X=);Y=zeros(n);Times=1; %置迭代次数为初始值%创立节点导纳矩阵for i=1:n1 if B1(i,6)=0 %不含变压器的支路 p=B1(i,1); q=B1(i,2); Y(p,q)=Y(p,q)-1/B1(i,3); Y(q,p)=Y(p,q); Y(p,p)=Y(p,p)+1/B1(i,3)+0.5*B1(i,4); Y(q,q)=Y(q,q)+1/B1(i,3)+0.5*B1(i,4); else %具有变压器的支路 p=B1(i,1); q=B1(i,2); Y(p,q)=Y(p,q)-1/(B1(i,3)*B1(i,5); Y(q,p)=Y(p,q); Y(p,p)=Y(p,p)+1/B1(i,3); Y(q,q)=Y(q,q)+1/(B1(i,5)2*B1(i,3); endendYOrgS=zeros(2*n-2,1);DetaS=zeros(2*n-2,1); %将OrgS、DetaS初始化%创立OrgS，用于存储初始功率参数h=0;j=0;for i=1:n %对PQ节点的解决 if i=isb&B2(i,6)=2 h=h+1; for j=1:n OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3)*(real(Y(i,j)*real(B2(j,3)-imag(Y(i,j)*imag(B2(j,3)+imag(B2(i,3)*(real(Y(i,j)*imag(B2(j,3)+imag(Y(i,j)*real(B2(j,3); OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3)*(real(Y(i,j)*real(B2(j,3)-imag(Y(i,j)*imag(B2(j,3)-real(B2(i,3)*(real(Y(i,j)*imag(B2(j,3)+imag(Y(i,j)*real(B2(j,3); end endendfor i=1:n %对PV节点的解决，注意这时不可再将h初始化为0 if i=isb&B2(i,6)=3 h=h+1; for j=1:n OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3)*(real(Y(i,j)*real(B2(j,3)-imag(Y(i,j)*imag(B2(j,3)+imag(B2(i,3)*(real(Y(i,j)*imag(B2(j,3)+imag(Y(i,j)*real(B2(j,3); OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3)*(real(Y(i,j)*real(B2(j,3)-imag(Y(i,j)*imag(B2(j,3)-real(B2(i,3)*(real(Y(i,j)*imag(B2(j,3)+imag(Y(i,j)*real(B2(j,3); end endendOrgS%创立PVU 用于存储PV节点的初始电压PVU=zeros(n-h-1,1);t=0;for i=1:n if B2(i,6)=3 t=t+1; PVU(t,1)=B2(i,3); endendPVU%创立DetaS，用于存储有功功率、无功功率和电压幅值的不平衡量h=0;for i=1:n %对PQ节点的解决 if i=isb&B2(i,6)=2 h=h+1; DetaS(2*h-1,1)=real(B2(i,2)-OrgS(2*h-1,1); DetaS(2*h,1)=imag(B2(i,2)-OrgS(2*h,1); endendt=0;for i=1:n %对PV节点的解决，注意这时不可再将h初始化为0 if i=isb&B2(i,6)=3 h=h+1; t=t+1; DetaS(2*h-1,1)=real(B2(i,2)-OrgS(2*h-1,1); DetaS(2*h,1)=real(PVU(t,1)2+imag(PVU(t,1)2-real(B2(i,3)2-imag(B2(i,3)2; endendDetaS%创立I，用于存储节点电流参数i=zeros(n-1,1);h=0;for i=1:n if i=isb h=h+1; I(h,1)=(OrgS(2*h-1,1)-OrgS(2*h,1)*sqrt(-1)/conj(B2(i,3); endendI%创立Jacbi(雅可比矩阵)Jacbi=zeros(2*n-2);h=0;k=0;for i=1:n %对PQ节点的解决 if B2(i,6)=2 h=h+1; for j=1:n if j=isb k=k+1; if i=j %对角元素的解决 Jacbi(2*h-1,2*k-1)=-imag(Y(i,j)*real(B2(i,3)+real(Y(i,j)*imag(B2(i,3)+imag(I(h,1); Jacbi(2*h-1,2*k)=real(Y(i,j)*real(B2(i,3)+imag(Y(i,j)*imag(B2(i,3)+real(I(h,1); Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k)+2*real(I(h,1); Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1)-2*imag(I(h,1); else %非对角元素的解决 Jacbi(2*h-1,2*k-1)=-imag(Y(i,j)*real(B2(i,3)+real(Y(i,j)*imag(B2(i,3); Jacbi(2*h-1,2*k)=real(Y(i,j)*real(B2(i,3)+imag(Y(i,j)*imag(B2(i,3); Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k); Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1); end if k=(n-1) %将用于内循环的指针置于初始值，以保证雅可比矩阵换行 k=0; end end end endendk=0;for i=1:n %对PV节点的解决 if B2(i,6)=3 h=h+1; for j=1:n if j=isb k=k+1; if i=j %对角元素的解决 Jacbi(2*h-1,2*k-1)=-imag(Y(i,j)*real(B2(i,3)+real(Y(i,j)*imag(B2(i,3)+imag(I(h,1); Jacbi(2*h-1,2*k)=real(Y(i,j)*real(B2(i,3)+imag(Y(i,j)*imag(B2(i,3)+real(I(h,1); Jacbi(2*h,2*k-1)=2*imag(B2(i,3); Jacbi(2*h,2*k)=2*real(B2(i,3); else %非对角元素的解决 Jacbi(2*h-1,2*k-1)=-imag(Y(i,j)*real(B2(i,3)+real(Y(i,j)*imag(B2(i,3); Jacbi(2*h-1,2*k)=real(Y(i,j)*real(B2(i,3)+imag(Y(i,j)*imag(B2(i,3); Jacbi(2*h,2*k-1)=0; Jacbi(2*h,2*k)=0; end if k=(n-1) %将用于内循环的指针置于初始值，以保证雅可比矩阵换行 k=0; end end end endendJacbi%求解修正方程，获取节点电压的不平衡量DetaU=zeros(2*n-2,1);DetaU=inv(Jacbi)*DetaS;DetaU%修正节点电压j=0;for i=1:n %对PQ节点解决 if B2(i,6)=2 j=j+1; B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1); endendfor i=1:n %对PV节点的解决 if B2(i,6)=3 j=j+1; B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1); endendB2%开始循环*while abs(max(DetaU)prOrgS=zeros(2*n-2,1); %!初始功率参数在迭代过程中是不累加的，因此在这里必须将其初始化为零矩阵h=0;j=0;for i=1:n if i=isb&B2(i,6)=2 h=h+1; for j=1:n OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3)*(real(Y(i,j)*real(B2(j,3)-imag(Y(i,j)*imag(B2(j,3)+imag(B2(i,3)*(real(Y(i,j)*imag(B2(j,3)+imag(Y(i,j)*real(B2(j,3); OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3)*(real(Y(i,j)*real(B2(j,3)-imag(Y(i,j)*imag(B2(j,3)-real(B2(i,3)*(real(Y(i,j)*imag(B2(j,3)+imag(Y(i,j)*real(B2(j,3); end endendfor i=1:n if i=isb&B2(i,6)=3 h=h+1; for j=1:n OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3)*(real(Y(i,j)*real(B2(j,3)-imag(Y(i,j)*imag(B2(j,3)+imag(B2(i,3)*(real(Y(i,j)*imag(B2(j,3)+imag(Y(i,j)*real(B2(j,3); OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3)*(real(Y(i,j)*real(B2(j,3)-imag(Y(i,j)*imag(B2(j,3)-real(B2(i,3)*(real(Y(i,j)*imag(B2(j,3)+imag(Y(i,j)*real(B2(j,3); end endendOrgS%创立DetaSh=0;for i=1:n if i=isb&B2(i,6)=2 h=h+1; DetaS(2*h-1,1)=real(B2(i,2)-OrgS(2*h-1,1); DetaS(2*h,1)=imag(B2(i,2)-OrgS(2*h,1); endendt=0;for i=1:n if i=isb&B2(i,6)=3 h=h+1; t=t+1; DetaS(2*h-1,1)=real(B2(i,2)-OrgS(2*h-1,1); DetaS(2*h,1)=real(PVU(t,1)2+imag(PVU(t,1)2-real(B2(i,3)2-imag(B2(i,3)2; endendDetaS%创立Ii=zeros(n-1,1);h=0;for i=1:n if i=isb h=h+1; I(h,1)=(OrgS(2*h-1,1)-OrgS(2*h,1)*sqrt(-1)/conj(B2(i,3); endendI%创立JacbiJacbi=zeros(2*n-2);h=0;k=0;for i=1:n if B2(i,6)=2 h=h+1; for j=1:n if j=isb k=k+1; if i=j Jacbi(2*h-1,2*k-1)=-imag(Y(i,j)*real(B2(i,3)+real(Y(i,j)*imag(B2(i,3)+imag(I(h,1); Jacbi(2*h-1,2*k)=real(Y(i,j)*real(B2(i,3)+imag(Y(i,j)*imag(B2(i,3)+real(I(h,1); Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k)+2*real(I(h,1); Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1)-2*imag(I(h,1); else Jacbi(2*h-1,2*k-1)=-imag(Y(i,j)*real(B2(i,3)+real(Y(i,j)*imag(B2(i,3); Jacbi(2*h-1,2*k)=real(Y(i,j)*real(B2(i,3)+imag(Y(i,j)*imag(B2(i,3); Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k); Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1); end if k=(n-1) k=0; end end end endendk=0;for i=1:n if B2(i,6)=3 h=h+1; for j=1:n if j=isb k=k+1; if i=j Jacbi(2*h-1,2*k-1)=-imag(Y(i,j)*real(B2(i,3)+real(Y(i,j)*imag(B2(i,3)+imag(I(h,1); Jacbi(2*h-1,2*k)=real(Y(i,j)*real(B2(i,3)+imag(Y(i,j)*imag(B2(i,3)+real(I(h,1); Jacbi(2*h,2*k-1)=2*imag(B2(i,3); Jacbi(2*h,2*k)=2*real(B2(i,3); else Jacbi(2*h-1,2*k-1)=-imag(Y(i,j)*real(B2(i,3)+real(Y(i,j)*imag(B2(i,3); Jacbi(2*h-1,2*k)=real(Y(i,j)*real(B2(i,3)+imag(Y(i,j)*imag(B2(i,3); Jacbi(2*h,2*k-1)=0; Jacbi(2*h,2*k)=0; end if k=(n-1) k=0; end end end endendJacbiDetaU=zeros(2*n-2,1);DetaU=inv(Jacbi)*DetaS;DetaU%修正节点电压j=0;for i=1:n if B2(i,6)=2 j=j+1; B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1); endendfor i=1:n if B2(i,6)=3 j=j+1; B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1); endendB2Times=Times+1; %迭代次数加1endTimes5 总结通过本次课程设计让我有复习了一次潮流计算的有关知识，跟家清晰了什么事潮流计算以及潮流计算的在电力系统的重要性。电力系统的稳定运营状况即是正常运营状况，是指电力系统在稳定运营条件下电压、功率的分布，也称为潮流分布。电力系统分析的潮流计算是电力系统分析的一种重要的部分。通过对电力系统潮流分布的分析和计算，可进一步对系统运营的安全性，经济性进行分析、评估，提出改善措施。同步潮流分布也是电力系统规划设计的一项基本工作。整个计算过程的模型建立并不是十分复杂，但计算过程十分繁琐、计算量相称的大，并且由于枝节太多很容易算错。但是在计算潮流计算的过程中却对以往学过的电力系统分析的有关知识进行了一次较为进一步的复习。并且整个计算对计算量的规定很大，锻炼了我们的计算能力。并且对细节的把握也得到了锻炼，做题的精细限度得到了提高。参照文献1何仰赞, 温增银电力系统分析（第三版）M. 华中科技大学，2.3王守相，刘玉田 电力系统潮流计算研究现状-山东电力技术1996年05期4何仰赞，温增银. 电力系统分析（上册）第三版M. 湖北：华中科技大学出版社，5 刘同娟. MATLAB在电路分析中的应用.电气电子教学学报.6 西安交通大学等.电力系记录算M.北京：水利电力出版社，1993.127 李光琦.电力系统暂态分析M.北京： 水利电力出版社，.58何仰赞，温增银. 电力系统分析（下册）第三版M. 湖北：华中科技大学出版社，9韦化，李滨，杭乃善，等大规模水一火电力系统最优潮流的现代内点算法实现J中国电机工程学报，23(6)：13一l810Chen Luonan，Suzuki Hideki，Katou KazuoMean fieldtheory for optimal power flowJIEEE Transactions OilPower Systems，1997，12(4)：14811486