CN114626579A - 基于改进布谷鸟算法的集装箱船舶配载优化方法 - Google Patents
基于改进布谷鸟算法的集装箱船舶配载优化方法 Download PDFInfo
- Publication number
- CN114626579A CN114626579A CN202210162159.8A CN202210162159A CN114626579A CN 114626579 A CN114626579 A CN 114626579A CN 202210162159 A CN202210162159 A CN 202210162159A CN 114626579 A CN114626579 A CN 114626579A
- Authority
- CN
- China
- Prior art keywords
- container
- ship
- formula
- stowage
- time
- 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.)
- Pending
Links
- 241000544061 Cuculus canorus Species 0.000 title claims abstract description 34
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000005457 optimization Methods 0.000 title claims abstract description 21
- 235000005770 birds nest Nutrition 0.000 claims abstract description 41
- 235000005765 wild carrot Nutrition 0.000 claims abstract description 41
- 238000013178 mathematical model Methods 0.000 claims abstract description 9
- 230000008278 dynamic mechanism Effects 0.000 claims abstract description 6
- 238000003860 storage Methods 0.000 claims description 19
- 239000011159 matrix material Substances 0.000 claims description 15
- 238000010586 diagram Methods 0.000 claims description 4
- 235000015170 shellfish Nutrition 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 238000009826 distribution Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 description 12
- 238000011160 research Methods 0.000 description 8
- 230000002452 interceptive effect Effects 0.000 description 6
- 230000003993 interaction Effects 0.000 description 4
- 230000007246 mechanism Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000006872 improvement Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000010845 search algorithm Methods 0.000 description 2
- 108091026890 Coding region Proteins 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000007306 turnover Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial 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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/40—Business processes related to the transportation industry
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- General Health & Medical Sciences (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Educational Administration (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Primary Health Care (AREA)
- Ship Loading And Unloading (AREA)
Abstract
本发明涉及一种基于改进布谷鸟算法的集装箱船舶配载优化方法,其在满足船舶适航性的基础上,建立了以集装箱水平运输时间、岸桥作业时间最短且倒箱量最少为目标的多目标多约束的数学模型。运用改进的布谷鸟算法进行求解,首先在莱维飞行阶段,分别从种群整体、当代最优、自身个体中抽取鸟巢信息,建立多种位置更新模式;其次,采用动态机制控制发现概率;最后通过计算每个鸟巢的目标函数值,来找到最优解。本发明基于码头方考虑的目标函数模型和改进的布谷鸟算法进行集装箱船舶配载优化,得到的配载结果可以有效提高码头的作业效率,降低码头作业成本。
Description
技术领域
本发明涉及船舶配载技术领域,具体涉及一种基于改进布谷鸟算法的集装箱船舶配载优化方法。
背景技术
船舶配载计划是一个NP-Hard问题,现有的文献主要由船舶的预配计划与码头的实配计划研究组成。集装箱船舶配载在实际的操作过程中是一个多目标、多约束的的组合优化问题。配载计划的优劣程度关系到集装箱船舶在港口的装卸效率、船舶的适航性等。近些年来,我国港口集装箱吞吐量成倍增长,这给集装箱船舶配载带来了很多压力。
对于集装箱配载方面的研究,现在主要的研究是从船方的角度,而从码头角度对集装箱船舶实配研究的文献很少,且考虑堆场翻箱和作业过程冲突等因素较少。本发明基于以上研究的不足之处从集装箱船舶配载整体均衡出发,对集装箱码头配载优化问题进行研究。
发明内容
本发明的目的在于提供一种基于改进布谷鸟算法的集装箱船舶配载优化方法,其从码头方的角度综合考虑实际情况,建立了以最小化堆场翻箱量、最小的集装箱水平运输距离和最小的岸桥作业时间为配载计划的决策目标,并结合改进的布谷鸟算法进行码头配载优化,以提高优化的可行性和有效性。
为实现上述目的,本发明采用的技术方案是:
一种基于改进布谷鸟算法的集装箱船舶配载优化方法,其在满足船舶适航性的基础上,建立了以集装箱水平运输时间、岸桥作业时间最短且倒箱量最少为目标的多目标多约束的数学模型;并运用改进的布谷鸟算法进行求解:首先在莱维飞行阶段,分别从种群整体、当代最优、自身个体中抽取鸟巢信息,建立多种位置更新模式;其次,采用动态机制控制发现概率;最后通过计算每个鸟巢的目标函数值,来找到最优解。
所述数学模型包括目标函数及约束条件;
目标函数如下:
约束条件如下:
xip,zip=0或1 (12)
约束(2)、(3)为箱位限制约束,其中,式(2)表示每个集装箱只能被配载到一个船箱位上;式(3)表示每个船箱位只能配载一个集装箱;式(4)表示待装集装箱离开箱区必须符合桥吊作业的时间安排;式(5)表示箱区集装箱轻不压重的堆垛要求;式(6)表示位于同一堆栈不同层数的待装集装箱离开箱区的先后顺序和翻箱的关系;
式(7)表示待装集装箱在船舶箱位紧挨着的下方有集装箱即不能悬空;式(8)表示待装集装箱与船舶贝位的匹配满足预配要求;式(9) 表示船舶贝位同列集装箱的装载重量不能超过单列额定重量;式(10) 表示确保离场集装箱数量与待装集装箱数量相同;式(11)表示箱区作业量限制,防止作业量过大;式(12)决策变量类型;式(13)表示相关变量的取值范围;
c为列号,c=1,2,…C;r为层号,r=1,2,…R;b1为堆场箱区编号, b1=1,2,…B1;b2为堆场箱区区块编号,b2=1,2,…B2;b3为堆场小箱区编号,b3=1,2,…B3;b为堆场箱区具体编号,b=(b1,b2,b3); i,j为码头堆场待装集装箱编号,i=1,2…,I,j=1,2,…J;m为贝位号, m=1,2,…M;p为待装船集装箱的船箱位编号,p=1,2,…,P;εp为堆场箱位p的开始装箱时间;U(p)为位于船舶箱位p(k,c,r) 同贝同列上一层p(k,c,r+1)的箱位;t表示t-1到t的时间间隔,t=1, 2,…T,Δi为与集装箱i同一堆栈的上层集装箱集合;σ1为堆场翻一次箱所需要的时间;Wi为集装箱i的重量;
nbt为t时间段内箱区b已完成场到船任务的集装箱数量;
为岸桥p完成集装箱i到相应船箱位所用的时间,vp表示岸桥p 作业范围内待装船集装箱的数量;W为船舶贝位中单列允许的集装箱最大载重量;τip为集装箱从堆场场箱位到船箱位的水平运输时间;BLTip为场箱位为i装载到船箱位p的离场时间,BLTip=εp-τip,其中p∈P;为箱区b在时间段t的箱区作业量阈值;γip为二维0-1矩阵,若待装船集装箱i符合预配图要求从而配载到船舶所对应的船舶箱位p(k,c, r),γip为1,否则为0;θip为二维0-1矩阵,若待装船集装箱i位于箱区b,则矩阵对应的元素为1,否则为0;Q为一个极大的正数。
采用改进的布谷鸟算法对目标函数进行求解具体包括以下步骤:
1)设置种群的规模,初始的位置,迭代的次数等参数;
2)计算种群中所有个体的目标函数值,并保留最优解;
3)根据式(14)进行位置的更新,在位置更新之后计算每个个体的目标函数值,通过精英策略将最优解保留至下一代;
cacuhy的数学表达式为:
其中,g=1为尺度参数;
5)通过式(17)计算发现概率,比较pa与随机数H的大小,如果 H>Pa,则丢弃劣质鸟巢,继续更新;否则,将当前的鸟巢位置保留下来;
5)计算经过步骤4)变化后每个鸟巢的目标函数值,并将最优鸟巢保留下来进行下一轮的迭代;
6)如果输出的结果满足算法的终止条件;如果不是则返回步骤3) 继续进行迭代,直到迭代的结果满足结束条件为止。
采用上述方案后,在满足船舶适航性的基础上,建立了以集装箱水平运输时间、岸桥作业时间最短且倒箱量最少为目标的多目标多约束的数学模型。运用改进的布谷鸟算法进行求解,首先在莱维飞行阶段,分别从种群整体、当代最优、自身个体中抽取鸟巢信息,建立多种位置更新模式;其次,采用动态机制控制发现概率;最后通过计算每个鸟巢的目标函数值,来找到最优解。本发明基于码头方考虑的目标函数模型和改进的布谷鸟算法进行集装箱船舶配载优化,得到的配载结果可以有效提高码头的作业效率,降低码头作业成本。
附图说明
图1为改进布谷鸟算法的流程图;
图2为本发明实施例的预配载示意图。
具体实施方式
本发明揭示了一种基于改进布谷鸟算法的集装箱船舶配载优化方法,在满足船舶适航性的基础上,建立了以集装箱水平运输时间、岸桥作业时间最短且倒箱量最少为目标的多目标多约束的数学模型。运用改进的布谷鸟算法进行求解,首先在莱维飞行阶段,分别从种群整体、当代最优、自身个体中抽取鸟巢信息,建立多种位置更新模式;其次,采用动态机制控制发现概率;最后通过计算每个鸟巢的目标函数值,来找到最优解。具体如下:
1、模型建立
码头方在进行集装箱船舶配载的过程中,其主要考虑的是要保证集装箱船舶能够依照靠泊计划按时离港,所以装卸速度和效率是关键因素。本文在综合考虑研究意义及现实意义的情况下,决定建立一个先卸后装,并且把倒箱量最少、装卸效率最大作为目标的数学模型,并用改进的布谷鸟算法进行求解。
1.1假设条件
(1)在码头进行实际配载之前,预配船图需由船方提供给码头方,保证待配载集装箱的各种信息,如尺寸大小都已知。且本文是从码头的视角出发,所以主要考虑的是码头方装船配载方案的优化。
2)本文考虑的是单船配载作业,且为了避免装卸同时进行时设备冲突的增加以及翻箱数的上升,不考虑边装边卸。考虑到建模的方便,同时考虑问题求解的实现条件,做出以下合理的假设:
(a)危险品箱,冷藏箱等特殊箱在船上的配载位置已设定,且不考虑踮脚箱;
(b)桥吊和龙门吊等港口机械设备数量充足,配置合理,有足够的码头作业机械保证配载的效率不受大的影响,作业过程中不受机械故障停休的影响;
(c)假定不同工艺的集装箱的作业时间相同,不受集装箱作业工艺的影响;
(d)预配待出口集装箱的数量不超过计划可装载的船箱位数量;
1.2参数定义
c为列号,c=1,2,…C;
r为层号,r=1,2,…R;
b1为堆场箱区编号,b1=1,2,…B1;
b2为堆场箱区区块编号,b2=1,2,…B2;
b3为堆场小箱区编号,b3=1,2,…B3;
b为堆场箱区具体编号,b=(b1,b2,b3);
i,j为码头堆场待装集装箱编号,i=1,2…,I,j=1,2,…J;
m为贝位号,m=1,2,…M;
p为待装船集装箱的船箱位编号,p=1,2,…,P;
εp为堆场箱位p的开始装箱时间;
U(p)为位于船舶箱位p(k,c,r)同贝同列上一层p(k,c,r+1) 的箱位;
Δi为与集装箱i同一堆栈的上层集装箱集合;
σ1为堆场翻一次箱所需要的时间;
Wi为集装箱i的重量;
nbt为t时间段内箱区b已完成场到船任务的集装箱数量;
W为船舶贝位中单列允许的集装箱最大载重量;
τip为集装箱从堆场场箱位到船箱位的水平运输时间;
γip为二维0-1矩阵,若待装船集装箱i符合预配图要求从而配载到船舶所对应的船舶箱位p(k,c,r),γip为1,否则为0;
θip为二维0-1矩阵,若待装船集装箱i位于箱区b,则矩阵对应的元素为1,否则为0;
Q为一个极大的正数。
1.3决策变量
xip为0-1变量,当且仅当集装箱i配载到船箱位p上时为1,否则为 0。
zij为0-1变量,当且仅当集装箱i在集装箱j之前配载时为1,否则为0。其中j∈Δ(i)。
vbpt为第t时间段箱区b完成对应船箱位p配载任务的作业时间。
2.4目标函数及约束条件
约束条件:
xip,zip=0或1 (12)
约束(2)、(3)为箱位限制约束,其中,
式(2)表示每个集装箱只能被配载到一个船箱位上;
式(3)表示每个船箱位只能配载一个集装箱;
式(4)表示待装集装箱离开箱区必须符合桥吊作业的时间安排;
式(5)表示箱区集装箱轻不压重的堆垛要求;
式(6)表示位于同一堆栈不同层数的待装集装箱离开箱区的先后顺序和翻箱的关系;
式(7)表示待装集装箱在船舶箱位紧挨着的下方有集装箱即不能悬空;
式(8)表示待装集装箱与船舶贝位的匹配满足预配要求;
式(9)表示船舶贝位同列集装箱的装载重量不能超过单列额定重量;
式(10)表示确保离场集装箱数量与待装集装箱数量相同;
式(11)表示箱区作业量限制,防止作业量过大;
式(12)决策变量类型;
式(13)表示相关变量的取值范围。
2、改进的布谷鸟算法(ICS)设计
2.1、莱维飞行机制的改进
在传统的CS算法中,莱维飞行虽然具有短步长与长步长相互交替,能有效探索全局空间,但其中长步长将短步长组成的集聚分隔开来,短步长出现的频率较高,长步长出现的频率较低,且在迭代的中后期,个体间相似度较高,算法陷入局部极值的概率增大。并且,在莱维飞行的机制中,鸟巢信息位置的更新单一,采集的个体信息片面,实时性、综合性的种群信息在寻优过程中会被忽略。所以,在保留莱维飞行的特性的基础上,每一步迭代的过程中引进新的个体信息,将多类别交互信息融入位置更新策略中,以此在搜索过程中扩展信息的多样性,探索新的搜索领域。交互信息反映种群内个体之间适应度的差异,根据差异可及时调整算法进化的方向,提高搜索的速度,加速算法的收敛进度[13]。改进的莱维机制表达式如下所示:
cacuhy的数学表达式为:
其中,g=1为尺度参数。
式(14)中,构建3组位置的更新计算策略,通过减法构建种群信息交互机制,分别从种群整体(N1),当代最优(N2),当前个体(N3) 获取鸟巢的交互信息,将此信息嵌入位置更新策略中,及时掌握鸟巢内适应度的优劣情况况,再利用随机方式选择位置更新策略。为了保证交互信息的全面性和总体性,采用Cauchy算子对随机鸟巢进行扰动,柯西算子能无重复地遍历种群中每一个个体,使算法在最大的限度上收集种群个体的信息,拓展交互信息的领域。
采用更新方式1时,在柯西算子的扰动下,通过计算种群内两个随机鸟巢的之间的差值来衡量所有个体的差异度(构建N1与N1的交互信息),在N1&N1(种群整体)范围进行全局寻优。在更新策略2 中,依据柯西算子随机选择一个鸟巢,并评价N1与N2的适应度之间的差异(构建N1与N2的交互信息),在N1&N2的维度上更新鸟巢位置(进行以N2为核心的局部搜索);更新策略3中,构建N1与 N3的交互信息,考量到N1与N3之间存在差异,在N1&N3内寻优搜索(进行以N3为核心的局部搜索)。这样的改进策略,依据N1、N2、 N3之间的交互信息进行多维度的位置更新,更全面的评估了种群差异与寻优情况,充分的扩展了算法的寻找范围;随机的位置更新策略,使算法在全局与局部之间动态变化寻优,避免了勘测开发不平衡问题。
2.2改进的发现概率
算法研究发现,如果发现概率pa越大,其收敛速度就越快,搜索能力越强;pa值越小,越不利于收敛速度及搜索精度。因此,本发明采用动态机制:
改进的布谷鸟算法的主要步骤为:
1)设置种群的规模,初始的位置,迭代的次数等参数;
2)计算种群中所有个体的目标函数值,并保留最优解;
3)根据式(14)进行位置的更新,在位置更新之后计算每个个体的目标函数值,通过精英策略将最优解保留至下一代。
4)通过式(17)计算发现概率,比较pa与随机数H的大小,如果 H>Pa,则丢弃劣质鸟巢,继续更新;否则,将当前的鸟巢位置保留下来。
5)计算经过步骤4)变化后每个鸟巢的目标函数值,并将最优鸟巢保留下来进行下一轮的迭代。
6)如果输出的结果满足算法的终止条件;如果不是则返回步骤3) 继续进行迭代,直到迭代的结果满足结束条件为止。
2.3编码方式
集装箱码头船舶配载的核心问题在于在满足约束条件,使目标函数达到最优的情况下确定待装集装箱在船上的位置,使得最初的配载计划既满足船方的配载要求也能满足码头方的配载目标和要求,所以怎样确定集装箱在堆场的位置怎么安排到船舶舱位的位置,所以本节采用图2 的编码方式来确定配载的位置信息。
本文的假设条件中,假设码头方的各作业计划已知,其中包括桥吊作业和场桥作业的计划已知,所以相对应作业船贝的作业顺序和作业开始时间也可以得到确定。所以对应堆场安排到具体的船位置也是可以确定的,船箱位的配载次序也是知道的。需要解决的问题是怎么确定场箱位与船箱位的匹配关系。这里设计1~n的自然数对场箱位进行了实数编码,用1~m代表船舶仓箱位编码顺序。编码的上层代表场箱位,下层每一位代表船舶舱箱位。表1展示了n个场箱位匹配船箱位的匹配图。
4 | 7 | 9 | 19 | 58 | … | 66 | 89 | 99 | n |
1 | 3 | 6 | 11 | 17 | … | 44 | 55 | 64 | m |
表1
如上表1所示表示的是场箱位编号为4的待装集装箱匹配船箱位为 1的位置,第19位场箱位装载船箱位为11位置。按照这样的编码方式编码出的任意序列在满足相关约束的条件下代表了配载计划的一个初始解(即对应的待装集装箱场箱位对应的船箱位位置),这样的编码方式是为了后续更好的利用模型进行求解,得到最佳的配载结果。
3、算例求解
3.1、实验数据
为验证本文所提出的算法模型的有效性,本文根据厦门远海集装箱码头的实际配载数据,选取了ALS金星(ALS Venus),航次为9525948 的集装箱船的09贝位为实验研究对象进行模型的算例分析。其中包括集装箱船舶09贝位的船箱位信息,待装在场集装箱信息以及船舶09贝位的预配计划图的相关数据如下。
09贝位的85个预配船箱位信息如下表2所示。其中主要包括:船箱位、箱型(高箱、平箱)、卸货港、尺寸(20ft、40ft)、空箱和重箱;其中在箱型一列中GP代表平箱、HC代表高箱;空/重一列,F代表重箱, E代表空箱。
表2
待配载到09贝位的85个在场集装箱的信息如下表3
表3
通过实际调查码头的运作情况,设定相关参数,箱区一个待装集装箱不考虑翻箱的场桥平均作业时间即平均场桥就位和取箱的时间为 2min,由拖车的平均移动速度,根据集装箱所在箱区到船贝之间的距离,可以得到配载过程中每个集装箱的水平运输所需的时间。堆场箱区平均一次翻箱所用的时间为3.3min,式(1)中的权重系数设置为:α1=0.4,α2=0.4,α3=0.3,对改进的布谷鸟算法参数设置如下:迭代次数种群规模大小为150,最大迭代次数Tmax=250,CS的发现概率为Pa=0.25。
3.2结果分析
基于以上数据带入到算法中,来对配载方案进行优化求解,得到改进前后布谷鸟算法优化配载结果的实验结果及收敛曲线见表4。
表4
为了检验改进的布谷鸟搜索算法应用于本文模型的收敛性,将本文采用09贝位的85个待配载的在场集装箱,设立了10组不同的算例,进行实验并对结果进行数据统计,可以得到改进的布谷鸟搜索算法中目标函数值与迭代次数之间的关系,如图所示。
在Tmax=250的条件下,目标函数值先随着迭代次数的增加先不断减少,随着迭代次数的增加,目标函数值逐渐趋于平稳状态。并且随着迭代次数的增加,布谷鸟搜索的速度加快,目标函数值逐渐趋于精确,最终找到最优解。从图可以看出当迭代次数到达80次时,目标函数值达到平稳状态,不再改变,证明了改进的布谷鸟算法可以在有 限的迭代次数之后找到目标函数的最优解,从而说明了本文所提出的 该进的布谷鸟算法解决提出的集装箱配载模型在解决集装箱配载问 题是有效的。
综上,改进的布谷鸟算法在求解集装箱配载方面,本文所构建的基于码头方考虑的目标函数模型和改进的布谷鸟算法,得到的配载结果可以有效提高码头的作业效率,降低码头作业成本,表明模型和算法是可靠的。
以上所述,仅是本发明实施例而已,并非对本发明的技术范围作任何限制,故凡是依据本发明的技术实质对以上实施例所作的任何细微修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (3)
1.一种基于改进布谷鸟算法的集装箱船舶配载优化方法,其特征在于:所述方法在满足船舶适航性的基础上,建立了以集装箱水平运输时间、岸桥作业时间最短且倒箱量最少为目标的多目标多约束的数学模型;并运用改进的布谷鸟算法进行求解:首先在莱维飞行阶段,分别从种群整体、当代最优、自身个体中抽取鸟巢信息,建立多种位置更新模式;其次,采用动态机制控制发现概率;最后通过计算每个鸟巢的目标函数值,来找到最优解。
2.根据权利要求1所述的一种基于改进布谷鸟算法的集装箱船舶配载优化方法,其特征在于:所述数学模型包括目标函数及约束条件;
目标函数如下:
约束条件如下:
xip,zip=0或1 (12)
约束(2)、(3)为箱位限制约束,其中,式(2)表示每个集装箱只能被配载到一个船箱位上;式(3)表示每个船箱位只能配载一个集装箱;式(4)表示待装集装箱离开箱区必须符合桥吊作业的时间安排;式(5)表示箱区集装箱轻不压重的堆垛要求;式(6)表示位于同一堆栈不同层数的待装集装箱离开箱区的先后顺序和翻箱的关系;
式(7)表示待装集装箱在船舶箱位紧挨着的下方有集装箱即不能悬空;式(8)表示待装集装箱与船舶贝位的匹配满足预配要求;式(9)表示船舶贝位同列集装箱的装载重量不能超过单列额定重量;式(10)表示确保离场集装箱数量与待装集装箱数量相同;式(11)表示箱区作业量限制,防止作业量过大;式(12)决策变量类型;式(13)表示相关变量的取值范围;
c为列号,c=1,2,…C;r为层号,r=1,2,…R;b1为堆场箱区编号,b1=1,2,…B1;b2为堆场箱区区块编号,b2=1,2,…B2;b3为堆场小箱区编号,b3=1,2,…B3;b为堆场箱区具体编号,b=(b1,b2,b3);i,j为码头堆场待装集装箱编号,i=1,2…,I,j=1,2,…J;m为贝位号,m=1,2,…M;p为待装船集装箱的船箱位编号,p=1,2,…,P;εp为堆场箱位p的开始装箱时间;U(p)为位于船舶箱位p(k,c,r)同贝同列上一层p(k,c,r+1)的箱位;t表示t-1到t的时间间隔,t=1,2,…T,Δi为与集装箱i同一堆栈的上层集装箱集合;σ1为堆场翻一次箱所需要的时间;Wi为集装箱i的重量;
nbt为t时间段内箱区b已完成场到船任务的集装箱数量;
3.根据权利要求2所述的一种基于改进布谷鸟算法的集装箱船舶配载优化方法,其特征在于:采用改进的布谷鸟算法对目标函数进行求解具体包括以下步骤:
1)设置种群的规模,初始的位置,迭代的次数等参数;
2)计算种群中所有个体的目标函数值,并保留最优解;
3)根据式(14)进行位置的更新,在位置更新之后计算每个个体的目标函数值,通过精英策略将最优解保留至下一代;
cacuhy的数学表达式为:
其中,g=1为尺度参数;
4)通过式(17)计算发现概率,比较pa与随机数H的大小,如果H>Pa,则丢弃劣质鸟巢,继续更新;否则,将当前的鸟巢位置保留下来;
5)计算经过步骤4)变化后每个鸟巢的目标函数值,并将最优鸟巢保留下来进行下一轮的迭代;
6)如果输出的结果满足算法的终止条件;如果不是则返回步骤3)继续进行迭代,直到迭代的结果满足结束条件为止。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210162159.8A CN114626579A (zh) | 2022-02-22 | 2022-02-22 | 基于改进布谷鸟算法的集装箱船舶配载优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210162159.8A CN114626579A (zh) | 2022-02-22 | 2022-02-22 | 基于改进布谷鸟算法的集装箱船舶配载优化方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114626579A true CN114626579A (zh) | 2022-06-14 |
Family
ID=81900277
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210162159.8A Pending CN114626579A (zh) | 2022-02-22 | 2022-02-22 | 基于改进布谷鸟算法的集装箱船舶配载优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114626579A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115936323A (zh) * | 2022-08-15 | 2023-04-07 | 苏州诀智科技有限公司 | 基于集装箱船舶的智能舱内配载方法、***及计算机介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103544586A (zh) * | 2013-10-25 | 2014-01-29 | 东北大学 | 一种提高集装箱船舶稳定性和岸桥作业效率的配载方法 |
CN109388848A (zh) * | 2018-08-30 | 2019-02-26 | 大连海事大学 | 一种基于双四十英尺岸桥的集装箱船舶配载方法 |
-
2022
- 2022-02-22 CN CN202210162159.8A patent/CN114626579A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103544586A (zh) * | 2013-10-25 | 2014-01-29 | 东北大学 | 一种提高集装箱船舶稳定性和岸桥作业效率的配载方法 |
CN109388848A (zh) * | 2018-08-30 | 2019-02-26 | 大连海事大学 | 一种基于双四十英尺岸桥的集装箱船舶配载方法 |
Non-Patent Citations (3)
Title |
---|
MARICHELVAM, MK; PRABAHARAN, T AND YANG, XS: "Improved cuckoo search algorithm for hybrid flow shop scheduling problems to minimize makespan", APPLIED SOFT COMPUTING, no. 19, 30 June 2014 (2014-06-30), pages 93 - 101 * |
王永金;黄鹏飞;陈麒龙;储雨峰;张宁;: "基于改进布谷鸟算法的集装箱码头配载优化研究", 中国航海, vol. 46, no. 3, 30 September 2023 (2023-09-30), pages 65 - 71 * |
申经伟: "基于改进布谷鸟算法的港口群海铁联运网络优化", 中国优秀硕士学位论文全文数据库工程科技Ⅱ辑, no. 2021, 15 January 2021 (2021-01-15), pages 034 - 4 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115936323A (zh) * | 2022-08-15 | 2023-04-07 | 苏州诀智科技有限公司 | 基于集装箱船舶的智能舱内配载方法、***及计算机介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103544586B (zh) | 一种提高集装箱船舶稳定性和岸桥作业效率的配载方法 | |
Parreño et al. | A GRASP algorithm for the container stowage slot planning problem | |
CN111861205B (zh) | 一种自动化集装箱码头出口箱空间分配方法 | |
US8352404B2 (en) | Container vessel stowage planning | |
CN111498696B (zh) | 一种双吊具场桥调度优化方法及装置 | |
CN107451720B (zh) | 面向自动化集装箱码头的船次作业计划评价方法 | |
CN112434870A (zh) | 一种集装箱区垂直布置的双自动化场桥动态调度方法 | |
CN110599000A (zh) | 自动化码头翻箱评估方法、箱位分配方法及相关装置 | |
Wang et al. | Tree based searching approaches for integrated vehicle dispatching and container allocation in a transshipment hub | |
CN114626579A (zh) | 基于改进布谷鸟算法的集装箱船舶配载优化方法 | |
Xue et al. | An ant colony algorithm for yard truck scheduling and yard location assignment problems with precedence constraints | |
CN116629424A (zh) | 新型集装箱堆场作业模式的自动化双场桥智能调度方法 | |
CN109019056B (zh) | 一种集装箱堆场双层高架轨道垂直运输设备调度方法 | |
CN115271140A (zh) | 基于动态规划的集装箱船舶配载箱位预留方法 | |
Li et al. | An integer linear programming for container stowage problem | |
Nguyen et al. | Heuristic algorithms for constructing transporter pools in container terminals | |
KR100841682B1 (ko) | 안벽크레인의 양적하 작업 계획 산출 방법 | |
Park et al. | Optimizing stacking policies using an MOEA for an automated container terminal | |
CN108256769B (zh) | 基于nsga-ii改进的场桥调度方法 | |
CN112926873A (zh) | 集装箱码头船舶配载深度学习混合训练*** | |
Hamdi et al. | A heuristic for the container stacking problem in automated maritime ports | |
Xia et al. | Loading sequencing problem in container terminal with deep Q-learning | |
CN106516807A (zh) | 一种提高多种钢铁产品混装海运稳定性的配载方法 | |
CN101893854B (zh) | 集装箱码头装船倍内发箱顺序决策*** | |
CN112053046A (zh) | 一种带有时间窗的自动化集装箱码头agv重进重出路径规划方法 |
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 |