码头货轮集装箱装卸的优化问题升级版

1、-. z.码头货轮集装箱装卸的优化问题-. z.码头货轮集装箱装卸的优化问题摘要集装箱货币化已成为开展趋势，而港口开展渐渐滞后于集装箱的吞吐量，研究集装箱装卸的优化问题能有效扩大港口生产力，提高港口经济效益。本文将建立集卡线路规划模型和岸桥、集卡与龙门吊协同优化模型，通过禁忌搜索算法进展求解，并通过*港的数据对模型进展实证分析。对于提高装卸效率，降低装卸本钱这一问题，我们将其分解为线路规划、协同优化和模型检验三个子问题进展分析。针对问题一，我们建立了集卡线路规划模型。通过对*港前湾港集装箱码头QQCT的航拍图和雷达图进展分析，画出了码头泊位到堆场的平面图，按照相应的比例尺，得到实际码头与堆场间

2、的距离、各堆场间的相互距离。通过集卡行驶的速度，计算得到集卡从码头到堆场的时间、集卡在各堆场之间行驶的相互时间和集卡从堆场返回码头的时间。集卡在运输过程中，要尽量减少空集卡的行驶，即运送集装箱返回的途中携带需要装运到船上的集装箱。利用第一阶段的禁忌搜索算法，当所需装卸集装箱位置确定后，最短的行驶路线也就计算出来。针对问题二，我们建立了桥吊、集卡和龙门吊的协同优化模型。问题一计算的集卡最正确线路分配结果，继续作为桥吊、集卡和龙门吊协同优化的条件。第二阶段的禁忌搜素算法分析出最适宜的桥吊、集卡与龙门吊的比例，桥吊在不等待集卡的情况下效率高。通过协同优化，得到最高效率的设备分配比例。针对问题三，我们

3、汇总了附件中所有集装箱的装卸数据，对模型进展检验分析。以*前湾港区为例，通过带入实际数据，得到如下比例关系，即桥吊：集卡：龙门吊为2:10:5。2辆桥吊工作时配备10辆集卡，5辆轮式龙门吊；3辆桥吊工作时配备15辆集卡，7辆龙门吊；如此分配使相对本钱与效率到达最大化。本文的亮点在于：利用港口的雷达图和航拍图，绘制了港口的分布平面图，分析更贴近实际；以集卡线路规划为突破口，并以此为条件，建立了以集装箱类型为依据的集卡一站式效劳岸桥到堆场的线路标准化；对数据的分类处理，使计算简洁；协同了集卡、桥吊、龙门吊，采用两个阶段的禁忌搜索算法，将集装箱的装与卸混合在一起计算，比原来对集卡、桥吊，集卡、龙门吊

4、等局部优化更加贴近实际，大大提升了港口的运行效率，并且降低的了本钱。关键词：集卡，桥吊，龙门吊，线路规划，协同优化，禁忌搜索算法-. z.目录TOC o 1-3 h u HYPERLINK l _Toc5396 摘要 PAGEREF _Toc5396 1 HYPERLINK l _Toc13604 一、 问题重述 PAGEREF _Toc13604 3 HYPERLINK l _Toc17405 1.1问题的背景 PAGEREF _Toc17405 3 HYPERLINK l _Toc11905 1.2要解决的问题 PAGEREF _Toc11905 3 HYPERLINK l _Toc156

5、76 二、 问题分析 PAGEREF _Toc15676 3 HYPERLINK l _Toc11242 2.1概论 PAGEREF _Toc11242 3 HYPERLINK l _Toc15922 2.2问题一的分析 PAGEREF _Toc15922 5 HYPERLINK l _Toc8082 2.3问题二的分析 PAGEREF _Toc8082 5 HYPERLINK l _Toc22845 2.4问题三的分析 PAGEREF _Toc22845 5 HYPERLINK l _Toc24458 三、模型假设 PAGEREF _Toc24458 5 HYPERLINK l _Toc32

6、40 四、符号说明 PAGEREF _Toc3240 6 HYPERLINK l _Toc30572 五、模型建立与求解 PAGEREF _Toc30572 7 HYPERLINK l _Toc2621 5.1集卡线路优化模型 PAGEREF _Toc2621 7 HYPERLINK l _Toc7201 模型分析 PAGEREF _Toc7201 7 HYPERLINK l _Toc31835 模型建立 PAGEREF _Toc31835 7 HYPERLINK l _Toc31530 模型求解 PAGEREF _Toc31530 9 HYPERLINK l _Toc15543 5.2岸桥、

7、集卡和龙门吊的协同优化模型 PAGEREF _Toc15543 10 HYPERLINK l _Toc10449 模型的分析 PAGEREF _Toc10449 10 HYPERLINK l _Toc21044 模型建立 PAGEREF _Toc21044 10 HYPERLINK l _Toc19747 模型求解 PAGEREF _Toc19747 10 HYPERLINK l _Toc30331 5.3模型的检验分析 PAGEREF _Toc30331 13 HYPERLINK l _Toc28086 5.3.1模型分析 PAGEREF _Toc28086 13 HYPERLINK l _

8、Toc10542 5.3.2模型建立 PAGEREF _Toc10542 15 HYPERLINK l _Toc25208 5.3.3模型求解 PAGEREF _Toc25208 15 HYPERLINK l _Toc17910 六、 模型评价 PAGEREF _Toc17910 15 HYPERLINK l _Toc13852 6.1模型优点 PAGEREF _Toc13852 15 HYPERLINK l _Toc31673 6.2模型缺点 PAGEREF _Toc31673 16 HYPERLINK l _Toc30441 6.3模型改良 PAGEREF _Toc30441 16 HYP

9、ERLINK l _Toc14004 七、参考文献 PAGEREF _Toc14004 17 HYPERLINK l _Toc25327 附录 PAGEREF _Toc25327 18-. z.问题重述1.1问题的背景集装箱码头是海陆联运的枢纽站，在各个经济体的贸易中都占据着举足轻重的地位。港口的装卸货能力在一定程度上代表着一个港口的生产力，在集装箱吞吐量不断增大而港口开展渐渐滞后的现状下，研究港口集装箱装卸的优化问题就显得尤为重要。影响集装箱装卸效率的因素主要有装卸设备的硬件配备和在现有设备根底上对各方资源的优化协调程度。本文将以*港前湾港集装箱码头QQCT为现实背景，来考虑码头货轮集装箱装

10、卸的优化问题。QQCT坐落于*胶州湾的前港港区，具有水深域阔、不淤不冻、避风浪的优良码头条件，可全天候停靠第六代及以上集装箱船舶。泊位长度3400米，泊位水深-17.5米，巷道水深-15米，堆场面积225平方米。该港区配备了世界上最先进、最大型的新型桥吊，可装卸目前世界上最大型的超巴拿马型集装箱船。1.2要解决的问题本问将以QQCT为原型，利用卸船箱和装船箱的数据，简化港口条件，仅在投入2个和3个QC时，构造数学模型，分别使相应本钱尽量小，效率尽量高。为方便解决问题，我们将问题进展了分解，主要包括以下几个问题：基于集装箱类型的集卡线路一站式优化线路标准化；岸桥、集卡和龙门吊的协同调度；以*港的

11、实际数据对模型进展检验分析。问题分析2.1概论针对以*港为原型的码头集装箱装卸优化问题，我们首先通过对*港实地情况的分析，将目标港进一步细化，选择了四个港中最专业的集装箱装运码头，即*前湾集装箱码头。通过航拍图图1和雷达图图2，我们初步拿到QQCT的港口分布图，并基于此，我们绘制了港口分布立体图图3和平面图图4，来使分析更加具体和形象。 图1 航拍图 图2 雷达图来源：必应地图截图图3 立体图来源：百度图片 图4 平面图泊位确定2.2问题一的分析这是一个线路规划问题，通过对港口平面图、集卡运行速度、桥吊和龙门吊的工作效率分析，得到了集卡的最优工作路线。问题的特点在于集卡在运输进口、过境等集装箱

12、时是联合搭配的。问题的难点在于堆场与堆场之间可以相互通行，集卡返回码头时携带需要装船的集装箱，而且需要卸的集装箱往往和需要装的集装箱数目不一致，加大了理解和计算的难度。2.3问题二的分析这是一个资源优化问题，通过对桥吊、集卡、龙门吊三个因素进展综合考虑，得到三种资源的最优配比。问题的特点在于问题一的结论影响问题二的分析，桥吊、集卡、龙门吊之间的相互调配相互制约，但存在整体最优的情况。问题的难点是在各因素相互制约的条件下寻找最优解，合理有效利用第一问的数据。2.4问题三的分析这是一个检验分析问题，通过对QQCT集装箱装卸的数据进展分析整理，对模型进展实证检验，以证明模型的可靠性和实用性。问题的特

