CN111525549A - 发电机组对直驱风电机组并网次同步振荡特性的分析方法 - Google Patents

发电机组对直驱风电机组并网次同步振荡特性的分析方法 Download PDF

Info

Publication number
CN111525549A
CN111525549A CN202010272580.5A CN202010272580A CN111525549A CN 111525549 A CN111525549 A CN 111525549A CN 202010272580 A CN202010272580 A CN 202010272580A CN 111525549 A CN111525549 A CN 111525549A
Authority
CN
China
Prior art keywords
direct
generator set
grid
frequency
wind turbine
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202010272580.5A
Other languages
English (en)
Other versions
CN111525549B (zh
Inventor
肖仕武
张馨元
高陈
王利超
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
North China Electric Power University
Electric Power Research Institute of State Grid Xinjiang Electric Power Co Ltd
Original Assignee
North China Electric Power University
Electric Power Research Institute of State Grid Xinjiang Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by North China Electric Power University, Electric Power Research Institute of State Grid Xinjiang Electric Power Co Ltd filed Critical North China Electric Power University
Priority to CN202010272580.5A priority Critical patent/CN111525549B/zh
Publication of CN111525549A publication Critical patent/CN111525549A/zh
Application granted granted Critical
Publication of CN111525549B publication Critical patent/CN111525549B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/24Arrangements for preventing or reducing oscillations of power in networks
    • H02J3/241The oscillation concerning frequency
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/381Dispersed generators
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2300/00Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20The dispersed energy generation being of renewable origin
    • H02J2300/28The renewable source being wind energy
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/76Power conversion electric or electronic aspects

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Eletrric Generators (AREA)
  • Wind Motors (AREA)

Abstract

本发明公开了属于电力系统建模及分析领域的一种发电机组对直驱风电机组并网次同步振荡特性的分析方法,在建立准确的考虑多质块轴系发电机组dq轴阻抗模型的基础上,将多质块轴系发电机组机端dq轴阻抗特性以及直驱风电机组和外部交流电网dq轴阻抗特性并联,依据波特图判据,判断发电机组多质块轴系扭振模态频率与直驱风电机组并网电气谐振频率之间的关系;并基于阻抗分析法,应用传递函数的阻抗极点判据,从而准确量化直驱风机对火电机组并网SSO影响,得到发电机组对直驱风机并网次同步振荡影响特性;此时,从而在机理上对含直驱风机的风火组合不同送出系统SSO做出稳定性判断,具有较高准确性。

Description

发电机组对直驱风电机组并网次同步振荡特性的分析方法
技术领域
本发明属于电力系统建模及分析领域,特别涉及一种发电机组对直驱风电机组并网次同步振荡特性的分析方法;具体来说是一种针对考虑多质块汽轮机轴系建模,研究火电机组对直驱风电机组并网次同步振荡相互作用影响,从直驱风电场与发电机组并网发生的振荡频率远近的差异提出的分析方法。
背景技术
次同步谐振(SSR)的本质是电力系统特殊的机电耦合作用引起的汽轮发电机转轴的扭振作用。汽轮发电机转子可以被视为多个集中质量块弹性连接。在次同步振荡中,多质块汽轮发电机轴系在受扰后,各个质量块在同步旋转的同时,还会发生相互间的相对扭转振荡。由于系统对该振荡所呈现的弱阻尼、无阻尼、甚至负阻尼特性,使这种振荡的振幅呈逐渐增大的趋势。随着新能源接入电力系统,电力系统发生一些显著变化,近年来,国内外有些地区,直驱风电机组在不含串补的工况下,发生了持续性的次同步振荡问题,导致电能质量恶化,由此增加了系统损耗,危及电网的安全稳定运行,对人民生命和财产安全构成严重威胁。因此有必要就弱电网条件下在建立准确的考虑多质块的发电机数学模型基础上,对直驱风机组对火电机组并网SSO影响展开稳定性分析。
本发明提出了一种基于模态频率差异的发电机组对直驱风电场并网次同步振荡特性的分析方法,在建立考虑多质块发电机阻抗模型的基础上,将直驱风电场并网振荡频率分为相距轴系扭振频率较远和较近两类来进行影响分析,并通过阻抗分析法,应用传递函数的阻抗极点判据,从而准确量化直驱风机对火电机组并网SSO影响,得到发电机组对直驱风机并网次同步振荡影响特性。
发明内容
本发明提出一种发电机组对直驱风电机组并网次同步振荡特性的分析方法,其特征在于,该方法是一种基于模态频率差异的发电机组对直驱风电机组并网次同步振荡特性的分析方法;包括如下步骤:
S1:建立多质块轴系发电机组的机端dq轴阻抗模型;
S2:基于波特图判据,判断发电机组多质块轴系扭振频率与直驱风电机组并网电气谐振频率之间的关系;
S3:如果发电机组轴系模态频率与直驱风电机组并网电气谐振频率相差大于2HZ,则认为直驱风电机组并网振荡频率相距轴系扭振频率远,于是将直驱风电机组等效为dq轴下一固定RL阻抗;
S4:如果发电机组轴系模态频率与直驱风电机组并网电气谐振频率相差小于2HZ,则认为直驱风电机组并网振荡频率相距轴系扭振频率近,于是将风电、火电组合系统划分为两个组合子系统,并将多质块轴系发电机组机端dq轴阻抗特性以及直驱风电机组和外部交流电网dq轴阻抗特性并联,应用阻抗函数的极点判断整个系统的稳定性。
所述步骤1建立多质块轴系发电机组的机端dq轴阻抗模型;首先建立发电机轴系模型,对发电机组其他各个部分线性化数学模型,做一定简化,忽略汽轮机调速器影响、励磁调节,并认为励磁电压恒定;其多质块轴系的运动方程在运行点线性化后,得到:
Figure BDA0002443650500000031
其中,ωb为基准角速度。
Δω=[Δω1,Δω2,Δω3,Δω4,Δω5,Δω6]Τ,ωi(i=1,…,6)表示轴系各质量块的角速度;Δδ=[Δδ1,Δδ2,Δδ3,Δδ4,Δδ5,Δδ6]Τ,δi(i=1,…,6)表示轴系各质量块的角位移;Μ=diag[M1,M2,M3,M4,M5,M6],Mi(i=1,…,6)表示各质量块惯性时间常数;
ΔΤ=[ΔT1,ΔT2,ΔT3,ΔT4,ΔT5,ΔT6]Τ,Ti(i=1,…,6)表示各质量块上机械力矩,其中,质量块5对应发电机转子,ΔT5=-ΔTe,ΔTe为发电机电磁转矩,质量块6对应励磁机转子,取ΔT6=0;T代表各质量块的惯性时间常数,p为微分算子。
其中,Dii(i=1,…,6)为轴系各质量块自阻尼系数,Dij(i=1,…,6)为相邻质量块i和j之间互阻尼系数;Kij(i=1,…,6)为质量块i和j之间的弹性系数
所述步骤2基于波特图判据,基于风电机组阻抗和电网阻抗频率特性曲线的波特图判据与Nyquist判据的基本原理相似。通过解析法或者实测法获得风电机组和电网的阻抗-频率波特图,找出直驱风电机组并网发生的振荡频率;通过分析风电机组阻抗幅频曲线与电网阻抗幅频曲线交点频率对应的二者相位差
Figure BDA0002443650500000032
判断系统稳定性,若
Figure BDA0002443650500000033
系统稳定,反之,系统不稳定。
所述如果发电机组轴系模态频率与直驱风电机组并网电气谐振频率相差小于2HZ,则认为直驱风机并网振荡频率相距轴系扭振频率较近,于是将风火组合系统划分为两个组合子系统,并将多质块轴系发电机组机端dq轴阻抗特性以及直驱风电机组和外部交流电网dq轴阻抗特性并联,应用阻抗函数的极点判断整个系统的稳定性。
所述发电机的的次同步谐振实质上是与自然振荡频率的互补频率,接近于发电机轴系的模态2频率(26.56Hz),这一振荡趋势增大,激发了模态2的强烈的轴系扭振,从而形成机械系统与电气系统之间的相互激励,在系统总阻尼为负或者很小的情况下,相互激励引起的振荡难以维持导致系统失稳。
与现有技术相比,本发明的有益效果是
一种基于模态频率差异的发电机组对直驱风电机组并网次同步振荡特性的分析方法。在建立准确的考虑多质块轴系发电机组dq轴阻抗模型的基础上,基于波特图判据,将直驱风电机组并网振荡频率分为相距轴系扭振频率较远和较近两类来进行影响分析。在相距轴系扭振频率较远时,将直驱风电机组等效为dq轴下一固定RL阻抗,量化分析直驱风电机组对火电机组次同步振荡的影响;在相距轴系扭振频率较近时,通过已建立的直驱风电机组和火电机组的阻抗模型,应用传递函数的极点判断整个系统的稳定性。从而本发明在一定程度有效分析了直驱风电机组对火电机组并网次同步特性,具有较高准确性。
附图说明
图1为发电机组对直驱风电机组并网次同步振荡特性的分析流程图。
图2为风电机组并网系统阻抗-频率波特图;其中a,直驱风机振荡频率示意图;b,可能发生振荡区域示意图。
图3为具有补偿线路的简单辐射状系统;
图4为含有等效风机的补偿线路的简单辐射状系统;
图5为系统判稳传递框图;
图6为直驱风机有功输出波形;
图7为直驱风机有功频谱图;
图8为发电机的各模态分量曲线图,其中a为模态s21;b为模态s26。
具体实施方式
本发明是提出一种发电机组对直驱风电机组并网次同步振荡特性的分析方法,具体是一种基于模态频率差异的发电机组对直驱风电机组并网次同步振荡特性的分析方法。包括以下步骤:
S1:建立多质块轴系发电机组的机端dq轴阻抗模型;
S2:基于波特图判据,判断发电机组多质块轴系扭振模态频率与直驱风电机组并网电气谐振频率之间的关系;
S3:如果发电机组轴系模态频率与直驱风电机组并网电气谐振频率相差大于2HZ,则认为直驱风机并网振荡频率相距轴系扭振频率较远,于是可将直驱风电机组等效为dq轴下一固定RL阻抗;
S4:如果发电机组轴系模态频率与直驱风电机组并网电气谐振频率相差小于2HZ,则认为直驱风电机组并网振荡频率相距轴系扭振频率较近,于是将风火组合系统划分为两个组合子系统,并将多质块轴系发电机组机端dq轴阻抗特性以及直驱风电机组和外部交流电网dq轴阻抗特性并联,应用阻抗函数的极点判断整个系统的稳定性。
以下结合附图和实施例对本发明做进一步的说明。
如图1所示发电机组对直驱风电机组并网次同步振荡特性的分析流程图,包括步骤:
(1)建立多质块轴系发电机组的机端dq轴阻抗模型,
首先建立发电机轴系模型,轴系的运动方程在运行点线性化后,得到:
Figure BDA0002443650500000061
其中,ωb为基准角速度。
Δω=[Δω1,Δω2,Δω3,Δω4,Δω5,Δω6]Τ,ωi(i=1,…,6)表示轴系各质量块的角速度;Δδ=[Δδ1,Δδ2,Δδ3,Δδ4,Δδ5,Δδ6]Τ,δi(i=1,…,6)表示轴系各质量块的角位移;Μ=diag[M1,M2,M3,M4,M5,M6],Mi(i=1,…,6)表示各质量块惯性时间常数;ΔΤ=[ΔT1,ΔT2,ΔT3,ΔT4,ΔT5,ΔT6]Τ,Ti(i=1,…,6)表示各质量块上机械力矩,其中,质量块5对应发电机转子,ΔT5=-ΔTe,ΔTe为发电机电磁转矩,质量块6对应励磁机转子,一般取ΔT6=0;
Figure BDA0002443650500000071
为轴系机械阻尼矩阵,其中Dii(i=1,…,6)为轴系各质量块自阻尼系数,Dij(i=1,…,6)为相邻质量块i和j之间互阻尼系数,保守计算中,自阻尼系数和互阻尼系数一般取为0;
Figure BDA0002443650500000072
为轴系质量块弹性系数矩阵,其中Kij(i=1,…,6)为质量块i和j之间的弹性系数。
对发电机组其他各个部分线性化数学模型,做一定简化,忽略汽轮机调速器影响、励磁调节,并认为励磁电压恒定。
将Δω5=Δω代入同步发电机模型,Δuf=0,整理得
Figure BDA0002443650500000081
而轴系模型的状态空间模型如下:
Figure BDA0002443650500000082
忽略励磁机上产生电枢反应转矩ΔT6,同时忽略ΔT1~ΔT4,保留Δδ5=Δδ和ΔT5=-ΔTe两个状态变量,消去其他变量,整理得:
其中,
Figure BDA0002443650500000083
继续整理得:
Figure BDA0002443650500000084
将以上表达式由时域转换到频域,即用s代替微分算子p即可,那么有:
Figure BDA0002443650500000094
对于Δω5有:
Figure BDA0002443650500000091
不妨令Q=A第11行(sI-A)-1B,则有:
Figure BDA0002443650500000092
而电磁磁转矩和电流的关系存在
Figure BDA0002443650500000093
于是变换到频域里有:
Figure BDA0002443650500000101
将上式代入上述同步发电机模型,得到
Figure BDA0002443650500000102
由此,根据以上各式,进一步整理,可得:
Figure BDA0002443650500000111
Figure BDA0002443650500000112
式子(14)、(15)组成了考虑多质块的同步发电机状态空间方程。
Figure BDA0002443650500000113
在以上状态方程中,状态变量X3,输出量Y3,输入量u3以及各系数矩阵(A3、B3、C3、D3)分别为:
Figure BDA0002443650500000114
由于以发电机的端口电流为idq输出变量,以端口电压udq为输入变量,则发电机的导纳变为:
Figure BDA0002443650500000121
其阻抗为:
Figure BDA0002443650500000122
(2)基于波特图找到直驱风电机组并网发生的振荡区域
基于风电机组阻抗和电网阻抗频率特性曲线的波特图判据与Nyquist判据的基本原理相似。通过解析法或者实测法获得风电机组和电网的阻抗-频率波特图,通过分析风电机组阻抗幅频曲线与电网阻抗幅频曲线交点频率对应的二者相位差
Figure BDA0002443650500000123
判断系统稳定性,若
Figure BDA0002443650500000124
系统稳定,反之,系统不稳定。
波特图判据不仅能够定性地判断系统稳定性,还能通过阻抗-相频曲线找出系统可能发生的振荡区域(如图2所示其中a,直驱风机振荡频率示意图;b,可能发生振荡区域示意图),这也是波特图判据相较于Nyquist判据的优势所在。
根据该波特图判据,可大致找出直驱风电机组并网发生的振荡频率λm,对于火电机组而言,更多关心的是轴系扭振频率λp。
在直驱风机振荡频率λm相距轴系频率λp大于2Hz时,即可认为振荡频率相距轴系扭振频率较远,此时有:
pm|>2Hz (19)
而在直驱风电机组振荡频率λm相距轴系频率λp小于2Hz时,即可认为振荡频率相距轴系扭振频率较近,此时有:
pm|<2Hz (20)
故此,分为以下两部分,即直驱风电机组并网振荡频率λm相距轴系扭振频率λp大于2Hz和小于2Hz两类来进行分析。
(3)判断当振荡频率相距轴系扭振频率大于2Hz时,
在直驱风机振荡频率λm相距轴系频率λp大于2Hz时,认为直驱风机振荡频率相距轴系扭振频率较远;对于图2所示简单的辐射状电力系统,自然振荡频率为:
Figure BDA0002443650500000131
在直驱风机并网前:
发电机次暂态电抗X″可用发电机参数X″d代替,网侧集中参数设置:R=3.4562Ω、L=1.34H、C=26.36uF
变压器电抗:
Figure BDA0002443650500000132
那么该系统的自然振荡频率为:
Figure BDA0002443650500000133
f0-fer=50-23.0629=26.937Hz。
可以看出,发电机的的次同步谐振实质上是与自然振荡频率的互补频率,接近于发电机轴系的模态2频率(26.56Hz),这一振荡趋势增大,激发了模态2的强烈的轴系扭振,从而形成机械系统与电气系统之间的相互激励,在系统总阻尼为负或者很小的情况下,相互激励引起的振荡难以维持导致系统失稳。
而理论推导的直驱风机dq阻抗为:
Figure BDA0002443650500000141
不妨将工频f0,代入上式,即s=j2πf0替换上述式子中的s,可得直驱风机在工频f0处的阻抗为:
Figure BDA0002443650500000142
对于多质块发电机,此振荡频率较扭振频率距离远,因此在工频下,直驱风机的阻抗等效为RL串联,如图4所示。
本发明用d轴和q轴阻抗的平均值来表示风机的在工频处的最终阻抗,即:
Z(s)PMSG=[(Z(s)dd+Z(s)qq)+j(Z(s)qd-Z(s)dq)]/2 (24)
代入得,
Z(s=j2πfg)=Req+jωLeq (25)
那么工频下直驱风机的等效阻抗为:
Z(s=j2πf0)=Req+jXeq=-550.93+j12971.26=-550.93+j2πf0·41.31
与直驱风电机组连接的变压器阻抗为:
Figure BDA0002443650500000143
工频下直驱风机组电机组(500台)的等效阻抗为:
Req=-550.93/500=-1.10
Leq=(23.174+41.31)/500=0.1290Η
那么直驱风电机组接入后,由于与网侧的Rg、Lg、C并联,共同组成了“新的网侧”,必然改变原有的自然振荡频率,经过新的网侧的改变,等效风电机组的R-L并联之路,与原有网侧电路进行并联组合:
Z1=Rg+j(XT-火电侧+XL+XC)
Z2=R+j(XT-风机侧+X风机) (26)
并联之后有:
Z=Z1//Z2 (27)
再将火电机组的次暂态电抗串联,可得等效的总阻抗如下:
Ztatal=Z+jX″ (28)
令该网络的总阻抗虚部为零,可得该工况下的自然扭振频率为:
fer=22.81Hz
f0-fer=50-22.81=27.19Hz
可以看到,27.19Hz与轴系模态2频率26.5相差0.69Hz,即直驱风机接入后,必然会导致整个网络的自然振荡频率发生变化,与其互补的频率与发电机轴系扭振频率发生偏移,因此在直驱风电机组并网SSO振荡频率远离发电机轴系扭振频率的工况下,直驱风电机组接入,使得火电机组发生SSO的风险变小。
在风电机组并网的SSO振荡频率距离轴系扭振较远的情况下,以上分析了风机并网对多质块火电机组的发生次同步振荡的风险是起削弱作用。其机理解释:风电机组并网SSO频率距离轴系较远,忽略了其振荡或者基于理论推导可认为,在轴系扭振频率处,风电机组阻抗可用RL阻抗来代替,由于这个等效的RL,结合网侧的RL,改变了串补度,从而影响输电系统的自然振荡频率,从而避免了可能激发火电机组的次同步振荡的危险。(4)判断当振荡频率相距轴系扭振频率大于2Hz时
对于风火组合送出系统而言,当风电机组振荡频率相距轴系扭振频率大于2Hz时,可认为风电机组并网SSO频率距离火电机组轴系扭振频率较远;此时可将风火组合电力系统划分为两个组合子系统传递函数形式,如图5,子系统1为直驱风电场经弱交流系统,子系统2为多质量块同步发电机系统。
本发明所关注的模式为直驱风电机组并网振荡频率与发电机某一轴系模态频率相近的函数极点,子系统1直驱风电机组经弱交流系统的函数极点计算结果如表1,子系统2多质量块同步发电机系统的函数极点计算结果如表2。由表1、表2可知,直驱风电机组并网的振荡模式λp2和发电机的轴系模态λm2和相互靠近,λp2≈λm2,且当各自的模式阻尼为正时,那么,其中一个子系统模式所对应的并联模式阻尼被弱化,从而可能影响系统稳定性。
基于阻抗特性的函数极点法,将多质块轴系发电机组机端dq轴阻抗特性以及直驱风电机组和外部交流电网dq轴阻抗特性并联,可以求出并联后风火组合系统的函数极点。在众多函数极点中选出本发明所关注的模式,同时考虑模式阻尼为正的情况,并联后风火组合系统结果如表3所示。
表1子系统1部分极点结果
Figure BDA0002443650500000161
Figure BDA0002443650500000171
表2子系统2部分特征值结果
Figure BDA0002443650500000172
表3风火组合系统次同步振荡模式计算
Figure BDA0002443650500000173
本发明对此时进行时域仿真:
风火组合系统中,模拟风电机组单台出力0.1MW(1000台模拟风电机组),发电机出力200MW。此时风电机组接入弱交流系统后有功曲线如图6。可以看到,在1s时刻发生了收敛的持续性振荡,通过对其功率曲线的FFT进行分析,分析图如图7。FFT分析可得,有功功率中出现了幅值较大的频率为26.57Hz的次同步振荡分量,并且该振荡频率大小和本发明发电机的模态2轴系扭振频率十分接近,依据强耦合理论,那么必然会导致其中一个子系统模式阻尼严重弱化。
通过观察此时发电机的模态分量曲线图8(图8中a为模态s21;b为模态s26),可以看出,模态1收敛性振荡,该模态下系统稳定;但系统模态2发生了不收敛的次同步振荡,激发了火电机组轴系模态2的强烈扭振现象。因此在该工况下,直驱风电机组并入后,引发了火电机组的轴系强烈扭振。
此时计算λm2下风火组合系统相对子系统模式的偏移量及子系统模式的实部如表4:
表4模式判据计算结果
Figure BDA0002443650500000181
计算出风火组合系统模式相对子系统模式的偏移量
Figure BDA0002443650500000182
Figure BDA0002443650500000183
由此可得出振荡频率相距轴系扭振频率较小时的模式判稳方法的稳定判据,当满足:
real(A)>B2 (29)
Figure BDA0002443650500000184
处于复平面右半平面,系统处于不稳定状态。
那么,带入验证
Figure BDA0002443650500000185
B2=(-1.5014)2=2.254
此时real(A)=6.125>B2
此时
Figure BDA0002443650500000186
处于复平面右半平面,系统处于不稳定状态,与上述仿真结果一致,即在风电机组并网SSO频率距离火电机组轴系扭振频率较近时工况时,λp2和发电机轴系扭振模式λm2接近重合,两子系统间交互作用急剧增强,与之对应的风火组合模式
Figure BDA0002443650500000191
Figure BDA0002443650500000192
向相反方向大幅移动,并联组合模式
Figure BDA0002443650500000193
阻尼严重弱化,风电场接入使得
Figure BDA0002443650500000194
的模式阻尼变负,从而影响系统稳定性,发生了所仿真的次同步振荡4。
因此本发明给出的一种基于模态频率差异的发电机组对直驱风电机组并网次同步振荡特性的分析方法,在建立准确的考虑多质块轴系发电机组dq轴阻抗模型的基础上,通过将直驱风电机组并网振荡频率划分为相距轴系扭振频率较远和较近两类,并基于阻抗分析法,应用传递函数的阻抗极点判据,从而能够准确量化直驱风电机组对火电机组并网SSO影响,得到发电机组对直驱风电机组并网次同步振荡影响特性。

Claims (6)

1.一种发电机组对直驱风电机组并网次同步振荡特性的分析方法,其特征在于,该方法是一种基于模态频率差异的发电机组对直驱风电机组并网次同步振荡特性的分析方法;包括如下步骤:
S1:建立多质块轴系发电机组的机端dq轴阻抗模型;
S2:基于波特图判据,判断发电机组多质块轴系扭振频率与直驱风电机组并网电气谐振频率之间的关系;
S3:如果发电机组轴系扭振频率与直驱风电机组并网电气谐振频率λm相差大于2HZ,则认为直驱风电机组并网振荡频率λm相距轴系扭振频率λp远,于是将直驱风电机组等效为dq轴下一固定RL阻抗;
S4:如果发电机组轴系扭振频率λp与直驱风电机组并网电气谐振频率相差小于2HZ,则认为直驱风电机组并网振荡频率λm相距轴系扭振频率λp近,于是将风电、火电组合系统划分为两个组合子系统,并将多质块轴系发电机组机端dq轴阻抗特性以及直驱风电机组和外部交流电网dq轴阻抗特性并联,应用阻抗函数的极点判断整个系统的稳定性。
2.根据权利要求1所述发电机组对直驱风电机组并网次同步振荡特性的分析方法,其特征在于,所述步骤1建立多质块轴系发电机组的机端dq轴阻抗模型;首先建立发电机轴系模型,对发电机组其他各个部分线性化数学模型,做一定简化,忽略汽轮机调速器影响、励磁调节,并认为励磁电压恒定;其多质块轴系的运动方程在运行点线性化后,得到:
Figure FDA0002443650490000011
其中,ωb为基准角速度;
Δω=[Δω1,Δω2,Δω3,Δω4,Δω5,Δω6]Τ,ωi(i=1,…,6)表示轴系各质量块的角速度;
Δδ=[Δδ1,Δδ2,Δδ3,Δδ4,Δδ5,Δδ6]Τ,δi(i=1,…,6)表示轴系各质量块的角位移;
Μ=diag[M1,M2,M3,M4,M5,M6],Mi(i=1,…,6)表示各质量块惯性时间常数;
ΔΤ=[ΔT1,ΔT2,ΔT3,ΔT4,ΔT5,ΔT6]Τ,Ti(i=1,…,6)表示各质量块上机械力矩,其中,质量块5对应发电机转子,ΔT5=-ΔTe,ΔTe为发电机电磁转矩,质量块6对应励磁机转子,取ΔT6=0;T表示各质量块上机械力矩,p为微分算子;
其中,Dii(i=1,…,6)为轴系各质量块自阻尼系数,Dij(i=1,…,6)为相邻质量块i和j之间互阻尼系数;Kij(i=1,…,6)为质量块i和j之间的弹性系数。
3.根据权利要求1所述发电机组对直驱风电机组并网次同步振荡特性的分析方法,其特征在于,所述步骤2基于波特图判据,基于风电机组阻抗和电网阻抗频率特性曲线的波特图判据,通过解析法或者实测法获得风电机组和电网的阻抗-频率波特图,大致找出直驱风机并网发生的振荡频率,从而通过分析风电机组阻抗幅频曲线与电网阻抗幅频曲线交点频率对应的二者相位差
Figure FDA0002443650490000021
判断系统稳定性,从而定位系统可能发生的振荡区域,并通过判断直驱风机并网振荡频率差是否大于2Hz:进一步将直驱风机并网振荡频率λm划分为相距轴系扭振频率λp远和近两类来进行影响分析,具有较高稳定性;若
Figure FDA0002443650490000022
系统稳定,反之,系统不稳定。
4.根据权利要求1所述发电机组对直驱风电机组并网次同步振荡特性的分析方法,其特征在于,所述如果发电机组轴系扭振频率与直驱风电机组并网电气谐振频率相差小于2HZ,则认为直驱风机并网振荡频率λm相距轴系扭振频率λp近,于是将风火组合系统划分为两个组合子系统,并将多质块轴系发电机组机端dq轴阻抗特性以及直驱风电机组和外部交流电网dq轴阻抗特性并联,应用阻抗函数的极点判断整个系统的稳定性。
5.根据权利要求1所述发电机组对直驱风电机组并网次同步振荡特性的分析方法,其特征在于,发电机的次同步谐振实质上是与自然振荡频率的互补频率,接近于发电机轴系的模态2频率(26.56Hz),这一振荡趋势增大,激发了模态2的强烈的轴系扭振,从而形成机械系统与电气系统之间的相互激励,在系统总阻尼为负或者很小的情况下,相互激励引起的振荡难以维持导致系统失稳。
6.根据权利要求1所述一种基于模态频率差异的发电机组对直驱风电场并网次同步振荡特性的分析方法,其特征在于,所述步骤3通过将直驱风机等效为dq轴下一固定RL阻抗:
Z(s)PMSG=[(Z(s)dd+Z(s)qq)+j(Z(s)qd-Z(s)dq)]/2
从而量化分析了发电机组对直驱风电场并网次同步振荡特性。
CN202010272580.5A 2020-04-09 2020-04-09 发电机组对直驱风电机组并网次同步振荡特性的分析方法 Active CN111525549B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010272580.5A CN111525549B (zh) 2020-04-09 2020-04-09 发电机组对直驱风电机组并网次同步振荡特性的分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010272580.5A CN111525549B (zh) 2020-04-09 2020-04-09 发电机组对直驱风电机组并网次同步振荡特性的分析方法

Publications (2)

Publication Number Publication Date
CN111525549A true CN111525549A (zh) 2020-08-11
CN111525549B CN111525549B (zh) 2023-11-14

Family

ID=71902490

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010272580.5A Active CN111525549B (zh) 2020-04-09 2020-04-09 发电机组对直驱风电机组并网次同步振荡特性的分析方法

Country Status (1)

Country Link
CN (1) CN111525549B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112186759A (zh) * 2020-09-28 2021-01-05 西安热工研究院有限公司 自适应捕捉频点的双馈风电场次同步振荡抑制方法
CN113364005A (zh) * 2021-05-20 2021-09-07 国网冀北电力有限公司电力科学研究院 风电振荡激发汽轮机组轴系扭振风险的监测方法及装置
CN113612238A (zh) * 2021-07-20 2021-11-05 国网江苏省电力有限公司电力科学研究院 一种分析风电并网引发火电机组轴系小干扰振荡的改进阻尼转矩方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101615791A (zh) * 2009-07-29 2009-12-30 清华大学 有效抑制次同步谐振和振荡的次同步阻尼控制系统
CN104638675A (zh) * 2015-02-10 2015-05-20 北京四方继保自动化股份有限公司 发电机机端次同步阻尼控制器参数整定方法
US20160218650A1 (en) * 2015-01-23 2016-07-28 Rolls-Royce Plc Method and system for damping torsional oscillations
CN106599526A (zh) * 2017-02-17 2017-04-26 华北电力大学(保定) 一种火电机组次同步振荡的阻抗分析法
CN110148967A (zh) * 2019-06-18 2019-08-20 华北电力大学 一种基于导纳分析直驱风机次同步振荡特性的研究方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101615791A (zh) * 2009-07-29 2009-12-30 清华大学 有效抑制次同步谐振和振荡的次同步阻尼控制系统
US20160218650A1 (en) * 2015-01-23 2016-07-28 Rolls-Royce Plc Method and system for damping torsional oscillations
CN104638675A (zh) * 2015-02-10 2015-05-20 北京四方继保自动化股份有限公司 发电机机端次同步阻尼控制器参数整定方法
CN106599526A (zh) * 2017-02-17 2017-04-26 华北电力大学(保定) 一种火电机组次同步振荡的阻抗分析法
CN110148967A (zh) * 2019-06-18 2019-08-20 华北电力大学 一种基于导纳分析直驱风机次同步振荡特性的研究方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘丽萍;: "考虑轴系扭振的大型同步发电机误并列运行仿真分析", 新型工业化 *
董德勇;宿端鹏;: "次同步振荡误差评估与建模", 电气技术 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112186759A (zh) * 2020-09-28 2021-01-05 西安热工研究院有限公司 自适应捕捉频点的双馈风电场次同步振荡抑制方法
CN112186759B (zh) * 2020-09-28 2022-11-15 西安热工研究院有限公司 自适应捕捉频点的双馈风电场次同步振荡抑制方法
CN113364005A (zh) * 2021-05-20 2021-09-07 国网冀北电力有限公司电力科学研究院 风电振荡激发汽轮机组轴系扭振风险的监测方法及装置
CN113364005B (zh) * 2021-05-20 2023-09-26 国网冀北电力有限公司电力科学研究院 风电振荡激发汽轮机组轴系扭振风险的监测方法及装置
CN113612238A (zh) * 2021-07-20 2021-11-05 国网江苏省电力有限公司电力科学研究院 一种分析风电并网引发火电机组轴系小干扰振荡的改进阻尼转矩方法

Also Published As

Publication number Publication date
CN111525549B (zh) 2023-11-14

Similar Documents

Publication Publication Date Title
CN111525549A (zh) 发电机组对直驱风电机组并网次同步振荡特性的分析方法
CN102150356B (zh) 具有分量分离的直接功率控制
CN102157933B (zh) 一种电力系统稳定器设计方法
CN106599526B (zh) 一种火电机组次同步振荡的阻抗分析法
Hamane et al. Control for variable speed wind turbine driving a doubly fed induction generator using Fuzzy-PI control
CN111555312B (zh) 一种适用于电力系统超低频振荡稳定性评估的方法
Li et al. Impact of increased wind power generation on subsynchronous resonance of turbine-generator units
CN106533289B (zh) 一种非线性电压控制方法及系统
CN110611321B (zh) 一种补偿虚拟同步机负阻尼特性的虚拟电力系统稳定器设计方法
CN115882762A (zh) 构网型风电并网系统的频率优化控制方法
CN111740449B (zh) 风电机群对同步发电机组阻尼特性分析方法
CN107069697A (zh) 基于轨迹特征根的系统暂态稳定判断和扰动类型筛选方法
CN110417047B (zh) 基于复转矩系数分析双馈风机ssci阻尼特性的方法
CN113612238A (zh) 一种分析风电并网引发火电机组轴系小干扰振荡的改进阻尼转矩方法
CN112886611B (zh) 一种直驱风机并网系统的次同步振荡抑制方法
CN109713664B (zh) 直流孤岛频率稳定的网源协调控制策略计算方法及系统
CN114792055A (zh) 一种基于暂态电抗后电势的异步电动机等效惯量评估方法
CN112464391A (zh) 一种柴油发电机组哈密顿模型的建立方法
Zhang et al. Research of coordination control system between nonlinear robust excitation control and governor power system stabilizer in multi-machine power system
Cidrás et al. Regulation of synchronous generators by means of hydrostatic transmissions
Swe et al. Modelling and simulation for small signal stability of multi-machine power system under various disturbance conditions
Alomari et al. Modeling and simulation for bifurcations of SSR in large wind farms
CN110212582B (zh) 一种风电参与电力系统调频调压的仿真方法
Yuan et al. Impedance-based sub-synchronous torsional stability analysis in a HVDC system connected to multiple thermal power plants
CN115714401A (zh) 一种电力系统超低频率振荡风险分析方法、设备及介质

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant