CN110345006B - 一种风力发电机组低风速区最大功率跟踪优化控制方法 - Google Patents
一种风力发电机组低风速区最大功率跟踪优化控制方法 Download PDFInfo
- Publication number
- CN110345006B CN110345006B CN201910259041.5A CN201910259041A CN110345006B CN 110345006 B CN110345006 B CN 110345006B CN 201910259041 A CN201910259041 A CN 201910259041A CN 110345006 B CN110345006 B CN 110345006B
- Authority
- CN
- China
- Prior art keywords
- wind
- speed
- wind speed
- control
- algorithm
- 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
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000005457 optimization Methods 0.000 title claims abstract description 13
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 49
- 239000002245 particle Substances 0.000 claims abstract description 29
- 238000001914 filtration Methods 0.000 claims abstract description 15
- 238000005259 measurement Methods 0.000 claims abstract description 15
- 230000010355 oscillation Effects 0.000 claims abstract description 5
- 230000000694 effects Effects 0.000 claims description 12
- 230000005540 biological transmission Effects 0.000 claims description 8
- 230000006870 function Effects 0.000 claims description 8
- 230000003044 adaptive effect Effects 0.000 claims description 7
- 230000002349 favourable effect Effects 0.000 claims description 5
- 238000004088 simulation Methods 0.000 claims description 5
- 238000013016 damping Methods 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 4
- 230000004044 response Effects 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 230000013016 learning Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 238000013459 approach Methods 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 claims description 2
- 230000001149 cognitive effect Effects 0.000 claims description 2
- 238000012937 correction Methods 0.000 claims description 2
- 230000003247 decreasing effect Effects 0.000 claims description 2
- 230000009191 jumping Effects 0.000 claims description 2
- 230000007246 mechanism Effects 0.000 claims description 2
- 230000009326 social learning Effects 0.000 claims description 2
- 230000007704 transition Effects 0.000 claims description 2
- 238000013461 design Methods 0.000 abstract description 5
- 230000007547 defect Effects 0.000 abstract description 4
- 238000012545 processing Methods 0.000 abstract description 3
- 238000011217 control strategy Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000002028 premature Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000005251 gamma ray Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
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/04—Automatic control; Regulation
- F03D7/042—Automatic control; Regulation by means of an electrical or electronic controller
- F03D7/043—Automatic control; Regulation by means of an electrical or electronic controller characterised by the type of control logic
-
- 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/043—Automatic control; Regulation by means of an electrical or electronic controller characterised by the type of control logic
- F03D7/045—Automatic control; Regulation by means of an electrical or electronic controller characterised by the type of control logic with model-based controls
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
- F05B2260/00—Function
- F05B2260/96—Preventing, counteracting or reducing vibration or noise
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
- F05B2270/00—Control
- F05B2270/30—Control parameters, e.g. input parameters
- F05B2270/32—Wind speeds
-
- 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)
- Control Of Eletrric Generators (AREA)
Abstract
本发明提出了一种风力发电机组低风速区最大功率跟踪优化控制方法,主要包括风速信号的采集与滤波处理、改进控制器的设计、控制器参数的优化等三个部分。信号采集部分采用脉冲激光雷达测量风轮不同高度不同位置处的风速,并采用自适应卡尔曼滤波算法对风速信号进行滤波处理;改进控制器设计部分在传统最佳转矩控制的基础上引入超扭曲高阶滑模算法,以便提高转速跟踪性能;控制器参数寻优部分采用改进的粒子群算法,克服了常规粒子群算法容易早熟收敛以及算法后期易在全局最优解附近产生振荡现象的不足。本发明不但可以实现对实时变化风速的有效测量,而且可以提高风轮转速跟踪实时风速的性能,从而提高风能利用率,具有一定的实用价值。
Description
技术领域
本发明涉及风力发电机组的控制领域,特别涉及一种基于超扭曲优化算法的风机最大功率跟踪控制方法。
背景技术
众所周知,当今社会能源和环境问题日益突出,传统化石能源由于具有污染严重、不可再生等缺点,已经不能满足可持续发展的要求,因此,开发利用新能源显得尤为重要。而风力发电以其能源利用率高、技术相对成熟等优点渐渐成为最具有发展前景的新能源技术之一。
根据风速的高低,通常将风电机组的运行区域划分为四个部分:风速低于切入风速的部分、风速介于切入风速与额定风速之间的区域(低风速区)、风速接近额定风速的区域、风速高于额定风速的区域(高风速区)。其中,低风速区为系统运行的典型工作区之一,其主要控制目标是控制风轮转速跟踪由实时风速所决定的最佳转速,从而捕获更多风能,使风机输出功率最大化。因此,获得较优的控制效果既依赖于对实时风速的准确测量,又依赖于优化控制算法的合理设计。
在风速测量方面,传统的风速仪、风向标等传感器一般位于风机的逆风区,受到塔影效应和尾流效应的影响,与到达风轮处的真实风速有一定的差别。同时,由于湍流风具有随机性、时变性,风速在整个风轮扫掠面上不断变化,而机械式仪表具有转动惯量,无法确保风速测量的实时性、准确性,依据其设计的控制策略很难确保最优。此外,风电场恶劣的自然环境会缩短传感器的使用寿命。因此,现有的研究大多根据风电机组的输出特性对风速进行估计,具有较高的估计精度,但不可否认的是,估计的风速与实际风速相比具有滞后性,根据其采取的控制律不能实时控制风机,导致控制性能大打折扣。
控制策略方面,在工程中应用最为广泛的最佳转矩法具有简单、易于实现的特点,但是当风速较低或风电机组转动惯量较大时,风轮转速跟踪实时风速的响应速度会下降,导风能转换效率的降低。常规的滑模控制在处理非线性、不确定性扰动等方面具有很大的优势,能够使系统在一定特性下沿规定的状态轨迹运动,且结构简单、响应快速,但是会产生抖振。
发明内容
针对上述背景技术中存在的问题,本发明提出一种基于激光雷达测风速的、将超扭曲算法与最佳转矩控制相结合的最大功率跟踪改进控制策略,以5MW风力发电机组为对象,在对激光雷达测量风速进行处理的基础上进行控制策略的实施,同时改进采用粒子群算法在收敛域内依据工程要求对控制器参数进行进一步优化,以提升控制系统的综合性能。并将控制效果与最佳转矩法和常规滑模一最佳转矩控制的效果进行对比,通过仿真实验证明了该方法的有效性。
本发明的技术方案为:一种风力发电机组低风速区最大功率跟踪优化控制方法,实施步骤如下:
1)采用脉冲激光雷达测量风轮不同高度不同位置处的风速。
2)对测量的风速采用自适应卡尔曼滤波算法对所测量的风速信号进行滤波处理。
3)根据脉冲激光雷达所测量和滤波处理后的风速信号,计算风轮转速的参考值。
4)建立风力发电机组的空气动力系统模型、传动系统模型和风力发电机模型。
5)在传统最佳转矩控制法的基础上,引入超扭曲高阶滑模算法,设计风力发电机组的最大功率跟踪改进控制策略。
6)根据工程要求,在超扭曲高阶滑模算法的收敛域内采用改进粒子群算法对控制器参数进行进一步优化。
本发明的有益效果在于:本发明采用激光雷达测量风轮转子前方一定距离的风速信息,并采用自适应卡尔曼滤波算法在考虑测量噪声的情况下重建模型的缺失状态,大大提高了测量的实时性和准确性,根据其计算的风轮转速的参考值更接近实际情况、所设计的控制器更实用。将超扭曲高阶滑模算法与传统的最佳转矩法相结合,可提高风轮转速跟踪实时风速的性能,提高风能利用效率,同时可抑制发电机转矩的抖振。采用改进粒子群算法对控制器参数进行优化,可以取得更好的控制性能,同时克服了常规粒子群算法容易早熟收敛以及算法后期易在全局最优解附近产生振荡现象的不足。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图,其中:
图1为本发明的一较佳实施例的安装在风机上的脉冲激光雷达的示意图;
图2为本发明的一较佳实施例的Cp-λ关系曲线图;
图3为本发明的一较佳实施例的风机控制系统结构图;
图4为本发明的一较佳实施例的风速曲线图;
图5为本发明的一较佳实施例的风轮角速度对比图;
图6为本发明的一较佳实施例的发电机角速度对比图;
图7为本发明的一较佳实施例的发电机输出功率对比图;
图8为本发明的一较佳实施例的风能利用系数(功率系数)对比图;
图9为本发明的一较佳实施例的发电机转矩对比图;
表1为本发明的一较佳实施例的风电机组的主要参数表;
表2为本发明的一较佳实施例的性能对比表。
具体实施方式
下面将对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
请参阅图1至图9,本发明实施例包括:
一种风力发电机组低风速区最大功率跟踪优化控制方法,包括以下步骤:
步骤1:采用脉冲激光雷达测量风轮不同高度不同位置处的风速。
如图1所示,为安装在风机上的脉冲激光雷达的示意图。脉冲激光雷达在风轮转子前发送三个激光脉冲,提供三个75%叶片跨度位置的预览风速测量值。通过脉冲激光雷达的距离门距,可以获得风轮转子前方一定距离的预览风。风速矢量包含三个维度一水平分量u、轴向分量v和垂直分量w,由于水平分量和垂直分量较小,故可忽略不计,激光测量点几乎沿轴向。轴向分量v可近似用下式表示:
其中:θ为激光束的圆锥角,vL为风速。
平均风速vp可以通过激光雷达的最后扫描圆测量,其表达式为:
其中:T为脉冲激光雷达的扫描周期。
步骤2:对测量的风速采用自适应卡尔曼滤波算法对所测量的风速信号进行滤波处理。由于自然风行进到风轮转子时仍然具有一定的时间延迟,同时,测量信号中难免会夹杂一些噪声,从而降低激光雷达测量风速信号的准确性,从而影响后续控制器的设计。
为了解决上述问题,本发明采用Sage-Husa自适应卡尔曼滤波器对风轮转子预览中的风速信号进行滤波处理,利用观测输出对系统噪声和观测噪声的统计特性进行实时动态估计,估计的噪声参数用于状态估算,从而重建模型的缺失状态。状态和估计误差方差预测公式为:
其中:Pk,k-1为估计风速误差的协方差,为零均值、协方差为Qk-1的高斯白噪声序列。
协方差自适应调整机制为:
dk=(1-b)(1-bk+1)
其中:Φk为系统状态转移矩阵,H为观测矩阵,zk为系统的观测值,Vk为系统的观测噪声,0<b<1为遗忘因子。
将估计误差用于补偿激光雷达预览,从而获得更实时、准确的优化风速。
步骤3:根据脉冲激光雷达所测量和滤波处理后的风速信号,计算风轮转速的参考值,计算公式为:
其中:v为激光雷达测量并经滤波处理后的轴向风速,R为风轮半径,λopt为最佳叶尖速比。如图2所示,本发明所选风电机组的最佳叶尖速比λopt=9.031。
步骤4:建立风力发电机组的空气动力系统模型。空气动力系统的核心设备是风轮,当外部风速为v时,风轮捕获的有效功率为:
其中:ρ为空气密度;R为风轮的半径;v为风速,它与叶尖速比λ、桨距角β成非线性关系,如图2所示。
叶尖速比λ的定义为:其中:ωr为风轮角速度。
风轮所获得的气动转矩为:
其中:
步骤5:建立风力发电机组的传动系统模型。传动系统由高速轴和低速轴通过齿轮箱相互连接而组成,在忽略轴系刚度和阻尼系数情况下,可用如下单质量块模型加以描述:
B=Br+ng 2Bg
J=Jr+ng 2Jg
其中,Br和Bg分别为低速轴和高速轴的阻尼系数,Jr和Jg分别为低速轴和高速轴的转动惯量,ng为齿轮箱的传动比,Tg为发电机转矩。
步骤6:建立风力发电机模型。本发明采用双馈异步发电机,为方便起见,可将其电磁部分简化为一个一阶线性模型:
其中:τg中为转矩系数,Tg,ref为发电机转矩的参考值。
发电机输出功率可表示为:
Pe=Tgωg
其中:ωg为发电机角速度。
由于本发明的控制器设计以发电机转矩参考值Tg,ref为主要控制量,故可将发电机电磁动态的影响作为扰动ψ加到风机模型,即
Tg=Tg,ref+ψ
步骤7:定义滑模面s:
要实现风电机组在低风速区的控制目标,必须使风轮角速度ωr(t)实时跟踪参考角速度ωret(t),为了使风轮角速度与参考角速度的误差以及误差的变化率在有限时间内收敛到0,可定义滑模面为:
其中:c>0,e(t)为风轮角速度与参考角速度的误差。
步骤8:控制器设计。采用最佳转矩控制实现对发电机角速度的闭环反馈调节,使发电机角速度自行趋近最佳角速度。为了提高响应速度,考虑引入滑模变结构控制的思想,因高阶滑模控制既保留了常规滑模控制的优点又能有效抑制常规滑模控制普遍具有的抖振现象,其中的超扭曲算法更是得到了广泛应用。为此,本发明在最佳转矩控制的基础上引入超扭曲算法,得到改进最佳转矩控制的控制律:
u=Tg,ref=ueq+usw
其中:ueq为等效控制,对应于最佳转矩控制,usw为切换控制,对应于超扭曲高阶滑模算法,γ1和γ2为控制器参数,0<d<1。当d接近0时,控制器的滑模效果会增强,相当于常规的滑模控制器,虽能提高转速跟踪性能,但抖振比较明显;当d接近1时,控制器接近PI控制器,虽能消除抖振但超调会产生得更加频繁,影响转速跟踪效果。为了在提高转速跟踪性能的同时有效抑制发电机转矩抖振,本发明折衷选取d=1/2。
步骤9:选取Lyapunov函数V(s,y)=ηTMη
其中:
选取合适的参数,使Lyapunov函数的导数小于0,即可保证系统在有限时间内收敛到滑模面。从而可得收敛域:
只要控制器参数γ1和γ2处于上述范围,即可保证整个风机控制系统的稳定性和基本的控制性能。
步骤10:采用改进粒子群算法对控制器参数进行迭代寻优。为了在实际应用中获得最优的控制性能,本发明拟采用改进粒子群算法在上式给出的范围内依据工程要求对参数γ1和γ2进行进一步优化。定义如下目标函数作为粒子群算法的适应度函数:
其中:t1为完成单次仿真的时间,Pe为发电机输出功率,e(t)为风轮角速度跟踪误差,std(Tg)为发电机转矩的标准差,w1~w3为加权系数,可根据实际情况进行相应调整。
步骤11:在基本粒子群优化算法的基础上,引入惯性权重因子w,得到修正之后的粒子速度更新公式:
vid(t+1)=wvid(t)+c1r1(pid(t)-xid(t))+c2r2(Pgd(t)-xid(t))
粒子位置更新方程为:
xid(t+1)=xid(t)+vid(t+1)
其中:vid(t)、xid(t)为t时刻粒子i第d维的速度和位置;vid(t+1)、xid(t+1)为t+1时刻粒子i第d维的速度和位置;w称为惯性权重因子。c1和c2分别称为认知学习因子和社会学习因子,一般取值为0.5;r1和r2为相互独立的随机数,它们在[0,1]区间内服从均匀分布。
步骤12:由于较大的权重因子有利于跳出局部最小点,便于全局搜索,而较小的惯性因子则有利于对当前的搜索区域进行精确局部搜索,以利于算法收敛。因此,针对常规粒子群算法容易早熟收敛以及算法后期易在全局最优解附近产生振荡的现象,本发明采用线性变化的权重,让惯性权重从最大值ωmax线性减小到最小值ωmin。随算法迭代次数的公式为:
其中:ωmax,ωmin分别为ω的最大值和最小值,t为当前迭代步数,tmax为最大迭代步数,通常取ωmax=0.9,ωmin=0.4。
采用美国国家能源部可再生能源实验室开发的TURBSIM软件生成平均风速为7m/s、湍流强度为A类级别、时长为400s的风速序列,采用FAST软件模拟激光雷达,由图4不难看出,采用激光雷达测量的风速与实际风速基本吻合,精度较高,风电机组的基本参数如表1所示。在上述实验条件下进行仿真,图4-图9为仿真结果图,表2为性能对比表。将本发明所提方法(STC-OT)与传统最佳转矩控制(OT)和常规滑模-最佳转矩控制(SMC-OT)的效果进行对比,如图5-图9以及表2所示。本发明选取滑模面参数c=50,选择粒子群种群数为50,学习因子c1=c2=0.5,粒子群算法经过30余次循环的迭代使得寻优参数趋于稳定,最终得到最佳的控制器参数值为:γ1=-3.875e4;γ2=-605.335。
表1 风电机组的主要参数
表2 性能对比
图5和图6分别为风轮角速度曲线和发电机角速度曲线,可见,三种控制策略均能使风轮角速度和发电机角速度跟踪其参考曲线。OT、SMC-OT、STC-OT的风轮角速度平均跟踪误差(AVG(e))的数值分别为0.0125、0.0067、-0.0024,说明STC-OT的转速跟踪效果更好。
发电机输出功率曲线和风能利用系数曲线如图7和图8所示,性能对比数据如表2所示,OT、SMC-OT、STC-OT的发电机平均输出功率(AVG(Pe))分别为1.9491e6、1.9862e6、2.0349e6,可见,采用STC-OT时发电机平均功率比采用OT时提高了4.40%,比采用SMC-OT时提高了2.45%,平均风能利用系数比采用OT时提高了4.39%,比采用SMC-OT时提高了2.12%,说明常规的滑模控制和超扭曲控制均能提高风能利用效率,但超扭曲控制更胜一筹。
图9为发电机转矩曲线,可见,相对于SMC-OT而言,STC-OT使发电机转矩的标准差(STD(Tg)减小了10.8%,说明超扭曲滑模控制可有效抑制常规滑模控制普遍具有的抖振现象。
本发明风力发电机组低风速区最大功率跟踪优化控制方法的有益效果是:
1、本发明采用激光雷达测量风轮转子前方一定距离的风速信息,并采用自适应卡尔曼滤波算法在考虑测量噪声的情况下重建模型的缺失状态,大大提高了测量的实时性和准确性,根据其计算的风轮转速的参考值更接近实际情况、所设计的控制器更实用。
2、本发明将超扭曲高阶滑模算法与传统的最佳转矩法相结合,可提高风轮转速跟踪实时风速的性能,提高风能利用效率,同时可抑制发电机转矩的抖振。
3、本发明采用改进粒子群算法对控制器参数进行优化,可以取得更好的控制性能,同时克服了常规粒子群算法容易早熟收敛以及算法后期易在全局最优解附近产生振荡现象的不足。
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书内容所作的等效结构或等效流程变换,或直接或间接运用在其它相关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (9)
1.一种风力发电机组低风速区最大功率跟踪优化控制方法,其特征在于,包括以下步骤:(1)采用脉冲激光雷达测量风轮不同高度不同位置处的风速;
(2)对测量的风速采用自适应卡尔曼滤波算法对所测量的风速信号进行滤波处理;
(3)根据脉冲激光雷达所测量和滤波处理后的风速信号,计算风轮转速的参考值;
(4)建立风力发电机组的空气动力系统模型;
(5)建立风力发电机组的传动系统模型;
(6)建立风力发电机模型;
(7)定义滑模面s;
(8)控制器设计;
(9)采用改进粒子群算法对控制器参数进行迭代寻优;
所述脉冲激光雷达测量风速:脉冲激光雷达在风轮转子前发送三个激光脉冲,提供三个75%叶片跨度位置的预览风速测量值,通过脉冲激光雷达的距离门距,获得风轮转子前方一定距离的预览风,风速矢量包含三个维度-水平分量u、轴向分量v和垂直分量w,由于水平分量和垂直分量较小,故可忽略不计,激光测量点几乎沿轴向,轴向分量v可近似用下式表示:
其中:θ为激光束的圆锥角,vL为风速;
平均风速vp可以通过激光雷达的最后扫描圆测量,其表达式为:
其中:T为脉冲激光雷达的扫描周期。
2.根据权利要求1所述的风力发电机组低风速区最大功率跟踪优化控制方法,其特征在于,所述风速信号的滤波处理,采用Sage-Husa自适应卡尔曼滤波器对风轮转子预览中的风速信号进行滤波处理,利用观测输出对系统噪声和观测噪声的统计特性进行实时动态估计,估计的噪声参数用于状态估算,从而重建模型的缺失状态,状态和估计误差方差预测公式为:
其中:Pk,k-1为估计风速误差的协方差,为零均值、协方差为Qk-1的高斯白噪声序列,协方差自适应调整机制为:
dk=(l-b)(1-bk+1)
其中:Φk为系统状态转移矩阵,H为观测矩阵,zk为系统的观测值,Vk为系统的观测噪声,0<b<1为遗忘因子,
将估计误差用于补偿激光雷达预览,从而获得更实时、准确的优化风速。
3.根据权利要求1所述的风力发电机组低风速区最大功率跟踪优化控制方法,其特征在于,所述风轮转速的参考值的计算,计算公式为:
其中:v为激光雷达测量并经滤波处理后的轴向风速,R为风轮半径,λopt为最佳叶尖速比。
4.根据权利要求1所述的风力发电机组低风速区最大功率跟踪优化控制方法,其特征在于,所述空气动力系统模型:空气动力系统的核心设备是风轮,当外部风速为v时,风轮捕获的有效功率为:
其中:ρ为空气密度;R为风轮的半径;v为风速,Cp与叶尖速比λ、桨距角β成非线性关系;
叶尖速比λ的定义为:
其中:ωr为风轮角速度;
风轮所获得的气动转矩为:
其中:
5.根据权利要求1所述的风力发电机组低风速区最大功率跟踪优化控制方法,其特征在于,所述风力发电机组的传动系统模型:传动系统由高速轴和低速轴通过齿轮箱相互连接而组成,在忽略轴系刚度和阻尼系数情况下,可用如下单质量块模型加以描述:
B=Br+ng 2Bg
J=Jr+ng 2Jg
其中,Br和Bg分别为低速轴和高速轴的阻尼系数,Jr和Jg分别为低速轴和高速轴的转动惯量,ng为齿轮箱的传动比,Tg为发电机转矩。
6.根据权利要求1所述的风力发电机组低风速区最大功率跟踪优化控制方法,其特征在于,所述风力发电机模型:采用双馈异步发电机,为方便起见,可将其电磁部分简化为一个一阶线性模型:
其中:τg中为转矩系数,Tg,ref为发电机转矩的参考值;
发电机输出功率可表示为:
Pe=Tgωg
其中:ωg为发电机角速度;
由于控制器设计以发电机转矩参考值Tg,ref为主要控制量,故可将发电机电磁动态的影响作为扰动ψ加到风机模型,即Tg=Tg,ref+ψ。
7.根据权利要求1所述的风力发电机组低风速区最大功率跟踪优化控制方法,其特征在于,所述滑模面的定义:要实现风电机组在低风速区的控制目标,必须使风轮角速度ωr(t)实时跟踪参考角速度ωref(t),为了使风轮角速度与参考角速度的误差以及误差的变化率在有限时间内收敛到0,可定义滑模面为:
其中:c>0,e(t)为风轮角速度与参考角速度的误差。
8.根据权利要求1所述的风力发电机组低风速区最大功率跟踪优化控制方法,其特征在于,所述控制器的设计:采用最佳转矩控制实现对发电机角速度的闭环反馈调节,使发电机角速度自行趋近最佳角速度,为了提高响应速度,考虑引入滑模变结构控制的思想,因高阶滑模控制既保留了常规滑模控制的优点又能有效抑制常规滑模控制普遍具有的抖振现象,其中的超扭曲算法更是得到了广泛应用,为此,在最佳转矩控制的基础上引入超扭曲算法,得到改进最佳转矩控制的控制律:
u=Tg,ref=ueq+usw
其中:ueq为等效控制,对应于最佳转矩控制,usw为切换控制,对应于超扭曲高阶滑模算法,γ1和γ2为控制器参数,0<d<1;当d接近0时,控制器的滑模效果会增强,相当于常规的滑模控制器,虽能提高转速跟踪性能,但抖振比较明显;当d接近1时,控制器接近PI控制器,虽能消除抖振但超调会产生得更加频繁,影响转速跟踪效果;为了在提高转速跟踪性能的同时有效抑制发电机转矩抖振,选取d=1/2,选取Lyapunov函数V(s,y)=ηTMη,其中:
选取合适的参数,使Lyapunov函数的导数小于0,即可保证系统在有限时间内收敛到滑模面,从而可得收敛域:
只要控制器参数γ1和γ2处于上述范围,即可保证整个风机控制系统的稳定性和基本的控制性能。
9.根据权利要求1所述的风力发电机组低风速区最大功率跟踪优化控制方法,其特征在于,采用改进粒子群算法对控制器参数进行迭代寻优:为了在实际应用中获得最优的控制性能,采用改进粒子群算法在上式给出的范围内依据工程要求对参数γ1和γ2进行进一步优化;定义如下目标函数作为粒子群算法的适应度函数:
其中:t1为完成单次仿真的时间,Pe为发电机输出功率,e(t)为风轮角速度跟踪误差,std(Tg)为发电机转矩的标准差,w1~w3为加权系数,可根据实际情况进行相应调整;
在基本粒子群优化算法的基础上,引入惯性权重因子w,得到修正之后的粒子速度更新公式:
vid(t+1)=wvid(t)+c1r1(pid(t)-xid(t))+c2r2(pgd(t)-xid(t));
粒子位置更新方程为:
xid(t+1)=xid(t)+vid(t+1)
其中:vid(t)、xid(t)为t时刻粒子i第d维的速度和位置;vid(t+1)、xid(t+1)为t+1时刻粒子i第d维的速度和位置;w称为惯性权重因子;c1和c2分别称为认知学习因子和社会学习因子,一般取值为0.5;r1和r2为相互独立的随机数,它们在[0,1]区间内服从均匀分布;
由于较大的权重因子有利于跳出局部最小点,便于全局搜索,而较小的惯性因子则有利于对当前的搜索区域进行精确局部搜索,以利于算法收敛,因此,针对常规粒子群算法容易早熟收敛以及算法后期易在全局最优解附近产生振荡的现象,采用线性变化的权重,让惯性权重从最大值ωmax线性减小到最小值ωmin,随算法迭代次数的公式为:
其中:ωmax,ωmin分别为ω的最大值和最小值,t为当前迭代步数,tmax为最大迭代步数,通常取ωmax=0.9,ωmin=0.4。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910259041.5A CN110345006B (zh) | 2019-03-29 | 2019-03-29 | 一种风力发电机组低风速区最大功率跟踪优化控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910259041.5A CN110345006B (zh) | 2019-03-29 | 2019-03-29 | 一种风力发电机组低风速区最大功率跟踪优化控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110345006A CN110345006A (zh) | 2019-10-18 |
CN110345006B true CN110345006B (zh) | 2023-08-18 |
Family
ID=68174137
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910259041.5A Active CN110345006B (zh) | 2019-03-29 | 2019-03-29 | 一种风力发电机组低风速区最大功率跟踪优化控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110345006B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110940334B (zh) * | 2019-10-23 | 2023-04-18 | 南京华科广发通信科技有限公司 | 一种人体行走测速徽章及测速方法 |
CN110889781B (zh) * | 2019-12-04 | 2022-05-27 | 浙江大学 | 基于滑模控制的风电机组保性能最大功率跟踪方法 |
CN111336062B (zh) * | 2020-03-05 | 2021-11-09 | 中国大唐集团科学技术研究院有限公司华中电力试验研究院 | 一种基于测量风速的风力发电机组最大风能捕获方法 |
CN111541410B (zh) * | 2020-05-12 | 2021-08-24 | 郑州轻工业大学 | 一种开绕组无刷双馈发电机直接功率控制的方法 |
CN111894802B (zh) * | 2020-07-10 | 2021-05-04 | 南京理工大学 | 基于Jensen模型和双光束雷达的风机尾流区偏航误差修正方法 |
CN111859651A (zh) * | 2020-07-10 | 2020-10-30 | 浙江运达风电股份有限公司 | 一种低空气密度下风电机组发电性能优化方法 |
CN112145376B (zh) * | 2020-09-29 | 2021-06-22 | 沈阳航空航天大学 | 一种风力机全时效率测定方法 |
CN112486019B (zh) * | 2020-11-17 | 2022-06-07 | 九江学院 | 不确定风力发电机系统的最大功率跟踪模糊控制方法 |
CN113759708B (zh) * | 2021-02-09 | 2024-07-16 | 京东城市(北京)数字科技有限公司 | 系统的优化控制方法、装置和电子设备 |
CN114215688B (zh) * | 2021-12-15 | 2023-05-02 | 中国华能集团清洁能源技术研究院有限公司 | 风电机组最大功率跟踪控制方法、装置、设备和存储介质 |
CN116517791B (zh) * | 2023-05-30 | 2024-02-02 | 中电投新疆能源化工集团哈密有限公司 | 基于激光测风雷达的多环境适用风速综合监控系统及方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107524563A (zh) * | 2017-08-21 | 2017-12-29 | 华南理工大学 | 一种基于滑模极值搜索的控制方法 |
-
2019
- 2019-03-29 CN CN201910259041.5A patent/CN110345006B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107524563A (zh) * | 2017-08-21 | 2017-12-29 | 华南理工大学 | 一种基于滑模极值搜索的控制方法 |
Non-Patent Citations (1)
Title |
---|
双馈风力发电系统粒子群优化算法研究;侯丽珺;中国优秀硕士学位论文全文数据库工程科技II辑;C042-1401 * |
Also Published As
Publication number | Publication date |
---|---|
CN110345006A (zh) | 2019-10-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110345006B (zh) | 一种风力发电机组低风速区最大功率跟踪优化控制方法 | |
CN108334672B (zh) | 基于有效风速估计的变速风力发电机组最大风能捕获方法 | |
CN104976056B (zh) | 用于风力涡轮的推力速度控制的系统和方法 | |
Quarton | The evolution of wind turbine design analysis—a twenty year progress review | |
CN109989884B (zh) | 一种风机偏航控制方法 | |
CN105179164B (zh) | 基于t-s模糊模型的风能转换系统滑模控制方法及装置 | |
CN104141591B (zh) | 风力发电最大功率点跟踪的改进自适应转矩控制方法 | |
Foti et al. | Wake meandering of a model wind turbine operating in two different regimes | |
CN110454329B (zh) | 一种风电机组桨距角控制方法 | |
CN109416021A (zh) | 风力涡轮机的基于激光雷达的多变量前馈控制 | |
CN102797629A (zh) | 一种风电机组的控制方法、控制器及其控制系统 | |
CN108488035B (zh) | 永磁直驱风力发电机组失速和变桨混合控制方法 | |
KR101734757B1 (ko) | 윈드쉬어에 의한 풍력터빈 블레이드의 추력 모델링 방법 및 시스템 | |
US10451036B2 (en) | Adjustment factor for aerodynamic performance map | |
CN109296500A (zh) | 基于鲁棒控制理论的最大风能捕获方法 | |
EP3514373B1 (en) | Twist correction factor for aerodynamic performance map used in wind turbine control | |
Hafidi et al. | Wind speed estimation for wind turbine control | |
CN108547736A (zh) | 风速风向预测方法及风力发电机组的偏航控制方法 | |
CN106930898A (zh) | 一种基于自适应扰动补偿的漂浮式风力机功率控制方法 | |
CN113279907A (zh) | 用于操作风力涡轮机的方法 | |
CN115807734A (zh) | 一种基于尾流跟踪的海上风电场场级协同控制策略 | |
Pedersen et al. | Using wind speed from a blade-mounted flow sensor for power and load assessment on modern wind turbines | |
CN108537372B (zh) | 一种风向预测方法及风力发电机组的偏航控制方法 | |
CN101252334B (zh) | 变速恒频风力发电机动态最优能量捕获方法 | |
Bertašienė et al. | Synergies of Wind Turbine control techniques |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
DD01 | Delivery of document by public notice |
Addressee: Suzhou University of Science and Technology Document name: Notification to Make Rectification |
|
DD01 | Delivery of document by public notice | ||
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 |