1. 引言
在集装箱多式联运运输组织中,枢纽港作为联运网络的“制高点”,对整个集装箱多式联运效率具有重要影响,其生产作业的高效组织和集装箱快速通过等受到了广泛关注。其中,装卸设备作为枢纽港内部作业的核心资源,其有效组织和协同调度等与实现枢纽港高效作业、节能减排等目标直接相关且尤其重要。
目前已经有很多学者针对枢纽港设备调度问题展开研究,王志刚 [1] 等针对自动化码头的设备调度和AGV路径规划问题,构建双层调度模型,并且设计双层遗传算法进行求解;Lu Yiqin [2] 提出了自动化集装箱码头调度的三阶段集成优化模型,以自动化集装箱码头装卸任务作业时间最小为目标;Hop [3] 等在解决集卡和场吊作业调度优化问题时建立了时间最小的单目标数学模型;代江涛 [4] 等考虑了集装箱码头装卸设备在不同作业状态下的能耗,建立了完工时间最小和能耗最低的双目标模型;秦琴等 [5] 等使用柔性流水车间调度理论解决了缓冲支架、AGV、双小车岸桥以及自动化轨道吊的协同调度问题,利用NEH启发式算法得到初始解的遗传算法进行求解;张笑菊 [6] 基于码头岸桥同贝同步装卸,构建了枢纽港集卡和岸桥联合调度模型,以完工时间最小为目标。
从以上研究来看,已有的建模场景以集装箱码头或者铁水联运为主,且主要针对作业完工时间、装卸成本等目标进行优化。从节能减排现实需求出发考虑多式联运枢纽港设备能耗和完工时间开展的研究仍较为缺乏。
2. 问题描述

Figure 1. Schematic diagram for loading and unloading of export boxes
图1. 出口箱装卸示意图
多式联运枢纽港集装箱作业主要有三个环节,一是在铁水联运中由门吊承担的集装箱从火车上取下,在公水联运中由正面吊承担的集装箱从货车上取下,在水水转运过程中由岸桥承担的集装箱从船上取下环节;二是集卡承担的码头与场区的水平作业环节;三是岸桥承担的装船环节,如图1所示。集装箱在枢纽港中转过程中,各种作业设备相互配合完成中转。设备的作业序列是多设备间协同调度需要考虑的核心问题,决定了枢纽港的整体作业水平,不合理的作业序列会导致设备间的彼此等待,装卸设备的利用率过低,产生不必要的设备能耗。综上本文以龙门吊、正面吊、集卡和岸桥协同调度问题为研究对象,综合考虑不同设备作业能力等约束,以作业完工时间和能耗最小为目标展开研究。
3. 数学模型
本文模型在参考文献 [7] 的基础上,拓展考虑了公水转运、水水转运多种联运模式,并以岸桥(卸)、正面吊、龙门吊、集卡、岸桥(装)作业阶段的设备协同调度为对象来构建。
3.1. 符号定义
R定义为铁水联运集装箱任务集合
,L定义为公水联运集装箱任务集合
,S定义为水水转运集装箱任务集合
,
表示阶段j的作业设备集合,
时分别代表集卡,岸桥,正面吊,龙门吊的作业阶段;
表示阶段j中装卸设备的工作效率;
表示阶段j中设备单位时间的能耗指标,
;
表示阶段j中设备单位时间的能耗指标,
;
表示电力能源折算系数;
表示柴油能源折算系数;
表示设备在阶段j的作业时间效率;Z表示虚拟任务箱;
表示在阶段j设备mk对任务箱r的装卸时间,
;
表示在阶段j设备mk对任务箱l的装卸时间,
;
表示在阶段j设备mk对任务箱s的装卸时间,
;
表示0~1变量,铁水联运在阶段j,任务箱r由设备mk操作为
,否则为0,
;
表示0~1变量,公水联运在阶段j,任务箱l由设备mk操作为
,否则为0,
;
表示0~1变量,水水转运在阶段j,任务箱s由设备mk操作为
,否则为0,
;
是0~1变量,铁水联运时在阶段j,当同一设备mk完成任务箱i的操作之后接着对任务箱
进行作业时
,反之为0,
,
,
;
是0~1变量,公水联运时在阶段j,当同一设备mk完成任务箱i的操作之后接着对任务箱
进行作业时
,反之为0,
,
,
;
是0~1变量,水水联运时在阶段j,当同一设备mk完成任务箱i的操作之后接着对任务箱
进行作业时
,反之为0,
,
,
;
;
决策变量,铁水联运在阶段j,设备mk对任务箱r开始作业的时间;
决策变量,公水联运在阶段j,设备mk对任务箱l开始作业的时间;
决策变量,水水转运在阶段j,设备mk对任务箱s开始作业的时间;
决策变量,铁水联运在阶段j,设备mk对任务箱r完成作业的时间;
决策变量,公水联运在阶段j,设备mk对任务箱l完成作业的时间;
决策变量,水水转运在阶段j,设备mk对任务箱s完成作业的时间。
3.2. 模型构建
在多式联运条件下,枢纽港中转装卸的协同优化问题的主要优化目标是作业完工时间和能耗。其中在集装箱的中转装卸过程中消耗电能的是岸桥和龙门吊,消耗柴油的是集卡和正面吊。为了方便计算装卸过程中各设备消耗的总能耗,引入设备能耗折算系数α,单位是kg标准煤/kg或kg标准煤/kWh,然后将不同的设备能耗统一转换成标准煤能耗,结合设备的单位能耗指标μ,单位为kWh/TEU或kg/TEU,以及各设备单位时间的作业效率
单位为kWh/h通过计算得到装卸设备的单位时间能耗
,单位为kWh/h或kg/h,如式(1)示:
(1)
进而构建“铁水”“公水”“水水”出口集装箱中转的最小能耗公式:
1) 铁水联运集装箱中转的最小能耗计算公式为:
(2)
2) 公水联运集装箱中转的最小能耗为:
(3)
3) 水水转运集装箱中转的最小能耗为:
(4)
多式联运出口集装箱各设备协同作业的目标模型为:
(5)
(6)
式(5)表示所有集装箱中转的总完工时间,式(6)表示最小化集卡、岸桥、正面吊、龙门吊、协同作业过程的总能耗。模型的条件约束主要是:
1) 作业序列约束.
(7)
(8)
(9)
(10)
(11)
(12)
(13)
(14)
(15)
(16)
(17)
(18)
(19)
(20)
式(7)代表铁水联运出口集装箱在各个阶段作业的唯一性;式(8)代表公水联运出口集装箱在各个阶段作业的唯一性;(9)铁水联运出口集装箱在当前作业阶段有且仅有一个紧前箱或紧后箱;式(10)公水联运出口集装箱在当前作业阶段只有一个紧前箱或紧后箱;式(11)水水转运出口集装箱在当前作业阶段只有一个紧前箱或紧后箱;式(12) (13)代表在铁水联运各阶段的起始箱是虚拟箱0,终止箱是虚拟箱Z;式(14) (15)代表在公水联运各阶段起始箱是虚拟箱0,终止箱是虚拟箱Z;式(16) (17)代表在水水转运各阶段起始箱是虚拟箱0,终止箱是虚拟箱Z;式(18)确定铁水联运各阶段设备的作业顺序;式(19)确定公水联运各阶段设备的作业顺序;式(20)确定水水转运各阶段设备的作业顺序
2) 设备能力约束。
(21)
(22)
(23)
(21)表示铁水联运阶段各阶段设备在可利用范围内作业;式(22)表示公水联运阶段各阶段设备在可利用范围内作业;式(23)表示水水转运阶段各阶段设备在可利用范围内作业。
3) 变量约束.
(24)
(25)
(26)
式(24) (25) (26)为决策变量约束。
4. 算法求解
设备调度问题是NP难问题 [8] ,遗传算法比较适合此类问题的求解。代江涛 [4] 等在遗传算法迭代时改进的自适应交叉变异概率公式,可以有效解决变异概率或交叉概率过大适应度高的染色体被破坏或二者的值过小迭代初期无法产生新的个体陷入局部最优解的问题,将公式应用到设备调度问题中以找到最优设备协同调度方案,使完工时间和能耗最小。
公式如下:
(27)
(28)
式(27)和(28)中,
是种群个体的平均适应度值;
是种群个体的最大适应度值;
是进行交叉操作的两个个体之中较大的适应度值;
是进行变异操作个体的适应度值;
是最大交叉概率,
;
是最小交叉概率,
;
是最大变异概率,
;
是最小变异概率,
。
4.1. 染色体编码与解码
本文采用实数编码的方式。以铁水中转为例,在遗传算法的初始阶段构造一个
的矩阵,根据各阶段的设备数量确定mk,设定为
然后随机向下取整生成部分的初始解。得到初始解后对数据进行处理使其中一列数据转化为一行,得到另一部分染色的编码最终表达方式,分别代表了第一阶段设备龙门吊作业、第二阶段设备集卡作业、第三阶段设备岸桥作业的设备作业序列。
4.2. 种群初始化
在约束条件(7)-(26)的限制下生成一定数量的个体,每个个体的任务编号和作业设备序号利用随机数产生生成的所有个体形成初始种群。
4.3. 适应度评价和选择
本文采取的是分别给两个目标函数分别赋予权重,将双目标问题转化为单目标问题,将目标函数
的权重系数设置为
,
的权重系数设置为
。则归一化后的目标函数 [9] 为:
(29)
4.4. 遗传操作
1) 选择
使用本文将采取随机通用采样的方式来选取父代个体,这种方式每随机选择一次就能够选择N个父代的个体。通过选择好的个体适应力分配求出各个体被选择比例,得出各个体中可选择的期望个体,构造出一个轮盘,每转动一个轮盘,N条指针所指到的个体被选作父体进行次代演化。
2) 交叉
选择两个染色体的随机位置进行断点,然后将两个断点前后的基因进行交换。
3) 变异
随机选择一染色体的两个基因交换它们的位置,即完成变异。
5. 算例分析
结合Sha等 [10] 的实验参数和武汉港相关作业数据开展算例研究,改进遗传算法(Improved Genetic Algorithm, IGA)使用Matlab进行编程实现。采用的计算机参数为AMD Ryzen 7 5800H with Radeon Graphics @3.20 GHz 16.0 GB。
5.1. 算例设计
设定初始数据按照铁路集装箱:公路集装箱:水路集装箱 = 4:5:7的比例设置,集装箱160个其中铁路集装箱为40个,公路集装箱50个,水路集装箱70。同时设置集卡10辆;岸桥7个(其中前三个岸桥负责水路集装箱卸下,后4个岸桥负责集装箱装船);正面吊3个;龙门吊3个,基本参数如表1所示。改进遗传算法的参数设置为,迭代次数300,交叉概率最大为0.9,最小为0.6,变异概率最大为0.1最小为0.01 [4] 时间权重0.8,能耗权重0.2 [7] 。

