一种风力发电机组独立变桨控制方法和系统与流程

1.本发明属于风力发电领域，涉及一种风力发电机组独立变桨控制方法和系统。


背景技术：

2.随着风力发电机组设计开发向大型化发展，叶片长度及塔架高度在不断增加，在不同风况下各大部件的受力情况也变得复杂，导致极限载荷水平也在不断增加。通过先进的控制策略架构及使用新的传感器引入控制策略可以有效地解决面对复杂风况带来的极限载荷问题。
3.现有针对此问题的技术方案包括：在叶片根部安装载荷传感器，实时检测叶根载荷。通过坐标转换成d轴载荷分量及q轴载荷分量，使用两个pid控制器分别以两个载荷分量作为控制输入，输出三个不同的变桨角度给三个变桨执行机构。但由于两个pid控制回路会出现耦合现象，且控制系统输入噪声与测量噪声无法忽视。同时设置控制增益只能单向调节，不能兼顾其控制成本及控制效果。


技术实现要素：

4.本发明的目的在于克服上述现有技术中，传统单输入单输出的方式对三只叶片的叶根载荷进行分解成两个轴，分别进行控制，但多控制回路的方式易造成耦合，且控制增益未能兼顾控制成本及控制偏差两个相矛盾的指标的缺点，提供一种风力发电机组独立变桨控制方法和系统。
5.为了达到上述目的，本发明采用以下技术方案予以实现：
6.一种风力发电机组独立变桨控制方法，包括如下步骤：
7.步骤1)设定检测周期，在检测周期内检测风电机组叶片的叶根载荷及叶轮方位角，建立旋转坐标系；
8.步骤2)基于风电机组的叶片模态、发电机模态和塔架模态，建立线性化模型；对旋转坐标系下的叶根载荷进行坐标变换，得到非旋转坐标系下的分量，基于非旋转坐标系下的分量对线性化模型进行状态估计；
9.步骤3)通过估计后的风电机组状态计算全状态负反馈的控制增益，进一步计算得到非旋转坐标系下的控制指令；
10.步骤4)将非旋转坐标系下的控制指令通过坐标反变换计算旋转坐标系下的独立变桨角度指令；
11.步骤5)在检测周期内检测发电机转速，计算统一的变桨角度指令，结合独立变桨角度指令，计算得到三个叶片变桨角度指令，基于三个叶片变桨角度指令进行变桨动作。
12.优选地，步骤2)中，基于非旋转坐标系下的分量，结合线性化模型进行状态华估计，采用卡尔曼滤波的方式计算，使用ricatti方程计算卡尔曼滤波器的增益。
13.优选地，步骤3)中，计算全状态负反馈的控制增益具体为：
14.引入代价方程对控制成本进行均衡，之后使用ricatti方程来计算全状态负反馈
的控制增益。
15.优选地，步骤4)中，还包括对独立变桨角度指令进行限幅。
16.优选地，步骤5)中，计算统一的变桨角度指令过程中，以控制发电机转速为目标进行转速
‑
变桨回路的pid控制。
17.优选地，步骤5)中，三个叶片变桨角度指令的具体计算过程为：
18.将限幅后的独立变桨角度指令叠加统一的变桨角度指令，得到三个叶片变桨角度指令，将三个叶片变桨角度指令传递给变桨执行机构进行变桨动作，以进行独立变桨降载及发电机转速控制。
19.一种风力发电机组独立变桨控制系统，包括：
20.旋转坐标系建立单元，用于获取检测周期内风电机组叶片的叶根载荷及叶轮方位角，建立旋转坐标系；
21.状态估计单元，与旋转坐标系建立单元相交互，基于风电机组的叶片模态、发电机模态和塔架模态，建立线性化模型；对旋转坐标系下的叶根载荷进行坐标变换，得到非旋转坐标系下的分量，基于非旋转坐标系下的分量对线性化模型进行状态估计；
22.全状态负反馈处理单元，与状态估计单元相交互，用于计算全状态负反馈的控制增益，进一步计算得到非旋转坐标系下的控制指令；
23.变桨动作控制单元，与全状态负反馈处理单元相交互，在检测周期内检测发电机转速，计算统一的变桨角度指令，进一步计算得到三个叶片变桨角度指令，基于三个叶片变桨角度指令进行变桨动作。
24.优选地，变桨动作控制单元包括独立变桨角度指令模块，用于将非旋转坐标系下的控制指令通过坐标反变换计算旋转坐标系下的独立变桨角度指令；三个叶片变桨角度指令由统一的变桨角度指令和独立变桨角度指令结合计算得到。
25.与现有技术相比，本发明具有以下有益效果：
26.本发明提供了一种风力发电机组独立变桨控制方法，通过检测三支叶片在旋转坐标系下的载荷及叶轮方位角，对风力发电机组进行状态空间建模，对其进行简化，在简化模型的基础上对风力发电机组的状态进行有效估计，将估计的状态进行全状态反馈并使用线性二次调节器对控制增益进行计算，最大程度避免多个单输入单输出回路之间耦合的影响，并兼顾控制成本与控制效果，在复杂风况下通过独立变桨提高机组在极端风况下的运行安全。本发明对不同叶轮方位角下的叶片模态塔架模态进行了线性化建模，克服了传统方式中采用单输入单输出多回路的方式，并创新地采用状态估计并使用线性二次调节器进行全状态反馈的方式来进行独立变桨控制。
27.本发明还公开了一种风力发电机组独立变桨控制系统，替代已有独立变桨的方案，从而避免由于多回路之间耦合，未考虑输入噪声及测量噪声给控制效果带来的影响。同时依据已采集的叶根载荷信号，叶轮方位角信号，塔架位移信号等可估计机组线性化模型，依据此模型可较为精确的估计机组状态信息从而进行更为精准的反馈闭环控制。
附图说明
28.图1为本发明基于线性二次调节器的风力发电机组独立变桨控制方法的流程图；
29.图2为线性二次调节器 卡尔曼滤波器的控制架构图。
具体实施方式
30.下面结合附图对本发明做进一步详细描述：
31.实施例1
32.一种风力发电机组独立变桨控制方法，包括如下步骤：
33.步骤1)设定检测周期，在检测周期内检测风电机组叶片的叶根载荷及叶轮方位角，建立旋转坐标系；
34.步骤2)基于风电机组的叶片模态、发电机模态和塔架模态，建立线性化模型；对旋转坐标系下的叶根载荷进行坐标变换，得到非旋转坐标系下的分量，基于非旋转坐标系下的分量对线性化模型进行状态估计；
35.步骤3)通过估计后的风电机组状态计算全状态负反馈的控制增益，进一步计算得到非旋转坐标系下的控制指令；
36.步骤4)将非旋转坐标系下的控制指令通过坐标反变换计算旋转坐标系下的独立变桨角度指令；
37.步骤5)在检测周期内检测发电机转速，计算统一的变桨角度指令，结合独立变桨角度指令，计算得到三个叶片变桨角度指令，基于三个叶片变桨角度指令进行变桨动作。
38.实施例2
39.一种风力发电机组独立变桨控制方法，包括：
40.在当前的检测周期内检测三支叶片的叶根载荷及叶轮方位角，对整机发电机组进行线性化建模，并同时将旋转坐标系下的叶根载荷进行坐标变换将其转换成非旋转坐标系下的分量，由此分量结合已搭建的线性化模型进行状态估计，此状态估计采用卡尔曼滤波的方式计算，使用ricatti方程计算卡尔曼滤波器的增益。
41.在得到风力发电机组估计的状态后使用全状态负反馈的方式，并引入代价方程对控制成本及控制效果进行均衡，同样使用ricatti方程来计算全状态负反馈的控制增益来计算在非旋转坐标系下的控制指令。
42.将非旋转坐标系下的控制指令通过坐标反变换计算旋转坐标系下的独立变桨角度指令。对此独立变桨角度指令进行限幅。
43.在当前的检测周期内检测发电机转速，以控制发电机转速为目标进行转速
‑
变桨回路的pid控制，计算统一的变桨角度指令。
44.将限幅后的独立变桨角度指令叠加统一的变桨角度指令得到最终的三个叶片变桨角度指令传递给变桨执行机构进行变桨动作，以进行独立变桨降载及发电机转速控制。
45.实施例3
46.一种风力发电机组独立变桨控制方法，如图1所示，包括如下步骤：
47.1)风力发电机组的非线性气动运动方程可以描述为其中m为质量矩阵，f为非线性运动方程，分别为风力发电机组自由度变量的位移，速度及加速度，u为控制输入向量，u
d
为风输入干扰向量，t为时间。
48.2)通过对在运行点上进行微小的摄动(δ)，记为其中op记为运行点。
49.3)对第一步的公式以第二步在运行点上的微小摄动进行泰勒展开成二阶线性化
方程，并将以替代从而得到状态空间方程其中a
r
为状态矩阵，b
r
为输入矩阵，b
dr
为风输入干扰矩阵，d
dr
为风输入转换矩阵，x
r
为状态变量，u
r
为输入变量，y
r
为输出变量，w
r
为风输入干扰变量，v
r
为测量噪声，其中下标r代表是旋转坐标系，意味着其为随时间变化的。
50.4)由于上述状态空间方程为随时间变化，为应用线性时不变模型需要将其转化为非旋转坐标系。
51.5)第b个叶片叶轮方位角可表示为：ψ
b
＝ψ (b
‑
1)*(2*pi/3)。其中b为叶片个数，ψ为叶轮方位角。定义为ψ＝0为第一个叶片竖直朝上。
52.6)由于非旋转坐标系为旋转坐标系的投影，三个叶片的旋转坐标系可转换成两个轴的非旋转坐标系其关系为其中q0为统一模态，q
c
为cos循环模态，q
s
为sin循环模态。将此标记为转换矩阵，将第一步的状态空间方程中的x
r
＝t
s
(ψ)x
nr
，u
r
＝t
c
(ψ)u
nr
，y
r
＝t
o
(ψ)y
nr
，可以得到转换后的状态空间方程为其中nr代表非旋转坐标系。
53.7)由于其已经转换为弱叶轮方位角相关，因此将其对每个叶轮方位角进行平均可以在不缺失重要模态信息下得到线性时不变的状态空间方程记为
54.8)本实施例中，q代表的风力发电机组自由度可选择为叶片1，2，3的挥舞方向位移，叶轮方位角及塔架一阶前后位移。u输入变量选为三个叶片桨距角β1，β2，β3。
55.9)为解决多控制回路互相耦合的问题引入多输入多输出控制框架中的lqg算法。第7步得到的状态空间方程中w和v记为风输入干扰噪声及测量噪声且是零均值高斯白噪声，两者互相没有关联。
56.10)lqg最优控制运用最少的控制花费来获得最少的控制偏差，因此引入代价方程其中q代表状态权衡矩阵，r代表输入权衡矩阵，且两矩阵定义为半正定，即q＝q
t
≥0和r＝r
t
≥0。此两矩阵分别代表权衡快速控制及低控制成本的控制矩阵，以此代价方程最小为控制目标。
57.11)由第10步的代价方程引出需要确定状态变量与输入变量之间的关系k
c
，使用全状态负反馈关系也即线性二次调节器代表为u＝
‑
k
c
x。其中k
c
＝r
‑1b
t
p
c
，其中p
c
也为半正定矩阵即p
c
＝p
ct
≥0，p
c
由计算ricatti方程来得出，公式为a
t
p
c
p
c
a
‑
p
c
br
‑1b
t
p
c
m
t
qm＝0。
58.12)由于第11步需要全状态反馈，状态变量由卡尔曼滤波器来估计风力发电机组
的状态记为其卡尔曼增益由k
f
＝p
f
c
t
v
‑1，其目标为减小，其目标为减小其中p
f
也为半正定矩阵即p
f
＝p
ft
≥0，p
c
由计算ricatti方程来得出，公式为p
f
a
t
ap
f
‑
p
f
c
t
v
‑1cp
f
b
d
wb
dt
＝0。
59.13)可由第7步线性化得到的模型估计风力发电机组的状态，其公式为将线性二次调节器u＝
‑
k
c
x及y＝cx du d
a
v代替可得整个的lqg控制框架为其中，a为系统状态系数矩阵，b为系统控制系数矩阵，c为输出状态系数矩阵，d为输出控制系数矩阵，kf为状态估计增益，kc为全状态反馈增益，x为系统状态向量，y为系统输出向量，w为噪声输入向量。
60.14)通过线性二次调节器u＝
‑
k
c
x计算出来的控制指令依然为非旋转坐标系下，因此需要对非旋转坐标系下的控制指令通过坐标反变换得到旋转坐标系下的独立变桨角度指令。
61.15)将三个不同的独立变桨角度指令进行限幅后叠加相同的统一变桨角度指令得到最终的变桨角度指令并发送给变桨执行机构进行变桨动作来完成独立变桨控制。
62.结合图2，y代表风力发电机组的输出，u代表风力发电机组的控制输入，其由控制增益
‑
kc与估计的系统状态向量乘积来取得。此为全状态负反馈。通过a、b、c、d及kf结合第13步方程可得到估计的系统状态向量当得到估计的系统状态向量及控制增益
‑
kc可得到整个控制算法的控制输入u。
63.实施例4
64.一种风力发电机组独立变桨控制系统，包括：
65.线性化模型建立单元，用于获取检测周期内风电机组叶片的叶根载荷及叶轮方位角，建立线性化模型；
66.状态估计单元，对旋转坐标系下的叶根载荷进行坐标变换，得到非旋转坐标系下的分量，对线性化模型进行状态估计；
67.全状态负反馈处理单元，与状态估计单元相交互，用于计算全状态负反馈的控制增益，进一步计算得到非旋转坐标系下的控制指令；
68.变桨动作控制单元，在检测周期内检测发电机转速，计算统一的变桨角度指令，结合独立变桨角度指令，计算得到三个叶片变桨角度指令，基于三个叶片变桨角度指令进行变桨动作。
69.本发明创新地引入了多输入多输出的方式，对非线性强耦合时变的风力发电机组模型进行简化建模，通过坐标转换的方式将旋转坐标系下的模型转换成非旋转坐标系，并利用线性化的模型对机组内部的状态进行有效的估计，并对估计的状态进行全状态反馈，其控制增益采用代价方程的方式计算得出，可以有效解决控制成本与控制偏差相矛盾的问题，也即在能够解决降载问题的同时兼顾变桨执行器的寿命。限幅的操作，例如当计算的独立变桨指令超过统一变桨指令的50％以上，则最大输出指令值即为统一变桨指令的50％。
70.需要说明的是，由于三支叶片在不断旋转，测量到的叶根载荷是随叶轮方位角不
断变化，此旋转坐标系代表了随叶轮方位角变化的随时间变化的坐标系。对包含叶片模态，发电机模态，塔架模态等运动方程进行建模并线性化简化后，由于运动方程中已包含相关模态信息，因此可以对状态进行估计。计算统一的变桨角度指令具体是：采用pid或pi控制的方式，通过变桨来保持发电机转速在额定转速，因此控制输入为发电机转速，控制算法为pid或pi，控制输出为变桨角度指令(三支叶片一致)。
71.本发明在由旋转坐标系转换成非旋转坐标系过程中采用叶轮方位角平均的方式，对叶轮方位角划分越多则线性化过程计算量越大，但不限于此方式，其他如未采用叶轮方位角平均即不考虑360度内不同叶轮方位角下状态空间方程的不同。考虑不平衡载荷对三支不同的叶片变桨角度进行控制的方式，因此线性化模型仅包含变桨控制，但不限于此方式，其他如转矩控制的方式也可包含在内。在线性化过程中仅包含了叶片挥舞方向模态，塔架一阶模态，发电机旋转状态，因此状态变量为10个，但不限于此方式，其他如包含更高阶的叶片模态及塔架模态，采用更多的状态变量以使线性化模型包含更多的模态信息。在使用线性二次调节器过程中采用的全状态负反馈的方式，但不限于此方式，其他如部分状态反馈，正反馈等方式。在传统转速变桨单输入单输出控制回路给出的统一变桨指令叠加多输入多输出控制器给出独立变桨的方式，但不限于此方式，其他如仅采用多输入多输出的方式同时对发电机转速及叶根载荷进行控制。采用的是使用三支旋转叶片叶根载荷的方式，但不限于此方式，其他如直接采集非旋转坐标系下的轴承载荷，塔架载荷等方式。
72.综上所述，本发明提供了一种基于线性二次调节器的风力发电机组独立变桨控制的方法，替代已有独立变桨的方案，从而避免由于多回路之间耦合，未考虑输入噪声及测量噪声给控制效果带来的影响。同时依据已采集的叶根载荷信号，叶轮方位角信号，塔架位移信号等可估计机组线性化模型，依据此模型可较为精确的估计机组状态信息从而进行更为精准的反馈闭环控制。本发明设计了一种基于线性二次调节器的风力发电机组独立变桨控制的方法，从而依据已测量的叶根载荷，结合卡尔曼滤波的状态估计及线性二次调节的闭环反馈更为精准的进行独立变桨控制。本发明创新地采用了卡尔曼滤波进行状态估计并依据此运用线性二次调节器进行状态控制，在原有独立变桨控制基础上，结合多输入多输出控制框架，避免多回路耦合、控制成本与效果不兼顾、未考虑系统输入及测量噪声的问题。
73.以上内容仅为说明本发明的技术思想，不能以此限定本发明的保护范围，凡是按照本发明提出的技术思想，在技术方案基础上所做的任何改动，均落入本发明权利要求书的保护范围之内。