CN115807734A - 一种基于尾流跟踪的海上风电场场级协同控制策略 - Google Patents

一种基于尾流跟踪的海上风电场场级协同控制策略 Download PDF

Info

Publication number
CN115807734A
CN115807734A CN202211143225.3A CN202211143225A CN115807734A CN 115807734 A CN115807734 A CN 115807734A CN 202211143225 A CN202211143225 A CN 202211143225A CN 115807734 A CN115807734 A CN 115807734A
Authority
CN
China
Prior art keywords
wind
wake
wake flow
center position
power plant
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.)
Pending
Application number
CN202211143225.3A
Other languages
English (en)
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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN202211143225.3A priority Critical patent/CN115807734A/zh
Publication of CN115807734A publication Critical patent/CN115807734A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/72Wind turbines with rotation axis in wind direction

Landscapes

  • Wind Motors (AREA)

Abstract

本发明公开了一种基于尾流跟踪的海上风电场场级协同控制策略,包括机舱式激光雷达测风系统、尾流跟踪模块、优化器和场级控制器。首先利用机舱式激光雷达测风系统测得原始风信息,通过尾流跟踪模块进行风场反演完成环境输入风的参数评估和尾流特征参数辨识;基于海上风电场气动‑液压‑伺服‑弹性动力学仿真模型,通过智能优化算法求解不同环境工况下的最优尾流中心位置,建立多维智能决策数据库LUT和优化器;利用场级协同PI控制器实现尾流重定向和智能化控制。本发明的场级协同控制方法可以实现尾流的在线实时跟踪,动态改善尾流效应,实现海上风电场的增效降载;所述控制器具有控制效果好、执行速度快等特点。

Description

一种基于尾流跟踪的海上风电场场级协同控制策略
技术领域
本发明属于海上风电场协同控制技术领域,尤其是涉及一种基于尾流跟踪的海上风电场场级协同控制策略。
背景技术
随着经济社会的快速发展,电力需求大幅增长,需要大力发展可再生能源。海上风能作为可再生能源的重要一环,储量丰富,而且近海风电与用电密集的沿海城市距离较近,便于电能传输与消纳,可以有效节约陆地资源,降低视觉和噪音污染。因此,海上风力发电被认为是解决能源需求和环境问题的有效方式之一。
然而,现有的海上风电场往往由于海域限制排布紧密,从而形成尾流效应带来功率损耗和疲劳载荷。研究表明,尾流效应能够对下游风力机造成40%的功率亏损以及80%的疲劳载荷增长,最终降低风电场产能和疲劳寿命,影响整个风电场的经济效益和安全运行。
为了改善尾流效应,现有研究提出海上风电场场级协同控制策略,利用先进传感和控制方法,实现风力机的协同运行控制,降低尾流影响,实现风电场场级的功率提升和降载,具体控制方法包括尾流重定向控制和轴向诱导因子控制。
公开号为CN111980857A的中国专利文献公开了一种风电场的闭环控制方法,包括:通过风电场的工程尾流模型来计算得到风电场内每台风机不同来流下的初步最优偏航角;测量风电场的现场尾流流动数据;以及基于风电场的现场尾流流动数据作为反馈来修正每台风机不同来流下的初步最优偏航角。
公开号为CN108953060A的中国专利文献公开了一种基于激光雷达测风仪的风电场场级偏航控制方法,包括:通过设于风电场内预定位置和预定数量的激光雷达测量相应位置的风速风向数据;处理所述风速数据以得到风电场局部单机尾流模型;通过所述局部尾流模型,得出整个风电场的尾流模型;根据整个风电场的尾流模型,计算风电机组的当前偏航角度修正因子;通过前排激光雷达测量提前风向信息,作为第一排风电机组偏航的补偿信号;基于风速风向标的风向信号和所述当前偏航修正信号及所述补偿信号进行风电机组偏航控制。
然而,大气风场具有复杂性,尾流场结构多变,利用现有的技术方法难以实现海上风电场内精准的尾流评估,目前,现有海上风电场协同控制技术并没有形成成熟的尾流智能化控制能力。
发明内容
本发明提供了一种基于尾流跟踪的海上风电场场级协同控制策略,通过激光雷达测风技术实现尾流实时跟踪,并利用场级协同控制技术实现尾流重定向和智能化控制,从而提高海上风电场系统的发电量,实现海上风电场的增效降载。
一种基于尾流跟踪的海上风电场场级协同控制策略,包括:
(1)利用机舱式激光雷达测风系统测量得到每台海上风力机前侧流场风速分布数据;
(2)建立尾流跟踪模块,利用步骤(1)得到的风速分布数据通过三维风场反演完成来流风在无尾流影响下的风信息评估和尾流影响下的尾流特征参数辨识,尾流特征参数包括尾流深度、尾流宽度和尾流中心位置;
(3)根据实际海上风电场的排布特征和流场分布情况,建立海上风电场气动-液压-伺服-弹性动力学仿真模型;利用激光雷达测风系统所测的历史风速分布数据对仿真模型的参数进行校正,并将该仿真模型用于预测风电场在输入风况和控制策略下的运行状况;
(4)建立海上风电场基于尾流重定向的多维智能决策数据库,用于实现在不同环境风况下的最优控制策略在线实时搜索;
建立综合考虑风电场功率最大化和载荷最小化的多目标优化函数,通过智能优化算法得到风电场在不同环境风况输入下的最优尾流中心位置信息,并将上述参数信息记录到多维智能决策数据库;
(5)建立海上风电场场级协同PI控制器,利用查表法在线搜索多维智能决策数据库获得风电场在相应风况输入下的最优尾流中心位置信息,将其作为场级协同PI控制器的参考值,根据尾流跟踪模块测得的尾流信息测量值,通过场级协同PI控制器改变各机组与环境来流风的偏航不对中程度,实现尾流重定向和控制优化。
进一步地,尾流跟踪模块旨在利用海上风电场机舱式激光雷达测量获得前侧流场风速分布情况,通过风向对北校正、水平风速换算、三维风场反演及尾流参数辨识等数据处理完成来流风在无尾流影响下的风信息评估和尾流影响下的尾流特征参数辨识,获得尾流参数信息。
步骤(2)的具体过程为:
(2-1)三维风场反演
海上风电场机舱式激光雷达测风系统可精确地远程测量其安装位置前方的风况,获得风力机的风廓线和复杂的尾流状态,主要用于风电场性能优化和整体控制。在本发明所述的风电场场级协同控制过程中,海上风力机的机舱往往还面临着空间姿态位置随时间不断变化等问题,这给激光雷达在风向测量上带来难度,因此需要在三维风场反演时设计相应的风速算法来消除误差。
为了准确评估海上风电场内部风速分布信息,首先进行风向对北校正,通过坐标旋转转换矩阵算法来反演激光束在地球坐标系下的指向;激光束指向的坐标变换公式为:
LLOS=(H1H2H3)-1LLOS_Lidar (1)
其中,LLOS为激光束在地球坐标下的指向,H1、H2、H3为坐标旋转转换矩阵,LLOS_Lidar为激光束在机舱坐标下的指向;
定义机舱坐标系的方向:X0为转子轴向指向机舱头部方向,Y0为转子径向指向水平左侧,Z0与风力机塔筒轴同轴并指向机舱顶部方向;激光束指向通过方位角
Figure BDA0003854289950000041
和俯仰角θ0描述,其中方位角
Figure BDA0003854289950000042
为激光束在X0Y0面的投影与X0的夹角,俯仰角θ0为激光束与X0Y0面的夹角;风力机机舱在协同控制过程中形成姿态的变化,用翻滚角
Figure BDA0003854289950000043
俯仰角θ和偏航角ψ表示;因此,激光束在机舱坐标下的指向LLOS_Lidar及坐标旋转转换矩阵表示为如下:
Figure BDA0003854289950000044
Figure BDA0003854289950000045
Figure BDA0003854289950000046
Figure BDA0003854289950000047
(1)式可由(2)、(3)、(4)、(5)式带入直接得到激光束实际指向,进一步求解得到激光束视线风速vLOS
vLOS=vwind·LLOS (6)
在实际海上机舱式激光雷达扫描测量时,可获得任意高度上多个方位的视线风速信息。在相同高度和圆周扫描层,视线风速在理想情况下为三角函数,定义为:
f(vLOS)=a cos(θ-b)+c (7)
参数a、b和c由所测量高度上的风矢量信息决定,具体地:
Figure BDA0003854289950000051
其中,v1、v2和v3为测量风电场在地球坐标系下的南北方向、东西方向和垂直方向上的分速度,其决定了上述视线风速函数的参数大小;
通过最小二乘法拟合(7)式和(8)式所述的视线风速,通过最小化函数误差(9)式,求解得到曲线函数的参数a、b和c,完成三维风电场风速分布的反演;
Figure BDA0003854289950000052
其中,N为机舱式激光雷达在某个测量周期内的测试点个数;基于上述风场反演,获得风电场来流风在无尾流影响下的风信息,包括风速、风向和湍流强度,同时为尾流特征参数辨识提供数据基础;
(2-2)尾流特征参数辨识
基于海上风电场机舱式激光雷达测得原始数据,通过风向对北校正、水平风速换算完成三维风场反演,获得风力机前侧风速分布情况。利用风力机轮毂高度处的风速数据,完成基于反高斯分布的尾流参数辨识,获取尾流的特征参数。具体地,建立风力机前侧所受尾流速度分布的反高斯函数,公式如下:
Figure BDA0003854289950000053
其中,U为风力机轮毂高度的(x,zh)位置处在时间t时刻的风速,U0为时间t时刻的环境平均风速,Cs为尾流比例系数,用以量化尾流深度,e为自然常数,μ为轮毂高度处尾流中心位置,σ为尾流标准差,用以量化尾流宽度,π为圆周率;基于海上风电场的排布情况,本发明假设尾流中心位置约位于下游风力机轮毂高度处。
利用激光雷达测得的风力机前侧相应距离处的风速数据,利用最小二乘法完成反高斯函数的参数拟合,获得尾流特征参数,具体包括尾流深度、尾流宽度和尾流中心位置信息。
步骤(3)中,海上风电场气动-液压-伺服-弹性动力学仿真模型可用于预测海上风电场在不同环境工况和控制输入下的动态运行特征,具体包括流场分布和机组功率、载荷等动态响应等。
海上风电场气动-液压-伺服-弹性动力学仿真模型具体包括风电场气动仿真模型、海上风力机模型和尾流模型等。
所述的风电场气动仿真模型利用Taylor湍流冻结假设完成风电场三维流场在空间上的动态传播模拟;所述的海上风力机模型采用制动盘模型,根据实际机组的性能参数完成建模并获得在特定环境工况下机组功率、载荷动态响应;所述的尾流模型用于预测与风电场输出功率和机组负载相关的关键尾流特征,包括尾流风速赤字、尾流膨胀以及尾流蜿蜒等。
步骤(4)中,海上风电场基于尾流重定向的多维智能决策数据库LUT旨在记录海上风电场在不同环境工况下各机组前侧的最佳尾流中心位置信息。本发明所述尾流重定向是通过风电场主动偏航协同控制进行实现,通过上游风力机的不对风实现下游尾流轨迹的变化,从而改善下游风力机由于尾流效应带来的功率损耗和载荷增长。为了在实际海上风电场动态变化的环境工况下实现多目标优化,本发明建立的多维智能决策数据库LUT可用于实现最优尾流中心位置的在线实时搜索,从而利用协同偏航控制实现尾流重定向。
具体地,为了建立所述海上风电场基于尾流重定向的多维智能决策数据库LUT,利用步骤(3)所述海上风电场气动-液压-伺服-弹性动力学仿真模型获得海上风电场在不同环境工况和控制输入下的流场分布和机组功率、载荷等动态响应等。根据海上风电场系统在单位时间内功率和载荷的时域变化,分别计算总输出功率和机组关键部位疲劳载荷。
建立综合考虑风电场功率最大化和载荷最小化的多目标优化函数具体过程为:
(4-1)计算海上风力发电机组i的输出功率:
Figure BDA0003854289950000071
式中,Pi为风电场内机组i发电机的输出功率(W),Tq,i(t)为机组i发电机在t时刻的瞬时扭矩(N-m),ωi为机组i发电机在t时刻的瞬时转速(rpm),tk为当前时刻,ΔT为当前求解的时域周期。在本发明专利中可通过主动偏航控制改变环境风入流角,实现风电场运行性能的优化。
根据上述单台机组的功率计算公式,求解风电场系统的总发电功率即为风电场系统内所有风力发电机组的发电功率之和
Figure BDA0003854289950000072
为:
Figure BDA0003854289950000073
(4-2)计算风力发电机组i所受等效疲劳载荷:
利用风电场气动-液压-伺服-弹性动力学仿真模型获得风电场内风力发电机组i关键部位所受载荷时域变化情况,包括塔筒底部和叶片根部的载荷变化情况;通过雨流计数法,对风力发电机组i所受时域变化的载荷进行等效分析,获得关键部位等效疲劳载荷DEL:
Figure BDA0003854289950000074
式中,DEL为等效疲劳载荷,
Figure BDA0003854289950000075
为时间序列j内的等效循环次数,Nj,i为第i个工况在时间序列j内发生的次数,Lj,i为第i个工况在时间序列j内的载荷范围,m为材料S-N曲线的斜率;
(4-3)为实现风电场总输出功率最大化,并降低机组关键部位所受疲劳载荷,建立多目标优化函数如下:
Figure BDA0003854289950000076
式中,
Figure BDA0003854289950000077
为归一化后的风电场系统总发电功率,DELnorm为归一化后的各台风力发电机组关键部位所受最大疲劳载荷,α为权重系数。
建立多目标优化函数后,通过智能优化算法求解在特定环境工况下最大化目标函数J的最优偏航控制角,并通过风电场气动-液压-伺服-弹性动力学仿真模型记录各机组在该偏航控制输入下转子前侧的最优尾流中心位置;最后,将上述环境风况参数和对应的最优尾流中心位置信息保存为多维智能决策数据库,用以不同环境风况下各机组前侧最优尾流中心位置的在线实时搜索和场级协同控制。
具体地,步骤(4-3)中,所述的智能优化算法包括但不限于遗传算法、粒子群算法和博弈论算法。
步骤(5)中,所述的海上风电场场级协同PI控制器包含低通滤波器、多维智能决策数据库、尾流中心位置误差计算模块、PI控制器、偏航执行器和激光雷达测风系统。
将激光雷达测风系统获得的环境风况测量值作为输入,通过低通滤波器获得外界环境风输入的平均风速与风向,在低通滤波器后接多维智能决策数据库;利用查表法实现在线搜索,获得对应环境风况下的各机组前侧最优尾流中心位置信息;利用尾流中心位置误差计算模块完成各机组前侧实际尾流中心位置与期望尾流中心位置的误差估计;将尾流中心位置误差计算模块获得的中心位置偏差值输入PI控制器,对偏差值进行比例、积分运算,得到风电场在单位控制周期内进行主动协同偏航控制的调整值;通过偏航执行器改变各机组与环境来流风的偏航不对中程度,实现尾流重定向和控制优化。
利用尾流中心位置误差计算模块完成各机组前侧实际尾流中心位置与期望尾流中心位置的误差估计具体包括:
通过激光雷达测风系统的三维风场反演获得实测尾流中心位置信息参数,通过在线搜索多维智能决策数据库获取最优尾流中心位置信息参数,通过尾流中心位置误差计算模块对两者做差获得中心位置的偏差值。
与现有技术相比,本发明具有以下有益效果:
1、本发明专利建立尾流跟踪模块,可根据激光雷达测风系统所测实时原始数据完成三维风场反演,获得环境风况参数和各风力机组前侧尾流特征参数,从而完成实时尾流跟踪,可为后续的尾流智能化控制提供可靠的实测输入。
2、本发明专利采用多目标优化方法,综合考虑海上风电场系统总输出功率和机组关键部位所受疲劳载荷情况建立多目标优化函数,基于海上风电场气动-液压-伺服-弹性动力学仿真模型完成流场分布预测,并利用智能优化算法完成在不同环境风况下对最优尾流中心位置的求解,建立多维智能决策数据库LUT。该智能数据库可用于海上风电场的在线搜索和快速优化,实现风电场系统在变化环境风况下的尾流智能化控制。
3、本发明专利建立场级协同PI控制器,以激光雷达测风系统获得的实测尾流中心位置信息作为输入,以多维智能决策数据库LUT搜索获得的最优尾流中心位置信息作为目标值,求解得到单位控制周期内的协同偏航控制调整值,下发给各台机组实现场级协同控制和尾流重定向。本发明所述控制器具有控制效果好、执行速度快等特点。
附图说明
图1为本发明一种基于尾流跟踪的海上风电场场级协同控制策略的系统结构图;
图2为海上风电场场级协同PI控制器的结构图;
图3为PI控制器的结构图;
图4为多维智能决策数据库的建立流程图;
图5为本发明实施例中激光雷达测风系统的实测风况信息;
图6是本发明实施例中海上风电场气动-液压-伺服-弹性动力学仿真模型的WT2、WT3、WT5和WT6机组在优化控制前后其机组前侧所受尾流的中心位置变化情况;
图7是本发明实施例中海上风电场系统在优化控制前后运行性能时域变化情况。
具体实施方式
下面结合附图和实施例对本发明做进一步详细描述,需要指出的是,以下所述实施例旨在便于对本发明的理解,而对其不起任何限定作用。
如图1所示,基于尾流跟踪的海上风电场场级协同控制策略的控制系统,包括机舱式的激光雷达测风系统、尾流跟踪模块、优化器以及场级控制器。利用激光雷达测风系统完成对环境风况和尾流特征参数的辨识,通过优化器中建立的多维智能决策数据库获得最优尾流中心位置信息参考值,利用PI协同控制器通过主动偏航完成尾流重定向与智能化控制。
具体步骤如下:
步骤1,建立尾流跟踪模块,利用海上风电场机舱式激光雷达测风系统测量获得各台海上风力机前侧流场风速分布情况,根据测量的风速分布数据通过三维风场反演完成来流风在无尾流影响下的风信息评估和尾流影响下的尾流特征辨识,获得尾流参数信息和尾流中心位置,完成尾流实时跟踪。
步骤2,根据实际海上风电场的排布特征和流场分布情况,建立海上风电场气动-液压-伺服-弹性动力学仿真模型,具体包括风电场气动仿真模型、海上风力机模型和尾流模型等,并利用上述激光雷达所测历史风速信息对模型参数进行校正,并将该模型用于预测风电场在特定环境输入风况和控制策略下的运行状况。
步骤3,海上风电场基于尾流重定向的多维智能决策数据库LUT用于实现在不同环境风况下的最优控制策略在线实时搜索。建立综合考虑风电场功率最大化和载荷最小化的多目标优化函数,通过智能优化算法得到风电场在不同环境风况输入下的最优尾流中心位置信息,并将上述参数信息记录为多维智能决策数据库。
步骤4,海上风电场场级协同控制器,利用查表法在线搜索多维智能决策数据库LUT获得风电场在相应风况输入下的最优尾流中心位置信息,将其作为场级控制器的参考值,根据尾流跟踪模块测得的尾流信息测量值,通过场级协同PI控制器改变各机组与环境来流风的偏航不对中程度,实现尾流重定向和控制优化。
如图2所示,海上风电场场级协同PI控制器具体包括低通滤波器、多维智能决策数据库LUT、尾流中心位置误差计算模块、PI控制器、偏航执行器和基于测风激光雷达系统的尾流跟踪模块等。在单位控制周期内,通过激光雷达测风系统测量风况数据完成三维风场反演,获得环境输入风的平均风速和风向信息,经过低通滤波后输入多维智能决策数据库LUT在线搜索获得相应工况下的最优尾流中心位置信息,利用PI控制器通过偏航执行器完成尾流重定向,使得实测尾流中心位置与最优尾流中心位置信息一致,完成尾流智能化控制。
如图3所示,PI控制器利用上述激光雷达测风系统和多维智能决策数据库LUT求解得到实测尾流中心位置和最优尾流中心位置的偏差值,通过比例、积分运算得到场级协同偏航控制的调整量,完成尾流重定向。
如图4所示,基于海上风电场气动-液压-伺服-弹性动力学仿真模型预测实际海上风电场在相应环境工况下的动态运行特征,获得不同环境工况参数、机组控制参数与风电场输出功率和机组负载之间的映射关系模型。建立综合考虑降低功率损耗和疲劳载荷的多目标优化函数,利用智能优化算法完成对最优偏航控制策略的寻优,记录相应控制下的最优尾流中心位置信息,建立并保存为多维智能决策数据库LUT。
在具体的实施例中,上述海上风力机模型包括机械模型和控制模型等,可通过修改模型参数来模拟不同型号的海上风电机组。本实施例选择NREL 5MW风电机组进行仿真模拟。
上述流场模型和尾流模型根据海上风电场环境输入风况信息和尾流传播特征完成各机组所受风况的预测,其具体的模型参数由测风激光雷达系统获得的实测数据完成参数校正。本实施例选择实际机舱式激光雷达所测风况信息,该实测原始数据为2019年5月19日所测某风电场内一机组搭载的机舱式激光雷达测风系统在350分钟内的风况数据。基于实测数据,完成三维风场反演,获得不同激光束所测视线风速,并对机组前侧风速时序列进行求解,获得风电场环境风输入的平均风速和风向信息,风信息时序列如图5所示。由上述数据处理后可得,该测量时间内风电场的输入环境风平均风速为8.3892m/s,平均风向为268.6361°(选择正北方向为0°)。
在具体的实施例中,上述海上风电场仿真环境的建立是基于实际海上风电场的特征,具体包括风电机组的排布和环境风资源分布情况等,使得仿真结果与实际工况接近。基于上述输入风况信息,将测风激光雷达所测风信息作为风电场仿真模型的输入工况。建立2×3海上风电场气动-液压-伺服-弹性动力学仿真模型,各机组模型选择NREL 5MW风电机组,机组的转子直径为126m。通过仿真预测该风电场在800s内的运行特性,具体计算其输出功率时序变化和关键部位疲劳载荷时序变化。其中,在各机组轮毂高度处建立X-Y平面坐标系,则海上风电场系统内各机组排布位置信息如下表1所示,其中270°输入风向为OX轴方向(东西方向),Y轴方向为南北方向。
表1海上风电场仿真模型各机组排布位置
机组序号 坐标信息(m)
WT1 (0,-189)
WT2 (630,-189)
WT3 (1260,-189)
WT4 (0,189)
WT5 (630,189)
WT6 (1260,189)
根据获得运行特性的时域变化规律,求解单位控制周期的800s内风电场系统的总输出功率和关键部位疲劳载荷情况。建立综合考虑增效减载的多目标优化函数,如下:
Figure BDA0003854289950000131
式中,
Figure BDA0003854289950000132
为归一化后的风电场系统总发电功率,DELnorm为归一化后的各台风力发电机组关键部位所受最大疲劳载荷,α为权重系数。本实施例选择优化权重系数α为0.6。利用智能优化算法完成对最优偏航控制角的寻优,使得最大化该目标函数的结果。具体地,本实施例选择遗传算法完成最优解的寻优。观察海上风电场系统在相应的环境输入风况和主动偏航协同控制下,各机组前侧所受上游尾流影响情况,记录各机组所受尾流中心位置信息,将其保存为多维智能决策数据库LUT以实现在线实时搜索与优化。
具体地,本实施例在激光雷达测风系统所测风况下进行仿真模拟和多目标优化控制,在线搜索多维智能决策数据库LUT对应工况下的最优尾流中心位置信息,通过场级PI协同控制器完成尾流重定向和智能化控制。在闭环控制下,将在线搜索得到的最优尾流中心位置作为系统的参考值,三维风场反演得到的实测尾流中心位置作为观测值,风电场系统内各机组通过主动协同偏航控制实时跟踪最优尾流中心位置信息,实现尾流的智能化控制。
如图6所示,海上风电场气动-液压-伺服-弹性动力学仿真模型的WT2、WT3、WT5和WT6机组在优化控制前后其机组前侧所受尾流的中心位置变化情况,其纵坐标为各机组前侧所受尾流的中心位置在XY坐标系中Y轴的坐标值变化。各机组前侧的尾流中心位置在上游风力机转子下游处呈现动态蜿蜒变化状态,且机组下游距离越远,尾流蜿蜒现象越显著。
图7为海上风电场系统在优化控制前后运行性能时域变化情况,具体包括总输出功率、塔筒俯仰弯矩时域变化、叶片俯仰弯矩时域变化情况等。由于尾流延迟传播作用,海上风电场系统约在450s后达到稳定状态,因此本实施例主要关注450s-800s时间内的结果。通过计算可得,优化控制后风电场系统的总输出功率可提升10.33%,塔筒疲劳载荷可降低9.74%,而叶片疲劳载荷略有上升,总体而言,本发明所述优化控制策略可显著实现增效降载的多目标优化。
以上所述的实施例对本发明的技术方案和有益效果进行了详细说明,应理解的是以上所述仅为本发明的具体实施例,并不用于限制本发明,凡在本发明的原则范围内所做的任何修改、补充和等同替换,均应包含在本发明的保护范围之内。

