电力系统的稳定性研究分析毕业论文.docx
- 文档编号:5972133
- 上传时间:2023-01-02
- 格式:DOCX
- 页数:42
- 大小:469.83KB
电力系统的稳定性研究分析毕业论文.docx
《电力系统的稳定性研究分析毕业论文.docx》由会员分享,可在线阅读,更多相关《电力系统的稳定性研究分析毕业论文.docx(42页珍藏版)》请在冰豆网上搜索。
电力系统的稳定性研究分析毕业论文
电力系统的稳定性研究分析毕业论文
第一章概述
第1.1节稳定性概述
电力系统是由发电机、变压器、输电线路、用电设备组成的网络,它包括通过电的或机械的方式连接在网络中的所有设备。
电力系统的运行状态由运行参量来描述。
电力系统中同步发电机只有在同步运行状态下,其送出的电磁功率为定值,同时在电力系统中各节点的电压及各支路功率潮流也都是定值,这就是电力系统的稳定运行状态。
反之,如果电力系统中各发电机不能保持同步,则发电机送出的电磁功率和全系统各节点的电压及支路的功率将发生很大幅度的波动。
如果不能使电力系统中各发电机间恢复同步运行,电力系统将持续处于失步运行状态,即电力系统失去稳定状态。
保证电力系统稳定是电力系统正常运行的必要条件。
只有在保持电力系统稳定的条件下,电力系统才能不间断的向各类用户提供合乎质量要求的电能。
电力系统失去稳定的原因是在运行中不断受到内部和外部的干扰,小的负荷波动,大的如电力元件发生短路故障等,使电气连接在一起的各同步发电机的机械输入转矩与电磁转矩失去平衡。
电力系统稳定一般按电力系统承受干扰的大小分为静态和暂态稳定两大类。
在大的干扰下电力系统的运行参数将发生很大的偏移和振荡,所以必须考虑电力系统的非线性,从电力系统的机电暂态过程来判断系统的稳定性。
第1.2节电力系统暂态稳定
电力系统在某一运行方式下,受到外界大干扰后,经过一个机电暂态过程,能够恢复到原始稳定运行方式,则认为电力系统在这一运行方式下是暂态稳定的。
电力系统暂态稳定性与干扰的形式有关,一般有三种形式:
1)突然变化电力系统的结构特性,最常见的是短路,无故障断开线路也属于这一类干扰。
2)突然增加或减少发电机出力,如切除一台容量较大的发电机。
3)突然增加或减少大量负荷,如切除或投入一个大负荷。
在电力系统受到大的干扰后,其机电暂态过程是一组非线性状态方程式,不能进行线性化,所以一般采用数值积分的时域分析法,将计算结果绘出运行参数对时间的曲线,用以判断电力系统的暂态稳定性。
第1.3节电力系统稳定性的解决
根据不同的电力系统稳定问题及特点,可以采用不同的研究方法,。
目前主要的方法是:
1)干扰下的电力系统稳定问题。
可将电力系统的数学模型进行线性化处理,所以一般用频率法,即计算电力系统参数矩阵的特征值和特征相量,可以用来确定静态和动态稳定性,设计和整定各种提高电力系统稳定性的措施和自动调节装置等。
2)对大干扰下的稳定问题,由于要求解非线性方程组,目前几乎无例外地采用时域法即用各种数值积分的方法。
但是,随着电力系统的发展,运行方式的复杂多变新的机组设备及自动调节和控制装置的投入,各种控制系统和机械—-电气系统间的相互作用,以及对电力系统供电安全性要求的日益提高,对电力系统稳定性的研究提出了一系列新的问题。
当今电力系统稳定性的问题已是一个多元的问题,它涉及不同的网络结构、运行方式、控制方式及参数、故障条件、过程的时间跨度等,如加上工程和经济上的考虑,往往是一个十分复杂和得不到一个唯一解决的问题。
因此,要求进一步发展有关研究电力系统稳定性的理论及分析方法,以适应不断变化的需要。
应用任何一种方法来研究电力系统稳定性问题,都是一种数学模拟,即用数学模型来研究物理现象的过程。
所以,在研究电力系统稳定性问题时,首先要建立电力系统的数学模型,将表示电力系统特征的主要变量用合适的数学公式联系起来,一般用一组代数方程和一组微分方程的函数关系来表示。
建立数学模型不仅是简单的选择若干数学公式,而主要是对电力系统物理过程的正确和抽象化。
在以科学研究为目的时常常希望用数学模拟的方法尽可能准确地再现某一物理过程。
这时,在构成数学模型时不仅要考虑主要因素,也要计及次要因素,最大可能使计算结果与实际过程相符合。
但是精确的数学模型往往导致复杂化。
在实际的工作中,人们要求所选择的数学模型在一定的范围内再现系统过程的实际变化规律。
所以,要求在保持研究现象的主要方面和其重要规律的合理精确前提下,忽略一些不重要因素,经过合理的假设和简化,用合适的数学模型来描述。
在电力系统稳定性研究中,根据所研究的目的和要求,以及各个不同元件在整个过程中的作用,可以选择不同详尽程度的数学模型,同时考虑到获得这些模型相应的设计或实测参数的可能性。
在本次设计中,根据设计需要我们选择了简化模型,其更加详细的在后面各节中有详细的介绍。
而本次设计,主要采用了牛顿迭代法,数值积分法对简化模型进行了编程计算,并且对结果进行分析,得到提高暂态稳定性的能力的措施和方法。
第三章数学模型的建立
第3.1节发电机的数学模型
在本次的设计中,对发电机数学模型要求了两个方程式分别为电磁功率方程和转子运动方程式;在对两公式进行推导的过程中用到了如下一关系式:
q
(3)方程式的建立(以凸极机为例):
由上图可以得到:
IqXq=Usin
Idxd二Eq-Ucos、
当发电机等效电阻被略去时
P二巳=Ulcosj-UIcos(U)二UIqcos、Uldsin:
.
将(3.3)式代入得:
同样由图可得:
Xd
上式即为发电机电磁功率方程式。
然而在本次设计中多机系统,假设发电机以一个等值电抗和该电抗后的电动势来表
示E'和Xd。
假定除了发电机节点外,已消去了网络中其他的中间节点,贝U任一发电机的电磁功率为:
:
G
Pe尸ReEili)E2G古EDEYj色亦弋0门笊
j=i
jf
式中Y为发电机电动势节点I和j之间的互导纳,Yij|为Yj的模值。
G为发电机的台数;6为Ei和Ej相量之间的夹角,即「-门。
:
i和、:
j为电动势的相对于某一相对向量的夹角;Lj为tg‘Gj/Bj。
上式表明,任一发电机发出的电磁功率是该发电机电动势相对于其他发电机电动势向量的相角差的函数。
3.1.2•转子运动方程
根据旋转物体力学可以直接写出发电机的转子关系式:
(3.7)
j九-Te
式中:
J----转子的转动惯量
Tm――原动机转矩,
Te发电机电磁转矩
:
――转子机械角加速度,
以TB去除(3.10)式两边,得出以标么值表示的运动方程式为:
假设在转子旋转上施以一个单位标么值的转矩,・订*=1,则:
0
t「M(川=M:
-M(3.12)
0
M是以秒为单位的惯性常数。
它表示在发电机转子上施以一个单位标么值的转矩,将转子有静止(w=C)拖动到同步转速(w0)所需的时间MS为
式中:
GD2——电机转子转动惯量
n电机的额定转速,
于是,发电机转子运动方程可以表示为:
(3.13)
d2、.
dt2
式中:
「•的单位是电弧度;M值是以电机额定功率KVA为基准的。
如果计算中不以电机额定功率为基准,则应根据所取得基准进行换算。
对以上所推得的转子运动方程利用分段计算法可以求得计算机程序得以实现的方法为:
(3.15)
第3.2节负荷数学模型
在本次设计中,负荷采用恒定阻抗表示,分为正序阻抗和负序阻抗
3.2.1
,正序阻抗
因为
—SP
S=UI二1-
jQL
U
U
PL-jQL
所以
U二ZL1
U2
(3--16)
Zl+
Pl-jQL
3.2.2.
负序阻抗
负序阻抗用一个经验公式来表示:
ZL2=0.05*0.35*ZL1|+jO.35ZL1(3--17)
第3.3节网络的数学模型
3.3.1.导纳矩阵的形成
abc三相电力
电力网络是一个线性网络。
线性网络可以进行线形变换。
利用对称分量法,网络可以变换为120三序网络。
设三相网络导纳型网络方程为
I(abc)网络节点三相电流列向量
Y(abc)网络节点三相导纳列向量
设网络节点数为n。
每个节点都有abc三相,因此,U(abc)和I(abc)均为3n维向量,Y呦为3n3n维矩阵。
网络节点三相电压列向量u(abc)、注入电流列向量l(abc)分别与节点三序电压列向量
u(120),注入电流列向量
I(120)有如下关系:
(abc)(120)
U1Au
(abc)(120)
I=TAI
其中Ta是三相量和三序量的变换矩阵,维数是3n3n它是对角元素为T的对角矩阵。
111
T=a2a1
2
>a门
其中a二ej12°;a二ej24°。
导纳型三序网络方程有两种排列方式。
一是按节点块排列,即
v(120)
•11
v(120)
丫12...
y1201
U
(120)q
1
「I
(120)
1
1
、/(120)
、/(120)
120
U
(120)
(120)
丫21
丫22
丫2n
2
=
I
2
■
:
1
:
1
:
式中u(120)i节点三序电压列向量
1:
120)i节点三序电流列向量
£120)i节点三序互导纳
Yi(120)i节点三序自导纳
二是按序网块排列,即
式中:
u1——正序网节点电压列向量
u2——负序网节点电压列向量
u0――零序网节点电压列向量
I1正序网节点电流列向量
12——负序网节点电流列向量
I0――零序网节点电流列向量
Y――正序网络节点导纳矩阵
Y――负序网络节点导纳矩阵
匕一一零序网络节点导纳矩阵
其余的导纳为耦合导纳。
在正常的情况下,一般假定网络是对称的,耦合导纳为零矩阵。
只有故障时,网络结构才会出现不对称。
,即各序间出现耦合导纳。
即需要修改导纳矩阵。
3.3.2.导纳矩阵的修改
当电网发生故障时,每重故障都可以找到一个修改的导纳矩阵。
修改导纳矩阵
訂⑴取决于故障点和故障类型。
它是反映故障情况的导纳矩阵。
(1)短路故障时的卅
现以i节点的b相发生接地短路故障为例,阐述如何求得短路故障的。
短路故障相当于原网络中故障点增加一条对地支路,当发生三相短路时,该支路是对称的;
当发生非三相短路时,该支路是不对称的。
为阐述方便,称该支路为修改支路,意为修改
图中z是b相接地电阻。
当发生金属性短路时,丫=1修改支路的一相导纳。
根据图,可以写出
z
网络的支路。
短路修改支路的三相电路下图所示。
z取一个最小值,它不影响计算结果。
修改支路三相电压和三相电流的关系式:
三相电压、电流和三序电压,电流有如下关系:
其中a=ej120°,a2=ej240°。
上标1,2,0分别表示正,负,零序。
将式
(2)代入
(1),得修改支路对称分量方程:
或写成:
AyU:
20=
1120
=Ii
=U:
20
[
yn
y12
y131
_y
-
1
2a
al
其中:
Ay=
y21
y22
y23
a
1
2a
3
2
i
y01
y02
y03一
-
a
a
1
即为短路修改的导纳矩阵。
短路修改支路是一条对地支路。
在按节点块排列的网络节点三序导纳矩阵中,它将使i节
点的自导纳矩阵增加Y。
在按序网块排列的网络节点三序导纳矩阵中,其相应的修改导纳矩阵为:
i
△Y=n+j*xx
2n+i|_xx乂一
厶丫⑴中有九个非零元素,用叉表示,为标出的元素均为0。
非零的九个元素为:
Y(i,i)=y11;Y(ni,i)=:
Y(i,2ni)二力°;
Y(ni,i)二y21;Y(ni,ni)二y22;:
Y(ni,2ni)=y20;
Y(2ni,i)二yo1;:
Y(2ni,i)二y°1;•:
Y(2ni,2n'i)二y。
。
;
(2)断线故障时的.-:
Y(t)
下面以a相断线为例,说明如何求得断线故障的修改导纳矩阵.-■:
Y(t),如下图:
断线处用很小的一个阻抗z表示。
z很小,对计算结果不会有影响。
1
丫二一是断线支路的一相导纳。
断线支路由正常支路和修改支路迭加而成。
正常支路已Z
在形成导纳矩阵时记及。
所以修改导纳矩阵只反映修改支路的情况。
由图(C)可以写出修
改支路三相电压和三相电流之间的关系式
_-y001
Ui;!
■in
000
*Uijb
二lb
000
1儿
利用式
(2)可以得出断线支路对称分量方程
-1-1-1uiii
22
-1-1—1屮Uij=lij
.....0.0--1-1-1」2—
或写为弓u;20胡20
y是断线修改支路的导纳矩阵。
断线修改支路是一条非对地支路。
增加一条断线修改支路,不仅影响断线支路两端节点的自导纳矩阵,还影响节点之间的互导纳矩阵,也就是,两节点的自导纳矩阵各自增加-弓
按序网块排列,相应的修改导纳矩阵为:
[
i
XXXXXX
j
XXXXXX
,八n十i
XXXXXX
△丫⑴=ni
n+j
XXXXXX
2n+i
XXXXXX
2n+j
XXXXXX
1
.一
厶丫⑴矩阵中有三十六个非零元素,用叉表示,未标出元素均为零
:
Y(i,i)二yn;:
Y(i,ni)二:
Y(i,2ni)=如°;
:
Y(i,i)二丫门;:
Y(i,ni)二『22;:
Y(i,2ni)二y?
。
;
.:
Y(2ni,i)二y01;.:
Y(2ni,ni)二y02;
:
Y(j,j)二yii;:
Y(j,nj)®;
Y(nj,j)=y2i;:
Y(nj,nj)722;
Y(2nj,j)=y°i;:
Y(2nj,nj)=y°2;
:
Y(i,j)--yii;:
Y(i,nj)-5;
Y(ni,j)--y2i;:
Y(ni,nj)-_y22;
Y(2ni,j)--y°i,Y(2ni,nj)--y°2;
:
Y(j,i)=「yii;:
Y(j,ni)一%2;
Y(nj,i)=_y2i;.:
Y(nj,ni)--y22;
:
Y(2nj,i)二-yoi;.:
Y(2nj,ni)二-y°2;
Y(2ni,2ni)=y。
。
;
Y(j,2nj)=%。
;
Y(nj,2ni)=y2。
;
Y(2nj,2nj)~。
。
;
.:
丫(i,2nj)--yio;
Y(ni,2nj)--y?
。
;
.:
丫(2ni,2nj)--y。
。
;
:
Y(j,2ni)一yi。
;
:
Y(nj,2ni)=V。
;
Y(2nj,2ni)=-y。
。
;
3.3.3网络模型的建立
暂态稳定计算中,一般不计网络的暂态过程,即网络的电气特性用代数方程描述。
网络方程为:
Y(t)*u⑴=|⑴(3.18)
上标t表示随时间而变化。
n(t)
且在(3.18)式中:
Y(t).lY^(3.19)
K=1
其中n(t)表示t时刻故障重数。
2n
将方程(3.18)展开;l&(m)y(ms)*U&(s)(3.20)
S=2j=1
当m=01、2时,I严为if,h⑴,lf),它们是i节点的零序,正序,负序的注入电流。
当s=0,1,2时,U(s)为U(0),Uj
(1),U
(2),表示j节点的零序,正序,负序的电压。
当m=0
1,2,s=0,1,2时,Yj(ms)分别为Yj(11),Yjg,Yj(10),Yj(21),Yj(22),Yj(20),Yj(01),Yj(02),Y(00),分别是
Y,丫11,丫10,Ym1,丫11,丫110,Y)1,Y)11,Y。
导纳矩阵的第i行第j列元素
式(3.20)用直角坐标表示可写为:
2n
ixm)+jiym)=迟迟(G(ms)+jB罟)*(uXjs)+juys))(3.21)
s=Qj4
由于在暂态稳定计算中,发电机定子回路采用矩阵形式为了与注入电流相匹配,因此将
由(3.24)代入(3.22)可得到对应于正序网中发电机节点g的方程:
以上(3.22)和(3.25)式是本次设计中所用的网络方程,对应于正序网的发电机节点用
方程(3.25)。
对应于正序网非发电机节点的,负序网的所有节点的,采用方程(3.22)
第四章电力系统暂态稳定分析计算的方法和内容
第4.1节数值法和直接法的简介
电力系统受到大的干扰后能否保持稳定的运行状态,还需要进一步的研究才能得到回答,目前暂态稳定性分析的基本方法可分为两类:
一类是数值计算法,在列出描述系统暂态过程的微分方程和代数方程组后,应用各种数值积分法进行计算求解,然后根据发电机转子之间的相对角的变化情况来判别稳定性。
主要步骤是:
(1)、建立一表示电力系统及其各种元件动态行为的数学模型;
(2)、设置初始运行方式和要研究的干扰(类型、地点等);
(3)、通过建立电力系统各种状态下的时域解,确定电力系统的暂态响应;
(4)、分析所得的时域解,判断是否在干扰后可达一新的和可接受的稳定运行方式,或者判断是否失去稳定。
(5)、重复对别的初始运行条件和干扰下的电力系统动态行为进行计算分析。
具体的流程图见大图纸。
另一类是直接法,其中有些方法是关于李雅普诺夫直接法进行近似处理后发展而成的实用的方法,有的则是将简单系统中的稳定判别方法推广应用于多机系统。
根据暂态稳定性的定义,在遭受扰动后如果系统是稳定的,则它最终过渡到一个稳定运行状态,那时各发电机之间的转子角度、转速和其他所有状态变量将重新保持不变,即到达一个平衡状态。
这个平衡状态是静态稳定的,否则这种稳定运行情况是不可能存在的。
对于故障后稳态运行情况下各个变量的取值,可以用状态空间Xs来表示,并称为平衡点(SFP,Stable
EquilibriumPoint)。
直接法的基本方法是,在状态空间中找出一个包围稳定平衡点Xs的区域R/,凡是属于这一区域的任何扰动,系统以后的运行都趋于稳定运行平衡点。
然而,直接法所选取的数学模型比较粗略,其计算结果的精确性尚不另人满意。
但是数值计算法是目前广泛应用的分析方法,已发展的比较成熟,并基本上能满足电力系统的规划、设计和运行中进行的离线暂态稳定分析计算的速度和精确度的要求。
第4.1节插值法简介
插值的目的就是根据给定的数据表,寻找一个解析函数:
x,近似的代替fx。
用代数多相式作为研究插值的工具,就是所谓的代数插值,对代数插值而言问题的提法是这样的,当给出了nJ个点上的函数表如下图:
X
XoX1X2……Xn
y
yoy1y2……..齐
要构建一个多项式:
x,应满足以下两个条件:
1)「x是一个不超过n次的多项式;
2)在给定的点Xii=0.1.2……n上与fx取相同值,即「Xi=yi=0.1.2.•…n我们称x为fx的插值函数,点Xi插值点。
4.2.1线性插值
给出函数表
X
Xo
y
yo
y1
如何构造一个插值函数巴(x),使%(x)满足上面两个条件呢?
简单的就是过两点(Xo,y°),
X1,%作一条直线,把直线方程表示为
■1x;=axb(4.1)
为确定a、b,把两点代入方程(5.1)中,得
{
axob=y°
a^b=%
只要Xo不等于X1,即可解a、b。
由于用直线近似的代替函数fx,所以称这种插值为线形插值。
4.2.2二次插值
线形插值是用两个点xo,yo和X1,y1来构造y二f(x)的插值函数。
下面我们用三个点
Xo,yo,X1,y1,X2,y2来构造y二f(x)过三点的插值函数。
过三点可以做出一条抛物线,假设其方程为
2
;:
2x=a0a1xa2x
为确定系数将三点代入方程得:
yo
=ao
◎Xoa2Xo
*:
=a°
a1x1a2x「
y2
QX2a2X2
当Xo,Xi,X2互异时,方程组(4.2)的存在唯一解,且可以求出系数。
从而可以进行插值。
4.2.3.分段线性插值
1概念:
设在区间a,b1上,给定n+1个插值节点
a=x0x1x2...xn=b
和相应的函数值y0,%,y2,...,yn,求作一个插值函数(x)具有下列性质:
(1)「(Xj)=yj,j=0,1,2,...,n。
(2)(x)在每个小区间上是线性函数。
插值函数(x)叫做区间a,b1上对数据(Xi,yj(i=0,1,2,,,n)的分段线性插值函数。
2用分段线性插值法计算转子运动方程
这种方法计算步骤简单,尤其适用于简单系统的手算。
0*型=360"
在计算中,常用度数来表示「•,另外将换成转差(与同步角速度之差),则转子运动方程可改写为:
(4.3)
dt
d」__2dt_Tj
式中■为标幺值。
分段计算法的基本出发点是将转子运动方程分为一系列很小的时间段,并且假定:
(1)从一个时间段的中点至下一个时间段的中点的一段时间内,过剩功率厶P保持不变,
并等于下一个时间段开始时的过剩功率。
(2)每个时间段内的相对角速度沁不变,就等于这个时间段中点的相对角速度。
这种假设与实际情况是不一致的,把连续变化的量用阶梯变化的量代替了。
但如果时段
取得足够小,误差是不大的。
一般4可取0.05—0.10s。
以下介绍应用分段计算法计算、:
_t曲线的具体步骤。
若已知n-1时段结束时的角度「g),贝U这时的电磁功率P(n书和过剩功率沖
(2)均
可求得。
由式(4.3)的第二式可得相对角速度的变化量为:
-P(nd)
=「「(nJ)P(nj)
式中△6(n)、占为度数;心P为标幺值;Tj、加的单位均为秒(S);K为常数。
其
值为360f。
■:
t2
(4.8)
'(n)=“(n」)*"'(n)
而求得;.(n)后又可求・申、仁。
」)以及■■:
(n1)。
继续这样一点一点地计算下去。
最后可
做出6_t曲线。
递推计算公式即为式(4.7)和(4.8)。
在发生故障或切除故障的瞬间,由于运行点跃变。
计算这个瞬间相对角速度的变化量时,应当用跃变前后的两个过剩功率的平均值。
A%=塑严
式中■■■P(n)为突变前的过剩功率,■-■P('n)为突变后的过剩功率
第五章上机程序简介
本次设计为分析110千伏电力系统暂态稳定分析,根据前几章讨论的数学模型和计算方法,可以编制简化模型暂态稳定计算的原理框图。
框图见附录六:
框
(1)~(8)为暂态稳定计算做好准备,这些内容将在下面进一步说明。
框(7)计算发电机暂态电动势E'和转子位置角「•,它们是暂态稳定的初值。
在简化模型中E'保持不变。
框(9)中t是暂态时间,k是时段数,
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 电力系统 稳定性 研究 分析 毕业论文
![提示](https://static.bdocx.com/images/bang_tan.gif)