CN111092444A - 变速恒频风电机组惯性能量支撑控制方法及装置 - Google Patents

变速恒频风电机组惯性能量支撑控制方法及装置 Download PDF

Info

Publication number
CN111092444A
CN111092444A CN201911360494.3A CN201911360494A CN111092444A CN 111092444 A CN111092444 A CN 111092444A CN 201911360494 A CN201911360494 A CN 201911360494A CN 111092444 A CN111092444 A CN 111092444A
Authority
CN
China
Prior art keywords
wind turbine
variable
turbine generator
speed constant
frequency wind
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
CN201911360494.3A
Other languages
English (en)
Other versions
CN111092444B (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.)
China Electric Power Research Institute Co Ltd CEPRI
Zhejiang Windey Co Ltd
Original Assignee
China Electric Power Research Institute Co Ltd CEPRI
Zhejiang Windey 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 China Electric Power Research Institute Co Ltd CEPRI, Zhejiang Windey Co Ltd filed Critical China Electric Power Research Institute Co Ltd CEPRI
Priority to CN201911360494.3A priority Critical patent/CN111092444B/zh
Publication of CN111092444A publication Critical patent/CN111092444A/zh
Application granted granted Critical
Publication of CN111092444B publication Critical patent/CN111092444B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/24—Arrangements for preventing or reducing oscillations of power in networks
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
    • F03D—WIND MOTORS
    • F03D9/00—Adaptations of wind motors for special use; Combinations of wind motors with apparatus driven thereby; Wind motors specially adapted for installation in particular locations
    • F03D9/20—Wind motors characterised by the driven apparatus
    • F03D9/25—Wind motors characterised by the driven apparatus the apparatus being an electrical generator
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/70—Wind energy
    • Y02E10/72—Wind turbines with rotation axis in wind direction
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/70—Wind energy
    • Y02E10/76—Power conversion electric or electronic aspects

Abstract

本申请公开了一种变速恒频风电机组惯性能量控制方法、装置、设备及计算机可读存储介质，方法包括：根据机组的运行数据得到Cp‑λ散点分布阵列；得到机组在惯性能量支撑过程中的Cp‑λ曲线；根据Cp‑λ曲线及机组在惯性能量支撑过程中的气动功率得到ΔEaero；计算机组在惯性能量支撑过程中的ΔEe和ΔEloss，并根据转速约束得到ΔEkic，且根据ΔEaero‑ΔEe‑ΔEloss‑ΔEkic＝0得到机组在惯性能量支撑过程中的支撑时间；当机组的功率支撑时间达到支撑时间时，控制机组退出惯性能量支撑过程。本申请公开的上述技术方案，根据机组的实际运行情况来计算该机组在惯性能量支撑过程中的支撑时间，并根据支撑时间对该机组进行控制，以使得该机组能够发挥最大的支撑作用，从而达到较好的电网支撑效果，并防止机组意外停机。

Description

变速恒频风电机组惯性能量支撑控制方法及装置
技术领域
本申请涉及风电机组控制技术领域，更具体地说，涉及一种变速恒频风电机组惯性能量支撑控制方法、装置、设备及计算机可读存储介质。
背景技术
变速恒频风电机组的机械功率与电磁功率解耦、转子转速与电网频率解耦，因此，在电网频率波动时，无法表现出类似同步发电机的惯量特性，故当系统中风电接入比例增大、电网较弱时，势必会相对降低整个系统的总惯量和调频能力。为了更好地保证风电并网后系统的安全稳定运行，许多国家在其电网导则中已经对风电参与电力系统调频有了较为明确的要求和规定，其中，变速恒频风电机组提取惯性能量，参与系统频率的调节已成为一项重要而迫切的任务。
变速恒频风电机组以变转速载荷的方式运行，其旋转动能随时变化，因而在电网频率波动时能提供的支撑能量也在时刻变化。目前，变速恒频风电机组的惯性能量支撑控制一般是通过附加转矩控制方式实现，具体为变速恒频风电机组在感受到频率波动时，以固定功率增量进行一定时间的支撑，其在整个过程中支撑的能量比较大，因此，对于电网频率稳定可以起到更好的作用。但是，目前在采用该种控制方法进行控制时，无法获知变速恒频风电机组能够支撑的惯性能量的多少，且不知变速恒频风电机组具体应该在何时退出惯性能量支撑这一过程，若退出过早，则会使变速恒频风电机组无法发挥最好的支撑性能，无法获得最好的电网支撑效果，若退出过晚，则会导致变速恒频风电机组意外停机。
综上所述，如何对变速恒频风电机组惯性能量支撑过程进行控制，以达到较好的电网支撑效果，并防止变速恒频风电机组出现意外停机，是目前本领域技术人员亟待解决的技术问题。
发明内容
有鉴于此，本申请的目的是提供一种变速恒频风电机组惯性能量控制方法、装置、设备及计算机可读存储介质，用于对变速恒频风电机组惯性能量支撑过程进行控制，以达到较好的电网支撑效果，并防止变速恒频风电机组出现意外停机。
为了实现上述目的，本申请提供如下技术方案：
一种变速恒频风电机组惯性能量控制方法，包括：
采集变速恒频风电机组在并网运行过程中的运行数据，根据所述运行数据得到Cp-λ散点分布阵列；其中，Cp为风能利用系数，λ为叶尖速比；
对局部区域中的所述Cp-λ散点分布阵列进行拟合，得到所述变速恒频风电机组在惯性能量支撑过程中的Cp-λ曲线；
根据所述Cp-λ曲线及所述变速恒频风电机组在所述惯性能量支撑过程中的气动功率得到ΔEaero；其中，ΔEaero为叶轮吸收的风能；
计算所述变速恒频风电机组在所述惯性能量支撑过程中的ΔEe和ΔEloss，并根据转速约束得到ΔEkic，且根据ΔEaero-ΔEe-ΔEloss-ΔEkic＝0得到所述变速恒频风电机组在所述惯性能量支撑过程中的支撑时间；其中，ΔEe为所述变速恒频风电机组内除风能之外的能量转换成的输出电能，ΔEloss为机械和电气损耗，ΔEkic为惯性旋转动能；
当所述变速恒频风电机组在所述惯性能量支撑过程中的功率支撑时间达到所述支撑时间时，则控制所述变速恒频风电机组退出所述惯性能量支撑过程。
优选的，根据所述Cp-λ曲线及所述变速恒频风电机组在所述惯性能量支撑过程中的气动功率得到ΔEaero，包括：
根据所述Cp-λ曲线及
Figure BDA0002337053150000021
得到：
Figure BDA0002337053150000022
其中，Paero为所述变速恒频风电机组在所述惯性能量支撑过程中的气动功率，ρ为空气密度，A为叶轮扫掠面积，V为通过叶轮的等效风速，β为桨距角，R为叶轮半径，Tspt为所述支撑时间，ωT为叶轮转速，ωT＝Gωg，G为齿轮箱速比，ωg为发电机转速，
Figure BDA0002337053150000031
为所述变速恒频风电机组功率支撑终止状态时的叶轮转速，
Figure BDA0002337053150000032
为所述变速恒频风电机组功率支撑起始状态时的叶轮转速。
优选的，计算所述变速恒频风电机组在所述惯性能量支撑过程中的ΔEloss，包括：
根据
Figure BDA0002337053150000033
得到ΔEloss；
其中，
Figure BDA0002337053150000034
η为发电效率。
优选的，根据转速约束得到ΔEkic，包括：
根据
Figure BDA0002337053150000035
得到ΔEkic；
其中，J为转动惯量，ωg2为所述变速恒频风电机组功率支撑终止状态时的发电机转速，ωg1为所述变速恒频风电机组功率支撑起始状态时的发电机转速。
优选的，对局部区域中的所述Cp-λ散点分布阵列进行拟合，包括：
利用最小二乘拟合对局部区域中的所述Cp-λ散点分布阵列进行拟合。
优选的，利用最小二乘拟合对局部区域中的所述Cp-λ散点分布阵列进行拟合，包括：
利用最小二乘拟合且选取二阶多项式对局部区域中的所述Cp-λ散点分布阵列进行拟合。
一种变速恒频风电机组惯性能量控制装置，包括：
采集模块，用于采集变速恒频风电机组在并网运行过程中的运行数据，根据所述运行数据得到Cp-λ散点分布阵列；其中，Cp为风能利用系数，λ为叶尖速比；
拟合模块，用于对局部区域中的所述Cp-λ散点分布阵列进行拟合，得到所述变速恒频风电机组在惯性能量支撑过程中的Cp-λ曲线；
第一计算模块，用于根据所述Cp-λ曲线及所述变速恒频风电机组在所述惯性能量支撑过程中的气动功率得到ΔEaero；其中，ΔEaero为叶轮吸收的风能；
第二计算模块，用于计算所述变速恒频风电机组在所述惯性能量支撑过程中的ΔEe和ΔEloss，并根据转速约束得到ΔEkic，且根据ΔEaero-ΔEe-ΔEloss-ΔEkic＝0得到所述变速恒频风电机组在所述惯性能量支撑过程中的支撑时间；其中，ΔEe为所述变速恒频风电机组内除风能之外的能量转换成的输出电能，ΔEloss为机械和电气损耗，ΔEkic为惯性旋转动能；
控制模块，用于当所述变速恒频风电机组在所述惯性能量支撑过程中的功率支撑时间达到所述支撑时间时，则控制所述变速恒频风电机组退出所述惯性能量支撑过程。
优选的，所述第一计算模块包括：
第一计算单元，用于根据所述Cp-λ曲线及
Figure BDA0002337053150000041
得到：
Figure BDA0002337053150000042
其中，Paero为所述变速恒频风电机组在所述惯性能量支撑过程中的气动功率，ρ为空气密度，A为叶轮扫掠面积，V为通过叶轮的等效风速，β为桨距角，R为叶轮半径，Tspt为所述支撑时间，ωT为叶轮转速，ωT＝Gωg，G为齿轮箱速比，ωg为发电机转速，
Figure BDA0002337053150000043
为所述变速恒频风电机组功率支撑终止状态时的叶轮转速，
Figure BDA0002337053150000044
为所述变速恒频风电机组功率支撑起始状态时的叶轮转速。
一种变速恒频风电机组惯性能量控制设备，包括：
存储器，用于存储计算机程序；
处理器，用于执行所述计算机程序时实现如上述任一项所述的变速恒频风电机组惯性能量控制方法的步骤。
一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现如上述任一项所述的变速恒频风电机组惯性能量控制方法的步骤。
本申请提供了一种变速恒频风电机组惯性能量控制方法、装置、设备及计算机可读存储介质，其中，该方法包括：采集变速恒频风电机组在并网运行过程中的运行数据，根据运行数据得到Cp-λ散点分布阵列；其中，Cp为风能利用系数，λ为叶尖速比；对局部区域中的Cp-λ散点分布阵列进行拟合，得到变速恒频风电机组在惯性能量支撑过程中的Cp-λ曲线；根据Cp-λ曲线及变速恒频风电机组在惯性能量支撑过程中的气动功率得到ΔEaero；其中，ΔEaero为叶轮吸收的风能；计算变速恒频风电机组在惯性能量支撑过程中的ΔEe和ΔEloss，并根据转速约束得到ΔEkic，且根据ΔEaero-ΔEe-ΔEloss-ΔEkic＝0得到变速恒频风电机组在惯性能量支撑过程中的支撑时间；其中，ΔEe为变速恒频风电机组内除风能之外的能量转换成的输出电能，ΔEloss为机械和电气损耗，ΔEkic为惯性旋转动能；当变速恒频风电机组在惯性能量支撑过程中的功率支撑时间达到支撑时间时，则控制变速恒频风电机组退出惯性能量支撑过程。
本申请公开的上述技术方案，采集变速恒频风电机组在并网运行过程中的运行数据，并根据运行数据得到Cp-λ散点分布阵列，并通过Cp-λ散点分布阵列得到变速恒频风电机组在惯性能量支撑过程中的Cp-λ曲线，然后，利用Cp-λ曲线及变速恒频风电机组在惯性能量支撑过程中的气动功率得到叶轮吸收的风能ΔEaero，且计算变速恒频风电机组在所述惯性能量支撑过程中的ΔEe和ΔEloss，并根据转速约束得到ΔEkic，然后，根据ΔEaero-ΔEe-ΔEloss-ΔEkic＝0得到变速恒频风电机组在所述惯性能量支撑过程中的支撑时间，并根据支撑时间对变速恒频风电机组惯性能量支撑过程进行控制，即实现根据变速恒频风电机组的实际运行情况来计算得到该变速恒频风电机组在惯性能量支撑过程中的支撑时间，并根据所计算出的支撑时间对该变速恒频风电机组进行控制，以使得变速恒频风电机组能够发挥最大的支撑作用，从而达到较好的电网支撑效果，并防止变速恒频风电机组出现意外停机。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图。
图1为本申请实施例提供的一种变速恒频风电机组惯性能量控制方法的流程图；
图2为本申请实施例提供的变速恒频风电机组惯性能量支撑过程的工作过程示意图；
图3为本申请实施例提供的惯性能量支撑调节曲线的选择示意图；
图4为本申请实施例提供的惯性能量支撑测试效果-功率曲线图；
图5为本申请实施例提供的惯性能量支撑测试效果-转速曲线图；
图6为本申请实施例提供的Cp-λ曲线拟合结果示意图；
图7为本申请实施例提供的一种变速恒频风电机组惯性能量控制装置的结构示意图；
图8为本申请实施例提供的一种变速恒频风电机组惯性能量控制设备的结构示意图。
具体实施方式
下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
参见图1至图5，其中，图1示出了本申请实施例提供的一种变速恒频风电机组惯性能量控制方法的流程图，图2示出了本申请实施例提供的变速恒频风电机组惯性能量支撑过程的工作过程示意图，图3示出了本申请实施例提供的惯性能量支撑调节曲线的选择示意图，图4示出了本申请实施例提供的惯性能量支撑测试效果-功率曲线图，图5示出了本申请实施例提供的惯性能量支撑测试效果-转速曲线图。本申请实施例提供的一种变速恒频风电机组惯性能量控制方法，可以包括：
S11：采集变速恒频风电机组在并网运行过程中的运行数据，根据运行数据得到Cp-λ散点分布阵列。
其中，Cp为风能利用系数，λ为叶尖速比。
在变速恒频风电机组并网运行时，可以实时采集变速恒频风电机组的运行数据，其中，该运行数据包括但不限于变速恒频风电机组输出功率、空气密度、叶轮转速、风速的秒级数据等，然后，可以根据变速恒频风电机组的运行数据推导其的空气动力学特性关系，具体地，可以根据所采集到的运行数据及变速恒频风电机组本身的特性参数得到Cp-λ散点分布阵列，在Cp-λ散点分布阵列中每个Cp-λ散点具体为：
Figure BDA0002337053150000061
其中，Cp为风能利用系数，λ为叶尖速比，Ploss为损耗功率，Pe为发电功率，ρ为空气密度，A为叶轮扫掠面积，V为通过叶轮的等效风速，ωT为叶轮转速，R为叶轮半径。
S12：对局部区域中的Cp-λ散点分布阵列进行拟合，得到变速恒频风电机组在惯性能量支撑过程中的Cp-λ曲线。
在得到Cp-λ散点分布阵列之后，可以对局部区域(具体为与变速恒频风电机组做惯性能量支撑相对应的区域)中的Cp-λ散点分布阵列进行拟合，并通过拟合得到变速恒频风电机组在惯性能量支撑过程中的Cp-λ曲线。
S13：根据Cp-λ曲线及变速恒频风电机组在惯性能量支撑过程中的气动功率得到ΔEaero。
其中，ΔEaero为叶轮吸收的风能。
然后，根据拟合得到的Cp-λ曲线及变速恒频风电机组在惯性能量支撑过程中的气动功率
Figure BDA0002337053150000071
得到
Figure BDA0002337053150000072
其中，ΔEaero为叶轮吸收的风能，CP(β,λ)为风能利用系数，β为桨距角，Tspt为支撑时间。
上述过程即是将对变速恒频风电机组的气动性能产生影响的地形、气流、温度、湿度、老化程度等多方面的因素考虑在内，以得到变速恒频风电机组的运行数据，并将变速恒频风电机组的运行数据参与到叶轮吸收的风能计算中，从而使其参与到变速恒频风电机组惯性能量支撑过程中，以最终提高支撑时间计算的准确性。
S14：计算变速恒频风电机组在惯性能量支撑过程中的ΔEe和ΔEloss，并根据转速约束得到ΔEkic，且根据ΔEaero-ΔEe-ΔEloss-ΔEkic＝0得到变速恒频风电机组在惯性能量支撑过程中的支撑时间。
其中，ΔEe为变速恒频风电机组内除风能之外的能量转换成的输出电能，ΔEloss为机械和电气损耗，ΔEkic为惯性旋转动能。
同时，计算变速恒频风电机组在惯性能量支撑过程中的ΔEe和ΔEloss，并根据转速约束得到ΔEkic。在整个惯性能量支撑过程中，ΔEaero、ΔEe、ΔEloss和ΔEkic满足动态能量平衡关系ΔEaero-ΔEe-ΔEloss-ΔEkic＝0，而且由于
Figure BDA0002337053150000073
则可以根据动态能量平衡关系得到变速恒频风电机组在惯性能量支撑过程中的支撑时间：
Figure BDA0002337053150000074
其中，ΔEe为变速恒频风电机组内除风能之外的能量转换成的输出电能，Pspt为支撑功率，ΔEloss为机械和电气损耗，ΔEkic为惯性旋转动能。
需要说明的是，上述所提及的转速约束条件是由实际采集到的叶轮转速的最大值而确定的，以便于计算出风速恒频风电机组在每一运行时刻可进行功率支撑的最大惯性旋转动能(即ΔEkic)，并使该惯性旋转动能参与到惯性能量支撑过程中，以便于确定变速恒频风电机组实现最大化能量支撑效果的支撑时间，从而使得变速恒频风电机组可以发挥最大的支撑作用。
S15：当变速恒频风电机组在惯性能量支撑过程中的功率支撑时间达到支撑时间时，则控制变速恒频风电机组退出惯性能量支撑过程。
结合图2可知，当系统频率发生突降时，变速恒频风电机组会超额发出固定电功率支持电网，当变速恒频风电机组在惯性能量支撑过程中的功率支撑时间达到支撑时间Tspt时，则控制变速恒频风电机组退出惯性能量支撑过程，以使得变速恒频风电机组可以在正确的时间退出惯性能量支撑这一过程，从而使得变速恒频风电机组可以发挥最大的支撑作用，达到较好的电网支撑效果，并防止变速恒频风电机组出现意外停机，从而提高变速恒频风电机组运行的可靠性。
具体地，结合图2可知，在经过Tspt(即图2中从点2运动至点3)后，转速下降，继而输出功率下降，使得变速恒频风电机组重新加速，回调到原有状态。
本申请公开的上述技术方案，采集变速恒频风电机组在并网运行过程中的运行数据，并根据运行数据得到Cp-λ散点分布阵列，并通过Cp-λ散点分布阵列得到变速恒频风电机组在惯性能量支撑过程中的Cp-λ曲线，然后，利用Cp-λ曲线及变速恒频风电机组在惯性能量支撑过程中的气动功率得到叶轮吸收的风能ΔEaero，且计算变速恒频风电机组在所述惯性能量支撑过程中的ΔEe和ΔEloss，并根据转速约束得到ΔEkic，然后，根据ΔEaero-ΔEe-ΔEloss-ΔEkic＝0得到变速恒频风电机组在所述惯性能量支撑过程中的支撑时间，并根据支撑时间对变速恒频风电机组惯性能量支撑过程进行控制，即实现根据变速恒频风电机组的实际运行情况来计算得到该变速恒频风电机组在惯性能量支撑过程中的支撑时间，并根据所计算出的支撑时间对该变速恒频风电机组进行控制，以使得变速恒频风电机组能够发挥最大的支撑作用，从而达到较好的电网支撑效果，并防止变速恒频风电机组出现意外停机。
本申请实施例提供的一种变速恒频风电机组惯性能量控制方法，根据Cp-λ曲线及变速恒频风电机组在惯性能量支撑过程中的气动功率得到ΔEaero，可以包括：
根据Cp-λ曲线及
Figure BDA0002337053150000091
得到：
Figure BDA0002337053150000092
其中，Paero为变速恒频风电机组在惯性能量支撑过程中的气动功率，ρ为空气密度，A为叶轮扫掠面积，V为通过叶轮的等效风速，β为桨距角，R为叶轮半径，Tspt为支撑时间，ωT为叶轮转速，ωT＝Gωg，G为齿轮箱速比，ωg为发电机转速，
Figure BDA0002337053150000093
为变速恒频风电机组功率支撑终止状态时的叶轮转速，
Figure BDA0002337053150000094
为变速恒频风电机组功率支撑起始状态时的叶轮转速。
根据Cp-λ曲线及变速恒频风电机组在惯性能量支撑过程中的气动功率得到ΔEaero的具体过程为：
在惯性能量支撑过程中，可以认为风速保持不变，则：
Figure BDA0002337053150000095
同时，由于变速恒频风电机组在变速运行段不存在β变化对CP(β,λ)的影响，且由于CP(β,λ)连续变化，则可以近似对Cp-λ曲线中的CP进行二阶泰勒级数展开：
Figure BDA0002337053150000096
进而
Figure BDA0002337053150000097
Figure BDA0002337053150000098
将使(3)和式(4)带入式(1)，同时近似认为在进行惯性能量支撑过程中转速线性下降，那么最终可得
Figure BDA0002337053150000099
其中，ωT＝Gωg，G为齿轮箱速比，ωg为发电机转速，
Figure BDA00023370531500000910
为变速恒频风电机组功率支撑终止状态时的叶轮转速，
Figure BDA0002337053150000101
为变速恒频风电机组功率支撑起始状态时的叶轮转速。
本申请实施例提供的一种变速恒频风电机组惯性能量控制方法，计算变速恒频风电机组在惯性能量支撑过程中的ΔEloss，可以包括：
根据
Figure BDA0002337053150000102
得到ΔEloss；
其中，
Figure BDA0002337053150000103
η为发电效率。
可以根据
Figure BDA0002337053150000104
得到ΔEloss，其中，η为发电效率，以提高ΔEloss计算的便利性。
本申请实施例提供的一种变速恒频风电机组惯性能量控制方法，根据转速约束得到ΔEkic，可以包括：
根据
Figure BDA0002337053150000105
得到ΔEkic；
其中，J为转动惯量，ωg2为变速恒频风电机组功率支撑终止状态时的发电机转速，ωg1为变速恒频风电机组功率支撑起始状态时的发电机转速。
在根据转速约束得到ΔEkic时，具体可以利用
Figure BDA0002337053150000106
得到ΔEkic，其中，J为转动惯量，ωg2为变速恒频风电机组功率支撑终止状态时的发电机转速，ωg1为变速恒频风电机组功率支撑起始状态时的发电机转速，以便于根据变速恒频风电机组在惯性能量支撑起始点的转速得到变速恒频风电机组的惯性旋转动能。
本申请实施例提供的一种变速恒频风电机组惯性能量控制方法，对局部区域中的Cp-λ散点分布阵列进行拟合，可以包括：
利用最小二乘拟合对局部区域中的Cp-λ散点分布阵列进行拟合。
在对局部区域中的Cp-λ散点分布阵列进行拟合时，可以利用最小二乘拟合对局部区域中的Cp-λ散点分布阵列进行拟合，其具有简单、便于计算等特点。
具体地，可以设定Cp-λ散点分布阵列在上述局部区域中的关系为f(x)＝a0+a1x+a2x2+…+aixi+Λ，为了使拟合出的Cp-λ曲线能够尽量反映所给数据的变化取值，则要求所有数据点上的残差|δi|＝|f(xi)-yi|都较小，其中，f(xi)为xi对应的拟合值，yi为xi对应的实际值。为了达到该目标，可以令残差的平方和最小，即：
Figure BDA0002337053150000111
通过上述关系得到参数估计值a0、a1、a2…ai，最终获取到Cp-λ曲线及其对应的多项式。
本申请实施例提供的一种变速恒频风电机组惯性能量控制方法，利用最小二乘拟合对局部区域中的Cp-λ散点分布阵列进行拟合，可以包括：
利用最小二乘拟合且选取二阶多项式对局部区域中的Cp-λ散点分布阵列进行拟合。
在利用最小二乘拟合对局部区域中的Cp-λ散点分布阵列进行拟合时，为了简化计算过程，则可以取二阶多项式进行拟合，对应结果可以参见图6，其示出了本申请实施例提供的Cp-λ曲线拟合结果示意图，当然，也可以采用三阶多项式等对Cp-λ散点分布阵列进行拟合，本申请对多项式项数不做任何限定。
本申请实施例还提供了一种变速恒频风电机组惯性能量控制装置，参见图7，其示出了本申请实施例提供的一种变速恒频风电机组惯性能量控制装置的结构示意图，可以包括：
采集模块71，用于采集变速恒频风电机组在并网运行过程中的运行数据，根据运行数据得到Cp-λ散点分布阵列；其中，Cp为风能利用系数，λ为叶尖速比；
拟合模块72，用于对局部区域中的Cp-λ散点分布阵列进行拟合，得到变速恒频风电机组在惯性能量支撑过程中的Cp-λ曲线；
第一计算模块73，用于根据Cp-λ曲线及变速恒频风电机组在惯性能量支撑过程中的气动功率得到ΔEaero；其中，ΔEaero为叶轮吸收的风能；
第二计算模块74，用于计算变速恒频风电机组在惯性能量支撑过程中的ΔEe和ΔEloss，并根据转速约束得到ΔEkic，且根据ΔEaero-ΔEe-ΔEloss-ΔEkic＝0得到变速恒频风电机组在惯性能量支撑过程中的支撑时间；其中，ΔEe为变速恒频风电机组内除风能之外的能量转换成的输出电能，ΔEloss为机械和电气损耗，ΔEkic为惯性旋转动能；
控制模块75，用于当变速恒频风电机组在惯性能量支撑过程中的功率支撑时间达到支撑时间时，则控制变速恒频风电机组退出惯性能量支撑过程。
本申请实施例提供的一种变速恒频风电机组惯性能量控制装置，第一计算模块73可以包括：
第一计算单元，用于根据Cp-λ曲线及
Figure BDA0002337053150000121
得到：
Figure BDA0002337053150000122
其中，Paero为变速恒频风电机组在惯性能量支撑过程中的气动功率，ρ为空气密度，A为叶轮扫掠面积，V为通过叶轮的等效风速，β为桨距角，R为叶轮半径，ωT为叶轮转速，Tspt为支撑时间，ωT＝Gωg，G为齿轮箱速比，ωg为发电机转速，
Figure BDA0002337053150000123
为变速恒频风电机组功率支撑终止状态时的叶轮转速，
Figure BDA0002337053150000124
为变速恒频风电机组功率支撑起始状态时的叶轮转速。
本申请实施例提供的一种变速恒频风电机组惯性能量控制装置，第二计算模块74可以包括：
第二计算单元，用于根据
Figure BDA0002337053150000125
得到ΔEloss；
其中，
Figure BDA0002337053150000126
η为发电效率。
本申请实施例提供的一种变速恒频风电机组惯性能量控制装置，第二计算模块74可以包括：
第三计算单元，用于根据
Figure BDA0002337053150000127
得到ΔEkic；
其中，J为转动惯量，ωg2为变速恒频风电机组功率支撑终止状态时的发电机转速，ωg1为变速恒频风电机组功率支撑起始状态时的发电机转速。
本申请实施例提供的一种变速恒频风电机组惯性能量控制装置，拟合模块72可以包括：
拟合单元，用于利用最小二乘拟合对Cp-λ散点分布阵列进行拟合。
本申请实施例提供的一种变速恒频风电机组惯性能量控制装置，拟合单元可以包括：
拟合子单元，用于利用最小二乘拟合且选取二阶多项式对Cp-λ散点分布阵列进行拟合。
本申请实施例还提供了一种变速恒频风电机组惯性能量控制设备，参见图8，其示出了本申请实施例提供的一种变速恒频风电机组惯性能量控制设备的结构示意图，可以包括：
存储器81，用于存储计算机程序；
处理器82，用于执行存储器81存储的计算机程序时可实现如下步骤：
采集变速恒频风电机组在并网运行过程中的运行数据，根据运行数据得到Cp-λ散点分布阵列；其中，Cp为风能利用系数，λ为叶尖速比；对局部区域中的Cp-λ散点分布阵列进行拟合，得到变速恒频风电机组在惯性能量支撑过程中的Cp-λ曲线；根据Cp-λ曲线及变速恒频风电机组在惯性能量支撑过程中的气动功率得到ΔEaero；其中，ΔEaero为叶轮吸收的风能；计算变速恒频风电机组在惯性能量支撑过程中的ΔEe和ΔEloss，并根据转速约束得到ΔEkic，且根据ΔEaero-ΔEe-ΔEloss-ΔEkic＝0得到变速恒频风电机组在惯性能量支撑过程中的支撑时间；其中，ΔEe为变速恒频风电机组内除风能之外的能量转换成的输出电能，ΔEloss为机械和电气损耗，ΔEkic为惯性旋转动能；当变速恒频风电机组在惯性能量支撑过程中的功率支撑时间达到支撑时间时，则控制变速恒频风电机组退出惯性能量支撑过程。
本申请实施例还提供了一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序，计算机程序被处理器执行时可实现如下步骤：
采集变速恒频风电机组在并网运行过程中的运行数据，根据运行数据得到Cp-λ散点分布阵列；其中，Cp为风能利用系数，λ为叶尖速比；对局部区域中的Cp-λ散点分布阵列进行拟合，得到变速恒频风电机组在惯性能量支撑过程中的Cp-λ曲线；根据Cp-λ曲线及变速恒频风电机组在惯性能量支撑过程中的气动功率得到ΔEaero；其中，ΔEaero为叶轮吸收的风能；计算变速恒频风电机组在惯性能量支撑过程中的ΔEe和ΔEloss，并根据转速约束得到ΔEkic，且根据ΔEaero-ΔEe-ΔEloss-ΔEkic＝0得到变速恒频风电机组在惯性能量支撑过程中的支撑时间；其中，ΔEe为变速恒频风电机组内除风能之外的能量转换成的输出电能，ΔEloss为机械和电气损耗，ΔEkic为惯性旋转动能；当变速恒频风电机组在惯性能量支撑过程中的功率支撑时间达到支撑时间时，则控制变速恒频风电机组退出惯性能量支撑过程。
该计算机可读存储介质包括：U盘、移动硬盘、只读存储器(Read-Only Memory，ROM)、随机存取存储器(Random Access Memory，RAM)、磁碟或者光盘等各种可以存储程序代码的介质。
本申请实施例提供的一种变速恒频风电机组惯性能量控制装置、设备及计算机可读存储介质中相关部分的说明可以参见本申请实施例提供的一种变速恒频风电机组惯性能量控制方法中对应部分的详细说明，在此不再赘述。
需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。另外，本申请实施例提供的上述技术方案中与现有技术中对应技术方案实现原理一致的部分并未详细说明，以免过多赘述。
对所公开的实施例的上述说明，使本领域技术人员能够实现或使用本申请。对这些实施例的多种修改对本领域技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本申请的精神或范围的情况下，在其它实施例中实现。因此，本申请将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (10)