Claims (8)

1.一种基于尾流跟踪的海上风电场场级协同控制策略,其特征在于,包括:
(1)利用机舱式激光雷达测风系统测量得到每台海上风力机前侧流场风速分布数据;
(2)建立尾流跟踪模块,利用步骤(1)得到的风速分布数据通过三维风场反演完成来流风在无尾流影响下的风信息评估和尾流影响下的尾流特征参数辨识,尾流特征参数包括尾流深度、尾流宽度和尾流中心位置;
(3)根据实际海上风电场的排布特征和流场分布情况,建立海上风电场气动-液压-伺服-弹性动力学仿真模型;利用激光雷达测风系统所测的历史风速分布数据对仿真模型的参数进行校正,并将该仿真模型用于预测风电场在输入风况和控制策略下的运行状况;
(4)建立海上风电场基于尾流重定向的多维智能决策数据库,用于实现在不同环境风况下的最优控制策略在线实时搜索;
建立综合考虑风电场功率最大化和载荷最小化的多目标优化函数,通过智能优化算法得到风电场在不同环境风况输入下的最优尾流中心位置信息,并将上述参数信息记录到多维智能决策数据库;
(5)建立海上风电场场级协同PI控制器,利用查表法在线搜索多维智能决策数据库获得风电场在相应风况输入下的最优尾流中心位置信息,将其作为场级协同PI控制器的参考值,根据尾流跟踪模块测得的尾流信息测量值,通过场级协同PI控制器改变各机组与环境来流风的偏航不对中程度,实现尾流重定向和控制优化。
2.根据权利要求1所述的基于尾流跟踪的海上风电场场级协同控制策略,其特征在于,步骤(2)的具体过程为:
(2-1)三维风场反演
为了准确评估海上风电场内部风速分布信息,首先进行风向对北校正,通过坐标旋转转换矩阵算法来反演激光束在地球坐标系下的指向;激光束指向的坐标变换公式为:
LLOS=(H1H2H3)-1LLOS_Lidar (1)
其中,LLOS为激光束在地球坐标下的指向,H1、H2、H3为坐标旋转转换矩阵,LLOS_Lidar为激光束在机舱坐标下的指向;
定义机舱坐标系的方向:X0为转子轴向指向机舱头部方向,Y0为转子径向指向水平左侧,Z0与风力机塔筒轴同轴并指向机舱顶部方向;激光束指向通过方位角
Figure FDA0003854289940000021
和俯仰角θ0描述,其中方位角
Figure FDA0003854289940000022
为激光束在X0Y0面的投影与X0的夹角,俯仰角θ0为激光束与X0Y0面的夹角;风力机机舱在协同控制过程中形成姿态的变化,用翻滚角
Figure FDA0003854289940000023
俯仰角θ和偏航角ψ表示;因此,激光束在机舱坐标下的指向LLOS_Lidar及坐标旋转转换矩阵表示为如下:
Figure FDA0003854289940000024
Figure FDA0003854289940000025
Figure FDA0003854289940000026
Figure FDA0003854289940000027
进一步求解得到激光束视线风速vLOS
vLOS=vwind·LLOS (6)
在实际海上机舱式激光雷达扫描测量时,在相同高度和圆周扫描层,视线风速在理想情况下为三角函数,定义为:
f(vLOS)=acos(θ-b)+c (7)
参数a、b和c由所测量高度上的风矢量信息决定,具体地:
Figure FDA0003854289940000031
其中,v1、v2和v3为测量风电场在地球坐标系下的南北方向、东西方向和垂直方向上的分速度,其决定了上述视线风速函数的参数大小;
通过最小二乘法拟合(7)式和(8)式所述的视线风速,通过最小化函数误差(9)式,求解得到曲线函数的参数a、b和c,完成三维风电场风速分布的反演;
Figure FDA0003854289940000032
其中,N为机舱式激光雷达在某个测量周期内的测试点个数;基于上述风场反演,获得风电场来流风在无尾流影响下的风信息,包括风速、风向和湍流强度,同时为尾流特征参数辨识提供数据基础;
(2-2)尾流特征参数辨识
建立风力机前侧所受尾流速度分布的反高斯函数,公式如下:
Figure FDA0003854289940000033
其中,U为风力机轮毂高度的(x,zh)位置处在时间t时刻的风速,U0为时间t时刻的环境平均风速,Cs为尾流比例系数,用以量化尾流深度,e为自然常数,μ为轮毂高度处尾流中心位置,σ为尾流标准差,用以量化尾流宽度,π为圆周率;
利用激光雷达测得的风力机前侧相应距离处的风速数据,利用最小二乘法完成反高斯函数的参数拟合,获得尾流特征参数,具体包括尾流深度、尾流宽度和尾流中心位置信息。
3.根据权利要求1所述的基于尾流跟踪的海上风电场场级协同控制策略,其特征在于,步骤(3)中,海上风电场气动-液压-伺服-弹性动力学仿真模型具体包括风电场气动仿真模型、海上风力机模型和尾流模型。
4.根据权利要求3所述的基于尾流跟踪的海上风电场场级协同控制策略,其特征在于,所述的风电场气动仿真模型利用Taylor湍流冻结假设完成风电场三维流场在空间上的动态传播模拟;所述的海上风力机模型采用制动盘模型,根据实际机组的性能参数完成建模并获得在特定环境工况下机组功率、载荷动态响应;所述的尾流模型用于预测与风电场输出功率和机组负载相关的关键尾流特征,包括尾流风速赤字、尾流膨胀以及尾流蜿蜒。
5.根据权利要求1所述的基于尾流跟踪的海上风电场场级协同控制策略,其特征在于,步骤(4)中,建立综合考虑风电场功率最大化和载荷最小化的多目标优化函数具体过程为:
(4-1)计算海上风力发电机组i的输出功率:
Figure FDA0003854289940000041
式中,Pi为风电场内机组i发电机的输出功率,Tq,i(t)为机组i发电机在t时刻的瞬时扭矩,ωi为机组i发电机在t时刻的瞬时转速,tk为当前时刻,ΔT为当前求解时域周期;
根据上述单台机组的功率计算公式,求解风电场系统的总发电功率即为风电场系统内所有风力发电机组的发电功率之和
Figure FDA0003854289940000043
为:
Figure FDA0003854289940000042
(4-2)计算风力发电机组i所受等效疲劳载荷:
利用风电场气动-液压-伺服-弹性动力学仿真模型获得风电场内风力发电机组i关键部位所受载荷时域变化情况,包括塔筒底部和叶片根部的载荷变化情况;通过雨流计数法,对风力发电机组i所受时域变化的载荷进行等效分析,获得关键部位等效疲劳载荷DEL:
Figure FDA0003854289940000051
式中,DEL为等效疲劳载荷,
Figure FDA0003854289940000052
为时间序列j内的等效循环次数,Nj,i为第i个工况在时间序列j内发生的次数,Lj,i为第i个工况在时间序列j内的载荷范围,m为材料S-N曲线的斜率;
(4-3)为实现风电场总输出功率最大化,并降低机组关键部位所受疲劳载荷,建立多目标优化函数如下:
Figure FDA0003854289940000053
式中,
Figure FDA0003854289940000054
为归一化后的风电场系统总发电功率,DELnorm为归一化后的各台风力发电机组关键部位所受最大疲劳载荷,α为权重系数。
6.根据权利要求5所述的基于尾流跟踪的海上风电场场级协同控制策略,其特征在于,建立多目标优化函数后,通过智能优化算法求解在特定环境工况下最大化目标函数J的最优偏航控制角,并通过风电场气动-液压-伺服-弹性动力学仿真模型记录各机组在该偏航控制输入下转子前侧的最优尾流中心位置;最后,将上述环境风况参数和对应的最优尾流中心位置信息保存为多维智能决策数据库,用以不同环境风况下各机组前侧最优尾流中心位置的在线实时搜索和场级协同控制。
7.根据权利要求1所述的基于尾流跟踪的海上风电场场级协同控制策略,其特征在于,步骤(5)中,所述的海上风电场场级协同PI控制器包含低通滤波器、多维智能决策数据库、尾流中心位置误差计算模块、PI控制器、偏航执行器和激光雷达测风系统;
将激光雷达测风系统获得的环境风况测量值作为输入,通过低通滤波器获得外界环境风输入的平均风速与风向,在低通滤波器后接多维智能决策数据库;利用查表法实现在线搜索,获得对应环境风况下的各机组前侧最优尾流中心位置信息;利用尾流中心位置误差计算模块完成各机组前侧实际尾流中心位置与期望尾流中心位置的误差估计;将尾流中心位置误差计算模块获得的中心位置偏差值输入PI控制器,对偏差值进行比例、积分运算,得到风电场在单位控制周期内进行主动协同偏航控制的调整值;通过偏航执行器改变各机组与环境来流风的偏航不对中程度,实现尾流重定向和控制优化。
8.根据权利要求7所述的基于尾流跟踪的海上风电场场级协同控制策略,其特征在于,利用尾流中心位置误差计算模块完成各机组前侧实际尾流中心位置与期望尾流中心位置的误差估计具体包括:
通过激光雷达测风系统的三维风场反演获得实测尾流中心位置信息参数,通过在线搜索多维智能决策数据库获取最优尾流中心位置信息参数,通过尾流中心位置误差计算模块对两者做差获得中心位置的偏差值。
CN202211143225.3A 2022-09-20 2022-09-20 一种基于尾流跟踪的海上风电场场级协同控制策略 Pending CN115807734A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211143225.3A CN115807734A (zh) 2022-09-20 2022-09-20 一种基于尾流跟踪的海上风电场场级协同控制策略

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211143225.3A CN115807734A (zh) 2022-09-20 2022-09-20 一种基于尾流跟踪的海上风电场场级协同控制策略

