CN113723010B - 一种基于lstm温度-位移相关模型的桥梁损伤预警方法 - Google Patents
一种基于lstm温度-位移相关模型的桥梁损伤预警方法 Download PDFInfo
- Publication number
- CN113723010B CN113723010B CN202111052885.6A CN202111052885A CN113723010B CN 113723010 B CN113723010 B CN 113723010B CN 202111052885 A CN202111052885 A CN 202111052885A CN 113723010 B CN113723010 B CN 113723010B
- Authority
- CN
- China
- Prior art keywords
- data
- temperature
- strain
- model
- lstm
- 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
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000006073 displacement reaction Methods 0.000 title claims abstract description 15
- 238000012549 training Methods 0.000 claims abstract description 61
- 238000012544 monitoring process Methods 0.000 claims abstract description 25
- 238000012360 testing method Methods 0.000 claims abstract description 19
- 238000003062 neural network model Methods 0.000 claims abstract description 9
- 238000007781 pre-processing Methods 0.000 claims abstract description 5
- 230000011218 segmentation Effects 0.000 claims abstract description 4
- 230000002159 abnormal effect Effects 0.000 claims description 11
- 238000010219 correlation analysis Methods 0.000 claims description 9
- ORILYTVJVMAKLC-UHFFFAOYSA-N Adamantane Natural products C1C(C2)CC3CC1CC2C3 ORILYTVJVMAKLC-UHFFFAOYSA-N 0.000 claims description 8
- 238000012216 screening Methods 0.000 claims description 7
- 230000009466 transformation Effects 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 5
- 230000004913 activation Effects 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 4
- 238000005457 optimization Methods 0.000 claims description 4
- 238000012163 sequencing technique Methods 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 2
- 230000006870 function Effects 0.000 description 18
- 230000000875 corresponding effect Effects 0.000 description 12
- 230000036541 health Effects 0.000 description 7
- 230000007774 longterm Effects 0.000 description 5
- 238000013528 artificial neural network Methods 0.000 description 3
- 238000010606 normalization Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000004927 fusion Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 241001123248 Arma Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008033 biological extinction Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000013135 deep learning Methods 0.000 description 1
- 238000013136 deep learning model Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
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
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/25—Fusion techniques
-
- 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/04—Architecture, e.g. interconnection topology
- G06N3/044—Recurrent networks, e.g. Hopfield networks
-
- 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
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- 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
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Software Systems (AREA)
- Computing Systems (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Computational Linguistics (AREA)
- Biophysics (AREA)
- Mathematical Physics (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Medical Informatics (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明涉及一种基于LSTM温度‑位移相关模型的桥梁损伤预警方法。包括:对桥梁监测数据进行预处理,提取出相关性最强的温度传感器数据和应变传感器数据,并遍历进行数据分割;融合温度传感器数据和应变传感器数据,提取两个传感器数据间的低阶信息和高阶信息,并得到融合后的温度‑应变传感器数据,而后划分为训练集和测试集,并打乱每组训练集中数据的顺序;搭建具有两个LSTM层和一个Dense层的神经网络模型,训练得到基于LSTM温度‑应变相关模型;根据当前时刻桥梁数据通过基于LSTM温度‑应变相关模型进行预测,当桥梁下一时刻的真实数据传回时,计算真实值与预测值之间的偏差,判断偏差是否大于设定的阈值,从而实现桥梁损伤预警。
Description
技术领域
本发明涉及桥梁监测技术领域,具体涉及一种基于LSTM温度-位移相关模型的桥梁损伤预警方法。
背景技术
随着铁路和高速公路的快速发展,桥梁的数量也在不断增多。然而桥梁结构易受外在和内在因素的影响发生变化,从而导致安全事故的发生,因此对桥梁健康状况的监测也就尤为重要。虽然当今国内外很多大跨度桥梁基本都安装了桥梁健康监测***,监测内容不仅是桥梁自身的状况与行为,更注重对环境条件的监测与记录分析,但是这些桥梁健康监测***极大部分仅对监测数据进行采集和保存,极少具备数据分析的功能。对桥梁监测数据进行数据分析的主要目的是为了实现桥梁损伤预警,提供桥梁健康评估功能,同时减少传统桥梁检测所需要投入的人力与物力。目前常用的桥梁预警方法包括神经网络预警方法、时间序列预警方法、灰色***预警方法等。
目前桥梁健康监测***在桥梁上的应用日益增多,损伤预警健康监测理论也取得相应的进展,如:趋势预测应用于山西风陵渡大桥上、时间序列ARMA模型应用于重庆高家花园嘉陵江大桥长期检测扰度数据预测、神经网络应用于重庆马桑溪大桥跨中扰度数据预测等,通过预测数据与真实数据的比较能够对桥梁损伤起到预警作用。
应变能够表征大跨桥梁由于局部损伤或构件破坏引起的受力性能变化,是评估大跨桥梁服役安全的有效指标。然而,不仅损伤会导致大跨桥梁应变变化,温度也会引起显著的应变变化。错误估计温度引起的应变不仅会降低桥梁安全评估结果的可靠性,甚至可能得出错误的结论。因此,有必要深入研究温度引起的桥梁应变响应。
然而应用于桥梁损伤预警的时间序列模型和深度学习中的序列模型大都仅考虑传感器数据的自相关性,即用该传感器的历史数据建立模型来预测下一时刻的数据,忽略了其它因素可能对该传感器数据产生的影响。本专利是一种基于LSTM温度-位移相关模型的桥梁损伤预警方法,不仅考虑单个应变传感器数据的历史信息,同时将与之相关性较强的温度传感器数据也纳入模型训练进行预测,从而实现桥梁损伤预警。
发明内容
本发明的目的在于提供一种基于LSTM温度-位移相关模型的桥梁损伤预警方法,该方法使用温度-应变多变量的LSTM模型预测应变,相比现有的使用单变量LSTM模型进行预测,本发明模型预测值与真实值之间的平均绝对误差更小,模型的预测更为精准。
为实现上述目的,本发明的技术方案是:一种基于LSTM温度-位移相关模型的桥梁损伤预警方法,包括如下步骤:
S1、对桥梁监测数据进行预处理;
S2、对预处理后的数据进行相关性分析,提取出相关性最强的温度传感器数据和应变传感器数据;
S3、遍历基于步骤S2提取的温度传感器数据和应变传感器数据,进行数据分割;
S4、融合温度传感器数据和应变传感器数据,提取两个传感器数据间的低阶信息和高阶信息,并得到融合后的温度-应变传感器数据;
S5、将融合后的温度-应变传感器数据划分为训练集和测试集,并打乱每组训练集中数据的顺序;
S6、搭建具有两个LSTM层和一个Dense层的神经网络模型,输入训练集数据进行训练,得到基于LSTM温度-应变相关模型;
S7、根据当前时刻桥梁数据通过基于LSTM温度-应变相关模型进行预测,当桥梁下一时刻的真实数据传回时,计算真实值与预测值之间的偏差,判断偏差是否大于设定的阈值,从而实现桥梁损伤预警。
在本发明一实施例中,所述步骤S1具体实现方式为:
S11、准备桥梁监测数据原始数据集,其数据结构为【时间、桥梁名称、传感器类型、传感器编号、对应时间的传感器数据平均值、传感器位置】;将桥梁监测数据原始数据集按时间属性进行排序;
S12、将桥梁监测数据原始数据集整理成时间为唯一索引,列名为传感器编号的数据表,其数据为对应时间的传感器数据单位时间的平均值;将数据整理成数据表后,对单位时间间隔的时间缺失进行填补,将该时间所对应的缺失数据使用空值来进行替代;
S13、对所有接收到的传感器数据进行筛选,筛选出不在量程、精度误差范围内的所有异常值并剔除,剔除完成后在对应位置使用空值进行替代;
S14、对于数据中存在的空值,将空值用对应列的均值填补;
S15、对温度传感器和应变传感器的数据做归一化的处理,将数据映射到0和1之间:
其中x为待归一化的数据,x_min和x_max分别表示待归一化数据中的最小值和最大值。
在本发明一实施例中,所述步骤S2具体实现方式为:
S21、对数据表中各个温度传感器和应变传感器进行相关性分析,计算两两传感器之间的皮尔逊相关系数;根据相关系数的强弱,找出相关性最强的温度传感器和应变传感器;
S22、筛选出相关性最强的温度传感器和应变传感器后,将两个传感器的数据按时间进行内连接,其数据结构为【时间、温度传感器、应变传感器】,其对应数据为传感器单位时间的平均值。
在本发明一实施例中,所述步骤S3具体实现方式为:
将步骤S2提取的相关性最强的温度传感器数据和应变传感器数据中温度传感器的数据记为:
其中表示温度传感器数据中第i个数据,i∈{1,2…n};
应变传感器的数据记为:
其中表示应变传感器数据中第i个数据,i∈{1,2…n};
遍历温度传感器数据和应变传感器数据,设置滑动窗口的长度为w,提取连续w个单位时间的温度数据和应变数据作为输入特征,第w+1到w+f单位时间的应变数据作为标签对数据进行分割,窗口滑动步长为1,其中f为预测长度,表示向后预测f个单位时间的应变数据。
在本发明一实施例中,所述步骤S4具体实现方式为:
将连续w个小时的温度数据和应变数据进行融合,提取温度与应变数据间的低阶信息和高阶信息,其中:
低阶信息为温度与应变数据的线性变换,表示为:
xl=Wl[xt,xs]+bl
其中Wl为可训练参数矩阵,xt和xs分别为温度数据与应变数据,bl为偏置量;
高阶信息为温度与应变数据的线性变换再加上交叉项信息,表示为:
xm=Wm[xt,xs]+Wtsxtxs+bm
其中Wm为可训练参数矩阵,xt和xs分别为温度数据与应变数据,bm为偏置量;
低阶信息和高阶信息经过sigmoid函数激活后便得到最终的特征:
x=sigmoid(xl+xm)。
xl和xm为计算出的低阶信息和高阶信息。
在本发明一实施例中,所述步骤S5具体实现方式为:
将融合后的温度-应变传感器数据划分为训练集和测试集,其中80%的数据用于训练,20%的数据用于测试,并对训练集数据进行shuffle,打乱每组训练集中数据的顺序。
在本发明一实施例中,所述步骤S6具体实现方式为:
搭建两个LSTM层,和一个全连接Dense层的神经网络模型,设置激活函数为tanh,同时,为避免模型过拟合,每层LSTM设置Dropout率为0.2;在训练过程中以平均绝对误差MAE作为损失函数,以平均绝对误差最小化为优化目标,采用Adam优化器,设置每批训练的样本数BatchSize为128,训练集最大迭代次数Epoch为500,模型训练过程如下:
S61、数据输入
将步骤S5的数据作为神经网络模型的输入;
S62、LSTM层和Dense层
输入的数据经过两个LSTM层和一个Dense层后输出预测应变数据;
S63、计算损失函数
采用MAE作为模型的损失函数,计算Dense层输出的预测应变与真实应变间的损失;
S64、梯度下降
使用Adam优化器梯度下降更新参数,使得模型的损失函数最小化,直到模型收敛或者模型到达最大之迭代次数,得到训练后的基于LSTM温度-应变相关模型。
在本发明一实施例中,所述MAE公式如下:
其中X为样本空间,h为训练的模型,m为样本总数,h(xi)表示第i个样本经过模型h输出的预测值,yi表示第i个样本的真实值。
在本发明一实施例中,步骤S63中,为更加直观的体现模型在训练集和测试集上的误差,绘制损失曲线,其中横轴为训练集迭代次数Epoch,纵轴为模型的损失函数值Loss。
相较于现有技术,本发明具有以下有益效果:
1、本发明提出的基于温度-应变相关模型的LSTM桥梁损伤预警方法,使用温度-应变多特征融合的LSTM模型预测应变,相比现有的使用单变量LSTM模型进行预测,我们的模型预测值与真实值之间的平均绝对误差更小,模型的预测更为精准;
2、本发明通过设置超参数w、f对数据分隔进行控制,设置滑动窗口的长度为w,提取训练集中连续f个单位时间的温度和应变数据作为输入特征,第w到w+f个单位时间的应变数据作为标签对训练集进行分割,采用同样的方式处理测试集。其中f为预测长度,表示向后预测f个单位时间,提升模型的泛化能力。
3、本发明采用没有长期依赖问题的LSTM算法,是因为桥梁的数据都为时间序列数据,时间序列数据是在不同时间点上统计同一指标,并按照时间先后排列成的一个集合。时间序列的主要作用是了解一个指标的长期趋势和预测未来,由于传统的RNN存在长期依赖问题,长期依赖的根本问题在于,多阶段的反向传播会导致梯度消失和梯度***。因此相比传统的RNN,使用LSTM能够有更好的表现。
4、本发明使用Adam优化器可以使计算更加高效,使得模型更有效更稳定,适合应用于桥梁大规模的数据和参数的场景;使用SHUFFLE方法避免了数据投入的顺序对网络训练造成影响,增加随机性,提高模型的泛化能力;激活函数使用tanh可以提高模型的训练效率;使用Dropout防止模型过拟合,提升模型泛化能力。
附图说明
图1为本发明方法流程图。
图2为LSTM原理图。
图3为本发明一实例损失曲线图。
图4为本发明一实例预测值与真实值对比图。
具体实施方式
下面结合附图,对本发明的技术方案进行具体说明。
本发明一种基于LSTM温度-位移相关模型的桥梁损伤预警方法,包括如下步骤:
S1、对桥梁监测数据进行预处理;
S2、对预处理后的数据进行相关性分析,提取出相关性最强的温度传感器数据和应变传感器数据;
S3、遍历基于步骤S2提取的温度传感器数据和应变传感器数据,进行数据分割;
S4、融合温度传感器数据和应变传感器数据,提取两个传感器数据间的低阶信息和高阶信息,并得到融合后的温度-应变传感器数据;
S5、将融合后的温度-应变传感器数据划分为训练集和测试集,并打乱每组训练集中数据的顺序;
S6、搭建具有两个LSTM层和一个Dense层的神经网络模型,输入训练集数据进行训练,得到基于LSTM温度-应变相关模型;
S7、根据当前时刻桥梁数据通过基于LSTM温度-应变相关模型进行预测,当桥梁下一时刻的真实数据传回时,计算真实值与预测值之间的偏差,判断偏差是否大于设定的阈值,从而实现桥梁损伤预警。
以下为本发明具体实施实例。
如图1所示,本发明一种基于LSTM温度-位移相关模型的桥梁损伤预警方法,利用桥梁传感器传回的历史数据,首先对其进行预处理,主要包括去除异常值、填补缺失值以及归一化处理。然后进行相关性分析,提取出相关性较强的温度和应变传感器数据。进一步遍历提取的温度和应变传感器数据,设置滑动窗口的长度为w,提取连续w个单位时间的温度和应变数据作为输入特征,第w+1到w+f单位时间的应变数据作为标签对数据进行分割。同时融合温度和应变数据,提取两个传感器数据间的低阶信息和高阶信息。划分训练集和测试集,其中80%的数据用于训练模型,20%的数据用于测试。最终搭建具有两个LSTM层和一个Dense层的神经网络结构,以前w个单位时间的温度和应变数据作为输入,以w+1到w+f个单位时间的应变数据作为标签纳入模型训练,在训练过程中以平均绝对误差(MeanAbsolute Error,MAE)作为损失函数,以平均绝对误差最小化为优化目标,采用Adam优化器,设置每批训练的样本数(BatchSize)为128,训练集最大迭代次数(Epoch)500,最终输出预测应变数据。当下一时刻的真实数据传回时,计算真实值与预测值之间的偏差,判断偏差是否大于设定的阈值,从而实现桥梁损伤预警。具体流程如下:
步骤1、桥梁监测数据
步骤1.1、数据输入
准备桥梁监测数据原始数据集,其数据结构为:【时间、桥梁名称、传感器类型、传感器编号、对应时间的传感器数据平均值(时间间隔可为小时、分钟或秒,根据计算能力调整)、传感器位置】。将桥梁监测数据原始数据集按时间属性进行排序。
步骤1.2、数据整理
将桥梁监测数据原始数据集整理成时间为唯一索引,列名为传感器编号的数据表,其数据为对应时间的传感器数据单位时间的平均值。将数据整理成数据表后,对单位时间间隔的时间缺失进行填补,将该时间所对应的缺失数据使用空值来进行替代。以这样的格式进行输入,方便后续的步骤2进行处理。
步骤2、数据预处理
步骤2.1、异常值处理
步骤1.2完成之后,对所有接收到的传感器数据进行筛选,筛选出不在量程、精度误差范围内的所有异常值并剔除,剔除完成后在对应位置使用空值进行替代。
步骤2.2、缺失值填补
步骤2.1完成之后,对于数据中存在的空值,将空值用对应列的均值填补。
步骤2.3、数据归一化
为了提升模型的精度,同时加快模型的收敛速度,剔除量纲对模型的影响。对温度传感器和应变传感器的数据做归一化(Min-Max Normalization)的处理,将数据映射到0和1之间。
步骤3、相关性分析
步骤3.1
对数据表中各个温度和应变传感器进行相关性分析,计算两两传感器之间的皮尔逊相关系数。根据相关系数的强弱,找出相关性最强的温度和应变传感器。
步骤3.2
筛选出最相关的温度和应变传感器,将两个传感器的数据按时间进行内连接,其数据结构为【时间、温度传感器、应变传感器】,其对应数据为传感器单位时间的平均值。
步骤4、数据分隔
由步骤3的相关性分析已经得出了相关性最强的温度和应变传感器数据,其中温度传感器的数据记为:
应变传感器的数据记为:
遍历提取的温度和应变传感器数据,设置滑动窗口的长度为w,提取连续w个单位时间的温度和应变数据作为输入特征,第w+1到w+f单位时间的应变数据作为标签对数据进行分割,窗口滑动步长为1。其中f为预测长度,表示向后预测f个单位时间的应变数据。
步骤5、特征融合
由于以往的时间序列模型仅考虑单个数据的自相关性,忽略了其他信息可能对数据产生的影响,因此将步骤4的输入特征(连续w个小时的温度和应变数据)进行融合,提取温度与应变数据间的低阶信息和高阶信息。
低阶信息为温度与应变数据的线性变换,表示为:
xl=Wl[xt,xs]+bl
高阶信息为温度与应变数据的线性变换再加上交叉项信息,表示为:
xm=Wm[xt,xs]+Wtsxtxs+bm
低阶信息和高阶信息经过sigmoid函数激活后便得到最终的特征:
x=sigmoid(xl+xm)
步骤6、训练集测试集划分
将融合后的温度和应变传感器的数据划分为训练集和测试集,其中80%的数据用于训练模型,20%的数据用于测试。
步骤7、SHUFFLE
为了避免数据投入的顺序对网络训练造成影响,增加随机性,提高模型的泛化能力,导致权重更新时的梯度过于极端,避免模型最终过拟合或者欠拟合,对分割后的训练集数据进行shuffle(即打乱每组训练数据的顺序)。
步骤8、模型训练
搭建两个LSTM层(图2为LSTM原理图),和一个全连接(Dense)层的神经网络模型,设置激活函数为tanh,同时,为了避免模型过拟合,每层LSTM设置Dropout率为0.2。在训练过程中以平均绝对误差(Mean Absolute Error)作为损失函数,以平均绝对误差最小化为优化目标,采用Adam优化器,设置每批训练的样本数(BatchSize)为128,训练集最大迭代次数(Epoch)500。
步骤8.1数据输入
将步骤7整合好的数据作为深度学习模型的输入。
步骤8.2 LSTM层和Dense层
输入的数据经过两个LSTM层和一个Dense层后输出预测应变数据。
步骤8.3计算损失函数
采用MAE作为模型的损失函数,计算Dense层输出的预测应变与真实应变间的损失。
步骤8.4梯度下降
使用Adam优化器梯度下降更新参数,使得模型的损失函数最小化。直到模型收敛或者模型到达最大之迭代次数。
步骤9、模型评估
为评估模型的预测准确度,使用绝对平均误差(MAE)对模型进行打分评价。平均绝对误差是所有单个观测值与算术平均值的偏差的绝对值的平均。平均绝对误差可以避免误差相互抵消的问题,因而可以准确反映实际预测误差的大小。其计算方法如下:
同时为了更加直观的体现模型在训练集和测试集上的误差,绘制损失曲线,其中横轴为训练集迭代次数(Epoch),纵轴为模型的损失函数值(Loss)。当训练集迭代次数达到100次时,绝对平均误差在训练集和验证集上得分都趋于稳定。
进一步地,为了检验基于LSTM温度-应变相关模型的有效性,与仅用应变的LSTM模型和RNN模型作对比,结果表明在测试集上基于LSTM温度-应变相关模型的均方误差、均方根误差和平均绝对误差均优于仅用应变的LSTM模型和RNN模型。
步骤10、桥梁损伤预警
当桥梁下一时刻的真实数据传回时,计算真实值与预测值之间的偏差,若偏差大于设定的阈值m,则进行预警。
模型案例
1、损失曲线
如损失曲线图3所示,横轴是训练集迭代次数(Epoch),可见循环次数越多,绝对平均误差越小。当迭代次数为100次时,绝对平均误差在训练集和验证集上得分都表现的很好。因此我们将迭代次数设置为100次。
2、预测值与真实值对比
设置超参数w为60,f为1的情形下,绘制了140个单位时间(小时)的预测与真实值的对比图,如图4所示。可见,预测值和真实值吻合较好。表1为前10个单位时间的预测值与真实值的具体数值。说明基于LSTM温度-应变相关模型的桥梁损伤预警方法的预测较为准确。
表1前10个单位时间的预测值与真实值的具体数值
3、与其它模型的对比
如表2,所示,为了检验基于LSTM温度-应变相关模型的有效性,与仅用应变的LSTM模型和RNN模型作对比,结果表明在测试集上基于LSTM温度-应变相关模型的均方误差、均方根误差和平均绝对误差均优于仅用应变的LSTM模型和RNN模型。
表2本发明模型与其他模型对比
模型 | mse | rmse | mae |
RNN | 10.24 | 3.20 | 2.58 |
LSTM | 13.73 | 3.70 | 2.69 |
MF-LSTM(ours) | 7.51 | 2.74 | 2.20 |
以上是本发明的较佳实施例,凡依本发明技术方案所作的改变,所产生的功能作用未超出本发明技术方案的范围时,均属于本发明的保护范围。
Claims (6)
1.一种基于LSTM温度-位移相关模型的桥梁损伤预警方法,其特征在于,包括如下步骤:
S1、对桥梁监测数据进行预处理;
S2、对预处理后的数据进行相关性分析,提取出相关性最强的温度传感器数据和应变传感器数据;
S3、遍历基于步骤S2提取的温度传感器数据和应变传感器数据,进行数据分割;
S4、融合温度传感器数据和应变传感器数据,提取两个传感器数据间的低阶信息和高阶信息,并得到融合后的温度-应变传感器数据;
S5、将融合后的温度-应变传感器数据划分为训练集和测试集,并打乱每组训练集中数据的顺序;
S6、搭建具有两个LSTM层和一个Dense层的神经网络模型,输入训练集数据进行训练,得到基于LSTM温度-应变相关模型;
S7、根据当前时刻桥梁数据通过基于LSTM温度-应变相关模型进行预测,当桥梁下一时刻的真实数据传回时,计算真实值与预测值之间的偏差,判断偏差是否大于设定的阈值,从而实现桥梁损伤预警;
所述步骤S2具体实现方式为:
S21、对数据表中各个温度传感器和应变传感器进行相关性分析,计算两两传感器之间的皮尔逊相关系数;根据相关系数的强弱,找出相关性最强的温度传感器和应变传感器;
S22、筛选出相关性最强的温度传感器和应变传感器后,将两个传感器的数据按时间进行内连接,其数据结构为【时间、温度传感器、应变传感器】,其对应数据为传感器单位时间的平均值;
所述步骤S3具体实现方式为:
将步骤S2提取的相关性最强的温度传感器数据和应变传感器数据中温度传感器的数据记为:
其中表示温度传感器数据中第i个数据,i∈{1,2…n};
应变传感器的数据记为:
其中表示应变传感器数据中第i个数据,i∈{1,2…n};
遍历温度传感器数据和应变传感器数据,设置滑动窗口的长度为w,提取连续w个单位时间的温度数据和应变数据作为输入特征,第w+1到w+f单位时间的应变数据作为标签对数据进行分割,窗口滑动步长为1,其中f为预测长度,表示向后预测f个单位时间的应变数据;
所述步骤S4具体实现方式为:
将连续w个小时的温度数据和应变数据进行融合,提取温度与应变数据间的低阶信息和高阶信息,其中:
低阶信息为温度与应变数据的线性变换,表示为:
xl=Wl[xt,xs]+bl
其中Wl为可训练参数矩阵,xt和xs分别为温度数据与应变数据,bl为偏置量;
高阶信息为温度与应变数据的线性变换再加上交叉项信息,表示为:
xm=Wm[xt,xs]+Wtsxtxs+bm
其中Wm为可训练参数矩阵,xt和xs分别为温度数据与应变数据,bm为偏置量;
低阶信息和高阶信息经过sigmoid函数激活后便得到最终的特征:
x=sigmoid(xl+xm)
xl和xm为计算出的低阶信息和高阶信息。
2.根据权利要求1所述的一种基于LSTM温度-位移相关模型的桥梁损伤预警方法,其特征在于,所述步骤S1具体实现方式为:
S11、准备桥梁监测数据原始数据集,其数据结构为【时间、桥梁名称、传感器类型、传感器编号、对应时间的传感器数据平均值、传感器位置】;将桥梁监测数据原始数据集按时间属性进行排序;
S12、将桥梁监测数据原始数据集整理成时间为唯一索引,列名为传感器编号的数据表,其数据为对应时间的传感器数据单位时间的平均值;将数据整理成数据表后,对单位时间间隔的时间缺失进行填补,将该时间所对应的缺失数据使用空值来进行替代;
S13、对所有接收到的传感器数据进行筛选,筛选出不在量程、精度误差范围内的所有异常值并剔除,剔除完成后在对应位置使用空值进行替代;
S14、对于数据中存在的空值,将空值用对应列的均值填补;
S15、对温度传感器和应变传感器的数据做归一化的处理,将数据映射到0和1之间:
其中x为待归一化的数据,x_min和x_max分别表示待归一化数据中的最小值和最大值。
3.根据权利要求1所述的一种基于LSTM温度-位移相关模型的桥梁损伤预警方法,其特征在于,所述步骤S5具体实现方式为:
将融合后的温度-应变传感器数据划分为训练集和测试集,其中80%的数据用于训练,20%的数据用于测试,并对训练集数据进行shuffle,打乱每组训练集中数据的顺序。
4.根据权利要求1所述的一种基于LSTM温度-位移相关模型的桥梁损伤预警方法,其特征在于,所述步骤S6具体实现方式为:
搭建两个LSTM层,和一个全连接Dense层的神经网络模型,设置激活函数为tanh,同时,为避免模型过拟合,每层LSTM设置Dropout率为0.2;在训练过程中以平均绝对误差MAE作为损失函数,以平均绝对误差最小化为优化目标,采用Adam优化器,设置每批训练的样本数BatchSize为128,训练集最大迭代次数Epoch为500,模型训练过程如下:
S61、数据输入
将步骤S5的数据作为神经网络模型的输入;
S62、LSTM层和Dense层
输入的数据经过两个LSTM层和一个Dense层后输出预测应变数据;
S63、计算损失函数
采用MAE作为模型的损失函数,计算Dense层输出的预测应变与真实应变间的损失;
S64、梯度下降
使用Adam优化器梯度下降更新参数,使得模型的损失函数最小化,直到模型收敛或者模型到达最大之迭代次数,得到训练后的基于LSTM温度-应变相关模型。
5.根据权利要求4所述的一种基于LSTM温度-位移相关模型的桥梁损伤预警方法,其特征在于,所述MAE公式如下:
其中X为样本空间,h为训练的模型,m为样本总数,h(xi)表示第i个样本经过模型h输出的预测值,yi表示第i个样本的真实值。
6.根据权利要求4所述的一种基于LSTM温度-位移相关模型的桥梁损伤预警方法,其特征在于,步骤S63中,为更加直观的体现模型在训练集和测试集上的误差,绘制损失曲线,其中横轴为训练集迭代次数Epoch,纵轴为模型的损失函数值Loss。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111052885.6A CN113723010B (zh) | 2021-09-08 | 2021-09-08 | 一种基于lstm温度-位移相关模型的桥梁损伤预警方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111052885.6A CN113723010B (zh) | 2021-09-08 | 2021-09-08 | 一种基于lstm温度-位移相关模型的桥梁损伤预警方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113723010A CN113723010A (zh) | 2021-11-30 |
CN113723010B true CN113723010B (zh) | 2024-01-26 |
Family
ID=78682750
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111052885.6A Active CN113723010B (zh) | 2021-09-08 | 2021-09-08 | 一种基于lstm温度-位移相关模型的桥梁损伤预警方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113723010B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114418071A (zh) * | 2022-01-24 | 2022-04-29 | 中国光大银行股份有限公司 | 循环神经网络训练方法 |
CN114925518B (zh) * | 2022-05-17 | 2024-04-26 | 重庆大学 | 一种基于时间序列和多传感器融合的桥梁健康预警方法 |
CN115031620B (zh) * | 2022-06-07 | 2023-05-30 | 山东高速工程检测有限公司 | 基于无线低功耗多通道采集技术的桥梁监测方法及装置 |
CN115329812B (zh) * | 2022-08-10 | 2023-07-21 | 贵州桥梁建设集团有限责任公司 | 一种基于人工智能的桥梁基础设施异常监测方法 |
CN115659243B (zh) * | 2022-12-22 | 2023-04-28 | 四川九通智路科技有限公司 | 基于mems的基础设施风险监测方法及监测*** |
CN116542146A (zh) * | 2023-01-06 | 2023-08-04 | 中路高科交通检测检验认证有限公司 | 桥梁监测温度场-应变场时空相关模型及健康诊断方法 |
CN116792155B (zh) * | 2023-06-26 | 2024-06-07 | 华南理工大学 | 一种基于分布式光纤传感的隧道健康状态监测预警方法 |
CN117216844B (zh) * | 2023-09-12 | 2024-03-26 | 汕头大学 | 一种桥梁结构损伤检测方法、***和存储介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107858883A (zh) * | 2017-11-29 | 2018-03-30 | 北京交通大学 | 一种轨道***安全状态综合监测及智能分析方法 |
CN109033509A (zh) * | 2018-06-14 | 2018-12-18 | 东南大学 | 一种基于大坝变形分析的融合模型构建方法 |
CN110909762A (zh) * | 2019-10-15 | 2020-03-24 | 北京航空航天大学 | 基于多传感器融合的机器人姿态识别方法及装置 |
CN111860839A (zh) * | 2020-07-28 | 2020-10-30 | 上海海事大学 | 基于多信号融合及Adam优化算法的岸桥故障监测方法 |
CN113011763A (zh) * | 2021-03-29 | 2021-06-22 | 华南理工大学 | 一种基于时空图卷积注意力的桥梁损伤识别方法 |
-
2021
- 2021-09-08 CN CN202111052885.6A patent/CN113723010B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107858883A (zh) * | 2017-11-29 | 2018-03-30 | 北京交通大学 | 一种轨道***安全状态综合监测及智能分析方法 |
CN109033509A (zh) * | 2018-06-14 | 2018-12-18 | 东南大学 | 一种基于大坝变形分析的融合模型构建方法 |
CN110909762A (zh) * | 2019-10-15 | 2020-03-24 | 北京航空航天大学 | 基于多传感器融合的机器人姿态识别方法及装置 |
CN111860839A (zh) * | 2020-07-28 | 2020-10-30 | 上海海事大学 | 基于多信号融合及Adam优化算法的岸桥故障监测方法 |
CN113011763A (zh) * | 2021-03-29 | 2021-06-22 | 华南理工大学 | 一种基于时空图卷积注意力的桥梁损伤识别方法 |
Non-Patent Citations (4)
Title |
---|
基于桥梁健康监测数据的异常检测与数据质量评估方法研究;陈臻;重庆大学硕士学位论文;第3.3.1节 * |
基于长短时记忆神经网络的大跨拱桥温度-位移相关模型建立方法;郑秋怡等;《工程力学》;第1-2节 * |
郑秋怡等.基于长短时记忆神经网络的大跨拱桥温度-位移相关模型建立方法.《工程力学》.2021,第1-2节. * |
面向桥梁长期健康监测的大数据处理技术及应用;涂成枫;刘泽佳;张舸;周立成;陈溢涛;程楠;辜家伟;董守斌;邓志华;王勇;汤立群;;实验力学(第05期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113723010A (zh) | 2021-11-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113723010B (zh) | 一种基于lstm温度-位移相关模型的桥梁损伤预警方法 | |
CN109308522B (zh) | 一种基于循环神经网络的gis故障预测方法 | |
CN116757534B (zh) | 一种基于神经训练网络的智能冰箱可靠性分析方法 | |
CN111460728B (zh) | 一种工业设备剩余寿命预测方法、装置、存储介质及设备 | |
CN107886168A (zh) | 一种使用多层感知器神经网络进行电梯故障识别的方法 | |
CN111737909B (zh) | 基于时空图卷积网络的结构健康监测数据异常识别方法 | |
CN110210169B (zh) | 一种基于lstm的盾构机故障预测方法 | |
CN110837888A (zh) | 一种基于双向循环神经网络的交通缺失数据补全方法 | |
CN111046961A (zh) | 基于双向长短时记忆单元和胶囊网络的故障分类方法 | |
CN116737510B (zh) | 一种基于数据分析的键盘智能监测方法及*** | |
CN110580213A (zh) | 一种基于循环标记时间点过程的数据库异常检测方法 | |
Ye et al. | A deep learning-based method for automatic abnormal data detection: Case study for bridge structural health monitoring | |
CN117032165A (zh) | 一种工业设备故障诊断方法 | |
CN114548494A (zh) | 一种可视化造价数据预测智能分析*** | |
CN113988210A (zh) | 结构监测传感网失真数据修复方法、装置及存储介质 | |
Li et al. | Structural health monitoring data anomaly detection by transformer enhanced densely connected neural networks | |
CN116894180B (zh) | 一种基于异构图注意力网络的产品制造质量预测方法 | |
CN117809164A (zh) | 基于多模态融合的变电站设备故障检测方法及*** | |
CN115761654B (zh) | 一种车辆重识别方法 | |
CN116663126A (zh) | 一种基于通道注意力BiLSTM模型的桥梁温度效应预测方法 | |
CN116842358A (zh) | 一种基于多尺度卷积和自适应特征融合的软测量建模方法 | |
CN116720095A (zh) | 一种基于遗传算法优化模糊c均值的电特性信号聚类方法 | |
CN116245212A (zh) | 基于pca-lstm的电力数据异常检测与预测方法及*** | |
CN114372640A (zh) | 一种基于波动序列分类校正的风电功率预测方法 | |
CN114722721A (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 |