CN109871571B - 一种码头资源配置及船舶队长精准计算方法 - Google Patents
一种码头资源配置及船舶队长精准计算方法 Download PDFInfo
- Publication number
- CN109871571B CN109871571B CN201811612287.8A CN201811612287A CN109871571B CN 109871571 B CN109871571 B CN 109871571B CN 201811612287 A CN201811612287 A CN 201811612287A CN 109871571 B CN109871571 B CN 109871571B
- Authority
- CN
- China
- Prior art keywords
- wharf
- ship
- cost
- berth
- captain
- 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.)
- Active
Links
- 238000013468 resource allocation Methods 0.000 title claims abstract description 38
- 238000004364 calculation method Methods 0.000 title claims abstract description 22
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 20
- 238000012937 correction Methods 0.000 claims abstract description 20
- 238000004088 simulation Methods 0.000 claims abstract description 18
- 238000000034 method Methods 0.000 claims description 40
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 claims description 9
- 238000005457 optimization Methods 0.000 claims description 6
- 238000011160 research Methods 0.000 claims description 4
- 230000002354 daily effect Effects 0.000 claims 4
- 230000003203 everyday effect Effects 0.000 claims 1
- 238000011161 development Methods 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 8
- 230000000694 effects Effects 0.000 description 6
- 238000013459 approach Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 230000002860 competitive effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
Images
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
- Ship Loading And Unloading (AREA)
Abstract
一种码头资源配置及船舶队长精准计算方法,包括如下步骤:确定码头的基本信息,包括历史日均船舶到达数量,泊位数量,岸桥数量,岸桥作业时间分布;结合设备成本构建泊位岸桥资源配置模型;设计考虑岸桥移动的马尔科夫链计算码头中船舶队长期望值,通过历史数据设计误差修正函数并通过修正函数修正误差,将修正后的队长期望值引入泊位岸桥资源配置模型;在Matlab软件中开发模拟迭代算法求解资源配置模型;将所求最优解数据指导码头发展。
Description
技术领域
本发明属于码头资源及运营状况评价领域,具体涉及一种码头资源配置及船舶队长精准计算方法。
背景技术
随着世界经济一体化的快速发展,集装箱航运在国际物流中扮演着极为重要的角色。集装箱运输发展至今近30年,除2009年全球经济危机时期外,集装箱吞吐量一直保持稳定增长,其中前20年增长率为10%左右,近10年虽然增长有所减缓,但仍保持5%左右。集装箱运输的快速发展带动了集装箱码头之间的激烈竞争,只有缩短码头作业时间和降低码头集装箱船舶的运营成本,才能有竞争优势。
码头运营前期规划由可预测范围内的每日船舶到达数来确定投入使用的码头泊位、岸桥数量,而随着集装箱吞吐量不断增加,泊位、岸桥资源将会成为码头发展瓶颈。因此,根据集装箱船舶到达码头的时间分布规律,精准系统化地对集装箱码头泊位、岸桥资源重新配置,同时精准计算船舶队长以评价码头运营状况显得非常重要。
发明内容
为了解决集装箱码头泊位、岸桥资源配置及船舶队长等运营指标计算不够精准的问题,本发明提供了一种码头资源配置及船舶队长精准计算方法。
本发明通过以下技术方案实现:
一种码头资源配置及船舶队长精准计算方法,包括如下步骤:
步骤1:以码头设备的数量优化作为研究对象,结合历史日均船舶到达数量,泊位数量,岸桥数量,岸桥作业时间分布;将设备资源通过折旧法得到每日均摊成本,同时根据船舶在泊成本,结合设备成本构建码头的泊位和岸桥资源配置模型;
步骤2:根据码头船舶及泊位尺寸确定一艘船舶平均可同时进行作业的岸桥数量,并构建马尔科夫链得到码头稳态时的船舶队长;
步骤3:将步骤2中得到的船舶队长与历史真实数据对比,设计队长修正函数,并通过修正函数修正计算得到的船舶队长;
步骤4:将修正后的船舶队长公式带入码头泊位和岸桥资源配置模型;
步骤5:设计模拟迭代算法求解资源配置模型;
步骤6:将求解结果应用于码头实际,并指导码头设备数量优化。
进一步地,在上述技术方案中,所述所述的成本考虑四方面的运营成本,包括泊位固定投资成本、岸桥固定投资成本、货物的时间价值及船舶无法及时靠泊产生的惩罚成本,目标函数为:
其中,Cv为船舶在码头的运营费用;Cb为单个泊位每天的运营成本;Cq为单座岸桥每天的运营成本;V为每艘船舶所运载集装箱价值的平均值;r为现金年利率;Em为单位船舶在锚地落锚产生的额外费用;Eg为单位船舶无法及时落锚产生的额外费用;Cp为码头每天的运营成本;L为码头服务系统中船舶数量的平均值;Lq为码头服务系统中船舶等待队长的平均值;Lm为码头服务系统中船舶在锚地等待数量的平均值;Lg为船舶无法驶入锚地落锚数量的平均值;B为泊位的数量;Q为岸桥的数量;CbB表示建设泊位分摊的每日固定投资成本;CqQ表示建设岸桥分摊的每日固定投资成本;表示一艘船舶在码头一天的运营费用及货物积压费用;EmLm表示船舶在锚地产生的额外费用;EgLg表示船舶无法进入码头产生的额外费用;
进一步地,在上述技术方案中,所述的马尔科夫链考虑了岸桥移动对泊位服务能力产生的影响,得到系统服务能力与船舶数满足如下关系:
其中μk表示码头内船舶数为k时系统总服务能力;μ为单座岸桥的装卸能力。
进一步地,在上述技术方案中,所述队长修正函数根据年份使用权重法,其公式为
进一步地,在上述技术方案中,所述模拟迭代算法中的泊位数与岸桥数的比例系数η的取值范围为1<η≤3;
进一步地,在上述技术方案中,所述的马尔科夫链根据岸桥作业时间分布构建,当岸桥作业时间呈poisson分布时,各队长指标近似解如下:
Lg=Lq-Lm (6)
当岸桥作业时间呈Erlang分布时,各队长指标近似解如下:
Lg=Lq-Lm (10)
其中,Pni为码头中有n艘船舶在进行第i道装卸阶段的概率;
0状态:λP0=2μ*P12,λP0+λP12=2μ*P11,P1=P11+P12;
1状态:λP11+λP12=4μ*P22,λP11+λP22=4μ*P21,P2=P21+P22;
…
j状态:λPj1+λPj2=(2jμ*+μ)P(j+1)2,λPj1+λP(j+1)2=(2jμ*+μ)P(j+1)1,Pj+1=P(j+1)1+P(j+1)2;
j+1状态:λP(j+1)1+λP(j+1)2=(2jμ*+2μ)P(j+2)2,λP(j+1)1+λP(j+2)2=(2jμ*+2μ)P(j+2)1,Pj+2=P(j+2)1+P(j+2)2;
…
B状态:λPB1+λPB2=2QμP(B+1)2,λPB1+λP(B+1)2=2QμP(B+1)1,PB+1=P(B+1)1+P(B+1)2;
B+1状态:λP(B+1)1+λP(B+1)2=2QμP(B+2)2,λP(B+1)1+λP(B+2)2=2QμP(B+2)1,PB+2=P(B+2)1+P(B+2)2;
…
进一步地,在上述技术方案中,所述的模拟迭代算法具体包括如下步骤:
步骤(1):根据现有的码头资源配置现状,确定模拟迭代算法的初始泊位数B和初始岸桥数Q;
步骤(2):根据码头岸桥分布规律,将初始泊位数B和初始岸桥数Q等原始数据代入各队长指标公式后得到具体队长数值;
步骤(3):将各队长数值及泊位数B和岸桥数Q代入码头资源配置模型,求得目标函数,即码头日均成本Cp(B,Q);
步骤(4):根据比例系数η的取值范围1<η≤3,将岸桥数Q迭代为岸桥数Q+1,并进入步骤(2),得到新的队长数值,并得到新的日均成本Cp(B,Q+1);当岸桥数Q迭代超过1<η≤3时,将泊位数B迭代为泊位数B+1,同时岸桥数更新为初始岸桥数Q,进入步骤(2),得到新的队长数值并得到新的日均成本Cp(B+1,Q);当Cp(B,Q-1)<Cp(B,Q)<Cp(B,Q+1)时,得到码头资源配置最优解,即泊位数B,岸桥数Q。
本发明的有益效果为:
(1)本发明所述码头资源配置及船舶队长精准计算方法,提出了考虑岸桥移动对马尔科夫链的影响,相比传统的经典排队论模型计算队长,其所求的队长值更接近码头实际情况,而后通过队长修正函数,使所求解的码头运营指标更为精准。
(2)本发明所述码头资源配置及船舶队长精准计算方法,构建的码头资源配置模型考虑了码头实际作业中所产生的多种成本,使模型所反映的情况更趋近码头实际运营成本;设计的模拟迭代算法结合码头泊位与岸桥的相互关系,减少了迭代次数,在计算大规模码头时能较为明显得降低计算时间;最后求解得到的码头泊位和岸桥资源配置数量能够有效指导货运量增加后的码头发展规划。
附图说明
下面结合附图和具体实施方法对本发明作进一步详细的说明。
图1为码头资源配置及船舶队长精准计算方法流程示意图;
图2为码头内各设备运营示意图;
图3为岸桥作业时间服从Poisson分布时的马尔科夫链图;
图4为岸桥作业时间服从Erlang分布时的马尔科夫链图;
图5为模拟迭代算法示意图;
图6(a)~图6(b)为本发明计算方法与经典排队论方法的优劣势对比图,其中图6(a)是实施例1的对比图,图6(b)是实施例2的对比图;
图7(a)~图7(b)为本发明计算方法更为精准的原因图,其中图7(a)是码头服务能力与码头内船舶数量关系图,图7(b)是码头中船舶数量的概率图;
图8(a)~图8(b)为本发明计算方法指导下的码头运营指标图,其中图8(a)是本发明计算方法优化后的每年的码头队长指标图,图8(b)是每年的船舶及时靠泊率图。
具体实施方法
下面结合附图和实施案例对本发明的技术方案进行完整的描述。
实施例1
如图1所示,一种码头资源配置及船舶队长精准计算方法,包括如下步骤:步骤1:以码头设备的数量优化作为研究对象,如图2所示,结合历史日均船舶到达数量λ,泊位数量B,岸桥数量Q,岸桥作业时间分布μ;将设备资源通过折旧法得到每日均摊成本,同时根据船舶在泊成本,结合设备成本构建码头的泊位和岸桥资源配置模型;
步骤2:根据码头船舶及泊位尺寸确定一艘船舶平均可同时进行作业的岸桥数量为3个,并构建马尔科夫链得到码头稳态时的船舶队长L,Lq,Lm,Lg;
步骤3:将步骤2中得到的船舶队长与历史真实数据对比,设计队长修正函数,并通过修正函数修正计算得到的船舶队长;
步骤4:将修正后的船舶队长公式带入码头泊位和岸桥资源配置模型;
步骤5:设计模拟迭代算法求解资源配置模型;
步骤6:将求解结果应用于码头实际,并指导码头设备数量优化。
进一步地,在上述技术方案中,所述所述的成本考虑四方面的运营成本,包括泊位固定投资成本、岸桥固定投资成本、货物的时间价值及船舶无法及时靠泊产生的惩罚成本,目标函数为:
其中,Cv为船舶在码头的运营费用;Cb为单个泊位每天的运营成本;Cq为单座岸桥每天的运营成本;V为每艘船舶所运载集装箱价值的平均值;r为现金年利率;Em为单位船舶在锚地落锚产生的额外费用;Eg为单位船舶无法及时落锚产生的额外费用;Cp为码头每天的运营成本;L为码头服务系统中船舶数量的平均值;Lq为码头服务系统中船舶等待队长的平均值;Lm为码头服务系统中船舶在锚地等待数量的平均值;Lg为船舶无法驶入锚地落锚数量的平均值;B为泊位的数量;Q为岸桥的数量;CbB表示建设泊位分摊的每日固定投资成本;CqQ表示建设岸桥分摊的每日固定投资成本;表示一艘船舶在码头一天的运营费用及货物积压费用;EmLm表示船舶在锚地产生的额外费用;EgLg表示船舶无法进入码头产生的额外费用;
进一步地,在上述技术方案中,所述的马尔科夫链考虑了岸桥移动对泊位服务能力产生的影响,如图3所示,得到系统服务能力与船舶数满足如下关系:
其中μk表示码头内船舶数为k时系统总服务能力;μ为单座岸桥的装卸能力。进一步地,在上述技术方案中,所述队长修正函数根据年份使用权重法,其公式为
进一步地,在上述技术方案中,所述模拟迭代算法中的泊位数与岸桥数的比例系数η的取值范围为1<η≤3;
进一步地,在上述技术方案中,所述的马尔科夫链根据岸桥作业时间分布构建,当岸桥作业时间呈poisson分布时,各队长指标近似解如下:
Lg=Lq-Lm (6)
进一步地,在上述技术方案中,所述的模拟迭代算法如图5所示,具体包括如下步骤:
步骤(1):根据现有的码头资源配置现状,确定模拟迭代算法的初始泊位数B和初始岸桥数Q;
步骤(2):根据码头岸桥分布规律,将初始泊位数B和初始岸桥数Q等原始数据代入各队长指标公式后得到具体队长数值;
步骤(3):将各队长数值及泊位数B和岸桥数Q代入码头资源配置模型,求得目标函数,即码头日均成本Cp(B,Q);
步骤(4):根据比例系数η的取值范围1<η≤3,将岸桥数Q迭代为岸桥数Q+1,并进入步骤(2),得到新的队长数值,并得到新的日均成本Cp(B,Q+1);当岸桥数Q迭代超过1<η≤3时,将泊位数B迭代为泊位数B+1,同时岸桥数更新为初始岸桥数Q,进入步骤(2),得到新的队长数值并得到新的日均成本Cp(B+1,Q);当Cp(B,Q-1)<Cp(B,Q)<Cp(B,Q+1)时,得到码头资源配置最优解,即泊位数B,岸桥数Q。
本例1规定码头内的船舶遵循FCFS规则;且岸桥或其他设备不出现故障。
设定码头的配置为岸桥数量Q=13,泊位数量B=5,μ=1.6,其余参数Cv=32万元/天·艘,V=57600万元/艘,r=15%,Cb=6.1万元/天,Cq=1.1万元/天,Em=0.3万元/天·艘,Eg=0.5万元/天·艘。
船舶到达数量数据如表1所示。
表1船舶到达数量数据
将表1数据通过本申请方法计算得到日均船舶队长为Ls,其与历史数据显示的实际队长的偏差如图6(b)所示;在图6(b)中显示了通过传统计算方法的队长与实际值的偏差。通过表1的原始参数及图6(b)的结果,可以看出本发明的考虑岸桥移动的马尔科夫链计算船舶队长方法有着良好的效果,相比传统经典排队论方法,其解的效果能提高10%左右,自身与实际偏差仅为1%左右,因此,本申请设计的求解算法是有效的。
通过指数平滑法预测未来船舶到达数据后,根据本申请设计的码头资源配置方法,得到最优解如表2所示。
表2资源优化表
由表2可知,本申请的方法能够有效合理配置码头资源数量,并降低码头日均成本。证明了本申请方法的有效性。
实施例2
如图1所示,一种码头资源配置及船舶队长精准计算方法,包括如下步骤:
步骤1:以码头设备的数量优化作为研究对象,如图2所示,结合历史日均船舶到达数量λ,泊位数量B,岸桥数量Q,岸桥作业时间分布μ;将设备资源通过折旧法得到每日均摊成本,同时根据船舶在泊成本,结合设备成本构建码头的泊位和岸桥资源配置模型;
步骤2:根据码头船舶及泊位尺寸确定一艘船舶平均可同时进行作业的岸桥数量为3个,并构建马尔科夫链得到码头稳态时的船舶队长L,Lq,Lm,Lg;
步骤3:将步骤2中得到的船舶队长与历史真实数据对比,设计队长修正函数,并通过修正函数修正计算得到的船舶队长;
步骤4:将修正后的船舶队长公式带入码头泊位和岸桥资源配置模型;
步骤5:设计模拟迭代算法求解资源配置模型;
步骤6:将求解结果应用于码头实际,并指导码头设备数量优化。
进一步地,在上述技术方案中,所述所述的成本考虑四方面的运营成本,包括泊位固定投资成本、岸桥固定投资成本、货物的时间价值及船舶无法及时靠泊产生的惩罚成本,目标函数为:
其中,Cv为船舶在码头的运营费用;Cb为单个泊位每天的运营成本;Cq为单座岸桥每天的运营成本;V为每艘船舶所运载集装箱价值的平均值;r为现金年利率;Em为单位船舶在锚地落锚产生的额外费用;Eg为单位船舶无法及时落锚产生的额外费用;Cp为码头每天的运营成本;L为码头服务系统中船舶数量的平均值;Lq为码头服务系统中船舶等待队长的平均值;Lm为码头服务系统中船舶在锚地等待数量的平均值;Lg为船舶无法驶入锚地落锚数量的平均值;B为泊位的数量;Q为岸桥的数量;CbB表示建设泊位分摊的每日固定投资成本;CqQ表示建设岸桥分摊的每日固定投资成本;表示一艘船舶在码头一天的运营费用及货物积压费用;EmLm表示船舶在锚地产生的额外费用;EgLg表示船舶无法进入码头产生的额外费用。
进一步地,在上述技术方案中,所述的马尔科夫链考虑了岸桥移动对泊位服务能力产生的影响,如图4所示,得到系统服务能力与船舶数满足如下关系:
其中μk表示码头内船舶数为k时系统总服务能力;μ为单座岸桥的装卸能力。
进一步地,在上述技术方案中,所述队长修正函数根据年份使用权重法,其公式为
进一步地,在上述技术方案中,所述模拟迭代算法中的泊位数与岸桥数的比例系数η的取值范围为1<η≤3;
进一步地,在上述技术方案中,所述的马尔科夫链根据岸桥作业时间分布构建,当岸桥作业时间呈Erlang分布时,各队长指标近似解如下:
Lg=Lq-Lm(10)
其中,Pni为码头中有n艘船舶在进行第i道装卸阶段的概率;
0状态:λP0=2μ*P12,λP0+λP12=2μ*P11,P1=P11+P12;
1状态:λP11+λP12=4μ*P22,λP11+λP22=4μ*P21,P2=P21+P22;
…
j状态:λPj1+λPj2=(2jμ*+μ)P(j+1)2,λPj1+λP(j+1)2=(2jμ*+μ)P(j+1)1,Pj+1=P(j+1)1+P(j+1)2;
j+1状态:λP(j+1)1+λP(j+1)2=(2jμ*+2μ)P(j+2)2,λP(j+1)1+λP(j+2)2=(2jμ*+2μ)P(j+2)1,Pj+2=P(j+2)1+P(j+2)2;
…
B状态:λPB1+λPB2=2QμP(B+1)2,λPB1+λP(B+1)2=2QμP(B+1)1,PB+1=P(B+1)1+P(B+1)2;
B+1状态:λP(B+1)1+λP(B+1)2=2QμP(B+2)2,λP(B+1)1+λP(B+2)2=2QμP(B+2)1,PB+2=P(B+2)1+P(B+2)2;
…
进一步地,在上述技术方案中,所述的模拟迭代算法如图5所示,具体包括如下步骤:
步骤(1):根据现有的码头资源配置现状,确定模拟迭代算法的初始泊位数B和初始岸桥数Q;
步骤(2):根据码头岸桥分布规律,将初始泊位数B和初始岸桥数Q等原始数据代入各队长指标公式后得到具体队长数值;
步骤(3):将各队长数值及泊位数B和岸桥数Q代入码头资源配置模型,求得目标函数,即码头日均成本Cp(B,Q);
步骤(4):根据比例系数η的取值范围1<η≤3,将岸桥数Q迭代为岸桥数Q+1,并进入步骤(2),得到新的队长数值,并得到新的日均成本Cp(B,Q+1);当岸桥数Q迭代超过1<η≤3时,将泊位数B迭代为泊位数B+1,同时岸桥数更新为初始岸桥数Q,进入步骤(2),得到新的队长数值并得到新的日均成本Cp(B+1,Q);当Cp(B,Q-1)<Cp(B,Q)<Cp(B,Q+1)时,得到码头资源配置最优解,即泊位数B,岸桥数Q。
本例2规定码头内的船舶遵循FCFS规则;且岸桥或其他设备不出现故障。
设定码头的配置为岸桥数量Q=11,泊位数量B=4,μ=1.33,其余参数Cv=32万元/天·艘,V=57600万元/艘,r=15%,Cb=6.1万元/天,Cq=1.1万元/天,Em=0.3万元/天·艘,Eg=0.5万元/天·艘。
船舶到达数量数据如表3所示。
表3船舶到达数量数据
将表3数据通过本申请方法计算得到日均船舶队长为Ls,其与历史数据显示的实际队长的偏差如图6(a)所示;在图6(a)中显示了通过传统计算方法的队长与实际值的偏差。通过表3的原始参数及图6(a)的结果,可以看出本发明的考虑岸桥移动的马尔科夫链计算船舶队长方法有着良好的效果,相比传统经典排队论方法,其解的效果能提高5%左右,自身与实际偏差仅为0.5%左右,因此,本申请设计的求解算法是有效的。
通过指数平滑法预测未来船到达数据后,根据本申请设计的码头资源配置方法,得到最优解如表4所示。
表4资源优化表
由表4可知,本申请的方法能够有效合理配置码头资源数量,并降低码头日均成本。证明了本申请方法的有效性。
对本申请设计的方法与传统的经典排队论方法对比,得到在上述两例情况下,两种方法的差异图,如图7(a)和图7(b)所示。图7(a)中的TA和TB点表示两种方法差别的极值点。可知,是由于两种方法的规律差异范围及码头服务能力差异范围高度重合引起。
经本申请设计的方法所配置的码头的运营指标与传统方法相比,其优势非常明显,如图8(a)和图8(b)所示,在每一年的各项指标中,本申请设计的方法均较优。
本发明的有益效果为:
(1)本发明所述码头资源配置及船舶队长精准计算方法,提出了考虑岸桥移动对马尔科夫链的影响,相比传统的经典排队论模型计算队长,其所求的队长值更接近码头实际情况,而后通过队长修正函数,使所求解的码头运营指标更为精准。
(2)本发明所述码头资源配置及船舶队长精准计算方法,构建的码头资源配置模型考虑了码头实际作业中所产生的多种成本,使模型所反映的情况更趋近码头实际运营成本;设计的模拟迭代算法结合码头泊位与岸桥的相互关系,减少了迭代次数,在计算大规模码头时能较为明显得降低计算时间;最后求解得到的码头泊位和岸桥资源配置数量能够有效指导货运量增加后的码头发展规划。
以上所述,仅为本发明较佳的具体实施方法,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。
Claims (1)
1.一种码头资源配置及船舶队长精准计算方法,其特征在于:包括如下步骤:
步骤1:以码头设备的数量优化作为研究对象,结合历史日均船舶到达数量,泊位数量,岸桥数量,岸桥作业时间分布;将设备资源通过折旧法得到每日均摊成本,同时根据船舶在泊成本,结合设备成本构建码头的泊位和岸桥资源配置模型,其中成本考虑四方面的运营成本,包括泊位固定投资成本、岸桥固定投资成本、货物的时间价值及船舶无法及时靠泊产生的惩罚成本;
目标函数为:
其中,Cv为船舶在码头的运营费用;Cb为单个泊位每天的运营成本;Cq为单座岸桥每天的运营成本;V为每艘船舶所运载集装箱价值的平均值;r为现金年利率;Em为单位船舶在锚地落锚产生的额外费用;Eg为单位船舶无法及时落锚产生的额外费用;Cp为码头每天的运营成本;L为码头服务系统中船舶数量的平均值;Lq为码头服务系统中船舶等待队长的平均值;Lm为码头服务系统中船舶在锚地等待数量的平均值;Lg为船舶无法驶入锚地落锚数量的平均值;B为泊位的数量;Q为岸桥的数量;CbB表示建设泊位分摊的每日固定投资成本;CqQ表示建设岸桥分摊的每日固定投资成本;表示一艘船舶在码头一天的运营费用及货物积压费用;EmLm表示船舶在锚地产生的额外费用;EgLg表示船舶无法进入码头产生的额外费用;
步骤2:根据码头船舶及泊位尺寸确定一艘船舶平均可同时进行作业的岸桥数量,并构建马尔科夫链得到码头稳态时的船舶队长,马尔科夫链中考虑了岸桥移动对泊位服务能力产生的影响,系统服务能力与船舶数满足如下关系:
其中μk表示码头内船舶数为k时系统总服务能力;μ为单座岸桥的装卸能力;
马尔科夫链根据岸桥作业时间分布构建,当岸桥作业时间呈poisson分布时,各队长指标近似解如下:
Lg=Lq-Lm (6)
当岸桥作业时间呈Erlang分布时,各队长指标近似解如下:
Lg=Lq-Lm (10)
其中,Pni为码头中有n艘船舶在进行第i道装卸阶段的概率;
0状态:λP0=2μ*P12,λP0+λP12=2μ*P11,P1=P11+P12;
1状态:λP11+λP12=4μ*P22,λP11+λP22=4μ*P21,P2=P21+P22;
…
j状态:λPj1+λPj2=(2jμ*+μ)P(j+1)2,λPj1+λP(j+1)2=(2jμ*+μ)P(j+1)1,Pj+1=P(j+1)1+P(j+1)2;
j+1状态:λP(j+1)1+λP(j+1)2=(2jμ*+2μ)P(j+2)2,λP(j+1)1+λP(j+2)2=(2jμ*+2μ)P(j+2)1,Pj+2=P(j+2)1+P(j+2)2;
…
B状态:λPB1+λPB2=2QμP(B+1)2,λPB1+λP(B+1)2=2QμP(B+1)1,PB+1=P(B+1)1+P(B+1)2;
B+1状态:λP(B+1)1+λP(B+1)2=2QμP(B+2)2,λP(B+1)1+λP(B+2)2=2QμP(B+2)1,PB+2=P(B+2)1+P(B+2)2;
…
步骤3:将步骤2中得到的船舶队长与历史真实数据对比,设计队长修正函数,并通过修正函数修正计算得到的船舶队长,队长修正函数根据年份使用权重法,其公式为
步骤4:将修正后的船舶队长公式带入码头泊位和岸桥资源配置模型;
步骤5:设计模拟迭代算法求解资源配置模型,其中模拟迭代算法中的泊位数与岸桥数的比例系数η的取值范围为1<η≤3;
步骤6:将求解结果应用于码头实际,并指导码头设备数量优化。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811612287.8A CN109871571B (zh) | 2018-12-27 | 2018-12-27 | 一种码头资源配置及船舶队长精准计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811612287.8A CN109871571B (zh) | 2018-12-27 | 2018-12-27 | 一种码头资源配置及船舶队长精准计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109871571A CN109871571A (zh) | 2019-06-11 |
CN109871571B true CN109871571B (zh) | 2023-04-07 |
Family
ID=66917235
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811612287.8A Active CN109871571B (zh) | 2018-12-27 | 2018-12-27 | 一种码头资源配置及船舶队长精准计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109871571B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112249583B (zh) * | 2020-09-30 | 2022-10-14 | 华东师范大学 | 一种散货码头堆场垛位的智能选择方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101458788A (zh) * | 2008-12-19 | 2009-06-17 | 天津市市政工程设计研究院 | 集装箱码头物流仿真优化系统 |
CN101789093A (zh) * | 2009-07-21 | 2010-07-28 | 上海海事大学 | 一种集装箱码头泊位与岸桥分配方法 |
CN101986313A (zh) * | 2010-11-03 | 2011-03-16 | 上海海事大学 | 基于知识的集装箱码头泊位与岸桥调度方法 |
CN103606124A (zh) * | 2013-11-25 | 2014-02-26 | 天津市市政工程设计研究院 | 集装箱码头装卸交通系统布局与交通组织微观仿真方法 |
CN108090599A (zh) * | 2016-11-21 | 2018-05-29 | 上海海勃物流软件有限公司 | 集装箱码头运营成本控制和客户价值分析方法 |
-
2018
- 2018-12-27 CN CN201811612287.8A patent/CN109871571B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101458788A (zh) * | 2008-12-19 | 2009-06-17 | 天津市市政工程设计研究院 | 集装箱码头物流仿真优化系统 |
CN101789093A (zh) * | 2009-07-21 | 2010-07-28 | 上海海事大学 | 一种集装箱码头泊位与岸桥分配方法 |
CN101986313A (zh) * | 2010-11-03 | 2011-03-16 | 上海海事大学 | 基于知识的集装箱码头泊位与岸桥调度方法 |
CN103606124A (zh) * | 2013-11-25 | 2014-02-26 | 天津市市政工程设计研究院 | 集装箱码头装卸交通系统布局与交通组织微观仿真方法 |
CN108090599A (zh) * | 2016-11-21 | 2018-05-29 | 上海海勃物流软件有限公司 | 集装箱码头运营成本控制和客户价值分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109871571A (zh) | 2019-06-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Soner et al. | Application of AHP and VIKOR methods under interval type 2 fuzzy environment in maritime transportation | |
He et al. | Modeling berth allocation and quay crane assignment considering QC driver cost and operating efficiency | |
CN101789093A (zh) | 一种集装箱码头泊位与岸桥分配方法 | |
CN106651280B (zh) | 一种集装箱轮船物流运输调度方法及系统 | |
CN102799733B (zh) | 集装箱牵引车资源整合仿真系统及其仿真方法 | |
Wang et al. | Impact of navigational safety level on seaport fairway capacity | |
Alcalde et al. | Optimal space for storage yard considering yard inventory forecasts and terminal performance | |
CN107480830A (zh) | 基于港口‑航道‑锚地复合系统仿真的锚位数优化方法 | |
CN109871571B (zh) | 一种码头资源配置及船舶队长精准计算方法 | |
CN115330227A (zh) | 集装箱堆场的箱位分配方法、装置、设备及介质 | |
CN113822518A (zh) | 一种ais大数据驱动的集装箱港口装卸效率计算方法 | |
Rouky et al. | A robust metaheuristic for the rail shuttle routing problem with uncertainty: a real case study in the Le Havre Port | |
US20230252395A1 (en) | A quay crane operation method | |
Chen et al. | Monitoring and evaluation of ship operation congestion status at container ports based on AIS data | |
Yang et al. | Applications of chaotic quantum adaptive satin bower bird optimizer algorithm in berth-tugboat-quay crane allocation optimization | |
JP5569199B2 (ja) | 配船計画作成装置、方法及びプログラム | |
CN117252509A (zh) | 基于商品库存压力的水运排船方法及装置 | |
CN112700086A (zh) | 一种基于分层约束条件的炼厂泊位调度方法及装置 | |
Maglić et al. | Effectiveness of a barge-based ballast water treatment system for multi-terminal ports | |
Zhang et al. | Optimising feeder routing for container ships through an electronic chart display and information system | |
Zheng et al. | Hybrid dynamic modeling and receding horizon speed optimization for liner shipping operations from schedule reliability and energy efficiency perspectives | |
Ruscă et al. | Risk assessment for dangerous goods in maritime transport | |
Buhari et al. | An assessment of the relationship among cargo-throughput, vessel turnaround time and port-revenue in Nigeria.(A study of Lagos port complex) | |
Park et al. | A Predictive Discrete Event Simulation for Predicting Operation times in Container Terminal | |
JP4498878B2 (ja) | 原料船配船計画装置及び原料船配船計画方法、並びにプログラム |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |