CN110397553B - 一种不基于模型的风电场尾流管理方法及系统 - Google Patents
一种不基于模型的风电场尾流管理方法及系统 Download PDFInfo
- Publication number
- CN110397553B CN110397553B CN201910682535.4A CN201910682535A CN110397553B CN 110397553 B CN110397553 B CN 110397553B CN 201910682535 A CN201910682535 A CN 201910682535A CN 110397553 B CN110397553 B CN 110397553B
- Authority
- CN
- China
- Prior art keywords
- unit
- wind
- value
- fatigue load
- standard deviation
- 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.)
- Active
Links
- 238000007726 management method Methods 0.000 title claims abstract description 32
- 238000005457 optimization Methods 0.000 claims abstract description 140
- 238000000034 method Methods 0.000 claims abstract description 52
- 230000008569 process Effects 0.000 claims abstract description 32
- 238000004422 calculation algorithm Methods 0.000 claims description 17
- 238000004364 calculation method Methods 0.000 claims description 17
- 238000010606 normalization Methods 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 6
- 230000006870 function Effects 0.000 claims description 2
- 238000010248 power generation Methods 0.000 abstract description 6
- 230000008859 change Effects 0.000 description 4
- 238000001514 detection method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 231100000668 minimum lethal dose Toxicity 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000013135 deep learning Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
Images
Classifications
-
- 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
- F03D7/00—Controlling wind motors
- F03D7/02—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor
- F03D7/0204—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor for orientation in relation to wind direction
-
- 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
- F03D7/00—Controlling wind motors
- F03D7/02—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor
- F03D7/022—Adjusting aerodynamic properties of the blades
- F03D7/0224—Adjusting blade pitch
-
- 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
- F03D7/00—Controlling wind motors
- F03D7/02—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor
- F03D7/028—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor controlling wind motor output power
-
- 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
- F03D7/00—Controlling wind motors
- F03D7/02—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor
- F03D7/028—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor controlling wind motor output power
- F03D7/0292—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor controlling wind motor output power to reduce fatigue
-
- 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
- F03D7/00—Controlling wind motors
- F03D7/02—Controlling wind motors the wind motors having rotation axis substantially parallel to the air flow entering the rotor
- F03D7/04—Automatic control; Regulation
- F03D7/042—Automatic control; Regulation by means of an electrical or electronic controller
- F03D7/048—Automatic control; Regulation by means of an electrical or electronic controller controlling wind farms
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Sustainable Development (AREA)
- Sustainable Energy (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- Wind Motors (AREA)
Abstract
本发明公开了一种不基于模型的风电场尾流管理方及系统,实现对尾流的精确控制,优化机组运行状态,提升整个风场的发电效率。该方法包括以下步骤:获取当前工况下风电场中各机组信息,计算当前工况下的风场总功率和机组疲劳载荷标准差;计算每台机组的偏航角度偏置值或者变桨角度偏置值,将得到的每台机组的偏航角度偏置值或者变桨角度偏置值与当前工况下对应机组的风向偏差或变桨角度相叠加,得到当前工况下各机组最终的风向偏差或变桨角度;基于当前工况下的风场总功率和机组疲劳载荷标准差,计算当前工况下的优化目标,将得到的优化目标与优化过程中的目标最优值进行对比,确定控制参数。
Description
技术领域
本公开涉及新能源风力发电领域,具体涉及一种不基于模型的风电场尾流管理方法及系统。
背景技术
风力机从风中获取能量的同时在其下游形成风速下降的尾流区,若下游风机位于尾流区域,则下游机组的出力会受到影响,从而导致整个风电场的总功率受到影响。
目前针对风电场出力优化的研究多以单台机组的控制为主,风电机组普遍采用最大风能捕获策略来提高出力。由于实际风电场中存在尾流遮挡区,前一排机组会对后一排甚至后面几排机组产生遮挡,使得下游区域机组处于风速下降的尾流区,损失部分发电量。有数据显示,根据地形地貌、机组间的距离和风的湍流强度的不同,由于尾流效应导致的发电量损失最小是2%,最大可达30%。因此,有专家提出一些风电场最大出力优化控制方案,在优化尾流效应的背景下协调分配所有机组的有功功率参考值,满足风电场运行和电力系统调度的要求。发明人在研发过程中发现,目前的研究多是基于Jensen尾流模型或在该模型的基础上进行的改进,若要实现对尾流的精确控制,还要建立与实际风场高一致性的风电场模型,实现难度较大,由于地形、环境条件等的影响,模型的精度也很难把控。
发明内容
为了克服上述现有技术的不足,本公开提供了一种不基于模型的风电场尾流管理方法及系统,实现对尾流的精确控制,优化机组运行状态,提升整个风场的发电效率。
本公开一方面提供的一种不基于模型的风电场尾流管理方法的技术方案是:
一种不基于模型的风电场尾流管理方法,该方法包括以下步骤:
获取当前工况下风电场中各机组信息,包括机组的疲劳载荷、实发功率、变桨角度和风向偏差,计算当前工况下的风场总功率和机组疲劳载荷标准差;
计算每台机组的偏航角度偏置值或者变桨角度偏置值,将得到的每台机组的偏航角度偏置值或者变桨角度偏置值与当前工况下对应机组的风向偏差或变桨角度相叠加,得到当前工况下各机组最终的风向偏差或变桨角度;
基于当前工况下的风场总功率和机组疲劳载荷标准差,计算当前工况下的优化目标,将得到的优化目标与优化过程中的目标最优值进行对比,若优化目标大于优化过程中的目标最优值,则将当前工况下风向偏差或变桨角度值作为控制参数,若优化目标小于优化过程中的目标最优值,则将上一次循环得到的风向偏差或变桨角度值作为控制参数;
循环执行上述步骤,直到满足结束条件,将最后得到的控制参数作为最优参数进行机组控制。
本公开另一方面提供的一种不基于模型的风电场尾流管理系统的技术方案是:
一种不基于模型的风电场尾流管理系统,该系统包括:
数据获取模块,用于获取当前工况下风电场中各机组信息,包括机组的疲劳载荷、实发功率、变桨角度和风向偏差,计算当前工况下的风场总功率和机组疲劳载荷标准差;
角度计算模块,用于计算每台机组的偏航角度偏置值或者变桨角度偏置值;将得到的每台机组的偏航角度偏置值或者变桨角度偏置值与当前工况下对应机组的风向偏差或变桨角度相叠加,得到当前工况下各机组最终的风向偏差或变桨角度;
优化目标计算模块,用于基于当前工况下的风场总功率和机组疲劳载荷标准差,计算当前工况下的优化目标;
控制参数确定模块,用于将得到的优化目标与优化过程中的目标最优值进行对比,若优化目标大于优化过程中的目标最优值,则将当前工况下风向偏差或变桨角度值作为控制参数,若优化目标小于优化过程中的目标最优值,则将上一次循环得到的风向偏差或变桨角度值作为控制参数。
本公开另一方面提供的一种计算机可读存储介质的技术方案是:
一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现如下步骤:
获取当前工况下风电场中各机组信息,包括机组的疲劳载荷、实发功率、变桨角度和风向偏差,计算当前工况下的风场总功率和机组疲劳载荷标准差;
计算每台机组的偏航角度偏置值或者变桨角度偏置值,将得到的每台机组的偏航角度偏置值或者变桨角度偏置值与当前工况下对应机组的风向偏差或变桨角度相叠加,得到当前工况下各机组最终的风向偏差或变桨角度;
基于当前工况下的风场总功率和机组疲劳载荷标准差,计算当前工况下的优化目标,将得到的优化目标与优化过程中的目标最优值进行对比,若优化目标大于优化过程中的目标最优值,则将当前工况下风向偏差或变桨角度值作为控制参数,若优化目标小于优化过程中的目标最优值,则将上一次循环得到的风向偏差或变桨角度值作为控制参数;
循环执行上述步骤,直到满足结束条件,将最后得到的控制参数作为最优参数进行机组控制。
本公开另一方面提供的一种处理装置的技术方案是:
一种处理装置,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征是,所述处理器执行所述程序时实现如下步骤:
获取当前工况下风电场中各机组信息,包括机组的疲劳载荷、实发功率、变桨角度和风向偏差,计算当前工况下的风场总功率和机组疲劳载荷标准差;
计算每台机组的偏航角度偏置值或者变桨角度偏置值,将得到的每台机组的偏航角度偏置值或者变桨角度偏置值与当前工况下对应机组的风向偏差或变桨角度相叠加,得到当前工况下各机组最终的风向偏差或变桨角度;
基于当前工况下的风场总功率和机组疲劳载荷标准差,计算当前工况下的优化目标,将得到的优化目标与优化过程中的目标最优值进行对比,若优化目标大于优化过程中的目标最优值,则将当前工况下风向偏差或变桨角度值作为控制参数,若优化目标小于优化过程中的目标最优值,则将上一次循环得到的风向偏差或变桨角度值作为控制参数;
循环执行上述步骤,直到满足结束条件,将最后得到的控制参数作为最优参数进行机组控制。
通过上述技术方案,本公开的有益效果是:
(1)本公开通过风电场的协同尾流管理控制,可以改善机组之间的尾流影响,从而提高风电场整体发电性能;
(2)本公开采用机组的变桨角度作为控制参数,改变尾流状态;采用机组的偏航角度,来改变尾流尾迹,从而减小尾流对下游风机的影响;可选择上述控制参数中的任意一个或者对两个控制参数同时进行处理,实现对风电场进行尾流管理;
(3)本公开的优化目标为风电场总功率最大化和机组疲劳载荷均衡性最优化等;优化目标也可以为功率和疲劳载荷的协同优化,也可以为功率或疲劳载荷的单独优化,能够得到最优的控制参数,实现对机组优化控制。
附图说明
构成本公开的一部分的说明书附图用来提供对本公开的进一步理解,本公开的示意性实施例及其说明用于解释本申请,并不构成对本公开的不当限定。
图1是实施例一风电场尾流管理方法的流程图;
图2是实施例二风电场尾流管理方法的流程图;
图3是实施例三风电场尾流管理系统的结构图。
具体实施方式
下面结合附图与实施例对本公开作进一步说明。
应该指出,以下详细说明都是例示性的,旨在对本公开提供进一步的说明。除非另有指明,本公开使用的所有技术和科学术语具有与本公开所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本申请的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
实施例一
正如背景技术所介绍的,为了优化风电场整体的运行状态,本实施例提出了一种不基于模型的风电场尾流管理方法,实现了对风电场尾流控制,可提升风电场整体发电量、优化机组运行状态,避免基于模型的方法由于建模误差导致的运算精度偏低等问题。
请参阅附图1,该风电场尾流管理方法包括以下步骤:
S101,获取当前工况下风电场中各机组的信息和风电场风况信息,计算当前工况下的风场总功率和机组疲劳载荷均衡性,并将风电场中各机组信息、风场总功率和机组疲劳载荷均衡性等参数保存在数据库中。
在本实施例中,所述机组信息包括机组的疲劳载荷、实发功率、变桨角度、偏航角度等,其中,机组疲劳载荷从机组载荷估计器或利用应变片现场实测获得,其他参数均可从风电场Scada系统和风电机组控制器中获得。
在本实施例中,所述风电场风况信息包括风速、风向、湍流度等。
具体地,所述风场总功率的计算方法为:
将风电场中每台机组的功率相加,得到风场总功率。风场总功率的表达式为:
其中,P为风场总功率;pi为风场中机组i的功率,N为风场中机组总台数。
在本实施例中,机组疲劳载荷可以为塔筒、叶根、轴承等部位的载荷,机组疲劳载荷均衡性是指风场中所有机组的特定结构或部件的疲劳载荷的标准差值,是衡量所有机组使用寿命差异的指标。
具体地,所述机组疲劳载荷均衡性(所有机组的疲劳载荷的标准差值)的计算方法为:
利用标准差计算函数(STD)计算所有机组的疲劳载荷的标准偏差,其计算公式为:
LD=STD(l1,l2,l3…lN)
其中,LD为风电场所有机组疲劳载荷的标准差值,该值可以体现风场中不同机组疲劳载荷分布的均匀程度,LD值越小说明各机组疲劳载荷之间疲劳载荷分布越均匀;N为风场中机组总台数,li为机组i的疲劳载荷。
具体地,控制参数为机组的变桨角度或偏航角度,可以单独对变桨或者偏航角度进行操作,也可以同时对两者进行操作,根据实际情况进行选择;可以优化的控制参数包括但不限于以上两者。变桨角度,其控制变量是变桨设定值,通过控制变桨角度来改变轴向诱导因子,进而控制风机下游尾流状态;偏航角度,其控制变量是风向偏差,通过给定偏航角度偏置来改变风机的对风角度,进而控制风机下游的尾流尾迹。
在本实施例中,所述机组的偏航角度偏置值或者变桨角度偏置值的计算方法为:
设定随机数,判断随机数与偏航偏置角度更新概率的大小;
若随机数小于偏航偏置角度更新概率,则更新随机数,并利用新的随机数、偏航偏置角度调整最小步长以及机组的最小风向偏差或变桨角度,计算机组的偏航角度偏置值或者变桨角度偏置值。
S103,将步骤102计算得到的每台机组的偏航角度偏置值或者变桨角度偏置值与当前工况下对应机组的风向标检测风向或变桨角度相叠加,得到当前工况下各机组最终的风向偏差或变桨角度。
S104,分别设置风电场总功率和机组疲劳载荷均衡性的权值,计算优化目标。
优化目标为风电场总功率最大化或机组疲劳载荷均衡性最优化,可以进行风电场总功率最大化或机组疲劳载荷均衡性最优化的单目标优化,也可以进行风电场总功率最大化和机组疲劳载荷均衡性最优化的协同优化。
在本实施例中,采用风电场总功率最大化和机组疲劳载荷均衡性最优化的协同优化。采用尾流管理算法计算机组不同运行状态下的整体最优发电效率以及机组疲劳载荷均衡性。具体为:
归一化。将总功率和疲劳载荷标准差分别归一化到(0,1)区间内,避免因两个优化目标的数量级或数值相差太大而导致算法偏向某一个目标的优化而导致的多目标优化无意义。归一化方法如下:
约定单台机组的最大疲劳载荷标准差为LDmax,最小疲劳载荷标准差为LDmin。那么风场机组疲劳载荷标准差的归一化如下:
LDNorm=(LD-LDmin)/(LDmax-LDmin)
将归一化后的风电场总功率和机组疲劳载荷标准差分别乘以权值后求差值,得到优化目标Φ。优化目标Φ为:
Φ=α*PNorm-β*LDNorm
其中,α为风电场总功率的权值;β为机组疲劳载荷标准差的权值,PNorm为机组功率的归一化值;LDNorm为风场机组疲劳载荷标准差的归一化值。
根据风电场实际需要决定权值的大小关系,若以优化功率为主,则设置权值α>0.5,若以优化机组疲劳载荷为主则设置权值β>0.5。α、β的具体大小还需要根据对优化目标的重视程度来确定。
所述尾流管理算法可以为博弈论算法、遗传算法、粒子群算法、蚁群算法、机器学习算法、深度学习算法等不基于被控对象模型的优化算法,根据实际情况进行选择。
在本实施例中,采用机组疲劳载荷均衡性最优化进行单独优化时,则将机组疲劳载荷标准差的负值-LD作为优化目标Φ。
S105,将步骤104的计算得到的优化目标Φ与设置的优化过程中的目标最优值进行对比,若优化目标提升,即则记录叠加后的风向偏差或变桨角度值作为控制参数,若优化目标没有提升,即则记录上一次循环中的各机组的风向偏差或桨距角作为控制参数。
S106,按照上述步骤进行寻优,直到满足结束条件,结束条件可为达到循环次数或者优化目标收敛。将最后记录的控制参数作为最优参数进行机组控制。
本实施例提出的不基于模型的风电场尾流管理方法。采用机组的变桨角度作为控制参数,改变尾流状态;采用机组的偏航角度,来改变尾流尾迹,从而减小尾流对下游风机的影响;可选择上述控制参数中的任意一个或者对两个控制参数同时进行处理,实现对风电场进行尾流管理;计算的优化目标为风电场总功率最大化和机组疲劳载荷均衡性最优化等;优化目标也可以为功率和疲劳载荷的协同优化,也可以为功率或疲劳载荷的单独优化,能够得到最优的控制参数,实现对机组优化控制。
实施例二
为了使本领域的技术人员更好地了解本公开的技术方案,下面列举一个具体的实施例。
本实施例提供一种不基于模型的风电场尾流管理方法,如图2所示,该方法利用博弈论算法进行尾流的管理,该方法的主要步骤如下:
S201:以偏航角度为控制系数,进行如下定义:
2)k=0为优化迭代次数,E∈[0,1]偏航偏置角度更新概率;
4)Zn∈[0,1,2...n]为正整数,Δr为偏航偏置角度调整最小步长;
S202,优化目标处理,对风电场总功率和机组疲劳载荷进行协同优化,令为风电场总功率,LD=STD(l1,l2,l3…lN)为风电场所有机组疲劳载荷的标准差值,该值可以体现风场中不同机组疲劳载荷分布的均匀程度,LD值越小说明各机组疲劳载荷之间疲劳载荷分布越均匀。
进行多目标优化,将两个优化目标进行归一化处理,将总功率和疲劳载荷标准差分别归一化到(0,1)区间内,避免因两个优化目标的数量级或数值相差太大而导致算法偏向某一个目标的优化而导致多目标优化无意义。归一化方法如下:
约定单台机组的最大疲劳载荷标准差为LDmax,最小疲劳载荷标准差为LDmin。那么风场机组疲劳载荷标准差的归一化如下:
LDNorm=(LD-LDmin)/(LDmax-LDmin)
将归一化后的风电场总功率和机组疲劳载荷标准差分别乘以权值后求差值,得到优化目标Φ。优化目标Φ为:
Φ=α*PNorm-β*LDNorm
其中,α为风电场总功率的权值;β为机组疲劳载荷标准差的权值。
根据风电场实际需要决定权值的大小关系,若以优化功率为主,则设置权值α>0.5,若以优化机组疲劳载荷为主则设置权值β>0.5。α、β的具体大小还需要根据对优化目标的重视程度来确定。
S203,获取当前工况下各机组的偏航角度、桨距角、实发功率、机组疲劳载荷等信息,并计算出风场总功率、机组疲劳载荷均衡性,保存在数据库中;
S204,初始化:
S205,令k=k+1,针对所有的机组进行以下操作:
定义R1=random[0,1],为随机数;
若R1<E,则更新
R2=random[Zn],
ri=rmin+R2*Δr
否则,
即,每一次循环,是否对机组施加偏航偏置值由随机数R1和概率E来决定,即在一次循环中并不是每台机组的偏航指令都会发生变化。
S206,将步骤205计算出的偏航偏置值与对应机组的当前风向标检测风向叠加后作为机组最终的偏航角度;
S207,计算当前偏航角度和工况下的优化目标Φ,并将计算结果保存;
S209,判断如果k达到设置次数或者优化目标收敛到某个小的区间内,则停止算法循环,用计算出的控制参数作为最优参数用以进行机组控制;否则返回205。
在步骤201中,控制参数也可以是变桨角度,或者对偏航角度和变桨角度同时进行操作。
在步骤207中,若要单独进行机组疲劳载荷优化,则将步骤207中的优化目标Φ替换为机组疲劳载荷标准差的负值,-LD。
本实施例提出的不基于模型的风电场尾流管理方法。采用机组的偏航角度作为控制参数,改变尾流尾迹,实现对风电场进行尾流管理;优化目标为功率和疲劳载荷的协同优化,能够得到最优的控制参数,实现对机组优化控制。
实施例三
本实施例提供一种不基于模型的风电场尾流管理系统,该系统包括:
数据获取模块,用于获取当前工况下各机组的风电场机组信息,包括风电机组的疲劳载荷、实发功率、变桨角度和风向标检测风向,计算当前工况下的风场总功率和机组疲劳载荷标准差;若要进行风场总功率和疲劳载荷的多目标优化,还需要对两者分别进行归一化处理。
数据存储模块,用于存储当前工况下各机组的风电场机组信息、风场总功率和机组疲劳载荷标准差;
角度计算模块,用于计算每台机组的偏航角度偏置值或者变桨角度偏置值,并设定偏航偏置角度或变桨角度的初始值以及优化过程中的目标最优值;将得到的每台机组的偏航角度偏置值或者变桨角度偏置值与当前工况下对应机组的风向偏差或变桨角度相叠加,得到当前工况下各机组最终的风向偏差或变桨角度;
优化目标计算模块,用于基于当前工况下的风场总功率和机组疲劳载荷标准差,计算当前工况下的优化目标;
控制参数确定模块,用于将当前得到的优化目标与优化过程中的目标最优值进行对比,若优化目标大于优化过程中的目标最优值,则将当前(叠加后的)风向或变桨角度值的作为控制参数,若优化目标小于优化过程中的目标最优值,则将上一次循环中的风向偏差或桨距角作为控制参数;利用控制参数控制机组运行。
本实施例提出的不基于模型的风电场尾流管理系统,采用风电场控制器,不基于尾流模型和风电场模型,利用先进的优化算法,用于实际风电场进行尾流管理,优化风场整体运行状态。
实施例四
本实施例提供一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现如下步骤:
获取当前工况下风电场中各机组信息,包括机组的疲劳载荷、实发功率、变桨角度和风向偏差,计算当前工况下的风场总功率和机组疲劳载荷标准差;
计算每台机组的偏航角度偏置值或者变桨角度偏置值,将得到的每台机组的偏航角度偏置值或者变桨角度偏置值与当前工况下对应机组的风向偏差或变桨角度相叠加,得到当前工况下各机组最终的风向偏差或变桨角度;
基于当前工况下的风场总功率和机组疲劳载荷标准差,计算当前工况下的优化目标,将得到的优化目标与优化过程中的目标最优值进行对比,若优化目标大于优化过程中的目标最优值,则将当前工况下风向偏差或变桨角度值作为控制参数,若优化目标小于优化过程中的目标最优值,则将上一次循环得到的风向偏差或变桨角度值作为控制参数;
循环执行上述步骤,直到满足结束条件,将最后得到的控制参数作为最优参数进行机组控制。
实施例五
本实施例提供一种处理装置,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征是,所述处理器执行所述程序时实现如下步骤:
获取当前工况下风电场中各机组信息,包括机组的疲劳载荷、实发功率、变桨角度和风向偏差,计算当前工况下的风场总功率和机组疲劳载荷标准差;
计算每台机组的偏航角度偏置值或者变桨角度偏置值,将得到的每台机组的偏航角度偏置值或者变桨角度偏置值与当前工况下对应机组的风向偏差或变桨角度相叠加,得到当前工况下各机组最终的风向偏差或变桨角度;
基于当前工况下的风场总功率和机组疲劳载荷标准差,计算当前工况下的优化目标,将得到的优化目标与优化过程中的目标最优值进行对比,若优化目标大于优化过程中的目标最优值,则将当前工况下风向偏差或变桨角度值作为控制参数,若优化目标小于优化过程中的目标最优值,则将上一次循环得到的风向偏差或变桨角度值作为控制参数;
循环执行上述步骤,直到满足结束条件,将最后得到的控制参数作为最优参数进行机组控制。
上述虽然结合附图对本公开的具体实施方式进行了描述,但并非对本公开保护范围的限制,所属领域技术人员应该明白,在本公开的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本公开的保护范围以内。
Claims (8)
1.一种不基于模型的风电场尾流管理方法,其特征是,包括以下步骤:
获取当前工况下风电场中各机组信息,包括机组的疲劳载荷、实发功率、变桨角度和风向偏差,计算当前工况下的风电场总功率和机组疲劳载荷标准差;
计算每台机组的偏航角度偏置值或者变桨角度偏置值,将得到的每台机组的偏航角度偏置值或者变桨角度偏置值与当前工况下对应机组的风向偏差或变桨角度相叠加,得到当前工况下各机组最终的风向偏差或变桨角度;
基于当前工况下的风电场总功率和机组疲劳载荷标准差,计算当前工况下的优化目标,将得到的优化目标与优化过程中的目标最优值进行对比,若优化目标大于优化过程中的目标最优值,则将当前工况下风向偏差或变桨角度值作为控制参数,若优化目标小于优化过程中的目标最优值,则将上一次循环得到的风向偏差或变桨角度值作为控制参数;
循环执行上述步骤,直到满足结束条件,将最后得到的控制参数作为最优参数进行机组控制;
所述机组的偏航角度偏置值或者变桨角度偏置值的计算方法为:
设定随机数,判断随机数与偏航偏置角度更新概率的大小;
若随机数小于偏航偏置角度更新概率,则更新随机数,并利用新的随机数、偏航偏置角度调整最小步长以及机组的最小风向偏差或变桨角度,计算机组的偏航角度偏置值或者变桨角度偏置值;
所述当前工况下的优化目标的计算方法为:
将当前工况下的风电场总功率和疲劳载荷标准差分别归一化到(0,1)区间内,避免因两个优化目标的数量级或数值相差太大而导致算法偏向某一个目标的优化而导致的多目标优化无意义;归一化方法如下:
约定单台机组的最大疲劳载荷标准差为LDmax,最小疲劳载荷标准差为LDmin,则风场机组疲劳载荷标准差的归一化如下:
LDNorm=(LD-LDmin)/(LDmax-LDmin)
将归一化后的风电场总功率和机组疲劳载荷标准差分别乘以权值后求差值,得到优化目标Φ;优化目标Φ为:
Φ=α*PNorm-β*LDNorm
其中,α为风电场总功率的权值;β为机组疲劳载荷标准差的权值,PNorm为机组功率的归一化值;LDNorm为风场机组疲劳载荷标准差的归一化值;
根据风电场实际需要决定权值的大小关系,若以优化功率为主,则设置权值α>0.5,若以优化机组疲劳载荷为主则设置权值β>0.5,α、β的具体大小还需要根据对优化目标的重视程度来确定。
2.根据权利要求1所述的不基于模型的风电场尾流管理方法,其特征是,所述风电场总功率的计算方法为:
将风电场中每台机组的功率相加,得到风电场总功率。
3.根据权利要求1所述的不基于模型的风电场尾流管理方法,其特征是,所述机组疲劳载荷标准差的计算方法为:
利用标准差计算函数计算风场中所有机组的疲劳载荷的标准偏差。
4.根据权利要求1所述的不基于模型的风电场尾流管理方法,其特征是,对当前工况下的风电场总功率和机组疲劳载荷标准差进行归一化处理的步骤包括:
设定单台机组最大功率和最小功率,计算风电场总功率的最大值和风电场总功率的最小值,将当前工况下的风电场总功率和风电场总功率的最小值的差值与风电场总功率的最大值和最小值的差值相比,得到归一化后的风电场总功率;
设定单台机组的最大疲劳载荷标准差和最小疲劳载荷标准差,将当前工况下的机组疲劳载荷标准差和最小疲劳载荷标准差的差值与最大疲劳载荷标准差和最小疲劳载荷标准差的差值相比,得到归一化后的机组疲劳载荷标准差。
5.根据权利要求1所述的不基于模型的风电场尾流管理方法,其特征是,所述结束条件为达到循环次数或者优化目标收敛。
6.一种不基于模型的风电场尾流管理系统,其特征是,包括:
数据获取模块,用于获取当前工况下风电场中各机组信息,包括机组的疲劳载荷、实发功率、变桨角度和风向偏差,计算当前工况下的风电场总功率和机组疲劳载荷标准差;
角度计算模块,用于计算每台机组的偏航角度偏置值或者变桨角度偏置值,将得到的每台机组的偏航角度偏置值或者变桨角度偏置值与当前工况下对应机组的风向偏差或变桨角度相叠加,得到当前工况下各机组最终的风向偏差或变桨角度;
具体为:
设定随机数,判断随机数与偏航偏置角度更新概率的大小;
若随机数小于偏航偏置角度更新概率,则更新随机数,并利用新的随机数、偏航偏置角度调整最小步长以及机组的最小风向偏差或变桨角度,计算机组的偏航角度偏置值或者变桨角度偏置值;
优化目标计算模块,用于基于当前工况下的风电场总功率和机组疲劳载荷标准差,计算当前工况下的优化目标;
具体为:
将当前工况下的风电场总功率和疲劳载荷标准差分别归一化到(0,1)区间内,避免因两个优化目标的数量级或数值相差太大而导致算法偏向某一个目标的优化而导致的多目标优化无意义;归一化方法如下:
约定单台机组的最大疲劳载荷标准差为LDmax,最小疲劳载荷标准差为LDmin,则风场机组疲劳载荷标准差的归一化如下:
LDNorm=(LD-LDmin)/(LDmax-LDmin)
将归一化后的风电场总功率和机组疲劳载荷标准差分别乘以权值后求差值,得到优化目标Φ;优化目标Φ为:
Φ=α*PNorm-β*LDNorm
其中,α为风电场总功率的权值;β为机组疲劳载荷标准差的权值,PNorm为机组功率的归一化值;LDNorm为风场机组疲劳载荷标准差的归一化值;
根据风电场实际需要决定权值的大小关系,若以优化功率为主,则设置权值α>0.5,若以优化机组疲劳载荷为主则设置权值β>0.5,α、β的具体大小还需要根据对优化目标的重视程度来确定;
控制参数确定模块,用于将得到的优化目标与优化过程中的目标最优值进行对比,若优化目标大于优化过程中的目标最优值,则将当前工况下风向偏差或变桨角度值作为控制参数,若优化目标小于优化过程中的目标最优值,则将上一次循环得到的风向偏差或变桨角度值作为控制参数。
7.一种计算机可读存储介质,其上存储有计算机程序,其特征是,该程序被处理器执行时实现如权利要求1-5中任一项所述的不基于模型的风电场尾流管理方法中的步骤。
8.一种处理装置,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征是,所述处理器执行所述程序时实现如权利要求1-5中任一项所述的不基于模型的风电场尾流管理方法中的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910682535.4A CN110397553B (zh) | 2019-07-26 | 2019-07-26 | 一种不基于模型的风电场尾流管理方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910682535.4A CN110397553B (zh) | 2019-07-26 | 2019-07-26 | 一种不基于模型的风电场尾流管理方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110397553A CN110397553A (zh) | 2019-11-01 |
CN110397553B true CN110397553B (zh) | 2020-09-25 |
Family
ID=68325158
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910682535.4A Active CN110397553B (zh) | 2019-07-26 | 2019-07-26 | 一种不基于模型的风电场尾流管理方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110397553B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111980857A (zh) * | 2020-08-26 | 2020-11-24 | 上海电气风电集团股份有限公司 | 风电场的闭环控制方法及其装置及计算机可读存储介质 |
CN112459965B (zh) * | 2020-11-18 | 2022-03-01 | 上海电气风电集团股份有限公司 | 考虑风电场尾流的偏航优化控制方法、装置、设备及介质 |
CN112879220B (zh) * | 2021-03-16 | 2022-11-01 | 上海电气风电集团股份有限公司 | 风机控制方法、系统和可读存储介质 |
CN119532106A (zh) * | 2024-11-25 | 2025-02-28 | 北京宏大天恒科技有限公司 | 多普勒激光测风的风电机组阵列尾流管理方法及系统 |
CN119994986A (zh) * | 2025-04-11 | 2025-05-13 | 山东大学 | 一种考虑载荷均衡的风电场控制优化方法及系统 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1856645A (zh) * | 2003-06-14 | 2006-11-01 | 荷兰能源研究中心基金会 | 用于从流动的流体吸取能量的方法和设备 |
CN101413483A (zh) * | 2007-10-16 | 2009-04-22 | 通用电气公司 | 用于优化风力涡轮之间的尾流交互作用的系统和方法 |
CN101876289A (zh) * | 2009-04-30 | 2010-11-03 | 通用电气公司 | 用于增强具有多个风力涡轮机的风电厂布局的方法 |
CN105556117A (zh) * | 2013-09-17 | 2016-05-04 | 维斯塔斯风力系统集团公司 | 用于风力涡轮机的控制方法 |
EP3096006A1 (en) * | 2015-05-21 | 2016-11-23 | General Electric Company | System and methods for controlling noise propagation of wind turbines |
CN107850048A (zh) * | 2015-06-30 | 2018-03-27 | 维斯塔斯风力系统集团公司 | 用于生成风力涡轮机控制安排的方法以及系统 |
CN108708825A (zh) * | 2018-05-15 | 2018-10-26 | 浙江运达风电股份有限公司 | 一种考虑机组尾流的风电场偏航控制方法 |
CN109340062A (zh) * | 2018-12-18 | 2019-02-15 | 国电联合动力技术有限公司 | 一种低风速风电机组数字双胞胎式的疲劳损伤预测方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9574546B2 (en) * | 2012-06-14 | 2017-02-21 | General Electric Company | Wind turbine rotor control |
US10465655B2 (en) * | 2016-07-05 | 2019-11-05 | Inventus Holdings, Llc | Wind turbine wake steering apparatus |
-
2019
- 2019-07-26 CN CN201910682535.4A patent/CN110397553B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1856645A (zh) * | 2003-06-14 | 2006-11-01 | 荷兰能源研究中心基金会 | 用于从流动的流体吸取能量的方法和设备 |
CN101413483A (zh) * | 2007-10-16 | 2009-04-22 | 通用电气公司 | 用于优化风力涡轮之间的尾流交互作用的系统和方法 |
CN101876289A (zh) * | 2009-04-30 | 2010-11-03 | 通用电气公司 | 用于增强具有多个风力涡轮机的风电厂布局的方法 |
CN105556117A (zh) * | 2013-09-17 | 2016-05-04 | 维斯塔斯风力系统集团公司 | 用于风力涡轮机的控制方法 |
EP3096006A1 (en) * | 2015-05-21 | 2016-11-23 | General Electric Company | System and methods for controlling noise propagation of wind turbines |
CN107850048A (zh) * | 2015-06-30 | 2018-03-27 | 维斯塔斯风力系统集团公司 | 用于生成风力涡轮机控制安排的方法以及系统 |
CN108708825A (zh) * | 2018-05-15 | 2018-10-26 | 浙江运达风电股份有限公司 | 一种考虑机组尾流的风电场偏航控制方法 |
CN109340062A (zh) * | 2018-12-18 | 2019-02-15 | 国电联合动力技术有限公司 | 一种低风速风电机组数字双胞胎式的疲劳损伤预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110397553A (zh) | 2019-11-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110397553B (zh) | 一种不基于模型的风电场尾流管理方法及系统 | |
EP3121442B2 (en) | Operating wind turbines | |
US9201410B2 (en) | Methods and systems for optimizing farm-level metrics in a wind farm | |
US11713747B2 (en) | Wind turbine control method | |
CN102709939B (zh) | 一种提高风电场发电效率的风电场有功功率控制方法 | |
CN113054688B (zh) | 考虑不确定性的可再生能源集群出力分布式协调控制方法 | |
CN114465287A (zh) | 一种风电场快速有功功率优化方法及系统 | |
CN113394813B (zh) | 海上风电场的机组功率指令值计算方法和分布式调度方法 | |
CN113011994A (zh) | 风电厂的优化 | |
CN119844296A (zh) | 基于环境数据的风机运行实时动态调节控制方法及系统 | |
CN112879220B (zh) | 风机控制方法、系统和可读存储介质 | |
CN114294157A (zh) | 风电风机偏航偏差辨识方法、纠偏方法、电子设备和介质 | |
CN110970936B (zh) | 一种深度调峰机组一次调频性能计算方法 | |
CN110645145A (zh) | 风力发电机组的控制方法和控制设备 | |
CN118630827A (zh) | 风电场与共享储能有功调频裕度的综合估计方法及系统 | |
CN111931967A (zh) | 一种风电场短期功率预测方法 | |
CN117251995A (zh) | 基于可变遗忘因子最小二乘法的双馈风机惯量评估方法 | |
CN113250917B (zh) | 海上风机阵列输出指令控制方法、系统、装置及存储介质 | |
CN116169713A (zh) | 一种风电/抽蓄组合并网电力系统的低频振荡抑制方法 | |
CN114649831A (zh) | 一种基于分布式控制的风电场功率优化控制方法及系统 | |
CN107069767B (zh) | 一种提高无功补偿装置控制精度的方法及系统 | |
CN109672224B (zh) | 一种风电场有功功率控制方法及系统 | |
CN119026342B (zh) | 风电机组叶片气动优化方法、设备及介质 | |
CN118842100B (zh) | 一种风电集群控制方法、系统及存储介质 | |
CN118412935B (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 |