考虑场桥效率的集卡失约优化仿真_赵旭.pdf

收稿日期:20210414修回日期:20210515第 40 卷第 2 期计算机仿真2023 年 2 月文章编号:10069348(2023)02019007考虑场桥效率的集卡失约优化仿真赵旭，张宝国，索浩(大连海事大学交通运输工程学院，辽宁大连 116026)摘要:码头外集卡预约系统能够有效的减少港口拥堵，提高港口作业效率，但是失约率是影响预约系统的重要因素。为研究失约对预约系统的影响，降低外集卡在闸口和堆场排队时间，建立考虑场桥资源配置的具有失约现像的集卡预约优化模型。模型利用马尔可夫和排队论分别计算外集卡在闸口和堆场的等待时间，应用 matlab 仿真，对模型进行深入探讨并给出详细的数值性能分析，同时量化了失约现象对预约系统损耗的影响，得出了不同失约率成本最低值。关键词:预约排队;失约现象;马尔可夫;排队论;仿真中图分类号:U691.33文献标识码:BSimulation of Set Card Failure Optimization ConsideringField Bridge EfficiencyZHAO Xu，ZHANG Baoguo，SUO Hao(College of Transportation Engineering，Dalian Maritime University，Dalian Liaoning116026，China)ABSTACT:The reservation system of collecting cards outside the wharf can effectively reduce port congestion andimprove the efficiency of port operation，but the cancellation rate is an important factor affecting the reservation sys-tem In order to study the influence of missed appointment on the reservation system and reduce the queuing time ofexternal set card at the gate and yard，an optimization model of set card reservation with missed appointment imageconsidering the resource allocation of field bridge was established In this model，Markov and queuing theory wereused to calculate the waiting time of the outer set card at the gate and the yard respectively Matlab simulation wasused to conduct indepth discussion on the model and give detailed numerical performance analysis Meanwhile，theinfluence of nonappointment phenomenon on the loss of the reservation system was quantized，and the minimum costvalue of different nonappointment rate was obtainedKEYWODS:Appointment queue;Miss phenomenon;Markov;Queuing theory;Simulation1引言近年来，由于我国日益完善的集装箱港口硬件设备，集装箱吞吐量快速增加，港口堆场集装箱进出口流动大，这对集卡的运输效率提出了更高的要求，集卡预约系统能够降低集卡港内周转时间和排队等待时间，但是受客观因素影响，有部分集卡难以在预约时间段到达，这就造成了集卡的失约现象。现实情况中，集卡预约系统的实施效果是受到现实因素干扰的，比如交通堵塞、天气变化等带来的不确定性因素，增加了集卡的失约率，干扰了以降低排队时间为目标的集卡预约系统，并且一旦失约集卡的数量过多，甚至会造成原本预约模式的不可行。在集卡预约系统中，集卡的早到和晚到或未预约集卡的进入港口都会对港口的服务水平产生影响。这种干扰会造成港口设备的重叠使用，并且增加集卡的港内周转时间，并且造成集卡公司的营运成本增加。国内外学者在集卡预约方面己开展了众多前沿性的研究，取得了重大进展。一些研究致力于评价集卡预约系统的效果1，典型研究如 Giulianoa 和 Erera2 分析了集卡预约系统对闸口排队和碳排放的影响。Namboothiria 和 Erera3 以最小运输成本为目标构建了集卡调度模型，研究了集卡预约对集卡车队的影响。Zhao 和 Goodchild4 则考虑堆场翻箱，得出集卡预约可以降低堆场翻箱率。为提高集卡预约系统效率，一些学者针对预约系统的设计与优化问题开展了研究。如，Huynh 和 Walton5 通过仿真集卡在集装箱码头的作业过程，以码头和集卡车队效益最大为目标，通过设定实施目标，利用数学模型获得每个时间段091可接受的最大集卡数量。由于集卡到达的随机性，刻画集卡到达闸口的过程较为复杂，一些学者针对码头集卡服务排队问题开展了研究，Guan 和 Liu6 利用 M/EK/c 排队模型计算集卡的排队的平均队长，在此基础上构建了预约优化模型。集卡预约系统属于随机服务系统，准确把握集卡排队规律是预约决策的基础，由于集卡到达具有随时间变化的非稳态特征，传统的稳态排队模型虽然可得到集卡排队的平均情况，但不能准确反映集卡的到达规律，为解决此问题，非稳态排队论逐渐得到应用，如 Chen7 建立集卡预约优化模型，利用非稳态排队模型计算集卡在闸口的等待时间，以此为依据优化集卡预约份额。这些研究为集卡预约实践提供了理论基础，但仍然存在以下问题:1)集卡预约系统没有将失约现象作为重要影响因素;2)失约现象下集卡在闸口和堆场等待时间如何计算;3)能否调整闸口和堆场效率来降低失约带来的损耗。基于上述问题，本文以揭示预约机制下的集卡到达与服务规律、考虑集卡失约的情况最小化集卡等待时间为目标，运用马尔可夫、排队论和干扰管理等理论，将失约作为主要的影响因素，以集卡等待时间最小化为目标，利用 python 进行仿真进行深入探讨并详细的数值分析。2模型建立2.1问题描述与基本假设外集卡预约系统旨在降低闸口和堆场拥堵问题，提高港口工作效率，充分利用场桥等资源。预约系统的执行在于集卡公司和码头运营者进行良好沟通，确定每一个时间段集卡到达数，也就是每个时间段的预约份额。预约份额的确定一方面在于货量的多少，另一方面在于码头有多少资源来提供服务。本文不考虑货量，只考虑闸口效率和堆场效率以及场桥的利用率。本文中如无特别指定，集卡统一指外集卡。本文研究中，外集卡港内工作流程如下所示。如图 1。图 1外集卡港口工作流程2.2参数定义和模型建立B:为预约时段，将一天分成 B 个时段，即 b=1，2，3，b每个时间段长为 1440/BminT:反映集卡到达具有随时间变化的特征，更能准确地计算排队队长，将每天进一步细分为 T 个时间段，用来计算队长的最小时间单位，(Tb)N:一个工作周期内集卡数集合Mb:第 b 时段最优化预约的集卡数目，b=1，2，3，Z:闸口数(Z=1)P:集卡准时按照预约时间到达闸口的概率CZ:每辆集卡闸口处单位时间的等待损耗CY:每辆集卡堆场处单位时间的等待损耗CS:场桥为集卡服务而在系统中逗留的单位时间损耗Y:堆场集合g:闸口服务效率y:堆场服务效率k:为闸口和堆场服务时间服从的爱尔郎分布的阶数8 xi:第 i 个集卡和第 i+1 个集卡的预约时间间隔ti:第 i 个集卡的预约时间，显然有 ti=t1+i1j=1xjwzi:第 i 个集卡没有失约而是按时到达的情况下的闸口处平均等待时间wyi:第 i 个集卡没有失约而是按时到达的情况下的堆场里平均等待时间Nb(ti):b 时间段第 i 个集卡预约时间到达之前系统中的集卡数y:单位时间集卡到达堆场的数量C2a:P 指标比计算堆场等待时间的一个参数(C2a=1)9:计算集卡队长等待成本和场桥在系统中逗留成本的等待系数p(Nb(ti)=j):b 时间段第 i 个集卡预约时间到达之前系统中集卡数为 j 的概率S(n，p)/M/Z，已知有 n 个相互独立的预约的集卡，每个集卡以概率 p(0，1按照预约计划 S(np)到达闸口，且闸口和堆场对集卡提供 Markov 服务对于预约排队模型 S(n，p)/M/Z，预约的集卡以概率 p按时到达而以概率 1p 失约。由于闸口不能预知哪些集卡会按时到达而哪些会失约，因此闸口和场桥始终在系统中，期间不进行休假，直到最后一个集卡的预约时间 tn为止。若时刻 tn时系统仍有集卡未完成任务，则需要服务完所有剩余顾客方可离开系统。集卡按照预约顺序接受服务且允许预约同一时间。由于允许预约同一时间，所以定义一个规则:预约时间早的顾客如果没有失约则先对其服务。失约集卡相对于没有预约的集卡则先对其服务。优化具有失约现象的预约排队系统 S(n，p)/M/Z(n=2，3)的目的是要使所有按时到达集卡在闸口和堆场的等待损耗和场桥逗留损耗最小，目标是(?x)=czpni=1wzi+cyni=1wyi+cst1+n1i=1xi+wzn+wyn+py(1)最终目的是使(1)最小化。由于 t1=0 和 csw1=0 恒成191立，且 cspy为恒量，可以省略。上述函数式(1)可以改写为(?x)=csn1i=1xi+czpn1i=2wzi+(czp+cs)wzn+cyn1i=2wyi+(cy+cs)wyn(2)定义广义相对服务效率?=Cs(czp+cs)，用(czp+cs)除(2)，目标函数变为(?x)=?n1i=1xi+(1?)n1i=2wzi+wzn+cy(czp+cs)n1i=2wyi+(cy+cs)(czp+cs)wyn(3)2.3闸口处等待时间计算由于在考虑失约情况下的集卡预约系统异常复杂，本文中只考虑每个时间段的集卡预约情况，即每天的第 b 个时间段，时间长度为 1440/B min，集卡预约数 n=2，3 的 S(n，p)/M/1 系统。基于 Hassin 和 Mendel10 给出的 wzi，2in，的迭代表达式，为此建立一个 33 阶的概率矩阵。P(x1，x2)=P Nb(t1)=000P Nb(t2)=0 P Nb(t2)=10P Nb(t3)=0 P Nb(t3)=1 P Nb(t3)=2(4)第一行P Nb(t1)=0 00=100(5)第二行 P Nb(t2)=0 P Nb(t2)=1 0=1 pegx1pegx10(6)第三行P Nb(t3)=0 P Nb(t3)=1 P Nb(t3)=2=(1 pegx1)(1 pegx2)+pegx1(1 egx2 pgx2egx2)*(1 pegx1)pegx2+pegx1(egx2(pgx2+1 p)pegx1pegx2(7)闸口处集卡等待时间wzi=i1j=1(jgp Nb(ti)=j)(8)2.4堆场处等待时间计算闸口和堆场的服务时间服从阶数为 k 的爱尔郎分布，作业遵循先到先服务(First come first serve，FCFS)，服从参数 ky指数分布。根据吴登磊9 P 指标比法变异系数平方:C2S1=C2S2=1k闸口服务时间:S1=1g堆场服务时间:S2=1y堆场平均速率:yi=ni=1Nb(ti)T定义堆场相对利用率:2=yiy当服务时间服从 k 阶的爱尔郎分布，其指标比 P 近似和堆场等待时间分别为P 1 S21S221 21 2S1S2()(1 C2S1)S2S1(1+C2