13、点是将实际数据带入模型确定2个桥吊和3个桥吊最正确的资源配比。问题的难点在于6000多个数据的带入存在一定的技术困难和理解误差，会影响分析结果的合理性。三、模型假设1、船一旦靠岸只选择一个泊位，且泊位距各堆场和最近；2、每个集装箱的装或卸只进展单次作业，作业完成后不再考虑；3、各堆场位置，由集装箱的种类确定；4、集装箱装卸次序；5、只考虑一条船的情况，装箱与卸箱都是同一艘货轮；6、假设2个20寸的集装箱视为1个40寸的集装箱；7、将需装运的集装箱种类简化为过境箱、进口箱和出口箱。四、符号说明符号符号说明轮船进口、过境集装箱集合 出口集装箱集合 进口、过境集装箱数量出口集装箱的数量进口、过境集装

14、箱所存堆场的集合 出口集装箱所存堆场的集合 进口、过境集装箱所存堆场位置 出口集装箱所存堆场位置 集卡从船到堆场行驶所用的时间 集卡从堆场行驶到所用的时间 集卡从堆场到船所行驶的时间 集卡是否经过路线桥吊集合卸箱桥吊装箱桥吊集卡集合集卡桥吊完成一个集装箱装卸的时间龙门吊完成一个集装箱装卸的时间符号符号说明集装箱在船的位置集装箱从到所用的时间卸箱工作完成的时间装箱工作完成的时间无限大的数箱由桥吊运送后，与由桥吊完成，是前面的任务，则箱被集卡送，则箱与箱均由集卡运送，在前，五、模型建立与求解5.1集卡线路优化模型模型分析通过查阅资料，我们将集卡的工作流程抽象如下图5船堆场堆场集卡集卡集卡 图5 集

15、卡工作流程图模型建立集卡的空驾驶距离之和最小：(1)堆场对应具体货物，如下为集卡调度模型的表示：(2) (3) (4)约束条件：时， (5)现实生活中往往进口箱与出口箱数量不一致，参加虚拟进、出口箱集合(virtual)时，(6)(7)时， (8) (9)其中，(1)式为集卡行驶距离最短的的表达式；(2)式为行驶时间最小的目标函数；(3)、(4)式为进口、过境箱与出口箱数目相等的一一搭配；(5)式为变量的约束条件；(6)、(7)式进口、过境箱数大于出口箱数的虚拟搭配；(8)、(9)式进口、过境箱数小于出口箱数的虚拟搭配。模型求解算法求解逻辑：桥吊的工作效率为1min/箱，集卡的行驶速度为20k

16、m/h，桥吊与集卡的比例为1：5，桥吊到堆场的距离在2001500之间随机产生。计算五个不同脚本的数据。利用MATLAB 2014a 计算，两阶段禁忌搜索算法计算结果如下表1：卸箱数桥吊数集卡行驶距离km作业时间min计算时间初始解禁忌算法初始解禁忌算法100/1179.88161.301981810.6400/1711.37630.367767458.11000/12900.672839.803010297732.22000/14117.834001.563771370543.32500/14779.104712.514710465153.0 表1 两阶段禁忌搜索算法从结果上看，每个桥吊工作

17、1000到2000个效率最高。假设低于1000箱每个桥吊会造成桥吊等待集卡，导致工作效率大大下降。假设箱数高于2000箱每个桥吊，则造成集卡等待桥吊，造成码头拥挤，次序混乱，结果不令人满意。5.2岸桥、集卡和龙门吊的协同优化模型模型的分析以模型一集卡路线的最优化结果作为模型二的条件，综合考虑桥吊、集卡、龙门吊的情况，协同考虑混合装卸，利用禁忌搜索第二阶段的桥吊调度阶段最优解。模型建立(10)约束条件： (11)(12) (13) (14)(15)(16)(17)(18) (19) (20) (21) (22) (23) (24)(25)其中，(11)(14)式说明每个集装箱都有一辆集卡、桥吊进

18、展操作；(18)式为桥吊与龙门吊对每个集装箱操作先后次序一定且唯一； (19)(20)式为桥吊完成一次集装箱作业的时间关系； (21)(22)式为集卡完成一次集装箱作业的时间关系； (23)式为集装箱约束条件，如i,j集装箱由同一辆集装箱运输； (24)式为集卡装卸集装箱的顺序；(25)式为变量取值*围的约束。模型求解求解算法逻辑：桥吊的工作效率为1min/箱，集卡的行驶速度为20km/h，桥吊与集卡的比例为1：5，桥吊到堆场的距离在2001500之间随机产生。5000个集装箱等待装卸。利用MATLAB 2014a 计算，混合装卸调度模型的结果如表2：桥吊/集卡穿插作业集卡行驶距离(km)总装

19、卸时间(min)6/18是129503815否828036356/24是128653655否818535106/30是124853615否81253485 表2 混合装卸调度计算结果5.3模型的检验分析5.3.1数据分析根据码头集装箱装卸的数据表，我们得到如下信息表3工程编号123456789卸箱量61344922331585200装箱量006103870303813种类FHFRGPHCHHOTRFRH0H表3装、卸船箱在堆场的分布数据来源：两个表格数据汇总*前湾港区设配数据如表4：设备效率设备种类装集装箱个/min卸集装箱个/min桥吊11龙门吊0.50.5表4 *前湾港设配数据集卡从桥吊运

20、送集装箱到各堆场的行驶时间如下表5：堆场号123456789送箱时间单位：min/次4.213.662.511.661.662.513.364.215.06表5 集卡从桥吊运送集装箱到各堆场的时间集卡从i堆场到j堆场的单次行驶时间如表6(单位：min/次)：堆场123456789103.854.705.558.409.2510.1010.9511.8023.8503.854.75.558.49.2510.1010.9534.703.8503.854.705.558.409.2510.1045.554.703.8503.854.705.558.409.2558.405.554.703.8503.

21、854.705.558.4069.258.405.554.703.8503.854.705.55710.109.258.405.554.703.8503.854.70810.9510.109.258.405.554.703.8503.85911.8010.9510.09.258.405.554.703.850表6集卡从i堆场到j堆场的单次行驶时间每个集卡可以装1个40寸集装箱，2个20寸集装箱。5.3.2模型建立利用模型一、二的算法，带入进展数据进展算法合理性检验。5.3.3模型求解利用MATLAB 2014a 计算结果如下表7桥吊/集卡穿插作业集卡行驶距离km总装卸时间min6/18是186

22、485493否1192352346/24是154384386否988242146/30是149824338否97504182表7 模型求解结果当题目中所给轮船拥有2辆桥吊的时候，10辆集卡，5辆轮式龙门吊可以到达预期效率，拥有3辆桥吊的时候，15辆集卡，7辆轮式龙门吊可以高效完成装卸任务。题目中所给轮船巨大，完成任务时2辆桥吊最正确效率需要3.6天可以完成装卸任务，3辆桥吊时需要2.9天完成装卸任务，假设加上运输时间，此次轮船的作业时间会太长。建议增加桥吊、集卡、龙门吊的数量以便48小时内完成装卸和运输任务，方便港口的运行。模型评价6.1模型优点本模型采用两阶段的禁忌搜索算法，搜索*围以及精度

23、都有很大提高，计算过程清晰，题目数据处理采用分类汇总的方法，大大简化了计算难度。第一阶段禁忌搜索从集卡调度入手，准确分析得出集卡在当集装箱装卸次序的情况下的最正确行驶路线，即一旦装卸的集装箱确定，就可以程序化进展路线分配，大大减少了集卡运送集装箱的总路程。第一阶段集卡调度的计算结果作为第二阶段禁忌搜索的条件，使模型全局考虑桥吊、集卡、龙门吊的协同调度，到达了在装卸运输本钱一定的情况下，效率最高的设备搭配方式。采取混合装卸模型，更加贴近实际，算得最正确的桥吊、集卡、龙门吊之间的比例，方便工作人员进展设备分配的决策。模型最后计算出了装卸、运输总体最小时间，为在集装箱数目确定的情况下，投入设备的多少

24、取得最正确效益的决策提供了参考。6.2模型缺点本模型的集卡调度为静态调度模型，即每辆集卡满载且只可携带1个40寸，2个20寸的集装箱，实际生活中可根据需要，集卡既可以满载，也可不满载，集卡的调度静态调度存在缺陷。本模型只考虑了一条船的装卸作业，实际可能有多艘轮船需要完成装卸作业，轮船选择泊位时候过于简化。6.3模型改良将静态集卡调度改为集卡可加一节车厢，可以满载或不满载，根据实际情况确定的动态集卡调度模型。考虑多艘轮船的装卸，根据多艘轮船上的集装箱装卸次序，可用设备数量分配，综合考虑定下假设轮船选择不同泊位，确定同一个泊位不同船只的先后装卸次序。这样模型会更加贴近实际生产生活，方便决策者的决策

