双馈风电机组对电网电压振荡的载荷响应特性*

马骏超，杨 靖，彭 琰，孙 勇，王晨旭，花 斌†，孟军峰

双馈风电机组对电网电压振荡的载荷响应特性*

马骏超1，杨 靖2,3，彭 琰1，孙 勇2,3，王晨旭1，花 斌2,3†，孟军峰2

（1. 国网浙江省电力公司电力科学研究院，杭州 310014；2. 浙江运达风电股份有限公司，杭州 311106；3. 浙江省风力发电技术重点实验室，杭州 311106）

电网电压振荡会引起发电机电磁力矩振荡，增加传动链轴系和塔架左右方向载荷，给传动链和塔架带来一定的疲劳损伤。通过分析发电机电磁力矩振荡频率与电网电压振荡频率的关系，考察传动链轴系载荷、塔架左右方向载荷对发电机电磁力矩振荡的响应特性，从而建立联合仿真模型，对比分析不同电磁力矩振荡频率对传动链轴系载荷、塔架左右方向载荷的影响。结果表明，电网基波频率附近电压振荡会引起发电机电磁力矩低频振荡，进而显著增加传动链轴系和塔架左右方向载荷。

电压振荡；电磁力矩振荡；传动链轴系；塔架；载荷响应；联合仿真

0 引 言

风电是重要的可再生能源，加快风电发展对实现双碳目标和构建以新能源为主体的新型电力系统具有十分重要的意义[1]。风电一般采用电力电子变流器并网，电力电子设备与电网间的相互作用，易引发电力系统振荡[2]，如美国得克萨斯州某风电场双馈风电集群与串补输电线路发生频率约为25 Hz的次同步振荡[3]；我国新疆哈密市直驱风电机群发生频率为10 ～ 40 Hz的次同步振荡[4-5]，引起临近汽轮机组轴系扭振[4]；德国北海风电−柔性直流输电发生中高频段振荡现象等[6]。可见，大规模风电并网引发的宽频振荡问题大多为风电机组之间或风电机组与串补线路/直流输电线路间相互作用引起，振荡频率范围涉及几赫兹至几千赫兹。宽频振荡问题涉及多类型设备，多时间尺度，振荡信号呈现强时变特性、强非线性、多模态以及时空分布特性的特点[7]。

目前，针对高比例新能源并网引起的电力系统宽频振荡问题，国内外学者大多聚焦于振荡机理分析、振荡模态分析及监测和振荡抑制等方面。文献[8-10]综述了国内外因高比例新能源并网引起的宽频振荡机理，文献[2]着重分析国内实际振荡案例，指明下阶段研究方向。实际工程中频谱分析一般采用傅里叶变换[11]，近年来，基于希尔伯特−黄变换[12]、小波分析[13]、人工智能算法[14]等模态分析方法都取得较快发展，但复杂的计算过程限制了其在实际工程中的应用。文献[15-16]提出宽频振荡广域监测方法，但短期内无法实现产业化推广。振荡抑制可在早期电网规划阶段通过理论分析和半物理仿真来避免，此外风电机组附加阻尼控制[17-18]、风电场动态无功补偿装置及静止无功补偿器等装置附加阻尼控制[19-20]、附加滤波器[21]等措施均可以实现宽频振荡抑制。

可见，目前对于新能源并网后引起的宽频振荡问题的研究大多从电网的角度入手，且短时期内无法完成宽频振荡的快速监测和振荡抑制的产业化推广，因此从风电机组的角度出发，研究风电机组在电力系统振荡时机械部件的载荷响应特性对保障机组安全运行具有十分重要的现实意义。

针对上述问题，本文首先分析发电机电磁力矩振荡频率与电网电压振荡频率的关系，然后分析大型双馈异步风电机组传动链轴系载荷、塔架左右方向载荷对发电机电磁力矩的响应特性，最后搭建GH Bladed与Matlab/Simulink联合仿真模型，研究不同电磁力矩振荡频率下传动链轴系载荷和塔架左右方向载荷变化趋势。

1 发电机电磁力矩振荡频率与系统电压振荡频率的关系

对于双馈风电机组而言，在任意速旋转的坐标系中，电网电压轴分量与定子磁链轴分量的关系为：

