电力系统分析课程设计-基于Matlab的电力系统潮流计算文章资料讲解.doc
- 文档编号:607548
- 上传时间:2023-04-29
- 格式:DOC
- 页数:26
- 大小:1.27MB
电力系统分析课程设计-基于Matlab的电力系统潮流计算文章资料讲解.doc
《电力系统分析课程设计-基于Matlab的电力系统潮流计算文章资料讲解.doc》由会员分享,可在线阅读,更多相关《电力系统分析课程设计-基于Matlab的电力系统潮流计算文章资料讲解.doc(26页珍藏版)》请在冰点文库上搜索。
《电力系统分析》课程设计
基于Matlab的电力系统潮流计算
专业电气工程及其自动化
班级1402班
姓名郭娟
学号14039208
同组成员孙帆郭娟陈朝辉
苏新波吕应发
指导教师楚冰清
完成时间2016.12.30
1
目录
1概述 1
1.1设计目的 1
1.2设计题目 1
1.3设计内容 1
2电力系统潮流计算概述 2
2.1电力系统简介 2
2.2潮流计算简介 2
2.3潮流计算意义及其发展 2
3Matlab概述 3
3.1Matlab简介 3
3.2Matlab应用 4
4潮流计算过程 4
4.1系统图及参数 5
4.1.1系统图 5
4.1.2各元件参数 5
4.2电网潮流计算思路 5
4.3潮流计算过程 5
4.3.1各元件参数计算 5
4.3.2绘制等效电路图 7
4.3.3功率分布计算 7
5Matlab程序及运行结果 10
5.1Matlab程序 10
5.2运行结果 20
6设计心得与体会 23
参考文献 25
1概述
1.1设计目的
在如今的社会,电力已经成为人们必不可少的需求,而建立结构合理的大型电力系统不仅便于电能生产与消费的集中管理、统一调度和分配,减少总装机容量节省动力设施投资,且有利于地区能源资源的合理开发利用,更大限度地满足地区国民经济日益增长的用电需要。
电力系统建设往往是国家及地区国民经济发展规划的重要组成部分。
电力系统的出现,使高效、无污染、使用方便、易于调控的电能得到广泛应用,推动了社会生产各个领域的变化,开创了电力时代,发生了第二次技术革命。
电力系统的规模和技术水准已成为一个国家经济发展水平的标志之一。
电力系统稳态分析包括潮流计算(或潮流分析)和静态安全分析。
潮流计算针对电力革统各正常运行方式,而静态安全分析则要研究各种运行方式下个别系统元件退出运行后系统的状况。
其目的是校验系统是否能安全运行,即是否有过负荷的元件或电压过低的母线等。
原则上讲,静态安全分析也可U用潮流计算来代替。
但是一般静态安全分析需要校验的状态数非常多,用严格的潮流计算来分析这些状态往往计算量过大,因此不得不寻求一些特殊的算法以满足要求。
牛顿法是数学中解决非线性方程式的典型方法,有较好的收敛性。
解决电力系统潮流计算问题是以导纳距阵为基础的,因此,只要在迭代过程中尽可能保持方程式系数距阵的稀疏性,就可以大大提高牛顿法潮流程序的放率。
自从20世纪60年代中期利用了最佳顺序消去法以后,牛顿法在收敛性、内存要求、速度方面都超过了阻抗法,成为直到目前仍在广泛采用的优秀方法。
1.2设计题目
电力系统潮流计算(牛顿-拉夫逊法)
1.3设计内容
(1)根据电力系统网络推导电力网络数学模型,写出节点导纳矩阵;
(2)赋予各节点电压变量(直角坐标系形式)初值后,求解不平衡量;
(3)形成雅可比矩阵;
(4)求解的电压变量达到所要求的精度时,再计算各支路功率分布、功率损耗和平衡节点功率;
(5)上机编程调试;
(6)计算分析给定系统潮流分析并与手工计算结果做比较分析;
(7)书写课程设计报告。
2电力系统潮流计算概述
2.1电力系统简介
电力系统潮流计算是研究电力系统稳态运行情况的一种基本电气计算。
它的任务是根据给定的运行条件和网路结构确定整个系统的运行状态,如各母线上的电压(幅值及相角)、网络中的功率分布以及功率损耗等。
电力系统潮流计算是电力系统稳定计算和故障分析的基础。
潮流计算可以用传统的手工方式进行,也可以计算机软件完成。
两种方法各有优缺点。
手工方式可用来计算一些接线较简单的电力网,但若将其用于接线复杂的电力网则计算量过大,难于保证计算准确性;计算机方式从数学上看可归结为用数值方法解非线性代数方程,数学逻辑简单完整,可快速精确地完成计算,但其缺点是物理概念不明显,物理规律被埋没在循环往复的数值求解过程中,基本原理不太明显。
2.2潮流计算简介
潮流计算,电力学名词,指在给定电力系统网络拓扑、元件参数和发电、负荷参量条件下,计算有功功率、无功功率及电压在电力网中的分布。
潮流计算是根据给定的电网结构、参数和发电机、负荷等元件的运行条件,确定电力系统各部分稳态运行状态参数的计算。
通常给定的运行条件有系统中各电源和负荷点的功率、枢纽点电压、平衡点的电压和相位角。
待求的运行状态参量包括电网各母线节点的电压幅值和相角,以及各支路的功率分布、网络的功率损耗等。
2.3潮流计算的意义及其发展
电力系统潮流计算是电力系统分析中的一种最基本的计算,是对复杂电力系统正常和故障条件下稳态运行状态的计算。
潮流计算的目标是求取电力系统在给定运行状态的计算,即节点电压和功率分布,用以检查系统各元件是否过负荷。
各点电压是否满足要求,功率的分布和分配是否合理以及功率损耗等。
对现有的电力系统的运行和扩建,对新的电力系统进行规划设计以及对电力系统进行静态和稳态分析都是以潮流计算为基础。
潮流计算结果可用如电力系统稳态研究,安全估计或最优潮流等对潮流计算的模型和方法有直接影响。
实际电力系统的潮流技术那主要采用牛顿—拉夫逊法。
运行方式管理中,潮流是确定电网运行方式的基本出发点;在规划领域,需要进行潮流分析验证规划方案的合理性;在实时运行环境,调度员潮流提供了多个在预想操作情况下电网的潮流分布以及校验运行可靠性。
在电力系统调度运行的多个领域问题是研究电力系统稳态问题的基础和前提。
在用数字见算计算机解电力系统潮流问题的开始阶段,普遍采取以节点导纳矩阵为基础的逐次代入法。
这个方法的原理比较简单,要求的数字计算机内存量比较差下,适应50年代电子计算机制造水平和当时电力系统理论水平,但它的收敛性较差,当系统规模变大时,迭代次数急剧上升,在计算中往往出现迭代不收敛的情况。
这就迫使电力系统的计算人员转向以阻抗矩阵为基础的逐次代入法。
阻抗法改善了系统潮流计算问题的收敛性,解决了导纳无法求解的一些系统的潮流计算,在60年代获得了广泛的应用,阻抗法德主要缺点是占用计算机内存大,每次迭代的计算量大。
当系统不断扩大时,这些缺点就更加突出,为了克服这些缺点,60年代中期发展了以阻抗矩阵为基础的分块阻抗法。
这个方法把一个大系统分割为几个小的地区系统,在计算机内只需要存储各个地区系统的阻抗矩阵及它们之间联络的阻抗,这样不仅大幅度的节省了内存容量,同时也提高了计算速度。
克服阻抗法缺点是另一个途径是采用牛顿-拉夫逊法。
这是数学中解决非线性方程式的典型方法,有较好的收敛性。
在解决电力系统潮流计算问题时,是以导纳矩阵为基础的,因此,只要我们能在迭代过程中尽可能保持方程式系数矩阵的稀疏性,就可以大大提高牛顿法潮流程序的效率。
自从60年代中期,牛顿法中利用了最佳顺序消去法以后,牛顿法在收敛性。
内存要求。
速度方面都超过了阻抗法,成为了60年代末期以后广泛采用的优秀方法。
3Matlab概述
3.1Matlab简介
MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。
它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中,为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案,并在很大程度上摆脱了传统非交互式程序设计语言(如C、Fortran)的编辑模式,代表了当今国际科学计算软件的先进水平。
MATLAB和Mathematica、Maple并称为三大数学软件。
它在数学类科技应用软件中在数值计算方面首屈一指。
MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等,主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。
3.2Matlab的应用
MATLAB的基本数据单位是矩阵,它的指令表达式与数学、工程中常用的形式十分相似,故用MATLAB来解算问题要比用C,FORTRAN等语言完成相同的事情简捷得多,并且mathwork也吸收了像Maple等软件的优点,使MATLAB成为一个强大的数学软件。
在新的版本中也加入了对C,FORTRAN,C++,JAVA的支持。
可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用,此外许多的MATLAB爱好者都编写了一些经典的程序,用户可以直接进行下载就可以用。
MATLAB产品族可以用来进行以下各种工作:
数值分析
数值和符号计算
工程与科学绘图
控制系统的设计与仿真
数字图像处理技术
数字信号处理技术
通讯系统设计与仿真
财务与金融工程
MATLAB的应用范围非常广,包括信号和图像处理、通讯、控制系统设计、测试和测量、财务建模和分析以及计算生物学等众多应用领域。
附加的工具箱(单独提供的专用MATLAB函数集)扩展了MATLAB环境,以解决这些应用领域内特定类型的问题。
4潮流计算过程
4.1系统图及参数
4.1.1系统图
图1潮流计算设计图
4.1.2各元件参数
变压器T1、T2:
SFL1-16000/110,(121±2×2.5﹪)/6.3,ΔPs=110kW,ΔP0=10.5kW,U0﹪=10.5,I0﹪=0.9;
变压器T3:
SFL1-8000/110,(110±5﹪)/6.6,ΔPs=52kW,ΔP0=12.76kW,
Us﹪=10.5,I0﹪=1.1;
变压器T4:
2×SFL1-16000/110,(110±2×2.5﹪)/10.5,ΔPs=62kW,ΔP0=11.6kW,Us﹪=10.5,I0﹪=1.10;
导线型号均为LGJ-150,参数r0=0.21Ω/km,x0=0.4Ω/km,b0=2.8×10-6S/km。
4.2电网潮流计算思路
(1)计算各元件参数,画出等值电路;
(2)进行网络潮流计算;
(3)不满足供电要求,进行调压计算。
4.3潮流计算过程
4.3.1各元件参数计算
①120Km线路
(1)
②100Km线路
(2)
③70Km线路
(3)
④变压器T1,T2
(4)
⑤变压器T3
(5)
⑥变压器T4
(6)
4.3.2绘制等效电路
4.3.3功率分布计算
图2环形网络的等效图
(1)各元件功率损耗
①两台T4变压器并联损耗
(7)②T3变压器损耗
(8)③100Km与70Km线路交点4末端功率损耗
(9)④120Km与100Km线路交点3末端功率损耗
(10)
⑤1.4间100Km线路损耗
(11)
⑥1.3间120Km线路损耗
(12)
⑦2.4间70Km线路损耗
(13)
⑧2.3间100Km线路损耗
(14)
⑨位置1点总损耗:
(15)⑩位置2点总损耗:
(16)
3.3.4调压计算
计算1.4线路上的电压值
(17)
位置4由G1提供的电压为
(18)
由于119.84的输入电压大于110额定值,所以调压关系不满足。
(经验证其他3路均不满足关系),所以需要降低121KV端的输出电压和提高110KV的输额定值。
的变压器取输出
的变压器取输入
调压后:
(19)
所以满足调压关系。
经验证其他三路均满足,调压成功。
5Matlab程序及运行结果
5.1Matlab程序
clear;
n=8;
nl=8;
isb=1;
pr=0.00001;
B1=[128.5+20.1i0.000556i10;
1413.6+32.16i0.0002224i10;
1613.6+32.16i0.0002224i10;
231.495+40.335i01.051;
451.78+53.885i01.0251;
4610.2+24.12i0.0001668i10;
671.495+40.335i01.0251;
686.8+16.08i0.0004448i10;
891.78+53.885i01.0251;
8108.5+20.1i0.000556i10];
B2=[0022922901;
00220002;
050+30.987i220002;
00220002;
040+27.79i220002;
00220002;
050+30.987i220002;
00220002;
060+37.18i220002;
200022922903];
Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl);
SB=100;UB=220;
ifym~=0
%SB=input('请输入功率基准值:
SB=');
%UB=input('请输入电压基准值:
UB=');
YB=SB./UB./UB;
BB1=B1;
BB2=B2;
fori=1:
nl
B1(i,3)=B1(i,3)*YB;
B1(i,4)=B1(i,4)./YB;
end
disp('B1矩阵B1=');
disp(B1)
fori=1:
n
B2(i,1)=B2(i,1)./SB;
B2(i,2)=B2(i,2)./SB;
B2(i,3)=B2(i,3)./UB;
B2(i,4)=B2(i,4)./UB;
B2(i,5)=B2(i,5)./SB;
end
disp('B2矩阵B2=');
disp(B2)
end
%%%---------------------------------------------------
fori=1:
nl %支路数
ifB1(i,6)==0 %左节点处于低压侧
p=B1(i,1);q=B1(i,2);
else
p=B1(i,2);q=B1(i,1);
end
Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5)); %非对角元
Y(q,p)=Y(p,q);
Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2; %对角元K侧
Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2; %对角元1侧
end
%求导纳矩阵
disp('导纳矩阵Y=');
disp(Y)
%----------------------------------------------------------
G=real(Y);B=imag(Y); %分解出导纳阵的实部和虚部
fori=1:
n %给定各节点初始电压的实部和虚部
e(i)=real(B2(i,3));
f(i)=imag(B2(i,3));
V(i)=B2(i,4); %PV节点电压给定模值
end
fori=1:
n %给定各节点注入功率
S(i)=B2(i,1)-B2(i,2); %i节点注入功率SG-SL
B(i,i)=B(i,i)+B2(i,5); %i节点无功补偿量
end
%===================================================================
P=real(S);Q=imag(S);
ICT1=0;IT2=1;N0=2*n;N=N0+1;a=0;
whileIT2~=0
IT2=0;a=a+1;
fori=1:
n
ifi~=isb %非平衡节点
C(i)=0;D(i)=0;
forj1=1:
n
C(i)=C(i)+G(i,j1)*e(j1)-B(i,j1)*f(j1);%Σ(Gij*ej-Bij*fj)
D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1);%Σ(Gij*fj+Bij*ej)
end
P1=C(i)*e(i)+f(i)*D(i);%节点功率P计算eiΣ(Gij*ej-Bij*fj)+fiΣ(Gij*fj+Bij*ej)
Q1=C(i)*f(i)-e(i)*D(i);%节点功率Q计算fiΣ(Gij*ej-Bij*fj)-eiΣ(Gij*fj+Bij*ej)
%求P',Q'
V2=e(i)^2+f(i)^2; %电压模平方
%=========以下针对非PV节点来求取功率差及Jacobi矩阵元素=========
ifB2(i,6)~=3 %非PV节点
DP=P(i)-P1; %节点有功功率差
DQ=Q(i)-Q1;%节点无功功率差
%===============以上为除平衡节点外其它节点的功率计算=================
%=================求取Jacobi矩阵===================
forj1=1:
n
ifj1~=isb&j1~=i %非平衡节点&非对角元
X1=-G(i,j1)*e(i)-B(i,j1)*f(i); %dP/de=-dQ/df
X2=B(i,j1)*e(i)-G(i,j1)*f(i); %dP/df=dQ/de
X3=X2;%X2=dp/dfX3=dQ/de
X4=-X1;%X1=dP/deX4=dQ/df
p=2*i-1;q=2*j1-1;
J(p,q)=X3;J(p,N)=DQ;m=p+1;
J(m,q)=X1;J(m,N)=DP;q=q+1;
J(p,q)=X4;J(m,q)=X2;
elseifj1==i&j1~=isb %非平衡节点&对角元
X1=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);%dP/de
X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);%dP/df
X3=D(i)+B(i,i)*e(i)-G(i,i)*f(i);%dQ/de
X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);%dQ/df
p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;%扩展列△Q
m=p+1;
J(m,q)=X1;q=q+1;J(p,q)=X4;J(m,N)=DP;%扩展列△P
J(m,q)=X2;
end
end
else
%===============下面是针对PV节点来求取Jacobi矩阵的元素===========
DP=P(i)-P1; %PV节点有功误差
DV=V(i)^2-V2; %PV节点电压误差
forj1=1:
n
ifj1~=isb&j1~=i %非平衡节点&非对角元
X1=-G(i,j1)*e(i)-B(i,j1)*f(i);%dP/de
X2=B(i,j1)*e(i)-G(i,j1)*f(i);%dP/df
X5=0;X6=0;
p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV;
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 电力系统 分析 课程设计 基于 Matlab 潮流 计算 文章 资料 讲解