25、。七、参考文献1钱继锋. 集装箱码头岸桥集卡堆场作业方案的优化D.交通大学,2014.2冯春焕. 集装箱码头泊位岸桥集卡调度优化研究D.*海事大学,2011.3单浩. 集装箱码头泊位、岸桥和集卡协同调度优化研究D.*海事大学,2013.4曾庆成,杨忠振,陆靖. 集装箱码头同贝同步装卸调度模型与算法J. 交通运输工程学报,2010,01:88-93. 5曾庆成,胡祥培,杨忠振. 集装箱码头泊位分配-装卸桥调度干扰管理模型J. 系统工程理论与实践,2010,11:2026-2035. 6曾庆成. 集装箱码头装卸作业集成调度模型与方法D.*海事大学,2008. 7韩晓龙,牟少莉. 基于CHC算法的集

26、卡与岸桥协调调度优化问题J. *理工大学学报(信息与管理工程版),2014,02:233-236+245.附录% 禁忌搜索算法解决装卸混合调度问题 function BestShortcut,theMinDistance=TabuSearch clear; clc; Clist=装卸集装箱的贝位号; ContinaerNum=size(Clist,1);%问题的规模,即集装箱数目 dislist=zeros(ContinaerNum); for i=1:ContinaerNum for j=1:ContinaerNum dislist(i,j)=(Clist(i,1)-Clist(j,1)2+

27、(Clist(i,2)-Clist(j,2)2)0.5; end end TabuList=zeros(ContinaerNum);% (tabu list) TabuLength=round(ContinaerNum*(ContinaerNum-1)/2)0.5);%禁忌长度(tabu length) Candidates=200;%候选集的个数 (全部领域解个数) CandidateNum=zeros(Candidates,ContinaerNum);%候选解集合 S0=randperm(ContinaerNum);%随机产生初始解 BSF=S0; BestL=Inf; clf; figu

28、re(1); stop = uicontrol(style,toggle,string ,stop,background,white); tic; p=1; StopL=80*ContinaerNum; while pContinaerNum*(ContinaerNum-1)/2 disp(候选解个数不大于n*(n-1)/2!); break;end ALong(p)=Fun(dislist,S0); i=1; A=zeros(Candidates,2); while i=Candidates M=ContinaerNum*rand(1,2); M=ceil(M); if M(1)=M(2)

29、A(i,1)=ma*(M(1),M(2); A(i,2)=min(M(1),M(2); if i=1 isa=0; else for j=1:i-1 if A(i,1)=A(j,1) & A(i,2)=A(j,2) isa=1; break; else isa=0; end end end if isa i=i+1; else end else end end BestCandidateNum=100;%保存前BestCandidateNum个最好候选解作为第二阶段的条件 BestCandidate=Inf*ones(BestCandidateNum,4); F=zeros(1,Candida

30、tes); for i=1:Candidates CandidateNum(i,:)=S0; CandidateNum(i,A(i,2),A(i,1)=S0(A(i,1),A(i,2);F(i)=Fun(dislist,CandidateNum(i,:); if i=BestCandidateNum BestCandidate(i,2)=F(i); BestCandidate(i,1)=i; BestCandidate(i,3)=S0(A(i,1); BestCandidate(i,4)=S0(A(i,2); else for j=1:BestCandidateNum if F(i)BestC

31、andidate(j,2) BestCandidate(j,2)=F(i); BestCandidate(j,1)=i; BestCandidate(j,3)=S0(A(i,1); BestCandidate(j,4)=S0(A(i,2); break; end end end end %对BestCandidate JL,Inde*=sort(BestCandidate(:,2); SBest=BestCandidate(Inde*,:); BestCandidate=SBest; if BestCandidate(1,2)BestL BestL=BestCandidate(1,2); S0=CandidateNum(BestCandidate(1,1),:); BSF=S0; for m=1:ContinaerNum for n=1:ContinaerNum if TabuList(m,n)=0 TabuList(m,n)=TabuList(m,n)-1; end end end TabuList(BestCandidate(1,3),BestCandidate(1,4)=TabuLengt