数学建模——码头货轮集装箱装卸的优

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

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

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

4、的装与卸混合在一起计算，比原来对集卡、桥吊，集卡、龙门吊等部分 优化更加贴近实际，大大提升了港口的运行效率，并且降低的了成本。关键词：集卡，桥吊，龙门吊，线路规划，协同优化，禁忌搜索算法2目录摘要0.1、 问题重述31.1 问题的背景31.2 要解决的问题32、 问题分析32.1 概论32.2 问题一的分析52.3 问题二的分析52.4 问题三的分析53、 模型假设5四、符号说明64、 模型建立与求解85.1 集卡线路优化模型 .85.1.1 模型分析85.1.2 模型建立.85.1.3 模型求解.95.2 岸桥、集卡和龙门吊的协同优化模型 115.2.1 模型的分析115.2.2 模型建立

5、125.2.3 模型求解 125.3 模型的检验分析 155.3.1 模型分析155.3.2 模型建立 185.3.3 模型求解 18六、模型评价1 96.1 模型优点1 96.2 模型缺点1 96.3 模型改进20七、参考文献2 1附录2 214、问题重述1.1 问题的背景集装箱码头是海陆联运的枢纽站，在各个经济体的贸易中都占据着举足轻 重的地位。港口的装卸货能力在一定程度上代表着一个港口的生产力，在集装 箱吞吐量不断增大而港口发展渐渐滞后的现状下，研究港口集装箱装卸的优化 问题就显得尤为重要。影响集装箱装卸效率的因素主要有装卸设备的硬件配备 和在现有设备基础上对各方资源的优化协调程度。本文

6、将以青岛港前湾港集装 箱码头（QQCT）为现实背景，来考虑码头货轮集装箱装卸的优化问题。QQCT坐落于青岛胶州湾的前港港区，具有水深域阔、不淤不冻、避风浪 的优良码头条件，可全天候停靠第六代及以上集装箱船舶。泊位长度 3400米, 泊位水深-17.5米，巷道水深-15米，堆场面积225平方米。该港区配备了世界上 最先进、最大型的新型桥吊，可装卸目前世界上最大型的超巴拿马型集装箱船。1.2 要解决的问题本问将以QQCT为原型，利用卸船箱和装船箱的数据，简化港口条件，仅 在投入2个和3个QC时，构造数学模型，分别使相应成本尽量小，效率尽量高 为方便解决问题，我们将问题进行了分解，主要包括以下几个问

7、题：（1）基于集装箱类型的集卡线路一站式优化（线路标准化）；（2）岸桥、集卡和龙门吊的协同调度；（3）以青岛港的实际数据对模型进行检验分析。、问题分析2.1 概论针对以青岛港为原型的码头集装箱装卸优化问题，我们首先通过对青岛港 实地情况的分析，将目标港进一步细化，选择了四个港中最专业的集装箱装运 码头，即青岛前湾集装箱码头。通过航拍图（图 1）和雷达图（图2）,我们初步 拿到QQCT的港口分布图，并基于此，我们绘制了港口分布立体图（图 3）和 平面图（图4）,来使分析更加具体和形象。图1航拍图图2雷达图（来源：必应地图截图）图3立体图（来源：百度图片）图4平面图（泊位确定）2.2 问题一的分析

8、这是一个线路规划问题，通过对港口平面图、集卡运行速度、桥吊和龙门 吊的工作效率分析，得到了集卡的最优工作路线。问题的特点在于集卡在运输 进口、过境等集装箱时是联合搭配的。问题的难点在于堆场与堆场之间可以相 互通行，集卡返回码头时携带需要装船的集装箱，而且需要卸的集装箱往往和 需要装的集装箱数目不一致，加大了理解和计算的难度。2.3 问题二的分析这是一个资源优化问题，通过对桥吊、集卡、龙门吊三个因素进行综合考 虑，得到三种资源的最优配比。问题的特点在于问题一的结论影响问题二的分 析，桥吊、集卡、龙门吊之间的相互调配相互制约，但存在整体最优的情况。 问题的难点是在各因素相互制约的条件下寻找最优解，

9、合理有效利用第一问的 数据。2.4 问题三的分析这是一个检验分析问题，通过对 QQCT集装箱装卸的数据进行分析整理， 对模型进行实证检验，以证明模型的可靠性和实用性。问题的特点是将实际数 据带入模型确定2个桥吊和3个桥吊最佳的资源配比。问题的难点在于 6000多 个数据的带入存在一定的技术困难和理解误差，会影响分析结果的合理性。三、模型假设1、船一旦靠岸只选择一个泊位，且泊位距各堆场和最近；2、每个集装箱的装或卸只进行单次作业，作业完成后不再考虑;3、各堆场位置已知，由集装箱的种类确定；4、集装箱装卸次序已知；5、只考虑一条船的情况，装箱与卸箱都是同一艘货轮；6、假设2个20寸的集装箱视为1个

10、40寸的集装箱；7、将需装运的集装箱种类简化为过境箱、进口箱和出口箱。四、符号说明付方符号说明轮船进口、过境集装箱集合 Nr出口集装箱集合 NaDr进口、过境集装箱数量出口集装箱的数量一进口、过境集装箱所存堆场的集合r nDadir出口集装箱所存堆场的集合na进口、过境集装箱所存堆场位置 lid：出口集装箱所存堆场位置Ijtsdir集卡从船到堆场d：行驶所用的时 间17 r 7 a di da集卡从堆场dir行驶到da所用的时间tdas集卡从堆场d：到船所行驶的时间Rij1or0集卡是否经过路线ij 丫 桥吊集合qr卸箱桥吊装箱桥吊集卡集合集卡桥吊完成一个集装箱装卸的时间 龙门吊完成一个集装箱

11、装卸的时 间付方符号说明Mi,i i集装箱在船的位置集装箱从i到i 1所用的时间卸箱工作完成的时间装箱工作完成的时间Yq无限大的数箱i由桥吊q运送后,Yiq 1or0Yiq1与1由桥吊q完成,i是i前面的任务）则 1or0gic箱i被集卡c送）则gg 1or0gijc箱i与箱j均由集卡c运送）i在j前,gjc 1,or,gjic五、模型建立与求解5.1 集卡线路优化模型5.1.1 模型分析通过查阅资料，我们将集卡白工作流程抽象如下（图 5）3T 7集堆场I堆场w图5集卡工作流程图5.1.2模型建立tdasRdaS jjj ma(1)集卡的空驾驶距离之和最小:Min ijtsdrRsdrRddt

12、ji mr i i i dr j dai j堆场对应具体货物，如下为集卡调度模型的表示:Mintj Rdjda a(2)约束条件:Rdrda1 j dad d i jmr ma 时,i dr_rRrda1 i du dj(5：j daRd. 0or1现实生活中往往进口箱与出口箱数量不一致，加入虚拟进、出口箱集合(virtual)mr时,R Vidjvirtualmr时,Rjadidji virtualRHr1 i drdi iRd：dja i dri virtualRidjaj da(6)(8)Rd i ii drvirtual其中，(1)式为集卡行驶距离最短的的表达式；式为行驶时间最小的目标

13、函数；(3)、(4)式为进口、过境箱与出口箱数目相等的一一搭配;(5)式为变量的约束条件；(6)、(7)式进口、过境箱数大于出口箱数的虚拟搭配；(8)、(9)式进口、过境箱数小于出口箱数的虚拟搭配。5.1.3模型求解算法求解逻辑:桥吊的工彳效率为1min/箱，集卡的行驶速度为20km/h ,桥吊与集卡的比 例为1: 5,桥吊到堆场的距离在 2001500之间随机产生。计算五个不同脚本 的数据。利用MATLAB 2014a计算，两阶段禁忌搜索算法计算结果如下表 1:卸箱 数集卡行驶距 离(km)作业时间(min)计算 时间桥吊 数初始 解禁忌 算法100/1179.88161.30400/171

14、1.37630.361000/12900.672839.802000/14117.834001.562500/14779.104712.51初始 解禁忌 算法1981810.6776745823771370543.34710465153.0表1两阶段禁忌搜索算法从结果上看，每个桥吊工作1000到2000个效率最高。若低于1000箱每个桥吊会造成桥吊等待集卡，导致工作效率大大下降。若箱数高于2000箱每个桥吊，则造成集卡等待桥吊，造成码头拥挤，次序混乱，结果不令人满意5.2岸桥、集卡和龙门吊的协同优化模型5.2.1模型的分析以模型一集卡路线的最优化结果作为模型二的条件，

15、综合考虑桥吊、集卡、 龙门吊的情况，协同考虑混合装卸，利用禁忌搜索第二阶段的桥吊调度阶段最 优解。165.2.2模型建立Min max min max ：, a(10)约束条件:1,(11)Yjq Qa1,(12)gicC1,(13)g jc c C1,(14)Yi ri M1,Mr,Qr(15)Yjj qa j M a1,Qa(16)gijcri M1,Mr,(17)gijc j M a1,Na,(18)MHA1Yiqri M r,Qr(19)Mjjt iastsstsdirA1A1Yjj qagjc , ij Ma,A1 gjc ,Qa(20)Mr,j(21)i Mr,j(22)g ijc

16、 c CRj, iM r,(23)gic0.5 0.5 gjc gjicgic, iMr,j Ma,(24)Yj,gic,Yii g,gjj c 00r1(25)其中，(11卜(14)式说明每个集装箱都有一辆集卡、桥吊进行操作;21(15卜(18)式为桥吊与龙门吊对每个集装箱操作先后次序一定且唯(19卜(20)式为桥吊完成一次集装箱作业的时间关系；(21卜(22)式为集卡完成一次集装箱作业的时间关系；(23)式为集装箱约束条件，如i,j集装箱由同一辆集装箱运输；(24)式为集卡装卸集装箱的顺序；(25)式为变量取值范围的约束。5.2.3模型求解求解算法逻辑：开曲桥吊的工彳效率为1min/箱，集

17、卡的行驶速度为20km/h ,桥吊与集卡的比 例为1: 5,桥吊到堆场的距离在2001500之间随机产生。5000个集装箱等待 装卸。利用MATLAB 2014a计算，混合装卸调度模型的结果如表 2:桥吊/集 卡交叉作业集卡行驶距离(km)总装卸时间(min)6/18是129503815否828036356/24是128653655否818535106/30是124853615否81253485表2混合装卸调度计算结果5.3模型的检验分析5.3.1数据分析根据码头集装箱装卸的数据表，我们得到如下信息(表3)项目编号123456789卸箱 量61344922331585200装箱00613803

18、0383量071种类FHFRGPHCHHOTRFRH0H表3装、卸船箱在堆场的分布 （数据来源：两个表格数据汇总）青岛前湾港区设配数据如表4:、设氤效率 设备翁类、装集装箱（个/ min）卸集装箱（个/ min）桥吊11龙门吊0.50.5表4青岛前湾港设配数据集卡从桥吊运送集装箱到各堆场的行驶时间如下表5:堆场号123456789送箱时 间（单4.23.2.1.1.2.53.34.25.0位：min/次）1665166661616表5集卡从桥吊运送集装箱到各堆场的时间集卡从i堆场到j堆场的单次行驶时间如表 6（单位：min/次）:堆 场123456789103.854.7 05.558.4 0

19、9.2510.1010.9511.8023.8503.854.75.558.49.2510.1010.9534.703.8503.854.7 05.558.4 09.2510.1045.554.703.8503.854.705.558.409.2558.4 05.554.7 03.8503.854.7 05.558.4 069.28.45.54.73.803.84.75.5235050550510.9.28.45.54.73.83.84.77a r010505055010.10.9.28.45.54.73.83.889510505050511.10.10.9.28.45.54.73.89809

20、50505050表6集卡从i堆场到j堆场的单次行驶时间每个集卡可以装1个40寸集装箱,2个20寸集装箱。5.3.2模型建立利用模型一、二的算法，带入进行数据进行算法合理性检验。5.3.3模型求解利用MATLAB 2014a计算结果如下表7桥吊/集 卡交叉作 业集卡行驶距离(km)总装卸时间(min)6/18是186485493否1192352346/24是154384386否988242146/30是149824338否97504182表7模型求解结果当题目中所给轮船拥有2辆桥吊的时候，10辆集卡，5辆轮式龙门吊可以达 到预期效率，拥有3辆桥吊的时候，15辆集卡，7辆轮式龙门吊可以高效完成装

21、卸任务。题目中所给轮船巨大，完成任务时 2辆桥吊最佳效率需要3.6天可以完 成装卸任务，3辆桥吊时需要2.9天完成装卸任务，若加上运输时间，此次轮船 的作业时间会太长。建议增加桥吊、集卡、龙门吊的数量以便48小时内完成装卸和运输任务，方便港口的运行。六、模型评价6.1 模型优点本模型采用两阶段的禁忌搜索算法，搜索范围以及精度都有很大提高，计 算过程清晰，题目数据处理采用分类汇总的方法，大大简化了计算难度。第一阶段禁忌搜索从集卡调度入手，精确分析得出集卡在当集装箱装卸次 序已知的情况下的最佳行驶路线，即一旦装卸的集装箱确定，就可以程序化进 行路线分配，大大减少了集卡运送集装箱的总路程。第一阶段集

22、卡调度的计算结果作为第二阶段禁忌搜索的条件，使模型全局 考虑桥吊、集卡、龙门吊的协同调度，达到了在装卸运输成本一定的情况下， 效率最高的设备搭配方式。采取混合装卸模型，更加贴近实际，算得最佳的桥 吊、集卡、龙门吊之间的比例，方便工作人员进行设备分配的决策。模型最后 计算出了装卸、运输总体最小时间，为在集装箱数目确定的情况下，投入设备 的多少取得最佳效益的决策提供了参考。6.2 模型缺点本模型的集卡调度为静态调度模型，即每辆集卡满载且只可携带1个40寸, 2个20寸的集装箱，实际生活中可根据需要，集卡既可以满载，也可不满载，集卡的调度静态调度存在缺陷。本模型只考虑了一条船的装卸作业，实际可能有多

23、艘轮船需要完成装卸作 业，轮船选择泊位时候过于简化。6.3 模型改进将静态集卡调度改为集卡可加一节车厢，可以满载或不满载，根据实际情 况确定的动态集卡调度模型。考虑多艘轮船的装卸，根据多艘轮船上的集装箱装卸次序，可用设备数量 分配，综合考虑定下若轮船选择不同泊位，确定同一个泊位不同船只的先后装 卸次序。这样模型会更加贴近实际生产生活，方便决策者的决策。25七、参考文献1钱继锋.集装箱码头“岸桥一集卡一堆场”作业计划的优化D.北京交通大学,2014.2冯春焕.集装箱码头泊位一岸桥一集卡调度优化研究D.大连海事大学,2011.3单浩.集装箱码头泊位、岸桥和集卡协同调度优化研究D.大连海事大学,20

24、13.4曾庆成，杨忠振，陆靖.集装箱码头同贝同步装卸调度模型与算法J.交通运输工程学报，2010,01:88-93.5曾庆成，胡祥培，杨忠振.集装箱码头泊位分配-装卸桥调度干扰管理模型J.系统工程理论与实践，2010,11:2026-2035.6曾庆成.集装箱码头装卸作业集成调度模型与方法D.大连海事大 学,2008.7韩晓龙，牟少莉.基于CHCT法的集卡与岸桥协调调度优化问题J.武汉理 工大学学报（信息与管理工程版）,2014,02:233-236+245.附录%禁忌搜索算法解决装卸混合调度问题function BestShortcut,theMinDistance=TabuSearch c

25、lear;clc;Clist=装卸集装箱的贝位号 ;ContinaerNum=size(Clist,1);%问题的规模，即集装箱数目 dislist=zeros(ContinaerNum);for i=1:ContinaerNumfor j=1:ContinaerNumdislist(i,j)=(Clist(i,1)-ClistG1)A2+(Clist(i,2)-ClistG2)A2)A0.5;endendTabuList=zeros(ContinaerNum);% (tabu list)TabuLength=round(ContinaerNum*(ContinaerNum-1)/2)A0.5

26、);% 禁忌长度 (tabu length)Candidates=200;%候选集白个数(全部领域解个数)CandidateNum=zeros(Candidates,ContinaerNum);% 候选解集合S0=randperm(ContinaerNum);%随机产生初始解BSF=S0;BestL=Inf; clf;figure(1);stop = uicontrol('style','toggle','string','stop','background','white');tic;p=1；St

27、opL=80*ContinaerNum;while p<StopLif Candidates>ContinaerNum*(ContinaerNum-1)/2disp('候选解个数不大于n*(n-1)')；break;endALong(p)=Fun(dislist,S0);i=1;A=zeros(Candidates,2);while i<=CandidatesM=ContinaerNum*rand(1,2);M=ceil(M);if M(1)=M(2)A(i,1)=max(M(1),M(2);A(i,2)=min(M(1),M(2);if i=1 isa=0;

28、 elsefor j=1:i-1if A(i,1)=A(j,1) && A(i,2)=A(j,2)isa=1;break;elseisa=0;endendendif isai=i+1;else endelseendendBestCandidateNum=100;% 保留前 BestCandidateNum 个最好候选解作为第 二阶段的条件BestCandidate=Inf*ones(BestCandidateNum,4);F=zeros(1,Candidates);for i=1:CandidatesCandidateNum(i,:)=S0;CandidateNum(i,A(i

29、,2),A(i,1)=S0(A(i,1),A(i,2); F(i)=Fun(dislist,CandidateNum(i,:);if i<=BestCandidateNumBestCandidate(i,2)=F(i);BestCandidate(i,1)=i;BestCandidate(i,3)=S0(A(i,1);BestCandidate(i,4)=S0(A(i,2);elsefor j=1:BestCandidateNumif F(i)<BestCandidate(j,2)BestCandidate(j,2)=F(i);BestCandidate(j,1)=i;BestCa

30、ndidate(j,3)=S0(A(i,1);BestCandidate(j,4)=S0(A(i,2); break; end end end end %对 BestCandidate JL,Index=sort(BestCandidate(:,2);SBest=BestCandidate(Index,:); BestCandidate=SBest;if BestCandidate(1,2)<BestLBestL=BestCandidate(1,2);S0=CandidateNum(BestCandidate(1,1),:);BSF=S0; for m=1:ContinaerNum for n=1:ContinaerNum if TabuList(m,n)=0TabuList(m,n)=TabuList(m,n)-1; end end end TabuList(BestCandidate(1,3),BestCandidate(1,4)=TabuLength; else for i=1:BestCandidateNumif TabuList(BestCandidate(i,3),BestCandidate(i,4)=0 S0=CandidateNum(BestCandidate(i,1),:);for m=1:ContinaerNum for n=1:Contina