式中：sr、si为电网电压轴分量；sr、si为定子磁链轴分量；sr、si为定子电流轴分量；s为定子电阻。

稳态时，定子磁链的变化率为0，再忽略定子电阻，则式（1）可简化为：

在不考虑扰动时，发电机电磁转矩表达式为：

式中：p为发电机极对数；m为定转子同轴等效绕组间的互感；sd、sq为定子电流轴分量；rd、rq为转子电流轴分量。

假设电网存在频率为d的电压扰动Δ，且扰动在坐标系上的分量分别为Δsd和Δsq，其相量图如图1所示。图中，()s坐标系为两相静止坐标系；坐标系为同步速s旋转坐标系。

图1 扰动电压相量图

扰动分量在坐标轴上的分量分别为：

式中：|Δ|为扰动电压的幅值；为Δ与轴的夹角，0为其初始夹角，与扰动前电网运行状态相关。

由式（2）可知，在定子磁链暂态分量衰减结束后，由电压扰动Δsd、Δsq引起的磁链扰动在轴上的分量Δsd、Δsq为：

由式（4）和式（5）可以看出，频率为d的扰动电压，将在两相同步速旋转的轴坐标系上感应出频率为|d−s|的电压扰动和磁链扰动分量。

由双馈异步风电机组的定子磁链方程可知：

式中：s为坐标系中定子等效两相绕组自感。

因此，在轴有频率为|d−s|的定子电流扰动Δsd、Δsq和频率为|d−s|的转子电流扰动Δrd、Δrq，假设：

式中：1和2的频率均为|d−s|，初始角度与扰动前电网运行状态相关。1、2分别为Δs、Δr的幅值。

考虑定转子电流的扰动分量时，电磁转矩为：

将式（7）、式（8）代入式（9），可知：

因此，电磁力矩含有三种分量，第一种分量为pm(sqrd−sdrq)，这是正常运行时的直流分量；第二种分量为pm(rdΔsq−sdΔrq+sqΔrd−rqΔsd)，这是电压扰动导致的交流振荡分量，振荡频率为|d−s|；第三种分量为pm12sin(1−2)，这是电压扰动导致的直流分量，这部分幅值较小，一般可忽略。

综上所述，频率为d的电压振荡会产生频率为|d−s|的电磁力矩振荡。

2 风电机组传动链和塔架对电磁力矩扰动响应特性

电力系统振荡会导致发电机电磁力矩振荡，而发电机电磁力矩振荡会激励传动链轴系和塔架左右方向某阶模态发生共振，进而产生较大载荷，对传动链和塔架造成一定的疲劳损伤。本节重点分析传动链轴系载荷、塔架左右方向载荷对发电机电磁力矩振荡的响应特性。

2.1 传动链轴系响应特性

大型双馈风电机组传动链轴系通常采用等效质量法结合动力学方程来进行建模，不同等效方法可将传动链系统等效为不同的质量块模型，目前常见的有六质量块模型、四质量块模型、三质量块模型、两质量块模型和单质量块模型。由于两质量块模型能反映传动链轴系扭振动态特性且精度能满足研究需要[21-22]，因此选取两质量块模型为研究对象，各物理量均折算到低速轴侧，可得传动链轴系两质量块模型的动力学方程：

其中：r、g分别为风轮与发电机的转动惯量；s为传动链轴系等效阻尼；s为传动链轴系等效刚度；r、g、s分别为风轮、发电机与传动链轴系的角位移；r、g、s分别为风轮、发电机与传动链轴系的转矩。

消除式（12）中的r、g可知：

根据式（13）可以得出传动链轴系固有振荡频率和阻尼系数为：

传动链轴系载荷对发电机电磁力矩频率响应如图2所示，结合WD 107-2500型双馈风电机组关键参数[14]可以得出c= 9.41 rad/s，c= 0.0058。可以看出传动链轴系响应是一个典型的二阶欠阻尼系统，电磁力矩扰动易激发传动链轴系扭振。

图2 传动链轴系载荷对发电机电磁力矩频率响应

2.2 塔架左右方向响应特性

发电机电磁力矩是气动力矩的反力矩，同时也作为弯矩作用于塔架左右振动模态。根据拉格朗日方程可以得到：