Publications (1)

Publication Number Publication Date
CN115807734A true CN115807734A (zh) 2023-03-17

Family

ID=85482668

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211143225.3A Pending CN115807734A (zh) 2022-09-20 2022-09-20 一种基于尾流跟踪的海上风电场场级协同控制策略

Country Status (1)

Country Link
CN (1) CN115807734A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116771596A (zh) * 2023-06-30 2023-09-19 渤海石油航务建筑工程有限责任公司 海上风电场尾流转向控制方法及相关设备
CN117231425A (zh) * 2023-11-13 2023-12-15 中国海洋大学 基于平台前倾的漂浮式风机的尾流控制系统及方法
CN117872725A (zh) * 2024-03-11 2024-04-12 天津瀚福精密液压技术有限公司 一种液压数据控制方法

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116771596A (zh) * 2023-06-30 2023-09-19 渤海石油航务建筑工程有限责任公司 海上风电场尾流转向控制方法及相关设备
CN116771596B (zh) * 2023-06-30 2024-06-04 渤海石油航务建筑工程有限责任公司 海上风电场尾流转向控制方法及相关设备
CN117231425A (zh) * 2023-11-13 2023-12-15 中国海洋大学 基于平台前倾的漂浮式风机的尾流控制系统及方法
CN117231425B (zh) * 2023-11-13 2024-02-06 中国海洋大学 基于平台前倾的漂浮式风机的尾流控制系统及方法
CN117872725A (zh) * 2024-03-11 2024-04-12 天津瀚福精密液压技术有限公司 一种液压数据控制方法
CN117872725B (zh) * 2024-03-11 2024-05-14 天津瀚福精密液压技术有限公司 一种液压数据控制方法

