基于神经网络的集装箱船港口作业时间预测模型

第38卷第2期   计算机应用与软件
Vol 38No.22021年2月 
ComputerApplicationsandSoftware
Feb.2021
基于神经网络的集装箱船港口作业时间预测模型
韩宗垒 徐 斌 陈 佳
(大连海事大学辽宁省物流航运管理系统工程重点实验室 辽宁大连116026)
收稿日期:2019-07-05。韩宗垒,硕士生,主研领域:信息管理与信息系统。徐斌,副教授。陈佳,教授。
摘 要  集装箱船港口作业时间是制作泊位计划的一个重要依据,而集装箱船港口作业时间获取的主要来源是预测。传统预测方法是用装卸集装箱量除以岸桥装卸效率,预测精度较低,且受多种因素的影响,具有复杂的非线性特点。而神经网络在解决复杂的非线性问题方面具有很强的建模能力,所以选
取神经网络建立集装箱船港口作业时间预测模型。通过真实数据对预测模型进行训练学习,用测试数据集对模型进行验证,并且与传统预测方法相对比,结果表明了该预测模型在某集装箱港口预测应用的有效性。关键词  水路运输 集装箱船港口作业时间 神经网络 预测 模型
中图分类号 TP391    文献标志码 A    DOI:10.3969/j.issn.1000 386x.2021.02.014
CONTAINERSHIPPORTOPERATIONTIMEPREDICTIONMODELBASED
ONNEURALNETWORK
HanZonglei XuBin ChenJia
(LiaoningProvincialKeyLaboratoryofLogisticsandShippingManagementSystemEngineering,
DalianMaritimeUniversity,Dalian116026,Liaoning,China)
Abstract  Containershipportoperationtimeisanimportantbasisfortheproductionofberthplans,andthemainsourceofcontainershipportoperationtimeisprediction.Thetraditionalpredictionmethodistheloadingandunloadingcontainervolumedividedbytheshorebridgeloadingandunloadingefficiency.Thepredictionaccuracyofthismethodislow,anditisaffectedbymanyfactorsandhascomplexnonlinearcharacteristics.Theneuralnetworkhasstrongmodelingabilityinsolvingcomplexnonlinearproblems.Therefore,theneuralnetworkwasselectedtoestablishacontainershipportoperationtimepredictionmodel.Thepredictionmodelwastrainedbyrealdata,andthemodelwasvalidatedbythetestdataset.Comparedwiththetraditionalpredictionmethod,theeffectivenessofthepredictionmodelinacontainerportpredictionapplicationisdemonstrated.
Keywords  Waterwaytransportation Containershipportoperationtime Neuralnetworks Prediction Model
0 引 言
泊位是港口的重要资源,合理的泊位计划能有效提高港口生产效率,减少船舶在港时间,从而提高港口
的竞争力。Imai等[1]
?用连续区位空间的方法对连
续型泊位计划进行研究,建立了船舶等待时间和作业时间最小的数学模型并采用拉格朗日松驰系数算法进
行求解。张煜等[
2]
考虑在泊位计划中岸桥的分配影响集装箱船的装卸作业效率,依据规则建立确定船舶集装箱装卸作业时间和分配岸边起重机的算法。秦进
等[3]
提出基于时间窗约束的离散型泊位计划模型,建
立了模拟退火算法模型进行优化求解。韩晓龙等[4]
港口船舶的目标函数中加入了船舶作业时间,对船舶的服务时间和作业时间提出了时间窗约束,建立了以最小化卸船完工时间为优化目标的混合整数规划模型
和约束规划模型。曾庆成等[5]用船舶装卸集装箱量除
以岸桥装卸效率来预测集装箱船港口作业时间。
综上所述,大量的研究表明集装箱船港口作业时间对于制作科学高效的泊位计划是非常重要的,研究中用的是传统装卸集装箱量除以岸桥装卸效率的预测方法,预测精度较低,因此寻一种科学的方法来预测
第2期   韩宗垒,等:基于神经网络的集装箱船港口作业时间预测模型79
集装箱船港口作业时间是非常有意义的。本文基于神经网络针对集装箱船港口作业时间预测的问题建立模型,通过与传统预测方法的对比验证了模型的可行性。
1 问题描述
在集装箱船实际到达港口之前,港口需要根据船期表制作泊位计划确定集装箱船靠泊的时间和位置。集装箱船港口作业时间(船舶开始装卸第一个集装箱到完成装卸最后一个集装箱的这段时间)是制作泊位计划的主要依据,而集装箱船港口作业时间的主要获取方法是预测,所以预测出精确的集装箱船港口作业时间可以提高泊位计划的效率。
传统的集装箱船港口作业时间预测方法是船舶待装卸的集装箱量除以岸桥的装卸效率。这种预测方法不灵活并且预测精度较低,集装箱船港口作业时间受多种因素的影响[6],如船舶类型、岸桥数、装卸集装箱量、天气等,并且存在复杂的非线性关系。考虑到集装箱船港口作业时间与影响因素之间复杂的非线性关系,本文选取了比较适用的BP神经网络建立集装箱船港口作业时间预测模型。模型的目标是预测出更加精确的集装箱船港口作业时间,从而保证制作的泊位计划更加科学高效。
2 模型建立
BP(BackPropagation)神经网络是一种多层前馈神经网络,由输入层、隐含层、输出层组成,可以以任意的精确度逼近任意一个连续的函数,所以经常被用于非线性建模、函数逼近和模式分类等方面。BP神经网络的主要特点:信号是前向传播的,误差是反向传播的。
2.1 网络的层数
1989年RobertHecht Nielson证明了对于任何一个闭区间内的函数,都可以用有一个隐含层的BP神经网络来逼近,所以一个三层(含一个隐含层)的BP神经网络可以完成任意的n维向m维的映射。在多数的实际应用中一般都取三层的BP神经网络来解决问题,所以本模型选取了三层的神经网络。
2.2 输入和输出层神经元个数
BP神经网络的输入和输出层的神经元个数完全根据使用者的要求进行设计。影响船舶港口作业时间的因素有很多,通过分析确定了船舶类型、岸桥数量、卸20尺箱量、装20尺箱量、卸40尺箱量、装40尺箱量、卸特种箱量、装特种箱量和天气作为模型的输入。本模型根据输入样本的维度,将输入层设置为9个神经元。集装箱船港口作业时间预测模型最后输出的是集装箱船港口作业时间,所以输出层神经元的个数为1。
2.3 隐含层神经元个数
隐含层对于整个神经网络的精度起着至关重要的作用,隐含层神经元个数的选取与输入、输出层神经元个数都有直接的关系。神经网络精度的提高可以通过增加隐含层神经元的个数来实现,隐含层神经元的个数较少时,神经网络就不能较好地学习,导致模型不能较好地拟合数据,训练精度和预测精度都不高,出现欠拟合的情况;隐含层神经元个数较多时,模型结构可能过于复杂,导致过度的拟合训练数据,而对测试数据的预测能力较差,出现过拟合的情况,所以选取准确的隐含层神经元个数是很重要的。隐含层神经元个数选取的方法很多,结合开发的神经网络生成器,本模型采用式(1)选取隐含层神经元个数[7]。
p=m+n+a(1)式中:p为隐含层神经元个数;n为输入层神经元个数;m为输出层神经元个数;a为1~10之间的常数。经过试算最终确定隐含层神经元个数为10。
2.4 激活函数
BP神经网络的激活函数有多种。其中Sigmoid函数对于网络不同输入,将其输出范围控制在(0,1)之间,公式如下:
f(x)=
1+e-x
(2)Tanh(双曲正切)激活函数对于不同范围的输入,将其输出值范围控制在(-1,1)之间,公式如下:
f(x)=
ex-e-x
ex+e-x
(3)线性激活函数Purelin的输入与输出值可取任意值,公式如下:
f(x)=x(4)本模型选取Sigmoid作为隐含层的激活函数,Purelin作为输出层的激活函数。
2.5 LM BP算法
标准BP算法是根据梯度下降法来调整权值的:
Δw=-ηg(5)式中:Δw为权值阈值更新量;η为学习速率;g为梯度。
80
    计算机应用与软件2021年权值沿着与误差相反的方向移动,使得误差函数