根据式（16）可以得出塔架左右方向固有振荡频率和阻尼系数为：

塔架左右方向载荷对发电机电磁力矩频率响应如图3所示，结合WD 107-2500型双馈风电机组关键参数[22]可以得出ωt= 1.92 rad/s，ξt= 0.0064。可以看出塔架左右方向响应是一个典型的二阶欠阻尼系统，电磁力矩扰动易激发塔架左右方向晃动。

3 仿真分析

现有商业仿真软件很难同时满足计及载荷特性和电网特性的精细化建模（体现机电交互耦合特性）需求。GH Bladed软件是目前国内外整机制造商应用最广的商业化风机设计软件，其优势在于整机建模仿真、载荷计算、振动分析等，但要建立比较精细化的发电机、变流器、电网模型是比较困难的。Maltab软件是研究双馈风电机组的并网特性和电力系统仿真计算方面应用最广泛的仿真软件，其优势在于发电机、变流器、电网模型等并网特性仿真研究，但是在Matlab软件中构建风轮、塔架详细模型、载荷计算、振动分析等比较困难。因此，计及载荷特性和电网特性的精细化模型需要采用联合仿真技术，联合仿真模型结构如图4所示。

图4 联合仿真模型

3.1 电压振荡仿真分析

电力系统正常运行时，并网点（point of common coupling, PCC）电压本身存在不同频率的振荡分量，但幅值很小，可以忽略。但是负荷投切、电源出力的变化等扰动会加剧并网点电压的振荡。本节以投入新负荷、风电场参与电力系统调频为例来分析并网点电压基波附近振荡分量的变化趋势。仿真模型的电气部分采用如图5所示含风电场的四机两区模型[23]。图中，G1～ G4表示常规同步机；、分别表示负荷等效电感、电容；WF表示风电场，容量为300 MW；dist表示扰动负荷，位于节点9，容量为180 MW。

图5 含风电场的四机两区模型

仿真选取平均风速为14 m/s湍流风况，对扰动负荷投入前后的并网点电压利用离散傅里叶变换进行频谱分析，频谱分析结果如图6所示。扰动前后并网点电压振荡分量如表1所示。

图6 并网点电压频谱分析

表1 扰动前后并网点电压振荡分量

表1中频率偏差表示偏离基波的程度，如±2 Hz表示48 Hz与52 Hz的振荡分量，百分比表示该振荡分量与基波的比值。仿真结果表明，扰动后基波附近的振荡分量占比较扰动前均有明显上升，此部分振荡分量将会引起发电机电磁力矩的低频振荡，破坏轴系平衡，进而对风电机组关键部件的载荷造成一定的影响。

3.2 不同电磁力矩振荡频率风电机组载荷响应特性

仿真选取平均风速为14 m/s湍流风况，发电机电磁力矩振荡是由于电网电压振荡引起的，振荡幅值为5%n（n为额定电磁力矩），持续时间20 s，振荡频率分别是塔架左右一阶模态频率、传动链一阶模态频率、桨叶面内三阶模态频率、塔架左右二阶模态频率和高频（25 Hz、70 Hz）。疲劳载荷为用雨流计数法[24]处理之后的当量等效疲劳载荷（其中S-N斜率曲线取10）。

双馈风电机组仿真机型为WD107-2500[14]，该机型关键模态的频率和阻尼如表2所示。

表2 WD107-2500机型关键模态频率和阻尼系数

（1）振荡频率为塔架左右一阶模态频率

当电磁力矩振荡频率为塔架左右一阶模态频率（0.30 Hz）时，仿真结果如图7所示。

当振荡频率为塔架左右一阶模态频率时，塔架左右方向最大载荷增加6%，疲劳载荷增加34%，传动链轴系最大载荷不变，疲劳载荷增加23%。可见塔架左右方向振动较明显，传动链轴系也出现一定程度的扭振。

图7 振荡频率为塔架左右一阶模态频率时载荷响应特性

（2）振荡频率为传动链一阶模态频率

当电磁力矩振荡频率为传动链一阶模态频率（1.34 Hz）时，仿真结果如图8所示。当振荡频率为传动链一阶模态频率时，传动链轴系最大载荷增加8%，疲劳载荷增加172%；塔架左右方向最大载荷不变，疲劳载荷增加3%。可见传动链轴系疲劳载荷有明显增加，可能对传动荷造成较大疲劳损伤，应当引起足够重视。此外，发电机转速也有一定程度的振荡。

图8 振荡频率为传动链一阶模态频率时载荷响应特性

（3）振荡频率为桨叶面内三阶模态频率

当电磁力矩振荡频率为桨叶面内三阶模态频率（2.36 Hz）时，仿真结果如图9所示。振荡频率为桨叶面内三阶模态频率时，传动链轴系最大载荷不变，疲劳载荷增加36%；塔架左右方向最大载荷不变，疲劳载荷增加6%。可见传动链轴系出现一定程度的扭振。

（4）振荡频率为塔架左右二阶模态频率

当电磁力矩振荡频率为塔架左右二阶模态频率（2.71 Hz）时，仿真结果如图10所示。振荡频率为塔架左右二阶模态频率（2.71 Hz）时，传动链轴系最大载荷增加5%，疲劳载荷增加98%；塔架左右方向最大载荷增加28%，疲劳载荷增加103%。可见传动链轴系以及塔架左右方向疲劳载荷均有明显增加，传动链以及塔架左右方向振动明显，容易造成疲劳损伤，应当采取相应措施进行载荷优化。

（5）振荡频率为高频（25 Hz）

当电磁力矩振荡频率为高频（25 Hz）时，仿真结果如图11所示。振荡频率为高频（25 Hz）时，发电机电磁力矩产生高频振荡，但对风电机组传动链轴系和塔架左右方向的载荷影响很小，几乎没有变化。对于机械振动而言，25 Hz、70 Hz都属于高频的范畴，限于篇幅原因，本工况下仅给出振荡频率为25 Hz时发电机电磁力矩与传动链轴系载荷的仿真图。

对比各种电磁力矩振荡频率时，传动链轴系以及塔架左右方向最大载荷和疲劳载荷，其结果如表3所示。

表3 不同电磁力矩振荡频率最大载荷比与疲劳载荷比

从表3可以看出，对传动链轴系载荷影响较大的模态频率主要有传动链一阶模态，塔架左右方向二阶模态和桨叶面内三阶模态；对塔架左右方向载荷影响较大的模态频率主要有塔架左右方向一阶模态和塔架左右方向二阶模态。这与第2节中的分析一致。

对比表3中的数据可知，当电磁力矩振荡频率为塔架左右方向二阶模态频率（2.71 Hz）时，传动链轴系以及塔架左右方向载荷均有显著增加，给传动链、塔架造成的疲劳损伤最大，应当采取相应的措施抑制传动链轴系、塔架左右方向扭振，实现载荷优化。此外，当电磁力矩振荡频率为传动链一阶模态频率（1.34 Hz）时，传动链轴系载荷增加最大，会造成传动链轴系扭振，增加其疲劳损伤，同样需要引起重视。

4 结 论

通过研究双馈风电机组不同电磁力矩振荡频率载荷响应特性，得出以下结论：

（1）双馈风电机组定子侧直接并网，电网工频附近的电压振荡会引起发电机电磁力矩低频振荡。

（2）发电机电磁力矩低频振荡会增加风电机组机械部件载荷，如当发电机电磁力矩振荡频率为传动链一阶模态频率（1.34 Hz）时，传动链轴系载荷增加172%；发电机电磁力矩振荡频率为塔架左右方向二阶模态频率时，传动链轴系载荷增加98%，塔架左右方向载荷增加103%。

（3）本文仅选取WD107-2500型双馈风电机组的几个关键模态频率作为发电机电磁力矩振荡频率，但这并不意味着只有这几个模态频率才会激发风电机组关键部件的载荷变化，任意由电网电压扰动引起的发电机低频电磁力矩振荡均会导致风电机组关键部件载荷变化。

（4）发电机电磁力矩振荡频率为高频时，对风电机组传动链、塔架的载荷几乎没有影响。

本文仅从理论分析和联合仿真分析阐明不同的发电机电磁振荡频率对风电机组传动链、塔架载荷的影响，目前已有的基于附加阻尼控制的载荷优化方案大多处于理论研究阶段，后期应注重开展样机测试，完成风电机组的改造，实现风电机组振荡抑制功能。

[1] 舒印彪,张丽英,张运洲,等.我国电力碳达峰、碳中和路径研究[J].中国工程科学,2021,23(06):1-14. DOI:10.15302/J-SSCAE-2021.06.001.

[2] 李明节, 于钊, 许涛, 等. 新能源并网系统引发的复杂振荡问题及其对策研究[J]. 电网技术, 2017, 41(4): 1035-1042. DOI: 10.13335/j.1000-3673.pst.2016.3049.

[3] ADAMS J, PAPPU V A, DIXIT A. ERCOT experience screening for sub-synchronous control interaction in the vicinity of series capacitor banks[C]//Proceedings of 2012 IEEE Power and Energy SocietyGeneral Meeting. San Diego, CA: IEEE, 2012: 1-5. DOI: 10.1109/PESGM.2012.6345409.

[4] 谢小荣, 刘华坤, 贺静波, 等. 直驱风机风电场与交流电网相互作用引发次同步振荡的机理与特性分析[J]. 中国电机工程学报, 2016, 36(9): 2366-2372. DOI: 10.13334/j.0258-8013.pcsee.2016.09.007.

[5] LIU H K, XIE X R, HE J B, et al. Subsynchronous interaction between direct-drive PMSG based wind farms and weak AC networks[J]. IEEE transactions on power systems, 2017, 32(6): 4708-4720. DOI: 10.1109/tpwrs.2017.2682197.

[6] BUCHHAGEN C, RAUSCHER C, MENZE A, et al. BorWin1-first experiences with harmonic interactions in converter dominated grids[C]//Proceedings of International ETG Congress 2015, Die Energiewende-Blueprints for the New Energy Age. Bonn, Germany: VDE, 2015: 1-7.

[7] 赵妍, 唐文石, 聂永辉, 等. 基于格拉姆角差场和卷积神经网络的宽频振荡分类方法[J/OL]. 电网技术: 1-12[2022-04-14]. 10.13335/j.1000-3673.pst.2021.2079.https://doi.org/10.13335/j.1000-3673.pst.2021.2079.

[8] 刘芳, 刘威, 汪浩东, 等. 高比例新能源电力系统振荡机理及其分析方法研究综述[J]. 高电压技术, 2022, 48(1): 95-113. DOI: 10.13336/j.1003-6520.hve.20211187.

[9] 马宁宁, 谢小荣, 贺静波, 等. 高比例新能源和电力电子设备电力系统的宽频振荡研究综述[J]. 中国电机工程学报, 2020, 40(15): 4719-4731. DOI: 10.13334/j.0258-8013.pcsee.191968.

[10] 王伟胜, 张冲, 何国庆, 等. 大规模风电场并网系统次同步振荡研究综述[J]. 电网技术, 2017, 41(4): 1050-1060. DOI: 10.13335/j.1000-3673.pst.2017.0069.

[11] 马超, 王慧铮, 蒙永苹, 等. 低频振荡综合监测算法研究[J]. 电网技术, 2008, 32(S2): 31-35. DOI:10.13335/j.1000-3673.pst.2008.s2.072.

[12] 韩松, 何利铨, 孙斌, 等. 基于希尔伯特-黄变换的电力系统低频振荡的非线性非平稳分析及其应用[J]. 电网技术, 2008, 32(4): 56-60. DOI:10.13335/j.1000-3673.pst.2008.04.013.

[13] 冯双, 崔昊, 陈佳宁, 等. 人工智能在电力系统宽频振荡中的应用与挑战[J]. 中国电机工程学报, 2021, 41(23): 7889-7904. DOI: 10.13334/j.0258-8013.pcsee.210608.

[14] 马俊杰, 刘芳, 吴敏, 等. 基于类噪声小波分解的风电场次同步振荡辨识[J]. 电网技术, 2019, 43(4): 1294-1300. DOI: 10.13335/j.1000-3673.pst.2018.1468.

[15] 吴琛, 马宁宁, 刘旭斐, 等. 高比例新能源电力系统多模态振荡监测方法及装置设计[J]. 电网技术, 2021, 45(9): 3496-3504. DOI: 10.13335/j.1000-3673.pst.2020.1902.