Similar Documents

Publication Publication Date Title
CN108334672B (zh) 基于有效风速估计的变速风力发电机组最大风能捕获方法
CN115807734A (zh) 一种基于尾流跟踪的海上风电场场级协同控制策略
CN107559143B (zh) 一种大型风力机尾缘襟翼结构参数寻优及多目标襟翼优化控制方法
Li et al. Using neural networks to estimate wind turbine power generation
CN111881572B (zh) 一种基于多目标优化的风电场协同偏航智能控制方法
CN103850876B (zh) 一种适用于无载荷测量的风电机组独立变桨控制方法
CN103061980A (zh) 基于激光测风雷达的风力发电机组的前馈控制系统及其控制方法
CN108953060A (zh) 基于激光雷达测风仪的风电场场级偏航控制方法
CN112260325A (zh) 一种考虑风力发电机疲劳载荷的风电场有功功率控制方法
CN111476293A (zh) 一种计及风速风向的风电场分群方法
US20240309843A1 (en) Cooperative operation optimization control method for wind turbine groups
CN111342499B (zh) 一种基于风功率预测数据的风电场实时调度方法
CN111502914B (zh) 基于线性变参数系统的风机变桨控制器设计方法
CN113011994A (zh) 风电厂的优化
CN110889781B (zh) 基于滑模控制的风电机组保性能最大功率跟踪方法
CN115578016A (zh) 一种模型不完备风电场调频能力在线评估方法
CN114186407A (zh) 一种参数可自适应调整的风电场尾流速度场计算方法及系统
CN108223274B (zh) 基于优化rbf神经网络的大型风力机变桨距系统辨识方法
CN109256814A (zh) 一种风电机组功率曲线的智能修正方法
CN117610449A (zh) 一种基于GRA-XGBoost的风电机组入流风速估计方法
Chhor et al. Operation and control strategies for wind energy conversion systems: Review and simulation study
Xingjia et al. Individual pitch control for variable speed turbine blade load mitigation
CN116658359A (zh) 一种风电场实时协同偏航控制方法
CN116961021A (zh) 考虑风速空间分布和波动性的风电场调频能量评估方法
CN109753759A (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