电力系统潮流上机答辩报告

1《电力系统综合仿真》课程设计报告院系：电气与电子工程学院班级：实验电10班学号：1101180711学生姓名：李慧勇指导教师：姜彤设计周数：两周成绩：日期：2012年1月10日2《电力系统综合仿真》课程设计任务书一、目的与要求培养学生编程实现电力系统仿真计算的能力，掌握潮流计算、短路计算、暂态稳定计算的相关知识二、主要内容1．编写潮流计算程序，要求如下：2.1据给定的潮流计算任务书整理潮流计算的基础数据：节点的分类，线路模型，等值变压器模型，电压等级的归算，标幺值的计算；2.2基础数据的计算机存储：节点数据，支路数据（包括变压器）；2.3用牛顿-拉夫逊法计算；2.4根据所选潮流计算方法画流程图，划分出功能模块，有数据输入模块，导纳阵形成模块，解线性方程组模块，计算不平衡功率模块，形成雅可比矩阵模块，解修正方程模块，计算线路潮流，网损，PV节点无功功率和平衡节点功率，数据输出模块；2.5据上述模块编制程序并上机调试程序，得出潮流计算结果；2.6源程序及其程序中的符号说明集、程序流图GGG28793146532．利用潮流计算结果编写故障计算程序，要求如下：2.1发电机参数一律为PN=200MW，cosφ=0.85,19.0dx,125.0xx2.2变压器一律为Yd11接线，其中三角形侧接发电机，星形侧都直接接地2.3线路零序参数一律为103xx2.4负荷都采用Y/Y型变压器接入，中性点与系统相连侧部接地。负荷按照对称负荷考虑，纯阻抗形式。2.5计算5节点发生三相对称短路、单相接地短路、两相短路和两相接地短路时的故障电流，以及各支路电流和节点电压。3．思考题3.1潮流计算和短路电流计算的元件模型有哪些不同点？节点导纳矩阵有哪些不同？计算过程有哪些不同点？3.2短路电流计算中进行了哪些假设？3.3如果交给你一个任务，请你用已有的潮流计算软件计算北京城市电网的潮流，你应该做哪些工作？（收集哪些数据，如何整理，计算结果如何分析）3.4比较过去的计算机语言上机学习和本课程上机有哪些不同？设计中遇到的问题和解决的办法。三、进度计划序号设计内容完成时间备注1教师讲解本次设计的任务及相关知识课程设计开始前2整理潮流计算的原始数据，完成数据的读入读出模块一天3编写及调试形成导纳矩阵模块一天4计算不平衡功率，形成雅可比矩阵一天5解修正方程一天6计算线路潮流，网损，PV节点无功功率和平衡节点功率一天7读入故障计算需要的数据，计算三相短路一天8计算不对称短路一天9写设计报告一天10验收一天4四、设计成果要求1．提交任务书2．计算机计算潮流程序流程图；潮流计算开始赋电压初值形成节点导纳矩阵读入节点和支路数据得最终的节点电压解方程求电压修正量和电压新值形成Jacobi矩阵计算不平衡量∆P、∆Q判断是否符合迭代要求（用迭代次数k表示）计算平衡节点的功率、PV节点的无功、各支路功率和网损潮流结束数据输出53.计算机计算短路电流程序流程图；短路计算开始短路计算结束输出结果形成正序增广网络（各种类型的短路故障不同），并求出短路处电流）解方程得短路节点对应的各序阻抗列向量计算正、负、零序节点导纳矩阵读入节点导纳矩阵和节点电压计算各序电压并求出三相各节点电压64.完整的仿真程序(包括变量的定义，子程序的调用等，计算结果)，并有程序注释。（1）源程序#includestdafx.h#includemath.h#includestdio.h#includecomplex#includeiostream#includefstream#includeiomanip//输出格式控制头文件usingnamespacestd;usingnamespaceSystem;#includeNEquation.h//调用解方程文件#includeCNEquation.h//解虚数方程structBusData{intID;intType;doubleU,Angle,PL,QL,PG,QG,U0,G,B;}Bus[9];structBranchData{intBusA,BusB,Type;doubleR,X,B,K;}Branch[9];#definePI3.1415926intmain(arraySystem::String^^args){cout.setf(ios::left);ifstreamfile1(G:\\9data.txt);//导入节点和支路数据doublee[9]={0},f[9]={0};for(inti=0;i9;i++){file1Bus[i].IDBus[i].TypeBus[i].UBus[i].AngleBus[i].PLBus[i].QLBus[i].PGBus[i].QGBus[i].U0Bus[i].GBus[i].B;Bus[i].PL=Bus[i].PL/100;//功率化为标幺值,SB=100MVABus[i].QL=Bus[i].QL/100;Bus[i].PG=Bus[i].PG/100;Bus[i].QG=Bus[i].QG/100;e[i]=Bus[i].U0;f[i]=0;}for(inti=0;i9;i++){file1Branch[i].BusABranch[i].BusBBranch[i].TypeBranch[i].RBranch[i].XBranch[i].BBranch[i].K;}file1.close();doubleYG[9][9]={0},YB[9][9]={0};for(inti=0;i9;i++){doubleR=Branch[i].R;7doubleX=Branch[i].X;doubleB=Branch[i].B/2;doubleYLG=R/(R*R+X*X);doubleYLB=-X/(R*R+X*X);intm=Branch[i].BusA-1;intn=Branch[i].BusB-1;if(Branch[i].Type==1)//普通节点{YG[m][m]+=YLG;YG[n][n]+=YLG;YG[m][n]-=YLG;YG[n][m]-=YLG;YB[m][m]+=YLB+B+Bus[i].B;YB[n][n]+=YLB+B+Bus[i].B;YB[m][n]-=YLB;YB[n][m]-=YLB;}else//变压器节点{YG[m][m]+=YLG/(Branch[i].K*Branch[i].K);YG[m][n]-=YLG/Branch[i].K;YG[n][n]+=YLG;YG[n][m]-=YLG/Branch[i].K;YB[m][m]+=YLB/(Branch[i].K*Branch[i].K);YB[m][n]-=YLB/Branch[i].K;YB[n][n]+=YLB;YB[n][m]-=YLB/Branch[i].K;}}ofstreamfile2(G:\\节点导纳矩阵.txt);//形成节点导纳矩阵file2----------节点导纳矩阵:endl;for(inti=0;i9;i++)for(intj=0;j9;j++)file2Yi+1j+1=(YG[i][j],setw(8)YB[i][j])endl;file2.close();for(intk=1;k=10;k++)//k为求电压的迭代次数{doublea[8]={0},b[8]={0};inti;for(inti=1;i9;i++)for(intj=0;j9;j++){a[i-1]+=YG[i][j]*e[j]-YB[i][j]*f[j];b[i-1]+=YG[i][j]*f[j]+YB[i][j]*e[j];}doubleJacobi[16][16]={0};doubleH[8][8]={0},N[8][8]={0},J[8][8]={0},L[8][8]={0},R[8][8]={0},S[8][8]={0};for(inti=1;i9;i++)for(intj=1;j9;j++){H[i-1][j-1]=-YB[i][j]*e[i]+YG[i][j]*f[i];N[i-1][j-1]=YG[i][j]*e[i]+YB[i][j]*f[i];8J[i-1][j-1]=-N[i-1][j-1];L[i-1][j-1]=H[i-1][j-1];}for(inti=0;i8;i++){switch(Bus[i+1].Type){case(1)://对PQ节点求H、N、J、L阵{H[i][i]+=b[i];N[i][i]+=a[i];J[i][i]+=a[i];L[i][i]-=b[i];break;}case(2)://对PV节点求H、N、R、S阵{H[i][i]+=b[i];N[i][i]+=a[i];R[i][i]+=2*f[i+1];S[i][i]+=2*e[i+1];break;}}}for(inti=0;i8;i++){for(intj=0;j8;j++)switch(Bus[i+1].Type){case(1)://按节点类型填入Jacobi矩阵{Jacobi[2*i][2*j]=H[i][j];Jacobi[2*i][2*j+1]=N[i][j];Jacobi[2*i+1][2*j]=J[i][j];Jacobi[2*i+1][2*j+1]=L[i][j];break;}case(2):{Jacobi[2*i][2*j]=H[i][j];Jacobi[2*i][2*j+1]=N[i][j];Jacobi[2*i+1][2*j]=R[i][j];Jacobi[2*i+1][2*j+1]=S[i][j];break;}}}ofstreamfile3(G:\\雅克比矩阵.txt,ios::app);//形成Jacobi矩阵file3----------第k次迭代后雅克比矩阵:endl;for(inti=0;i16;i++){9for(intj=0;j16;j++)file3setw(10)Jacobi[i][j];file3endl;}file3.close();doubleDelta[16];//定义不平衡量for(inti=1;i9;i++){switch(Bus[i].Type){case(1)://PQ节点的不平衡量{Delta[2*(i-1)]=Bus[i].PG-Bus[i].PL-a[i-1]*e[i]-b[i-1]*f[i];Delta[2*i-1]=Bus[i].QG-Bus[i].QL-a[i-1]*f[i]+b[i-1]*e[i];break;}case(2)://PV节点的不平衡量{Delta[2*(i-1)]=Bus[i].PG-Bus[i].PL-a[i-1]*e[i]-b[i-1]*f[i];Delta[2*i-1]=Bus[i].U*Bus[i].U-f[i]*f[i]-e[i]*e[i];break;}}}ofstreamfile4(G:\\不平衡量.txt,ios::app);file4----------第k次迭代后不平衡量:endl;for(inti=0;i16;i++)file4setw(15)Delta[i];file4endl;file4.close();NEquationob;//解修正方程ob.SetSize(16);doublex[16]={0};for(inti=0;i16;i++){for(intj=0;j16;j++)ob.Data(i,j)=Jacobi[i][j];ob.Value(i)=Delta[i];}ob.Run();for(inti=0;i16;i++)x[i]=ob.Value(i);//修正量for(inti=0;i8;i++){f[i+1]+=x[2*i];//电压新值