Table 1. Basic parameter information
表1. 基本参数信息
5.2. 算例求解
上述算例用于验证IGA的可行性,求解结果为:算法约在第190次收敛,求解得到的最小完工时间为6791 s,最小能耗为628.97 kg,算法收敛过程如图2所示。得到的出口箱设备作业序列如表2所示,其中铁路出口箱第一个集装箱由第1个龙门吊,第3辆集卡,第4个岸桥完成出口箱卸下与装船;公路出口箱第一个集装箱则是由第3个正面吊,第6辆集卡,第9个岸桥完成出口箱卸下与装船;水路出口箱第一个集装箱则是由第2个岸桥,第10辆集卡,第6个岸桥完成出口箱中转。

Table 2. Operation sequence of container equipment
表2. 集装箱设备作业顺序
5.3. 算法对比
为进一步验证算法改进的有效性,设置多组不同集装箱数量的算例,将GA和IGA两种算法进行对比,结果如表3和图3所示。其中,图3表示IGA与GA求解结果之间的偏差值,用Gap表示。表3中结果表明,IGA求解的最小完工时间和能耗结果都优于GA,改进遗传算法求解效率更高,缩短10%左右求解时间。图3中结果表明,随着集装箱数量变化带来的问题规模增加,在求解结果方面,IGA与GA之间的偏差呈现出逐渐增大趋势,表明IGA在求解更大规模问题时具有更好地求解性能。

Table 3. Comparison of Results for Different Container Volumes
表3. 不同集装箱量结果对比
注:
;
。
Figure 3. Gap Comparison Chart
图3. Gap对比图
6. 结语
本文以集装箱多式联运枢纽港设备协同调度问题为研究对象,研究了枢纽港集装箱中转作业涉及的集卡、岸桥、龙门吊、正面吊等设备协同调度问题,建立了以作业完工时间最小和能耗最小的双目标模型,并改进遗传算法实现问题求解,得到较优设备作业序列;扩大算例研究结果表明集装箱数量增加,IGA最小完工时间和能耗结果都优于GA,随着作业箱量的增加,改进算法对优化目标的提升程度更加显著。其中,箱量为1120时IGA求解得到最小完工时间比GA缩短了4173 s,Gap1优化了8.68%;IGA求解的能耗比GA减少了152.14 kg,Gap2优化了3.44%。后续的研究中将考虑多式联运过程中列车到达的不确定性等约束,将各阶段设备数量作为决策变量,深入研究出口箱量不同的情况下各阶段作业设备的最佳配比等问题。
基金项目
湖北省教育厅科学技术研究计划项目(编号:Q20211110)。
NOTES
*第一作者。
#通讯作者。