一种双馈风机等效惯量评估方法及系统与流程

1.本发明涉及评估领域，特别是涉及一种双馈风机等效惯量评估方法及系统。


背景技术：

2.在碳达峰、碳中和的时代背景下，加快构建以新能源为主体的新型电力系统是电力行业的必由之路。高比例可再生能源和高比例电力电子设备正在成为电力系统发展的重要趋势和关键特征。作为风电市场的主要机型，双馈风机灵活快速的控制特性使其转子动能与输出有功功率相对解耦，基本不具备惯量支撑能力。随着风电机组渗透率的提高，电力系统整体惯量水平势必会被削弱，严重恶化系统抗扰能力和稳定特性。为保证电网的安全稳定运行，虚拟惯量控制技术已成为提高电力系统惯量水平的重要手段，在双馈风机传统控制策略的基础上将电网频率偏差相关量引入变流器有功控制环节，通过吸收或释放转子动能来抵抗电网频率波动，进而向电网提供惯量支撑。
3.一般情况下，同步机组的等效惯性时间常数已经给定，且其数值基本不会发生变化。但对于双馈风机来说，其惯量响应特性与变流器类型、控制方式和机组运行工况息息相关。双馈风机出力的强随机性、间歇性和波动性使系统运行方式复杂多变，系统惯量在时间尺度将表现出较大的波动特性，所以含虚拟惯量控制的双馈风机接入电网后，系统惯量不再是无法改变的固有特性，其时空变化特征将愈加突出，扰动后的惯量响应过程将会更加复杂多变。惯量评估可为电网规划和调度运行提供科学参考，有助于电网运行调度人员掌握当前系统的等效惯量水平，辨识出惯量薄弱时段、节点或区域，进而有针对性地采取局部惯量补偿措施，为精细化惯量调控提供基础。
4.系统暂态稳定以及暂态失稳模式不仅与系统总惯量相关，与整个系统的惯量分布更为相关。然而，现有惯量评估方法大多聚焦于整个系统的等效惯量，关于双馈风机自身惯量及响应特性的研究，当前缺少相关理论方法。惯量评估最常见的计算方法为利用摇摆方程计算惯性时间常数，此方法需要根据频率曲线精确测定计算频率变化率，计算量大，复杂程度高。无论是频率变化率的滤波取值，还是扰动发生时间的确定都会给计算带来难度，计算出错时会影响惯量评估的准确性。除此之外，此方法依赖于电网扰动下的时间序列量测数据，无法实现常态化下系统惯量的实时连续监测。


技术实现要素：

