CN113688773B - 一种基于深度学习的储罐穹顶位移数据修复方法及其装置 - Google Patents

一种基于深度学习的储罐穹顶位移数据修复方法及其装置 Download PDF

Info

Publication number
CN113688773B
CN113688773B CN202111030325.0A CN202111030325A CN113688773B CN 113688773 B CN113688773 B CN 113688773B CN 202111030325 A CN202111030325 A CN 202111030325A CN 113688773 B CN113688773 B CN 113688773B
Authority
CN
China
Prior art keywords
displacement data
measuring point
data
1dcnn
storage tank
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
Application number
CN202111030325.0A
Other languages
English (en)
Other versions
CN113688773A (zh
Inventor
陈增顺
华建民
黄乐鹏
薛暄译
高霖
王思宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing University
Original Assignee
Chongqing University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chongqing University filed Critical Chongqing University
Priority to CN202111030325.0A priority Critical patent/CN113688773B/zh
Publication of CN113688773A publication Critical patent/CN113688773A/zh
Application granted granted Critical
Publication of CN113688773B publication Critical patent/CN113688773B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/02—Preprocessing
    • 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/04—Architecture, e.g. interconnection topology
    • G06N3/045—Combinations of 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
    • G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/08—Feature extraction

Abstract

本发明涉及一种基于深度学习的储罐穹顶位移数据修复方法，属于位移数据修复技术领域，所述修复方法包括：获取待修复测点和周围关联测点的历史位移数据，通过集合经验模态分解算法分解为多个本征模态函数分量对1DCNN‑BiLSTM模型进行训练得到EEMD‑1DCNN‑BiLSTM模型，并通过EEMD‑1DCNN‑BiLSTM模型预测出缺失的位移数据，完成数据修复。本发明中，将EEMD、1DCNN和BiLSTM组合成一个新的模型，非常适合处理有空间关联的复杂长期时序动态信息，能够大大提高预测精度，非常适合用于LNG储罐穹顶缺失位移数据的修复。

Description

一种基于深度学习的储罐穹顶位移数据修复方法及其装置
技术领域
本发明属于储罐穹顶位移数据修复技术领域，涉及一种基于深度学习的储罐穹顶位移数据修复方法及其装置。
背景技术
采用位移传感器对于评估LNG储罐结构动力响应具有重要意义。但是在振动台实验中，某些位移传感器会发生失效或者异常从而导致数据丢失，这些数据很难恢复。
现有的基于人工智能方法对LNG储罐结构位移进行预测方法主要分为两种。一种是“浅层”机器学习方法，加速度传感数据具有高度非线性和非高斯性，“浅层”的模型对位移响应的长期预测具有一定的局限性，无法处理海量的监测数据且准确率较低。另一种方法是传统的深度神经网络模型，具有普遍性、效率高等特点，但准确性有待进一步的提高。因此，现有的预测方法无法用于修复传感器的位移数据。
发明内容
有鉴于此，本发明的目的在于提供一种基于EEMD-1DCNN-BiLSTM模型深度学习的储罐穹顶位移数据修复方法及装置。
为达到上述目的，本发明提供如下技术方案：
一种基于深度学习的储罐穹顶位移数据修复方法，用于LNG储罐穹顶位移数据缺失时的修复，包括以下步骤：
步骤S1、将LNG储罐穹顶位移数据缺失的测点作为待修复测点，在待修复测点周围选取多个测点作为关联测点，并获取待修复测点在数据缺失时段前某一时段的历史位移数据以及关联测点在对应时段的历史位移数据，通过集合经验模态分解算法将上述测点的历史位移数据分别分解为多个本征模态函数分量；
步骤S2、将分解得到的本征模态函数分量作为1DCNN-BiLSTM模型的输入特征，通过1DCNN模型提取关联测点的位移与待修复测点的位移的空间相关性特征，再把提取的空间相关性特征送入到BiLSTM模型，获取时间上的依赖度特性；
步骤S3、定义1DCNN-BiLSTM模型的损失函数，当1DCNN-BiLSTM模型的损失函数的值收敛为一个固定值保持不变时，结束训练，得到EEMD-1DCNN-BiLSTM模型；
步骤S4、获取待修复测点在数据缺失时段前的位移数据以及关联测点在数据缺失时段及数据缺失时段前的位移数据，将上述位移数据输入EEMD-1DCNN-BiLSTM模型预测出待修复测点在数据缺失时段的位移数据；将预测出的位移数据作为待修复测点在数据缺失时段的位移数据，完成缺失数据的修复。
进一步的，集合经验模态分解算法通过以下步骤实现：
步骤S11、选定原始信号的处理次数m；
步骤S12、选择m个不同幅值的随机白噪声，将原始信号分别与各个白噪声结合，得到m个新的信号；
步骤S13、对m个新的信号分别进行经验模态分解，得到一系列本征模态函数分量；
步骤S14、对相应模态的本征模态函数分量分别求均值，得到集合经验模态分解结果。
进一步的，m个所述白噪声的尺度呈现均匀分布状态，且其能量在频谱上也呈现均匀分布状态。
进一步的，经验模态分解算法通过以下步骤实现：
步骤S131、将原始信号与白噪声结合后得到的新的信号作为待分解信号x(t)，t表示时间；
步骤S132、对待分解信号x(t)进行筛分；具体为：找出待分解信号x(t)所有的极大值点，拟合成待分解信号的上包络线；找出待分解信号x(t)所有的极小值点，拟合成待分解信号的下包络线；计算上包络线和下包络线的均值得到待分解信号x(t)的平均包络线m1(t)；将待分解信号x(t)减去m1(t)后即可得到一个中间分量函数d1,1(t)；
步骤S133、判断d1,1(t)是否满足本征模态函数分量的条件，如果不满足，则以d1,1(t)代替待分解信号x(t)，继续按步骤S12对d1,1(t)进行筛分，经过K次筛分后的信号记为d1,k(t)，直到d1,k(t)满足本征模态函数分量的条件时，记为待分解信号x(t)的第一个IMF分量IMF1(t)；
步骤S134、从待分解信号x(t)中减去第一个IMF分量IMF1(t)，得到剩余分量r1(t)，对r1(t)继续按步骤S12和步骤S13进行分解；经过n次分解后，求得残余信号rn(t)；当rn(t)为单调函数时，停止分解，将剩余分量函数rn(t)作为残余量RES。
进一步的，本征模态函数分量满足以下条件：
函数在整个时间范围内，极值点个数与过零点个数相等或相差1；
在任意时刻点，上包络线和下包络线的均值均为0。
进一步的，一个测点的位移数据的时间序列形成一维数据；所述一维数据通过EEMD分解成多个IMF序列，形成二维数据；获取多个关联测点的数据后，形成三维数据，并将所述三维数据作为1DCNN模型的输入特征映射组。
进一步的，1DCNN-BiLSTM模型的损失函数l(x，y)定义为：
其中，N表示样本的个数，xi表示第i个样本的实际值，yi表示第i个样本的预测值。
一种基于深度学习的储罐穹顶位移数据修复装置，包括：
位移数据采集模块，用于实时采集LNG储罐穹顶各测点的位移数据并传输到计算分析模块；
计算分析模块，用于实时监测LNG储罐穹顶各测点的位移数据是否有缺失，当监测到某一测点的位移数据发生缺失时，将位移数据缺失的测点作为待修复测点，并在待修复测点周围选取多个测点作为关联测点，将待修复测点的位移数据和关联测点的位移数据输入EEMD-1DCNN-BiLSTM模型，对缺失的位移数据进行预测，使用预测的位移数据对缺失的位移数据进行补全，完成数据修复；以及
输出模块，用于输出位移数据采集模块采集的各测点的位移数据，以及待修复测点的修复的位移数据。
进一步的，所述计算分析模块包括数据读取单元、监测单元、EEMD-1DCNN-BiLSTM模型和存储单元；
所述数据读取单元用于读取位移数据采集模块采集的LNG储罐穹顶各测点的位移数据；
所述监测单元用于实时监测LNG储罐穹顶各测点的位移数据是否有缺失；
所述EEMD-1DCNN-BiLSTM模型包括集合经验模态分解单元和1DCNN-BiLSTM模型，所述集合经验模态分解单元用于通过集合经验模态分解算法将待修复测点和各关联测点的位移数据分别分解为多个本征模态函数分量，并将各本征模态函数分量形成的向量作为1DCNN-BiLSTM模型的输入特征；所述1DCNN-BiLSTM模型用于根据输入特征预测出待修复测点缺失的位移数据；
所述存储单元用于存储位移数据采集模块采集的LNG储罐穹顶各测点的位移数据，以及EEMD-1DCNN-BiLSTM模型预测的位移数据。
进一步的，还包括检测单元和报警单元，所述检测单元用于指定一位移数据正常的测点为待修复测点，并将EEMD-1DCNN-BiLSTM模型预测的该测点的位移数据与位移数据采集模块采集的该测点的真实位移数据进行比较，当两者的差值超出预先设置的阈值时，使报警单元输出报警信号，提示预测数据偏移过大。
本发明中，将EEMD、1DCNN和BiLSTM组合成一个新的模型，其中，EEMD能使复杂非线性位移数据分解为有限个频率由高到低的本征模态函数的线性组合，并且分解出来的各IMF分量包含了原信号的不同时间尺度的局部特征信号；1DCNN模型具有局部连接、权值共享等特点，能够保留、提取IMF之间的空间相关性特征；BiLSTM模型能够充分挖掘变量之间的非线性关系，自适应地感知上下时间序列特性信息，并能在任何时间点保存过去和将来的信息，从而具有捕获前后信息特征的能力，非常适合处理复杂的问题。因此，组合的新模型非常适合处理有空间关联的复杂长期时序动态信息，能够大大提高预测精度，非常适合用于LNG储罐穹顶缺失位移数据的修复。另外，EEMD算法与1DCNN-BiLSTM模型对硬件的要求不高，实现成本低。
附图说明
为了使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明作优选的详细描述，其中：
图1为本发明一种基于深度学习的储罐穹顶位移数据修复方法的一个优选实施例的流程图；
图2为获取待修复测点和关联测点的数据进行预测的示意图；
图3为通过EEMD算法对数据进行分解的示意图；
图4为1DCNN-BiLSTM模型的结构示意图；
图5为一维卷积神经网络计算过程的示意图；
图6为LSTM的单个神经元体系的结构示意图；
图7为本发明一种基于深度学习的储罐穹顶位移数据修复装置的一个优选实施例的结构框图。
具体实施方式
以下通过特定的具体实例说明本发明的实施方式，本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用，本说明书中的各项细节也可以基于不同观点与应用，在没有背离本发明的精神下进行各种修饰或改变。需要说明的是，以下实施例中所提供的图示仅以示意方式说明本发明的基本构想，在不冲突的情况下，以下实施例及实施例中的特征可以相互组合。
如图1所示，本发明一种基于深度学习的储罐穹顶位移数据修复方法的一个优选实施例包括以下步骤：
步骤S1、将LNG储罐穹顶位移数据缺失的测点作为待修复测点，在待修复测点周围选取多个测点作为关联测点，并获取待修复测点在数据缺失时段前某一时段的历史位移数据，以及关联测点在对应时段的历史位移数据。如图2所示，假设待修复测点缺失了时间步T到时间步(T+1)之间的位移数据，则获取待修复测点在时间步T之前K个时间步的位移数据，以及关联测点从时间步(T-K)到时间步(T+1)之间的位移数据。
之后，通过EEMD(EnsembleEmpirical Mode Decomposition，集合经验模态分解)算法，将上述测点的历史位移数据分别分解为多个IMF(Intrinsic Mode Functions，本征模态函数)分量。IMF分量满足以下两个条件：
(1)函数在整个时间范围内，极值点个数与过零点个数相等或相差1；
(2)在任意时刻点，上包络线和下包络线的均值均为0。
如图3所示，集合经验模态分解算法通过以下步骤实现：
步骤S11、选定原始信号的处理次数m；
步骤S12、选择m个不同幅值的随机白噪声，将原始信号分别与各个白噪声结合，得到m个新的信号n1(t)，n2(t)，……，nj(t)，……，nm(t)；m个所述白噪声的尺度呈现均匀分布状态，且其能量在频谱上也呈现均匀分布状态。
步骤S13、对m个新的信号分别进行EMD(Empirical Mode Decomposition，集合经验模态分解)，得到一系列IMF分量。下面以对n1(t)进行EMD为例对EMD算法的具体过程进行说明：
步骤S131、将n1(t)作为EMD的待分解信号x(t)。
步骤S132、对待分解信号x(t)进行筛分。筛分的过程就是从待分解信号中减去它的平均包络线函数，得到一个新的函数；具体为：找出待分解信号x(t)所有的极大值点，用三次样条函数拟合成待分解信号x(t)的上包络线；找出待分解信号x(t)所有的极小值点，用三次样条函数拟合成待分解信号x(t)的下包络线；计算上包络线和下包络线的均值得到第一个平均包络线函数m1(t)；将待分解信号x(t)减去第一个平均包络线函数m1(t)，得到第一个中间分量函数d1,1(t)。
步骤S133、判断中间分量函数d1,1(t)是否满足IMF分量的两个条件，如果满足则将d1,1(t)记为待分解信号的第一个IMF分量IMF1(t)；如果不满足则继续按步骤S12对d1,1(t)进行筛分，直至中间分量函数满足IMF分量的条件。假设K次筛分后得到的中间分量函数d1,k(t)满足IMF分量的条件，则将d1,k(t)记为待分解信号的第一个IMF分量，对于n1(t)，其第一个IMF分量记为a1,1。
步骤S134、从待分解信号x(t)中减去第一个IMF分量，得到第一个剩余分量函数r1(t)；将第一个剩余分量函数r1(t)继续按步骤S12和步骤S13进行分解(分解就是通过反复筛分的方式从信号中分解出IMF分量)，得到第二个IMF分量，对于n1(t)，其第二个IMF分量记为a2,1；使用第一个剩余分量函数r1(t)减去第二个IMF分量，得到第二个剩余分量函数r2(t)。继续按步骤S12和步骤S13对第二个剩余分量函数r2(t)进行分解；假设经过N次分解后，得到的第N个剩余分量函数rn(t)为单调函数，则停止分解，将剩余分量函数rn(t)作为残余量RES。如图2所示，此时，将信号n1(t)分解成了n个IMF分量(即a1,1、a2,1、……、ai,1、……、aN,1)和一个残余量RES。
按照以上方法，将n2(t)通过EMD分解为即a1,2、a2,2、……、ai,2、……、aN,2；
……；
将nj(t)通过EMD分解为即a1,j、a2,j、……、ai,j、……、aN,j；
……；
将nm(t)通过EMD分解为即a1,m、a2,m、……、ai,m、……、aN,m。
步骤S14、对相应模态的IMF分量分别求均值，得到a1、a2、……、ai、……、aN作为最终的IMF分量，即为集合经验模态分解结果。求平均值的公式为：
从EMD的分解过程中可以看出，与傅里叶变换和小波分解相比较，EMD不需要设定基函数，具有自适应性，因此适用范围更广泛。将待分解信号x(t)分解后，第一个IMF分量包含待分解信号x(t)中的时间尺度最小(频率最高)的成分，随着IMF分量阶数的增加，其对应的频率成分逐渐降低，rn(t)(即本实施例中的残余量RES)的频率成分最低。根据EMD分解的收敛条件，分解得到的残余量rn(t)为单调函数时，其时间周期将大于信号的记录长度，因此可以将残余量rn(t)作为待分解信号x(t)的趋势项。
步骤S2、将分解得到的本征模态函数分量(即a1、a2、……、ai、……、aN)作为1DCNN-BiLSTM模型的输入特征，通过1DCNN(一维卷积神经网络)模型提取关联测点的位移与待修复测点的位移的空间相关性特征，再把提取的空间相关性特征送入到BiLSTM(bidirectional long-short term memory；双向长短期记忆网络)模型，获取时间上的依赖度特性。如图4所示，所述1DCNN-BiLSTM模型由1DCNN模型和BiLSTM模型拼接而成。
CNN被广泛应用于图像处理领域，通常来说，图像是三维数据即X·Y·Z，本实施例中使用CNN处理时间序列的信号问题，因此需要将时间序列的信号数据变成三维数据。一个传感器(即测点)的位移数据的时间序列是1D，通过EEMD分解成多个IMF序列，变成2D数据，同时处理多个传感器，将会再增加一个维度，因此最终可以视为3D数据处理，从面以类似于图像数据的方式进行处理。1DCNN包括卷积层和池化层，其工作原理如下。
卷积层的作用是提取一个局部区域的特征，不同的卷积核相当于不同的特征提取器。卷积层的神经元和全连接网络一样都是一维结构。由于卷积网络主要应用在图像处理上，而图像为二维结构，因此为了更充分地利用图像的局部信息，通常将神经元组织为三维结构的神经层，其大小为高度q×宽度p×深度D，看成是由D个q×p大小的二维结构的特征映射构成。
特征映射(Feature Map)为经过卷积提取到的特征，每个特征映射可以作为一类抽取的特征。为了提高卷积网络的表示能力，可以在每一层使用多个不同的特征映射，以更好地表示特征。
不失一般性，假设一个卷积层的结构如下：
(1)输入特征映射组：X∈Rq×p×D为三维张量(Tensor)；其中，q表示一个传感器的位移数据包括的时间步个数，p表示传感器位移数据经过一次EEMD分解得到的IMF分量个数，D表示传感器的个数。每个切片(Slice)矩阵Xd∈Rq×p为一个输入特征映射，1≤d≤D。
(2)输出特征映射组：Y∈Rq′×p′×L为三维张量，其中每个切片矩阵CLt∈Rq′×p′为一个输出特征映射，1≤t≤L；Y＝[CL1，CL2，……，CLL]。
(3)卷积核：W∈RU×V×D×L为四维张量；其中，U表示卷积核的行数，V表示卷积核的列数，例如：U×V可取值3×5。每个切片矩阵Wt,d∈RU×V为一个二维卷积核，1≤d≤D；1≤t≤L。
为了计算输出特征映射Y，用卷积核Wt,1、Wt,2、……、Wt,D分别对输入特征映射X1、X2、……、XD进行卷积，然后将卷积结果相加，并加上一个标量偏置b得到卷积层的净输入Zt，这里净输入是指没有经过非线性激活函数的净活性值(Net Activation)。
再经过非线性激活函数后得到输出特征映射Yt。
CLt＝f(Zt)
其中为三维卷积核；bt表示偏置矩阵；f()为非线性激活函数，一般用ReLU函数。计算过程如图5所示，图中的虚线框表示卷积核。如果希望卷积层输出L个特征映射，可以将上述计算过程重复L次，得到L个特征映射CL1，CL2，……，CLL，L个特征映射经过池化后得到L个输出特征SL1，SL2，……，SLL。在输入为X∈Rq×p×D，输出为Y∈Rq′×p′×L的卷积层中，每一个输出特征映射都需要D个滤波器以及一个偏置。假设每个滤波器的大小为U×V，那么共需要L×D×(U×V)+L个参数。
从上述计算过程可以看出，1DCNN模型具有局部连接、权值共享等特点，能够保留、提取IMF之间的空间相关性特征。
1DCNN的输出特征被送入BiLSTM。BiLSTM即双向LSTM，由两个单独的LSTM(即前向LSTM和后向LSTM)组合而成，以两种方式对输入特征进行处理，一种方式是从过去到将来，另一种方式是从将来到过去，这种方法与单向LSTM的不同之处在于，在向后运行的LSTM中，保留了未来的信息并结合使用两个隐藏状态，可以在任何时间点保存过去和将来的信息，从而具有捕获前后信息特征的能力，能够处理非常复杂的问题。其计算公式如下所示：
其中，xt表示BiLSTM在时间t时刻的输入特征，即1DCNN对应时刻的输出特征SLt；表示t时刻的向前传播隐层状态；/>表示(t-1)时刻的向前传播隐层状态；/>表示t时刻的向后传播隐层状态；/>表示(t+1)时刻的向后传播隐层状态；Ot表示t时刻的隐层状态；αt为t时刻前向传播LSTM单元隐层输出权重；βt为t时刻后向传播LSTM单元隐层输出的权重；bt为t时刻隐层状态所对应的偏置量。
如图6所示，LSTM的单个神经单元的体系结构包括输入门、遗忘门、输出门以及记忆单元，用于实现信息的输入和输出，其运算过程如下：
Γi＝σ(Wi,xxt+Wi,hht-1+bi)
Γf＝σ(Wf,xxt+Wf,hht-1+bf)
Γo＝σ(Wo,xxt+Wo,hht-1+bo)
ht＝Γo*tanh(Ct)
其中，Wi,x、Wi,h、Wf,x、Wf,h、Wo,x、Wo,h、Wc,x、Wc,h表示权重矩阵；bi、bf、bc、bo表示偏置矩阵；xt表示时间t时刻的输入特征，即1DCNN对应时刻的输出特征SLt；ct-1表示更新前的神经元；ct表示更新后的神经元；ht-1表示上一时刻(前向传输时为t-1时刻，后向传输时为t+1时刻)的输出特征；ht表示当前时刻(即t时刻)的输出特征；Γi表示输入门；Γf表示遗忘门；Γo表示输出门；为候选神经元；σ为Sigmoid函数；tanh为双曲正切函数。
步骤S3、定义1DCNN-BiLSTM模型的损失函数，当1DCNN-BiLSTM模型的损失函数的值收敛为一个固定值保持不变时，结束训练，得到EEMD-1DCNN-BiLSTM模型。
其中，1DCNN-BiLSTM模型的损失函数l(x，y)可定义为：
其中，N表示样本的个数，xi表示第i个样本的实际值(即真实值)，yi表示第i个样本的预测值。
当损失函数的值收敛为一个固定值保持不变时，认为此时1DCNN-BiLSTM模型的参数为最优的模型参数，停止模型训练。
步骤S4、获取待修复测点在数据缺失时段前的位移数据，以及关联测点在数据缺失时段及数据缺失时段前的位移数据，将上述位移数据输入EEMD-1DCNN-BiLSTM模型预测出待修复测点在数据缺失时段的位移数据；将预测出的位移数据作为待修复测点在数据缺失时段的位移数据，完成缺失数据的修复。
本发明还提供一种基于深度学习的储罐穹顶位移数据修复装置，如附图7所示，本发明一种基于深度学习的储罐穹顶位移数据修复装置的一个优选实施例包括位移数据采集模块、计算分析模块和输出模块。
所述位移数据采集模块用于实时采集LNG储罐穹顶各测点的位移数据并传输到计算分析模块；
所述计算分析模块用于实时监测LNG储罐穹顶各测点的位移数据是否有缺失，当监测到某一测点的位移数据发生缺失时，将位移数据缺失的测点作为待修复测点，并在待修复测点周围选取多个测点作为关联测点，将待修复测点的位移数据和关联测点的位移数据输入EEMD-1DCNN-BiLSTM模型，对缺失的位移数据进行预测，使用预测的位移数据对缺失的位移数据进行补全，完成数据修复。
所述计算分析模块包括数据读取单元、监测单元、EEMD-1DCNN-BiLSTM模型和存储单元。
所述数据读取单元用于读取位移数据采集模块采集的LNG储罐穹顶各测点的位移数据；优选为采用包括GPS数据采集单元和/或北斗定位数据采集单元的模块。
所述监测单元用于实时监测LNG储罐穹顶各测点的位移数据是否有缺失；
所述EEMD-1DCNN-BiLSTM模型包括集合经验模态分解单元和1DCNN-BiLSTM模型，所述集合经验模态分解单元用于通过集合经验模态分解算法将待修复测点和各关联测点的位移数据分别分解为多个本征模态函数分量，并将各本征模态函数分量形成的向量作为1DCNN-BiLSTM模型的输入特征；所述1DCNN-BiLSTM模型用于根据输入特征预测出待修复测点缺失的位移数据。
所述存储单元用于存储位移数据采集模块采集的LNG储罐穹顶各测点的位移数据，以及EEMD-1DCNN-BiLSTM模型预测的位移数据。
所述输出模块用于输出位移数据采集模块采集的各测点的位移数据，以及待修复测点的修复的位移数据。所述预测数据输出模块优选为采用可视化模块，例如显示器，以可视的方式输出历史数据和预测数据。
为检测模型的预测偏差是否过大，还可包括检测单元和报警单元，所述检测单元用于指定一位移数据正常的测点为待修复测点，通过EEMD-1DCNN-BiLSTM模型预测出该测点的位移数据，并将预测的位移数据与位移数据采集模块采集的该测点的真实位移数据进行比较，当两者的差值超出预先设置的阈值时，使报警单元输出报警信号，提示预测数据偏移过大；从而提醒操作人员可能需要重新对模型进行训练，以提高模型的预测准确度。当然，操作人员不对报警进行处理也不会影响模型工作。
本实施例中，分析模块采用EEMD算法与1DCNN-BiLSTM模型，对计算和存储能力要求不高，对硬件的要求不高，实现成本低。
最后说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本技术方案的宗旨和范围，其均应涵盖在本发明的权利要求范围当中。

Claims (10)

1.一种基于深度学习的储罐穹顶位移数据修复方法，用于LNG储罐穹顶位移数据缺失时的修复，其特征在于，包括以下步骤：
步骤S1、将LNG储罐穹顶位移数据缺失的测点作为待修复测点，在待修复测点周围选取多个测点作为关联测点，并获取待修复测点在数据缺失时段前某一时段的历史位移数据以及关联测点在对应时段的历史位移数据，通过集合经验模态分解算法将上述测点的历史位移数据分别分解为多个本征模态函数分量；
步骤S2、将分解得到的本征模态函数分量作为1DCNN-BiLSTM模型的输入特征，通过1DCNN模型提取关联测点的位移与待修复测点的位移的空间相关性特征，再把提取的空间相关性特征送入到BiLSTM模型，获取时间上的依赖度特性；
步骤S3、定义1DCNN-BiLSTM模型的损失函数，当1DCNN-BiLSTM模型的损失函数的值收敛为一个固定值保持不变时，结束训练，得到EEMD-1DCNN-BiLSTM模型；
步骤S4、获取待修复测点在数据缺失时段前的位移数据以及关联测点在数据缺失时段及数据缺失时段前的位移数据，将上述位移数据输入EEMD-1DCNN-BiLSTM模型预测出待修复测点在数据缺失时段的位移数据；将预测出的位移数据作为待修复测点在数据缺失时段的位移数据，完成缺失数据的修复。
2.根据权利要求1所述的一种基于深度学习的储罐穹顶位移数据修复方法，其特征在于，集合经验模态分解算法通过以下步骤实现：
步骤S11、选定原始信号的处理次数m；
步骤S12、选择m个不同幅值的随机白噪声，将原始信号分别与各个白噪声结合，得到m个新的信号；
步骤S13、对m个新的信号分别进行经验模态分解，得到一系列本征模态函数分量；
步骤S14、对相应模态的本征模态函数分量分别求均值，得到集合经验模态分解结果。
3.根据权利要求2所述的一种基于深度学习的储罐穹顶位移数据修复方法，其特征在于，m个所述白噪声的尺度呈现均匀分布状态，且其能量在频谱上也呈现均匀分布状态。
4.根据权利要求2所述的一种基于深度学习的储罐穹顶位移数据修复方法，其特征在于，经验模态分解算法通过以下步骤实现：
步骤S131、将原始信号与白噪声结合后得到的新的信号作为待分解信号x(t)，t表示时间；
步骤S132、对待分解信号x(t)进行筛分；具体为：找出待分解信号x(t)所有的极大值点，拟合成待分解信号的上包络线；找出待分解信号x(t)所有的极小值点，拟合成待分解信号的下包络线；计算上包络线和下包络线的均值得到待分解信号x(t)的平均包络线m1(t)；将待分解信号x(t)减去m1(t)后即可得到一个中间分量函数d1,1(t)；
步骤S133、判断d1,1(t)是否满足本征模态函数分量的条件，如果不满足，则以d1,1(t)代替待分解信号x(t)，继续按步骤S12对d1,1(t)进行筛分，经过K次筛分后的信号记为d1,k(t)，直到d1,k(t)满足本征模态函数分量的条件时，记为待分解信号x(t)的第一个IMF分量IMF1(t)；
步骤S134、从待分解信号x(t)中减去第一个IMF分量IMF1(t)，得到剩余分量r1(t)，对r1(t)继续按步骤S12和步骤S13进行分解；经过n次分解后，求得残余信号rn(t)；当rn(t)为单调函数时，停止分解，将剩余分量函数rn(t)作为残余量RES。
5.根据权利要求4所述的一种基于深度学习的储罐穹顶位移数据修复方法，其特征在于，本征模态函数分量满足以下条件：
函数在整个时间范围内，极值点个数与过零点个数相等或相差1；
在任意时刻点，上包络线和下包络线的均值均为0。
6.根据权利要求1所述的一种基于深度学习的储罐穹顶位移数据修复方法，其特征在于，一个测点的位移数据的时间序列形成一维数据；所述一维数据通过EEMD分解成多个IMF序列，形成二维数据；获取多个关联测点的数据后，形成三维数据，并将所述三维数据作为1DCNN模型的输入特征映射组。
7.根据权利要求1所述的一种基于深度学习的储罐穹顶位移数据修复方法，其特征在于，1DCNN-BiLSTM模型的损失函数l(x，y)定义为：
其中，N表示样本的个数，xi表示第i个样本的实际值，yi表示第i个样本的预测值。
8.一种用于实现如权利要求1至7任一项所述的一种基于深度学习的储罐穹顶位移数据修复方法的基于深度学习的储罐穹顶位移数据修复装置，其特征在于，包括：
位移数据采集模块，用于实时采集LNG储罐穹顶各测点的位移数据并传输到计算分析模块；
计算分析模块，用于实时监测LNG储罐穹顶各测点的位移数据是否有缺失，当监测到某一测点的位移数据发生缺失时，将位移数据缺失的测点作为待修复测点，并在待修复测点周围选取多个测点作为关联测点，将待修复测点的位移数据和关联测点的位移数据输入EEMD-1DCNN-BiLSTM模型，对缺失的位移数据进行预测，使用预测的位移数据对缺失的位移数据进行补全，完成数据修复；以及
输出模块，用于输出位移数据采集模块采集的各测点的位移数据，以及待修复测点的修复的位移数据。
9.根据权利要求8所述的一种基于深度学习的储罐穹顶位移数据修复装置，其特征在于，所述计算分析模块包括数据读取单元、监测单元、EEMD-1DCNN-BiLSTM模型和存储单元；
所述数据读取单元用于读取位移数据采集模块采集的LNG储罐穹顶各测点的位移数据；
所述监测单元用于实时监测LNG储罐穹顶各测点的位移数据是否有缺失；
所述EEMD-1DCNN-BiLSTM模型包括集合经验模态分解单元和1DCNN-BiLSTM模型，所述集合经验模态分解单元用于通过集合经验模态分解算法将待修复测点和各关联测点的位移数据分别分解为多个本征模态函数分量，并将各本征模态函数分量形成的向量作为1DCNN-BiLSTM模型的输入特征；所述1DCNN-BiLSTM模型用于根据输入特征预测出待修复测点缺失的位移数据；
所述存储单元用于存储位移数据采集模块采集的LNG储罐穹顶各测点的位移数据，以及EEMD-1DCNN-BiLSTM模型预测的位移数据。
10.根据权利要求8所述的一种基于深度学习的储罐穹顶位移数据修复装置，其特征在于，还包括检测单元和报警单元，所述检测单元用于指定一位移数据正常的测点为待修复测点，并将EEMD-1DCNN-BiLSTM模型预测的该测点的位移数据与位移数据采集模块采集的该测点的真实位移数据进行比较，当两者的差值超出预先设置的阈值时，使报警单元输出报警信号，提示预测数据偏移过大。
CN202111030325.0A 2021-09-03 2021-09-03 一种基于深度学习的储罐穹顶位移数据修复方法及其装置 Active CN113688773B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111030325.0A CN113688773B (zh) 2021-09-03 2021-09-03 一种基于深度学习的储罐穹顶位移数据修复方法及其装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111030325.0A CN113688773B (zh) 2021-09-03 2021-09-03 一种基于深度学习的储罐穹顶位移数据修复方法及其装置

Publications (2)

Publication Number Publication Date
CN113688773A CN113688773A (zh) 2021-11-23
CN113688773B true CN113688773B (zh) 2023-09-26

Family

ID=78585144

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111030325.0A Active CN113688773B (zh) 2021-09-03 2021-09-03 一种基于深度学习的储罐穹顶位移数据修复方法及其装置

Country Status (1)

Country Link
CN (1) CN113688773B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0756746A1 (fr) * 1995-02-16 1997-02-05 Gemplus Procede de mise a jour securisee de memoire eeprom
CN103206614A (zh) * 2013-04-01 2013-07-17 中山大学 用于承压储罐的bleve事故监测的预警装置及方法
CN105760839A (zh) * 2016-02-22 2016-07-13 重庆大学 基于多特征流形学习与支持向量机的轴承故障诊断方法
CN109615014A (zh) * 2018-12-17 2019-04-12 清华大学 一种基于kl散度优化的数据分类系统与方法
CN110333076A (zh) * 2019-06-19 2019-10-15 电子科技大学 基于CNN-Stacking的轴承故障诊断方法
CN110909931A (zh) * 2019-11-20 2020-03-24 成都理工大学 基于模态分解重构及深度lstm-rnn模型的测井曲线预测方法
CN111240268A (zh) * 2020-01-14 2020-06-05 重庆大学 基于slstm神经网络的轴系统热误差建模方法和热误差补偿系统
CN112478483A (zh) * 2020-12-01 2021-03-12 华设设计集团股份有限公司 一种危化品常压储罐自动化监测预警方法及系统

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160070276A1 (en) * 2014-09-08 2016-03-10 Leeo, Inc. Ecosystem with dynamically aggregated combinations of components
US20180284745A1 (en) * 2016-05-09 2018-10-04 StrongForce IoT Portfolio 2016, LLC Methods and systems for self-organization of collected data using 3rd party data from a data marketplace in an industrial internet of things environment

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0756746A1 (fr) * 1995-02-16 1997-02-05 Gemplus Procede de mise a jour securisee de memoire eeprom
CN103206614A (zh) * 2013-04-01 2013-07-17 中山大学 用于承压储罐的bleve事故监测的预警装置及方法
CN105760839A (zh) * 2016-02-22 2016-07-13 重庆大学 基于多特征流形学习与支持向量机的轴承故障诊断方法
CN109615014A (zh) * 2018-12-17 2019-04-12 清华大学 一种基于kl散度优化的数据分类系统与方法
CN110333076A (zh) * 2019-06-19 2019-10-15 电子科技大学 基于CNN-Stacking的轴承故障诊断方法
CN110909931A (zh) * 2019-11-20 2020-03-24 成都理工大学 基于模态分解重构及深度lstm-rnn模型的测井曲线预测方法
CN111240268A (zh) * 2020-01-14 2020-06-05 重庆大学 基于slstm神经网络的轴系统热误差建模方法和热误差补偿系统
CN112478483A (zh) * 2020-12-01 2021-03-12 华设设计集团股份有限公司 一种危化品常压储罐自动化监测预警方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"An improved method based on EEMD-LSTM to predict missing measured data of structural sensors";Zengshun Chen等;《Applied sciences》;第12卷(第8期);第1-26页 *

Also Published As

Publication number Publication date
CN113688773A (zh) 2021-11-23

Similar Documents

Publication Publication Date Title
CN109947086B (zh) 基于对抗学习的机械故障迁移诊断方法及系统
CN109766930B (zh) 一种基于dcnn模型的矿井机械设备的剩余寿命预测方法
Famili et al. Data preprocessing and intelligent data analysis
CN110361778B (zh) 一种基于生成对抗网络的地震数据重建方法
CN111767930A (zh) 物联网时序数据异常检测方法及其相关设备
Aydemir et al. Image-based prognostics using deep learning approach
CN114429150A (zh) 一种基于改进深度子域适应网络的变工况下滚动轴承故障诊断方法及系统
CN115062272A (zh) 一种水质监测数据异常识别及预警方法
CN111861028A (zh) 一种基于时空数据融合的犯罪数量预测方法
CN113988357A (zh) 基于深度学习的高层建筑风致响应预测方法及装置
CN115858609A (zh) 电动汽车充电桩状态监测方法、故障鉴定方法及电子设备
CN113688773B (zh) 一种基于深度学习的储罐穹顶位移数据修复方法及其装置
CN116933152A (zh) 一种基于多维emd-pso-lstm神经网络的海浪信息预测方法及系统
CN113722308B (zh) 基于EEMD-MultiCNN-LSTM的加速度响应数据补全方法及装置
CN113743297A (zh) 基于深度学习的储罐穹顶位移数据修复方法及装置
CN113688770A (zh) 高层建筑长期风压缺失数据补全方法及装置
Colace et al. Unsupervised Learning Techniques for Vibration-Based Structural Health Monitoring Systems Driven by Data: A General Overview
CN113688774B (zh) 基于深度学习的高层建筑风致响应预测、训练方法及装置
CN116882538B (zh) 一种海洋环境预测模型的训练方法及相关装置
CN113762203B (zh) 基于仿真数据的跨域自适应sar图像分类方法、装置及设备
Nurmalasari et al. Classification for Papaya Fruit Maturity Level With Convolutional Neural Network
CN117933082A (zh) 一种基于特征分解和深度学习的大坝变形预测方法及系统
Kumar et al. Artificial intelligence and image processing approaches in damage assessment and material evaluation
Mondal et al. Forest-based Underwater Temperature Prediction Model for IoT Devices
Zhou et al. Identification of Strong Motion Record Baseline Drift Based on Bayesian Optimized Transformer Network

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