减小,缺点是神经网络收敛较慢,且学习速率不容易被
确定。LM(LevenbergMarquardt)算法是一种利用标准
的数值优化技术的快速算法,是梯度下降法和高斯 牛
顿法的结合,既有高斯 牛顿法的局部收敛性,又有梯
度下降法的全局特性,具有收敛速度快、鲁棒性好的特
点。下面对LM算法做简要阐述:
设网络的误差函数为:
E(w)=1
∑p
i=1
e2
(w)(6)
式中:w是权值和阈值组成的向量;p是样本数;e2
(w)是误差的平方。
用wk表示在第k次迭代的权值和阈值,迭代完成后的权值和阈值组成的向量为wk+1,Δw是权值和阈值的改变量,则有:
wk+1=wk+Δw(7)牛顿法是通过最小二乘法求解误差函数E(w):Δw=-[ E2(w)]-1 E(w)(8)式中: E2(w)是误差E(w)的Hessian矩阵, E(w)表示梯度,对Hessian矩阵进行近似计算,可以表明:
E(w)=JT(w)e(w)(9) E2(w)=JT(w)e(w)+S(w)(10)式中:J(w)是e(w)的Jacobian矩阵;S(w)是误差矩阵。
J(w)=
e1(w)
w1
e1(w)
w2
e1(w)神经网络预测
wn e2(w)
w1
e2(w)
w2
e2(w)
wn
eN(w)
w1
eN(w)
w2
eN(w)
在靠近极值点时S(w)≈0,牛顿法可以修正为高斯 牛顿法,经过改进得到修正权值阈值的公式:
Δw=-[JT(w)J(w)]-1J(w)e(w)(11)LM算法将高斯 牛顿法经过改进得到修正权值阈值的公式:
Δw=[JT(w)J(w)+μI]-1J(w)e(w)(12)式中:I为单位矩阵;μ为大于0的常数。
系数μ的值很小时,LM算法就近似等于高斯 牛顿法,当μ的值很大时,就近似等于梯度下降法。每迭代成功一次μ就会除以比例系数β(β>1),这样在接近目标误差的时候就基本与高斯 牛顿法相等,计算速度快,精确度也高,否则μ乘比例系数β,LM算法利用近似二阶导数信息,比梯度下降法快得多。在实际应用中μ是一个试探性的参数,对于一个给定值,如果求得Δw能使E(w)降低,则μ降低,反之μ增加。2.6 相关参数
(1)权值和阈值。选取处于(-1,1)之间的随机数作为权值和阈值的初始值。
(2)学习速率。神经网络权值每次的变化量都取决于学习速率的大小,如果学习速率选取较大,系统可能因此而动荡不稳定;学习速率选取较小则收敛速度慢,训练时间长,网络误差值与误差最小值更趋近的目标无法保障。实际应用中常选取较小学习速率给系统提供稳定性保障,所以学习速率的选取区间是[0.01,0.9]。本模型选取的学习速率为0.01。
(3)其他。最大训练次数为1000次,训练要求精度为0.0001,极小值认定次数为50,μ的初始值为0.00001,比例系数β为10。
3 算例分析
某港口近年集装箱船港口作业数据中包含船舶类型、分配岸桥数、装卸集装箱量、天气、计划作业时间、实际作业时间等。集装箱船的第一代和第六代相当少,这里只考虑第二至第五代的船型。因为原始数据中船舶类型和天气都是文字形式,所以需要对船舶类型和天气情况进行编码数字化。船舶类型和天气情况编码后的结果如表1和表2所示。
表1 编码后的船舶类型
类型
船长/m船宽/m
吃水
深度/m
载箱量
/TEU
编码
二代
175~
225
25~
30
9.5~
10.5
1000~
2000
三代
240~
275
30~
32
10.5~
12.0
2000~
3000
四代
275~
295
32~
35
11.5~
12.5
3000~
4500
五代
280~
300
32.2~
39.4
11.5~
13.5
4500~
6000
第2期   韩宗垒,等:基于神经网络的集装箱船港口作业时间预测模型81
 表2 编码后的天气
风速天气编码
1~6级晴天、多云、阴天1小雨2
中雨3
小雪4
中雪5
6~8级晴天、多云、阴天6小雨7
中雨8
小雪9
中雪10
随机选取1000条集装箱船港口作业信息数据作为模型的训练集和测试集数据,由于数据量较大,这里只展示部分编码后的数据,如表3所示。
表3 编码后的船舶信息
船舶类型岸
20
20
40
40
/h
2324129888122453399322381759599011541242229741296147732104213311898142854721022871852191899328913.533273126229101267969.533173159907811214048.3
1)传统的集装箱船装卸集装箱量除以岸桥装卸效率的预测方法可以根据式(13)计算:
h=
v×n
(13)
式中:h为集装箱船港口作业时间;Y为装卸集装箱总量;v为单个岸桥装卸效率;n为分配岸桥数。
港口单个岸桥装卸效率为35箱/h,用传统预测方法计算集装箱船港口作业时间结果如表4所示。
表4 传统方法预测作业时间
装卸集装箱量分配岸桥数预计作业时间/h82737.9
722210.3
61328.8
62328.9
801211.4
83437.9
75237.2
2)用神经网络模型预测集装箱船港口作业时间:为了消除各参数由于单位等的影响,并且样本不一定包含极大和极小值,所以对数据用式(14)做规范化处理,使规范化后的数据范围为[-1,1]。
y=
(y
max
-y
min
)×(x-x
min
max
-x
min
+y
min
(14)
式中:y为规范化后的数值;y
max
=1;y
min
=-1;x
max
每一属性中的最大值;x
min
为每一属性中的最小值;x为需要规范化的数据值。
模型训练学习过程如下:选取样本中的3条数据如表5所示,规范化后的样本数据如表6所示,模型的
输入为x
,x
,…,x
,输出为y。
表5 样本数据
20
20
40
40
/h
32621372692146118310
22176101859399128712
431141011141072514226.5
表6 规范化后的样本数据
0-1-11-1-11-1-0.60.27-1-11-10.34-0.870.22-0.1711
11-0.09-111-11-1-1(1)设定好网络的期望误差值ε=0.0001,系数β=10,μ=0.1,学习速率η=0.01,以及权值和阈值的向量:
82
    计算机应用与软件2021年
wk=0.84,-0.86,0.38,0.51,-0.18,-0.60,0.67,0.86,0.56,0.75,0.87,-0.86,0.47,-0.21,0.34,-0.46,-0.45,-0.75,-0.83,0.99,0.78,-0.90,0.60,0.35,0.28,0.70,0.85,0.82,
-0.03,0.87,-0.79,0.59,-0.67,0.39,-0.50,0.37,0.82,-0.03,0.87,-0.79,0.59,-0.67,0.39,-0.50,0.37,0.25,0.59,0.36,0.69,0.61,-0.75,0.77,0.8
5,-0.29,-0.86,-0.76,0.54,0.69,-0.05,-0.03,0.06,0.36,0.85,
-0.63,-0.24,0.69,0.56,-0.12,0.29,-0.10,0.80,0.64,0.03,0.05,-0.20,-0.41,0.49,
-0.34,-0.85,0.59,0.32,0.72,0.67,0.31,0.97,0.43,0.16,0.75,0.23,0.72,0.87,
0.86,-0.61,-0.87,0.50,-0.51,0.47,-0.03,0.30,0.68,0.68,0.72,-0.93,0.20,-0.56,0.66,-0.17,0.82,-0.64,0.
58
(2)计算神经网络的输出y^=0.4377。(3)计算Jacobian矩阵:
J(w)=
-0.00528,0.00540,-0.00239,-0.00320,0.00113,0.00377,0.00377,-0.00421,
-0.00540,-0.00352,-0.00662,-0.00768,0.00759,-0.00415,0.00185,-0.00300,
-0.00300,0.00406,0.00397,0.00662
-0.00206,0.00246,0.00194,-0.00223,0.00149,0.00087,0.00087,0.00069,
0.00174,0.00211,-0.00119,0.00004,
-0.00126,0.00115,-0.00086,0.00097,0.00097,-0.00057,0.00073,-0.00054,0.00047,0.00112,0.00068,0.00131,
0.00116,-0.00142,-0.00142,0.00146,0.00161,-0.00055,0.00333,0.00294,
-0.00209,-0.00267,0.00019,0.00012,0.00012,-0.00023,-0.00139,-0.003290.00098,0.00037,-0.00107,-0.00087,0.00019,-0.00045,-0.00045,0.0
0016,
-0.00124,-0.00099,-0.00007,-0.00011,0.00045,0.00092,-0.00110,0.00077,0.00077,0.00191,-0.00133,-0.000720.00104,0.00097,0.00045,0.00140,
0.00062,0.00023,0.00023,0.00108,
0.00033,0.00104,-0.00085,-0.00084,0.00060,0.00085,-0.00049,0.00050,0.00050,-0.00046,0.00003,-0.00029,
-0.02719,-0.03543,0.02966,-0.00577,0.03214,-0.01772,-0.01483,-0.03461,0.02307,-0.01689,0.
02389
(4)通过式(12)、式(6)计算出Δw:Δw=
0,0.00006,0.00006,-0.00006,0.00006,0.00006,
-0.00006,0.00006,0.00004,0,0.00009,0.00009,
-0.00009,0.00009,0.00009,-0.00009,0.00009,0.00005,0,-0.00002,-0.00002,0.00002,-0.00002,-0.00002,0.00002,-0.00002,-0.00001,0.00001,0.00001,-0.00001,0.00001,0.00001,-0.00001,0.00001,0.00001,0,-0.00002,-0.00002,0.00002,
-0.00002,-0.00002,0.00002,-0.00002,-0.000010,0.00004,0.00004,-0.00004,0.00004,0.00004,
-0.00004,0.00004,0.00002,0,0.00002,0.00002,
-0.00002,0.00002,0.00002,-0.00002,0.00002,0.00001,0,0.00002,0.00002,-0.00002,0.00002,0.00002,-0.00002,0.00002,0.00001,0,-0.00001,
-0.00001,0.00001,-0.00001,-0.00001,0.00001,
-0.00001,-0.00001,0,0.00001,0.00001,-0.00001,0.00001,0.00001,-0.00001,0.00001,0.00001,
-0.00015,-0.00020,-0.00004,-0.00020,
-0.00003,-0.00013,-0.00005,-0.00003,
-0.00003,-0.00003,-0.00006,-0.00009,
0.00002,-0.00001,0.00002,-0.00004,-0.00002,
-0.00002,0.00001,-0.00001,0.
00041
E(wk)=0.01406
(5)如果E(wk)<ε,转到式(10),否则,用wk+1为权值和阈值计算误差E(wk+1)=0.01137。
wk+1=
0.84000,-0.85994,0.38006,0.50994,-0.17994,
-0.59994,0.66994,0.86006,0.56004,0.75000,
0.87009,-0.85991,0.46991,-0.20991,0.34009,
-0.46009,-0.44991,-0.74995,-0.83000,0.98998,0.77998,-0.89998,0.59998,0.34998,0.28002,
0.69998,0.84999,0.82000,-0.02999,0.87001,
-0.79001,0.59001,-0.66999,0.38999,-0.49999,0.37001,0.25000,0.58998,0.35998,0.69002,
0.60998,-0.75002,0.77002,0.84998,-0.29001,
-0.86000,-0.75996,0.54004,0.68996,-0.04996,
-0.02996,0.05996,0.36004,0.85002,-0.63000,
-0.23998,0.69002,0.55998,-0.11998,0.29002,
-0.10002,0.80002,0.64001,0.03000,0.05002,
-0.19998,-0.41002,0.49002,-0.33998,-0.85002,0.59002,0.32001,0.72000,0.66999,0.30999,
0.97001,0.42999,0.15999,0.75001,0.22999,
0.71999,0.87000,0.86001,-0.60999,-0.87001,0.50001,-0.50999,0.46999,-0.02999,0.30001,0.65985,0.85980,-0.72004,0.13980,-0.78003,0.42987,0.35995,0.83997,-0.56003,0.40997,
0.67994,0.67991,0.72002,-0.93001,0.20002,
-0.56004,0.65998,-0.17002,0.82001,-0.64001,0.
57959

本文发布于:2024-09-21 01:30:29,感谢您对本站的认可!

本文链接:https://www.17tex.com/tex/1/361726.html

版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系,我们将在24小时内删除。

标签:时间   作业   港口   集装箱船   预测
留言与评论(共有 0 条评论)
   
验证码:
Copyright ©2019-2024 Comsenz Inc.Powered by © 易纺专利技术学习网 豫ICP备2022007602号 豫公网安备41160202000603 站长QQ:729038198 关于我们 投诉建议