[16] 鞠平, 谢欢, 孟远景, 等. 基于广域测量信息在线辨识低频振荡[J]. 中国电机工程学报, 2005, 25(22): 56-60. DOI: 10.3321/j.issn:0258-8013.2005.22.010.

[17] VIETO I, SUN J. Damping of subsynchronous resonance involving Type-III wind turbines[C]//Proceedings of the IEEE 16th Workshop on Control and Modeling for Power Electronics (COMPEL). Vancouver, BC, Canada: IEEE, 2015: 1-8. DOI: 10.1109/COMPEL.2015.7236515.

[18] 谢震, 刘坤, 张兴, 等. 双馈风力发电机在电网电压不对称骤升下无功功率优化控制[J]. 中国电机工程学报, 2015, 35(13): 3211-3220. DOI:10.13334/j.0258-8013. pcsee.2015.13.003.

[19] EL-MOURSI M S, BAK-JENSEN B, ABDEL-RAHMAN M H. Novel STATCOM controller for mitigating SSR and damping power system oscillations in a series compensated wind park[J]. IEEE transactions on power electronics, 2010, 25(2): 429-441. DOI: 10.1109/TPEL.2009.2026650.

[20] 孙贤大, 刘国辉, 吕孝国, 等. 基于SVG的次同步振荡抑制方法研究与应用[J]. 电力电子技术, 2020, 54(3): 75-78.

[21] LIU H K, XIE X R, LI Y, et al. Damping subsynchronous resonance in series-compensated wind farms by adding notch filters to DFIG controllers[C]//Proceedings of 2015 IEEE Innovative Smart Grid Technologies - Asia (ISGT ASIA). Bangkok, Thailand: IEEE, 2016. DOI: 10.1109/ISGT-Asia.2015.7386978.

[22] 应有, 孙勇, 杨靖, 等. 大型双馈风电机组电网故障穿越过程载荷特性分析[J]. 电力系统自动化, 2020, 44(12): 131-138. DOI: 10.7500/AEPS20190911006.

[23] KUNDUR P. Power system stability and control[M]. New York: McGraw-Hill, 1994.

[24] 徐敏敏, 杨春兰, 段月星, 等. 基于雨流计数法的疲劳寿命预测[J]. 机械设计与研究, 2016, 32(5): 184-187. DOI: 10.13952/j.cnki.jofmdr.2016.0217.

Load Response Characteristics of Doubly-Fed Induction Generator to Power System Voltage Oscillation

MA Jun-chao1, YANG Jing2,3, PENG Yan1, SUN Yong2,3,WANG Chen-xu1, HUA Bin2,3, MENG Jun-feng2

(1. State Grid Zhejiang Electric Power Research Institute, Hangzhou 310014, China;2. Zhejiang Windey Co., Ltd., Hangzhou 311106, China;3. Key Laboratory of Wind Power Technology of Zhejiang Province, Hangzhou 311106, China)

Grid voltage oscillates will cause the electromagnetic torque of the generator to oscillate, which will increase the load of the transmission chain shaft and the tower load in the left and right directions, and bring fatigue damage to the transmission chain and the tower. In this work, the relationship between the generator electromagnetic torque oscillation frequency and the grid voltage oscillation frequency was analyzed, and the response characteristics of the transmission chain shaft load and the tower load in the left and right directions to the generator electromagnetic torque oscillation was investigated. Then a co-simulation model was established to compare the effects of different electromagnetic torque oscillation frequencies on the load of the transmission chain shaft and the load in the left and right directions of the tower. Results showed that the voltage oscillation near the fundamental frequency of the power grid would cause the low-frequency oscillation of the generator electromagnetic torque, which would significantly increase the transmission chain shaft load and tower load in the left and right directions.

voltage oscillation; electromagnetic torque oscillation; transmission chain shaft; tower; load response; co-simulation

2095-560X（2022）03-0271-09

TK83；TM614

A

10.3969/j.issn.2095-560X.2022.03.012

2022-03-29

2022-04-18

国网浙江省电力有限公司科技项目（B311DS21000G）

花 斌，E-mail：huab@chinawindey.com

马骏超（1989-），男，博士，高级工程师，主要从事新能源并网运行与消纳、直流配电网运行调控关键技术研究。

花 斌（1994-），男，硕士，主要从事新能源并网运行与消纳研究。