一种基于KDE-FA的冷水机组故障特征刻画方法与流程

一种基于kde
‑
fa的冷水机组故障特征刻画方法
技术领域
1.本发明属于空调系统中冷水机组故障检测领域，具体涉及一种基于非参数核密度估计和因子分析(kde
‑
fa)的冷水机组故障特征刻画方法。


背景技术：

2.冷水机组是暖通空调系统中的主要耗能设备，通过冷水机组的故障检测(fd)，及时发现和识别故障，并作出相应决策，对实现空调系统节能，保持空调系统高能效运行和维持室内环境舒适度具有重要的工程意义和经济意义。
3.故障特征刻画主要包括故障数据特征选择和特征降维。选择合适的故障特征是fd的前提和关键，特征选择的好坏很大程度上决定了故障检测的性能。故障特征维度较高及特征之间具有信息冗余性是现有fd方法研究所面临的主要问题之一，这会导致维度灾难，增加工作量，需要更多的时间和空间，从而降低故障检测效率。
4.大多数现有fd方法通常以其检测性能最佳为原则选择表征故障的特征参数，这往往导致获取这些特征参数的传感器数量多、成本高。冷水机组现场传感器安装现状的调研结果显示：温度传感器是现场存在最多且成本较低的传感器，压力和流量传感器成本较高且现场大量缺乏。因此选择大量压力或流量特征进行fd工作不具备经济意义，基于成本进行特征刻画成为fd方法的出发点。
5.目前，在冷水机组故障检测技术中，故障特征刻画的传统技术途径有：1)考虑现场成本和特征敏感性进行特征选择；2)基于pca等方法进行特征降维。然而，上述两条传统途径常常存在两大问题：1)不能完全从样本数据本身出发，反映故障特征的真实分布，从而导致对fd有用的敏感特征可能被遗漏；2)原始变量和公共因子之间的关系无法得到较好解释，提取出来的主成分无法清晰地解释其代表的含义。


技术实现要素：

6.为解决现有技术中存在的上述缺陷，本发明提出了一种基于kde
‑
fa的冷水机组故障特征刻画方法，根据该方法所确定的故障特征以及匹配的故障检测方法，在现场冷水机组故障检测应用中，既可以不遗漏有用的敏感特征，真实刻画故障数据特征，又可以降低特征维度，减少工作量，提高检测效率及精度，从而有效克服目前故障特征刻画存在的主要问题。
7.实现本发明的技术路径是：首先基于kde进行敏感特征选择，然后基于fa进行特征降维，从而达到在不影响检测精度的条件下，以最少特征表征冷水机组故障的目的；最后基于kde
‑
fa方法建立特征刻画模型，并将此模型用于bn和svdd的fd技术。其中，融入kde
‑
fa特征刻画方法的fd包括三个过程：1)数据预处理；2)特征刻画；3)故障检测过程。
8.本发明具体技术方案如下。
9.一种基于kde
‑
fa模型的冷水机组故障特征刻画方法，包括下述步骤：
10.步骤1：选用移动平均滤波器法对冷水机组原始故障数据进行平滑处理，然后选用
计算几何加权运行平均值和几何加权运行方差的方法对平滑数据进行稳态筛选；
11.步骤2：采用kde方法对稳态筛选后的故障数据进行故障特征选择：
12.步骤3：采用fa方法对kde选择出的冷水机组故障特征进行特征降维，完成故障特征刻画工作。采用fa方法对kde选择出的故障特征进行特征降维时，包括因子描述、因子提取、因子旋转和因子得分四个环节。
13.步骤4：构建kde
‑
bn模型、fa
‑
bn模型及kde
‑
fa
‑
bn模型，进行fd故障检测，对冷水机组的故障特征刻画的结果进行评价；
14.步骤5：构建kde
‑
svdd模型、fa
‑
svdd模型及kde
‑
fa
‑
svdd模型，进行fd故障检测，对冷水机组的故障特征刻画的结果进行评价；
15.步骤6：分别将构建的kde
‑
bn模型、fa
‑
bn模型及kde
‑
fa
‑
bn模型与传统bn模型的检测性能进行比较，若正常状态被误判为故障状态的概率越低，则判定该fd方法效果越好；
16.步骤7：分别将构建的kde
‑
svdd模型、fa
‑
svdd模型及kde
‑
fa
‑
svdd模型与传统svdd模型的检测性能进行比较，若正常状态被误判为故障状态的概率越低，则判定该fd方法效果越好。
17.本发明的有益效果是：
18.引入能反映故障特征真实分布的kde方法从原始特征集中尽可能选择出所有敏感特征；引入能对原始变量和公共因子之间进行解释的fa方法大幅降低特征维度，且去除特征间冗余度；提出一种基于非参数核密度估计和因子分析(kde
‑
fa)的冷水机组故障特征刻画方法，以解决以上存在的技术问题。
19.本发明先基于kde进行敏感特征选择，然后基于fa进行特征降维。在不影响检测精度的条件下，以最少特征表征冷水机组故障；然后基于贝叶斯网络(bn)和支持向量数据描述(svdd)两个模型的fd方法对特征选择结果进行评价。
20.本发明提出的一种基于kde
‑
fa的冷水机组故障特征刻画方法，kde方法基于敏感最大化的原则，尽量选出原始故障特征集中的所有敏感特征，相对于其他特征选择方法，所选出的敏感特征更多，保留信息更完整；fa方法去除了特征之间的冗余性，在保证不损失特征信息的基础上，大幅度降低特征维度，简化后续fd方法模型的结构复杂度；kde
‑
fa方法在实现特征维度降低的同时，保证fd方法模型的检测性能提高。
21.其中，kde是一种不对故障数据总体做任何假设，完全从故障样本数据本身出发，估计故障数据总体概率密度函数的方法。将kde应用于冷水机组fd的特征选择，能反映故障特征的真实分布，提高特征对故障类型的敏感分辨率；能刻画特征与故障类型之间的映射差异性。因此，该方法能够从原始故障特征集中尽可能选择出所有敏感特征。
22.fa是从故障特征中提取共性因子，根据数据的内在逻辑性将多个变量归纳为少数几个因子，从而进行降维处理，避免信息重叠。通过数学变换得到的新特征能够反映更多的原始信息；新特征和原始故障特征之间的可解释性较强，能够解释每个原始故障特征对每个新特征的贡献度。因此，该方法大幅降低特征维度，且去除特征间冗余度。
附图说明
23.图1为融入kde
‑
fa特征刻画方法的fd框架图；
24.图2为bn的结构图；
25.图3为svdd的原理图；
26.图4为碎石图；
27.图5为三维因子空间的载荷散点图；
28.图6为四种bn模型的far对比；
29.图7为传统svdd模型正常样本检测结果；
30.图8为kde
‑
svdd模型正常样本检测结果；
31.图9为fa
‑
svdd模型正常样本检测结果；
32.图10为kde
‑
fa
‑
svdd模型正常样本检测结果；
33.图11为四种svdd模型的far对比。
具体实施方式
34.下面结合附图和实施例对发明作进一步的详细说明，但并不作为对发明做任何限制的依据。
35.如图1所示，对本发明提出的一种基于kde
‑
fa的冷水机组故障特征刻画的具体实施步骤说明如下：
36.步骤1：本发明首先选用移动平均滤波器法对冷水机组原始故障数据进行平滑处理，然后选用计算几何加权运行平均值和几何加权运行方差的方法对平滑数据进行稳态筛选。
37.步骤2：采用kde方法对稳态筛选后的故障数据进行特征选择。
38.采用kde方法对冷水机组故障的稳态数据进行特征选择：
39.对于某一故障，设x1，x1，......x
n
是取自故障总体x的样本，在任意点x处的总体密度函数f(x)的kde采用下式计算：
[0040][0041]
其中，k()称为某一故障的核函数，h称为某一故障窗宽或带宽，n为故障数据的样本容量。
[0042]
核函数选择gaussian核函数，故障的最优带宽h
o
采用下式计算：
[0043][0044]
其中，s为故障数据的样本标准差。
[0045]
利用上述方法生成故障数据与正常数据的核密度对比图，然后基于尽量保证数据信息完整、尽量选择温度和压力参数以节约传感器成本两个原则，对核密度对比图中正常数据曲线与故障数据曲线的偏离度进行判断，最终选出故障的敏感特征，即为利用kde方法所选出来的特征。
[0046]
假设冷水机组有k类故障，每类故障有w个劣化等级，加上正常状况，则冷水机组共有wk 1种状况，假设每种状况可获得的特征参数有q个，均来自传感器。利用kde方法对冷水机组的wk 1种状况的q个特征参数进行估计，共得到q(wk 1)个概率密度函数(pdf)估计图和累积分布函数(cdf)估计图。然后针对每一类故障的每一个参数，将w个劣化等级数据的
pdf曲线与正常数据的pdf曲线进行对比，观察各劣化等级状况和正常状况之间的pdf曲线偏离度，进而判断此参数对此故障的敏感度，并作出决策是否选择此参数作为此故障的特征。
[0047]
步骤3：采用fa方法对kde选择出的冷水机组故障特征进行特征降维，完成故障特征工作。采用fa方法对kde选择出的故障特征进行特征降维时，包括因子描述、因子提取、因子旋转和因子得分四个环节：
[0048]
3a)在因子描述环节，进行bartlett球形检验和kmo统计量检验；当bartlett球形检验结果中的卡方统计量足够大、显著性值小于0.05、kmo统计量检验结果中的kmo大于0.5三个条件满足其一，则认为故障特征之间的冗余性较强，此时应采用fa方法进行特征降维；
[0049]
3b)在因子提取环节，采用方差贡献率法和碎石图法确定公共因子的提取个数；因子提取后输出的公共因子能够反映全部数据信息的95％时，选取此时对应的故障特征的公共因子个数；碎石图中以曲线由陡峭变平缓的转折点为判断点，选取陡峭曲线对应的故障特征的公共因子个数；
[0050]
3c)在因子旋转环节，当故障特征原始变量和提取出的公共因子之间的关系不确定时，采用正交旋转法对因子载荷矩阵进行旋转，使故障特征原始变量和公共因子两者之间的关系明确；对公共因子命名；
[0051]
3d)在因子得分环节，确定每一样本在公共因子上的具体数值，将其作为新的变量运用于下一步模型构建；选择回归法确定正常状态和各故障状态在因子变量的数值因子变量的数值采用下式计算：
[0052][0053]
其中，a为因子载荷矩阵，∑为标准化故障特征原始变量的协方差矩阵，x为故障特征原始变量矩阵。
[0054]
步骤4：构建kde
‑
bn模型、fa
‑
bn模型及kde
‑
fa
‑
bn模型，进行fd故障检测，对冷水机组的故障特征刻画的结果进行评价。
[0055]
基于kde
‑
bn模型进行fd故障检测包括下述过程：
[0056]
41a)确定kde
‑
bn模型结构：
[0057]
如图2所示，用于fd的kde
‑
bn的结构包括：第1层检测层和第2层特征层；第1层中的节点代表冷水机组检测结果，即正常和故障状态；第2层中的节点由n个特征组成；
[0058]
41b)确定kde
‑
bn模型参数：
[0059]
包括每个根节点的先验概率和每个子节点的条件概率；第一层检测层的正常状态和故障状态的先验概率由专家经验或通过对维修服务的历史记录统计其发生频率确定；
[0060]
第二层特征层中的特征节点的条件概率通过对历史数据采用kde方法进行故障特征选择后，带入5维核密度函数估计得到，由下式确定：
[0061][0062]
其中，k()表示高斯核函数；h1,h2,h3,h4,h5分别表示因子变量1、2、3、4、5在特征节点各状态下的数据带宽，x
i1
,x
i2
,x
i3
,x
i4
,x
i5
分别表示一维至五维的故障总体x的样本，x1,
x2,x3,x4,x5分别表示一维至五维的任意点处的数值；
[0063]
41c)基于kde
‑
bn模型进行fd故障检测：
[0064]
进行故障检测的规则如下：
[0065][0066]
其中，s为故障特征，f
i
为冷水机组运行状态，p(f
i
)为类f
i
的先验概率，p(s|f
i
)为给定f
i
时s的条件概率。
[0067]
基于fa
‑
bn模型进行fd故障检测包括下述过程：
[0068]
42a)按照步骤41a)的方法确定fa
‑
bn模型结构；
[0069]
42b)确定fa
‑
bn模型参数：
[0070]
包括每个根节点的先验概率和每个子节点的条件概率；第一层检测层的正常状态和故障状态的先验概率由专家经验或通过对维修服务的历史记录统计其发生频率确定；
[0071]
第二层特征层中的特征节点的条件概率通过对历史数据采用fa方法进行故障特征降维后，带入5维核密度函数估计得到；
[0072]
42c)按照步骤41c)的方法基于fa
‑
bn模型进行fd故障检测。
[0073]
基于kde
‑
fa
‑
bn模型进行fd故障检测包括下述过程：
[0074]
43a)按照步骤41a)的方法确定kde
‑
fa
‑
bn模型结构；
[0075]
43b)确定kde
‑
fa
‑
bn模型参数：
[0076]
包括每个根节点的先验概率和每个子节点的条件概率；第一层检测层的正常状态和故障状态的先验概率由专家经验或通过对维修服务的历史记录统计其发生频率确定。
[0077]
第二层特征层中的特征节点的条件概率通过对历史数据先采用kde方法进行故障特征选择，再采用fa方法进行故障特征降维后，最后带入5维核密度函数估计得到；
[0078]
43c)按照步骤41c)的方法基于kde
‑
fa
‑
bn模型进行fd故障检测。
[0079]
步骤5：构建kde
‑
svdd模型、fa
‑
svdd模型及kde
‑
fa
‑
svdd模型，进行fd故障检测，对冷水机组的故障特征刻画的结果进行评价：
[0080]
基于kde
‑
svdd模型进行fd故障检测包括下述过程：
[0081]
51a)确定kde
‑
svdd模型结构：
[0082]
kde
‑
svdd模型的结构是：将目标数据映射到高维空间，在高维空间构建出一个以a为中心r为半径且体积最小的超球体；
[0083]
51b)确定kde
‑
svdd模型参数：
[0084]
对所有历史数据采用kde方法进行故障特征选择后确定kde
‑
svdd模型的参数；
[0085]
51c)基于kde
‑
svdd模型进行fd故障检测：
[0086]
svdd的原理如图3所示，kde
‑
svdd模型的目标数据为冷水机组正常运行数据，svdd判别方法即故障检测边界内为正常运行区域，边界外为故障区域；比较冷水机组的观测数据z到超球体中心距离d(z)与超球体半径r的大小来判断是否发生故障；若d(z)≤r，观测数据z在超球体内，判定为正常；反之，判定为故障。
[0087]
基于fa
‑
svdd模型进行fd故障检测包括下述过程：
[0088]
52a)按照步骤51a)的方法确定fa
‑
svdd模型结构；
[0089]
52b)确定fa
‑
svdd模型参数：
[0090]
对所有历史数据采用fa方法进行故障特征降维后确定kde
‑
svdd模型的参数；
[0091]
52c)按照步骤51c)的方法基于fa
‑
svdd模型进行fd故障检测。
[0092]
基于kde
‑
fa
‑
svdd模型进行fd故障检测包括下述过程：
[0093]
53a)按照步骤51a)的方法确定kde
‑
fa
‑
svdd模型结构；
[0094]
53b)确定kde
‑
fa
‑
svdd模型参数：
[0095]
对所有历史数据先采用kde方法进行故障特征选择，再采用fa方法进行故障特征降维后确定kde
‑
fa
‑
svdd模型的参数；
[0096]
53c)按照步骤51c)的方法基于kde
‑
fa
‑
svdd模型进行fd故障检测。
[0097]
步骤6：分别将构建的kde
‑
bn模型、fa
‑
bn模型及kde
‑
fa
‑
bn模型与传统bn模型的检测性能进行比较，若正常状态被误判为故障状态的概率越低，则判定该fd方法效果越好。
[0098]
步骤7：分别将构建的kde
‑
svdd模型、fa
‑
svdd模型及kde
‑
fa
‑
svdd模型与传统svdd模型的检测性能进行比较，若正常状态被误判为故障状态的概率越低，则判定该fd方法效果越好。
[0099]
下面通过具体实施例进一步阐述本发明提出的一种基于kde
‑
fa的冷水机组故障特征刻画方法的检测步骤和验证本发明的有益效果：
[0100]
实施例：
[0101]
本实施例使用的故障历史数据来源于ashrae rp
‑
1043故障实验，为一台90冷吨(约316kw)的离心式冷水机组，热交换器(冷凝器、蒸发器和其他换热器)均为满液式壳管式换热器，膨胀阀为热力膨胀阀，制冷剂选用r134a。实验通过改变参数模拟了27个运行工况；通过实验模拟了7类典型软故障，且每类故障划分为4个劣化等级，从低到高依次为sl1，sl2，sl3，sl4，7类典型故障见表1；通过传感器以10秒的间隔进行数据采样，共需要采集29组数据集，其中包括一组正常数据集和28组故障数据集(7类故障，每类故障分4个劣化等级)，每一组数据包含64个特征参数，每一个特征参数包含5191个数据。
[0102]
表1 7类典型软故障
[0103][0104]
步骤1：选用移动平均滤波器法对冷水机组原始故障数据进行平滑处理，然后选用计算几何加权运行平均值和几何加权运行方差的方法对平滑数据进行稳态筛选；
[0105]
步骤2：采用kde方法对稳态筛选后的故障数据进行特征选择：
[0106]
基于尽量选出所有敏感特征和尽量选择温度参数以节约传感器成本两个原则，对核密度对比图中正常数据曲线与故障数据曲线的偏离度进行主观判断，分别选出7类软故障的敏感特征，结果见表2。7类故障所有的敏感特征即为本研究利用kde方法所选出来的特征，特征选择结果见表3。
[0107]
表2各故障的敏感特征
[0108][0109]
表3特征选择结果
[0110][0111]
步骤3：采用fa方法对kde选择出的21个故障特征进行特征降维：
[0112]
3a)在因子描述环节，进行bartlett球形检验和kmo统计量检验，当bartlett球形检验结果中的卡方统计量足够大、显著性值小于0.05，kmo统计量检验结果中的kmo大于0.5三个条件满足其一，则认为特征之间的冗余性较强，此时应进行fa特征降维。
[0113]
fa方法的bartlett球形检验和kmo统计量检验结果见表4，由结果可知应进行fa特征降维。
[0114]
表4 bartlett球形检验和kmo统计量检验结果
[0115]
[0116]
3b)在因子提取环节，采用方差贡献率法和碎石图法确定公共因子的提取个数；因子提取后输出的公共因子能够反映全部数据信息的95％时，选取此时对应的故障特征的公共因子个数；碎石图中以曲线由陡峭变平缓的转折点为判断点，选取陡峭曲线对应的故障特征的公共因子个数。
[0117]
因子提取结果见表5和图4由结果可知应提取5个公共因子。
[0118]
表5公共因子的累积方差贡献率
[0119][0120][0121]
3c)在因子旋转环节，当故障特征原始变量和提取出的公共因子之间的关系无法得到较好解释时，采用正交旋转法对因子载荷矩阵进行旋转，使两者之间的关系更明确；然后对公共因子的物理意义给出合理解释，并对公共因子命名。首先计算因子载荷矩阵，然后采用最大方差法对因子载荷矩阵进行旋转，输出旋转后的因子载荷矩阵、三维因子空间的载荷散点图(若公共因子数目多于3个，只能输出前3个公共因子的因子空间载荷散点图)。最后根据输出结果对目标变量和公共因子之间的关系进行解释。
[0122]
以上输出结果分别见表6、表7和图5。
[0123]
表6因子载荷矩阵
[0124][0125][0126]
表7旋转后的因子载荷矩阵
[0127][0128]
由表7可知，目标变量teo、tei、tre、pre、tea、t_suc、tsh_suc在公共因子1上的载荷较大，所以公共因子1主要解释了蒸发器特征参数和压缩机吸气口处特征参数的信息；目标变量tco、tci、prc、trc在公共因子2上的载荷较大，所以公共因子2主要解释了冷凝器特征参数的信息；目标变量tr_dis、tsh_dis、kw、tcwd、tewd在公共因子3上的载荷较大，所以公共因子3主要解释了压缩机排气口处特征参数和热交换器温差特征参数的信息；目标变量trc_sub、tca在公共因子4上的载荷较大，所以公共因子4主要解释了制冷剂本身温度和过冷度特征参数的信息；目标变量po_feed、to_sump、to_feed在公共因子5上的载荷较大，所以公共因子5主要解释了润滑油特征参数的信息，这也正好说明了为什么公共因子5解释目标变量信息的累积贡献相对其他4个公共因子较少的原因。
[0129]
从以上数据可知，目标变量和公共因子在因子旋转后能够更好的解释之间的关系，公共因子的物理意义也更加明确。
[0130]
3d)在因子得分环节，确定每一样本在公共因子上的具体数值，将其作为新的变量运用于下一步的模型构建工作；选择回归法确定正常状态和各故障状态在因子变量上的具
体数值。
[0131]
通过fa方法，本发明将21个特征参数的维度降到5维，将新生成的5个因子变量用于后续的fd研究工作。本发明从新变量的正常样本和7类故障的所有劣化等级样本中共随机选取1800个样本，并划分为1200个样本的训练集和600个样本的测试集。训练集的样本数据用于确定fd模型的参数，测试集的样本数据用于对模型的验证。
[0132]
步骤4：构建kde
‑
bn模型、fa
‑
bn模型及kde
‑
fa
‑
bn模型，进行fd故障检测，对冷水机组的故障特征刻画的结果进行评价。
[0133]
分别使用1200个训练样本确定传统bn模型、kde
‑
bn模型、fa
‑
bn模型和kde
‑
fa
‑
bn模型的参数，即各故障的先验概率和各类故障发生时各特征参数的条件概率。然后基于四个训练好的模型对600个测试样本进行判定，并对比fd结果。
[0134]
传统bn模型分别对冷水机组8个状态(正常状态和7类故障状态)下的600个测试样本进行fd判定，结果见表8。
[0135]
表8传统bn模型的fd结果矩阵
[0136][0137]
kde
‑
bn模型分别对冷水机组8个状态(正常状态和7类故障状态)下的600个测试样本进行fd判定，结果见表9。
[0138]
表9 kde
‑
bn模型的fd结果矩阵
[0139][0140]
fa
‑
bn模型分别对冷水机组8个状态(正常状态和7类故障状态)下的600个测试样本进行fd判定，结果见表10。
[0141]
表10 fa
‑
bn模型的fd结果矩阵
[0142][0143]
kde
‑
fa
‑
bn模型分别对冷水机组8个状态(正常状态和7类故障状态)下的600个测试样本进行fd判定，结果见表11。
[0144]
表11 kde
‑
fa
‑
bn模型的fd结果矩阵
[0145][0146]
由表8至表11可计算四种模型的far分别为25.3(152/600)、22.7％(136/600)、23.3％(140/600)、20.2％(121/600)；如图6，kde
‑
fa
‑
bn模型相对于其他三类bn模型，检测结果变化不大，且far相对降低。
[0147]
步骤5：构建kde
‑
svdd模型、fa
‑
svdd模型及kde
‑
fa
‑
svdd模型，进行fd故障检测，对冷水机组的故障特征刻画的结果进行评价。
[0148]
分别使用1200个训练样本确定传统svdd模型、kde
‑
svdd模型、fa
‑
svdd模型、kde
‑
fa
‑
svdd模型的参数，即惩罚因子和松弛变量，然后基于四个训练好的模型对600个测试样本进行判定，并对比fd结果。
[0149]
传统svdd模型分别对冷水机组8个状态(正常状态和7类故障状态)下的600个测试样本进行fd判定，结果见图7。
[0150]
kde
‑
svdd模型分别对冷水机组8个状态(正常状态和7类故障状态)下的600个测试样本进行fd判定，结果见图8。
[0151]
fa
‑
svdd模型分别对冷水机组8个状态(正常状态和7类故障状态)下的600个测试样本进行fd判定，结果见图9。
[0152]
kde
‑
fa
‑
svdd模型分别对冷水机组8个状态(正常状态和7类故障状态)下的600个测试样本进行fd判定，结果见图10。
[0153]
由图7至图10可计算四种方法的far分别为23.8％(143/600)、18.7％(112/600)、22％(132/600)、16.2％(97/600)。如图11，kde
‑
fa
‑
svdd相对于三类svdd模型，检测结果变化不大，且far相对降低。
[0154]
步骤6：分别将构建的kde
‑
bn模型、fa
‑
bn模型及kde
‑
fa
‑
bn模型与传统bn模型的检测性能进行比较，构建的3种bn模型的far均低于传统bn模型，且kde
‑
fa
‑
bn模型的far最低，为20.2％。
[0155]
步骤7：分别将构建的kde
‑
svdd模型、fa
‑
svdd模型及kde
‑
fa
‑
svdd模型与传统svdd模型的检测性能进行比较，构建的3种svdd模型的far均低于传统svdd模型，且kde
‑
fa
‑
svdd模型的far最低，为16.2％。
[0156]
结果同时显示：本发明提出的一种基于kde
‑
fa的冷水机组故障特征刻画方法，kde
方法基于敏感最大化的原则，尽量选出原始特征集中的所有敏感特征，相对于其他特征选择方法，所选出的敏感特征更多，保留信息更完整；fa方法去除了特征之间的冗余性，在保证不损失特征信息的基础上，大幅度降低特征维度，简化后续fd方法的结构复杂度；kde
‑
fa方法在实现特征维度降低的同时，保证fd方法的检测性能提高。
[0157]
本发明并不局限于上述实施例，在本发明公开的技术方案的基础上，本领域的技术人员根据所公开的技术内容，不需要创造性的劳动就可以对其中的一些技术特征作出一些替换和变形，这些替换和变形均在本发明的保护范围内。

技术特征：
1.一种基于kde
‑
fa的冷水机组故障特征刻画方法，其特征在于，包括下述步骤：步骤1：选用移动平均滤波器法对冷水机组原始故障数据进行平滑处理，然后选用计算几何加权运行平均值和几何加权运行方差的方法对平滑数据进行稳态筛选；步骤2：采用kde方法对稳态筛选后的故障数据进行故障特征选择；步骤3：采用fa方法对kde方法选择出的冷水机组故障特征进行特征降维，完成故障特征刻画工作；步骤4：构建kde
‑
bn模型、fa
‑
bn模型及kde
‑
fa
‑
bn模型，进行fd故障检测，对冷水机组的故障特征刻画的结果进行评价；步骤5：构建kde
‑
svdd模型、fa
‑
svdd模型及kde
‑
fa
‑
svdd模型，进行fd故障检测，对冷水机组的故障特征刻画的结果进行评价；步骤6：分别将构建的kde
‑
bn模型、fa
‑
bn模型及kde
‑
fa
‑
bn模型与传统bn模型的检测性能进行比较，若正常状态被误判为故障状态的概率越低，则判定该fd方法效果越好；步骤7：分别将构建的kde
‑
svdd模型、fa
‑
svdd模型及kde
‑
fa
‑
svdd模型与传统svdd模型的检测性能进行比较，若正常状态被误判为故障状态的概率越低，则判定该fd方法效果越好。2.根据权利要求1所述的一种基于kde
‑
fa的冷水机组故障特征刻画方法，其特征在于，所述步骤2中，采用kde方法对冷水机组故障的稳态数据进行特征选择：对于某一故障，设x1，x1，......x
n
是取自故障总体x的样本，在任意点x处的总体密度函数f(x)的kde采用下式计算：其中，k()称为某一故障的核函数，h称为某一故障窗宽或带宽，n为故障数据的样本容量；核函数选择gaussian核函数，故障的最优带宽h
o
采用下式计算：其中，s为故障数据的样本标准差。3.根据权利要求1所述的一种基于kde
‑
fa的冷水机组故障特征刻画方法，其特征在于，所述步骤3中，采用fa方法对kde选择出的故障特征进行特征降维时，包括因子描述、因子提取、因子旋转和因子得分环节。4.根据权利要求3所述的一种基于kde
‑
fa的冷水机组故障特征刻画方法，其特征在于，所述步骤3中，采用fa方法对kde选择出的故障特征进行特征降维时，具体通过四个环节实现：3a)在因子描述环节，进行bartlett球形检验和kmo统计量检验；当bartlett球形检验结果中的卡方统计量足够大、显著性值小于0.05、kmo统计量检验结果中的kmo大于0.5三个条件满足其一，则认为故障特征之间的冗余性较强，此时应采用fa方法进行特征降维；3b)在因子提取环节，采用方差贡献率法和碎石图法确定公共因子的提取个数；因子提取后输出的公共因子能够反映全部数据信息的95％时，选取此时对应的故障特征的公共因
子个数；碎石图中以曲线由陡峭变平缓的转折点为判断点，选取陡峭曲线对应的故障特征的公共因子个数；3c)在因子旋转环节，当故障特征原始变量和提取出的公共因子之间的关系不确定时，采用正交旋转法对因子载荷矩阵进行旋转，使故障特征原始变量和公共因子两者之间的关系明确；对公共因子命名；3d)在因子得分环节，确定每一样本在公共因子上的数值，将其作为新的变量运用于下一步模型构建；选择回归法确定正常状态和各故障状态在因子变量的数值因子变量的数值采用下式计算：其中，a为因子载荷矩阵，∑为标准化故障特征原始变量的协方差矩阵，x为故障特征原始变量矩阵。5.根据权利要求1所述的一种基于kde
‑
fa的冷水机组故障特征刻画方法，其特征在于，所述步骤4中，构建kde
‑
bn模型，进行fd故障检测包括下述过程：41a)确定kde
‑
bn模型结构：用于fd的kde
‑
bn的结构包括：第1层检测层和第2层特征层；第1层中的节点代表冷水机组检测结果，即正常和故障状态；第2层中的节点由n个特征组成；41b)确定kde
‑
bn模型参数：包括每个根节点的先验概率和每个子节点的条件概率；第一层检测层的正常状态和故障状态的先验概率由专家经验或通过对维修服务的历史记录统计其发生频率确定；第二层特征层中的特征节点的条件概率通过对历史数据采用kde方法进行故障特征选择后，带入5维核密度函数估计得到，由下式确定：其中，k()表示高斯核函数；h1,h2,h3,h4,h5分别表示因子变量1、2、3、4、5在特征节点各状态下的数据带宽，x
i1
,x
i2
,x
i3
,x
i4
,x
i5
分别表示一维至五维的故障总体x的样本，x1,x2,x3,x4,x5分别表示一维至五维的任意点处的数值；41c)基于kde
‑
bn模型进行fd故障检测：进行故障检测的规则如下：其中，s为故障特征，f
i
为冷水机组运行状态，p(f
i
)为类f
i
的先验概率，p(s|f
i
)为给定f
i
时s的条件概率。6.根据权利要求5所述的一种基于kde
‑
fa的冷水机组故障特征刻画方法，其特征在于，构建fa
‑
bn模型，进行fd故障检测包括下述过程：42a)按照步骤41a)的方法确定fa
‑
bn模型结构；42b)确定fa
‑
bn模型参数：包括每个根节点的先验概率和每个子节点的条件概率；第一层检测层的正常状态和故障状态的先验概率由专家经验或通过对维修服务的历史记录统计其发生频率确定；
第二层特征层中的特征节点的条件概率通过对历史数据采用fa方法进行故障特征降维后，带入5维核密度函数估计得到；42c)按照步骤41c)的方法基于fa
‑
bn模型进行fd故障检测。7.根据权利要求5所述的一种基于kde
‑
fa的冷水机组故障特征刻画方法，其特征在于，构建kde
‑
fa
‑
bn模型，进行fd故障检测包括下述过程：43a)按照步骤41a)的方法确定kde
‑
fa
‑
bn模型结构；43b)确定kde
‑
fa
‑
bn模型参数：包括每个根节点的先验概率和每个子节点的条件概率；第一层检测层的正常状态和故障状态的先验概率由专家经验或通过对维修服务的历史记录统计其发生频率确定；第二层特征层中的特征节点的条件概率通过对历史数据先采用kde方法进行故障特征选择，再采用fa方法进行故障特征降维后，最后带入5维核密度函数估计得到；43c)按照步骤41c)的方法基于kde
‑
fa
‑
bn模型进行fd故障检测。8.根据权利要求1所述的一种基于kde
‑
fa的冷水机组故障特征刻画方法，其特征在于，所述步骤5中，构建kde
‑
svdd模型，进行fd故障检测包括下述过程：51a)确定kde
‑
svdd模型结构：kde
‑
svdd模型的结构是：将目标数据映射到高维空间，在高维空间构建出一个以a为中心r为半径且体积最小的超球体；51b)确定kde
‑
svdd模型参数：对所有历史数据采用kde方法进行故障特征选择后确定kde
‑
svdd模型的参数；51c)基于kde
‑
svdd模型进行fd故障检测：kde
‑
svdd模型的目标数据为冷水机组正常运行数据，svdd判别方法即故障检测边界内为正常运行区域，边界外为故障区域；比较冷水机组的观测数据z到超球体中心距离d(z)与超球体半径r的大小来判断是否发生故障；若d(z)≤r，观测数据z在超球体内，判定为正常；反之，判定为故障。9.根据权利要求8所述的一种基于kde
‑
fa的冷水机组故障特征刻画方法，其特征在于，构建fa
‑
svdd模型，进行fd故障检测包括下述过程：52a)按照步骤51a)的方法确定fa
‑
svdd模型结构；52b)确定fa
‑
svdd模型参数：对所有历史数据采用fa方法进行故障特征降维后确定kde
‑
svdd模型的参数；52c)按照步骤51c)的方法基于fa
‑
svdd模型进行fd故障检测。10.根据权利要求8所述的一种基于kde
‑
fa的冷水机组故障特征刻画方法，其特征在于，构建kde
‑
fa
‑
svdd模型，进行fd故障检测包括下述过程：53a)按照步骤51a)的方法确定kde
‑
fa
‑
svdd模型结构；53b)确定kde
‑
fa
‑
svdd模型参数：对所有历史数据先采用kde方法进行故障特征选择，再采用fa方法进行故障特征降维后确定kde
‑
fa
‑
svdd模型的参数；53c)按照步骤51c)的方法基于kde
‑
fa
‑
svdd模型进行fd故障检测。
技术总结
本发明公开了一种基于KDE


技术研发人员：王智伟 丁书久 李玉娇
受保护的技术使用者：西安建筑科技大学
技术研发日：2021.04.02
技术公布日：2021/6/29