CN113705094A - 一种基于pso-gru的船舶燃油管道故障预测方法 - Google Patents

一种基于pso-gru的船舶燃油管道故障预测方法 Download PDF

Info

Publication number
CN113705094A
CN113705094A CN202110995049.5A CN202110995049A CN113705094A CN 113705094 A CN113705094 A CN 113705094A CN 202110995049 A CN202110995049 A CN 202110995049A CN 113705094 A CN113705094 A CN 113705094A
Authority
CN
China
Prior art keywords
gru
value
model
pso
gate
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.)
Granted
Application number
CN202110995049.5A
Other languages
English (en)
Other versions
CN113705094B (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.)
Jiangsu Jierui Information Technology Co Ltd
Original Assignee
Jiangsu Jierui Information Technology Co Ltd
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 Jiangsu Jierui Information Technology Co Ltd filed Critical Jiangsu Jierui Information Technology Co Ltd
Priority to CN202110995049.5A priority Critical patent/CN113705094B/zh
Publication of CN113705094A publication Critical patent/CN113705094A/zh
Application granted granted Critical
Publication of CN113705094B publication Critical patent/CN113705094B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/27Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/004Artificial life, i.e. computing arrangements simulating life
    • G06N3/006Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/044Recurrent networks, e.g. Hopfield networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/08Logistics, e.g. warehousing, loading or distribution; Inventory or stock management
    • G06Q10/083Shipping
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/14Pipes
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Artificial Intelligence (AREA)
  • Software Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Strategic Management (AREA)
  • Computing Systems (AREA)
  • General Health & Medical Sciences (AREA)
  • Computational Linguistics (AREA)
  • Mathematical Physics (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Molecular Biology (AREA)
  • Health & Medical Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Tourism & Hospitality (AREA)
  • Development Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Marketing (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • General Business, Economics & Management (AREA)
  • Game Theory and Decision Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Medical Informatics (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种基于PSO‑GRU的船舶燃油管道故障预测方法,通过准确预测船舶燃油设备供应泵压力值实现对燃油管道粗滤器和供压设备运行状态的有效评估;首先采集船舶燃油设备供应泵压力数据和船舶设备信息,经过归一化处理后划分训练集和测试集;然后构建门循环单元网络模型(GRU),并通过粒子群优化算法对GRU模型进行优化以获得最优预测模型参数,从而完成PSO‑GRU模型的搭建;最后将测试集输入训练后的PSO‑GRU模型进行预测,得到压力预测值。本发明提供的方法可以准确预测船舶燃油管道内压力值,对船舶设备故障诊断及运行状况提供理论依据,降低维护成本,保障船舶安全平稳运行。

Description

一种基于PSO-GRU的船舶燃油管道故障预测方法
技术领域
本发明涉及船舶燃油管道故障预测领域,尤其涉及一种基于PSO-GRU混合模型的船舶燃油管道故障预测方法。
背景技术
随着世界贸易的全面流通,海运已成为国家之间进出口货物的重要运输手段;随之而来的是船舶数量的与日俱增。其中以燃烧燃油为动力的船占绝大多数,故由船舶燃油管道故障引起的事故不在少数;又因为船舶燃油管道长度大、面积大、持续工作时间长、其中储存燃料有害且易燃易爆、高压高温下等特点,若是发生泄漏、炸裂等问题会产生巨大的灾害。因此,燃油管道的故障预测技术对于保证船舶安全性能与可靠性,降低维护成本而言十分重要。
在故障预测研究领域中,人工智能的技术与理论已经趋于成熟,目前得到广泛应用的算法大致可分为四种:一是决策树,其在已知故障发生概率的基础上,通过对输入数据信息的期望值进行增值,来进行预测与分类。二是贝叶斯网络,该方法是一种基于概率推理的有向无环图模型,该方法克服了决策树的缺点——特征模型间的相互联系,但是因为本身结构的原因,会导致模型训练时间过长,以至于模型算法做不到收敛状态,导致建模失败。三是支持向量机(SVM),该方法是通过找寻数据序列之间的最大间距来对数据进行划分,重点在于区分数据,拥有着强大的鲁棒性和可移植性,但是其解决多分类问题困难,对参数和核函数选择敏感。四是神经网络,因其特殊的网络结构和较强的模型拟合能力发展十分迅速,其中RNN其因特殊的循环记忆结构可以充分的作用历史数据对预测结果的影响,所以对处理时间序列预测问题时有着十分显著的效果。长短期记忆网络(LSTM)是一种具有很强信息捕获能力的特殊的RNN模型,而门循环单元网络(GRU)是LSTM的一种更为简单高效的变体。由于在使用GRU模型预测过程中许多参数通常需要人为依靠经验确定,具有不确定性,致使模型精准度降低。而粒子群算法(PSO)可以调节较少的参数就可以实现快速收敛,通用性强,近年来被广泛应用于确定模型参数方面。针对燃油管道内压力温度等因素会在管道老化过程中诱发故障的问题,本文提出了一种基于PSO算法和GRU网络相结合的模型,提升燃油管道故障预测准确率。
发明内容
本发明要解决的技术问题是提供一种基于粒子群优化算法的门循环单元网络(PSO-GRU)混合模型的船舶燃油管道故障预测方法,使用该方法可以快速确定GRU网络模型最优参数,提高训练效率和预测精度,为燃油管道故障诊断提供理论预测,对于保证船舶安全、降低事故发生率发挥着重要作用。
为解决上述技术问题,本发明提出了一种基于PSO-GRU的船舶燃油管道故障预测方法,包括以下步骤:
步骤1:收集船舶燃油管道的供应泵压力、增压泵压力、加热器温度、冷却器温度作为特征参量;
步骤2:获取步骤1中的特征参量在历史时间维度船舶燃油管道中的样本数据,并对样本数据进行缺失检查,进行归一化处理,划分训练集数据和测试集数据;
步骤3:构建包括更新门、重置门的门循环单元网络模型GRU;
步骤4:用步骤2中获得的训练集数据对步骤3中构建的门循环单元网络预测模型GRU进行粒子群算法优化,确定门循环单元网络预测模型GRU的最优参数,以此为基础重新建立优化后的门循环单元网络模型PSO-GRU;
步骤5:将步骤2得到的训练集输入到步骤4中的PSO-GRU网络模型进行训练,得到的预测结果反归一化处理后和传感器正常工作数值范围进行对比,超出范围则判定为故障预警。
上述步骤2中利用线性插值法对数据缺失值进行补齐,其计算公式如下:
yk=yp+(yn-yp)(k-p)/(n-p) (1)
其中yk为要补齐的值,yp为yk前一个已知数据,yn为yk后一个已知数据。
进一步地,步骤2对样本数据进行归一化处理指将样本数据限制在一定范围内,本发明采用MinMaxScaler进行数据归一化和还原,其归一化的公式如下:
Figure BDA0003233563040000021
Figure BDA0003233563040000022
其中,xi为原始数据样本值;X.Min为样本序列中的最小值;X.Max为样本序列中的最大值;min为缩放的最小值;max为缩放的最大值;yi为标准化后的值。
上述步骤3中构建门循环单元网络预测模型GRU,主要是通过更新门和重置门的计算来完成整个神经单元的计算,计算过程如下:
zt=σ(Wz·[ht-1,xt]) (4)
rt=σ(Wr·[ht-1,xt]) (5)
h't=tanh(Wh'·[rt*ht-1,xt]) (6)
ht=(1-zt)*ht-1+zt*h't (7)
其中xt表示当前时刻的输入;ht-1表示上一个神经元的输出,rt为重置门,zt为更新门,h’t为隐藏层状态,ht为当前输出。·代表矩阵乘积;*代表点积;σ代表sigmoid函数,tanh为双曲正切函数,Wz为更新门权重矩阵,Wr为重置门权重矩阵。
上述步骤4中使用粒子群优化算法PSO优化门循环网络模型参数的过程如下:
步骤4.1:以隐藏层数、时间序列大小、隐藏层节点数作为PSO算法的优化超参数,根据经验给予超参数各自的范围与粒子速度等条件,设置迭代的最大次数以及粒子最大飞行速度等。
步骤4.2:确定一个让粒子达到标准的函数,通过每个粒子的运动寻找该粒子的最优解,再通过每一个粒子的最优解推导出本次所有粒子的最优解,叫做本次全局最优解。与历史全局最优比较,进行更新。
步骤4.3:通过粒子不断地更新自己的速度与位置来寻求最优解,具体如下式所示:
vid k+1=ωVid k+z1r1(Pid k-Xid k)+z2r2(Pgd k-Xid k) (8)
Figure BDA0003233563040000031
其中:k表示当前迭代次数;Vid k、Xid k、Pid k、Pgd k分别表示粒子的速度、位置、本次单体最优解和本次全局最优解;z1与z2为环境加速因子∈[1,4];r1与r2为0到1之间的随机数;ω为惯性因子∈[0.5,1.5],其值非负,其值较大时全局寻优能力强,局部寻优能力较弱;其值较小时全局寻优能力弱,局部寻优能力强,所以可以通过调整ω的大小来调节寻优效果。
步骤4.4:然后将处理好的测试数据代入已经训练好的模型进行预测,以模型在测试数据集上的绝对平均百分比误差(MAPE)作为粒子适应度值,其公式如下:
Figure BDA0003233563040000032
其中m为实验预测次数,yi为真实值,
Figure BDA0003233563040000033
为模型预测值。
步骤4.5:当迭代次数达到设置的最大迭代次数,或者全局最优解到达设置的最小边界时算法终止。获得最优隐藏层数、时间序列大小、隐藏层节点数,将最优参数带入GRU单元网络模型进行训练,得到经过粒子群优化后的门循环单元网络模型PSO-GRU。
上述步骤5得到的预测结果使用均方根误差(RMSE)和平均绝对误差(MAE)作为评估指标对模型预测效果进行评价;RMSE越低,表示模型越稳定;MAE越低,表示模型精准度越高,具体公式定义如下:
Figure BDA0003233563040000041
Figure BDA0003233563040000042
其中m为实验预测次数,yi为真实值,
Figure BDA0003233563040000043
为模型预测值。
本发明提供的预测方法存在以下有益效果:本发明首先以船舶燃油管道供应泵压力值作为门循环单元网络模型的输入;然后通过粒子群优化算法迭代计算出GRU网络模型的最优参数(隐藏层数、时间序列大小、隐藏层节点数);最后建立了PSO-GRU网络预测模型,通过预测数据预判燃油管道故障发生;经过和常见的几种预测模型对比后得出:该优化后的模型解决了依据经验选取参数而导致的模型拟合能力不够、预测精度低的问题,可准确预测燃油管道内压力变化趋势,为后续燃油管道的故障诊断、状态评估提供理论指导,对保障船舶安全运行、降低管道维护成本起着重要作用。
附图说明
图1为本发明基于PSO-GRU网络模型对船舶燃油管道故障预测方法的整体流程图;
图2为本发明中数据预处理方法流程图;
图3为本发明中门循环单元网络模型GRU结构图;
图4为本发明中PSO-GRU神经网络模型结构图;
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明的一部分实施例,而不是全部的实施例;基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
参照图1-4,本发明提供了一种基于PSO-GRU混合模型的船舶燃油管道故障预测方法,该方法主要包括以下步骤:
步骤1:收集船舶燃油管道的供应泵压力、增压泵压力、加热器温度、冷却器温度作为特征参量,具体为:
步骤1.1:从航海船舶终端设备传感器采集燃油管道的供应泵压力、增压泵压力、加热器温度、冷却器温度,采集信息频率为1次每秒;
步骤1.2:在约定好的通信规则下将时间信息、船只代码、设备代码以及步骤1.1传感器数据上传至服务器;
步骤1.3:将步骤1.2全部原始数据保存到数据库中,并发布到RabbitMQ消息队列中;
步骤2:获取步骤1中的特征参量在历史时间维度船舶燃油管道中的样本数据,并对样本数据进行缺失检查,进行归一化处理,划分训练集数据和测试集数据,具体为:
步骤2.1:燃油管道供应泵压力值可能由于传感器故障或者传输失误产生数据缺失,为了得到干净完整的数据,提高故障预测准确性,本发明利用线性插值法对数据缺失值进行补齐,其计算公式如下:
yk=yp+(yn-yp)(k-p)/(n-p) (1)
其中yk为要补齐的值,yp为yk前一个已知数据,yn为yk后一个已知数据。
步骤2.2:在将数据输入网络模型之前,要进行归一化处理,归一化使得梯度始终朝着最小值的方向前进,在一定程度上提高模型精度,提升收敛速度。本发明采用MinMaxScaler进行数据归一化和还原,具体公式如下所示:
Figure BDA0003233563040000051
Figure BDA0003233563040000052
其中,xi为原始数据样本值;X.Min为样本序列中的最小值;X.Max为样本序列中的最大值;min为缩放的最小值;max为缩放的最大值;yi为标准化后的值。
步骤2.3:将经过步骤2.1、2.2处理的数据90%作为训练集数据,10%作为测试集数据;
步骤3:构建包括更新门、重置门的门循环单元网络模型GRU,具体为:
步骤3.1:构建门循环单元网络模型GRU,主要是通过更新门和重置门的计算来完成整个神经单元的计算,计算过程如下:
zt=σ(Wz·[ht-1,xt]) (4)
rt=σ(Wr·[ht-1,xt]) (5)
h't=tanh(Wh'·[rt*ht-1,xt]) (6)
ht=(1-zt)*ht-1+zt*h't (7)
其中xt表示当前时刻的输入;ht-1表示上一个神经元的输出,rt为重置门,zt为更新门,h’t为隐藏层状态,ht为当前输出。·代表矩阵乘积;*代表点积;σ代表sigmoid函数,tanh为双曲正切函数,Wz为更新门权重矩阵,Wr为重置门权重矩阵。
步骤3.2:选择Leaky ReLUs激活函数作为模型的激活函数,其计算公式如下式所示:
Figure BDA0003233563040000061
其中,a通常是通过过往经验而人为给予的值。在取值区间上的值不为0,且斜率、梯度都不为0,这使得在梯度下降过程中,这一部分也将得到更新,保留了网络模型复杂度,同时因为斜率一直存在也可以减缓梯度消失的现象。
步骤3.3:本发明在设计网络模型时选取均方误差法作为损失函数,其计算公式如下式所示:
Figure BDA0003233563040000062
其中zi代表实际值,z代表预测值。
步骤3.4:本发明在模型训练时应用的参数优化方法为梯度下降法,具体公式如下式所示:
Figure BDA0003233563040000063
其中η代表学习速率,即每次下降的步长。
步骤4:用步骤2中获得的训练集数据对步骤3中构建的门循环单元网络预测模型GRU进行粒子群算法优化,确定门循环单元网络预测模型GRU的最优参数,以此为基础重新建立优化后的门循环单元网络预测模型PSO-GRU,具体为:
步骤4.1:以隐藏层数、时间序列大小、隐藏层节点数作为PSO算法的优化超参数,根据经验给予超参数各自的范围与粒子速度等条件,设置迭代的最大次数为1000、以及粒子最大飞行速度等。
步骤4.2:确定一个让粒子达到标准的函数,通过每个粒子的运动寻找该粒子的最优解,再通过每一个粒子的最优解推导出本次所有粒子的最优解,叫做本次全局最优解。与历史全局最优比较,进行更新。
步骤4.3:通过粒子不断地更新自己的速度与位置来寻求最优解,具体如下式所示:
vid k+1=ωVid k+z1r1(Pid k-Xid k)+z2r2(Pgd k-Xid k) (11)
Figure BDA0003233563040000071
其中:k表示当前迭代次数;Vid k、Xid k、Pid k、Pgd k分别表示粒子的速度、位置、本次单体最优解和本次全局最优解;z1与z2为环境加速因子∈[1,4];r1与r2为0到1之间的随机数;ω为惯性因子∈[0.5,1.5],其值非负,其值较大时全局寻优能力强,局部寻优能力较弱;其值较小时全局寻优能力弱,局部寻优能力强,所以可以通过调整ω的大小来调节寻优效果。
步骤4.4:然后将处理好的测试数据代入已经训练好的模型进行预测,以模型在测试数据集上的绝对平均百分比误差(MAPE)作为粒子适应度值,其公式如下:
Figure BDA0003233563040000072
其中m为实验预测次数,yi为真实值,
Figure BDA0003233563040000073
为模型预测值。
步骤4.5:当迭代次数达到设置的最大迭代次数,或者全局最优解到达设置的最小边界时算法终止。获得最优隐藏层数为2、时间序列大小为20、隐藏层节点数为30,将最优参数带入GRU单元网络模型进行训练,得到经过粒子群优化后的门循环单元网络模型PSO-GRU。
步骤5:将步骤2得到的训练集输入到步骤4中的PSO-GRU网络模型进行训练,得到的预测结果反归一化处理后和传感器正常工作数值范围进行对比,超出范围则判定为故障预警,具体为:
步骤5.1:本发明得到的预测结果使用均方根误差(RMSE)和平均绝对误差(MAE)作为评估指标对模型预测效果进行评价。RMSE越低,模型越稳定;MAE越低,模型精准度越高,具体公式定义如下:
Figure BDA0003233563040000074
Figure BDA0003233563040000075
其中m为实验预测次数,yi为真实值,
Figure BDA0003233563040000076
为模型预测值。
步骤5.2:将步骤2中的训练集带入普通LSTM模型、GRU模型、PSO-GRU模型进行训练,预测结果同样根据步骤5.1的评估指标进行计算,与本发明基于PSO-GRU模型预测结果进行比较,试验结果表明PSO-GRU预测模型的均方根误差(RMSE)为0.0055,平均绝对误差(MAE)为0.0689,与LSTM模型、GRU模型和PSO-LSTM模型相比分别降低了39.6%、50.4%,39.5%、49.6%和9.8%、13.4%,对船舶燃油管道具有良好预测效果。
以上对本发明所提出的一种基于粒子群算法和门循环单元网络混合模型(PSO-GRU)的船舶燃油管道故障预测方法,进行了详细介绍,本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。综上所述,本说明书内容不应理解为对本发明的限制。

Claims (5)

1.一种基于PSO-GRU的船舶燃油管道故障预测方法,其特征在于,包括以下步骤:
步骤1:收集船舶燃油管道的供应泵压力、增压泵压力、加热器温度、冷却器温度作为特征参量;
步骤2:获取步骤1中的特征参量在历史时间维度船舶燃油管道中的样本数据,并对样本数据进行缺失检查,进行归一化处理,划分训练集数据和测试集数据;
步骤3:构建包括更新门、重置门的门循环单元网络模型GRU;
步骤4:用步骤2中获得的训练集数据对步骤3中构建的门循环单元网络模型GRU进行粒子群算法优化,确定门循环单元网络模型GRU的最优参数,以此为基础重新建立优化后的门循环单元网络预测模型PSO-GRU;
步骤5:将步骤2得到的训练集输入到步骤4中的PSO-GRU网络模型进行训练,得到的预测结果反归一化处理后和传感器正常工作数值范围进行对比,超出范围则判定为故障预警。
2.根据权利要求1所述的一种基于PSO-GRU的船舶燃油管道故障预测方法,其特征在于,步骤2中所述的对数据进行缺失检查,采用线性插值方法对缺失值进行补齐,其计算公式如下:
yk=yp+(yn-yp)(k-p)/(n-p) (1)
其中yk为要补齐的值,yp为yk前一个已知数据,yn为yk后一个已知数据;
步骤2中对样本数据采用MinMaxScaler进行数据归一化和还原,其归一化的公式如下:
Figure FDA0003233563030000011
Figure FDA0003233563030000012
其中,xi为原始数据样本值;X.Min为样本序列中的最小值;X.Max为样本序列中的最大值;min为缩放的最小值;max为缩放的最大值;yi为标准化后的值。
3.根据权利要求1所述的一种基于PSO-GRU的船舶燃油管道故障预测方法,其特征在于,步骤3中所述的构建门循环单元网络模型GRU,主要过程如下:
步骤3.1:构建门循环单元网络模型GRU,主要是通过更新门和重置门的计算来完成整个神经单元的计算,计算过程如下:
zt=σ(Wz·[ht-1,xt]) (4)
rt=σ(Wr·[ht-1,xt]) (5)
h't=tanh(Wh'·[rt*ht-1,xt]) (6)
ht=(1-zt)*ht-1+zt*h't (7)
其中xt表示当前时刻的输入;ht-1表示上一个神经元的输出,rt为重置门,zt为更新门,h’t为隐藏层状态,ht为当前输出,·代表矩阵乘积;*代表点积;σ代表sigmoid函数,tanh为双曲正切函数,Wz为更新门权重矩阵,Wr为重置门权重矩阵;
步骤3.2:选择Leaky ReLUs激活函数作为模型的激活函数,其计算公式如下式所示:
Figure FDA0003233563030000021
其中,a通常是通过过往经验而人为给予的值,在取值区间上的值不为0,且斜率、梯度都不为0;
步骤3.3:在设计网络模型时选取均方误差法作为损失函数,其计算公式如下式所示:
Figure FDA0003233563030000022
其中zi代表实际值,z代表预测值;
步骤3.4:在模型训练时应用的参数优化方法为梯度下降法,具体公式如下式所示:
Figure FDA0003233563030000023
其中η代表学习速率,即每次下降的步长。
4.根据权利要求1所述的一种基于PSO-GRU的船舶燃油管道故障预测方法,其特征在于,步骤4中所述的使用粒子群优化算法PSO优化门循环单元网络模型参数的过程如下:
步骤4.1:以隐藏层数、时间序列大小、隐藏层节点数作为PSO算法的优化超参数,根据经验给予超参数各自的范围与粒子速度,设置迭代的最大次数以及粒子最大飞行速度;
步骤4.2:确定一个让粒子达到标准的函数,通过每个粒子的运动寻找该粒子的最优解,再通过每一个粒子的最优解推导出本次所有粒子的最优解,叫做本次全局最优解;与历史全局最优比较,进行更新;
步骤4.3:通过粒子不断地更新自己的速度与位置来寻求最优解,具体如下式所示:
vid k+1=ωVid k+z1r1(Pid k-Xid k)+z2r2(Pgd k-Xid k) (11)
Figure FDA0003233563030000031
其中:k表示当前迭代次数;Vid k、Xid k、Pid k、Pgd k分别表示粒子的速度、位置、本次单体最优解和本次全局最优解;z1与z2为环境加速因子∈[1,4];r1与r2为0到1之间的随机数;ω为惯性因子∈[0.5,1.5],其值非负,其值较大时全局寻优能力强,局部寻优能力较弱;其值较小时全局寻优能力弱,局部寻优能力强,通过调整ω的大小来调节寻优效果;
步骤4.4:然后将处理好的测试数据代入已经训练好的模型进行预测,以模型在测试数据集上的绝对平均百分比误差(MAPE)作为粒子适应度值,其公式如下:
Figure FDA0003233563030000032
其中m为实验预测次数,yi为真实值,
Figure FDA0003233563030000033
为模型预测值;
步骤4.5:当迭代次数达到设置的最大迭代次数,或者全局最优解到达设置的最小边界时算法终止;获得最优隐藏层数、时间序列大小、隐藏层节点数,将最优参数带入门循环单元网络模型GRU进行训练,得到经过粒子群优化后的门循环单元网络模型PSO-GRU。
5.根据权利要求1所述的一种基于PSO-GRU的船舶燃油管道故障预测方法,其特征在于,步骤5中得到的预测结果使用均方根误差(RMSE)和平均绝对误差(MAE)作为评估指标对模型预测效果进行评价,具体公式定义如下:
Figure FDA0003233563030000041
Figure FDA0003233563030000042
其中m为实验预测次数,yi为真实值,
Figure FDA0003233563030000043
为模型预测值;RMSE越低,表示模型越稳定,MAE越低,表示模型精准度越高。
CN202110995049.5A 2021-08-27 2021-08-27 一种基于pso-gru的船舶燃油管道故障预测方法 Active CN113705094B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110995049.5A CN113705094B (zh) 2021-08-27 2021-08-27 一种基于pso-gru的船舶燃油管道故障预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110995049.5A CN113705094B (zh) 2021-08-27 2021-08-27 一种基于pso-gru的船舶燃油管道故障预测方法

Publications (2)

Publication Number Publication Date
CN113705094A true CN113705094A (zh) 2021-11-26
CN113705094B CN113705094B (zh) 2024-01-23

Family

ID=78655920

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110995049.5A Active CN113705094B (zh) 2021-08-27 2021-08-27 一种基于pso-gru的船舶燃油管道故障预测方法

Country Status (1)

Country Link
CN (1) CN113705094B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114580250A (zh) * 2022-01-26 2022-06-03 江苏康缘药业股份有限公司 基于量子粒子群算法的gru网络的能耗预测方法、装置
CN115860718A (zh) * 2022-12-02 2023-03-28 深圳市黑金工业制造有限公司 一种基于大数据的触摸一体机综合维护管理系统及方法
CN116242368A (zh) * 2023-05-11 2023-06-09 北京航空航天大学 一种基于偏振光场信息预测的偏振/惯性容错导航方法
CN116384257A (zh) * 2023-05-29 2023-07-04 浙江大学 一种空分整装冷箱装配误差预测与公差优化方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110263907A (zh) * 2019-07-01 2019-09-20 上海海事大学 基于改进ga-pso-bp的船舶短路故障诊断方法
CN110320886A (zh) * 2019-06-04 2019-10-11 武汉理工大学 一种船舶辅助系统船岸一体化监控和故障诊断的方法
CN111426816A (zh) * 2020-04-10 2020-07-17 昆明理工大学 一种基于pso-lstm的变压器油中溶解气体浓度预测方法
CN111738512A (zh) * 2020-06-22 2020-10-02 昆明理工大学 一种基于cnn-ipso-gru混合模型的短期电力负荷预测方法
CN112348271A (zh) * 2020-11-12 2021-02-09 华北电力大学 基于vmd-ipso-gru的短期光伏功率预测方法
US20210117796A1 (en) * 2020-12-01 2021-04-22 Harbin Engineering University Ship Motion Prediction Method Based on Long Short-Term Memory Network and Gaussian Process Regression

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110320886A (zh) * 2019-06-04 2019-10-11 武汉理工大学 一种船舶辅助系统船岸一体化监控和故障诊断的方法
CN110263907A (zh) * 2019-07-01 2019-09-20 上海海事大学 基于改进ga-pso-bp的船舶短路故障诊断方法
CN111426816A (zh) * 2020-04-10 2020-07-17 昆明理工大学 一种基于pso-lstm的变压器油中溶解气体浓度预测方法
CN111738512A (zh) * 2020-06-22 2020-10-02 昆明理工大学 一种基于cnn-ipso-gru混合模型的短期电力负荷预测方法
CN112348271A (zh) * 2020-11-12 2021-02-09 华北电力大学 基于vmd-ipso-gru的短期光伏功率预测方法
US20210117796A1 (en) * 2020-12-01 2021-04-22 Harbin Engineering University Ship Motion Prediction Method Based on Long Short-Term Memory Network and Gaussian Process Regression

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
SHAO-PEI JI等: "GRU-corr Neural Network Optimized by Improved PSO Algorithm for Time Series Prediction", 《INTERNATIONAL JOURNAL ON ARTIFICIAL INTELLIGENCE TOOLS》, pages 302 - 306 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114580250A (zh) * 2022-01-26 2022-06-03 江苏康缘药业股份有限公司 基于量子粒子群算法的gru网络的能耗预测方法、装置
CN115860718A (zh) * 2022-12-02 2023-03-28 深圳市黑金工业制造有限公司 一种基于大数据的触摸一体机综合维护管理系统及方法
CN116242368A (zh) * 2023-05-11 2023-06-09 北京航空航天大学 一种基于偏振光场信息预测的偏振/惯性容错导航方法
CN116242368B (zh) * 2023-05-11 2023-07-25 北京航空航天大学 一种基于偏振光场信息预测的偏振/惯性容错导航方法
CN116384257A (zh) * 2023-05-29 2023-07-04 浙江大学 一种空分整装冷箱装配误差预测与公差优化方法
CN116384257B (zh) * 2023-05-29 2023-09-29 浙江大学 一种空分整装冷箱装配误差预测与公差优化方法

Also Published As

Publication number Publication date
CN113705094B (zh) 2024-01-23

Similar Documents

Publication Publication Date Title
CN113705094B (zh) 一种基于pso-gru的船舶燃油管道故障预测方法
CN110321603B (zh) 一种用于航空发动机气路故障诊断的深度计算模型
CN111815033A (zh) 一种基于rcnn和气象时序特征的海上风电功率预测方法
Dao et al. Fault diagnosis of hydro-turbine via the incorporation of bayesian algorithm optimized CNN-LSTM neural network
CN111242351A (zh) 基于自编码器和gru神经网络的热带气旋轨迹预测方法
CN107944648A (zh) 一种大型船舶航速油耗率精准预测方法
CN113895272A (zh) 基于深度学习的电动汽车交流充电状态监测和故障预警方法
CN113947182B (zh) 基于双阶段堆叠图卷积网络的交通流预测模型构建方法
CN109471049B (zh) 一种基于改进堆叠自编码器的卫星电源系统异常检测方法
Yao et al. Wave height forecast method with multi-step training set extension LSTM neural network
CN115877483A (zh) 一种基于随机森林和gru的台风路径预报方法
CN116248392A (zh) 一种基于多头注意力机制的网络恶意流量检测系统及方法
CN112241802A (zh) 一种风电功率的区间预测方法
CN114492217A (zh) 一种基于混合时间序列感知机模型的台风浪高预测方法
Pérez-Pérez et al. Fault detection and isolation in wind turbines based on neuro-fuzzy qLPV zonotopic observers
CN113673774A (zh) 基于自编码器和时序卷积网络的航空发动机余寿预测方法
CN112465136A (zh) 一种水电机组劣化预测方法和系统
WO2014173271A1 (zh) 数字化人机界面监视单元数量优化方法及系统
CN116663414A (zh) 一种用于电力变压器的故障诊断方法及系统
Yao et al. Extreme motion prediction and early-warning assessment of semisubmersible platform based on deep learning method
CN115051834A (zh) 一种基于STSA-transformer算法的新型电力系统APT攻击检测方法
CN115208619A (zh) 基于STSA-transformer算法的新型电力系统APT攻击检测方法
CN115081483A (zh) 基于特征选择和gwo-bp的水轮发电机转子故障诊断方法
CN113761795A (zh) 一种飞机发动机故障检测方法和系统
CN113837443A (zh) 基于深度BiLSTM的变电站线路负载预测方法

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