法律状态公告日
法律状态信息
法律状态
2020-03-31
授权
授权
2019-08-09
著录事项变更 IPC(主分类):G06F17/50 变更前: 变更后: 申请日:20180622
著录事项变更
2018-09-18
实质审查的生效 IPC(主分类):G06F17/50 申请日:20180622
实质审查的生效
2018-08-24
公开
公开
技术领域
本发明涉及非高斯随机系统故障诊断领域,具体为基于广义互熵—DPCA算法的有机朗肯循环系统故障诊断方法。
背景技术
节能减排,提高能源的利用率成为保持可持续发展的重要措施。低温余热发电系统将锅炉尾部烟道中的烟气余热通过有机朗肯循环(ORC,Organic Rankine Cycle)转换为机械能,进而转化为高品位电能,在节能、节水、减少有害气体排放等方面具有重要意义。随着系统长时间的运行,系统中的各个部件可能发生故障,导致系统效率下降、性能恶化,造成巨大的经济损失,甚至严重威胁人身安全。因此,建立ORC过程的状态监测与故障诊断系统来保证系统的安全、稳定、有效运行越来越得到重视。
在有机朗肯循环系统中,随机干扰是不可避免的且随机噪声未必服从高斯分布,仅仅采用均值和方差为不能全面反映系统输出或跟踪误差的高阶统计特性。因此采用非高斯随机控制理论来研究基于有机朗肯循环系统的故障诊断将具有更具有一般性的研究意义。
由于ORC系统不服从高斯分布且具有较多的变量,采用传统的DPCA方法(动态主元分析Dynamic PCA,DPCA)会出现误报和漏报的现象,本发明针对这一缺点,在原有主元分析的基础上,采用广义互熵准则对传统的DPCA进行改进,以解决传统DPCA算法的不足。
发明内容
本发明为了解决由于ORC系统不服从高斯分布且具有较多的变量,因此采取传统的DPCA方法会出现误报和漏报的问题,提供了一种基于广义互熵—DPCA算法的有机朗肯循环系统故障诊断方法。
本发明是通过如下技术方案来实现的:一种基于广义互熵—DPCA算法的有机朗肯循环系统的故障诊断方法,包括如下步骤:
步骤一:采集正常工况下的数据X∈RN×m作为正常数据,并且对其进行归一化处理,其中N表示样本数,m表示变量个数;
步骤二:确定时间滞后长度l,构建具有动态信息的増广矩阵,具体过程如下:
①设定l=0;
②通过执行广义互熵—DPCA算法,实现数据降维,计算所有的主元得分,即依据重构误差建立性能指标并通过优化该性能指标求得最优的方向矩阵,实现故障诊断,包括以下步骤:
a.计算平均重构误差
即
其中xi为任一样本点,
b.构建广义互熵—DPCA性能指标
互熵描述两个随机变量的随机性,互熵的具体表达式为:
其中:X与
广义高斯核函数的表达式为:
将广义高斯核函数代替高斯核函数后,重新定义后的互熵的表达式为:
其中α>0,β>0分别为核参数,Γ为伽玛函数;
互熵为正并且有边界,其边界表达式为:
当互熵达到最大值时,X与
最终得出:
其中
c.采用梯度下降法对性能指标进行优化,得到所有方向的方向矩阵P:
PCA的本质可以看作是一种坐标变换,由上述建立广义互熵—DPCA算法,当重构误差最小时说明我们已经把原始的监测数据不失真的转换到新的坐标系中。广义互熵是两个随机变量的相似程度的度量,即当广义互熵最小时,检测数据X与重构数据X具有较强的相似性。因此,通过优化性能指标即可得到最优的方向矩阵P。由于性能指标(8)是带约束的优化问题,因此在本发明中引入拉格朗日因子使得问题转换成无约束的参数优化问题,如式(9)所示:
优化性能指标有很多方法,在本发明中采用梯度下降法对性能指标进行优化;梯度下降法分为两步实现:
c1.分别对方向矩阵P和拉格朗日因子λ求偏倒;
c2.依据递推公式求得未知参数P和λ:
其中μ和η为学习率,且μ∈[0,1],η∈[0,1];其作用为平衡偏倒的部分在幅值上对所求最小值的影响(即平衡
d.采用交叉验证方法确定主元个数,实现数据降维:
主元个数的确定有很多方法,如主元累计方差贡献百分比法,未重构方差,交叉验证等。由于百分比法和未重构方差方法都含有较多的主观因素,交叉验证主要依据数据本身的特性提取主元,因此本发明采用交叉验证方法求取主元,具体步骤如下:
d1.分裂正常数据XN×m即X∈RN×m为两部分,其中一部分为
d2.对训练数据
d3.计算测试数据
d4.计算重构误差e=xi-ti(P-i)T;重复执行d1—d4,当得到所有的残差后停止循环;
d5.计算
③设定j=m×(l+1)和r(l)=0,判断第j个主元得分是否小于设定值,如果是,执行④;如果否,执行⑤;
④j=j-1,r(l)=r(l)+1,重复执行③直到j=0;
⑤计算新的动态关系:
⑥如果rnew(l)≤0,停止;否则l=l+1,跳转到②;
⑦对于正常数据X∈RN×m,根据上述步骤得到的时间延迟长度l构建具有动态信息的增广矩阵,其具体模型为:
其中xT(i)表示第i时刻的观测值;
步骤三:针对增广数据矩阵,运用步骤二中的a—d步骤计算得到方向矩阵P;
步骤四:设置置信度,利用核密度估计方法估计SPE统计量以及T2统计量的概率密度函数,通过对两个统计量的概率密度函数积分得到故障检测控制限,包括如下步骤:
1)SPE表示每个采样数据与统计模型之间的误差,具体的计算为:
SPE=||(I-PPT)x||2>
2)T2统计量表示得分向量的平方和,其具体计算为:
T2=xTPPTx>
3)由于数据是非高斯分布,而传统求控制限的方法是在高斯假设的前提下,并不适应用非高斯情况,因此利用核密度估计方法估计SPE和T2的概率密度函数:
其中λ∈[0,1],
4):设定置信度为α,对概率密度函数积分求得T2的控制限
步骤五:采集故障工况下的数据作为故障数据,并且对其进行归一化处理,随后将故障数据利用方向矩阵计算得到故障数据模型的SPE统计量以及T2统计量,通过与已计算出的控制限对比,检测出系统是否发生故障。
与现有技术相比本发明具有以下有益效果:(1)工业过程中不可避免受到非高斯噪声的影响,本发明针对这一现象设计了广义互熵—DPCA算法对ORC系统进行故障诊断,且广义互熵采用广义高斯核,更具有一般性。
(2)与传统的求取控制限的方法不同,本发明充分考虑到非高斯噪声的影响,在计算控制限的过程中采用核密度估计方法,有效解决了传统DPCA算法在高斯情况的假设,提高了故障检测的精度。
附图说明
图1为广义互熵—DPCA算法的故障诊断结构框图。
图2为具体数值实施例加入故障时的仿真图。
具体实施方式
以下结合具体实施例对本发明作进一步说明。
步骤一:采集正常工况下的数据X∈RN×m作为正常数据,并且对其进行归一化处理,其中N表示样本数,m表示变量个数。
ORC系统主要由蒸发器,冷凝器,工质泵,膨胀机,阀门这六个主要部件构成。蒸发器模型中入口质量流量与焓值由工质泵提供,出口质量流量是膨胀机入口的质量流量。类似的,冷凝器入口质量流量与焓值由膨胀机提供,出口质量流量是泵的入口质量流量,进而将各个部件的模型串联起来,得到整个ORC系统的非线性模型如下式所示:
其中
由于系统中各变量的量纲不同导致数据差别很大,因此需要采用归一化处理训练数据,使得数据位于同一数量级,去均值处理后的监测数据为
表一为有机朗肯系统故障描述表,此表中列举了多种朗肯系统的故障分类,体现了实际应用中的故障类型。
表一:有机朗肯循环系统故障描述表
步骤二:确定时间滞后长度l,构建具有动态信息的増广矩阵。对于监测数据
①设定l=0;
②通过执行广义互熵—DPCA算法,实现数据降维,计算所有的主元得分,即依据重构误差建立性能指标并通过优化该性能指标求得最优的方向矩阵,实现故障诊断,包括以下步骤:
a.计算平均重构误差
即
其中xi为任一样本点,
b.根据广义互熵准则构建性能指标
c.利用梯度下降法优化性能指标:
d.依据交叉验证方法确定主元个数,num(num<11),则数据矩阵Xtrain∈R500×11→Xtrain∈R500×num,实现此时完成了广义高斯DPCA算法的构建。
③设定j=m×(l+1)和r(l)=0,利用a—d计算得到的方向矩阵P计算主元得分且判断第j个主元得分是否小于设定值,如果是,执行④;如果否,执行⑤;
④j=j-1,r(l)=r(l)+1,重复执行③直到j=0;
⑤计算新的动态关系:
⑥如果rnew(l)≤0,停止;否则l=l+1,跳转到②;
⑦对于正常数据X∈R500×11,根据上述步骤得到的时间延迟长度l,构建具有动态信息的增广矩阵,其具体模型为:
其中xT(i)表示第i时刻的观测值。
步骤三:针对增广数据矩阵X,运用步骤二中的a—d步骤计算得到方向矩阵P。
步骤四:设置置信度,利用核密度估计方法估计SPE统计量以及T2统计量的概率密度函数,通过对两个统计量的概率密度函数积分得到故障检测控制限。
1)SPE表示每个采样数据与统计模型之间的误差,其具体的计算为:
SPE=||(I-PPT)x||2>
2)T2统计量表示得分向量的平方和,其具体计算为:
T2=xTPPTx>
3):利用核密度估计方法估计SPE和T2的概率密度函数。
4):设定置信度为α,对概率密度函数积分求得T2的控制限
步骤五:采集故障工况下的数据作为故障数据,并且对其进行归一化处理,随后将故障数据利用方向矩阵计算得到故障数据模型的SPE统计量以及T2统计量,通过与已计算出的控制限对比,检测出系统是否发生故障。
表二:有机郎肯循环系统故障引入时间表
以第一种故障传感器恒增益故障为例:首先对采集到的具有传感器故障的数据进行归一化处理,随后确定时延参数l=3,因此构建的训练数据矩阵为
最后,随后将故障数据利用方向矩阵计算得到故障数据模型的SPE统计量以及T2统计量,比较SPE与ClSPE、T2与如果超过控制限则说明传感器发生故障,没有超过控制限则说明系统无故障。
例如,如图2所示,其中一个数据表达式如下:
X(t)=-R(t)3+3R(t)2+E(t)>
X∈R500×7,其中R(t)和E(t)分别为服从重尾分布的数组,在250s时刻,数据X的第二个变量以每秒0.5倍的速率降低作为故障数据,其故障诊断结果如图2所示,从SPE的结果可以看到,在250s超过了控制限,说明系统出现了故障。
机译: 有机朗肯循环系统,具有包括用于利用风扇驱动的有机朗肯循环工作流体回路的膨胀发动机,以及在有机朗肯循环过程中由工作单元驱动的给油泵和清除剂泵。
机译: 一个扩展器,一种有机朗肯循环系统,包括这种膨胀机和制造包括这种膨胀机的有机朗肯循环系统的方法
机译: 扩展器,一种有机朗肯循环系统,包括这种膨胀机和制造包括这种膨胀机的有机朗肯循环系统的方法