摘要
为准确描述构网型直驱风机的动态特性,并分析其在不同电网强度下的稳定性,考虑风机源端特性对系统进行小信号建模。首先,基于状态空间建模方法,详细构建了基于虚拟同步机控制的构网型直驱风机全系统小信号模型,并通过仿真测试验证了所建小信号模型的准确性。其次,通过特征模态分析,研究了不同电网强度对系统稳定性的影响,评估其在弱电网环境下的适应性。同时,分析风速阶跃扰动对系统低频振荡的诱发机制,探讨系统结构参数(如直流侧电容)及控制器参数对系统动态性能的影响。最后,基于 RT-LAB 半实物仿真平台搭建系统模型,验证了模型在实时仿真环境下的响应能力和准确性。结果表明,构网型直驱风机在弱电网条件下具有较强的适应性,但在风速阶跃扰动下易诱发低频振荡。此外,合理调整系统结构参数和优化控制器参数可有效抑制低频振荡,提高系统在弱电网环境下的稳定性。本研究为构网型风电机组的动态特性分析提供了一种系统性的建模方法,并为进一步优化其稳定控制策略提供了理论支持。
Abstract
To accurately describe the dynamic characteristics of grid-forming direct-drive wind turbines and analyze their stability under different grid strengths, this study considers the source-side characteristics of wind turbines and conducts small-signal modeling of the system. First, based on the state-space modeling method, a comprehensive small-signal model of the entire grid-forming direct-drive wind turbine system with virtual synchronous generator (VSG) control is constructed, and its accuracy is verified through simulation tests. Second, eigenvalue-based modal analysis is employed to investigate the impact of different grid strengths on system stability and assess its adaptability in weak grid environments. Additionally, the mechanism of low-frequency oscillations induced by wind speed step disturbances is analyzed, along with the effects of system structural parameters ( such as DC-side capacitance) and controller parameters on the system’s dynamic performance. Finally, an RT-LAB hardware-inthe-loop simulation platform is used to build the system model, validating its response capability and accuracy in real-time simulation environments. The results indicate that grid-forming direct-drive wind turbines exhibit strong adaptability in weak grid conditions but may prompt low-frequency oscillations under wind speed step disturbances. Moreover, properly adjusting system structural parameters and optimizing controller parameters can effectively suppress low-frequency oscillations and enhance system stability in weak grid environments. This study provides a systematic modeling approach for analyzing the dynamic characteristics of grid-forming wind turbines and offers theoretical support for further optimizing their stability control strategies.
随着电网结构由以同步发电机为主向新能源发电和电力电子设备为核心转变,电网逐渐呈现出高新能源渗透率和高电力电子设备比例的“双高” 特性[1-3]。这一转型对系统的稳定性提出了更高要求。构网型控制技术,特别是虚拟同步机( virtual synchronous generator,VSG)控制,以其显著的功率同步和惯性同步特性,为满足现代电力系统中新能源并网的需求提供了重要支持[4-8]。然而,传统的构网型模型通常将直流侧简化为理想电压源,难以全面表征风电、光伏等动态复杂的分布式电源特性。
作为主流新能源机组,直驱风机( direct-drive wind turbine,DDWT)在采用虚拟同步控制实现并网后,其稳定性问题成为关键研究课题之一[9-10]。针对构网型逆变器并网系统,文献[11-15]通过小信号建模和特征模态分析揭示了 VSG 控制的频率调节特性,并分析了参数变化对稳定性的影响。文献[16]进一步研究了 VSG 多机并网的低频振荡特性。然而,这些研究多假设源端为理想直流电压源,忽略了新能源源端动态特性对系统稳定性的影响。
针对构网型 DDWT 的研究较少,文献[17]通过频域分析建立了构网型 DDWT 的输出阻抗模型,验证了模型的正确性,但未从时域角度深入探讨源端扰动的动态影响。文献[18-19]将 VSG 技术引入 DDWT 系统,构建了小信号模型,并分析了不同参数对系统稳定性的影响,但未涉及风速扰动时参数动态变化对直流母线电压波动的影响。文献[20] 研究了 VSG 直驱风电机组并联运行时的低频特征模态稳定性,但未充分考虑源端动态特性对系统耦合振荡的作用。
综上所述,当新能源设备作为动态复杂的源端与电网耦合时,其运行状态对系统动态稳定性具有显著影响。因此,亟须构建考虑源端特性的构网型 DDWT 并网系统小信号模型,以揭示关键参数对系统稳定性的影响,并为系统优化设计提供理论支持。
1 构网型 DDWT 的数学模型
基于 VSG 的构网型直驱风机的结构如图1所示,主要由风机、永磁同步发电机、全功率变换器,以及控制系统构成。
图1构网型直驱风机结构框图
Fig.1Schematic diagram of the grid-forming direct-drive wind turbine structure
1.1 DDWT 的数学模型
电流和电压在 DDWT 内部的相互作用通过引入交流源和串联阻抗来描述,对应的数学表达式如下:
(1)
式中:ωm 为机械角速度,umd 、umq分别为定子电压的 d、q 轴分量,imd 、imq分别为定子电流的 d、q 轴分量,lmd、lmq分别为定子电感的 d、q 轴分量,rs 为定子电阻,Φf 为转子磁通,P 为极对数。
DDWT 转子运动方程为
(2)
式中:Jw为风机自身转动惯量,ρ 为空气密度,R 为风机叶片半径,Cp为风能利用系数,v 为风速。
1.2 机侧变流器控制的数学模型
构网型 DDWT 机侧变流器控制目标为维持直流电压稳定,控制框图见图2。
引入中间变量 z1、z2、z3后,得到机侧变流器控制方程如下:
(3)
式中:Udcref为直流侧电压参考值,Udc为直流侧电压值,imdqref为输入至内环的电流参考值的 d、q 轴分量,umdqref为机侧换流器直流侧电压参考值的 d、q 轴分量,Kp1、Ki1为直流电压外环控制参数,Kp2、Ki2 、Kp3、Ki3为机侧电流内环控制参数。
图2机侧变流器控制框图
Fig.2Block diagram of Machine-side converter control
1.3 直流侧电容的数学模型
在风力发电系统中,直流侧电容器不仅连接发电机侧和电网侧,还起到平稳系统功率波动的作用,其模型可以表示为
(4)
式中:Cdc为直流侧电容值,pm 为定子侧发出的有功功率,pg为网侧输入的有功功率,i gdq为滤波后网侧电流的 d、q 轴分量,ugdq为滤波后网侧电压的 d、q 轴分量。
1.4 网侧变流器控制的数学模型
采集滤波器输出电压与电流经瞬时功率计算可以得到 VSG 输出功率如下:
(5)
控制框图如图3所示,根据 VSG 控制的功频和励磁控制器的框图可得 VSG 控制算法如下:
(6)
式中:w0 为电网额定频率,wg 为网侧同步角频率,θ 为网侧虚拟同步角度,Pe 为网侧输出有功功率,Pref为风机输出的功率指令参考值,Kp为有功调频系数,J 为虚拟惯量系数,Dp为虚拟阻尼系数,Q 为网侧输出无功功率,Qref为网侧无功参考值,Uref为机端电压幅值参考值,UN为空载电压,Dq为无功电压下垂系数,Um为网侧输出电压幅值,e 为电压的偏差量。
图3网侧变流器控制框图
Fig.3Block diagram of grid-side converter control
引入中间变量 z4 、z5 、z6、z7 后,得到网侧变流器控制方程如下:
(7)
式中:ugdqref为输入至电压外环的电压参考值的 d、q 轴分量,i gdqref为输入至电流内环的电流参考值的 d、q 轴分量,udqref为网侧换流器输出电压参考值的 d、q 轴分量,Kp4、Ki4、Kp7、Ki7为网侧电压外环控制参数,Kp5、Ki5、Kp6、Ki6为网侧电流内环控制参数,Lf为滤波电感值。在电压外环控制中引入虚拟阻抗,其中 Rv、Lv分别为虚拟电阻与虚拟电感。
1.5 滤波器与线路的数学模型
假设逆变器输出电压能够准确跟踪给定值,忽略电力电子器件开关延迟和损耗,结合基本的电路定律可得滤波器的数学表达式如下:
(8)
其中,电网电压 up的 d、q 轴分量可表示为:
式中:i odq为逆变器输出电流的 d、q 轴分量,uodq为逆变器输出电压的 d、q 轴分量,updq为汇流母线处电压的 d、q 轴分量,Rf为线路电阻,Lg为线路电感,Cf为滤波电容。
2 系统小信号建模与稳定性分析
联立构网型 DDWT 各模块的数学模型,对系统进行小信号建模,可得风电机组系统小信号模型为
(9)
其中:
风电机组的系数矩阵为
(10)
其中, ABCD 的表达式如下:
输入端的系数矩阵可表示为 BF,其表达式为
(11)
2.1 小信号模型准确性验证
针对本文建立的直驱风机状态空间小信号模型验证其精确有效性,本文在 MATLAB / Simulink 软件中搭建了构网型直驱风机并网系统仿真模型。同时依据式(9)建立了对应的线性小信号模型,设初始风速为 10 m / s,对两种模型施加相同的风速扰动(3 s时风速由 10 m / s 降低为 9 m / s)。对比两者在有功功率、无功功率,以及直流母线电压的动态响应曲线,可以观察其呈现出高度一致性,见图4。结果表明,所建立的 DDWT 状态空间小信号模型能够充分反映系统的动态行为,证明了该模型的准确性和有效性。
2.2 稳定性分析
利用 MATLAB 仿真平台,通过输入系统状态矩阵 AF,计算并求得 20 对特征根。再通过绘制主导特征值的根轨迹图,探讨主要控制参数对系统动态稳定性的影响。 DDWT 的特征根及主要状态参与因子见表1。
图4线性化的状态空间小信号模型对比验证
Fig.4Comparative validation of linearized state-space smallsignal models
表1DDWT 的参与因子
Tab.1Participating factors of DDWT
表1(续)
图5为主导模态 λ12、λ13的参与因子分析结果。
图5主导模态 λ12、λ13参与因子分析结果
Fig.5Participation factor analysis results of the dominant modes λ12 and λ13
由图5可以看出,在主导模态 λ12、λ13 中,与虚拟同步控制相关的状态变量具有较高的参与因子,说明该模态的动态特性由虚拟同步控制主导。图6为主导模态 λ14、λ15 的参与因子分析结果。由图6可以看出,在主导模态 λ14、λ15则由与直流母线电压控制相关的状态变量主导,反映出直流母线电压控制对其动态特性的关键影响。上述分析明确了虚拟同步控制与直流母线电压控制对不同模态的主导作用,为构网型控制优化设计提供了理论依据。
图6主导模态 λ14、λ15参与因子分析结果
Fig.6Participation factor analysis results of the dominant modes λ14 and λ15
为分析电网强度短路比( short circuit ratio,SCR)对构网型直驱风机并网稳定性的影响,保持其余参数恒定,计算短路比 SCR 从 3 变化到 5 时虚轴附近模态的根轨迹,结果见图7。由图7可以看出: 随着 SCR 增大,主导模态的实部向右移动,虚部变化较小,表明系统阻尼特性减弱,稳定性降低,而振荡频率基本保持稳定。
图7主导模态随 SCR 变化的根轨迹
Fig.7Root locus of dominant modes with the change of SCR
通过 MATLAB 建立时域仿真,对特征模态分析结果进一步验证。在 t = 3 s 时施加风速扰动,观察不同 SCR 条件下的有功功率动态响应,结果见图8。仿真显示,SCR 增大会显著增加功率波形的振荡幅度,振荡频率变化较小,与模态分析结果一致。整体分析表明,电网强度 SCR 的提升削弱了系统的阻尼性能,降低了系统的动态稳定性。
图8不同 SCR 下构网型直驱风机有功功率波形图
Fig.8Waveforms of grid-forming DDWT’s active power under different SCRs
通过特征模态分析和时域仿真可知,构网型直驱风机在低短路比条件下具备良好的弱电网适应性,能够实现稳定并网。然而,风速阶跃变化易引发低频模态振荡,表明阻尼特性减弱,稳定性下降。需进一步优化关键控制参数,以增强低频模态稳定性,提高系统运行可靠性。
3 关键参数对低频振荡的影响
构网型直驱风机虽具备良好的弱电网适应性,但在风速阶跃扰动下,其动态响应可能表现为低频振荡特性(振荡频率通常在 0.1~2.0 Hz),主要由阻尼特性不足引起。为深入研究风速扰动对系统动态稳定性的影响,本文分析关键参数对低频模态稳定性的作用机制,旨在通过优化参数设计提升低频模态的阻尼特性,从而确保系统在风速变化时实现稳定运行。
3.1 直流侧电容参数的影响
为研究直流侧电容参数 Cdc对于构网型直驱风机并网稳定性的影响,取短路比 SCR = 3,保持其他参数不变,计算 Cdc从 0.3~0.9 变化时虚轴附近模态的根轨迹,结果见图9。由图9可以看出:随着 Cdc增大,主导模态的实部向左移动,表明系统阻尼性能增强,稳定性提高; 同时虚部增大,对应振荡频率增加。
通过 MATLAB 仿真验证,在 t = 3 s 时施加风速扰动,观察不同 Cdc下直流母线电压的响应,结果见图10。由图10可以看出:随着 Cdc增大,电压的振荡幅度逐渐减小,衰减速度加快,动态恢复更平稳。结果表明,较大的直流侧电容有助于提升系统的阻尼性能和动态稳定性,与特征模态分析结果一致。
图9主导模态随 Cdc变化的根轨迹
Fig.9Root locus of dominant modes with the change of Cdc
图10不同 Cdc下构网型直驱风机直流母线电压波形图
Fig.10Waveforms of grid-forming DDWT’s DC bus voltage under different Cdc
3.2 虚拟同步控制参数的影响
为研究虚拟惯量 J 对于构网型直驱风机并网稳定性的影响,取短路比 SCR = 3,保持其他参数不变,计算 J 从 8~20 变化时虚轴附近模态的根轨迹,结果见图11。由图11可以看出:随着 J 的逐渐增大,主导模态的实部向右移动,表明系统阻尼性能减弱,整体稳定性下降; 虚部变化幅度较小,振荡频率略有降低。
通过 MATLAB 仿真对特征模态分析结果进行验证,保持其他参数不变,在 t = 3 s 时施加风速扰动,观察不同虚拟惯量下构网型直驱风机有功功率和直流母线电压的响应曲线。图12、13分别为不同虚拟惯量下构网型直驱风机有功功率与直流母线电压波形图。由图12、13可以看出:随着虚拟惯量 J 的增大,有功功率和直流母线电压的振荡幅度增大,衰减速度减慢,系统恢复时间延长,表明阻尼性能下降,系统稳定性减弱; 振荡频率略有降低,与特征模态分析结果一致。
图11主导模态随 J 变化的根轨迹
Fig.11Root locus of dominant modes with the change of J
图12不同 J 下构网型直驱风机有功功率波形图
Fig.12Waveforms of grid-forming DDWT’s active power under different J
图13不同 J 下构网型直驱风机直流母线电压波形图
Fig.13Waveforms of grid-forming DDWT’s DC bus voltage under different J
为研究虚拟阻尼 D 对于构网型直驱风机并网稳定性的影响,取短路比 SCR = 3,保持其他参数不变,计算 D 从 600~400 变化时虚轴附近模态的根轨迹,结果见图14。由图14可以看出:随着 D 的逐渐减小,主导模态的实部向右移动,虚部略有增加,表明系统阻尼特性减弱,振荡频率略微升高,整体稳定性下降。
图14主导模态随 D 变化的根轨迹
Fig.14Root locus of dominant modes with the change of D
通过 MATLAB 仿真对特征模态分析结果进行验证,保持其他参数不变,在 t = 3 s 时施加风速扰动,观察不同虚拟阻尼下构网型直驱风机有功功率和直流母线电压的响应曲线。图15、16分别为不同虚拟阻尼下构网型直驱风机有功功率与直流母线电压波形图。由图15、16可以看出,随着虚拟阻尼 D 的减小,有功功率和直流母线电压的波形振荡幅度增大,衰减速度明显减慢,系统恢复时间延长,表明阻尼性能下降,系统稳定性降低,与特征模态分析结果一致。
图15不同 D 下构网型直驱风机有功功率波形图
Fig.15Waveforms of grid-forming DDWT’s active power under different D
图16不同 D 下构网型直驱风机直流母线电压波形图
Fig.16Waveforms of grid-forming DDWT’s DC bus voltage under different D
3.3 虚拟阻抗参数的影响
为研究虚拟电阻 Rv对于构网型直驱风机并网稳定性的影响,取短路比 SCR = 3,保持其他参数不变,计算 Rv从 0.2~0.6 变化时虚轴附近模态的根轨迹,结果见图17。由图17可以看出:随着 Rv的增大,主导模态的实部向右移动,表明系统阻尼性能减弱,稳定性下降; 虚部变化不显著,振荡频率基本保持不变。
图17主导模态随 Rv变化的根轨迹
Fig.17Root locus of dominant modes with the change of Rv
通过 MATLAB 仿真对特征模态分析结果进行验证,保持其他参数不变,在 t = 3 s 时施加风速扰动,观察不同虚拟电阻下构网型直驱风机有功功率和直流母线电压的响应曲线。图18、19分别为不同虚拟电阻下构网型直驱风机有功功率与直流母线电压波形图。由图18、19可以看出,随着虚拟电阻 Rv 的增大,有功功率波形的振荡幅度略微增加,衰减速度变慢,直流母线电压的动态响应表现为振荡幅度增加和恢复时间延长的特征,表明系统阻尼性能下降,与特征模态分析结果一致。
图18不同 Rv下构网型直驱风机有功功率波形图
Fig.18Waveform of grid-forming DDWT’s active power under different Rv
图19不同 Rv下构网型直驱风机直流母线电压波形图
Fig.19Waveforms of grid-forming DDWT’s DC bus voltage under different Rv
为研究虚拟电抗 Lv对于构网型直驱风机并网稳定性的影响,取短路比 SCR = 3,保持其他参数不变,计算 Lv从 0. 004~0. 006 变化时虚轴附近模态的根轨迹,结果见图20。由图20可以看出:随着 Lv的增大,主导模态的实部向右移动,表明系统阻尼减弱,稳定性下降; 虚部变化不显著,振荡频率基本保持不变。
图20主导模态随 Lv变化的根轨迹
Fig.20Root locus of dominant modes with the change of Lv
通过 MATLAB 仿真对特征模态分析结果进行验证,保持其他参数不变,在 t = 3 s 时施加风速扰动,观察不同虚拟电抗下构网型直驱风机有功功率和直流母线电压的响应曲线。图21、22分别为不同虚拟电抗下构网型直驱风机有功功率与直流母线电压波形图。由图21、22可以看出,随着虚拟电抗 Lv 的增大,波形振荡幅度显著增大,衰减速度减慢,表明系统阻尼性能下降,与特征模态分析结果一致。
图21不同 Lv下构网型直驱风机有功功率波形图
Fig.21Waveforms of grid-forming DDWT’s active power under different Lv
图22不同 Lv下构网型直驱风机直流母线电压波形图
Fig.22Waveforms of grid-forming DDWT’s DC bus voltage under different Lv
综上分析,直流侧电容参数、虚拟同步控制参数,以及虚拟阻抗参数均显著影响构网型直驱风机的低频模态稳定性。通过合理调整这些控制参数,可有效改善低频模态的阻尼特性,从而提高系统稳定性,确保风机在复杂工况下的稳定运行。
4 仿真验证
为验证构网型 DDWT 并网系统的小信号模型的准确性,基于 RT-LAB 半实物仿真平台开展了动态仿真试验,具体仿真参数见表 2 [21]。
表2仿真参数
Tab.2Simulation parameters
4.1 基于 RT-LAB 的风速变化动态响应仿真
图23展示了风速发生阶跃变化时系统的动态响应行为。由图23( a)可以看出,当风速在 2 s 和 4 s 分别发生上升与下降的阶跃扰动时,系统频率随之波动,显示出明显的动态响应特性。与此对应,图23(b)、23(c)分别显示了有功功率和无功功率的低频振荡现象,这表明系统在风速扰动下的动态调节能力。
图23风速变化时系统响应仿真波形
Fig.23Simulated waveforms of system response when wind speed changes
图23(d)显示了直流母线电压的响应情况,其在风速扰动后出现短暂波动,但迅速恢复至稳态,体现了系统良好的调压特性。图23(e)反映了风机机械角速度的变化趋势,风速的阶跃扰动导致机械角速度发生调整,以释放或吸收动能。图23( f)为输出电流的动态响应波形,可见风速扰动引起的功率调整使输出电流随之变化,系统表现出类似同步发电机的惯量支撑特性。
4.2 基于 RT-LAB 的频率变化动态响应仿真
为验证构网型直驱风机的频率主动支撑能力,在不改变系统参数的情况下,设置电网频率于 2 s 时下降 0.05 Hz,4 s 时恢复。图24( a)显示了系统频率的扰动变化过程; 图24(b)展示了有功功率、无功功率,以及直流母线电压的动态响应; 图24(c)则反映了风机机械角速度的变化。
由图24(b)可以看出,频率下降时,系统通过调节有功功率增量进行频率支撑,无功功率也随之变化以稳定系统电压,同时直流母线电压的波动表明系统具备较好的抗扰动调压能力。图24( c)显示,风机控制器在检测到频率下降后,通过降低发电机转速释放动能,为系统提供短时能量支撑。整体结果验证了构网型直驱风机在频率扰动下的动态调节能力和稳定运行特性。
图24频率变化时系统响应仿真波形
Fig.24Simulated waveforms of system response when frequency changes
5 结论
本文针对构网型 DDWT 并网系统,基于源端动态特性建立小信号模型,经过稳定性分析与仿真验证,得出以下结论:
1)构网型直驱风机具备较强的弱电网适应性,即使短路比相对较低时,仍能够稳定并网运行,显示出其在弱电网环境下的良好性能。本文通过对风机并网动态特性的深入分析,验证了其在不同电网强度条件下的稳态运行能力,并量化了关键参数对系统稳定性的影响,为风电场在弱电网环境下的适应性评估提供了理论支撑。
2)研究采用改进的系统建模方法,将直流侧电容的动态作用纳入风电机组并网过程的稳定性分析中。结果表明,直流侧电容参数的变化对直流母线电压的动态响应具有显著影响,进而影响系统的低频振荡和功率输出稳定性。与传统忽略电容影响的模型相比,本文方法更精准地刻画了直流侧电压波动特性,提高了对风电机组动态行为的描述准确性。
3)网侧参数的调整会改变构网型直驱风机的低频模态特性,进而影响系统的阻尼性能,具体表现为直流母线电压的波动幅度和有功功率的动态响应特性。未来研究将进一步分析多风机并网的动态交互特性及参数优化策略,以提升弱电网环境下系统的整体稳定性和工程适用性。