5.本发明的目的是提供一种双馈风机等效惯量评估方法及系统，在计算等效惯性时间常数时过程简单、计算量小，能够提高惯量评估的准确性。
6.为实现上述目的，本发明提供了如下方案：
7.一种双馈风机等效惯量评估方法，包括：
8.获取稳态信息，所述稳态信息包括双馈风机的输出功率和电网频率；
9.根据所述双馈风机的输出功率计算初始转子角速度；
10.根据所述电网频率计算系统初始同步角速度；
11.基于所述系统初始同步角速度和所述初始转子角速度计算等效惯性时间常数，所述等效惯性时间常数用于评估双馈风机等效惯量。
12.可选的，所述基于所述系统初始同步角速度和所述初始转子角速度计算等效惯性时间常数，具体包括：
13.获取双馈风机参数，所述双馈风机参数包括：增益参数、频域微分算子、滤波器的时间常数、双馈风机固有惯性时间常数和双馈风机额定角速度；
14.根据所述双馈风机参数、所述系统初始同步角速度和所述初始转子角速度计算等效惯性时间常数。
15.可选的，所述根据所述双馈风机参数、所述系统初始同步角速度和所述初始转子角速度计算等效惯性时间常数，具体为：
16.根据求解公式计算等效惯性时间常数，其中，h
eq
表示等效惯性时间常数，k
dvic
表示虚拟惯量控制的微分增益，s表示频域微分算子，k
pvic
表示虚拟惯量控制的比例增益，t
f
表示滤波器的时间常数，k
p_pll
表示锁相环的比例增益，k
i_pll
表示锁相环的积分增益，ω
r0
表示初始转子角速度，ω
s0
表示系统初始同步角速度，h
d
表示双馈风机固有惯性时间常数，k
ps
表示速度控制器的比例增益，k
is
表示速度控制器的积分增益，ω
n
表示双馈风机额定角速度。
17.可选的，所述求解公式的确定方法为：
18.根据转速参考值与双馈风机输出的有功功率确定关系表达式；
19.以定子电压相角增量为输入以虚拟惯量控制提供的电磁转矩参考值增量为输出构建基于锁相环频率的比例微分虚拟惯量控制模型；
20.以双馈风机的电磁功率为输入以电磁转矩参考值增量为输出构建双馈风机的最大功率追踪模型；
21.以电磁转矩参考值为输入以双馈风机转子转速为输出构建双馈风机的发电机暂态模型；
22.根据第一动能变化量表达式和第二动能变化量表达式确定等效惯性时间常数表达式；所述第一动能变化量表达式为附加虚拟惯量控制前双馈风机惯量响应过程中的动能变化量表达式，所述第二动能变化量表达式为附加虚拟惯量控制后双馈风机惯量响应过程中的动能变化量；
23.根据三相同步锁相环动态方程、所述关系表达式、所述基于锁相环频率的比例微分虚拟惯量控制模型、所述双馈风机的最大功率追踪模型发电机暂态模型和所述等效惯性时间表达式确定所述求解公式。
24.可选的，所述根据第一动能变化量表达式和第二动能变化量表达式确定等效惯性时间常数表达式，具体包括：
25.根据所述第一动能变化量表达式和所述第二动能变化量表达式确定虚拟转动惯量表达式；
26.基于所述虚拟转动惯量表达式确定等效惯性时间常数表达式。
27.一种双馈风机等效惯量评估系统，包括：
28.获取模块，用于获取稳态信息，所述稳态信息包括双馈风机的输出功率和电网频率；
29.初始转子角速度确定模块，用于根据所述双馈风机的输出功率计算初始转子角速度；
30.系统初始同步角速度确定模块，用于根据所述电网频率计算系统初始同步角速度；
31.等效惯性时间常数确定模块，用于基于所述系统初始同步角速度和所述初始转子角速度计算等效惯性时间常数，所述等效惯性时间常数用于评估双馈风机等效惯量。
32.可选的，所述等效惯性时间常数确定模块，具体包括：
33.获取子模块，用于获取双馈风机参数，所述双馈风机参数包括：增益参数、频域微分算子、滤波器的时间常数、双馈风机固有惯性时间常数和双馈风机额定角速度；
34.等效惯性时间常数确定子模块，用于根据所述双馈风机参数、所述系统初始同步角速度和所述初始转子角速度计算等效惯性时间常数。
35.可选的，所述等效惯性时间常数确定子模块，具体包括：
36.求解单元，用于根据求解公式
37.计算等效惯性时间常数，其中，h
eq
表示等效惯性时间常数，k
dvic
表示虚拟惯量控制的微分增益，s表示频域微分算子，k
pvic
表示虚拟惯量控制的比例增益，t
f
表示滤波器的时间常数，k
p_pll
表示锁相环的比例增益，k
i_pll
表示锁相环的积分增益，ω
r0
表示初始转子角速度，ω
s0
表示系统初始同步角速度，h
d
表示双馈风机固有惯性时间常数，k
ps
表示速度控制器的比例增益，k
is
表示速度控制器的积分增益，ω
n
表示双馈风机额定角速度。
38.可选的，所述双馈风机等效惯量评估系统，还包括：
39.关系表达式确定模块，用于根据转速参考值与双馈风机输出的有功功率确定关系表达式；
40.第一模型构建模块，用于以定子电压相角增量为输入以虚拟惯量控制提供的电磁转矩参考值增量为输出构建基于锁相环频率的比例微分虚拟惯量控制模型；
41.第二模型构建模块，用于以双馈风机的电磁功率为输入以电磁转矩参考值增量为输出构建双馈风机的最大功率追踪模型；
42.第三模型构建模块，用于以电磁转矩参考值为输入以双馈风机转子转速为输出构建双馈风机的发电机暂态模型；
43.等效惯性时间常数表达式确定模块，用于根据第一动能变化量表达式和第二动能变化量表达式确定等效惯性时间常数表达式；所述第一动能变化量表达式为附加虚拟惯量控制前双馈风机惯量响应过程中的动能变化量表达式，所述第二动能变化量表达式为附加虚拟惯量控制后双馈风机惯量响应过程中的动能变化量；
44.求解公式确定模块，用于根据三相同步锁相环动态方程、所述关系表达式、所述基于锁相环频率的比例微分虚拟惯量控制模型、所述双馈风机的最大功率追踪模型发电机暂态模型和所述等效惯性时间表达式确定所述求解公式。
45.可选的，所述等效惯性时间常数表达式确定模块，具体包括：
46.虚拟转动惯量表达式确定子模块，用于根据所述第一动能变化量表达式和所述第二动能变化量表达式确定虚拟转动惯量表达式；
47.等效惯性时间常数表达式确定子模块，用于基于所述虚拟转动惯量表达式确定等效惯性时间常数表达式。
48.根据本发明提供的具体实施例，本发明公开了以下技术效果：本发明获取稳态信息，所述稳态信息包括双馈风机的输出功率和电网频率；根据所述双馈风机的输出功率计算初始转子角速度；根据所述电网频率计算系统初始同步角速度；基于所述系统初始同步角速度和所述初始转子角速度计算等效惯性时间常数，所述等效惯性时间常数用于评估双馈风机等效惯量，根据稳态信息就可以得到初始转子角速度和系统初始同步角速度，不需要计算频率变化率，相比于背景技术中利用摇摆方程计算惯性时间常数的计算方法本发明的计算方法计算量小没有复杂的计算过程，使得可以准确实时评估双馈风机等效惯量，可为电网规划和调度运行提供科学参考，有助于电网运行调度人员掌握当前系统的等效惯量水平，辨识出惯量薄弱时段、节点或区域，进而有针对性地采取局部惯量补偿措施，为精细化惯量调控提供基础。
附图说明
49.为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
50.图1为本发明实施例提供的的双馈风机等效惯量评估方法的流程图；
51.图2为本发明实施例提供的含虚拟惯量控制的双馈风机在转子转速时间尺度下的等效模型的结构图；
52.图3为本发明实施例提供的基于锁相环频率的比例微分虚拟惯量控制模型的结构图；
53.图4为本发明实施例提供的双馈风机的最大功率追踪模型的结构图；
54.图5为本发明实施例提供的双馈风机的发电机暂态模型的结构图；
55.图6为本发明实施例提供的改进的四机两区测试系统的结构图；
56.图7为采用本方法在双馈风机p
e
＝0.3时得到的等效惯性时间常数的评估值与真实值的比较结果图；
57.图8为采用本方法在双馈风机p
e
＝0.6时得到的等效惯性时间常数的评估值与真实值的比较结果图；
58.图9为本发明实施例提供的的双馈风机等效惯量评估系统的结构图。
具体实施方式
59.下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他
实施例，都属于本发明保护的范围。
60.为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
61.本实施例提供了一种双馈风机等效惯量评估方法，如图1所示，所述方法，包括：
62.步骤101：获取稳态信息。所述稳态信息包括双馈风机的输出功率和电网频率，是双馈风机稳态运行下的输出功率和电网频率。
63.步骤102：根据所述双馈风机的输出功率计算初始转子角速度ω
r0
。
64.步骤103：根据所述电网频率计算系统初始同步角速度ω
s0
。
65.步骤104：基于所述系统初始同步角速度和所述初始转子角速度计算等效惯性时间常数，所述等效惯性时间常数用于评估双馈风机等效惯量。
66.在实际应用中，ω
s0
可根据实时量测电网频率数据获取；不考虑风速的变化，稳态时风机的初始转子角速度ω
r0
等于转子转速参考值ω
r
*，ω
r
*的数值与风电场的实时风速有关，实际应用中可根据功率
‑
转速最优曲线即式(3)表达式获得。因此，利用双馈风机稳态运行下的输出功率和电网频率的量测信息分别获取ω
r0
和ω
s0
，进而实时评估双馈风机的等效惯量。
67.在实际应用中，步骤104，具体包括：
68.获取双馈风机参数，所述双馈风机参数包括：增益参数、频域微分算子、滤波器的时间常数、双馈风机固有惯性时间常数和双馈风机额定角速度；增益参数包括：虚拟惯量控制的微分增益、虚拟惯量控制的比例增益、锁相环的比例增益、锁相环的积分增益、速度控制器的比例增益和速度控制器的积分增益。
69.根据所述双馈风机参数、所述系统初始同步角速度和所述初始转子角速度计算等效惯性时间常数。
70.在实际应用中，双馈风机参数可根据铭牌数据或技术手册获得。
71.在实际应用中，所述根据所述双馈风机参数、所述系统初始同步角速度和所述初始转子角速度计算等效惯性时间常数，具体为：
72.根据求解公式(14)计算等效惯性时间常数。
73.在实际应用中，所述求解公式的确定方法为：
74.根据转速参考值与双馈风机输出的有功功率确定关系表达式。
75.以定子电压相角增量为输入以虚拟惯量控制提供的电磁转矩参考值增量为输出构建基于锁相环频率的比例微分虚拟惯量控制模型。
76.以双馈风机的电磁功率为输入以电磁转矩参考值增量为输出构建双馈风机的最大功率追踪模型。
77.以电磁转矩参考值为输入以双馈风机转子转速为输出构建双馈风机的发电机暂态模型。
78.根据第一动能变化量表达式和第二动能变化量表达式确定等效惯性时间常数表达式；所述第一动能变化量表达式为附加虚拟惯量控制前双馈风机惯量响应过程中的动能变化量表达式，所述第二动能变化量表达式为附加虚拟惯量控制后双馈风机惯量响应过程中的动能变化量。
79.根据三相同步锁相环动态方程、所述关系表达式、所述基于锁相环频率的比例微
分虚拟惯量控制模型、所述双馈风机的最大功率追踪模型发电机暂态模型和所述等效惯性时间表达式确定所述求解公式。
80.在实际应用中，所述根据第一动能变化量表达式和第二动能变化量表达式确定等效惯性时间常数表达式，具体包括：
81.根据所述第一动能变化量表达式和所述第二动能变化量表达式确定虚拟转动惯量表达式。
82.基于所述虚拟转动惯量表达式确定等效惯性时间常数表达式。
83.在实际应用中，基于锁相环频率的比例微分虚拟惯量控制模型、双馈风机的最大功率追踪模型和双馈风机的发电机暂态模型统称为含虚拟惯量控制的双馈风机在转子转速时间尺度下的等效模型，结构如图2所示。
84.其中，双馈风机借助于锁相环实现定子电压定向控制，本发明采用目前应用广泛的三相同步锁相环，三相同步锁相环动态方程为：
[0085][0086]
其中，x
pll
是跟踪定子d轴电压时的累积误差；是跟踪定子d轴电压时的累积误差的微分；u
ds
为定子电压u
s
的d轴分量；δ
pll
表示锁相环输出的相角；表示锁相环输出的相角的微分；ω
pll
为锁相环估算的dq坐标系旋转角速度；ω
s
为电网同步旋转角速度；k
i_pll
和k
p_pll
分别为锁相环的积分增益与锁相环的比例增益。
[0087]
比例微分虚拟惯量控制是在风电机组常规控制的基础上引入一个与系统频率偏差的比例和微分量相关的辅助功率。当系统频率波动时，双馈风机通过锁相环估算dq坐标系旋转角速度ω
pll
来测量系统频率，虚拟惯量控制提供的电磁转矩参考值增量δt
vic
的频域表达式(基于锁相环频率的比例微分虚拟惯量控制模型)为：
[0088][0089]
其中，k
dvic
和k
pvic
分别为虚拟惯量控制的微分增益和虚拟惯量控制的比例增益；t
f
为滤波器的时间常数；δω
pll
为锁相环估算的角速度变化，s表示频域微分算子，δδ
s
为定子电压相角增量，根据上述内容，可以得到基于锁相环频率的比例微分虚拟惯量控制模型如图3所示，其中输入为定子电压相角增量δδ
s
，输出为虚拟惯量控制提供的电磁转矩参考值增量δt
vic
。
[0090]
其中，最大功率追踪模型包含最大功率追踪(maximum power point tracking，mppt)控制和速度控制器。实际中由于叶轮处平均风速不可直接测量，因此通常根据机组的输出功率大小来计算相应的转速指令值大小。风机输出功率p
e
计算转速参考值ω
r
*的近似算法(关系表达式)为：
[0091][0092]
当实际检测到的双馈风机转子转速ω
r
与mppt获得的转速参考值ω
r
*间出现转速偏差时，以转速偏差作为输入量经速度控制器调节发电机电磁转矩参考值来帮助转子转速恢复到最佳运行状态。最大功率追踪模型下的电磁转矩参考值增量δt*的传递函数表达式(双馈风机的最大功率追踪模型)为：
[0093][0094]
其中，k
ps
和k
is
分别为速度控制器的比例增益和速度控制器的积分增益；ω
r
为风机实际转速。根据上述内容，可以得到双馈风机的最大功率追踪模型如图4所示，其中输入为双馈风机的电磁功率(输出功率)p
e
，输出为最大功率追踪模型下的电磁转矩参考值增量δt*。
[0095]
其中，虚拟惯量控制提供的电磁转矩参考值增量δt
vic
的频域表达式(基于锁相环频率的比例微分虚拟惯量控制模型)和最大功率追踪模型下的电磁转矩参考值增量δt*的传递函数表达式(双馈风机的最大功率追踪模型)，可以得到转子电流控制的电磁转矩参考值t
e*
，转子侧变流器控制通过定子电压定向矢量控制实现，定子磁链方程和电压方程为：
[0096][0097]
其中，ψ
ds
和ψ
qs
分别为d轴定子磁链和q轴定子磁链；u
ds
和u
qs
分别为定子电压u
s
在d轴、q轴上的分量；x
s
为定子电感；x
m
为定转子互感；i
ds
和i
dr
分别为定子绕组和转子绕组d轴电流；i
qs
和i
qr
分别为分别为定子绕组和转子绕组q轴电流，ψ
s
为定子磁链、u
s
为定子电压。
[0098]
转子电流控制的电磁转矩方程为：
[0099][0100]
其中，t
e
为电磁转矩；p为定子绕组极对数，右上角加“*”代表风机各参数的指令值。
[0101]
发电机转子运动方程为：
[0102]
t
m
‑
t
e
＝2h
d
ω
r
·
s
ꢀꢀꢀ
(7)
[0103]
其中，t
m
为双馈风机机械转矩；h
d
为双馈风机固有惯性时间常数；ω
r
为双馈风机转
子转速。
[0104]
双馈风机的发电机暂态模型(公式6和公式7)如图5所示，其中输入为转子电流控制的电磁转矩参考值t
e*
，输出为双馈风机转子转速ω
r
，(根据公式6的第二个公式得到i
qr*
，根据i
qr*
得到i
qr
，根据i
qr
和公式6的第一个公式得到t
e
，根据t
e
得到ω
r
)，k1＝
‑
2x
s
/(3px
m
)，gq(s)为电流内环控制环节，电流内环控制响应速度远快于转子转速尺度，认为转子电流快速跟踪指令值可认为gq(s)≈1。
[0105]
在实际应用中，第一动能变化量表达式为双馈风机惯量响应过程中的动能变化量可用转子转速变化量和固有转动惯量表示：
[0106][0107]
其中，p，s
n
，j
dfig
，ω
r0
，δω
r
分别为双馈风机极对数、额定容量、固有转动惯量、初始转子角速度和转子角速度增量。
[0108]
在实际应用中，第二动能变化量表达式为在附件加虚拟惯量控制后，双馈风机转速与系统角速度耦合，惯量响应过程中的动能变化量可用系统同步角速度变化量与虚拟转动惯量表示：
[0109][0110]
其中，j
eq
，ω
s0
，δω
s
分别为双馈风机虚拟转动惯量、系统初始同步角速度和系统同步角速度增量。
[0111]
在实际应用中，虚拟转动惯量表达式
[0112][0113]
在实际应用中，基于所述虚拟转动惯量表达式确定等效惯性时间常数表达式：
[0114]
基于虚拟转动惯量的表达式，结合等效惯性时间常数的物理意义，可确定等效惯性时间常数表达式为：
[0115][0116]
其中，ω
n
为双馈风机额定角速度。
[0117]
在实际应用中，根据三相同步锁相环动态方程、所述关系表达式、所述基于锁相环频率的比例微分虚拟惯量控制模型、所述双馈风机的最大功率追踪模型、发电机暂态模型和所述等效惯性时间表达式确定所述求解公式，具体包括：
[0118]
根据双馈风机在转子转速时间尺度下的等效模型、三相同步锁相环动态方程和所述关系表达式，可以得到dfig电磁转矩增量传递函数关系式为：
[0119][0120]
不考虑机械转矩的变化，由式(7)得转子运动方程增量频域表达式为：
[0121]
2h
d
sδω
r
＝
‑
δt
e
ꢀꢀꢀ
(13)
[0122]
结合等效惯性时间表达式，根据公式(12)和公式(13)，可推导求解公式为。
[0123][0124]
其中，h
eq
表示等效惯性时间常数，k
dvic
表示虚拟惯量控制的微分增益，s表示频域微分算子，k
pvic
表示虚拟惯量控制的比例增益，t
f
表示滤波器的时间常数，k
p_pll
表示锁相环的比例增益，k
i_pll
表示锁相环的积分增益，ω
r0
表示初始转子角速度，ω
s0
表示系统初始同步角速度，h
d
表示双馈风机固有惯性时间常数，k
ps
表示速度控制器的比例增益，k
is
表示速度控制器的积分增益，ω
n
表示双馈风机额定角速度。
[0125]
本实施例应用上述方法对改进的四机两区测试系统如图6所示进行处理，双馈风机替换节点1的同步机，通过一台变压器的单元接线方式连接到节点5。系统频率事故通过节点7处的负荷突增来模拟，测试时间为5s，以负荷突增时刻作为起始时间。图7和图8分别为双馈风机在不同工况下(p
e
＝0.3和p
e
＝0.6)等效惯性时间常数评估值与真实值的比较结果。
[0126]
由图7和图8可以看出，双馈风机等效惯性时间常数在不同工况下的评估值与真实值随时间变化的趋势基本一致，验证了本发明所提出的一种基于转子动能虚拟惯量控制的双馈风机等效惯量评估方法的准确性。
[0127]
本实施例提供了一种与上述方法对应的双馈风机等效惯量评估系统，如图9所示，所述系统包括：
[0128]
获取模块a1，用于获取稳态信息，所述稳态信息包括双馈风机的输出功率和电网频率。
[0129]
初始转子角速度确定模块a2，用于根据所述双馈风机的输出功率计算初始转子角速度。
[0130]
系统初始同步角速度确定模块a3，用于根据所述电网频率计算系统初始同步角速度。
[0131]
等效惯性时间常数确定模块a4，用于基于所述系统初始同步角速度和所述初始转子角速度计算等效惯性时间常数，所述等效惯性时间常数用于评估双馈风机等效惯量。
[0132]
作为一种可选的实施方式，所述等效惯性时间常数确定模块，具体包括：
[0133]
获取子模块，用于获取双馈风机参数，所述双馈风机参数包括：增益参数、频域微分算子、滤波器的时间常数、双馈风机固有惯性时间常数和双馈风机额定角速度。
[0134]
等效惯性时间常数确定子模块，用于根据所述双馈风机参数、所述系统初始同步角速度和所述初始转子角速度计算等效惯性时间常数。
[0135]
作为一种可选的实施方式，所述等效惯性时间常数确定子模块，具体包括：
[0136]
求解单元，用于根据求解公式
[0137]
计算等效惯性时间常数，其中，h
eq
表示等效惯性时间常数，k
dvic
表示虚拟惯量控制的微分增益，s表示频域微分算子，k
pvic
表示虚拟惯量控制的比例增益，t
f
表示滤波器的时间常数，k
p_pll
表示锁相环的比例增益，k
i_pll
表示锁相环的积分增益，ω
r0
表示初始转子角速度，ω
s0
表示系统初始同
步角速度，h
d
表示双馈风机固有惯性时间常数，k
ps
表示速度控制器的比例增益，k
is
表示速度控制器的积分增益，ω
n
表示双馈风机额定角速度。
[0138]
作为一种可选的实施方式，双馈风机等效惯量评估系统，还包括：
[0139]
关系表达式确定模块，用于根据转速参考值与双馈风机输出的有功功率确定关系表达式。
[0140]
第一模型构建模块，用于以定子电压相角增量为输入以虚拟惯量控制提供的电磁转矩参考值增量为输出构建基于锁相环频率的比例微分虚拟惯量控制模型。
[0141]
第二模型构建模块，用于以双馈风机的电磁功率为输入以电磁转矩参考值增量为输出构建双馈风机的最大功率追踪模型。
[0142]
第三模型构建模块，用于以电磁转矩参考值为输入以双馈风机转子转速为输出构建双馈风机的发电机暂态模型。
[0143]
等效惯性时间常数表达式确定模块，用于根据第一动能变化量表达式和第二动能变化量表达式确定等效惯性时间常数表达式。所述第一动能变化量表达式为附加虚拟惯量控制前双馈风机惯量响应过程中的动能变化量表达式，所述第二动能变化量表达式为附加虚拟惯量控制后双馈风机惯量响应过程中的动能变化量。
[0144]
求解公式确定模块，用于根据三相同步锁相环动态方程、所述关系表达式、所述基于锁相环频率的比例微分虚拟惯量控制模型、所述双馈风机的最大功率追踪模型发电机暂态模型和所述等效惯性时间表达式确定所述求解公式。
[0145]
作为一种可选的实施方式，所述等效惯性时间常数表达式确定模块，具体包括：
[0146]
虚拟转动惯量表达式确定子模块，用于根据所述第一动能变化量表达式和所述第二动能变化量表达式确定虚拟转动惯量表达式。
[0147]
等效惯性时间常数表达式确定子模块，用于基于所述虚拟转动惯量表达式确定等效惯性时间常数表达式。
[0148]
本发明有以下技术效果：
[0149]
(1)本发明计算速度快、复杂程度低，仅需双馈风机稳态运行时的量测数据即可较为精准地实时在线评估双馈风机等效惯量，具有很好的应用价值。
[0150]
(2)结合系统稳态运行下的量测数据，准确评估含虚拟惯量控制的双馈风机等效惯量具有重大意义，其可以为风电机组控制参数整定和虚拟惯量分配提供科学依据，有利于调度运行人员掌握系统惯量的时空分布情况，进而为实际系统安全预警和实时紧急控制策略提供参考。
[0151]
(3)现有方法需要扰动信息，而此方法需要的是稳态信息，并且现有方法需要根据频率曲线精确测定计算频率变化率，计算量大，复杂程度高，而本方法不需要计算频率变化率。另外，从研究目的上说，现有惯量评估方法大多聚焦于整个系统的等效惯量，而本文是针对风机自身惯量，在经历上述两步改进后，本方法结果依旧较为精准。
[0152]
(4)转子转速控制时间尺度动态问题，主要是由风力机相关控制、发电机转速控制、锁相环控制等控制以及该时间尺度储能元件(转子)的动态所产生的稳定问题，适合于虚拟惯量响应研究。因此，含虚拟惯量控制的双馈风机在转子转速时间尺度下的等效模型可以有针对性地解决双馈风机等效惯量评估问题。
[0153]
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他
实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
[0154]
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。