一种可抑制多扰动因素的风机模糊自适应变桨距控制方法与流程

1.本发明涉及风力机组变桨距控制技术领域，具体的说是涉及一种可抑制多扰动因素的风机模糊自适应变桨距控制方法。


背景技术：

2.经济社会的持续快速发展，离不开有力的能源保障。风能作为一种绿色可再生能源，分布广泛，发电成本相对较低，而且安全性、可靠性较高。正是因为这些独特的优点，风力发电具有良好的社会效益和经济效益，因此风力发电技术越来越受到国家的重视。
3.由于风力机组在实际运行中不确定因素较多，传统风机的桨距控制策略仅考虑单一的风速扰动，并未综合考虑安装与制造时造成的初始参数误差、风力机组运行期间的环境变化(例如湍流风速、积雪、积霜等)以及风力机组叶片变化(例如叶片的弯曲、变形等)等多因素扰动，因此影响风力机的输出功率以及电能质量，另外还同时造成桨叶上的载荷不平衡，引起桨叶的拍打振动和扭矩波动。


技术实现要素：

4.本发明的目的在于提出一种可抑制多扰动因素的风机模糊自适应变桨距控制方法，建立计及多因素扰动的风电机组变桨控制系统闭环传递函数，设计一种基于模糊自适应pid算法的风力发电机组变桨距控制方法，实现对多种因素引起的扰动自适应快速响应，确保风电机组桨距角得到恰当调整，维持系统稳定，实现对风能的最优捕获。
5.为达到上述目的，本发明采用的技术方案如下：
6.本发明提供一种可抑制多扰动因素的风机模糊自适应变桨距控制方法，包括：
7.建立考虑多扰动因素的风力机组变桨距控制系统的闭环传递函数；
8.基于所述闭环传递函数和模糊规则确定模糊自适应pid控制参数；
9.采用模糊自适应pid控制对所述考虑多扰动因素的风力机组变桨距控制系统进行调节，输出风力机组桨距角。
10.进一步的，所述扰动因素包括风力机组运行期间叶片变化扰动，风力机组安装初始参数误差扰动和风力机组运行期间环境变化扰动；
11.建立考虑多扰动因素的风力机组变桨距控制系统的闭环传递函数，包括：
12.建立风力机组运行期间叶片变化扰动的传递函数：
[0013][0014]
其中，为叶根拍打弯矩m对桨距角β的线性化系数，为叶片的空气动力阻尼系数，ω
z
为叶片拍打振动模态的自然角频率，v(s)为外界风速；
[0015]
建立风力机组安装初始参数误差扰动的传递函数：
[0016][0017]
其中，k1，k2为扰动系数；
[0018]
建立风力机组运行期间环境变化扰动的传递函数：
[0019][0020]
其中，k3，k4，k5为扰动系数；
[0021]
综合上述传递函数，得到风力机组变桨距控制系统的闭环传递函数：
[0022][0023]
其中，d(s)＝d1(s) d2(s) d3(s)，k
p
、k
i
、k
d
为pid控制器的自适应参数，g
θ
为变桨控制器传递函数，g
c
为变桨执行器传递函数，t
β
为时间常数，g
p
为风力机的传递函数，j为转动惯量，δv(s)＝v(s)
‑
v0(s)，δβ(s)＝β(s)
‑
β0(s)，v0为初始风速，β(s)为桨距角，β0是风力机工作点处的浆距角，μ、ξ和γ均为线性系数。
[0024]
进一步的，基于所述闭环传递函数和模糊规则确定模糊自适应pid控制参数，包括：
[0025]
对所述闭环传递函数进行降阶处理，并基于降阶处理后的闭环传递函数计算风力机组变桨距控制系统的超调量和调节时间；
[0026]
根据超调量和调节时间整定模糊自适应pid控制的参数初值；
[0027]
基于模糊规则确定模糊自适应pid控制参数的调整量；
[0028]
将模糊自适应pid控制器参数的调整量与整定的参数初值相结合得到模糊自适应pi控制参数。
[0029]
进一步的，
[0030]
降阶后的闭环传递函数为：
[0031][0032]
其中，φ(s)为降阶后的闭环传递函数，ζ为阻尼比，ω
n
为自然频率，ζ和ω
n
根据主导极点得到，
[0033]
所述主导极点确定如下：
[0034]
分别令传递函数g(s)的分子、分母等于零得到零点z
i
和极点p
i
，在复平面内寻找出距离虚轴最近且附近无闭环零点的极点，即为主导极点。
[0035]
进一步的，所述计算风力机组变桨距控制系统的超调量和调节时间如下：
[0036]
当0＜ζ＜1时，超调量σ为：
[0037]
[0038]
以误差带δ＝0.05为标准，调节时间t
s
为：
[0039][0040]
当ζ＞1时，不存在超调量；以误差带δ＝0.05为标准，调节时间t
s
为：
[0041][0042][0043][0044]
进一步的，所述基于模糊规则确定模糊自适应pid控制参数的调整量，包括：
[0045]
将风力机组的转速理论值ω
ref
与实际值间的差值e以及误差变化量ec作为模糊控制器的输入量；
[0046]
将输入量模糊化，根据预设定的模糊规则运用mamdani模糊推理法得到模糊子集；
[0047]
将模糊子集解模糊得到模糊控制器的输出量作为pid控制器自适应参数的调整量δk
p
、δk
i
、δk
d
。
[0048]
进一步的，所述模糊控制器采用二维模糊控制器。
[0049]
进一步的，所述模糊控制器采用三角形隶属度函数。
[0050]
进一步的，采用重心法对模糊子集解模糊得到模糊控制器的输出量。
[0051]
本发明具有如下有益效果：
[0052]
本发明考虑到风力机组控制过程中各种扰动因素，建立计及多因素扰动的风电机组变桨控制系统闭环传递函数，确定模糊控制规则，实现对多种因素引起的扰动自适应快速响应，确保风电机组桨距角得到恰当调整，维持系统稳定，实现对风能的最优捕获。
[0053]
本发明将模糊控制应用在变桨距控制系统中，模糊控制具有鲁棒性和自适应能力强以及不需要建立精确的系统动态模型就能够实现很好的控制效果。
附图说明
[0054]
图1为本发明的一种可抑制多扰动因素的风机模糊自适应变桨距控制方法原理图；
[0055]
图2为本发明预设的模糊控制规则表。
具体实施方式
[0056]
下面对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制本发明的保护范围。
[0057]
风力机组在实际运行中不确定因素较多，包括风电机组在安装与制造时的初始参数误差、风力机组运行期间的环境变化(例如湍流风速、积雪、积霜等)，以及叶片变化(例如叶片的弯曲、变形等)等因素，本发明通过建立计及多因素扰动的风电机组变桨控制系统闭环传递函数，设计一种基于模糊自适应pid算法的风力发电机组变桨距控制策略，实现对多
种因素引起的扰动自适应快速响应，确保风电机组桨距角得到恰当调整，维持系统稳定，实现对风能的最优捕获。
[0058]
一种可抑制多扰动因素的风机模糊自适应变桨距控制方法，包括：
[0059]
建立考虑多扰动因素的风力机组变桨距控制系统的闭环传递函数；所述扰动因素包括风力机组运行期间叶片变化扰动，风力机组安装初始参数误差扰动和风力机组运行期间环境变化扰动；
[0060]
基于模糊自适应pid控制对考虑多扰动因素的风力机组变桨距控制系统进行调节，输出风力机组桨距角。
[0061]
作为一种优选的实施方式，本实施例中，建立考虑多扰动因素的风力机组变桨距控制系统的闭环传递函数参见图1，具体如下：
[0062]
建立风力机组运行期间叶片变化扰动的传递函数，
[0063]
叶片变化因素包括叶片的弯曲、变形等，可将该扰动用传递函数d1(s)近似表述为：
[0064][0065]
式中：为叶根拍打弯矩m对桨距角β的线性化系数，为叶片的空气动力阻尼系数，ω
z
为叶片拍打振动模态的自然角频率，和ω
z
均为已知量，v(s)为外界风速，s为传递函数参数。
[0066]
建立风力机组安装初始参数误差扰动的传递函数，
[0067]
风力机组在安装与制造时存在初始参数误差，可将该扰动用传递函数d2(s)近似表述为：
[0068][0069]
式中：k1，k2为扰动系数，其值根据风力机安装与制造时的误差取值，为已知量，v(s)为外界风速。
[0070]
建立风力机组运行期间环境变化扰动的传递函数，
[0071]
风力机组运行期间的环境变化包括湍流风速、积雪、积霜等，可将该扰动用传递函数d3(s)近似表述为：
[0072][0073]
式中：k3，k4，k5为扰动系数，其值根据运行时环境的变化取值，为已知量；v(s)为外界风速。
[0074]
综合以上各扰动因素，建立风力机组变桨距控制系统的闭环传递函数如下所示：
[0075][0076]
式中，d(s)＝d1(s) d2(s) d3(s)，k
p
、k
i
、k
d
为pid控制器的自适应参数，g
θ
、g
c
、g
p
均为组成风机系统各部分的传递函数，g
θ
为变桨控制器传递函数，g
c
为变
桨执行器传递函数，t
β
为时间常数，单位为s；g
p
为风力机的传递函数，j为转动惯量，δv(s)＝v(s)
‑
v0(s)，δβ(s)＝β(s)
‑
β0(s)，v0为初始风速，β(s)为桨距角，β0是风力机工作点处的浆距角，μ、ξ和γ均为线性系数。
[0077]
作为一种优选的实施方式，本实施例中，通过令传递函数g(s)的分子、分母等于零以求解出系统零、极点z
i
、p
i
(i＝1,2,3
…
)，p
i
是关于k
p
、k
i
、k
d
的量。
[0078]
在复平面内，寻找出距离虚轴最近且附近无闭环零点的极点，这些极点就是主导极点。从而进行等效简化实现对计及多因素扰动控制系统闭环传递函数的降阶，对于降阶后得到的二阶系统有如下表达式：
[0079][0080]
式中，ζ为阻尼比，ω
n
为自然频率，ζ和ω
n
需根据主导极点得到，因此ζ和ω
n
是关于k
p
、k
i
、k
d
的量。当0＜ζ＜1时，系统超调量σ为：
[0081][0082]
以误差带δ＝0.05为标准，调节时间t
s
为：
[0083][0084]
当ζ＞1时,不存在超调量；以误差带δ＝0.05为标准，调节时间t
s
为：
[0085][0086]
上式中，
[0087]
系统超调量和调节时间越小越好，因为ζ和ω
n
是关于k
p
、k
i
、k
d
的量，所以σ和t
s
也是关于k
p
、k
i
、k
d
的量，因此可根据超调量σ和调节时间t
s
的计算公式指导整定模糊自适应pid控制器的参数初值k
p0
、k
i0
、k
d0
使超调量和调节时间尽可能小。
[0088]
作为一种优选的实施方式，本实施例中，采用模糊控制器基于模糊规则确定pid控制器参数的调整量，将模糊控制器输出的pid控制器参数的调整量与整定的pid控制器初始参数结合得到pi控制器的自适应参数k
p
、k
i
、k
d
。
[0089]
参见图1，模糊控制器采用二维模糊控制器，将风力机组的转速理论值ω
ref
与实际值间的差值e及其误差变化量ec作为模糊控制器的输入量，ω
ref
为风力机转速给定值，一般设定为风机转速额定值，然后将输入量模糊化，再根据相应的模糊规则运用mamdani模糊推理法得到模糊子集，模糊控制规则参见图2，最后将模糊子集解模糊得到模糊控制器的输出量，该输出量为pid控制器参数的调整量δk
p
、δk
i
、δk
d
。
[0090]
本实施例中，模糊控制器采用三角形隶属度函数，并根据风机调节要求，确定模糊控制规则。
[0091]
本实施例中，将模糊控制器输入量转换为相应的模糊语言变量值以实现模糊化。
[0092]
本实施例中，采用重心法对模糊子集解模糊得到模糊控制器的输出量。
[0093]
综合模糊控制器输出量与pid控制器初始参数可以得到pi控制器的自适应参数k
p
、k
i
、k
d
，计算公式如下：
[0094]
k
p
＝k
p0
δk
p
[0095]
k
i
＝k
i0
δk
i
[0096]
k
d
＝k
d0
δk
d
[0097]
参见图1，ω
ref
和β
ref
分别为风力机转速的给定值和浆距角的给定值。基于得到pid控制器的自适应参数k
p
、k
i
、k
d
，采用pid控制器对考虑多扰动因素的风力机组变桨距控制系统进行调节，输出风力机组桨距角。
[0098]
本发明所提风机变桨距控制方法综合考虑到风力机组控制过程中各种扰动因素，实现了对多种因素引起的扰动自适应快速响应，确保风力机组桨距角得到恰当调整，维持系统稳定，实现对风能的最优捕获。
[0099]
本领域内的技术人员应明白，本技术的实施例可提供为方法、系统、或计算机程序产品。因此，本技术可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本技术可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd
‑
rom、光学存储器等)上实施的计算机程序产品的形式。
[0100]
本发明是参照根据本技术实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0101]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0102]
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0103]
最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解：依然可以对本发明的具体实施方式进行修改或者等同替换，而未脱离本发明精神和范围的任何修改或者等同替换，其均应涵盖在本发明的权利要求保护范围之内。