CN115081331A - 一种基于状态参数重构误差的风电机组运行状态异常检测方法 - Google Patents

一种基于状态参数重构误差的风电机组运行状态异常检测方法 Download PDF

Info

Publication number
CN115081331A
CN115081331A CN202210754815.3A CN202210754815A CN115081331A CN 115081331 A CN115081331 A CN 115081331A CN 202210754815 A CN202210754815 A CN 202210754815A CN 115081331 A CN115081331 A CN 115081331A
Authority
CN
China
Prior art keywords
data
state
monitoring index
reconstruction error
wind turbine
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
Application number
CN202210754815.3A
Other languages
English (en)
Inventor
陈俊生
孙荣利
刘明杰
石钧仁
黄健
严冬
马艺玮
李本川
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing University of Post and Telecommunications
Original Assignee
Chongqing University of Post and Telecommunications
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chongqing University of Post and Telecommunications filed Critical Chongqing University of Post and Telecommunications
Priority to CN202210754815.3A priority Critical patent/CN115081331A/zh
Publication of CN115081331A publication Critical patent/CN115081331A/zh
Pending legal-status Critical Current

Links

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
    • G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/08—Learning methods
    • 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/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
    • 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)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Physics (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Computational Mathematics (AREA)
  • Artificial Intelligence (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Databases & Information Systems (AREA)
  • Computational Linguistics (AREA)
  • Probability & Statistics with Applications (AREA)
  • Operations Research (AREA)
  • Evolutionary Biology (AREA)
  • Health & Medical Sciences (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Algebra (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Medical Informatics (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明属于新能源发电领域，具体涉及一种基于状态参数重构误差的风电机组运行状态异常检测方法，该方法包括：采用正常运行的风电机组数据设置监测指标阈值和连续越限时间阈值；将待检测数据经过滑动窗口处理获得增广状态数据矩阵，并对待检测数据进行重构；计算重构误差和监测指标；将计算出的重构误差和监测指标分别与设置的监测指标阈值和连续越限时间阈值进行对比，得到数据检测结果；本发明对风电机组运行数据进行重构，通过重构误差，以及连续越限数，便可通过辨识模型判断运行数据异常状态，计算方便快捷，能够获得较高异常辨识能力。

Description

一种基于状态参数重构误差的风电机组运行状态异常检测 方法
技术领域
本发明属于新能源发电领域，具体涉及一种基于状态参数重构误差的风电机组运行状态异常检测方法。
背景技术
风电场数据采集与监控系统提供了大量的风电机组运行数据，但其中绝大部分为机组健康运行数据，故障异常数据缺乏甚至难以获得，且当SCADA系统中风电机组状态参数发生越限报警时，机组的故障或异常往往已发展至一定程度。因此，如何有效地综合利用现有大量健康运行数据实现风电机组运行状态监测是所有风电研究者所要解决的难题。如果能够准确的辨识风电机组运行状态异常数据，对提高风电机组运行安全性、经济性和健康管理水平具有良好的实际应用价值。
目前基于SCADA数据的风电机组运行状态监测方法是通过对单一或少量机组状态参数的分析实现机组部件的运行状态监测和判断，难以全面反映机组整机的运行状态。
发明内容
针对现有的基于SCADA数据的风电机组运行状态监测方法中通过单一或少量数据对机组运行状态判断结果低的问题，本发明提出了一种基于状态参数重构误差的风电机组运行状态异常检测方法，该方法包括：
S1：获取待检测数据；
S2：采用滑动窗口对待检测的数据进行处理，得到增广状态数据矩阵；
S3：将待检测数据输入训练好的SDAE网络中，得到重构数据；
S4：根据增广状态数据矩阵和重构数据计算重构误差，根据重构误差计算监测指标；
S5：设置监测指标阈值和连续越限时间阈值；
S6：将计算出的监测指标与设置的监测指标阈值进行比较，若小大于设置的监测指标阈值，则执行步骤S7，否则返回步骤S1；
S7：计算重构数据的连续越限数，将连续越限数与设置的连续越限时间阈值进行对比，若大于连续越限时间阈值，则待检测数据为异常数据，否则待检测数据为噪声数据；所述连续越限数为：在监测的指标中出现连续多次超过设置的监测指标阈值的数量；
S8：根据异常数据得到风电机组运行状态的异常组件。
优选的，采用滑动窗口对待检测的数据进行处理的过程包括：
S21：设置滑动窗口的宽度和滑动窗口的数量；
S22：采用滑动窗口对待检测的数据进行采集；
S23：根据所有的滑动窗口采集的数据构建增广状态数据矩阵。
优选的，对SDAE网络进行训练的过程包括：SDAE网络包括第一DAE模型和第二DAE模型；
步骤1：获取原始SCADA数据，采用滑动窗口对原始SCADA数据进行处理，得到原始数据的增广矩阵；
步骤2：初始化SDAE网络参数，该参数包括初始加噪比c0、加噪比阈值cr
步骤3：根据加噪比C0对增广矩阵进行加噪处理；
步骤4：将经过加噪处理后的增广矩阵数据输入到第一DAE模型中口模型进行训练，得到第一权值，并更新加噪比；
步骤5：将更新后的加噪比与加噪比阈值进行对比，若不相等，则返回步骤3；若相等，则执行步骤6；
步骤6：将第一DAE模型的隐藏层输出和初始加噪比c0输入到第二DAE模型中进行训练，得到第二权值，并对加噪比进行更新；
步骤7：将更新后的加噪比与加噪比阈值进行对比，若小于加噪比阈值，则训练结束，输出训练参数，否则返回步骤6。
优选的，根据重构误差计算监测指标的公式为：
Figure BDA0003722055210000031
其中，tk表示监测指标，
Figure BDA0003722055210000032
表示第k个正常样本的重构误差，N表示重构后的数据数量，u表示均值向量，Cov-1表示协方差矩阵的逆矩阵。
优选的，设置监测指标阈值和连续越限时间阈值的过程包括：
S51：获取系统正常运行的原始SCADA数据，并对数据进行预处理；对预处理后的数据进行相关性分析，得到机组各状态参数时间序列数据的短时自相关性；
S52：采用滑动窗口对机组状态数据进行处理，得到系统正常运行的增广状态数据矩阵；
S53：根据机组各状态参数时间序列数据的短时自相关性将SCADA数据数据输入到训练好的SDAE网络中，得到第二重构数据；
S54：根据重构数据和正常运行的增广状态数据矩阵计算机组正常数据的重构误差和监测指标；
S55：根据机组正常数据的监测指标构确定监测指标阈值；
S56：根据滑动窗口大小和监测指标连续越限的最大值构建连续越限数指标阈值。
进一步的，对预处理后的数据进行相关性分析包括：采用综合相关性指标分析风电机组状态参数间的互相关性；采用自相关系数和偏自相关系数分析机组各状态参数时间序列数据的短时自相关性；根据SCADA系统采集的数据；计算时间序列延迟k阶的自相关系数；根据自相关系数计算数据的偏自相关系数。
进一步的，偏自相关系数的表达为：
Figure BDA0003722055210000033
其中，
Figure BDA0003722055210000034
表示第1个偏自相关系数，ρj表示第j个自相关系数，
Figure BDA0003722055210000035
表示第k+1个偏自相关系数，k表示序列的长度。
优选的，确定监测指标阈值的的过程包括：构建重构误差马氏距离的概率密度函数
Figure BDA0003722055210000041
设置置信度α；根据概率密度函数
Figure BDA0003722055210000042
计算机组状态监测指标小于监测指标阈值的概率值P(t<td)；根据设置的置信度α和概率值P(t<td)反向确定监测指标阈值td。
进一步的，计算机组状态监测指标小于监测指标阈值的概率值的公式为：
Figure BDA0003722055210000043
Figure BDA0003722055210000044
其中，P表示机组状态小于监测指标阈值的概率值，
Figure BDA0003722055210000045
表示重构误差马氏距离的概率密度函数，n表示样本个数数据，tk为样本数据中第k个样本数据点的监测指标，h为带宽参数，K(·)为核函数。
优选的，连续越限数指标阈值的表达式为：
Qd＝max(Qmax,β)
其中，Qmax为训练数据的监测指标连续越限最大值，β为滑动窗口宽度。
本发明的有益效果：
本发明采用滑动窗口的方法对待检测的数据进行处理，得到增广状态数据矩阵，并将SDAE网络重构后的数据与增广状态数据矩阵求取重构误差，将重构误差与设置的阈值进行比较，得到风电机组运行状态异常检测；将输入数据重构数据进行联系，使得最终的识别结果更准确；本发明设计了一种监测指标阈值和连续越限时间阈值的设置方法，通过该方式，使得设置阈值更精确，得到结果更准确；本发明对待检测的风电机组运行数据进行重构，通过重构误差和连续越限数判断运行数据异常状态，计算方便快捷，识别结果更准确。
附图说明
图1为本发明的基于状态参数重构误差的风电机组运行状态异常检测方法流程图；
图2为本发明的对对SDAE网络进行训练的流程图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
一种基于状态参数重构误差的风电机组运行状态异常检测方法的具体实施方式，如图1所示，该方法包括：
S1：获取待检测数据；
S2：采用滑动窗口对待检测的数据进行处理，得到增广状态数据矩阵；
S3：将待检测数据输入训练好的SDAE网络中，得到重构数据；
S4：根据增广状态数据矩阵和重构数据计算重构误差，根据重构误差计算监测指标；
S5：设置监测指标阈值和连续越限时间阈值；
S6：将计算出的监测指标与设置的监测指标阈值进行比较，若小大于设置的监测指标阈值，则执行步骤S7，否则返回步骤S1；
S7：计算重构数据的连续越限数，将连续越限数与设置的连续越限时间阈值进行对比，若大于连续越限时间阈值，则待检测数据为异常数据，否则待检测数据为噪声数据；所述连续越限数为：在监测的指标中出现连续多次超过设置的监测指标阈值的数量；连续多次为至少Qd次。
S8：根据异常数据得到风电机组运行状态的异常组件。
连续越限数：在采用SCADA系统对采集的原始运行数据进行不良数据识别和修正后得到的正常数据，在得到的数据中仍可能包含部分噪声数据，此类噪声数据会引起监测指标越限，但其连续越限时间较短，因此为了准确辨识机组运行状态，还需考虑监测指标的连续越限数指标Q。Q的阈值Qd是根据正常样本数据的监测指标统计分布获得。在经过滑动窗口对数据进行处理后，得到增广状态数据矩阵，增广状态数据矩阵的连续β个数据点均包含原始数据的部分数据信息，因此当原始数据中某个数据为噪声数据时，将引起增广状态数据矩阵连续个数据点的监测指标越限，因此监测指标连续越限数最小值应为滑动窗口的宽度。
采用滑动窗口对待检测的数据进行处理的过程包括：
步骤(21)、采用机组正常状态的数据矩阵XN,设置滑动窗口的宽度为β，m个样本数量共有m-β+1个滑动窗口，记
Figure BDA0003722055210000061
为第l个滑动窗口采集的第i个机组监测数据，则有公式：
Figure BDA0003722055210000062
其中，
Figure BDA0003722055210000063
表示第(l+β-1)个滑动窗口采集的第i个监测数据，T表示转置。
步骤(22)、得到的第l个滑动窗口数据Sl为
Figure BDA0003722055210000064
其中，
Figure BDA0003722055210000065
表示第l个滑动窗口采集的第n个机组监测数据；
步骤(23)、根据步骤(21)得到的第l个滑动窗口采集的第i个机组监测参量的数据和步骤(22)得到的第l个滑动窗口采集的数据，对状态参数数据采用滑动窗口处理，得到相应的增广状态数据矩阵为Y＝[S(1),S(2),…,S(m-β+1)]；其中，S(m-β+1)表示第m-β+1个滑动窗口采集的监测数据。
如图2所示，对SDAE网络进行训练的过程包括：SDAE网络包括第一DAE模型和第二DAE模型；
步骤1：获取原始SCADA数据，采用滑动窗口对原始SCADA数据进行处理，得到原始数据的增广矩阵；
步骤2：初始化SDAE网络参数，该参数包括初始加噪比c0、加噪比阈值cT；
步骤3：根据加噪比c0对增广矩阵进行加噪处理；
具体的，在得到增广状态数据矩阵YN后，根据初始加噪比c0，对矩阵YN加噪处理
Figure BDA0003722055210000071
得到输入数据
Figure BDA0003722055210000072
其中，c0表示初始加噪比，cT表示最终加噪比，c0>c1>…>cT，并且c0-c1＝…＝cT-1-cT＝Δc,Δc>0；q()表示加噪处理的过程，Δc表示两个加噪比的差值，是对加噪比的微调。
步骤4：将经过加噪处理后的增广矩阵数据输入到第一DAE模型中进行训练，得到第一权值，并更新加噪比；加噪比更新公式为：ct＝ct-1-Δc。
步骤5：将更新后的加噪比与加噪比阈值进行对比，若不相等，则返回步骤3；若相等，则执行步骤6；
步骤6：将第一DAE模型的隐藏层输出和初始加噪比c0输入到第二DAE模型中进行训练，得到第二权值，并对加噪比进行更新；
步骤7：将更新后的加噪比与加噪比阈值进行对比，若小于加噪比阈值，则训练结束，输出训练参数，否则返回步骤6。
训练结束后模型输出的最优参数包括W1、b1、W2和b2及h，其中，W1表示第一个权值，b1表示第一个偏置量，W2表示第二个权值，b2表示第二个偏置量，h表示隐藏层。
设置监测指标阈值和连续越限时间阈值的过程包括：
步骤1：获取系统正常运行的原始SCADA数据，并对数据进行预处理；对预处理后的数据进行相关性分析，得到机组各状态参数时间序列数据的短时自相关性；
步骤2：采用滑动窗口对机组状态数据进行处理，得到系统正常运行的增广状态数据矩阵；其中风电机组状态参数的自相关系数反映了状态参数先前时刻值与当前时刻值的线性相关程度；
步骤3：根据机组各状态参数时间序列数据的短时自相关性将SCADA数据数据输入到训练好的SDAE网络中，得到第二重构数据；
步骤4：根据重构数据和正常运行的增广状态数据矩阵计算机组正常数据的重构误差和监测指标；
步骤5：根据机组正常数据的监测指标构确定监测指标阈值；
步骤6：根据滑动窗口大小和监测指标连续越限的最大值构建连续越限数指标阈值。
对原始SCADA数据进行预处理和对预处理后的数据进行相关性分析的过程包括：
步骤(11)、对SCADA系统所采集到的原始机组状态数据进行不良数据识别与修正，能够得到机组正常状态的数据矩阵，为
Figure BDA0003722055210000081
Figure BDA0003722055210000082
其中，
Figure BDA0003722055210000083
表示第j个状态参数的第i个样本数据，n为采集的样本数据数量，m为机组状态参数的个数。
步骤(12)、采用综合相关性指标分析风电机组状态参数间的互相关性，其中的综合相关系数ICC的计算公式为：
Figure BDA0003722055210000084
其中，
Figure BDA0003722055210000085
为Pearson相关系数，
Figure BDA0003722055210000086
为Kendall相关系数，
Figure BDA0003722055210000087
为Spearman相关系数，x(k)表示第k个状态参数的样本数据，x(l)表示第l个状态参数的样本数据。
步骤(13)、采用自相关系数和偏自相关系数分析机组各状态参数时间序列数据的短时自相关性；根据SCADA系统采集的数据x1,x2,…,xn，定义时间序列{xt,t∈T}，时间序列{xt}延迟k阶的自相关系数为：
Figure BDA0003722055210000088
偏自相关系数：
Figure BDA0003722055210000089
其中，ρk表示时间序列延迟k阶的自相关系数，xt表示时间序列，cov表示协方差矩阵，xt-k表示延迟k阶的时间序列，var()表示协方差概率密度的二阶矩，
Figure BDA00037220552100000810
表示延迟k阶偏自相关系数为按j阶自回归模型对xt作线性最小方差估计的系数，
Figure BDA0003722055210000091
表示延迟k+1阶偏自相关系数为按j阶自回归模型对xt作线性最小方差估计的系数。
具体的，计算机组正常数据的重构误差EN和机组正常数据的监测指标t的过程包括：
步骤(41)、首先采集到原始数据XN后，滑动窗口处理得到增广状态数据矩阵YN，利用多重加噪比的SDAE模型训练，得到重构数据ZN，则得到重构误差为EN＝YN-ZN；
步骤(42)、因为SMDAE模型的输入为多元状态数据，所以EN为多元数据，为了实现机组异常状态，需从EN中提取出状态监测指标；采用马氏距离，来度量未知数据与已知道数据间的相似性；
步骤(43)、设置EN的均值向量为μ＝(μ1,μ2,…,μn)，选取重构误差的马氏距离为监测所需指标，得到监测指标公式为
Figure BDA0003722055210000092
其中，
Figure BDA0003722055210000093
表示第k个正常样本的重构误差，u表示均值向量，Cov-1表示协方差矩阵的逆矩阵。
确定监测指标的阈值td和监测指标连续越限时间阈值Qd的过程包括：
步骤(51)、确定机组正常运行数据的重构误差马氏距离概率密度函数
Figure BDA0003722055210000094
Figure BDA0003722055210000095
其中，tk为样本数据中第k个样本数据点的监测指标，n表示样本数据个数，h为带宽参数，K(·)为核函数，满足公式：
Figure BDA0003722055210000096
其中，K(u)表示u的核函数，u2(K)表示约束系数与核函数的乘积。
步骤(52)、确定监测指标的阈值td，其过程包括：设置置信度α；根据概率密度函数
Figure BDA0003722055210000101
计算机组状态监测指标小于监测指标阈值的概率值P(t<td)；根据设置的置信度α和概率值P(t<td)反向确定监测指标阈值td；计算机组状态监测指标小于监测指标阈值的概率值的公式为：
Figure BDA0003722055210000102
Figure BDA0003722055210000103
其中，P表示机组状态小于监测指标阈值的概率值，
Figure BDA0003722055210000104
表示重构误差马氏距离的概率密度函数，n表示样本个数数据，tk为样本数据中第k个样本数据点的监测指标，h为带宽参数，K(·)为核函数。
步骤(53)、噪声数据会造成监测指标越限，这类噪声的连续越限时间较短，为了能够准确辨识，因此需要进一步采用连续越限数指标Q；连续越限数指标Q的阈值Qd公式满足：
Qd＝max(Qmax,β)
其中，Qmax为训练数据的监测指标连续越限最大值；β为滑动窗口宽度。
一种基于状态参数重构误差的风电机组运行状态异常检测方法的具体实施方式，该方法包括：
第一步，采用原始正常运行的SCADA数据，记为XN,状态参数间的互相关性采用综合相关性指标分析，采用自相关系数和偏自相关系数分析机组各状态参数时间序列数据的短时自相关性；采用滑动窗口技术处理机组状态数据，得到计算参数短时自相关性的增广状态数据矩阵YN；采用深度网络结构的SDAE模型，SDAE模型具有高强度非线性投射能力；在此，引入加噪比，他是DAE网络十分重要的调优参数，因此，采用多重加噪比训练SDAE模型，并获得正常数据的重构输出ZN；
第二步，根据公式EN＝YN-ZN计算机组正常数据的重构误差EN,根据重构误差计算监测指标的公式计算机组正常数据的监测指标t；根据机组状态监测指标小于监测指标阈值的概率值确定监测指标的阈值td，根据连续越限数指标阈值的表达式确定监测指标连续越限时间阈值Qd；
第三步，输入待检测数据XM,经滑动窗口处理获得相应的增广状态数据矩阵YM,根据正常样本训练得到的SDAE网络计算待检测数据的重构输出ZM,根据EM＝YM-XM计算待检测数据的重构误差EM,根据
Figure BDA0003722055210000111
计算待检测数据的监测指标tM；判断tM>td，如果错误，则返回待检测数据，如果正确，则统计获得连续越限数QM,判断QM>Qd，如果错误，则认定为噪声数据，如果正确，则表明是异常数据，接着计算个参数的贡献度以确定与机组状态异常最相关的状态参数或者部件。
以上所举实施例，对本发明的目的、技术方案和优点进行了进一步的详细说明，所应理解的是，以上所举实施例仅为本发明的优选实施方式而已，并不用以限制本发明，凡在本发明的精神和原则之内对本发明所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种基于状态参数重构误差的风电机组运行状态异常检测方法，其特征在于，该方法包括：
S1：获取待检测数据；
S2：采用滑动窗口对待检测的数据进行处理，得到增广状态数据矩阵；
S3：将待检测数据输入训练好的SDAE网络中，得到重构数据；
S4：根据增广状态数据矩阵和重构数据计算重构误差，根据重构误差计算监测指标；
S5：设置监测指标阈值和连续越限时间阈值；
S6：将计算出的监测指标与设置的监测指标阈值进行比较，若小大于设置的监测指标阈值，则执行步骤S7，否则返回步骤S1；
S7：计算重构数据的连续越限数，将连续越限数与设置的连续越限时间阈值进行对比，若大于连续越限时间阈值，则待检测数据为异常数据，否则待检测数据为噪声数据；所述连续越限数为：在监测的指标中出现连续超过设置的监测指标阈值的数量；
S8：根据异常数据得到风电机组运行状态的异常组件。
2.根据权利要求1所述的一种基于状态参数重构误差的风电机组运行状态异常检测方法，其特征在于，采用滑动窗口对待检测的数据进行处理的过程包括：
S21：设置滑动窗口的宽度和滑动窗口的数量；
S22：采用滑动窗口对待检测的数据进行采集；
S23：根据所有的滑动窗口采集的数据构建增广状态数据矩阵。
3.根据权利要求1所述的一种基于状态参数重构误差的风电机组运行状态异常检测方法，其特征在于，对SDAE网络进行训练的过程包括：SDAE网络包括第一DAE模型和第二DAE模型；
步骤1：获取原始SCADA数据，采用滑动窗口对原始SCADA数据进行处理，得到原始数据的增广矩阵；
步骤2：初始化SDAE网络参数，该参数包括初始加噪比c0、加噪比阈值cT；
步骤3：根据加噪比c0对增广矩阵进行加噪处理；
步骤4：将经过加噪处理后的增广矩阵数据输入到第一DAE模型中进行训练，得到第一权值，并更新加噪比；
步骤5：将更新后的加噪比与加噪比阈值进行对比，若不相等，则返回步骤3；若相等，则执行步骤6；
步骤6：将第一DAE模型的隐藏层输出和初始加噪比c0输入到第二DAE模型中进行训练，得到第二权值，并对加噪比进行更新；
步骤7：将更新后的加噪比与加噪比阈值进行对比，若小于加噪比阈值，则训练结束，输出训练参数，否则返回步骤6。
4.根据权利要求1所述的一种基于状态参数重构误差的风电机组运行状态异常检测方法，其特征在于，根据重构误差计算监测指标的公式为：
Figure FDA0003722055200000021
其中，tk表示监测指标，
Figure FDA0003722055200000022
表示第k个正常样本的重构误差，N表示重构后的数据数量，u表示均值向量，Cov-1表示协方差矩阵的逆矩阵。
5.根据权利要求1所述的一种基于状态参数重构误差的风电机组运行状态异常检测方法，其特征在于，设置监测指标阈值和连续越限时间阈值的过程包括：
S51：获取系统正常运行的原始SCADA数据，并对数据进行预处理；对预处理后的数据进行相关性分析，得到机组各状态参数时间序列数据的短时自相关性；
S52：采用滑动窗口对机组状态数据进行处理，得到系统正常运行的增广状态数据矩阵；
S53：根据机组各状态参数时间序列数据的短时自相关性将SCADA数据数据输入到训练好的SDAE网络中，得到第二重构数据；
S54：根据重构数据和正常运行的增广状态数据矩阵计算机组正常数据的重构误差和监测指标；
S55：根据机组正常数据的监测指标构确定监测指标阈值；
S56：根据滑动窗口大小和监测指标连续越限的最大值构建连续越限数指标阈值。
6.根据权利要求5所述的一种基于状态参数重构误差的风电机组运行状态异常检测方法，其特征在于，对预处理后的数据进行相关性分析包括：采用综合相关性指标分析风电机组状态参数间的互相关性；采用自相关系数和偏自相关系数分析机组各状态参数时间序列数据的短时自相关性；根据SCADA系统采集的数据；计算时间序列延迟k阶的自相关系数；根据自相关系数计算数据的偏自相关系数。
7.根据权利要求6所述的一种基于状态参数重构误差的风电机组运行状态异常检测方法，其特征在于，偏自相关系数的表达为：
Figure FDA0003722055200000031
其中，
Figure FDA0003722055200000032
表示第1个偏自相关系数，ρj表示第j个自相关系数，
Figure FDA0003722055200000033
表示第k+1个偏自相关系数，k表示序列的长度。
8.根据权利要求5所述的一种基于状态参数重构误差的风电机组运行状态异常检测方法，其特征在于，确定监测指标阈值的的过程包括：构建重构误差马氏距离的概率密度函数
Figure FDA0003722055200000034
设置置信度α；根据概率密度函数
Figure FDA0003722055200000035
计算机组状态监测指标小于监测指标阈值的概率值P(t<td)；根据设置的置信度α和概率值P(t<td)确定监测指标阈值td。
9.根据权利要求8所述的一种基于状态参数重构误差的风电机组运行状态异常检测方法，计算机组状态监测指标小于监测指标阈值的概率值的公式为：
Figure FDA0003722055200000041
Figure FDA0003722055200000042
其中，P表示机组状态小于监测指标阈值的概率值，
Figure FDA0003722055200000043
表示重构误差马氏距离的概率密度函数，n表示样本个数数据，tk为样本数据中第k个样本数据点的监测指标，h为带宽参数，K(·)为核函数。
10.根据权利要求5所述的一种基于状态参数重构误差的风电机组运行状态异常检测方法，其特征在于，连续越限数指标阈值的表达式为：
Qd＝max(Qmax,β)
其中，Qmax为训练数据的监测指标连续越限最大值，β为滑动窗口宽度。
CN202210754815.3A 2022-06-30 2022-06-30 一种基于状态参数重构误差的风电机组运行状态异常检测方法 Pending CN115081331A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210754815.3A CN115081331A (zh) 2022-06-30 2022-06-30 一种基于状态参数重构误差的风电机组运行状态异常检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210754815.3A CN115081331A (zh) 2022-06-30 2022-06-30 一种基于状态参数重构误差的风电机组运行状态异常检测方法

Publications (1)

Publication Number Publication Date
CN115081331A true CN115081331A (zh) 2022-09-20

Family

ID=83256690

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210754815.3A Pending CN115081331A (zh) 2022-06-30 2022-06-30 一种基于状态参数重构误差的风电机组运行状态异常检测方法

Country Status (1)

Country Link
CN (1) CN115081331A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115496644A (zh) * 2022-11-18 2022-12-20 山东超华环保智能装备有限公司 基于数据识别的固体废弃物处理设备监测方法
CN116664000A (zh) * 2023-06-13 2023-08-29 无锡物联网创新中心有限公司 基于长短期记忆网络的工业设备异常检测方法及相关装置

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115496644A (zh) * 2022-11-18 2022-12-20 山东超华环保智能装备有限公司 基于数据识别的固体废弃物处理设备监测方法
CN115496644B (zh) * 2022-11-18 2023-09-26 南通万达能源动力科技有限公司 基于数据识别的固体废弃物处理设备监测方法
CN116664000A (zh) * 2023-06-13 2023-08-29 无锡物联网创新中心有限公司 基于长短期记忆网络的工业设备异常检测方法及相关装置

Similar Documents

Publication Publication Date Title
CN110610035B (zh) 一种基于gru神经网络的滚动轴承剩余寿命预测方法
CN111562996B (zh) 一种关键性能指标数据的时序异常检测方法及系统
CN110071913B (zh) 一种基于无监督学习的时间序列异常检测方法
CN115081331A (zh) 一种基于状态参数重构误差的风电机组运行状态异常检测方法
CN111275288A (zh) 基于XGBoost的多维数据异常检测方法与装置
CN112414694B (zh) 基于多元状态估计技术的设备多级异常状态识别方法及装置
CN108435819B (zh) 一种铝型材挤压机能耗异常检测方法
CN111426905B (zh) 一种配电网同母线变关系异常诊断方法、装置及系统
CN115412455A (zh) 一种基于时间序列的服务器多性能指标异常检测方法及装置
CN113076986B (zh) 过滤式与封装式评价策略相结合的光伏故障电弧特征选择方法
CN114167838A (zh) 一种伺服系统多尺度健康评估与故障预测方法
CN107220500A (zh) 基于逆高斯过程的性能退化试验贝叶斯可靠性评估方法
CN110220885B (zh) 一种机械设备磨损状态综合判断方法
CN117074961A (zh) 一种变分自编码器动力电池异常检测方法、系统、装置及存储介质
CN114429238A (zh) 一种基于时空特征提取的风电机组故障预警方法
CN116840765B (zh) 一种基于多元时序分析的电压互感器误差状态评估方法
CN110852906B (zh) 一种基于高维随机矩阵进行窃电嫌疑识别的方法和系统
CN116128690B (zh) 一种碳排放量成本值计算方法、装置、设备及介质
CN116641941A (zh) 一种基于典型变量分析的液压系统早期故障动态检测方法
CN108594156B (zh) 一种改进的电流互感器饱和特性识别方法
CN111273126A (zh) 一种配电网拓扑快速感知方法
CN114397569A (zh) 基于vmd参数优化、样本熵的断路器故障电弧检测方法
CN113598784A (zh) 心律失常检测方法及系统
CN113051809A (zh) 一种基于改进受限玻尔兹曼机的虚拟健康因子构建方法
CN117349781B (zh) 一种变压器故障智能诊断方法及系统

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