1.一种变速恒频风电机组惯性能量控制方法，其特征在于，包括：
采集变速恒频风电机组在并网运行过程中的运行数据，根据所述运行数据得到Cp-λ散点分布阵列；其中，Cp为风能利用系数，λ为叶尖速比；
对局部区域中的所述Cp-λ散点分布阵列进行拟合，得到所述变速恒频风电机组在惯性能量支撑过程中的Cp-λ曲线；
根据所述Cp-λ曲线及所述变速恒频风电机组在所述惯性能量支撑过程中的气动功率得到ΔEaero；其中，ΔEaero为叶轮吸收的风能；
计算所述变速恒频风电机组在所述惯性能量支撑过程中的ΔEe和ΔEloss，并根据转速约束得到ΔEkic，且根据ΔEaero-ΔEe-ΔEloss-ΔEkic＝0得到所述变速恒频风电机组在所述惯性能量支撑过程中的支撑时间；其中，ΔEe为所述变速恒频风电机组内除风能之外的能量转换成的输出电能，ΔEloss为机械和电气损耗，ΔEkic为惯性旋转动能；
当所述变速恒频风电机组在所述惯性能量支撑过程中的功率支撑时间达到所述支撑时间时，则控制所述变速恒频风电机组退出所述惯性能量支撑过程。
2.根据权利要求1所述的变速恒频风电机组惯性能量控制方法，其特征在于，根据所述Cp-λ曲线及所述变速恒频风电机组在所述惯性能量支撑过程中的气动功率得到ΔEaero，包括：
根据所述Cp-λ曲线及
Figure FDA0002337053140000011
得到：
Figure FDA0002337053140000012
其中，Paero为所述变速恒频风电机组在所述惯性能量支撑过程中的气动功率，ρ为空气密度，A为叶轮扫掠面积，V为通过叶轮的等效风速，β为桨距角，R为叶轮半径，Tspt为所述支撑时间，ωT为叶轮转速，ωT＝Gωg，G为齿轮箱速比，ωg为发电机转速，
Figure FDA0002337053140000014
为所述变速恒频风电机组功率支撑终止状态时的叶轮转速，
Figure FDA0002337053140000015
为所述变速恒频风电机组功率支撑起始状态时的叶轮转速。
3.根据权利要求2所述的变速恒频风电机组惯性能量控制方法，其特征在于，计算所述变速恒频风电机组在所述惯性能量支撑过程中的ΔEloss，包括：
根据
Figure FDA0002337053140000013
得到ΔEloss；
其中，
Figure FDA0002337053140000021
η为发电效率。
4.根据权利要求3所述的变速恒频风电机组惯性能量控制方法，其特征在于，根据转速约束得到ΔEkic，包括：
根据
Figure FDA0002337053140000022
得到ΔEkic；
其中，J为转动惯量，ωg2为所述变速恒频风电机组功率支撑终止状态时的发电机转速，ωg1为所述变速恒频风电机组功率支撑起始状态时的发电机转速。
5.根据权利要求1至4任一项所述的变速恒频风电机组惯性能量控制方法，其特征在于，对局部区域中的所述Cp-λ散点分布阵列进行拟合，包括：
利用最小二乘拟合对局部区域中的所述Cp-λ散点分布阵列进行拟合。
6.根据权利要求5所述的变速恒频风电机组惯性能量控制方法，其特征在于，利用最小二乘拟合对局部区域中的所述Cp-λ散点分布阵列进行拟合，包括：
利用最小二乘拟合且选取二阶多项式对局部区域中的所述Cp-λ散点分布阵列进行拟合。
7.一种变速恒频风电机组惯性能量控制装置，其特征在于，包括：
采集模块，用于采集变速恒频风电机组在并网运行过程中的运行数据，根据所述运行数据得到Cp-λ散点分布阵列；其中，Cp为风能利用系数，λ为叶尖速比；
拟合模块，用于对局部区域中的所述Cp-λ散点分布阵列进行拟合，得到所述变速恒频风电机组在惯性能量支撑过程中的Cp-λ曲线；
第一计算模块，用于根据所述Cp-λ曲线及所述变速恒频风电机组在所述惯性能量支撑过程中的气动功率得到ΔEaero；其中，ΔEaero为叶轮吸收的风能；
第二计算模块，用于计算所述变速恒频风电机组在所述惯性能量支撑过程中的ΔEe和ΔEloss，并根据转速约束得到ΔEkic，且根据ΔEaero-ΔEe-ΔEloss-ΔEkic＝0得到所述变速恒频风电机组在所述惯性能量支撑过程中的支撑时间；其中ΔEe为所述变速恒频风电机组内除风能之外的能量转换成的输出电能，ΔEloss为机械和电气损耗，ΔEkic为惯性旋转动能；
控制模块，用于当所述变速恒频风电机组在所述惯性能量支撑过程中的功率支撑时间达到所述支撑时间时，则控制所述变速恒频风电机组退出所述惯性能量支撑过程。
8.根据权利要求7所述的变速恒频风电机组惯性能量控制装置，其特征在于，所述第一计算模块包括：
第一计算单元，用于根据所述Cp-λ曲线及
Figure FDA0002337053140000031
得到：
Figure FDA0002337053140000032
其中，Paero为所述变速恒频风电机组在所述惯性能量支撑过程中的气动功率，ρ为空气密度，A为叶轮扫掠面积，V为通过叶轮的等效风速，β为桨距角，R为叶轮半径，Tspt为所述支撑时间，ωT为叶轮转速，ωT＝Gωg，G为齿轮箱速比，ωg为发电机转速，
Figure FDA0002337053140000033
为所述变速恒频风电机组功率支撑终止状态时的叶轮转速，
Figure FDA0002337053140000034
为所述变速恒频风电机组功率支撑起始状态时的叶轮转速。
9.一种变速恒频风电机组惯性能量控制设备，其特征在于，包括：
存储器，用于存储计算机程序；
处理器，用于执行所述计算机程序时实现如权利要求1至6任一项所述的变速恒频风电机组惯性能量控制方法的步骤。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现如权利要求1至6任一项所述的变速恒频风电机组惯性能量控制方法的步骤。
CN201911360494.3A 2019-12-25 2019-12-25 变速恒频风电机组惯性能量支撑控制方法及装置 Active CN111092444B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911360494.3A CN111092444B (zh) 2019-12-25 2019-12-25 变速恒频风电机组惯性能量支撑控制方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911360494.3A CN111092444B (zh) 2019-12-25 2019-12-25 变速恒频风电机组惯性能量支撑控制方法及装置

Publications (2)

Publication Number Publication Date
CN111092444A true CN111092444A (zh) 2020-05-01
CN111092444B CN111092444B (zh) 2021-11-09

Family

ID=70397453

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911360494.3A Active CN111092444B (zh) 2019-12-25 2019-12-25 变速恒频风电机组惯性能量支撑控制方法及装置

Country Status (1)

Country Link
CN (1) CN111092444B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112564132A (zh) * 2020-12-15 2021-03-26 深圳供电局有限公司 一种风电一次调频潜力不确定性建模方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110153099A1 (en) * 2008-06-30 2011-06-23 Vestas Wind Systems A/S Method and system for controlling a wind power plant comprising a number of wind turbine generators
CN102403732A (zh) * 2011-11-08 2012-04-04 华北电力大学（保定） 双馈风力发电机组的虚拟惯性控制方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110153099A1 (en) * 2008-06-30 2011-06-23 Vestas Wind Systems A/S Method and system for controlling a wind power plant comprising a number of wind turbine generators
CN102403732A (zh) * 2011-11-08 2012-04-04 华北电力大学（保定） 双馈风力发电机组的虚拟惯性控制方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
许国东等: "《风电机组惯性能量支撑过程量化分析》", 《浙江电力》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112564132A (zh) * 2020-12-15 2021-03-26 深圳供电局有限公司 一种风电一次调频潜力不确定性建模方法

Also Published As

Publication number Publication date
CN111092444B (zh) 2021-11-09

Similar Documents

Publication Publication Date Title
Zhang et al. Pitch angle control for variable speed wind turbines
CN104917201B (zh) 模拟惯性与超速相结合的双馈风机有功频率控制器及方法
CN105134487B (zh) 一种考虑湍流频率因素的风力机最大功率点跟踪控制方法
CN110098622A (zh) 风电机组的一次调频方法和设备
CN109560573B (zh) 一种变速风电机组频率控制器参数的优化方法及装置
CN102611132B (zh) 双馈变速风电机组的附加频率控制器参数整定方法
CN107742903B (zh) 风电机组参与一次调频后退出调频时的转速恢复方法
CN110417032A (zh) 一种双馈风机参与系统调频的多目标优化控制方法
CN111092444B (zh) 变速恒频风电机组惯性能量支撑控制方法及装置
CN106286128A (zh) 一种减载变速风电机组的系统频率控制方法及装置
CN111864769A (zh) 考虑风机和系统频率响应特性的调频参数确定方法及系统
Ciri et al. Evaluation of log‐of‐power extremum seeking control for wind turbines using large eddy simulations
CN114465287A (zh) 一种风电场快速有功功率优化方法及系统
Akbari Participation of DFIG based wind turbines in improving short term frequency regulation
CN117251995A (zh) 基于可变遗忘因子最小二乘法的双馈风机惯量评估方法
Xu et al. Optimal frequency controller parameters of wind turbines participating system frequency control
CN103603771B (zh) 一种风力发电机组最佳风能捕获方法
De Zutter et al. Modeling of active yaw systems for small and medium wind turbines
Senapati et al. Lagrange interpolating polynomial–based deloading control scheme for variable speed wind turbines
CN114597964A (zh) 一种风光电置换火电的调整方法及装置
Sarfejo et al. Mppt approach for fixed speed wind turbine by using proportional-integral controller
CN113036779B (zh) 一种风电参与电网调频的时变调频参数设定方法及系统
CN115241922A (zh) 风电场参与电网频率支撑的机组组合确定方法与系统
CN109888838B (zh) 风电机组参与电网紧急控制的可调功率计算方法及系统
CN114251235A (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