CN102588211B - 一种风力发电机组全工况模型预测控制方法及系统 - Google Patents
一种风力发电机组全工况模型预测控制方法及系统 Download PDFInfo
- Publication number
- CN102588211B CN102588211B CN2012100504184A CN201210050418A CN102588211B CN 102588211 B CN102588211 B CN 102588211B CN 2012100504184 A CN2012100504184 A CN 2012100504184A CN 201210050418 A CN201210050418 A CN 201210050418A CN 102588211 B CN102588211 B CN 102588211B
- Authority
- CN
- China
- Prior art keywords
- control
- output
- variable
- wind
- state
- 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 36
- 238000006073 displacement reaction Methods 0.000 claims abstract description 68
- 239000011159 matrix material Substances 0.000 claims description 62
- 238000005457 optimization Methods 0.000 claims description 33
- 238000005070 sampling Methods 0.000 claims description 17
- 230000006870 function Effects 0.000 claims description 14
- 230000008859 change Effects 0.000 claims description 13
- 238000012937 correction Methods 0.000 claims description 13
- 238000005096 rolling process Methods 0.000 claims description 13
- 230000008878 coupling Effects 0.000 claims description 12
- 238000010168 coupling process Methods 0.000 claims description 12
- 238000005859 coupling reaction Methods 0.000 claims description 12
- 238000004891 communication Methods 0.000 claims description 9
- 238000013016 damping Methods 0.000 claims description 9
- 230000005611 electricity Effects 0.000 claims description 9
- 230000010355 oscillation Effects 0.000 claims description 9
- 238000005259 measurement Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 7
- 230000000694 effects Effects 0.000 claims description 5
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 238000013459 approach Methods 0.000 claims description 3
- 239000003638 chemical reducing agent Substances 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 230000005764 inhibitory process Effects 0.000 claims description 3
- 230000001915 proofreading effect Effects 0.000 claims description 3
- 239000000523 sample Substances 0.000 claims description 3
- 230000004044 response Effects 0.000 description 12
- 230000001276 controlling effect Effects 0.000 description 8
- 238000012360 testing method Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000006243 chemical reaction Methods 0.000 description 5
- 230000009471 action Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 4
- 238000011217 control strategy Methods 0.000 description 4
- 230000007812 deficiency Effects 0.000 description 4
- 238000011161 development Methods 0.000 description 2
- 238000012806 monitoring device Methods 0.000 description 2
- 238000004886 process control Methods 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 241000703769 Culter Species 0.000 description 1
- 239000004698 Polyethylene Substances 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004523 catalytic cracking Methods 0.000 description 1
- 238000004939 coking Methods 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000004821 distillation Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000012467 final product Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- -1 polyethylene Polymers 0.000 description 1
- 229920000573 polyethylene Polymers 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000003860 storage 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/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
- F05B2270/00—Control
- F05B2270/40—Type of control system
- F05B2270/404—Type of control system active, predictive, or anticipative
-
- 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)
- Feedback Control In General (AREA)
Abstract
一种风力发电机组全工况模型预测控制方法及系统,系统包括MPC控制器、反馈信息测量器、风轮、传动链、塔筒、发电单元、变桨驱动器、变流器,风轮、传动链、塔筒和发电单元的状态变量通过反馈信息测量器检测得到,检测结果传递给MPC控制器,通过MPC控制器计算出叶片桨距角和发电机转矩的目标,通过变桨驱动器和变流器对叶片桨距角和风力发电机转矩进行调整。本方法采用变桨控制预测模型和转矩控制预测模型计算控制增量,考虑了传动链扭转位移、传动链扭转速度、叶片平面外一阶挥舞位移、叶片平面外一阶挥舞速度、塔筒前后一阶摆动位移、塔筒前后一阶摆动速度、机组的机械载荷等状态变量,两个预测模型可在不同的工况下进行自动切换,实现全工况运行。
Description
技术领域
本发明涉及风力发电机组控制技术领域,具体涉及一种风力发电机组全工况模型预测控制方法及系统。
背景技术
模型预测控制(MPC)又称基于模型的预测控制(MBPC),最初由Richalet和Culter等人提出,从1968年模型预测控制的四个基本原则提出发展至今,已经在很多工业控制领域有所应用。模型预测控制包含四项原则,预测模型、参考响应轨迹、滚动优化、反馈校正。
模型预测控制的第一个原则是预测模型。预测控制是一种基于模型的算法,该模型称为预测模型。这是模型预测控制的基础,因为模型的好坏直接影响控制的质量,另外它也是模型预测控制不同于普通PID控制的根本所在。预测控制根据系统的未来输入及其历史信息预测未来输出值,预测模型强调其预测功能而不论其形式根据预测模型展示的未来行为,依据某种参考轨迹及优化标准,结合反馈信息,及时调整,以求出未来输入序列作为有自适应机制的预测控制,预测模型不是一成不变的,它直接或间接地被改变预测控制的发展。
模型预测控制的第二个原则是参考响应轨迹。由于控制量的幅值受实际条件的限制,只能在一定的范围内变化,而且控制量的变化速度一般也在一定的范围之内,另外被控对象往往受到不可测干扰的影响,这些因素使得被控变量很难在一拍之内跟踪设定值,一般需要经历一段时间输出才能跟踪设定值,因此根据设定值建立参考响应轨迹,使被控变量沿着参考响应轨迹逐步趋向设定值,这样可防止直接跟踪可能带来的大幅超调而导致的长时间振荡。它是设定值经过柔化的产物。参考响应轨迹概念是Propoi于1963年提出的。参考响应轨迹在模型预测控制中的应用是控制领域的一大进步,也是它优于预估控制的一个重要方面,它使得输出预测的时间长度可以超过时滞长度。
模型预测控制的第三个原则是滚动优化。模型预测控制的优化是一种反馈控制律,其控制量的大小是根据预测输出偏差来调整的。预测输出偏差指输出预测值偏离参考输出值的大小,参考输出值根据参考响应轨迹确定。预测控制的优化不同于传统意义下的优化,它是一种变时段的滚动优化。控制量的预测同输出预测一样不只是下一个控制时刻的预测,在每一采样时刻,优化性能指标涉及从该时刻起到未来的某时段,这不同于预估控制只预测下一个控制时刻的控制量的大小。在下一采样时刻,优化时段同时向后推移,优化是反复在线进行的。模型预测控制在每一个控制时刻只把未来第一个控制时刻的预测控制值应用于控制,其它时刻的预测控制值用于目标函数的性能指标计算。滚动优化是预测控制区别于传统最优控制的根本点。传统的最优控制是用一个性能函数来判断全局最优化,预测控制的优化不是一次离线完成,而是反复进行的,因此适用于动态特性变化和存在不确定因素的复杂工业过程。
模型预测控制的第四个原则是反馈校正。预测控制是一种闭环控制方法,通过滚动优化,确定一系列未来输入后,预测控制仅把当前时刻的控制输入作用于对象,而不是逐一实施。各个输入序列到下一采样时刻,检测实际输出,并根据这一实时信息直接或间接地对原预测模型进行修正(比如,优化时加入预测值与实际值的误差分量或根据实时输出信息在线辨识修改模型),然后进行新的优化。输入控制序列不像传统的最优控制一次算出,而是反复在线求出提高了对模型失配及环境干扰的鲁棒性。
预测控制系统的另一个优越之处就在于通过对具有鲜明物理意义的参数调节(如控制时域、优化时域、参考轨迹时间常数、性能指标中的加权系数),极易调节成好的鲁棒性这给操作实现带来极大的方便,因此,模型预测控制也是极方便的控制方法学。
模型预测控制已经被广泛地应用于催化裂化、原油蒸馏、延迟焦化、聚乙烯反应器等大型工业过程控制中,但在风力发电机组的控制中很少有人使用。风力发电机组在每个运行阶段的控制量都为了是系统趋向于其控制目标,可以理解为快速的过程控制,因此,模型预测控制方法完全可以应用于风力发电机组。
预测控制的主要特征是:以预测模型为基础,采用二次在线滚动优化性能指标和反馈校正等控制策略来消除被控对象建模误差,克服结构、参数与环境等不确定性因素的影响,有效地弥补了现代控制理论对复杂受控对象本身无法避免的不足。针对风电的被控对象建模困难,结构与参数具有不确定,风资源大扰动等特点,更具备适用性。
目前风力发电机组的控制器存在的一个重要问题就是控制作用有一定量的延迟。现有控制器的控制作用大多是通过对前一时刻测量值的评估而做出的,当前的控制量直接基于上一时刻的测量值,因此,如果当前时刻外作用发生波动,控制系统在下一时刻才能响应,从而使得控制作用有一定量的滞后。
中国专利公开号CN101560950A,公开日是2009年10月21日,名称为“基于预测控制风力涡轮机的方法和装置”中公开了一种通过预测叶片周围风场来进行风轮桨距控制的方法,该专利以控制叶片的载荷为目标,通过预测时间序列得到风场预测模型,进而利用空气动力学模型计算得到空气产生的载荷值,从而给出用于叶片调整的控制值,将风力涡轮机上的气流影响和相应的叶片调整之间的延迟时间降低。
该专利给出的方法能够弥补以往控制方法有控制延迟的缺陷,相比以往的控制方法,有很大的进步。但该控制方法本质上不是模型预测控制。首先,该控制的预测是对风速的预测,而不是对系统未来输出状态的预测;其次,其反馈校正是对风俗预测模型的校正,而不是对控制量误差的校正;然后,该控制方法没有使用滚动优化和参考轨迹。
因此该方法本身存在不足,如没有考虑动态特性变化和不确定因素的影响,同时可能带来大幅超调,导致长时间振荡。
风力发电机组的控制变量主要是桨距角和转矩,在不同的运行阶段,其控制目标不同,采用的控制变量也不同。一台变速型风力发电机组典型的运行阶段和其运行轨迹如图1所示。
阶段1是并网转速之前,控制系统主要是监控风速状况是否满足启动条件,满足条件则执行启动过程,打开桨距角增加双馈风电机组转速到并网转速切入电网后转换到阶段12。
阶段12是阶段1和阶段2的过渡阶段,按一个设定的斜率增加转矩给定到阶段2的最优给定曲线上即可。
阶段2是低风速最优功率跟踪阶段,在阶段2控制发电机转矩从风中捕获最大的能量。
阶段23是阶段2和阶段3的过渡阶段,按一个设定的斜率增加转矩给定到阶段3。
阶段3是变桨控制限制功率阶段,控制变桨保证机组输出功率稳定在额定值。这个阶段双馈风电机组承受载荷比较大,对塔筒的晃动,叶片的挥舞以及主传动链的扭振都必须加以控制。
由上分析可以看出,风电机组控制器的控制目标是:1)风电机组在整个运行范围内稳定可靠地在上述5个阶段运行;2)优化机组的运行性能,提高机组的发电效率与质量,减小机组的机械载荷。
《华北电力大学学报》2008年第35期中“基于混合逻辑动态模型的变速风力机全工况运行的控制研究”中介绍了风力机全工况运行下的混杂系统模型,该文结合风力机在不同工况点的非线性模型,引入混合逻辑动态模型的建模方法,建立风力机全工况运行下的混杂系统模型,并在此模型基础上,引入广义模型预测控制进行最优控制,使在风速变化的情况下,风力机可以良好的跟踪设定的功率曲线。
其不足之处是其控制目标只考虑了使风电机组在整个运行范围内稳定可靠地运行,没有考虑风电机组的运行性能优化,如提高机组的发电质量,减小机组的机械载荷等。并且,其控制模型中没有考虑风速扰动的影响,无法保证控制策略在风速扰动时实现机组的稳定运行和快速跟随。
现有技术中,预测风速却没有考虑动态特性变化和不确定因素的影响;没有参考轨迹可能带来大幅超调,导致长时间振荡;未考虑风电机组的运行性能优化,如提高机组的发电质量,减小机组的机械载荷等;没有考虑风速扰动的影响,无法保证控制策略在风速扰动时实现机组的稳定运行和快速跟随。
发明内容
针对现有技术存在的不足,本发明提供一种风力发电机组全工况模型预测控制方法及系统。
一种风力发电机组全工况模型预测控制系统,包括MPC控制器、反馈信息测量器、风轮、传动链、塔筒、发电单元、变桨驱动器、变流器;
风轮、传动链、塔筒和发电单元为被控对象,他们的状态变量通过反馈信息测量器检测得到;检测结果传递给MPC控制器,通过MPC控制器计算出叶片桨距角和发电机转矩的调整值;变桨驱动器和变流器为执行装置,通过它们对风轮和发电单元进行调整。
发电单元的输出端接至电网,发电单元风轮配备的变桨驱动器,调整变桨角度,实现风能最大利用率;变桨驱动器输入端与MPC控制器输出端口相连,变流器输入端连接到MPC控制器的输出端口,变流器输出端连接到发电单元。
所述反馈信息测量器包括安装于风力发电机主轴的转速传感器、安装于滑环尾部的风轮绝对位置编码器、安装于风力发电机尾部的发电机编码器、安装于机舱底盘的振动分析仪、安装于发电机定子侧的电压测量传感器、安装于发电机定子侧的电流测量传感器、装设与机舱外部支架上的风速传感器和安装于变距减速机旁的叶片桨距角测量编码器,各器件的输出端与MPC控制器的输入端口相连。风轮、传动链、塔筒和发电单元的状态变量通过反馈信息测量器采集得到。
所述MPC控制器,MPC控制器选用PLC,PLC包括机舱子站和塔底子站,机舱子站和塔底子站之间通过FAST总线通信。内部嵌入状态估计模块、预测模块、在线校正模块、参考轨迹模块和滚动优化模块,预测模块中,定义两个预测模型:转矩控制预测模型和变桨控制预测模型。两个预测模型可在不同的工况下进行自动切换,在额定风速和额定风速以下时投入控制量为转矩的预测模型,额定风速以上时投入控制量为桨距角的预测模型;参考轨迹模块根据系统的工况选择对应的参考轨迹;所述预测模型为计入风速扰动变量的离散状态空间模型;
本系统的工作过程如下:风轮、传动链、塔筒和发电单元的状态变量通过监测装置检测得到,检测结果传递给MPC控制器,通过MPC控制器计算出叶片桨距角和发电机转矩的调整值,通过变桨驱动器和变流器对叶片和风力发电机进行调整。
采用上述控制系统进行风力发电机组全工况模型预测控制的方法,具体按如下步骤进行:
步骤1:系统的离线准备;
MPC控制器的离线准备工作,包括:理论计算状态矩阵A,控制输入增益矩阵B,扰动输入增益矩阵Γ;选择控制目标,确定测量状态输出矩阵Cy,为被控状态输出矩阵Cz;根据机型确定被控输入和被控输出的约束条件、被控输出设定值zr随风速变化曲线;选择合适的采样周期T、预测时域N、控制时域Nc、参考轨迹时间常数Tg、校正向量h、输出误差权重Q、控制量变化幅度权重R;将以上数据存入存储模块,以便实时调用。
步骤2:反馈信息测量器实时检测风力发电机组的状态变量,通过风力发电机组全工况预测模型计算控制增量;
风力发电机组的状态变量包括:风轮转速、传动链扭转位移、传动链扭转速度、叶片平面外一阶挥舞位移、叶片平面外一阶挥舞速度、塔筒前后一阶摆动位移和塔筒前后一阶摆动速度;
计算控制增量采用两个预测模型:变桨控制预测模型和转矩控制预测模型,在不同工况下进行自动切换,在额定风速和额定风速以下时采用转矩控制预测模型,控制量为转矩,状态变量包括:风轮转速、传动链扭转位移、传动链扭转速度、叶片平面外一阶挥舞位移、叶片平面外一阶挥舞速度、塔筒前后一阶摆动位移和塔筒前后一阶摆动速度,并将风速增量作为扰动输入变量;在额定风速以上时采用变桨控制预测模型,控制量为桨距角,状态变量包括:风轮转速、传动链扭转位移和传动链扭转速度,并将风速增量作为扰动输入变量。
所述变桨控制预测模型,在离散化之前的状态空间模型如下:
所述转矩控制预测模型,离散化之前的状态空间模型如下:
变桨控制预测模型和转矩控制预测模型的统一状态空间标准形式为:
y=Cx
公式(1)~(3)的变量说明如下:
其中q1为叶片1平面外一阶模态位移,q2为叶片2平面外一阶模态位移,q3为叶片3平面外一阶模态位移,q4为风轮方位角,q5为发电机方位角,q6为塔筒前后一阶模态位移,为q的微分。
为控制增量输入向量,Δυ∈RO为扰动增量输入向量,y∈RP为测量输出向量。
Irot为风轮转动惯量,Igen为发电机转动惯量(把高速轴、低速轴和齿轮箱都加入到发电机系统),Tgen为发电机电磁转矩。Kd是传动链扭振的扭转刚度,Cd是传动链扭振阻尼。Mij为第i个状态变量和第j个状态变量之间的耦合质量系数,Cij为第i个状态变量和第j个状态变量之间的耦合阻尼系数,Kij为第i个状态变量和第j个状态变量之间的耦合刚度系数,其中,0<<i<<j<<6。
α=(δT_aero)/δυ,γ=(δT_aero)/(δω_r),ζ=(δT_aero)/δβ,其中T_aero为风轮的启动转矩,δυ,δω_r,δβ分别是风速、风轮转速和桨距角相对平衡点的增量值。
ζ_t=(δThrust_aero)/δβ,α_t=(δThrust_aero)/δυ,Thrust_aero为塔筒的气动推力。
对公式(3)的状态空间标准形式离散化,得到的线性离散状态空间模型如下:
x(k+1)=Ax(k)+BΔu(k)+ΓΔυ(k) (4)
y(k)=Cyx(k) (5)
z(k)=Czx(k) (6)
其中,x为状态向量,Δu为控制输入向量,Δυ扰动输入向量,y为测量输出向量,z为被控输出变量,Cy∈RP×N为测量状态输出矩阵,Cz∈RP×N为被控状态输出矩阵。
通过风力发电机组全工况预测模型计算控制增量,具体按如下步骤进行:
步骤2.1:初始化;
模型预测的在线计算由初始化与实时控制组成。
初始化是通过反馈信息测量器检测预测模型的状态变量x(k),并把它设定为下一步的模型预测初值,从步骤2.2起即转入实时控制,在每一个采样时刻的在线计算流程如图3所示。
步骤2.2:采用卡尔曼滤波算法对不可测状态变量叶片一阶挥舞位移进行状态估计,将估计后的状态变量输出;
在每一个采样时刻,都需要预测未来一段时间的被控输出变量值。
假设当前时刻为k,信号Δu(k+j|k)表示输入Δu的未来在(k+j)时刻的值,信号x(k+j|k)、z(k+J|k)的含义同Δu(k+j|k)一样。
在每一时刻k,要确定从该时起的起Nc个控制增量,使被控对象在其作用下未来的N个时刻的输出预测值尽可能接近期望。这里,N为预测时域,Nc为控制时域,总是假设Nc<N。这意味着当j>Nc时Δu(k+j|k)=0,即u(k+j|k)=u(k+Nc-1|k),j>Nc,也就是说,当j>Nc时,一个零阶保持器作用在输入上。
利用公式(4)进行迭代计算,
当j≤Nc时有
当Nc≤j≤N时有
在每一时刻k的预测时域N内总的预测方程为:
定义P为状态变量的系数矩阵,Gz为控制输入增量的系数矩阵,Gυ为控制输入增量的系数矩阵。
则式(9)可写为:
控制输出z的预测方程(6)可以写作:
z(k+j|k)=Czx(k+j|k)j =1,...,N (11)
令
则公式(10)可以写作:
Z(k)=CzPx(k)+CzGzΔU(k)+CzGuΔυ(k) (12)
预测模型中的不可直接测量的状态变量:叶片平面外一阶挥舞位移,需要进行状态估计,状态估计采用卡尔曼滤波器算法,具体如下:
计算当前的估计状态:
M(k)=P(k|k-1)Cz T[CzP(k|k-1)Cz T+R]-1 (14)
其中P(k|k-1)是误差方差阵:P(k|k)=[I-M(k)Cz]P(k|k-1)
计算下一步更新状态和误差方差阵:
P(k+1|k)=AP(k|k)AT+BQBT (16)
其中R,Q是高斯白噪声的协方差阵。
步骤2.3:根据风速增量和估计后的状态变量,判断系统所处的运行工况,选择匹配的预测模型,即式(1)或式(2);
若风速在额定风速和额定风速以下时投入转矩控制预测模型,控制量为转矩,状态变量包括:风轮转速、传动链扭转位移、传动链扭转速度、叶片平面外一阶挥舞位移、叶片平面外一阶挥舞速度、塔筒前后一阶摆动位移和塔筒前后一阶摆动速度,并将风速增量作为扰动输入变量;若风速在额定风速以上时投入变桨控制预测模型,控制量为桨距角,状态变量包括:风轮转速、传动链扭转位移和传动链扭转速度,并将风速增量作为扰动输入变量。
步骤2.4:利用检测到的风力发电机组的状态变量和所选择的预测模型参数,根据式(12)计算控制时域内的控制输出变量;
步骤2.5:将预测值与上一时刻的实际输出进行比较,得出输出误差,如式(18);
在下一采样时刻首先检测对象的实际输出y(k+1),并把它与式(11)算出的模型预测输出z(k+1)相比较,构成输出误差
e(k+1)=y(k+1)-z(k+1|k) (18)
步骤2.6:将下一时刻的预测模型输出与本时刻的实际输出进行比较,得出输出误差,输出误差与本时刻的预测模型输出值叠加,作为在线校正后的预测值,如式(19);
采用对e(k+1)加权的方式修正对未来输出的预测:
其中 为校正后的输出预测向量,有权系数组成的N维向量 为校正向量。
步骤2.7:根据在线校正过的预测值,确定参考轨迹;
采用的参考轨迹是一阶指数变化形式
zr(k+j)=ajz(k)+(1-aj)zr (17)
目标函数定义如下:
约束条件为:
Δumin≤Δu(k+j|k)≤Δumax j=0,…,Nc-1 (21)
其中Q(j),R(j)是权系数,分别表示对跟踪误差和控制量变化的抑制。
将目标函数(20)写成向量形式:
其中
由权系数构成的对角阵Q和R分别称为误差权矩阵和控制权矩阵。
定义:
这个向量可以认为是轨迹误差,也就是目标轨迹和系统自由响应之间的偏差。这样,可以将目标函数改写为:
最终可以写成如下的标准形式:
其中H=2((CzGz)TQ(CzGz)+R),f=-2E(k)TQCzGz。
约束条件(18)可以写成一个矩阵的形式
其中b1,b2的维数为Nc×nu,分别包括了Nc个Δumax,Δumin。同样,z1,z2的维数为N,包含了N个zmax,zmin。
这样,问题转化为求解如下的有约束优化问题
并且满足公式(23),滚动优化策略要求在每个控制周期仅仅使用优化解得第一列(控制时域中的第一步),这样即时控制增量可写做:
步骤3:经风力发电机组全工况模型预测计算出的控制增量,作为下一时刻的桨距角增量或转矩增量,通过变桨驱动器和变流器对叶片和和发电机进行调整控制。
有益效果:
1.本发明控制系统采用的模型预测控制方法,由于具有滚动优化模块和参考轨迹模块,系统本身考虑了动态特性变化和不确定因素的影响;同时避免了系统大幅超调,导致长时间振荡的现象;
2.本发明的预测模型包含传动链扭转位移、传动链扭转速度、叶片平面外一阶挥舞位移、叶片平面外一阶挥舞速度、塔筒前后一阶摆动位移、塔筒前后一阶摆动速度,减小机组的机械载荷等状态变量,同时在优化模型中增加了各状态变量的约束条件,可以有效控制机组的机械载荷;
3.本发明的预测模型包含风速扰动变量,将每一采样时刻的风速扰动都计入了控制策略中,可以保证在风速扰动时实现机组的稳定运行和快速跟随;
4.本发明方法有两个预测模型,变桨控制预测模型和转矩控制预测模型,两个预测模型可在不同的工况下进行自动切换,实现了全工况运行。
附图说明
图1是一台变速型风力发电机组典型的运行阶段和其运行轨迹图;
图2是本发明实施例风力发电机组全工况模型预测控制原理图;
图3是本发明实施例模型MPC控制器原理框图;
图4为本发明实施例在线计算流程图;
图5为本发明实施例模型预测控制器PLC配置方案构成图;
图6为本发明实施例比较阶跃风输入下系统的响应曲线图;
图7为本发明实施例比较阶跃风输入下叶片挥舞位移曲线图;
图8为本发明实施例比较阶跃风输入下塔筒顶部前后摆动位移曲线图;
图9为本发明实施例试验湍流风输入曲线图;
图10为本发明实施例比较湍流风输入系统响应曲线图;
图11为本发明实施例比较湍流风输入下叶片挥舞位移曲线图;
图12为本发明实施例比较湍流风输入塔筒顶部前后摆动位移曲线图。
具体实施方式
下面结合附图和实施例对本发明做进一步说明。
本实施例采用以3MW双馈风力发电机组的全工况模型预测控制系统,包括MPC控制器1、反馈信息测量器6、风轮3、传动链4、塔筒5、发电单元8、变桨驱动器2、变流器7、控制原理如图2所示,发电单元的输出端接至电网,发电单元风轮配备的变桨驱动器,调整叶片桨距角,实现风能最大利用率;变桨驱动器输入端与MPC控制器输出端口相连,变流器输入端连接到MPC控制器的输出端口,变流器输出端连接到发电单元。
所述MPC控制器,内部嵌入状态估计模块、预测模块、在线校正模块、参考轨迹模块和滚动优化模块,原理如图3所示。预测模块中,定义两个预测模型:转矩控制预测模型和变桨控制预测模型。两个预测模型可在不同的工况下进行自动切换,在额定风速和额定风速以下时投入控制量为转矩的预测模型,额定风速以上时投入控制量为桨距角的预测模型;参考轨迹模块根据系统的工况选择对应的参考轨迹;所述预测模型为计入风速扰动变量的离散状态空间模型。
MPC控制器选用PLC,PLC包括机舱子站和塔底子站,分别置于机舱控制柜和塔底控制柜中,机舱子站和塔底子站之间通过FAST总线通信,如图5所示。
机舱子站配置如下:FS211/N为FAST总线子模块;DIO280为数字输入/输出模块;DI232为数字输入模块;PTAI216为温度记录模块;AIO288为模拟输入输出模块;CNT204为计数模块;RS204为串口通信模块;CM202为CAN总线模块;EM203为以太网主站模块;LM201为占位模块;WT310为显示屏。
图中塔底子站配置如下:SWI205为工业以太网交换机模块;NT250为电源模块;MPC240为处理器;PCC201为存储卡模块;FM211为FAST总线主模块;DIO248为数字输入/输出模块;PTAI216为温度记录模块;LM201为占位模块。
AIO288模块、RS204模块是MPC控制器与反馈信息测量器6的通信接口。反馈信息测量器6包括:安装于风力发电机主轴的转速传感器(普通接近开关)、安装于滑环尾部的绝对位置编码器(型号为HENGSTLER-AC58)、安装于风力发电机尾部的发电机编码器(型号为LeineLinde-RHI503)、安装于机舱底盘的振动分析仪(型号为PCH1024)、安装于发电机定子侧的电压测量传感器(型号为SIEMENS-3RT)、安装于发电机定子侧的电流测量传感器(型号为SIEMENS-ESM500-9920)、装设与机舱外部支架上的风速传感器(型号为FT702)和安装于变距减速机旁的叶片位置编码器(型号为Baumer-SSI)。
CM202模块是MPC控制器与变桨控制器2、变流器7的通信接口,它们之间通过CAN总线通信。
FS211模块与FM211模块是机舱柜和塔底柜通信的接口模块,通过FAST总线通信。
SWI205模块是MPC控制器与塔底编程/维护设备、风场SCADA、Mobile SMS的通信接口。
PCC201为存储卡模块,存储的数据有状态矩阵A、控制输入增益矩阵B、扰动输入增益矩阵Γ、测量状态输出矩阵Cy、被控状态输出矩阵Cz、被控输入和被控输出的约束条件、被控输出设定值zr随风速变化曲线、采样周期T、预测时域N、控制时域Nc、参考轨迹时间常数Tg、校正向量h、输出误差权重Q、控制量变化幅度权重R、历史风速、过程变量,以便处理器模块MPC240随时调用。
系统的工作过程如下:风轮、传动链、塔筒和发电单元的状态变量通过监测装置检测得到,检测结果传递给MPC控制器,通过MPC控制器计算出叶片桨距角和发电机转矩的调整值,通过变桨驱动器和变流器对叶片和风力发电机进行调整。
采用上述控制系统进行风力发电机组全工况模型预测控制的方法,具体按如下步骤进行:
步骤1:系统的离线准备;
MPC控制器的离线准备工作,包括:理论计算状态矩阵A,控制输入增益矩阵B,扰动输入增益矩阵Γ;选择控制目标,确定测量状态输出矩阵Cy,为被控状态输出矩阵Cz;根据机型确定被控输入和被控输出的约束条件、被控输出设定值zr随风速变化曲线;选择合适的采样周期T、预测时域N、控制时域Nc、参考轨迹时间常数Tg、校正向量h、输出误差权重Q、控制量变化幅度权重R;将以上数据存入存储模块,以便实时调用。
步骤2:反馈信息测量器实时检测风力发电机组的状态变量,通过风力发电机组进行全工况预测模型计算控制增量;
风力发电机组的状态变量包括:风轮转速、传动链扭转位移、传动链扭转速度、叶片平面外一阶挥舞位移、叶片平面外一阶挥舞速度、塔筒前后一阶摆动位移和塔筒前后一阶摆动速度;
两个预测模型:变桨控制预测模型和转矩控制预测模型,在不同工况下进行自动切换,在额定风速和额定风速以下时投入转矩控制预测模型,控制量为转矩,状态变量包括:风轮转速、传动链扭转位移、传动链扭转速度、叶片平面外一阶挥舞位移、叶片平面外一阶挥舞速度、塔筒前后一阶摆动位移和塔筒前后一阶摆动速度,并将风速增量作为扰动输入变量;在额定风速以上时投入变桨控制预测模型,控制量为桨距角,状态变量包括:风轮转速、传动链扭转位移和传动链扭转速度,并将风速增量作为扰动输入变量。
所述变桨控制预测模型,在离散化之前的状态空间模型如下:
所述转矩控制预测模型,离散化之前的状态空间模型如下:
变桨控制预测模型和转矩控制预测模型的统一状态空间标准形式为:
y=Cx
公式(1)~(3)的变量说明如下:
Irot为风轮转动惯量,Igen为发电机转动惯量(把高速轴、低速轴和齿轮箱都加入到发电机系统),Tgen为发电机电磁转矩。kd是传动链扭振的扭转刚度,Cd是传动链扭振阻尼。Mij为第i个状态变量和第j个状态变量之间的耦合质量系数,Cij为第i个状态变量和第j个状态变量之间的耦合阻尼系数,Kij为第i个状态变量和第j个状态变量之间的耦合刚度系数,其中,0<<i<<j<<6。
α=(δT_aero)/δυ,γ=(δT_aero)/(δω_r ),ζ=(δT_aero)/δβ,其中T_aero为风轮的启动转矩,δυ,δω_r,δβ分别是风速、风轮转速和桨距角相对平衡点的增量值。ζ_t=(δThrust_aero)/δβ,α_t=(δThrust=aero)/δυ,Thrust_aero为塔筒的气动推力。 其中δFaeroi为叶片i平面外一阶挥舞推理增量,βi,υu分别为叶片i的桨距角和叶片i处的风速,i=1,2,3.ζb,αb分别为ζbi和αbi的平均值。α,γ,ζ,α_t,ζt,αi,ζi,αbi,ζbi均为过程变量。
对公式(3)的状态空间标准形式离散化,得到的线性离散状态空间模型如下:
x(k+1)=Ax(k)+BΔu(k)+ΓΔυ(k) (4)
y(k)=Cyx(k) (5)
z(k)=Czx(k) (6)
其中,x为状态向量,Δu为控制输入向量,Δυ为扰动输入向量,y为测量输出向量,z为被控输出变量,Cy∈RP×N为测量状态输出矩阵,Cx∈RP×N为被控状态输出矩阵。
通过风力发电机组全工况预测模型计算控制增量,具体按如下步骤进行:
步骤2.1:初始化;
模型预测的在线计算由初始化与实时控制组成,流程如图4所示。
初始化是通过反馈信息测量器检测预测模型的状态变量x(k),并把它设定为下一步的模型预测初值,从步骤2.2起即转入实时控制,在每一个采样时刻的在线计算流程如图3所示。
步骤2.2:采用卡尔曼滤波算法对不可测状态变量叶片一阶挥舞位移进行状态估计,将估计后的状态变量输出;
在每一个采样时刻,都需要预测未来一段时间的被控输出变量值。
假设当前时刻为k,信号Δu(k+j|k)表示输入Δu的未来在(k+j)时刻的值,信号x(k+j|k)、z(k+j|k)的含义同Δ(k+j|k)一样。
在每一时刻k,要确定从该时起的起Nc个控制增量,使被控对象在其作用下未来的N个时刻的输出预测值尽可能接近期望。这里,N为预测时域,Nc为控制时域,总是假设Nc<N。这意味着当j>Nc时Δu(k+j|k)=0,即u(k+j|k)=u(k+Nc-1|k),j>Nc,也就是说,当j>Nc时,一个零阶保持器作用在输入上。
利用公式(4)进行迭代计算,
当j≤Nc时有
当Nc≤j≤N时有
在每一时刻k的预测时域N内总的预测方程为:
定义P为状态变量的系数矩阵,Cz为控制输入增量的系数矩阵,Cυ为控制输入增量的系数矩阵。
则式(9)可写为:
控制输出z的预测方程(6)可以写作:
z(k+j|k)=Czx(k+j|k)j=1,...,N (11)
令
则公式(10)可以写作:
Z(k)=CzPx(k)+CzGzΔU(k)+CzGuΔυ(k) (12)
预测模型中的不可直接测量的状态变量:叶片平面外一阶挥舞位移,需要进行状态估计,状态估计采用卡尔曼滤波器算法,具体如下:
计算当前的估计状态:
M(k)=P(k|k-1)Cz T[CzP(k|k-1)Cz T+R]-1 (14)
其中P(k|k-1)是误差方差阵:P(k|k)=[I-M(k)Cz]P(k|k-1)
计算下一步更新状态和误差方差阵:
P(k+1|k)=AP(k|k)AT+BQBT (16)
其中R,Q是高斯白噪声的协方差阵。
步骤2.3:根据风速增量和估计后的状态变量,判断系统所处的运行工况,选择匹配的预测模型,即式(1)或式(2);
若风速在额定风速和额定风速以下时采用转矩控制预测模型,控制量为转矩,状态变量包括:风轮转速、传动链扭转位移、传动链扭转速度、叶片平面外一阶挥舞位移、叶片平面外一阶挥舞速度、塔筒前后一阶摆动位移和塔筒前后一阶摆动速度,并将风速增量作为扰动输入变量;若风速在额定风速以上时采用变桨控制预测模型,控制量为桨距角,状态变量包括:风轮转速、传动链扭转位移和传动链扭转速度,并将风速增量作为扰动输入变量。
步骤2.4:利用检测到的风力发电机组的状态变量和所选择的预测模型参数,根据式(12)计算控制时域内的控制输出变量;
步骤2.5:将预测值与上一时刻的实际输出进行比较,得出输出误差,如式(18);
在下一采样时刻首先检测对象的实际输出y(k+1),并把它与式(11)算出的模型预测输出z(k+1)相比较,构成输出误差
e(k+1)=y(k+1)-z(k+1|k) (18)
步骤2.6:将下一时刻的预测模型输出与本时刻的实际输出进行比较,得出输出误差,输出误差与本时刻的预测模型输出值叠加,作为在线校正后的预测值,如式(19);
采用对e(k+1)加权的方式修正对未来输出的预测:
其中 为校正后的输出预测向量,有权系数组成的N维向量 为校正向量。
步骤2.7:根据在线校正过的预测值,确定参考轨迹;
采用的参考轨迹是一阶指数变化形式
zr(k+j)=ajz(k)+(1-aj)zr (17)
目标函数定义如下:
约束条件为:
Δumin≤Δu(k+j|k)≤Δumax j=0,…,Nc-1 (21)
其中Q(j),R(j)是权系数,分别表示对跟踪误差和控制量变化的抑制。
将目标函数(20)写成向量形式:
其中
由权系数构成的对角阵Q和R分别称为误差权矩阵和控制权矩阵。
定义:
这个向量可以认为是轨迹误差,也就是目标轨迹和系统自由响应之间的偏差。这样,可以将目标函数改写为:
最终可以写成如下的标准形式:
其中H=2((CzGz)TQ(CzGz)+R),f=-2E(k)TQCzGz。
约束条件(18)可以写成一个矩阵的形式
其中b1,b2的维数为Nc×nu,分别包括了Nc个Δumax,Δumin。同样,z1,z2的维数为N,包含了N个zmax,zmin。
这样,问题转化为求解如下的有约束优化问题
并且满足公式(23),滚动优化策略要求在每个控制周期仅仅使用优化解得第一列(控制时域中的第一步),这样即时控制增量可写作:
步骤3:经风力发电机组全工况模型预测计算出的控制增量,作为下一时刻的桨距角增量或转矩增量,通过变桨驱动器和变流器对叶片和和发电机进行调整控制。
本系统变桨控制预测模型和转矩控制预测模型切换风速为14m/s,即风速在14m/s以下时使用转矩控制预测模型,在14m/s以上时使用变桨控制预测模型。
其变桨控制线性离散预测模型是将变桨控制预测模型在阶段3的稳态工作点按控制周期0.02s离散化后获得,其阶段3的稳态工作点条件为:
风速υ=18(m/s);风轮转速ωr=14.25(rpm);叶片桨距角β=16.66(deg)。
其变桨控制状态线性离散预测模型参数如表1所示。
表13 MW风力发电机组变桨控制7状态线性离散预测模型参数
其转矩控制线性离散预测模型是将转矩控制预测模型在阶段2的稳态工作点按控制周期0.02s离散化后获得,其阶段2的稳态工作点条件为:
风速υ=13(m/s);风轮转速ωr=14.25(rpm);叶片桨距角β=6.941(deg)。
其转矩控制状态线性离散预测模型参数如表2所示。
表23 MW风力发电机组转矩控制3状态线性离散预测模型参数
试验阶跃风速输入下风电机组响应,风速输入为17m/s在250s处跳跃为18m/s。
控制目标(被控输出设定值zr)为:稳定风轮转速为14.25rpm、增加传动链模态阻尼、增加叶片挥舞模态阻尼和塔筒前后模态阻尼,传动链扭振、叶片一阶挥舞和塔筒前后一阶摆动偏移为零。
模型的输入输出约束为:
-9(deg/s)≤Δβ≤9(deg/m)
8.9(rpm)≤ωr≤16(rpm)
采样周期T为0.02秒,参考轨迹时间常数Tg为0.1,预测时域N、控制时域Nc选择为10和1,校正向量 输出误差权重 控制量变化幅度权重R为300。
在MATLAB中进行仿真,变增益PI控制和模型预测控制采用相同的输入,图6为阶跃风输入下风轮转速响应,图7为叶片平面外一阶挥舞位移,图8为塔筒顶部前后一阶摆动位移,可以看出模型预测控制在稳定传动链扭振、减小叶片挥舞载荷和减小塔筒前后摆动载荷方面效果明显。
对湍流风输入测试,结果如下:图9是测试湍流风,图10是湍流风输入下风轮转速响应,图11是湍流风输入下的叶片一阶挥舞位移,图12是湍流风输入下塔筒顶部前后一阶摆动位移,可以看出湍流风速输入下模型预测控制在载荷控制方面效果也是很明显的。
Claims (2)
1.一种风力发电机组全工况模型预测控制方法,该方法采用的风力发电机组全工况模型预测控制系统,包括风轮、传动链、塔筒、发电单元、变桨驱动器、变流器、反馈信息测量器和MPC控制器;
发电单元的输出端接至电网,变桨驱动器输入端与MPC控制器输出端口相连,变流器输入端连接到MPC的输出端口,变流器输出端连接到发电单元;
所述反馈信息测量器包括安装于风力发电机主轴的转速传感器、安装于滑环尾部的风轮绝对位置编码器、安装于风力发电机尾部的发电机编码器、安装于机舱底盘的振动分析仪、安装于发电机定子侧的电压测量传感器、安装于发电机定子侧的电流测量传感器、装设与机舱外部支架上的风速传感器和安装于变距减速机旁的叶片桨距角编码器,各反馈信息测量器的输出端与MPC控制器的输入端口相连;
所述MPC控制器选用PLC作为硬件系统,PLC包括机舱子站和塔底子站,机舱子站和塔底子站之间通过FAST总线通信;
其特征在于:该方法具体按如下步骤进行:
步骤1:系统的离线准备;
MPC控制器的离线准备工作,包括:理论计算:状态矩阵A,控制输入增益矩阵B,扰动输入增益矩阵Γ;选择控制目标,确定测量状态输出矩阵Cy,为被控状态输出矩阵Cz;根据机型确定被控输入和被控输出的约束条件、被控输出设定值zi随风速变化曲线;选择合适的采样周期T、预测时域N、控制时域Nc、参考轨迹时间常数Ts、校正向量h、输出误差权重Q、控制量变化幅度权重R;将以上数据存入存储模块,以便实时调用;
步骤2:反馈信息测量器实时检测风力发电机组的状态变量,通过风力发电机组全工况预测模型计算控制增量;
风力发电机组的状态变量包括:风轮转速、传动链扭转位移、传动链扭转速度、叶片平面外一阶挥舞位移、叶片平面外一阶挥舞速度、塔筒前后一阶摆动位移和塔筒前后一阶摆动速度;
计算控制增量采用两个预测模型:变桨控制预测模型和转矩控制预测模型,在不同工况下进行自动切换,在额定风速和额定风速以下时采用转矩控制预测模型,控制量为转矩,状态变量包括:风轮转速、传动链扭转位移、传动链扭转速度、叶片平面外一阶挥舞位移、叶片平面外一阶挥舞速度、塔筒前后一阶摆动位移和塔筒前后一阶摆动速度,并将风速增量作为扰动输入变量;在额定风速以上时采用变桨控制预测模型,控制量为桨距角,状态变量包 括:风轮转速、传动链扭转位移和传动链扭转速度,并将风速增量作为扰动输入变量;
所述变桨控制预测模型,在离散化之前的状态空间模型如下:
所述转矩控制预测模型,离散化之前的状态空间模型如下:
变桨控制预测模型和转矩控制预测模型的统一状态空间标准形式为:
y=Cx
公式(1)~(3)的变量说明如下:
Irot为风轮转动惯量,Igen为发电机转动惯量,Tgen为发电机电磁转矩;Kd是传动链扭振的扭转刚度,Cd是传动链扭振阻尼;Mij为第i个状态变量和第j个状态变量之间的耦合质量 系数,Cij为第i个状态变量和第j个状态变量之间的耦合阻尼系数,Kij为第i个状态变量和第j个状态变量之间的耦合刚度系数,其中,0<<i<<j<<6;
α=(δT_aero)/δυ,γ=(δT_aero)/(δω_r),ζ=(δT_aero)/δβ,其中T_aero为风轮的启动转矩,δυ,δω_r,δβ分别是风速、风轮转速和桨距角相对平衡点的增量值;
ζ_t=(δThrust_aero)/δβ,α_t=(δThrust_aero)/δυ,Thrust_aero为塔筒的气动推力; 其中δFaeroi为叶片i平面外一阶挥舞推理增量,βi,υi分别为叶片i的桨距角和叶片i处的风速,i=1,2,3.ζb,αb分别为ζbi和αbi的平均值;α,γ,ζ,α_t,ζτ,αi,ζi,αbi,ζbi均为过程变量;
对公式(3)的状态空间标准形式离散化,得到的线性离散状态空间模型如下:
x(k+1)=Ax(k)+BΔu(k)+ΓΔυ(k) (4)
y(k)=Cyx(k) (5)
z(k)=Czx(k) (6)
其中,x为状态向量,Δu为控制输入向量,Δυ为扰动输入向量,y为测量输出向量,z为被控输出变量,Cy∈Rp×N为测量状态输出矩阵,Cz∈Rp×N为被控状态输出矩阵;
步骤3:经风力发电机组全工况预测模型计算出的控制增量,作为下一时刻的桨距角增量或转矩增量,通过变桨驱动器和变流器对叶片和和发电机进行调整控制。
2.根据权利要求1所述的风力发电机组全工况模型预测控制方法,其特征在于:步骤2所述的通过风力发电机组全工况预测模型计算控制增量,具体按如下步骤进行:
步骤2.1:初始化;
初始化是通过反馈信息测量器检测预测模型的状态变量x(k),并把它设定为下一步的模 型预测初值,从步骤2.2起即转入实时控制;
步骤2.2:采用卡尔曼滤波算法对不可测状态变量叶片一阶挥舞位移进行状态估计,将估计后的状态变量输出;
在每一个采样时刻,都需要预测未来一段时间的被控输出变量值;
设当前时刻为k,信号Δu(k+j|k)表示输入Δu的未来在(k+j)时刻的值,信号x(k+j|k)、z(k+j|k)的含义同Δu(k+j|k)一样;
在每一时刻k,要确定从该时起的起Nc个控制增量,使被控对象在其作用下未来的N个时刻的输出预测值尽可能接近期望,N为预测时域,Nc为控制时域,总是假设Nc<N,即当j>Nc时Δu(k+j|k)=0,即u(k+k|k)=u(k+Nc-1|k),j>Nc,也就是说,当j>Nc时,一个零阶保持器作用在输入上;
利用公式(4)进行迭代计算,
当j≤Nc时有
当Nc≤j≤B时有
在每一时刻k的预测时域N内总的预测方程为:
定义P为状态变量的系数矩阵,Gz为控制输入增量的系数矩阵,Gυ为控制输入增量的系数矩阵,
则式(9)可写为:
控制输出z的预测方程(6)可以写作:
z(k+j|k)=Czx(k+k|k) j=1,...,N (11)
则公式(10)可以写作:
Z(k)=CzPx(k)+CzGzΔU(k)+CzGυΔυ(k) (12)
预测模型中的不可直接测量的状态变量:叶片平面外一阶挥舞位移,需要进行状态估计,状态估计采用卡尔曼滤波器算法,具体如下:
计算当前的估计状态:
M(k)=P(k|k-1)Cz T[CzP(k|k-1)Cz T+R]-1 (14)
其中P(k|k-1)是误差方差阵:P(k|k)=[1-M(k)Cz]P(k|k-1)
计算下一步更新状态和误差方差阵:
P(k+1|k)=AP(k|k)AT+BQBT (16)
其中R,Q是高斯白噪声的协方差阵;
步骤2.3:根据风速增量和估计后的状态变量,判断系统所处的运行工况,选择匹配的预测模型,即式(1)或式(2);
若风速在额定风速和额定风速以下时投入转矩控制预测模型,控制量为转矩,状态变量包括:风轮转速、传动链扭转位移、传动链扭转速度、叶片平面外一阶挥舞位移、叶片平面外一阶挥舞速度、塔筒前后一阶摆动位移和塔筒前后一阶摆动速度,并将风速增量作为扰动输入变量;若风速在额定风速以上时投入变桨控制预测模型,控制量为桨距角,状态变量包括:风轮转速、传动链扭转位移和传动链扭转速度,并将风速增量作为扰动输入变量;
步骤2.4:利用检测到的风力发电机组的状态变量和所选择的预测模型参数,根据式(12)计算控制时域内的控制输出变量;
步骤2·5:将预测值与上一时刻的实际输出进行比较,得出输出误差,如式(18);
在下一采样时刻首先检测对象的实际输出y(k+1),并把它与式(11)算出的模型预测输出z(k+1)相比较,构成输出误差
e(k+1)=y(k+1)-z(k+1|k) (18)
步骤2·6:将下一时刻的预测模型输出与本时刻的实际输出进行比较,得出输出误差,输出误差与本时刻的预测模型输出值叠加,作为在线校正后的预测值,如式(19);
采用对e(k+l)加权的方式修正对未来输出的预测:
步骤2.7:根据在线校正过的预测值,确定参考轨迹;
采用的参考轨迹是一阶指数变化形式
zr(k+j)=ajz(k)+(1-aj)zr (17)
目标函数定义如下:
约束条件为:
其中Q(j),R(j)是权系数,分别表示对跟踪误差和控制量变化的抑制;
将目标函数(20)写成向量形式:
其中
由权系数构成的对角阵Q和R分别称为误差权矩阵和控制权矩阵;
定义:
这个向量可以认为是轨迹误差,也就是目标轨迹和系统自由响应之间的偏差;这样,将目标函数改写为:
最终可以写成如下的标准形式:
其中H=2((CzGz)TQ(CzGz)+R),f=-2E(k)TQCzGz,
约束条件(18)可以写成一个矩阵的形式
其中bl,b2的维数为Nc×nυ,分别包括了Nc个Δumax,Δumin;同样,z1,z2的维数为N,包含了N个zmax,zmin;
这样,问题转化为求解如下的有约束优化问题
并且满足公式(23),滚动优化要求在每个控制周期仅仅使用优化解得第一列即控制时域中的第一步,这样即时控制增量可写做:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012100504184A CN102588211B (zh) | 2012-02-29 | 2012-02-29 | 一种风力发电机组全工况模型预测控制方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012100504184A CN102588211B (zh) | 2012-02-29 | 2012-02-29 | 一种风力发电机组全工况模型预测控制方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102588211A CN102588211A (zh) | 2012-07-18 |
CN102588211B true CN102588211B (zh) | 2013-08-14 |
Family
ID=46477329
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2012100504184A Active CN102588211B (zh) | 2012-02-29 | 2012-02-29 | 一种风力发电机组全工况模型预测控制方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102588211B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107869421A (zh) * | 2016-09-26 | 2018-04-03 | 北京金风科创风电设备有限公司 | 风力发电机变桨系统的控制方法和装置 |
CN110145436A (zh) * | 2019-04-28 | 2019-08-20 | 华北电力大学 | 应用于风机的非线性经济模型预测控制方法 |
CN111305999A (zh) * | 2020-04-01 | 2020-06-19 | 宁夏三聚阳光新能源科技有限公司 | 一种智能调整风力发电系统 |
US12106022B2 (en) * | 2023-02-02 | 2024-10-01 | North China Electric Power University | Method and device for simulating dynamic digital twin model of dominant operation of wind turbine generator assembly |
Families Citing this family (43)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102801380B (zh) * | 2012-08-24 | 2015-07-15 | 国电联合动力技术有限公司 | 风电机组传动链阻尼转矩校正的控制方法及装置 |
CN103092129B (zh) * | 2013-01-29 | 2014-10-22 | 电子科技大学 | 一种用于复杂数字控制系统的预测控制方法 |
CN103411659B (zh) * | 2013-08-12 | 2015-06-03 | 国电联合动力技术有限公司 | 一种风力发电机叶片与塔筒状态监测方法及系统 |
CN103455028B (zh) * | 2013-08-29 | 2016-04-20 | 国家电网公司 | 一种风电机组控制系统pid环节静态测试及校验方法 |
CN103758694B (zh) * | 2013-12-24 | 2016-05-25 | 青海能高新能源有限公司 | 一种风力发电机组状态监控系统数据压缩方法 |
CN103983454B (zh) * | 2014-05-26 | 2017-02-01 | 湖南科技大学 | 一种风电机组传动链全工况测试模拟装置 |
CN105320793A (zh) * | 2014-07-30 | 2016-02-10 | 南车株洲电力机车研究所有限公司 | 一种从动力学角度评估风电机组变浆控制模型的方法 |
CN104234934B (zh) * | 2014-08-06 | 2017-02-08 | 东北大学 | 一种风力发电机转速控制方法 |
ES2762519T3 (es) * | 2014-08-15 | 2020-05-25 | Vestas Wind Sys As | Control de un aerogenerador con una condición de fallo |
CN104242339B (zh) * | 2014-08-29 | 2017-02-15 | 清华大学 | 基于模型预测控制理论的风电场电压自动控制方法 |
CN104242768B (zh) * | 2014-09-11 | 2016-01-20 | 天津大学 | 一种应用于多电机控制系统的有限状态模型预测控制方法 |
US10691088B2 (en) | 2015-03-20 | 2020-06-23 | The Boeing Company | Programmable actuator simulation card |
CN105179181B (zh) * | 2015-09-15 | 2018-11-06 | 苏州格远电气有限公司 | 风电变桨驱动器的状态显示装置及显示方法 |
CN105569930B (zh) * | 2015-12-21 | 2018-06-12 | 南京理工大学 | 一种基于偏差抑制的风力机模拟器转动惯量补偿改进方法 |
CN105894391B (zh) * | 2016-01-26 | 2020-08-14 | 沈阳工业大学 | 基于scada运行数据提取的风电机组转矩控制性能评估方法 |
CN105758635A (zh) * | 2016-02-26 | 2016-07-13 | 内蒙古久和能源装备有限公司 | 风力发电机组传动链的转速震荡检测保护方法及装置 |
US10697431B2 (en) * | 2016-04-07 | 2020-06-30 | Vestas Wind Systems A/S | Control of a wind turbine taking noise into account |
CN105673357B (zh) * | 2016-04-14 | 2018-10-02 | 南京理工大学 | 一种适用于风力机模拟器的考虑时滞的转动惯量补偿方法 |
CN105863971A (zh) * | 2016-06-27 | 2016-08-17 | 国电联合动力技术有限公司 | 一种适用于风力发电机组塔筒避振的虚拟质量控制方法 |
CN106481507A (zh) * | 2016-10-28 | 2017-03-08 | 华北电力大学 | 考虑风湍流的风电机组模型预测控制器与控制方法 |
CN106545428A (zh) * | 2016-12-08 | 2017-03-29 | 中国北方发动机研究所(天津) | 一种军用柴油机工况预估的控制算法参数优化方法 |
CN106681151B (zh) * | 2017-01-20 | 2019-05-10 | 华北电力大学 | 基于时序分析的变桨距风机多模态运行优化控制装置及运行方法 |
CN106972804A (zh) * | 2017-03-22 | 2017-07-21 | 武汉工程大学 | 基于混合逻辑动态模型的无刷直流电机驱动系统建模方法 |
CN107882680B (zh) * | 2017-10-13 | 2019-10-25 | 上海致远绿色能源股份有限公司 | 一种用于风力发电机组的转速控制方法 |
JP6827992B2 (ja) * | 2017-10-16 | 2021-02-10 | 三菱重工業株式会社 | 風力発電装置及びその制御方法並びに制御プログラム |
CN108256693B (zh) * | 2018-02-11 | 2024-02-13 | 阳光电源股份有限公司 | 一种光伏发电功率预测方法、装置及系统 |
CN110210044A (zh) * | 2018-02-28 | 2019-09-06 | 北京金风科创风电设备有限公司 | 风力发电机组的载荷预测方法和装置 |
CN109030976B (zh) * | 2018-05-28 | 2020-11-24 | 上海大学 | 一种电力电子设备监测诊断系统和方法 |
CN108762089B (zh) * | 2018-06-15 | 2021-02-19 | 大连理工大学 | 一种基于模型预测的航空发动机在线优化及多变量控制设计方法 |
CN110875604B (zh) * | 2018-08-31 | 2021-05-11 | 北京金风科创风电设备有限公司 | 控制风力发电机组的方法及设备 |
CN109578205B (zh) * | 2018-11-27 | 2021-03-12 | 山东科技大学 | 一种风力机的智能快速刹车系统的控制方法 |
CN109340062B (zh) * | 2018-12-18 | 2020-01-31 | 国电联合动力技术有限公司 | 一种低风速风电机组数字双胞胎式的疲劳损伤预测方法 |
CN109667728B (zh) * | 2018-12-21 | 2020-09-08 | 北京金风科创风电设备有限公司 | 用于风力发电机组转速传感器的故障检测方法和装置 |
CN110259647B (zh) * | 2019-06-21 | 2020-05-22 | 江南大学 | 一种基于正多胞体滤波的风电机组缓变型故障诊断方法 |
CN110504880B (zh) * | 2019-07-24 | 2021-01-26 | 东南大学盐城新能源汽车研究院 | 一种磁通切换永磁直线电机干扰观测前馈补偿控制方法 |
CN110794367B (zh) * | 2019-10-12 | 2022-06-21 | 中国直升机设计研究所 | 一种桨涡干扰声源定位方法 |
CN111244968B (zh) * | 2020-02-05 | 2021-09-10 | 山东大学 | 计及电网电压支撑能力影响的风电场电压控制方法及系统 |
ES2875674B2 (es) * | 2020-05-07 | 2022-03-21 | Univ Del Pais Vasco / Euskal Herriko Unibertsitatea | Método de control predictivo basado en modelo para la reducción de carga estructural en turbinas eólicas, y turbina configurada para aplicar este método |
CN111817626A (zh) * | 2020-06-30 | 2020-10-23 | 浙江工业大学 | 一种风力发电机的分布式状态融合估计方法 |
CN112253406B (zh) * | 2020-09-29 | 2022-05-27 | 中国电建集团华东勘测设计研究院有限公司 | 一种海上风电机组的环境荷载预测方法与振动预控制系统 |
CN112682254B (zh) * | 2020-12-21 | 2022-09-16 | 北京华能新锐控制技术有限公司 | 一种基于动态多模型预测控制器的风机有功功率跟踪方法 |
CN112596489A (zh) * | 2020-12-29 | 2021-04-02 | 润电能源科学技术有限公司 | 工业控制系统的控制指令输出方法、装置及可读存储介质 |
CN117028145B (zh) * | 2023-10-08 | 2023-12-22 | 国网江苏省电力有限公司电力科学研究院 | 考虑塔架载荷抑制的风电机组有功功率控制方法及装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101560950A (zh) * | 2008-04-15 | 2009-10-21 | 西门子公司 | 基于预测控制风力涡轮机的方法和装置 |
CN102022264A (zh) * | 2009-08-31 | 2011-04-20 | 通用电气公司 | 用于风力涡轮机健康管理的系统和方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8021112B2 (en) * | 2010-06-29 | 2011-09-20 | General Electric Company | Methods and systems for monitoring operation of a wind turbine |
US8035242B2 (en) * | 2010-11-09 | 2011-10-11 | General Electric Company | Wind turbine farm and method of controlling at least one wind turbine |
-
2012
- 2012-02-29 CN CN2012100504184A patent/CN102588211B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101560950A (zh) * | 2008-04-15 | 2009-10-21 | 西门子公司 | 基于预测控制风力涡轮机的方法和装置 |
CN102022264A (zh) * | 2009-08-31 | 2011-04-20 | 通用电气公司 | 用于风力涡轮机健康管理的系统和方法 |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107869421A (zh) * | 2016-09-26 | 2018-04-03 | 北京金风科创风电设备有限公司 | 风力发电机变桨系统的控制方法和装置 |
CN107869421B (zh) * | 2016-09-26 | 2019-06-18 | 北京金风科创风电设备有限公司 | 风力发电机变桨系统的控制方法和装置 |
CN110145436A (zh) * | 2019-04-28 | 2019-08-20 | 华北电力大学 | 应用于风机的非线性经济模型预测控制方法 |
CN111305999A (zh) * | 2020-04-01 | 2020-06-19 | 宁夏三聚阳光新能源科技有限公司 | 一种智能调整风力发电系统 |
CN111305999B (zh) * | 2020-04-01 | 2020-12-01 | 宁夏三聚阳光新能源科技有限公司 | 一种智能调整风力发电系统 |
US12106022B2 (en) * | 2023-02-02 | 2024-10-01 | North China Electric Power University | Method and device for simulating dynamic digital twin model of dominant operation of wind turbine generator assembly |
Also Published As
Publication number | Publication date |
---|---|
CN102588211A (zh) | 2012-07-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102588211B (zh) | 一种风力发电机组全工况模型预测控制方法及系统 | |
Abdelbaky et al. | Design and implementation of partial offline fuzzy model-predictive pitch controller for large-scale wind-turbines | |
CN111801493B (zh) | 确定用于风力涡轮机的控制设置 | |
CN105179164B (zh) | 基于t-s模糊模型的风能转换系统滑模控制方法及装置 | |
CN110739893B (zh) | 一种改进自适应无轨迹卡尔曼滤波转动惯量辨识方法 | |
CN104932262B (zh) | 基于mpc和pi控制方法的带电动汽车的微网调频方法 | |
CN104595106B (zh) | 基于强化学习补偿的风力发电变桨距控制方法 | |
Sami et al. | Sensorless fractional order composite sliding mode control design for wind generation system | |
CN103758699B (zh) | 一种风力发电机组的桨距角控制方法及桨距角控制器 | |
CN105545595A (zh) | 基于径向基神经网络的风力机反馈线性化功率控制方法 | |
CN106481507A (zh) | 考虑风湍流的风电机组模型预测控制器与控制方法 | |
CN110397548A (zh) | 一种双馈式变速抽水蓄能机组的多模型预测控制方法 | |
CN105649877A (zh) | 一种大型风电机组的蚁群pid独立变桨控制方法 | |
CN105156270A (zh) | 风力发电机独立变桨控制系统及方法 | |
CN105626378A (zh) | 基于rbf神经网络的大型风电机组独立变桨控制方法 | |
CN115167140B (zh) | 风力发电机组多目标随机模型预测控制策略方法及系统 | |
Imran et al. | LQG controller design for pitch regulated variable speed wind turbine | |
WO2018193802A1 (ja) | 風力発電システム | |
Chakib et al. | Comparative study of active disturbance rejection control with RST control for variable wind speed turbine based on doubly fed induction generator connected to the grid | |
CN107979112A (zh) | 一种风机控制方法、系统、终端及可读存储介质 | |
CN116845886A (zh) | 一种基于模型预测的多端口自主光伏系统构网型控制方法 | |
Owen et al. | The augmented unscented H-infinity transform with H-infinity filtering for effective wind speed estimation in wind turbines | |
Shahravanmehr et al. | LQG controller based on fuzzy logic to control the power of wind turbine | |
Yao et al. | LQG controller for a variable speed pitch regulated wind turbine | |
Schlipf | Prospects of multivariable feedforward control of wind turbines using lidar |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |