CN104880217B - 一种基于测量值关联度的故障传感器信息重构方法 - Google Patents
一种基于测量值关联度的故障传感器信息重构方法 Download PDFInfo
- Publication number
- CN104880217B CN104880217B CN201510338144.2A CN201510338144A CN104880217B CN 104880217 B CN104880217 B CN 104880217B CN 201510338144 A CN201510338144 A CN 201510338144A CN 104880217 B CN104880217 B CN 104880217B
- Authority
- CN
- China
- Prior art keywords
- response
- measuring point
- variable
- reconstruct
- formula
- 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.)
- Active
Links
Landscapes
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
本发明提供一种基于测量值关联度的故障传感器信息重构方法该方法包括以下步骤:S1、针对故障传感器,计算该传感器正常工作时所在测点响应与其余测点响应的关联度;S2、通过对比关联度大小,确定关联模型建立所需的响应变量;S3、进而运用偏最小二乘法建立重构变量与响应变量的重构模型,运用结构健康监测实测数据,对故障传感器进行故障传感器响应信息重构。本发明提供的方法对于故障传感器信息重构效果很好,明显的降低了重构误差,同时保证了重构值的变化趋势与实际值的变化趋势也保持一致,显著的整体结构分析的可靠性。
Description
技术领域
本发明属于土木工程结构健康监测、信息处理领域,具体涉及一种基于测量值关联度的故障传感器信息重构方法。
背景技术
数据采集是结构健康监测***中非常重要的环节,结构安全性能的分析取决于数据采集的可靠与否。在土木结构健康监测的实际运营过程中,由于设备老化等因素不可避免有局部的传感器发生故障导致数据失真,同时由于一些传感器不易更换,进而造成数据缺失,降低整体数据的可靠性,难以实现对结构的整体分析,造成决策失误,失去了结构监测的意义;同时结构监测会产生大量的数据,要从海量数据中选取有用信息对故障传感器进行重构十分困难,而各个测点的响应具有关联性,有必要针对不同测点的响应的建立有效的关联模型对故障信息进行重构。
国内外学者对故障数据的恢复进行了相关的研究,有的学者为解决桥梁监测的数据失真问题,通过格兰杰因果关系分析各测点监测数据的关系,极限学习机输入变量选取与失真数据因果关系大的测点响应,实现失真数据的恢复,但这种方法的部分个体的预测会出现较大误差,对于数据的变化趋势也未能保持一致,训练结果不稳定,其训练结果有待进一步改进。有的学者提出了基于神经网络的不可靠数据恢复方法,建立数据的RBF神经网络模型,对挠度数据进行恢复,该方法需要大量数据进行建模,同时仍容易陷入局部最优。有的学者提出了基于卡尔曼滤波的响应预测方法,通过结构模型参数与外部激励建立预测模型,通过已知的加速度测量值预测位移值和速度值,并通过框架模型对方法进行验证,但该方法针对于未知参数的结构时,无法有效建立预测模型。有的学者针对无线传感器的数据问题探讨了数据压缩采样恢复的方法,通过随机矩阵将数据线性投影生成新的数据并进行传递,在基站通过数据压缩采样恢复的方法进行数据恢复,但当原始数据稀疏性较差时,加速度数据的恢复精度降低。有的学者提出基于最小二乘支持向量机预测期的数据恢复方法,利用传感器的时程数据建立基于该方法的预测模型,并通过在线学习诊断进行信号恢复。但该方法耗时少,且由于神经网络恢复方法,但主要针对传感器发生故障短时间传感器的数据重构,对于采样率高的传感器和实现长期故障预测效果不明显。
通过分析国内外学者对数据恢复和预测的研究,发现目前的方法对建立结构稳定的关联模型上仍有不足,对于长期的故障重构效果不明显,有必要研究建立基于监测数据关联模型的方法进行故障数据重构。
发明内容
有鉴于此,本发明的目的在于克服现有技术的不足,提供一种基于测量值关联度的故障传感器信息重构方法。
为实现以上目的,本发明采用如下技术方案:一种基于测量值关联度的故障传感器信息重构方法,该方法包括以下步骤:
S1、针对故障传感器,计算该传感器正常工作时所在测点响应与其余测点响应的关联度;
S2、通过对比关联度大小,确定关联模型建立所需的响应变量;
S3、进而运用偏最小二乘法建立重构变量与响应变量的重构模型,运用结构健康监测实测数据,对故障传感器进行故障传感器响应信息重构。
进一步地,所述步骤S1中采用Pearson相关系数衡量不同类型响应之间的关联程度,具体如下:
设X为m个测点在n个时间步的荷载作用下的结构响应矩阵;
xi=[xi1 xi2 … xin]T (2)
xi代表传感器测点i的响应值;
记y为故障传感器正常运行时的响应值;xi与y的相关系数ρi为:
式中σxi,σy为测点i和故障传感器测点的响应方差,Cov(xi,y)为其协方差,xi(k)和y(k)分别为测点i和故障传感器测点在第k时间步时的响应,μxi和μy为测点i和故障传感器测点在n个时间步的荷载作用下的平均值。
进一步地,所述的步骤S2中确定关联模型建立所需的响应变量的方法具体如下:
故障测点与其余测点响应的关联度矩阵P为
P=[ρ1 ρ2 … ρm] (7)
将P中相关系数ρi的绝对值从大到小排序,记为DP
DP={(ρ)1 …(ρ)i … (ρ)m} (8)
(ρ)i代表排序为第i位的相关系数,对应的测点响应为(x)i;则按照相关系数大小排序后的结构响应矩阵DX为
DX={(x)1 …(x)i … (x)m} (9)
选择p个测点监测信息用作响应变量XP,将故障传感器重构信息看作重构变量YP,即
XP={(x)1 (x)2 … (x)p} (10)
YP={y} (11)
进一步地,所述的步骤S3中建立重构变量与响应变量的重构模型的具体方法如下:
将响应变量和重构变量做标准化处理,
响应变量XP和重构变量YP经标准化处理后为
F0={y*} (15)
运用主成分分析方法,获取E0和F0的主成分集合T和U为
T={t1 t2 ... th ... tA} (16)
U={u1 u2 ... uh ... uA} (17)
运用偏最小二乘法和第一主成分,E0,F0的回归方程为
E0=t1p1 T+E1 (18)
F0=t1r1 T+F1 (19)
回归系数向量是:
取第h个主成分,Eh-1的回归方程为
Eh-1=thph T+Eh (22)
Fh-1=thrh T+Fh (23)
其中
E0的第h个主成分th和E0的第h个轴wh的关系为
th=Eh-1wh (26)
其中,wh是矩阵Eh-1 TFh-1Fh-1 TEh-1最大特征值θh 2的单位特征向量。
分析式(18),(19),(22),(43)可得
E0=t1p1 T+t2p2 T+…+thph T+Eh (27)
F0=t1r1 T+t2r2 T+…+thrh T+Fh (28)
式(26)代入式(22)得
Eh=Eh-1(I-whph T) (29)
式(26)和(29)代入式(27)得
式(30)代入式(26)得
式(31)代入式(28)得
F0=E0w1 *r1 T+…+E0wh *rh T+Fh (33)
取主成分满***叉有效性原则时的个数h,此时Fh趋近于0,可忽略不计。
则
F0=E0(w1 *r1 T+…+wh *rh T) (34)
式(14),(15)代入式(34)得
设
则y*=α1 *(x)1 *+…+αp *(x)p * (37)
式(12),(13)代入式(37)得重构变量与响应变量的重构方程为
y=α1(x)1+…+αp(x)p+α0 (38)。
本发明中所述的传感器为应力应变传感器、加速度传感器、位移传感器或风压传感器、风速仪等监测土木工程结构荷载与响应的传感器中的一种或两种以上的组合。
本发明采用以上技术方案,针对出现故障的传感器,计算该传感器正常工作时所在测点响应与其余测点响应的关联度,即重构变量与响应变量的关联度。通过对比关联度大小,确定关联模型建立所需的响应变量;进而运用偏最小二乘法建立重构变量与响应变量的重构模型,运用结构健康监测实测数据,对故障传感器进行故障传感器响应信息重构。本发明提供的方法对于故障传感器信息重构效果很好,明显的降低了重构误差,同时保证了重构值的变化趋势与实际值的变化趋势也保持一致,显著的整体结构分析的可靠性。
附图说明
图1是本发明基于测量值关联度的故障传感器信息重构方法流程图;
具体实施方式
下面将结合附图,对本发明实施例中的技术方案进行清楚、完整的描述。
如图1所示,本发明提供了一种基于测量值关联度的故障传感器信息重构方法,该方法包括以下步骤:
S1、针对故障传感器,计算该传感器正常工作时所在测点响应与其余测点响应的关联度;
S2、通过对比关联度大小,确定关联模型建立所需的响应变量;
S3、进而运用偏最小二乘法建立重构变量与响应变量的重构模型,运用结构健康监测实测数据,对故障传感器进行故障传感器响应信息重构。
结构响应间存在着关联性,Pearson相关系数是反映结构响应线性相关程度的指标,本实施例中的步骤S1中采用Pearson相关系数衡量不同类型响应之间的关联程度,具体如下:
设X为m个测点在n个时间步的荷载作用下的结构响应矩阵;
xi=[xi1 xi2 … xin]T (2)
xi代表传感器测点i的响应值;
记y为故障传感器正常运行时的响应值;xi与y的相关系数ρi为:
式中σxi,σy为测点i和故障传感器测点的响应方差,Cov(xi,y)为其协方差,xi(k)和y(k)分别为测点i和故障传感器测点在第k时间步时的响应,μxi和μy为测点i和故障传感器测点在n个时间步的荷载作用下的平均值。
作为一种优选的实施方式,所述的步骤S2中确定关联模型建立所需的响应变量的方法具体如下:
故障测点与其余测点响应的关联度矩阵P为
P=[ρ1 ρ2 … ρm] (7)
将P中相关系数ρi的绝对值从大到小排序,记为DP
DP={(ρ)1 …(ρ)i … (ρ)m} (8)
(ρ)i代表排序为第i位的相关系数,对应的测点响应为(x)i;则按照相关系数大小排序后的结构响应矩阵DX为
DX={(x)1 …(x)i …(x)m} (9)
选择p个测点监测信息用作响应变量XP,将故障传感器重构信息看作重构变量YP,即
XP={(x)1 (x)2 … (x)p} (10)
YP={y} (11)
现对步骤S3作进一步详述,关联模型是建立重构变量和响应变量的关系模型。偏最小二乘法是一种多元统计数据分析方法,是重构变量对多响应变量的建模方法[9],可以用来处理用普通多元回归无法解决的问题。多元回归中,如果选用的响应变量存在多重相关性,采用普通的最小二乘法会降低参数估计的准确性,增大模型误差;偏最小二乘法利用多元线性回归分析、主成分分析以及典型相关分析的原理对数据进行分解和筛选的方式,提取对重构变量解释性最强的综合变量,获取重构变量对选取的响应变量的回归模型。
所以本实施例中的步骤S3中建立重构变量与响应变量的重构模型的具体方法如下:
将响应变量和重构变量做标准化处理,
响应变量XP和重构变量YP经标准化处理后为
F0={y*} (15)
运用主成分分析方法,获取E0和F0的主成分集合T和U为
T={t1 t2 ...th ... tA} (16)
U={u1 u2 ... uh ... uA} (17)
运用偏最小二乘法和第一主成分,E0,F0的回归方程为
E0=t1p1 T+E1 (18)
F0=t1r1 T+F1 (19)
回归系数向量是:
取第h个主成分,Eh-1的回归方程为
Eh-1=thph T+Eh (22)
Fh-1=thrh T+Fh (23)
其中
E0的第h个主成分th和E0的第h个轴wh的关系为
th=Eh-1wh (26)
其中,wh是矩阵Eh-1 TFh-1Fh-1 TEh-1最大特征值θh 2的单位特征向量。
分析式(18),(19),(22),(43)可得
E0=t1p1 T+t2p2 T+…+thph T+Eh (27)
F0=t1r1 T+t2r2 T+…+thrh T+Fh (28)
式(26)代入式(22)得
Eh=Eh-1(I-whph T) (29)
式(26)和(29)代入式(27)得
式(30)代入式(26)得
式(31)代入式(28)得
F0=E0w1 *r1 T+…+E0wh *rh T+Fh (33)
取主成分满***叉有效性原则时的个数h,此时Fh趋近于0,可忽略不计。则
F0=E0(w1 *r1 T+…+wh *rh T) (34)
式(14),(15)代入式(34)得
设
则y*=α1 *(x)1 *+…+αp *(x)p * (37)
式(12),(13)代入式(37)得重构变量与响应变量的重构方程为
y=α1(x)1+…+αp(x)p+α0 (38)。
需要说明的是:本发明中所述的传感器为应力应变传感器、加速度传感器、位移传感器或风压传感器、风速仪等监测土木工程结构荷载与响应的传感器中的一种或两种以上的组合。
本发明不局限于上述最佳实施方式,任何人在本发明的启示下都可得出其他各种形式的产品,但不论在其形状或结构上作任何变化,凡是具有与本申请相同或相近似的技术方案,均落在本发明的保护范围之内。
Claims (3)
1.一种基于测量值关联度的故障传感器信息重构方法,其特征在于:该方法包括以下步骤:
S1、针对故障传感器,计算该传感器正常工作时所在测点响应与其余测点响应的关联度;
S2、通过对比关联度大小,确定关联模型建立所需的响应变量;
S3、进而运用偏最小二乘法建立重构变量与响应变量的重构模型,运用结构健康监测实测数据,对故障传感器进行故障传感器响应信息重构;
所述步骤S1中采用Pearson相关系数衡量不同类型响应之间的关联程度,具体如下:
设X为m个测点在n个时间步的荷载作用下的结构响应矩阵;
xi=[xi1 xi2 … xin]T (2)
xi代表传感器测点i的响应值;
记y为故障传感器正常运行时的响应值;xi与y的相关系数ρi为:
式中σxi,σy为测点i和故障传感器测点的响应方差,Cov(xi,y)为其协方差,xi(k)和y(k)分别为测点i和故障传感器测点在第k时间步时的响应,μxi和μy为测点i和故障传感器测点在n个时间步的荷载作用下的平均值;
所述的步骤S2中确定关联模型建立所需的响应变量的方法具体如下:
故障测点与其余测点响应的关联度矩阵P为
P=[ρ1 ρ2 … ρm] (7)
将P中相关系数ρi的绝对值从大到小排序,记为DP
DP={(ρ)1 … (ρ)i … (ρ)m} (8)
(ρ)i代表排序为第i位的相关系数,对应的测点响应为(x)i;则按照相关系数大小排序后的结构响应矩阵DX为
DX={(x)1 … (x)i … (x)m} (9)
选择p个测点监测信息用作响应变量XP,将故障传感器重构信息看作重构变量YP,即
XP={(x)1 (x)2 … (x)p} (10)
YP={y} (11)
所述的步骤S3中建立重构变量与响应变量的重构模型的具体方法如下:
将响应变量和重构变量做标准化处理,
响应变量XP和重构变量YP经标准化处理后为
F0={y*} (15)
运用主成分分析方法,获取E0和F0的主成分集合T和U为
T={t1 t2 ... th ... tA} (16)
U={u1 u2 ... uh ... uA} (17)
运用偏最小二乘法和第一主成分,E0,F0的回归方程为
E0=t1p1 T+E1 (18)
F0=t1r1 T+F1 (19)
回归系数向量是
取第h个主成分,Eh-1的回归方程为
Eh-1=thph T+Eh (22)
Fh-1=thrh T+Fh (23)
其中
E0的第h个主成分th和E0的第h个轴wh的关系为
th=Eh-1wh (26)
其中,wh是矩阵Eh-1 TFh-1Fh-1 TEh-1最大特征值θh 2的单位特征向量,分析式(18),(19),(22),(23)可得
E0=t1p1 T+t2p2 T+…+thph T+Eh (27)
F0=t1r1 T+t2r2 T+…+thrh T+Fh (28)
式(26)代入式(22)得
Eh=Eh-1(I-whph T) (29)
式(26)和(29)代入式(27)得
式(30)代入式(26)得
式(31)代入式(28)得
F0=E0w1 *r1 T+…+E0wh *rh T+Fh
(33)
取主成分满***叉有效性原则时的个数h,此时Fh趋近于0,可忽略不计,则
F0=E0(w1 *r1 T+…+wh *rh T) (34)
式(14),(15)代入式(34)得
设
则y*=α1 *(x)1 *+…+αp *(x)p * (37)
式(12),(13)代入式(37)得重构变量与响应变量的重构方程为
y=α1(x)1+…+αp(x)p+α0 (38)。
2.根据权利要求1所述的一种基于测量值关联度的故障传感器信息重构方法,其特征在于:所述的传感器为应力应变传感器、加速度传感器、位移传感器、风压传感器或风速仪中的一种或两种以上的组合。
3.根据权利要求1所述的一种基于测量值关联度的故障传感器信息重构方法,其特征在于:所述的传感器为监测土木工程结构荷载与响应的传感器。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510338144.2A CN104880217B (zh) | 2015-06-17 | 2015-06-17 | 一种基于测量值关联度的故障传感器信息重构方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510338144.2A CN104880217B (zh) | 2015-06-17 | 2015-06-17 | 一种基于测量值关联度的故障传感器信息重构方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104880217A CN104880217A (zh) | 2015-09-02 |
CN104880217B true CN104880217B (zh) | 2017-08-15 |
Family
ID=53947803
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510338144.2A Active CN104880217B (zh) | 2015-06-17 | 2015-06-17 | 一种基于测量值关联度的故障传感器信息重构方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104880217B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106372426A (zh) * | 2016-08-31 | 2017-02-01 | 郑州航空工业管理学院 | 基于主成分分析与神经网络的多响应参数优化方法 |
CN108416309B (zh) * | 2018-03-14 | 2022-02-18 | 揭阳职业技术学院 | 一种智能传感器多故障传感信号重构方法 |
CN109738220B (zh) * | 2019-01-07 | 2021-01-08 | 哈尔滨工业大学(深圳) | 一种基于多荷载工况结构响应关联的传感器优化布置方法 |
CN112179655B (zh) * | 2020-08-17 | 2021-07-09 | 中国农业大学 | 一种基于阈值分级的汽轮发电机故障预警方法 |
CN111982256B (zh) * | 2020-08-18 | 2021-11-16 | 上海陆杰电子科技有限公司 | 一种基于多通道cnn模型的汽车衡压力传感器故障预警方法 |
CN116484306B (zh) * | 2023-06-20 | 2023-09-26 | 蘑菇物联技术(深圳)有限公司 | 异常传感器的定位方法、装置、计算机设备和存储介质 |
CN116698323B (zh) * | 2023-08-07 | 2023-10-13 | 四川华腾公路试验检测有限责任公司 | 一种基于pca和扩展卡尔曼滤波的桥梁健康监测方法及*** |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101382439B (zh) * | 2008-10-27 | 2010-12-29 | 哈尔滨工业大学 | 多参数自确认传感器的状态自确认方法 |
CN101477375B (zh) * | 2009-01-05 | 2012-01-04 | 东南大学 | 基于矩阵奇异值关联规则挖掘的传感器数据校验方法 |
CN103198638B (zh) * | 2013-04-27 | 2015-11-25 | 上海交通大学 | 基于无线传感器执行器网络的铁路桥梁结构健康监测*** |
CN104462770B (zh) * | 2014-11-10 | 2017-04-12 | 福州大学 | 考虑传感器性能退化的信号重构方法 |
CN105318905B (zh) * | 2014-12-30 | 2018-05-25 | 国家电网公司 | 一种火电厂热力***传感器故障诊断方法 |
-
2015
- 2015-06-17 CN CN201510338144.2A patent/CN104880217B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104880217A (zh) | 2015-09-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104880217B (zh) | 一种基于测量值关联度的故障传感器信息重构方法 | |
CN112033463B (zh) | 一种核动力设备状态评估与预测一体化方法和*** | |
CN103776480B (zh) | 基于多次移动平均的微小故障检测方法和装置 | |
CN109489977B (zh) | 基于KNN-AdaBoost的轴承故障诊断方法 | |
CN109034191B (zh) | 基于elm的单维遥测数据异常判读方法 | |
CN111737909B (zh) | 基于时空图卷积网络的结构健康监测数据异常识别方法 | |
CN108399277B (zh) | 一种基于温度与应变关联性的桥梁损伤识别方法 | |
CN108873859B (zh) | 基于改进关联规则的桥式抓斗卸船机故障预测模型方法 | |
CN105629958B (zh) | 一种基于子时段mpca‑svm的间歇过程故障诊断方法 | |
CN111946559B (zh) | 一种风机基础和塔架结构检测方法 | |
CN105716814A (zh) | 一种评估桁架结构损伤的实时监测***及其方法 | |
CN110348150A (zh) | 一种基于相关概率模型的故障检测方法 | |
CN112132430B (zh) | 一种配电主设备分布式状态传感器可靠性评估方法及*** | |
CN112785091A (zh) | 一种对油田电潜泵进行故障预测与健康管理的方法 | |
CN112395684A (zh) | 一种高速列车走行部***智能故障诊断方法 | |
Son et al. | Deep learning-based anomaly detection to classify inaccurate data and damaged condition of a cable-stayed bridge | |
CN109543743A (zh) | 一种基于重建预测残差的制冷机组多传感器故障诊断方法 | |
Moeinifard et al. | Lost vibration test data recovery using convolutional neural network: a case study | |
CN114091600A (zh) | 一种数据驱动的卫星关联故障传播路径辨识方法及*** | |
CN107918381A (zh) | 一种基于组合核函数的类均值核主元故障诊断方法 | |
CN113988210A (zh) | 结构监测传感网失真数据修复方法、装置及存储介质 | |
CN207992717U (zh) | 一种水电站闸门在线状态监测*** | |
CN102938068B (zh) | 桥梁结构多体系损伤识别方法 | |
CN114266013A (zh) | 基于深度学习虚拟感知网络的传动***振动解耦方法 | |
CN107727392A (zh) | 基于信号检测和roc分析的状态指数评估和优化方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |