CN108021026B - 一种风力发电机组故障预警及控制参数在线优化方法 - Google Patents
一种风力发电机组故障预警及控制参数在线优化方法 Download PDFInfo
- Publication number
- CN108021026B CN108021026B CN201711101416.2A CN201711101416A CN108021026B CN 108021026 B CN108021026 B CN 108021026B CN 201711101416 A CN201711101416 A CN 201711101416A CN 108021026 B CN108021026 B CN 108021026B
- Authority
- CN
- China
- Prior art keywords
- generating set
- wind generating
- early warning
- fault early
- parameters
- 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
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Wind Motors (AREA)
Abstract
本发明公开了一种风力发电机组故障预警及控制参数在线优化方法,该方法是根据机组的振动、转速、温度、桨叶角度、变桨速率等易测的时域信号,采用频谱分析方法提取出机组各大部件的特征频率,对风力发电机组各个部件的运行状态进行监测。风力发电机组各部件正常运作时,其特征频率在一定阈值范围内;当某部件的特征频率超出阈值,则说明该部件损坏或出现异常。此时,机组报出相应的故障预警,并对塔架及传动链加阻、变桨与转矩滤波器和PID参数等主控参数进行在线优化调整,防微杜渐,保证机组的可靠运行,具有实际推广价值,值得推广。
Description
技术领域
本发明涉及风力发电机组故障预警及控制参数在线优化的技术领域,尤其是指一种风力发电机组故障预警及控制参数在线优化方法。
背景技术
随着风电行业的高速发展,风力发电机组装机规模不断增长,运行年限的增加,机组关键大部件(如传动系统、偏航系统、变桨系统等)的故障频出,关键部件故障的发生导致机组的运行和维护成本也越来越高,对风电设备制造商及风电开发商造成了重大经济损失。因此,提前预警故障,保证风力发电机组控制参数的准确性,降低机组大故障率及运维成本,成为整机制造商及风电开发商控制经济成本的最有效的途径。
目前,调整风力发电机组控制参数的方法有:
1)经典控制理论中的工程设计方法。该方法通过串联校正和反馈校正直接计算得到校正成典型系统后的控制器参数理论值,但该方法计算过程复杂,需要对模型进行近似化处理,且校正后还需进行大量的调整试凑,因此优化效果并不明显。
2)小信号建模分析法。该方法利用小信号模型在忽略控制器的条件下,对系统进行特征根分析,使用根轨迹分析控制器参数对系统性能的影响。虽然该方法物理概念清晰。但是小信号分析是在系统确定的运行点进行线性化处理并在此基础上进行零极点配置的,而当系统的运行状态改变或发生大干扰故障时,优化后的控制参数可能会使系统性能变差,需要重新进行优化等。
发明内容
本发明的目的在于克服现有技术的不足,提出了一种简单易行、收敛速度快的风力发电机组故障预警及控制参数在线优化方法,当机组需报出相应的故障预警后,对塔架及传动链加阻、变桨与转矩滤波器和PID参数等主控参数进行在线优化调整,防微杜渐,以保证机组的可靠运行。
为实现上述目的,本发明所提供的技术方案为:一种风力发电机组故障预警及控制参数在线优化方法,包括以下步骤:
1)故障预警
首先,实时采集风力发电机组的运行数据,包括振动、转速、温度、桨叶角度、变桨速率;然后对风力发电机组的关键部件进行频谱特性分析,最后再对风力发电机组各个部件的特征频率及对应幅值运行状态进行监测;其中,当风力发电机组各部件正常运作时,说明其特征频率不偏移且对应的幅值在预设阈值范围内;若某部件的特征频率发生大范围变化或对应的幅值超出理论阈值,则说明该部件损坏或出现异常,此时,机组报出相应的故障,为发生损坏的零件进行预警;
2)控制参数在线优化
当风力发电机组的某部件出现故障预警时,说明此刻的系统控制参数不利于当前机组的稳定运行,需要进行参数调整;此时,采用一种改进的粒子群参数优化算法,对控制参数进行在线调整优化,以满足风力发电机组的实际运行要求,保证机组的可靠运行;其中,改进的粒子群参数优化算法包括以下步骤:
2.1)假设在D维的目标空间中,有m个粒子组成一个群落,初始化第i个粒子的位置为Xi=(xi1,xi2,…,xiD),第i个粒子的速度为Vi=(vi1,vi2,…,viD),i=1,2,…,m;
2.2)输入第i个粒子参数,实时存储机组运行的数据;
2.3)若存储的数据量达标,则对存储的数据进行频谱分析,获得特征频率的能量幅值;若存储的数据量不达标,则返回步骤2.2)继续执行;
2.4)更新下一组粒子参数,并判断是否达到m个粒子;若达到m个粒子,则进入步骤2.5);否则,返回步骤2.2)继续执行;
2.5)根据频率幅值,找出局部最优解Pbesti和全局最优解Gbest;
2.6)根据下列公式不断更新粒子的速度及位置:
式中,t表示时刻,d表示目标空间维度d=1,2,…,D,vid(t)表示t时刻第i个粒子在d维的更新速度,vid(t+1)表示t+1时刻第i个粒子在d维的更新速度,xid(t)表示t时刻第i个粒子在d维的位置,xid(t+1)表示t+1时刻第i个粒子在d维的位置,c1和c2表示学习因子,rand()表示随机数,Pbesti表示第i个粒子的局部最优解,Gbest表示全局最优解;
2.7)若达到最大迭代次数或者小于或等于频率幅值的正常阈值,则输出最优参数;否则,返回步骤2.2)继续执行,直到符合终止条件。
在步骤1)中,对风力发电机组的关键部件进行频谱特性分析采用的是FFT、DFT、倒谱、包络谱这些频谱分析方法。
在步骤2)中,所述控制参数包括风力发电机组的塔架及传动链加阻、变桨与转矩滤波器和PID参数。
在步骤2)中,m的取值为10~20。
在步骤2)中,所述学习因子c1和c2通常等于2。
本发明与现有技术相比,具有如下优点与有益效果:
1、无需增加硬件传感器,且不影响机组的原有安全保护功能,就能对风力发电机组关键部件进行提前故障预警,防微杜渐,具体是在机组各部件出现异常的初期,进行故障预警,及时检修,防止部件出现严重的损坏,降低维修成本。
2、可实现在线调整优化控制参数(及时准确地在线调整合理的控制参数),防止突发事故,减少机组的停机时间,能有效提高风力发电机组的整体动态特性,节约机组的运维成本。
3、提升机组性能,进而提升发电量,增加风电场的直接经济效益。
附图说明
图1为本发明方法的框架图。
图2为经典粒子群优化算法流程图。
图3为本发明改进后的粒子群参数优化算法流程图。
图4为实际风机运行转速频谱图。
图5为实际风机运行Y轴振动谱图。
具体实施方式
下面结合具体实施例对本发明作进一步说明。
如图1所示,本实施例所提供的风力发电机组故障预警及控制参数在线优化方法,主要包含两部分:故障预警和控制参数在线优化,其具体情况如下:
1)故障预警
风力发电机组的状态监测技术是综合利用对风电设备振动、温度、压力状态信号进行检测分析结果,结合诊断对象的工作环境,来定量识别风力发电机组各个部件的运行状态。而本发明提出的故障预警判断依据是风力发电机组不同部件的频谱特性。
首先,实时采集风力发电机组运行数据,包括振动、转速、温度、桨叶角度、变桨速率等易测的时域信号;然后采用FFT、DFT、倒谱、包络谱等频谱分析方法(优选FFT频谱分析方法)对风力发电机组的关键部件进行频谱特性进行分析,最后对风力发电机组各个部件的特征频率及对应幅值运行状态进行监测。其中,当风力发电机组各部件正常运作时,说明其特征频率不偏移且对应的幅值在一定阈值范围内;若某部件的特征频率发生较大范围变化或对应的幅值超出理论阈值,则说明该部件损坏或出现异常。此时,机组则报出相应的故障,为发生损坏的零件进行预警。如图4、图5所示,我们举实例说明机组实时信号的频谱特性异常的故障预警,具体如下:
图4是某风场某台机组实际运行的转速信号频谱图。由图可知,该机组的传动链频率的能量幅值超出了正常能量幅值范围,说明传动链固有频率被激励,容易引起打滑或者传动链部件的损坏等问题,此时则机组会报出传动链频率能量过高的故障预警。现场运维人员也可第一时间进行硬件排查检修,避免由于发现较晚引起的更大故障。
图5是某风场某台机组实际运行的Y轴振动频谱图。正常运行机组的Y轴振动频谱图中,有塔架的固有频谱、3P频率等低频;但由图可知,该机组的Y轴振动频谱中出现了除了塔架的固有频谱、3P频率以外的低频信号0.18Hz。该低频信号的出现说明变桨环路参数出现了问题,此刻就需要对变桨增益参数进行调整,以避免引起变桨系统的损坏。
2)控制参数在线优化
当风力发电机组的某部件出现故障预警时,说明此刻的系统控制参数不利于当前机组的稳定运行,需要进行参数调整。本发明根据模拟鸟群捕食行为的粒子群算法,提出一种满足风力发电机组实际运行要求的改进的粒子群参数优化算法,对控制参数进行在线调整优化,防微杜渐,保证机组的可靠运行。其中,经典的粒子群优化算法具体算法流程参见图2所示。
由于机组出厂前,已经基于理论模型设定了最优的机组运行参数,故本发明所描述的控制参数在线优化,是在设定初始最优参数的基础上,在一定比例范围内,进行系数的微调。针对风力发电机组实际运行的特殊性,本发明对上述经典的粒子群优化算法进行相关步骤的改进优化,具体算法流程如图3所示。改进后的粒子群参数优化算法的基本步骤如下:
2.1)假设在D维的目标空间中,有m个粒子组成一个群落,初始化第i个粒子的位置为Xi=(xi1,xi2,…,xiD),第i个粒子的速度为Vi=(vi1,vi2,…,viD),i=1,2,…,m;
2.2)输入第i个粒子参数,实时存储机组运行的数据;
2.3)若存储的数据量达标,则对存储的数据进行频谱分析,获得特征频率的能量幅值;若存储的数据量不达标,则返回步骤2.2)继续执行;
2.4)更新下一组粒子参数,并判断是否达到m个粒子;若达到m个粒子,则进入步骤2.5);否则,返回步骤2.2)继续执行;
2.5)根据频率幅值,找出局部最优解Pbesti和全局最优解Gbest;
2.6)根据下列公式不断更新粒子的速度及位置:
式中,t表示时刻,d表示目标空间维度d=1,2,…,D,vid(t)表示t时刻第i个粒子在d维的更新速度,vid(t+1)表示t+1时刻第i个粒子在d维的更新速度,xid(t)表示t时刻第i个粒子在d维的位置,xid(t+1)表示t+1时刻第i个粒子在d维的位置,c1和c2表示学习因子,rand()表示随机数,Pbesti表示第i个粒子的局部最优解,Gbest表示全局最优解;
2.7)若达到最大迭代次数或者小于或等于频率幅值的正常阈值,则输出最优参数;否则,返回步骤2.2)继续执行,直到符合终止条件。
其中,上述算法的参数选择如下:
粒子数:一般取10~20,对于大部分的问题10个粒子已经足够可以取得好的结果。
粒子的维数:由优化问题决定,如变桨控制的目标参数有2个。
粒子的范围:由优化问题决定,每一维可是设定不同的范围。
学习因子:c1和c2通常等于2。
上述算法可在保证机组稳定运行的基础上,对有偏差的控制参数进行微调,防止机组发生较大的故障或部件损坏。
综上所述,本发明根据机组的振动、转速、温度、桨叶角度、变桨速率等易测的时域信号,采用频谱分析方法提取出机组各大部件的特征频率,对风力发电机组各个部件的运行状态进行监测。风力发电机组各部件正常运作时,其特征频率在一定阈值范围内;当某部件的特征频率超出阈值,则说明该部件损坏或出现异常。此时,机组报出相应的故障预警,并对塔架及传动链加阻、变桨与转矩滤波器和PID参数等主控参数进行在线优化调整,防微杜渐,保证机组的可靠运行,具有实际推广价值,值得推广。
以上所述实施例只为本发明之较佳实施例,并非以此限制本发明的实施范围,故凡依本发明之形状、原理所作的变化,均应涵盖在本发明的保护范围内。
Claims (5)
1.一种风力发电机组故障预警及控制参数在线优化方法,其特征在于,包括以下步骤:
1)故障预警
首先,实时采集风力发电机组的运行数据,包括振动、转速、温度、桨叶角度、变桨速率;然后对风力发电机组的关键部件进行频谱特性分析,最后再对风力发电机组各个部件的特征频率及对应幅值运行状态进行监测;其中,当风力发电机组各部件正常运作时,说明其特征频率不偏移且对应的幅值在预设阈值范围内;若某部件的特征频率发生大范围变化或对应的幅值超出理论阈值,则说明该部件损坏或出现异常,此时,机组报出相应的故障,为发生损坏的零件进行预警;
2)控制参数在线优化
当风力发电机组的某部件出现故障预警时,说明此刻的系统控制参数不利于当前机组的稳定运行,需要进行参数调整;此时,采用一种改进的粒子群参数优化算法,对控制参数进行在线调整优化,以满足风力发电机组的实际运行要求,保证机组的可靠运行;其中,改进的粒子群参数优化算法包括以下步骤:
2.1)假设在D维的目标空间中,有m个粒子组成一个群落,初始化第i个粒子的位置为Xi=(xi1,xi2,…,xiD),第i个粒子的速度为Vi=(vi1,vi2,…,viD),i=1,2,…,m;
2.2)输入第i个粒子参数,实时存储机组运行的数据;
2.3)若存储的数据量达标,则对存储的数据进行频谱分析,获得特征频率的能量幅值;若存储的数据量不达标,则返回步骤2.2)继续执行;
2.4)更新下一组粒子参数,并判断是否达到m个粒子;若达到m个粒子,则进入步骤2.5);否则,返回步骤2.2)继续执行;
2.5)根据频率幅值,找出局部最优解Pbesti和全局最优解Gbest;
2.6)根据下列公式不断更新粒子的速度及位置:
式中,t表示时刻,d表示目标空间维度d=1,2,…,D,vid(t)表示t时刻第i个粒子在d维的更新速度,vid(t+1)表示t+1时刻第i个粒子在d维的更新速度,xid(t)表示t时刻第i个粒子在d维的位置,xid(t+1)表示t+1时刻第i个粒子在d维的位置,c1和c2表示学习因子,rand()表示随机数,Pbesti表示第i个粒子的局部最优解,Gbest表示全局最优解;
2.7)若达到最大迭代次数或者小于或等于频率幅值的正常阈值,则输出最优参数;否则,返回步骤2.2)继续执行,直到符合终止条件。
2.根据权利要求1所述的一种风力发电机组故障预警及控制参数在线优化方法,其特征在于:在步骤1)中,对风力发电机组的关键部件进行频谱特性分析采用的是FFT、DFT、倒谱、包络谱这些频谱分析方法。
3.根据权利要求1所述的一种风力发电机组故障预警及控制参数在线优化方法,其特征在于:在步骤2)中,所述控制参数包括风力发电机组的塔架及传动链加阻、变桨与转矩滤波器和PID参数。
4.根据权利要求1所述的一种风力发电机组故障预警及控制参数在线优化方法,其特征在于:在步骤2)中,m的取值为10~20。
5.根据权利要求1所述的一种风力发电机组故障预警及控制参数在线优化方法,其特征在于:在步骤2)中,所述学习因子c1和c2等于2。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711101416.2A CN108021026B (zh) | 2017-11-10 | 2017-11-10 | 一种风力发电机组故障预警及控制参数在线优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711101416.2A CN108021026B (zh) | 2017-11-10 | 2017-11-10 | 一种风力发电机组故障预警及控制参数在线优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108021026A CN108021026A (zh) | 2018-05-11 |
CN108021026B true CN108021026B (zh) | 2020-11-06 |
Family
ID=62079770
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711101416.2A Active CN108021026B (zh) | 2017-11-10 | 2017-11-10 | 一种风力发电机组故障预警及控制参数在线优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108021026B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109100609A (zh) * | 2018-07-23 | 2018-12-28 | 华北电力大学(保定) | 一种基于智能优化的双馈风机定子匝间短路故障的诊断方法 |
JP6856591B2 (ja) * | 2018-09-11 | 2021-04-07 | ファナック株式会社 | 制御装置、cnc装置及び制御装置の制御方法 |
CN112434879A (zh) * | 2020-12-07 | 2021-03-02 | 河北工业大学 | 电机故障预测修正系统、方法、产品、介质、电子设备 |
CN113009457B (zh) * | 2021-02-26 | 2023-07-28 | 北京小马易行科技有限公司 | 激光雷达是否出现故障的确定方法、装置、处理器与车辆 |
CN116449688B (zh) * | 2023-06-19 | 2023-10-31 | 深圳市大族机器人有限公司 | 机器人电机pid参数确定方法、装置和计算机设备 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9864820B2 (en) * | 2012-10-03 | 2018-01-09 | Operation Technology, Inc. | Generator dynamic model parameter estimation and tuning using online data and subspace state space model |
CN103278326A (zh) * | 2013-06-14 | 2013-09-04 | 上海电机学院 | 风力发电机组齿轮箱的故障诊断方法 |
CN103401625B (zh) * | 2013-08-23 | 2015-04-08 | 西安电子科技大学 | 基于粒子群优化算法的协作频谱感知优化方法 |
US8933572B1 (en) * | 2013-09-04 | 2015-01-13 | King Fahd University Of Petroleum And Minerals | Adaptive superconductive magnetic energy storage (SMES) control method and system |
CN103887815B (zh) * | 2014-02-21 | 2016-03-02 | 华南理工大学 | 基于运行数据的风电场参数辨识及动态等值方法 |
CN104779638A (zh) * | 2015-02-06 | 2015-07-15 | 华北水利水电大学 | 一种风电场内机组优化调度方法及调度装置 |
CN106096170B (zh) * | 2016-06-21 | 2019-05-17 | 山西大学 | 基于数据驱动的风电机组多变量故障预测方法 |
CN106682814B (zh) * | 2016-11-28 | 2021-01-01 | 华北电力大学 | 一种基于故障知识库的风电机组故障智能诊断方法 |
-
2017
- 2017-11-10 CN CN201711101416.2A patent/CN108021026B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108021026A (zh) | 2018-05-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108021026B (zh) | 一种风力发电机组故障预警及控制参数在线优化方法 | |
EP3121442B2 (en) | Operating wind turbines | |
CN110761947A (zh) | 一种风电机组偏航校准方法及系统 | |
CA2704988C (en) | Wind-turbine-dynamic-characteristics monitoring apparatus and method therefor | |
EP1911968A1 (en) | Control system for a wind turbine and method of controlling said wind turbine | |
CN109296506B (zh) | 一种风电机组的振动检测方法、控制方法及装置 | |
CN106704103B (zh) | 一种基于叶片参数自学习的风电机组功率曲线获取方法 | |
US20230265832A1 (en) | Load control method and apparatus for wind turbine generator system | |
CN107045574B (zh) | 基于svr的风力发电机组低风速段有效风速估计方法 | |
CN112031996B (zh) | 一种优化风电机组偏航运动时机舱振动超限的方法及系统 | |
CN111287911A (zh) | 一种风电机组疲劳载荷的预警方法和系统 | |
CN113153633A (zh) | 一种风电机组风向仪静态偏差校准方法 | |
CN113847199B (zh) | 一种基于机载雷达在线偏航系统的偏航优化控制方法 | |
Habibi et al. | Power improvement of non-linear wind turbines during partial load operation using fuzzy inference control | |
CN108760037B (zh) | 一种基于频谱分析的风力发电机叶片结构损伤检测方法 | |
JP2018178900A (ja) | 風力発電システム | |
CN112483334B (zh) | 基于边缘计算的风电机组智能控制方法 | |
EP3828408A1 (en) | A method and an apparatus for computer-implemented monitoring of a wind turbine | |
EP4045791B1 (en) | Method and an apparatus for computer-implemented monitoring of a wind turbine | |
CN115288929A (zh) | 一种计及载荷计算的风电机组偏航控制方法及系统 | |
CN113048012B (zh) | 基于混合高斯模型的风电机组偏航角度识别方法和装置 | |
CN111794921B (zh) | 一种基于迁移成分分析的陆上风电机组叶片结冰诊断方法 | |
Liu et al. | Data-driven fault identification of ageing wind turbine | |
CN114186407A (zh) | 一种参数可自适应调整的风电场尾流速度场计算方法及系统 | |
Bao et al. | Identification of yaw error inherent misalignment for wind turbine based on SCADA data: A data mining approach |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |