CN113742932A - 一种用于风电机组齿轮箱故障预测诊断的方法 - Google Patents
一种用于风电机组齿轮箱故障预测诊断的方法 Download PDFInfo
- Publication number
- CN113742932A CN113742932A CN202111068824.9A CN202111068824A CN113742932A CN 113742932 A CN113742932 A CN 113742932A CN 202111068824 A CN202111068824 A CN 202111068824A CN 113742932 A CN113742932 A CN 113742932A
- Authority
- CN
- China
- Prior art keywords
- model
- sequence
- data
- prediction
- value
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- YHXISWVBGDMDLQ-UHFFFAOYSA-N moclobemide Chemical compound C1=CC(Cl)=CC=C1C(=O)NCCN1CCOCC1 YHXISWVBGDMDLQ-UHFFFAOYSA-N 0.000 claims abstract description 17
- 238000012544 monitoring process Methods 0.000 claims abstract description 11
- 238000004458 analytical method Methods 0.000 claims abstract description 10
- 238000005311 autocorrelation function Methods 0.000 claims description 16
- 238000003745 diagnosis Methods 0.000 claims description 10
- 238000007689 inspection Methods 0.000 claims description 10
- 239000000126 substance Substances 0.000 claims description 9
- 238000012360 testing method Methods 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 8
- 238000010998 test method Methods 0.000 claims description 5
- 241001123248 Arma Species 0.000 claims description 3
- 230000003247 decreasing effect Effects 0.000 claims description 3
- 238000002474 experimental method Methods 0.000 claims description 2
- 238000004519 manufacturing process Methods 0.000 claims description 2
- 238000005070 sampling Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 6
- 238000012423 maintenance Methods 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 238000013528 artificial neural network Methods 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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
- F03D17/00—Monitoring or testing of wind motors, e.g. diagnostics
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/02—Gearings; Transmission mechanisms
- G01M13/021—Gearings
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/06—Wind turbines or wind farms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/10—Noise analysis or noise optimisation
-
- 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
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P70/00—Climate change mitigation technologies in the production process for final industrial or consumer products
- Y02P70/50—Manufacturing or production processes characterised by the final manufactured product
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Sustainable Development (AREA)
- Sustainable Energy (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了一种用于风电机组齿轮箱故障预测诊断的方法,该方法通过对连续变化的状态监测数据建立基于ARIMA的时间序列模型揭示数据之间的某种显著关联性,依据时间序列模型对历史数据进行延拓,从而实现对状态趋势的预测来获知复杂机械系统的运行状况,对齿轮箱出口压力历史样本数据建立了ARIMA模型,能够对短期出口压力状态的主要变化趋势进行预测,具有比较满意的预测精度,对齿轮箱故障状态分析具有一定参考价值,其中,由于单个状态监测数据变化的缓急程度以及模型参数适宜程度都有可能使局部采样点的预测偏差较大,还需要修正相关参数和配合其他连续变化的指标数据一起进一步提高故障趋势预测的准确性和可靠性。
Description
技术领域
本发明涉及风机故障诊断领域,尤其涉及一种用于风电机组齿轮箱故障预测诊断的方法。
背景技术
风电机组的主轴轴承作为风力发电设备的关键机械部件,其维修成本高,维护时间长,在非维修时期产生的异常很有可能造成主轴轴承发生故障,传统的SCADA系统无法及时准确的定位故障隐患,影响风电机组的正常运行和风电并网的稳定性,因此有必要对风电机组主轴轴承的实时运行状态和故障诊断进行深入研究。
传统的主轴轴承故障诊断方法依据主轴轴承在运行状态下的振动信号来对已掌握的故障类型进行模式识别并借助神经网络或支持向量机对故障作出诊断。然而神经网络易陷入局部最小问题,会出现过适应现象。另一方面,粒子群算法的参数优化则决定了收敛精度和预测准度,现有的支持粒子群模型的参数优化算法收敛速度较差,极大地影响了故障诊断的精度和速度。
发明内容
为了克服现有技术存在的缺点与不足,本发明提供一种用于风电机组齿轮箱故障预测诊断的方法。
本发明所采用的技术方案如下:
一种用于风电机组齿轮箱故障预测诊断的方法,该方法利用ARIMA模型对风电机组齿轮箱故障趋势预测,该方法处理具有非线性和非平稳性特征的齿轮箱运行状态监测数据,用以时间序列的自相关分析为基础的模型预测状态监测时间序列数据的趋势变化;
ARIMA模型通过实验或生产现场采集获得的时间序列数据首先需要经过前提处理平稳性和白噪声检验才能用于模型分析,然后进行模型识别、参数估计以及模型检验,最后根据检验模型是否有效来决定是返回前述步骤重新识别模型还是可以用来预测未来趋势。
前提处理与模型识别,由于平稳非白噪声数据序列才适用于ARIMA模型分析,因此首先对时间序列数据进行平稳性检验,一般通过时序图和相关图来检验。如果观察时序图或相关图中的数据序列不是在坐标轴附近上下波动而是存在一定的增长或下降趋势,则说明是非平稳数据序列,需要对数据进行差分处理。接着进行非白噪声检验,可以用后文提到的LB统计量的检验法来检验。如果数据序列的各项数值之间是随机的,并没有任何相关关系,则说明序列没有多少信息可用;如果数据序列非白噪声序列,即过去的值对现在有影响,才能建立ARIMA模型进行预测。
时间序列数据经过上述前提处理后,需要进行模型识别,即求出序列样本的自相关函数和偏自相关函数的值,然后根据他们表现出来的性质确定模型。
自相关函数根据以下公式求得:
所述偏自相关函数的值可以利用样本自相关系数的值,根据以下公式求得:
确定模型类型后,对模型进行定阶处理,即选择最优模型,采用贝叶斯信息准则来从所有通过检验的模型中选择贝叶斯信息准则值最小的模型为最优模型,贝叶斯信息准则函数如下所示:
其中,n为样本数,这样就可以估计出所有的参数值。
其中,是样本k阶滞后的相关系数,n为样本个数,QLB近似地服从自由度为m的χ2分布;假设残差序列为白噪声序列,在给定的显著水平α下,当QLB的P值大于α时,则接受原假设,即原数据序列是白噪声序列,模型有效;否则,就拒绝原假设,模型无效需要重新建立模型。
模型预测,将数据代入进行实际预测时,期望是得到无误差或误差尽可能小的预测值,即根据式计算所得的均方差最小化原则得到最佳未来预测值。
该方法能够适应齿轮箱运行状态监测数据随时间的变化特征,反映出一定的运行状态变化趋势,具有较好的预测精度和较大的应用范围,对风电机组其他部件的故障趋势预测具有一定的应用参考价值。
附图说明
图1为本发明的ARIMA建模步骤流程图。
具体实施方式
需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互结合,下面结合附图和有具体实施例对本申请作进一步详细说明。
如图1所示,一种用于风电机组齿轮箱故障预测诊断的方法,该方法利用ARIMA模型对风电机组齿轮箱故障趋势预测,该方法处理具有非线性和非平稳性特征的齿轮箱运行状态监测数据,用以时间序列的自相关分析为基础的模型预测状态监测时间序列数据的趋势变化。
所述ARIMA模型通过实验或生产现场采集获得的时间序列数据首先需要经过平稳性和白噪声检验才能用于模型分析,然后进行模型识别、参数估计以及模型检验,最后根据检验模型是否有效来决定是返回前述步骤重新识别模型还是可以用来预测未来趋势。
前提处理与模型识别,平稳非白噪声数据序列用于ARIMA模型分析,首先对时间序列数据进行平稳性检验,通过时序图和相关图来检验,如果观察时序图或相关图中的数据序列不是在坐标轴附近上下波动而是存在一定的增长或下降趋势,则说明是非平稳数据序列,对数据进行差分处理;接着进行非白噪声检验,用LB统计量的检验法来检验;如果数据序列的各项数值之间是随机的,并没有任何相关关系,则说明序列没有多少信息可用;如果数据序列非白噪声序列,即过去的值对现在有影响,建立ARIMA模型进行预测;
时间序列数据经过所述前提处理后,进行模型识别,即求出序列样本的自相关函数和偏自相关函数的值,然后根据他们表现出来的性质确定模型。
样本自相关系数根据以下公式求得:
样本偏自相关系数可以利用样本自相关系数的值,根据以下公式求得:
确定模型类型后,还需要对模型进行定阶处理,即选择最优模型。这里采用贝叶斯信息准则(Bayesian Information Criterion,BIC)来从所有通过检验的模型中选择BIC值最小的模型为最优模型。BIC函数如下所示:
参数估计,在上一步模型识别的基础上,利用最小二乘法估计模型ARMA(p,q)的未知参数,即对θ1,θ2,…,θq以及白噪声方差σ2进行估计,分别得到估计值例如采用最小二乘法对上述式(1)计算残差平方和为
其中,n为样本数,这样就可以估计出所有的参数值。
模型检验,主要是验证模型的有效性,对于已经进行过模型识别和参数估计的模型,要通过残差序列检验,即判断数据序列是否为白噪声序列来决定检验是否通过。Ljung和Box提出了LB统计量的检验法,具体公式如下:
其中,是样本k阶滞后的相关系数,n为样本个数,QLB近似地服从自由度为m的χ2分布。假设残差序列为白噪声序列,那么在给定的显著水平α下,当QLB的P值大于α时,则接受原假设,即原数据序列是白噪声序列,模型有效;否则,就拒绝原假设,模型无效需要重新建立模型。
模型预测,将数据代入进行实际预测时,期望是得到无误差或误差尽可能小的预测值,即根据式(1)计算所得的均方差最小化原则得到最佳未来预测值。
该方法通过对连续变化的状态监测数据建立基于ARIMA的时间序列模型揭示数据之间的某种显著关联性,依据时间序列模型对历史数据进行延拓,从而实现对状态趋势的预测来获知复杂机械系统的运行状况,对齿轮箱出口压力历史样本数据建立了ARIMA模型,能够对短期出口压力状态的主要变化趋势进行预测,具有比较满意的预测精度,对齿轮箱故障状态分析具有一定参考价值,其中,由于单个状态监测数据变化的缓急程度以及模型参数适宜程度都有可能使局部采样点的预测偏差较大,还需要修正相关参数和配合其他连续变化的指标数据一起进一步提高故障趋势预测的准确性和可靠性。
尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解的是,在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种等效的变化、修改、替换和变型,本发明的范围由所附权利要求及其等同范围限定。
Claims (7)
1.一种用于风电机组齿轮箱故障预测诊断的方法,其特征在于,该方法利用ARIMA模型对风电机组齿轮箱故障趋势预测,该方法处理具有非线性和非平稳性特征的齿轮箱运行状态监测数据,用以时间序列的自相关分析为基础的模型预测状态监测时间序列数据的趋势变化;
所述ARIMA模型通过实验或生产现场采集获得的时间序列数据首先经过前提处理平稳性和白噪声检验才能用于模型分析,然后进行模型识别、参数估计以及模型检验,最后根据检验模型是否有效来决定是返回前述步骤重新识别模型还是用来预测未来趋势。
2.如权利要求1所述的一种用于风电机组齿轮箱故障预测诊断的方法,其特征在于,所述ARIMA模型,首先对时间序列数据进行平稳性检验,通过时序图和相关图来检验,如果观察时序图或相关图中的数据序列不是在坐标轴附近上下波动而是存在一定的增长或下降趋势,则说明是非平稳数据序列,对数据进行差分处理,接着进行非白噪声检验,利用LB统计量的检验法来检验,如果数据序列的各项数值之间是随机的,并没有任何相关关系,则说明序列没有多少信息可用;如果数据序列非白噪声序列,即过去的值对现在有影响,才建立ARIMA模型进行预测;
时间序列数据经过处理后,进行模型识别,即求出序列样本的自相关函数和偏自相关函数的值,然后根据他们表现出来的性质确定模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111068824.9A CN113742932A (zh) | 2021-09-13 | 2021-09-13 | 一种用于风电机组齿轮箱故障预测诊断的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111068824.9A CN113742932A (zh) | 2021-09-13 | 2021-09-13 | 一种用于风电机组齿轮箱故障预测诊断的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113742932A true CN113742932A (zh) | 2021-12-03 |
Family
ID=78738325
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111068824.9A Pending CN113742932A (zh) | 2021-09-13 | 2021-09-13 | 一种用于风电机组齿轮箱故障预测诊断的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113742932A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115131943A (zh) * | 2022-07-07 | 2022-09-30 | 杭州申昊科技股份有限公司 | 一种基于深度学习的声光联动预警方法 |
CN117312926A (zh) * | 2023-11-28 | 2023-12-29 | 沧州经济开发区武理工京津冀协同产业科技研究院 | 一种变压器故障预测的方法、装置、电子设备及介质 |
CN117666462A (zh) * | 2024-01-31 | 2024-03-08 | 成都苔岑智能设备有限公司 | 一种基于多参数独立调配的plc产品标准化控制系统 |
CN118088666A (zh) * | 2024-04-28 | 2024-05-28 | 南京讯联液压技术股份有限公司 | 一种风电齿轮箱分布式智能润滑监测与故障预警系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105158598A (zh) * | 2015-08-15 | 2015-12-16 | 国家电网公司 | 一种适用于电力设备的故障预测方法 |
-
2021
- 2021-09-13 CN CN202111068824.9A patent/CN113742932A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105158598A (zh) * | 2015-08-15 | 2015-12-16 | 国家电网公司 | 一种适用于电力设备的故障预测方法 |
Non-Patent Citations (1)
Title |
---|
杨艺等: "基于ARIMA风电机组齿轮箱故障趋势预测方法研究", 《重庆工商大学学报》, vol. 36, no. 3, pages 87 - 93 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115131943A (zh) * | 2022-07-07 | 2022-09-30 | 杭州申昊科技股份有限公司 | 一种基于深度学习的声光联动预警方法 |
CN115131943B (zh) * | 2022-07-07 | 2023-10-31 | 杭州申昊科技股份有限公司 | 一种基于深度学习的声光联动预警方法 |
CN117312926A (zh) * | 2023-11-28 | 2023-12-29 | 沧州经济开发区武理工京津冀协同产业科技研究院 | 一种变压器故障预测的方法、装置、电子设备及介质 |
CN117666462A (zh) * | 2024-01-31 | 2024-03-08 | 成都苔岑智能设备有限公司 | 一种基于多参数独立调配的plc产品标准化控制系统 |
CN117666462B (zh) * | 2024-01-31 | 2024-05-17 | 成都苔岑智能设备有限公司 | 一种基于多参数独立调配的plc产品标准化控制系统 |
CN118088666A (zh) * | 2024-04-28 | 2024-05-28 | 南京讯联液压技术股份有限公司 | 一种风电齿轮箱分布式智能润滑监测与故障预警系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113742932A (zh) | 一种用于风电机组齿轮箱故障预测诊断的方法 | |
Song et al. | Step-by-step fuzzy diagnosis method for equipment based on symptom extraction and trivalent logic fuzzy diagnosis theory | |
Huo et al. | Incipient fault diagnosis of roller bearing using optimized wavelet transform based multi-speed vibration signatures | |
Yu | Health condition monitoring of machines based on hidden Markov model and contribution analysis | |
US7275018B2 (en) | Method of condition monitoring | |
CN102033200B (zh) | 基于统计模型的交流电机在线监测和诊断方法 | |
CN112799898B (zh) | 基于分布式故障检测的互联系统故障节点定位方法及系统 | |
CN109710983B (zh) | 一种基于关键性能指标的柴油机气缸分层故障诊断方法 | |
Hayashi et al. | Study of machine fault diagnosis system using neural networks | |
Hwang et al. | SVM-RBM based predictive maintenance scheme for IoT-enabled smart factory | |
CN113947017A (zh) | 一种滚动轴承剩余使用寿命预测方法 | |
CN111504647A (zh) | 基于ar-mset的滚动轴承的性能退化评估方法 | |
CN117270514B (zh) | 基于工业物联网的生产过程全流程故障检测方法 | |
Barbieri et al. | Sensor-based degradation prediction and prognostics for remaining useful life estimation: Validation on experimental data of electric motors | |
CN115186701A (zh) | 轴承寿命预测方法、装置、电子设备及存储介质 | |
CN116976191A (zh) | 一种水轮发电机组轴系稳定性劣化趋势预测方法 | |
Xu et al. | Anomaly detection with gru based bi-autoencoder for industrial multimode process | |
CN113237619B (zh) | 变转速旋转机械振动的故障预警方法、装置、设备以及存储介质 | |
Yang et al. | Application research of arima model in wind turbine gearbox fault trend prediction | |
CN117786461A (zh) | 水泵的故障诊断方法、控制设备及其存储介质 | |
CN111474476B (zh) | 一种电机故障预测方法 | |
Cao et al. | Remaining useful life prediction of wind turbine generator bearing based on EMD with an indicator | |
CN116894475A (zh) | 一种基于神经网络的燃气轮机状态预警方法 | |
CN108629077B (zh) | 故障诊断装置、监视装置、故障诊断方法以及记录介质 | |
Farahat et al. | Similarity-based feature extraction from vibration data for prognostics |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20211203 |