CN110008565A - 一种基于运行参数关联性分析的工业过程异常工况预测方法 - Google Patents
一种基于运行参数关联性分析的工业过程异常工况预测方法 Download PDFInfo
- Publication number
- CN110008565A CN110008565A CN201910244872.5A CN201910244872A CN110008565A CN 110008565 A CN110008565 A CN 110008565A CN 201910244872 A CN201910244872 A CN 201910244872A CN 110008565 A CN110008565 A CN 110008565A
- Authority
- CN
- China
- Prior art keywords
- parameter
- relevance
- 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.)
- Granted
Links
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/20—Administration of product repair or maintenance
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- General Business, Economics & Management (AREA)
- Operations Research (AREA)
- Entrepreneurship & Innovation (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Testing And Monitoring For Control Systems (AREA)
Abstract
本发明公开了一种基于运行参数关联性分析的工业过程异常工况预测方法,可应用于工业过程的故障预测与健康管理。本发明从工业过程运行参数之间的关联性入手,基于运行参数的关联性分析进行异常工况预测。在单参数预测阶段,本发明根据已有传感器数据通过指数平滑方法对各个运行参数进行预测。在关联性分析阶段,本发明通过已知的各运行参数值以及参数预测值计算运行参数关联性,其中参数关联性由代表参数曲线的一系列指标的相似性表示。在关联性趋势预测阶段,本发明构造多元自回归模型对参数关联性进行预测。本发明提出的方法将运行参数关联性纳入考虑,能够获得更完备的设备异常信息以及更提前的预测结果,对于工业设备的故障预测有实际意义。
Description
技术领域
本发明属于可靠性工程技术领域,涉及一种基于运行参数关联性分析的工业过程异常工况预测方法。
背景技术
随着复杂***的不断出现以及工业过程实时监测的需求不断增加,现代工业设备在运行过程中往往配备多个传感器对其运行状态进行监测。同时,设备运行过程中可能会出现多种故障模式,某一故障可能对应若干征兆,在此情况下,单传感器信息已无法完全体现设备运行状态,基于多传感器信息的故障预测应运而生。基于多传感器信息的故障预测旨在利用全面的传感器信息分析设备的运行状态,从而进行更可靠的设备诊断和预测。随着传感技术的持续发展,使用多个传感器进行设备的状态监测、故障诊断和预测已经成为发展趋势。
对于设备运行过程中的多个传感器,其代表的运行参数并不是独立存在的,设备的运行过程中每一运行参数的变化实际上都是对设备当前运行状态的反应。在工业设备正常运行状态下,运行参数通常比较平稳,维持在相对稳定的水平,因而运行参数的关联性也比较稳定。但当异常工况出现时,各运行参数对于异常的响应是不同的,从而导致运行参数的关联性发生变化,其关联性以及关联性变化趋势必然隐含着设备异常乃至故障的信息。
发明内容
针对现有的技术状况,本发明的目的是针对设备在运行过程中各运行参数存在关联性的问题,通过运行参数的关联性对设备运行过程的状态进行判断,通过对运行参数关联性变化趋势的预测来进行设备的异常工况预测。
现将本发明的构思阐述如下:
本发明提出一种基于运行参数关联性变化趋势的工业过程异常工况预测方法,从工业过程运行参数之间的关联性入手,通过对运行参数关联性的趋势分析进行异常工况预测。在单参数预测阶段,本发明根据已有传感器数据通过指数平滑方法对各个运行参数进行预测。在关联性分析阶段,本发明通过已知的各运行参数值以及参数预测值计算运行参数关联性,其中参数关联性由代表参数曲线的一系列指标的相似性表示。在关联性趋势预测阶段,本发明构造多元自回归模型对参数关联性进行预测。本发明提出的方法将运行参数关联性纳入考虑,能够获得更完备的设备异常信息以及更提前的预测结果。
根据以上发明构思,本发明提出一种基于运行参数关联性分析的工业过程异常工况预测方法,包括如下步骤:
步骤1:采用Holt指数平滑模型,对工业过程各传感器所采集的测量值序列进行定步长预测;
步骤2:以滑动窗口的方式,构造以工业过程运行参数变化趋势为特征的三元组,对数据窗口内的测量值序列进行表示,并通过基于欧氏距离的关联性指标计算该窗口内任意两个运行参数的关联性;
步骤3:根据已计算获得的关联性数据构造关联性预测模型,所述的关联性预测模型为多元自回归模型,并通过偏最小二乘算法估计模型参数;
步骤4:对于新获得数据,根据关联性预测模型进行设备异常工况的预测,并在预测到异常工况发生前不断更新关联性预测模型以及模型参数。
基于上述方案,各步骤可具体采用如下实现方式:
作为优选,步骤1具体如下:
步骤1.1:针对具有多个传感器的工业设备,记传感器数量为N,当设备处于运行过程中,不断采集到表征设备运行状态的运行参数值即传感器数据,将各传感器测量序列记为其中K表示序列长度,表示传感器i在第k个采样时刻点的测量值;
步骤1.2:对于传感器i,采用Holt指数平滑模型对其测量值序列进行预测,给定测量值其平滑值可根据下式计算:
其中,表示的平滑值;是线性增长因子,代表平滑后的趋势;α和β是平滑系数,取值范围均为(0,1);Holt模型的初始条件如下:
步骤1.3:在获得测量数据的平滑值以及线性增长因子后,利用该结果进行预测,预测值为:
其中l代表预测步长,τ为设备传感器信号的采样间隔。
作为优选,步骤2包含以下子步骤:
步骤2.1:运行参数关联性分析分为时间序列分段拟合、三元组表示和关联性计算三个阶段;时间序列分段拟合阶段中,对于固定长度为L的数据窗口,记为窗口Wj,运行参数Xi,i=1,2,…,N,在该窗口内的L个数据为其中j为该窗口起点,对应的采样时刻为tj;若假设该窗口内其中一个时间段内的m个数据能够由一条线段进行拟合,对于其中的测量值其在拟合线段中所对应的值为则该线段的拟合误差ERR计算公式为:
在对于以j为起点的窗口数据进行分段线性化时,从开始对其进行线段拟合;在数据分段拟合的过程按以下步骤进行,其中记设定的拟合误差阈值为ωE:
步骤(1):设置拟合起点为拟合终点为h=2;其中,对于以j为起点的窗口Wj,初始的拟合起点为
步骤(2):对于数据采用直线回归的方式进行线段拟合,由此获得相应的线段数据根据所述的拟合误差ERR计算公式计算其拟合误差ERR;
步骤(3):若ERR≤ωE,则令h=h+1,并重复步骤(2);若ERR>ωE,则保存当前拟合终点(即数据分割点),重置h=2,并回到步骤1,以当前拟合终点为新的拟合起点进行下一部分数据的拟合;
重复以上步骤(1)~(3)直至该窗口Wj内所有数据均已线段化,即获得了分段线性化后的数据
步骤2.2:在线段三元组表示阶段,采用以下三元组的形式来描述一条线段sj,
其中,kj表示线段斜率,表示该线段在时间轴上的长度,rj表示该线段数值的增长率,即对于线段数据对于窗口Wj数据线段化后的新序列得到其三元组序列表示形式为{s1,s2,…,sn},其中n表示该窗口内数据分段后的线段数量;
步骤2.3:在关联性计算阶段,对于设备运行过程中的两个运行参数VA和运行参数VB,首先需对线段化后的数据进行分割:在窗口Wj内,记参数VA的分段点为参数VB的分段点为其中nA和nB分别表示运行参数A和B的测量数据在该窗口内分段后的线段数;对参数VA和VB的分段点进行合并,去除重复项并进行从小到大排列,得到两个参数的分割点序列为随后,根据该分割点序列以及三元组表示形式获得参数VA和VB新的三元组序列为和
获得窗口Wj内两个参数VA和VB的三元组序列后,通过基于欧氏距离的关联性指标dAB来计算该窗口内两个参数VA和VB的关联性:
式中:为参数VA第i条线段的斜率,为参数VB第i条线段的斜率,为参数VA第i条线段的数值增长率,为参数VB第i条线段的数值增长率。
作为优选,步骤3包含以下子步骤:
步骤3.1:构造关联性预测模型:设置预测步长为f,根据已知数据长度确定模型设计参数U和M使U+f+M-1=窗口Wj的总个数,并构造以下矩阵:
其中,{d1,d2,…,dU+f+M-1}表示参数VA和VB的关联性序列,dj表示窗口Wj内两个参数VA和VB的关联性指标dAB;
步骤3.2:对于每一关联性序列,使用其中U个关联性值预测f步长后的关联性值,构造关联性预测模型:
Fp=Dpθ
其中,参数θ=[θ1,θ2,…,θU]T,可由偏最小二乘算法获得。
作为优选,步骤4包含以下子步骤:
在预测时,利用传感器新获得的数据构造新的矩阵:
随后,使用已构造关联性预测模型对关联性值进行f步预测,即
当判断设备出现异常,其中dnormal为设备运行初期处于正常运行状态时参数VA和VB的关联性值,ωp为关联性值相对于初始正常的漂移量阈值;若则利用传感器获得的新数据重构模型以更新模型参数θ;
随着数据的更新不断预测一定预测步长的关联性值,从而预测设备异常工况发生时间。
本发明提出的基于运行参数关联性分析的工业过程异常工况预测方法,可用于具备多个传感器的复杂工业***。本发明从工业过程运行参数之间的关联性入手,基于运行参数的关联性分析进行异常工况预测。在单参数预测阶段,本发明根据已有传感器数据通过指数平滑方法对各个运行参数进行预测。在关联性分析阶段,本发明通过已知的各运行参数值以及参数预测值计算运行参数关联性,其中参数关联性由代表参数曲线的一系列指标的相似性表示。在关联性趋势预测阶段,本发明构造多元自回归模型对参数关联性进行预测。本发明提出的方法将运行参数关联性纳入考虑,能够获得更完备的设备异常信息以及更提前的预测结果。这将给后续的设备健康管理提供强有力的数据支撑,对于高可靠性的设备维护管理尤有价值,在实际工程应用方面具有广阔前景。
附图说明
图1汽轮机运行参数测量值与预测结果;
图2汽轮机真空A与其他参数关联性趋势预测结果与真值对比;
图3汽轮机异常发生时间预测结果。
具体实施方式
现结合附图对本发明的具体实施方式作进一步的说明,部分原理已在前面详细叙述,在此不再赘述。下面本例用一个基于汽轮机低真空保护跳机数据的真实案例来阐述具体操作步骤以及验证所提出方法的有效性。
该铣床数据记录了采用铣刀切削金属材料的运行退化过程。该汽轮机运行的初始工况为负荷250MW、凝汽器真空93kPa,以凝汽器真空值作为指示参数,真空A从第762个采样点开始指示异常,当真空值下降至81kPa时,该汽轮机跳机。工业过程异常工况预测方法包括以下步骤:
步骤1:采用Holt指数平滑模型,对工业过程各传感器所采集的测量值序列进行定步长预测。本步骤具体包含以下子步骤:
步骤1.1:针对具有多个传感器的工业设备,记传感器数量为N,当设备处于运行过程中,不断采集到表征设备运行状态的运行参数值即传感器数据,将各传感器测量序列记为其中K表示序列长度,表示传感器i在第k个采样时刻点的测量值;
步骤1.2:对于传感器i,采用Holt指数平滑模型对其测量值序列进行预测,给定测量值其平滑值可根据下式计算:
其中,表示的平滑值;是线性增长因子,代表平滑后的趋势;α和β是平滑系数,取值范围均为(0,1);Holt模型的初始条件如下:
步骤1.3:在获得测量数据的平滑值以及线性增长因子后,利用该结果进行预测,预测值为:
其中l代表预测步长,τ为设备传感器信号的采样间隔。在本例中,预设的预测步长为l=15。
根据步骤1,对于每一运行参数测量序列,进行定步长预测,结果在图1给出,与此同时还给出了实际的状态监测测量序列。
步骤2:以滑动窗口的方式,构造以工业过程运行参数变化趋势为特征的三元组,对数据窗口内的测量值序列进行表示,并通过基于欧氏距离的关联性指标计算该窗口内任意两个运行参数的关联性。本步骤具体包含以下子步骤:
步骤2.1:运行参数关联性分析分为时间序列分段拟合、三元组表示和关联性计算三个阶段;时间序列分段拟合阶段中,对于固定长度为L的数据窗口,记为窗口Wj,在本例中,滑动窗口长度为100。运行参数Xi,i=1,2,…,N,在该窗口内的L个数据为其中j为该窗口起点,对应的采样时刻为tj;若假设该窗口内其中一个时间段内的m个数据恰好能够由一条线段进行拟合,对于其中的测量值其在拟合线段中所对应的值为则该线段的拟合误差ERR计算公式为:
在对于以j为起点的窗口数据进行分段线性化时,从开始对其进行线段拟合;在数据分段拟合的过程按以下步骤进行(其中记设定的拟合误差阈值为ωE,本例中ωE=0.025):
步骤(1):设置拟合起点为拟合终点为h=2;其中,对于以j为起点的窗口Wj,初始的拟合起点为
步骤(2):对于数据采用直线回归的方式进行线段拟合,由此获得相应的线段数据根据所述的拟合误差ERR计算公式计算其拟合误差ERR;
步骤(3):若ERR≤ωE,则令h=h+1,并重复步骤(2);若ERR>ωE,则保存当前拟合终点(即数据分割点),重置h=2,并回到步骤1,以当前拟合终点为新的拟合起点进行下一部分数据的拟合;
重复以上步骤(1)~(3)直至该窗口Wj内所有数据均已线段化,即获得了分段线性化后的数据
步骤2.2:在线段三元组表示阶段,采用以下三元组的形式来描述一条线段sj,
其中,kj表示线段斜率,表示该线段在时间轴上的长度,rj表示该线段数值的增长率,即对于线段数据对于窗口Wj数据线段化后的新序列得到其三元组序列表示形式为{s1,s2,…,sn},其中n表示该窗口内数据分段后的线段数量;
步骤2.3:在关联性计算阶段,对于设备运行过程中的两个运行参数VA和运行参数VB,首先需对线段化后的数据进行分割:以运行参数VA和运行参数VB为例,在窗口Wj内,记参数VA的分段点为参数VB的分段点为其中nA和nB分别表示运行参数A和B的测量数据在该窗口内分段后的线段数;对参数VA和VB的分段点进行合并,去除重复项并进行从小到大排列,得到两个参数的分割点序列为随后,根据该分割点序列以及三元组表示形式获得参数VA和VB新的三元组序列为和
获得窗口Wj内两个参数VA和VB的三元组序列后,通过基于欧氏距离的关联性指标dAB来计算该窗口内两个参数VA和VB的关联性:
式中:为参数VA第i条线段的斜率,为参数VB第i条线段的斜率,为参数VA第i条线段的数值增长率,为参数VB第i条线段的数值增长率。
步骤3:根据已计算获得的关联性数据构造关联性预测模型,所述的关联性预测模型为多元自回归模型,并通过偏最小二乘算法估计模型参数。本步骤具体包含以下子步骤:
步骤3.1:构造关联性预测模型:设置预测步长为f,根据已知数据长度确定模型设计参数U和M使U+f+M-1=窗口Wj的总个数,并构造以下矩阵:
其中,{d1,d2,…,dU+f+M-1}表示参数VA和VB的关联性序列,dj表示窗口Wj内两个参数VA和VB的关联性指标dAB;
步骤3.2:对于每一关联性序列,使用其中U个关联性值预测f步长后的关联性值,构造关联性预测模型:
Fp=Dpθ
其中,参数θ=[θ1,θ2,…,θU]T,可由偏最小二乘算法获得。
步骤4:对于新获得数据,根据关联性预测模型进行设备异常工况的预测,并在预测到异常工况发生前不断更新关联性预测模型以及模型参数。本步骤具体包含以下子步骤:
在预测时,利用传感器新获得的数据构造新的矩阵:
随后,使用已构造关联性预测模型对关联性值进行f步预测,即
当判断设备出现异常,其中dnormal为设备运行初期处于正常运行状态时参数VA和VB的关联性值,ωp为关联性值相对于初始正常的漂移量阈值;若则利用传感器获得的新数据重构模型以更新模型参数θ;
随着数据的更新不断预测一定预测步长的关联性值,从而预测设备异常工况发生时间。在本例中,预测步长f=10,阈值ωp=0.1。
根据步骤2-步骤4,进行关联性计算以及关联性变化趋势预测,结果如图2所示。随后,根据设置的失效阈值ωp,进行异常工况预测,结果如图3所示,与此同时还给出了使用实际数据获得的关联性序列根据设定的阈值检测得到的失效时间。
图1给出了汽轮机运行参数测量值与预测结果。从图中可知,利用指数平滑预测便可对单个参数获得较好的预测效果。图2给出了汽轮机真空A与其他参数关联性趋势预测结果与真值对比。由图可见,得到了良好的预测效果。图3给出了汽轮机异常发生时间预测结果。图3横坐标的1-6分别代表真空A与其余6个参数的关联性,纵坐标表示预测的异常工况发生时间。由图3可知,利用几组参数关联性的预测都获得了比较准确的预测结果,即可以提前一定步长预测到异常工况的发生。另外,通过对原始数据的分析,已知真空A从第762个采样点开始指示异常,加速下降,而通过对关联性变化趋势的分析,从图3可以看出,最早预测到异常工况发生的时间为第590个采样点(ID6),即更早捕捉到了异常的发生。更重要的是,通过关联性变化趋势,找到最早发生异常所对应的参数,有助于判断异常发生的位置,对于排除异常起到重要的作用。
Claims (5)
1.一种基于运行参数关联性分析的工业过程异常工况预测方法,其特征在于,包括以下步骤:
步骤1:采用Holt指数平滑模型,对工业过程各传感器所采集的测量值序列进行定步长预测;
步骤2:以滑动窗口的方式,构造以工业过程运行参数变化趋势为特征的三元组,对数据窗口内的测量值序列进行表示,并通过基于欧氏距离的关联性指标计算该窗口内任意两个运行参数的关联性;
步骤3:根据已计算获得的关联性数据构造关联性预测模型,所述的关联性预测模型为多元自回归模型,并通过偏最小二乘算法估计模型参数;
步骤4:对于新获得数据,根据关联性预测模型进行设备异常工况的预测,并在预测到异常工况发生前不断更新关联性预测模型以及模型参数。
2.根据权利要求1所述的一种基于运行参数关联性分析的工业过程异常工况预测方法,其特征在于:步骤1包含以下子步骤:
步骤1.1:针对具有多个传感器的工业设备,记传感器数量为N,当设备处于运行过程中,不断采集到表征设备运行状态的运行参数值即传感器数据,将各传感器测量序列记为其中K表示序列长度,表示传感器i在第k个采样时刻点的测量值;
步骤1.2:对于传感器i,采用Holt指数平滑模型对其测量值序列进行预测,给定测量值其平滑值可根据下式计算:
其中,表示的平滑值;是线性增长因子,代表平滑后的趋势;α和β是平滑系数,取值范围均为(0,1);Holt模型的初始条件如下:
步骤1.3:在获得测量数据的平滑值以及线性增长因子后,利用该结果进行预测,预测值为:
其中l代表预测步长,τ为设备传感器信号的采样间隔。
3.根据权利要求2所述的一种基于运行参数关联性分析的工业过程异常工况预测方法,其特征在于:步骤2包含以下子步骤:
步骤2.1:运行参数关联性分析分为时间序列分段拟合、三元组表示和关联性计算三个阶段;时间序列分段拟合阶段中,对于固定长度为L的数据窗口,记为窗口Wj,运行参数Xi,i=1,2,…,N,在该窗口内的L个数据为其中j为该窗口起点,对应的采样时刻为tj;若假设该窗口内其中一个时间段内的m个数据能够由一条线段进行拟合,对于其中的测量值其在拟合线段中所对应的值为则该线段的拟合误差ERR计算公式为:
在对于以j为起点的窗口数据进行分段线性化时,从开始对其进行线段拟合;在数据分段拟合的过程按以下步骤进行,其中记设定的拟合误差阈值为ωE:
步骤(1):设置拟合起点为拟合终点为其中,对于以j为起点的窗口Wj,初始的拟合起点为
步骤(2):对于数据采用直线回归的方式进行线段拟合,由此获得相应的线段数据根据所述的拟合误差ERR计算公式计算其拟合误差ERR;
步骤(3):若ERR≤ωE,则令h=h+1,并重复步骤(2);若ERR>ωE,则保存当前拟合终点,重置h=2,并回到步骤1,以当前拟合终点为新的拟合起点进行下一部分数据的拟合;
重复以上步骤(1)~(3)直至该窗口Wj内所有数据均已线段化,即获得了分段线性化后的数据
步骤2.2:在线段三元组表示阶段,采用以下三元组的形式来描述一条线段sj,
其中,kj表示线段斜率,表示该线段在时间轴上的长度,rj表示该线段数值的增长率,即对于线段数据对于窗口Wj数据线段化后的新序列得到其三元组序列表示形式为{s1,s2,…,sn},其中n表示该窗口内数据分段后的线段数量;
步骤2.3:在关联性计算阶段,对于设备运行过程中的两个运行参数VA和运行参数VB,首先需对线段化后的数据进行分割:在窗口Wj内,记参数VA的分段点为参数VB的分段点为其中nA和nB分别表示运行参数A和B的测量数据在该窗口内分段后的线段数;对参数VA和VB的分段点进行合并,去除重复项并进行从小到大排列,得到两个参数的分割点序列为随后,根据该分割点序列以及三元组表示形式获得参数VA和VB新的三元组序列为和
获得窗口Wj内两个参数VA和VB的三元组序列后,通过基于欧氏距离的关联性指标dAB来计算该窗口内两个参数VA和VB的关联性:
式中:为参数VA第i条线段的斜率,为参数VB第i条线段的斜率,为参数VA第i条线段的数值增长率,为参数VB第i条线段的数值增长率。
4.根据权利要求3所述的一种基于运行参数关联性分析的工业过程异常工况预测方法,其特征在于:步骤3包含以下子步骤:
步骤3.1:构造关联性预测模型:设置预测步长为f,根据已知数据长度确定模型设计参数U和M使U+f+M-1=窗口Wj的总个数,并构造以下矩阵:
其中,{d1,d2,…,dU+f+M-1}表示参数VA和VB的关联性序列,dj表示窗口Wj内两个参数VA和VB的关联性指标dAB;
步骤3.2:对于每一关联性序列,使用其中U个关联性值预测f步长后的关联性值,构造关联性预测模型:
Fp=Dpθ
其中,参数θ=[θ1,θ2,…,θU]T,可由偏最小二乘算法获得。
5.根据权利要求4所述的一种基于运行参数关联性分析的工业过程异常工况预测方法,其特征在于:步骤4包含以下子步骤:
在预测时,利用传感器新获得的数据构造新的矩阵:
随后,使用已构造关联性预测模型对关联性值进行f步预测,即
当判断设备出现异常,其中dnormal为设备运行初期处于正常运行状态时参数VA和VB的关联性值,ωp为关联性值相对于初始正常的漂移量阈值;若则利用传感器获得的新数据重构模型以更新模型参数θ;
随着数据的更新不断预测一定预测步长的关联性值,从而预测设备异常工况发生时间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910244872.5A CN110008565B (zh) | 2019-03-28 | 2019-03-28 | 基于运行参数关联性分析的工业过程异常工况预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910244872.5A CN110008565B (zh) | 2019-03-28 | 2019-03-28 | 基于运行参数关联性分析的工业过程异常工况预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110008565A true CN110008565A (zh) | 2019-07-12 |
CN110008565B CN110008565B (zh) | 2020-12-11 |
Family
ID=67168700
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910244872.5A Active CN110008565B (zh) | 2019-03-28 | 2019-03-28 | 基于运行参数关联性分析的工业过程异常工况预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110008565B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111553606A (zh) * | 2020-05-06 | 2020-08-18 | 上海核工程研究设计院有限公司 | 一种对核电设备组件设计与设备指标关联性的评估方法 |
CN111913444A (zh) * | 2019-09-07 | 2020-11-10 | 宁波大学 | 一种基于时序多块建模策略的化工过程监测方法 |
CN111983478A (zh) * | 2020-07-07 | 2020-11-24 | 江苏大学 | 一种基于Holt线性趋势模型的电化学储能电站SOC异常检测方法 |
CN112464146A (zh) * | 2020-12-03 | 2021-03-09 | 北京航空航天大学 | 一种基于历史遥测数据的关键分***及单机关联性健康基线构建方法 |
CN112801426A (zh) * | 2021-04-06 | 2021-05-14 | 浙江浙能技术研究院有限公司 | 一种基于关联参数挖掘的工业过程故障融合预测方法 |
CN113624025A (zh) * | 2021-08-18 | 2021-11-09 | 浙江大学 | 基于运行参数相关性的凝汽器真空低跳机征兆捕捉方法 |
CN113792988A (zh) * | 2021-08-24 | 2021-12-14 | 河北先河环保科技股份有限公司 | 一种企业在线监测数据异常识别方法 |
CN113867262A (zh) * | 2020-06-30 | 2021-12-31 | 华晨宝马汽车有限公司 | 用于监测主轴设备的运行状态的装置、方法以及介质 |
CN115034094A (zh) * | 2022-08-10 | 2022-09-09 | 南通恒强轧辊有限公司 | 一种金属加工机床运行状态预测方法及*** |
CN115509187A (zh) * | 2022-09-20 | 2022-12-23 | 北京中佳瑞通科技有限公司 | 一种工业大数据处理方法和*** |
CN115689071A (zh) * | 2023-01-03 | 2023-02-03 | 南京工大金泓能源科技有限公司 | 一种基于关联参数挖掘的设备故障融合预测方法及*** |
CN115829422A (zh) * | 2023-02-21 | 2023-03-21 | 创银科技(南通)有限公司 | 基于大数据的工业设备运行异常状态识别方法 |
WO2023094916A1 (en) * | 2021-11-23 | 2023-06-01 | International Business Machines Corporation | Identifying persistent anomalies for failure prediction |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090265155A1 (en) * | 2008-04-21 | 2009-10-22 | Nec Electronics Corporation | Method of predicting reliability of semiconductor device, reliability prediction system using the same and storage medium storing program causing computer to execute the same |
CN108037387A (zh) * | 2017-12-05 | 2018-05-15 | 北京能源集团有限责任公司 | 基于聚类归集的设备故障分析方法及装置 |
CN108873859A (zh) * | 2018-05-31 | 2018-11-23 | 浙江工业大学 | 基于改进关联规则的桥式抓斗卸船机故障预测模型方法 |
CN109086804A (zh) * | 2018-07-12 | 2018-12-25 | 中石化石油机械股份有限公司 | 一种基于多源状态监测信息和可靠性特征融合的液压设备早期故障预测方法 |
CN109110608A (zh) * | 2018-10-25 | 2019-01-01 | 歌拉瑞电梯股份有限公司 | 一种基于大数据学习的电梯故障预测方法 |
-
2019
- 2019-03-28 CN CN201910244872.5A patent/CN110008565B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090265155A1 (en) * | 2008-04-21 | 2009-10-22 | Nec Electronics Corporation | Method of predicting reliability of semiconductor device, reliability prediction system using the same and storage medium storing program causing computer to execute the same |
CN108037387A (zh) * | 2017-12-05 | 2018-05-15 | 北京能源集团有限责任公司 | 基于聚类归集的设备故障分析方法及装置 |
CN108873859A (zh) * | 2018-05-31 | 2018-11-23 | 浙江工业大学 | 基于改进关联规则的桥式抓斗卸船机故障预测模型方法 |
CN109086804A (zh) * | 2018-07-12 | 2018-12-25 | 中石化石油机械股份有限公司 | 一种基于多源状态监测信息和可靠性特征融合的液压设备早期故障预测方法 |
CN109110608A (zh) * | 2018-10-25 | 2019-01-01 | 歌拉瑞电梯股份有限公司 | 一种基于大数据学习的电梯故障预测方法 |
Non-Patent Citations (3)
Title |
---|
CHUNHUI ZHAO 等: "Critical-to-Fault-Degradation Variable Analysis and Direction Extraction for Online Fault Prognostic", 《IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY》 * |
卢会国 等: "基于关联规则的故障诊断方法研究", 《成都信息工程学院学报》 * |
方名菊 等: "基于形式概念的汽轮机相似故障案例检索", 《山东科技大学学报(自然科学版)》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111913444A (zh) * | 2019-09-07 | 2020-11-10 | 宁波大学 | 一种基于时序多块建模策略的化工过程监测方法 |
CN111913444B (zh) * | 2019-09-07 | 2022-03-18 | 宁波大学 | 一种基于时序多块建模策略的化工过程监测方法 |
CN111553606A (zh) * | 2020-05-06 | 2020-08-18 | 上海核工程研究设计院有限公司 | 一种对核电设备组件设计与设备指标关联性的评估方法 |
CN113867262A (zh) * | 2020-06-30 | 2021-12-31 | 华晨宝马汽车有限公司 | 用于监测主轴设备的运行状态的装置、方法以及介质 |
CN111983478A (zh) * | 2020-07-07 | 2020-11-24 | 江苏大学 | 一种基于Holt线性趋势模型的电化学储能电站SOC异常检测方法 |
CN112464146A (zh) * | 2020-12-03 | 2021-03-09 | 北京航空航天大学 | 一种基于历史遥测数据的关键分***及单机关联性健康基线构建方法 |
CN112801426A (zh) * | 2021-04-06 | 2021-05-14 | 浙江浙能技术研究院有限公司 | 一种基于关联参数挖掘的工业过程故障融合预测方法 |
CN113624025A (zh) * | 2021-08-18 | 2021-11-09 | 浙江大学 | 基于运行参数相关性的凝汽器真空低跳机征兆捕捉方法 |
CN113792988A (zh) * | 2021-08-24 | 2021-12-14 | 河北先河环保科技股份有限公司 | 一种企业在线监测数据异常识别方法 |
WO2023094916A1 (en) * | 2021-11-23 | 2023-06-01 | International Business Machines Corporation | Identifying persistent anomalies for failure prediction |
CN115034094A (zh) * | 2022-08-10 | 2022-09-09 | 南通恒强轧辊有限公司 | 一种金属加工机床运行状态预测方法及*** |
CN115509187A (zh) * | 2022-09-20 | 2022-12-23 | 北京中佳瑞通科技有限公司 | 一种工业大数据处理方法和*** |
CN115509187B (zh) * | 2022-09-20 | 2023-04-18 | 北京中佳瑞通科技有限公司 | 一种工业大数据处理方法和*** |
CN115689071A (zh) * | 2023-01-03 | 2023-02-03 | 南京工大金泓能源科技有限公司 | 一种基于关联参数挖掘的设备故障融合预测方法及*** |
CN115689071B (zh) * | 2023-01-03 | 2023-05-02 | 南京工大金泓能源科技有限公司 | 一种基于关联参数挖掘的设备故障融合预测方法及*** |
CN115829422A (zh) * | 2023-02-21 | 2023-03-21 | 创银科技(南通)有限公司 | 基于大数据的工业设备运行异常状态识别方法 |
CN115829422B (zh) * | 2023-02-21 | 2024-01-02 | 北京瀚海蓝山智能科技有限公司 | 基于大数据的工业设备运行异常状态识别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110008565B (zh) | 2020-12-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110008565A (zh) | 一种基于运行参数关联性分析的工业过程异常工况预测方法 | |
CN104390657B (zh) | 一种发电机组运行参数测量传感器故障诊断方法及*** | |
CN110008253B (zh) | 一种工业数据关联规则挖掘及异常工况预测方法 | |
CN110414154B (zh) | 一种带有双测点的风机部件温度异常检测和报警方法 | |
CN104483575A (zh) | 用于非侵入式电力监测的自适应负荷事件检测方法 | |
EP3367198A1 (en) | Method for monitoring and diagnosing rotating machines based on signature analysis of electrical signals | |
CN102789676B (zh) | 一种基于报警证据融合的工业报警器设计方法 | |
CN113776794B (zh) | 一种嵌入式边缘计算的故障诊断方法、装置和*** | |
CN108038044A (zh) | 一种面向连续被监测对象的异常检测方法 | |
CN112418517B (zh) | 一种用于冻害短临预警的预测方法 | |
CN111860839A (zh) | 基于多信号融合及Adam优化算法的岸桥故障监测方法 | |
CN111665066A (zh) | 基于卷积神经网络的设备故障自适应上下预警界生成方法 | |
CN115455358A (zh) | 基于非线性回归模型的电气参数趋势预警及故障诊断方法 | |
CN111984705A (zh) | 一种电能质量监测装置的精度评估管理方法和装置 | |
CN117556366B (zh) | 基于数据筛选的数据异常检测***及方法 | |
WO2016136391A1 (ja) | 事故点標定装置及び方法、電力系統監視システム、設備計画支援システム | |
CN109933031B (zh) | 一种根据化验数据自动校正软测量仪表的***及方法 | |
US11339763B2 (en) | Method for windmill farm monitoring | |
CN109556861A (zh) | 一种基于案例推理的轴承实时故障诊断*** | |
KR102351538B1 (ko) | 패턴 학습 및 진단 알고리즘을 이용한 자동화 설비의 고장 예측 시스템 및 그 방법 | |
CN111986469A (zh) | 一种现场终端故障智能化诊断的方法 | |
WO2022071189A1 (ja) | 状態監視装置および状態監視方法 | |
CN109583124A (zh) | 一种基于adrc的hmm故障预测*** | |
CN118013468B (zh) | 一种风电机组部件健康度监测方法、装置、设备及介质 | |
CN115127605B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |