大型风电场并网运行系统暂态稳定性分析

第26卷第2期2008年4月

水　电　能　源　科　学Water Resources and Power Vol.26No.2Apr.2008

文章编号:100027709(2008)022*******

大型风电场并网运行系统暂态稳定性分析

杨　娜　张　焰

(上海交通大学电子信息与电气工程学院,上海200240)

摘要:采用时域仿真方法研究了大型风电场接入电网后的电力系统暂态稳定问题,建立了风力机传动部分两质量块的风力发电机数学模型,用PSA T 软件对大容量风电场接入的电力系统进行电网侧故障的暂态仿真,以IEEE 39节点算例进行分析。仿真结果表明,风电场并网后系统的暂态稳定性降低,通过加入TCSC 改善了系统的暂态稳定性。

关键词:电力系统;风电场;暂态稳定;PSA T 中图分类号:TM614;TM712

文献标志码:A

收稿日期:2007212212,修回日期:2008201228

作者简介:杨娜(19832),女,硕士研究生,研究方向为电力系统稳定分析和风力发电,E 2mail :nyang_sjtu @ 通讯作者:张焰(19582),女,教授、博导,研究方向为电力系统分析、规划、电力系统可靠性、配电系统自动化,E 2

mail :zhang_yan @

随着风力发电技术的迅速发展,大容量风电场并网运行势在必行。由于大型风力发电机并网

时具有和其他常规能源电厂不同的特点,可能会对电力系统稳定性产生一定影响[1～8]。建立风力发电机组及风电场的数学模型并进行仿真分析是重要的手段之一。近年来,大多侧重于风力发电系统本身的稳定性及其模型对稳定性影响的研究,而对风电并入后电网侧发生大扰动的系统暂态稳定性未作深入研究。文献[9,10]对包含风电场的电力系统小干扰稳定性进行了分析,但建模中忽略了风力机的传动部分。传动轴连接风轮和发电机两个惯量相差很大的器件必然是柔性的,在风力机模型中考虑传动轴的柔性可更准确地反映扰动后的动态过程。

本文以恒速恒频风力发电系统为研究对象,在风力发电机建模中考虑传动轴的柔性,采用风力机传动部分的两质量块模型,在PSA T 软件中搭建仿真模型,分析大型风电场并入电网侧发生大扰动时的系统暂态稳定性,并通过加入TCSC 改善了包含风电场的系统暂态稳定性。

1　风力发电机的数学模型

风电场由多台风力发电机按一定规则排列构

成。风力发电机主要由风力机和发电机等主要元

件组成,自然界的风推动风轮叶片,将风能转化为机械能。风力机的机械传动再将机械能传递到发电机的转子上。恒速恒频风力发电机由风轮、低速轴、增速齿轮箱、高速轴、发电机和无功补偿电容器组组成。1.1　风力机的数学模型

风力机机械传动部分如图1所示。

图1　实际的风力机传动模型

Fig.1　Ac t ual t rans mis s ion f or w ind t urbine

齿轮箱两侧的传动轴以两种转速旋转,连接了三个质量块,转速比为1∶n 。因可忽略齿轮箱的惯量,故本文将风轮和低速轴的惯量、低速轴的刚性系数和阻尼系数转换到高速侧(分别除以

n 2

)就可将图1的机械传动部分转化为两质量块

模型,如图2所示。 两质量块风机传动部分用s 函数表示为:

θk =s (ωt -ωg )

(1)T w -T k =(J t s +D t )ωt

(2)T k -T em =(J gen s +D g )ωg

(3)

图2　风力机传动的两质量块模型Fig.2　T wo mass model of wind t urbine t ransmission

T k=(k/s+D tg)(ωt-ωg)(4)式中,J t为风轮转换到高速侧的惯量;J gen为发电机转子惯量;θk为风轮和发电机转子间的转角差;ωt 为风轮的角速度;ωg为发电机的角速度;k为传动轴刚性系数;D tg为传动轴阻尼系数;D t为风轮的阻尼系数;D g为发电机的阻尼系数;T w为风推动风机产生的机械转矩,转换到高速侧应除以n。

其中

T w=

0V w

C p(λ)πρR2V3

2ωw T wbase V in

≤V w

0V out≤V w

(5)

式中,ρ为空气密度,kg/m3;ωw为风机的机械转速,r/s;R为风机的叶轮半径,m;C p(λ)为风力机的风能利用系数,即单位时间内风力机所吸收的风能与通过叶片旋转面的全部风能之比。按贝兹理论[11],C p,max(λ)取0.593,与风力机的叶尖速比λ(风力机叶片顶端线速度与风速之比)有关,λ=ω

w R/V;V in、V out分别为风力机的切入风速和切出风速。风速变化的时空模型由基本风V wa、阵风V wt、渐变风V wr、噪声风V wg四部分组成[10]。1.2　发电机的数学模型

异步发电机采用考虑转子机电暂态模型:

d E′d d t =

1

T′d0

[-E′d+(x-x′)i q-2πf0s T′d0E′q]

d E′q d t =

1

T′d0

[-E′q-(x-x′)i d+2πf0s T′d0E′d]

(6)

发电机机电转矩为:

T em=E′q i q+E′d i d(7)定子电压为:

U=-(r1+j x′)I+E′(8)式中,s为异步发电机滑差,s=(ωr-ω0)/ω0;ωr为发电机转子电角速度,ω0为同步转速;x′为发电机暂态电抗,x′=x1+x2x m/(x2+x m),x=x1+ x m,x为发电机电抗,x1为定子漏抗,x m为激磁电抗,x2为转子漏抗;r1、r2分别为定子电阻、转子电阻;T′d0为定子开路时转子回路的时间常数, T′d0=(x2+x m)/(2πf0r);E′=E′d+j E′q,E′为暂态电势,E′d为d轴暂态电势,E′q为q轴暂态电势。2　PSA T软件暂态仿真

PSAT为基于MATLAB环境的电力系统分析工具箱,可进行潮流、连续潮流、最优潮流的计算及分岔分析、小干扰稳定分析、时域仿真及PMU配置,并支持用户自定义模型。模型数据的输入有两种方式:①3.m文件格式输入;②PSAT 支持在SIMUL IN K环境下从模型库中拖拽元件搭接仿真算例(3.mdl文件)。SIMUL IN K环境下编译的文件最终是解释①的文本文件,由MATLAB 函数完成计算。因此,搭建暂态仿真模型时要符合元件使用的可行性规则[12],确保3.mdl文件编译成功并生成相应的3.m文件。

3　算例分析

在PSA T中搭建IEEE39节点仿真模型。IEEE39节点系统接线如图3所示,具体参数见文献[13]。发电机总有功/无功为6192.9MW/ 923.72Mvar,负荷总有功/无功为6098.1MW/ 1408.9Mvar。母线31取为松弛节点。同步发电机采用四阶模型。每台风力发电机出力2MW, J t=5.0,J gen=1

.0,k=0.3,D tg=0.02,D t= 0.02,D g=0.02,R=37.5,n=89,r1=0.01,x1= 0.1,r2=0.01,x2=0.08,x m=3.0,机端电容补偿无功使其功率因数达到0.98。

图3　IEEE39节点测试系统的PSA T仿真模型Fig.3　PSA T s imula t ion model of IEEE39bus s ys t e m 系统在1s时线路2—3靠近3处发生三相接地短路,故障在6个周波(1.1s)时由断开线路2—3而被消除。线路2—3在4s时重合闸。以下分系统无风电场、接入风电场、有风电场并接入TCSC三种情况,比较系统发生上述故障并断开线路、重合闸动作时的暂态稳定性。

・

3

9

1

・

第26卷第2期 杨　娜等:大型风电场并网运行系统暂态稳定性分析