游鑫夢(mèng) 梁承姬 張悅

摘要：為提高進(jìn)口集裝箱提取作業(yè)效率，制定合理的堆存策略和翻箱策略，考慮集裝箱堆存作業(yè)與提取作業(yè)的關(guān)聯(lián)性建立兩階段規(guī)劃模型。從降低期望翻箱率的角度，優(yōu)化進(jìn)口集裝箱的箱位分配，構(gòu)建第一階段箱區(qū)貝位分配模型，并采用遺傳算法求解?；诘谝浑A段的箱區(qū)貝位分配結(jié)果，在提箱順序已知的情況下，優(yōu)化障礙箱的落箱位，構(gòu)建第二階段提箱優(yōu)化模型，并設(shè)計(jì)啟發(fā)式算法求解。通過算例將該算法與已有算法進(jìn)行對(duì)比，結(jié)果表明設(shè)計(jì)的算法在降低貝內(nèi)二次翻箱率的效果上更顯著，進(jìn)而表明提出的兩階段規(guī)劃模型能有效降低翻箱率，提升進(jìn)口集裝箱提取作業(yè)效率。

關(guān)鍵詞： 進(jìn)口集裝箱; 堆存策略; 翻箱策略; 啟發(fā)式算法

Abstract： In order to enhance the efficiency of import container picking-up operation and formulate reasonable storage and re-handling strategies， a two-stage programming model is established with the consideration of the relationship between storage and picking-up operations of containers. From the point of view of reducing the expected re-handling rate， the position allocation of import containers is optimized， and the bay-in-block allocation model of the first stage is constructed and solved by the genetic algorithm. Based on the bay-in-block allocation results of the first stage and the known picking-up order of containers， the positions of obstacle containers are optimized， and an optimization model of container picking-up operation of the second stage is constructed and solved by the designed heuristic algorithm. The algorithm is compared with the existing algorithms through examples. The results show that the designed algorithm is more effective in reducing the double re-handling rate， which shows that the proposed two-stage programming model can effectively reduce the re-handling rate and improve the efficiency of import container picking-up operation.

Key words： import container; storage strategy; re-handling strategy; heuristic algorithm

0 引 言

近年來，我國(guó)水路運(yùn)輸業(yè)發(fā)展迅速，碼頭吞吐量急速上升，導(dǎo)致堆場(chǎng)堆放的集裝箱越來越多，堆場(chǎng)翻箱率居高不下。進(jìn)口集裝箱（以下簡(jiǎn)稱進(jìn)口箱）進(jìn)場(chǎng)堆存和提取過程均涉及翻箱問題，而翻箱作業(yè)制約了集裝箱作業(yè)效率，因此本文目的是制定合理的進(jìn)口箱堆存和翻箱策略，從降低翻箱率的角度提升進(jìn)口箱提取效率。

對(duì)于集裝箱堆存過程中的箱位分配問題：BAZZAZI等[1]在考慮集裝箱類型的條件下研究了每個(gè)時(shí)段箱區(qū)的集裝箱分配問題，并采用遺傳算法求解;YU等[2]提出非隔離、單周期隔離和多周期隔離3種堆存策略，并對(duì)這3種策略進(jìn)行了對(duì)比分析;周鵬飛等[3]提出通過統(tǒng)計(jì)箱組堆存時(shí)間獲取提箱優(yōu)先級(jí)的方法，對(duì)卸船箱組的箱位分配進(jìn)行優(yōu)化，并設(shè)計(jì)了啟發(fā)式算法進(jìn)行求解;嚴(yán)偉等[4]基于聚類算法設(shè)置出口集裝箱集港堆存規(guī)則，并通過仿真進(jìn)行了驗(yàn)證;周思方等[5]通過預(yù)估箱組提取時(shí)間范圍建立了箱位指派模型;梁承姬等[6]基于網(wǎng)絡(luò)流方法對(duì)箱區(qū)分配進(jìn)行了動(dòng)態(tài)研究，并采用禁忌搜索算法進(jìn)行求解;武慧榮等[7]以海鐵聯(lián)運(yùn)集裝箱碼頭進(jìn)出口集裝箱為研究對(duì)象，構(gòu)建了堆場(chǎng)箱區(qū)分配模型，并采用模擬退火算法進(jìn)行求解。

對(duì)于集裝箱提取過程中的翻箱問題：文獻(xiàn)[2]借鑒KIM[8]提出的計(jì)算期望翻箱量的方法，研究了箱區(qū)貝位堆存箱量和棧高與翻箱量之間的關(guān)系;LEE等[9]提出一種改進(jìn)模型來預(yù)估編組集裝箱的箱位，并設(shè)計(jì)了啟發(fā)式算法進(jìn)行求解;徐亞等[10]提出一種啟發(fā)式算法及其改進(jìn)算法，對(duì)提箱過程中的翻箱問題進(jìn)行了優(yōu)化;PETERING等[11]采用一種新的混合整數(shù)規(guī)劃模型研究了箱區(qū)集裝箱的重定位問題，并采用擴(kuò)展性的啟發(fā)式算法求解;鄭斯斯等[12]通過對(duì)多種翻箱規(guī)則下的優(yōu)先級(jí)進(jìn)行排序建立了倒箱路徑優(yōu)化模型，并采用啟發(fā)式算法求解;郭瑞智等[13]考慮提箱過程中場(chǎng)橋作業(yè)時(shí)間構(gòu)建了整數(shù)規(guī)劃模型，并設(shè)計(jì)了啟發(fā)式算法求解。

本文在前人研究的基礎(chǔ)上，從降低堆場(chǎng)翻箱率的角度對(duì)進(jìn)口箱的堆存策略和翻箱策略進(jìn)行研究。首先，以場(chǎng)橋期望提箱完工時(shí)間最短為優(yōu)化目標(biāo)，建立第一階段箱區(qū)貝位分配模型。其次，基于第一階段進(jìn)口箱分配優(yōu)化結(jié)果，以翻箱總次數(shù)最少為目標(biāo)，構(gòu)建第二階段提箱優(yōu)化模型。通過設(shè)計(jì)啟發(fā)式算法減少對(duì)障礙箱的二次操作。最后，借助算例將啟發(fā)式算法求解的結(jié)果與已有算法求解的結(jié)果進(jìn)行對(duì)比，驗(yàn)證本文提出算法的有效性。

1 問題描述

進(jìn)口箱從待卸船舶被卸載，并被送入堆場(chǎng)進(jìn)行堆存。根據(jù)文獻(xiàn)[8]的研究結(jié)論，箱區(qū)貝位堆存箱量和棧高與翻箱量之間存在定量關(guān)系，而且可以預(yù)估進(jìn)口箱在提取過程中的期望翻箱量。因此，進(jìn)口箱的堆存優(yōu)化需要在箱區(qū)貝位合理分配集裝箱堆存量，降低箱區(qū)首次翻箱率。

在上述研究的基礎(chǔ)上獲取貝內(nèi)集裝箱的堆存量。將每個(gè)階段需取走的集裝箱記作目標(biāo)箱，將目標(biāo)箱上方的集裝箱記作障礙箱，則優(yōu)化障礙箱的落箱位置是降低貝內(nèi)二次翻箱率的關(guān)鍵。

針對(duì)障礙箱的落箱位置選擇問題，設(shè)計(jì)提箱作業(yè)啟發(fā)式算法。設(shè)障礙箱的提箱優(yōu)先級(jí)為p0，障礙箱落箱候選堆棧的優(yōu)先級(jí)選擇順序如下：（1）若某堆棧內(nèi)堆放的集裝箱的最高優(yōu)先級(jí)低于p0，則該堆棧為障礙箱落箱最優(yōu)堆棧;（2）空棧為障礙箱落箱次優(yōu)堆棧;（3）比較每個(gè)堆?，F(xiàn)有集裝箱的最高優(yōu)先級(jí)，選擇集裝箱最高優(yōu)先級(jí)值最低的堆棧為障礙箱落箱末優(yōu)堆棧。若以上每種類型的堆棧存在多個(gè)，則優(yōu)先考慮鄰近堆棧上的空箱位為障礙箱落箱位置。本文中提箱優(yōu)先級(jí)指提箱先后次序，依該次序?yàn)榧b箱編號(hào)。

綜上所述，進(jìn)口箱堆存和提取均涉及翻箱問題，且兩個(gè)過程相互關(guān)聯(lián)。本文從堆存和提取方面進(jìn)行兩階段研究。第一階段以單箱區(qū)為研究對(duì)象，確定在初步降低箱區(qū)期望首次翻箱率的同時(shí)，優(yōu)化箱區(qū)各貝位的進(jìn)口箱量分配。第二階段在箱區(qū)各貝位集裝箱堆存狀態(tài)和提取順序已知的情況下研究貝內(nèi)翻箱問題。通過優(yōu)化障礙箱的落箱位置進(jìn)一步降低貝內(nèi)二次翻箱率。本文建立進(jìn)口箱堆存和提取兩階段規(guī)劃模型來實(shí)現(xiàn)上述過程。

2 進(jìn)口箱堆存和提取兩階段規(guī)劃模型

2.1 第一階段箱區(qū)貝位分配模型

2.1.1 模型假設(shè)

某時(shí)段被分配到箱區(qū)的集裝箱量已知;集裝箱進(jìn)入堆場(chǎng)時(shí)，分配到箱區(qū)各貝位的概率相同;貝位堆棧的最大層數(shù)為H時(shí)，為滿足所有情況下障礙箱的落箱需求，貝內(nèi)需留有H-1個(gè)空箱位;在無懸臂軌道吊的箱區(qū)布局中，只有一臺(tái)場(chǎng)橋進(jìn)行作業(yè)。

2.1.2 模型參數(shù)與變量

已知參數(shù)：Cu表示在某時(shí)段內(nèi)分配到箱區(qū)的卸船箱總量;n表示箱區(qū)的貝位總數(shù);i表示貝位編號(hào)（i=1，2，…，n）;S表示貝位最大堆棧數(shù)（s=1，2，…，S）;H表示貝內(nèi)堆棧最大層數(shù)（h=1，2，…，H）;C表示貝位堆放集裝箱的最大容量;Coi表示貝位i的初始集裝箱堆存量;to表示場(chǎng)橋單位翻箱作業(yè)時(shí)間;tc表示場(chǎng)橋單位提放箱作業(yè)時(shí)間;tm表示場(chǎng)橋單貝位移動(dòng)時(shí)間;R（C）表示當(dāng)貝內(nèi)集裝箱堆存量為C時(shí)，提取所有集裝箱的預(yù)估期望翻箱量。

決策變量：Ci表示從船舶卸載的集裝箱分配到箱區(qū)貝位i的箱量。

2.1.3 貝位期望翻箱量估計(jì)

借鑒文獻(xiàn)[8]估計(jì)貝位期望翻箱量的方法，研究得出箱區(qū)每個(gè)貝位期望翻箱量與堆存箱量和棧高有關(guān)。當(dāng)貝位最大堆棧數(shù)S已確定時(shí)，根據(jù)貝位的堆存箱量C計(jì)算提箱過程中產(chǎn)生的平均翻箱量R（C）。

2.1.4 箱區(qū)貝位分配模型

以場(chǎng)橋期望提箱完工時(shí)間最短為目標(biāo)，建立箱區(qū)貝位分配模型如下：

式（2）為目標(biāo)函數(shù)，表示目標(biāo)是箱區(qū)場(chǎng)橋期望提箱完工時(shí)間（由場(chǎng)橋期望翻箱作業(yè)時(shí)間、移動(dòng)時(shí)間和提放箱作業(yè)時(shí)間組成）最短，第一項(xiàng)中R（Coi+Ci）-RCoi表示貝位i新增的期望翻箱量;式（3）保證分配到各貝位的箱量之和等于卸船進(jìn)口箱總量;式（4）表示當(dāng)堆棧最底層的集裝箱出現(xiàn)翻箱作業(yè)時(shí)，確保貝內(nèi)有充足的空箱位提供，即限制貝內(nèi)集裝箱的分配數(shù)量不超過其最大容量;式（5）是對(duì)決策變量的取整約束。

2.2 第二階段提箱優(yōu)化模型

2.2.1 模型假設(shè)

貝內(nèi)堆存的集裝箱均為20英尺（1英尺≈0.304 8 m）的標(biāo)準(zhǔn)箱;考慮場(chǎng)橋作業(yè)的安全性，翻箱限制在同一貝內(nèi)進(jìn)行;貝內(nèi)集裝箱初始堆存狀態(tài)和提取順序已知;進(jìn)行提箱作業(yè)時(shí)，不考慮新的集裝箱入貝堆存。

2.2.2 模型參數(shù)與變量

已知參數(shù)：根據(jù)上述箱區(qū)貝位分配模型，可得到貝內(nèi)集裝箱堆存量C。在本階段研究中，共同參數(shù)部分與第一階段模型中的一致，其他具體參數(shù)有：集裝箱的提箱優(yōu)先級(jí)p（p=1，2，…，P）、障礙箱的提箱優(yōu)先級(jí)p0、提箱階段l（l=1，2，…，L）、貝內(nèi)集裝箱的初始堆存狀態(tài)Ipsh（表示提箱優(yōu)先級(jí)為p的集裝箱堆存在堆棧s的第h層）。

決策變量：xlpsh，若第l階段提箱優(yōu)先級(jí)為p的集裝箱堆存在堆棧s的第h層，則其值取1，否則為0;ylpshab，若第l階段提箱優(yōu)先級(jí)為p的集裝箱從堆棧s的第h層被翻至堆棧a的第b層，則其值取1，否則為0;zlpsh，若第l階段提箱優(yōu)先級(jí)為p的集裝箱被從堆棧s的第h層提走，則其值取1，否則為0。

2.2.3 提箱優(yōu)化模型

基于上述參數(shù)和決策變量，以貝內(nèi)翻箱總次數(shù)最少為目標(biāo)，建立進(jìn)口箱提箱優(yōu)化模型：

式（7）表示貝內(nèi)集裝箱初始堆存情況與第0階段的堆存情況一致;式（8）表示每次操作只作業(yè)一個(gè)集裝箱，且該集裝箱只能被提取或被翻倒至其他堆棧[14];式（9）表示若某個(gè)階段產(chǎn)生翻箱作業(yè)，則為障礙箱選擇的落箱位必須是空箱位;式（10）等號(hào)右邊第2項(xiàng)和第3項(xiàng)表示提箱過程中翻箱操作對(duì)位置變量xlpsh的影響，第4項(xiàng)表示提箱操作對(duì)xlpsh的影響[15];式（11）表示在堆存集裝箱或翻箱過程中不能出現(xiàn)懸空位置;式（12）確保貝內(nèi)所有集裝箱最終被提取完畢;式（13）表示任一提箱優(yōu)先級(jí)為p的集裝箱在貝內(nèi)均有確定的箱位，或已被提取離場(chǎng)，或仍堆存在貝內(nèi);式（14）表示貝內(nèi)任一箱位最多只能堆放一個(gè)集裝箱;式（15）表示按照既定的提箱順序依次提箱;式（16）表示決策變量的取值范圍。

3 模型求解

根據(jù)模型的特點(diǎn)，對(duì)第一階段箱區(qū)貝位分配模型采用遺傳算法進(jìn)行求解，對(duì)第二階段提箱優(yōu)化模型設(shè)計(jì)啟發(fā)式算法來求解。兩種求解方法均在MATLAB平臺(tái)上操作運(yùn)行。

3.1 遺傳算法主要求解步驟

（1）初始種群。在進(jìn)行編碼時(shí)采用整數(shù)編碼的形式，假設(shè)箱區(qū)貝位數(shù)為n，需要生成n個(gè)貝位的分配量，即染色體可以表示為：C1，C2，C3，…，Cn。其中，基因值C1，C2，C3，…，Cn分別對(duì)應(yīng)箱區(qū)各貝位進(jìn)口箱的箱量分配。圖1為染色體編碼示例。

（2）適應(yīng)度函數(shù)和選擇策略。目標(biāo)為場(chǎng)橋期望提箱完工時(shí)間最短，因此取目標(biāo)函數(shù)的倒數(shù)作為適應(yīng)度函數(shù)。

采取精英保留策略，即在父代種群中選擇適應(yīng)性強(qiáng)的個(gè)體插入子代種群中，從而保證子代中一定存在優(yōu)于上一代的個(gè)體。

（3）交叉。對(duì)每代種群以一定的交叉率pc進(jìn)行染色體交叉。由于染色體設(shè)計(jì)的特殊性，參考文獻(xiàn)[1]中的線性組合交叉策略，線性組合系數(shù)為k，在兩個(gè)個(gè)體P1，g、P2，g之間進(jìn)行算術(shù)交叉，產(chǎn)生新的個(gè)體。計(jì)算方式如下：

圖2為父代（第g代）染色體交叉示例：箱區(qū)設(shè)置10個(gè)貝位，分配到箱區(qū)的進(jìn)口箱總量為100 TEU，P1，g和P2，g為兩條父代染色體。當(dāng)k=0.85時(shí)，子代P1，g+1在1號(hào)貝位的數(shù)值為：0.85×9+0.15×8=8.85。按照此方法計(jì)算其余位置數(shù)值從而得到子代P1，g+1。子代P2，g+1在1號(hào)貝位的數(shù)值為：0.85×8+0.15×9=8.15。采取相同的方式得到子代P2，g+1。為滿足染色體的基因值為正整數(shù)，線性交叉后子代個(gè)體采用奇數(shù)貝位對(duì)應(yīng)的基因值向負(fù)無窮方向取整，偶數(shù)貝位對(duì)應(yīng)的基因值向正無窮方向取整的方式進(jìn)行保留?；蛐迯?fù)后子代P1，g+1在1號(hào)貝位的數(shù)值為8，子代P2，g+1在1號(hào)貝位的數(shù)值為8。若基因總值不等于100，則在超過或小于額定值的染色體中尋找最大或最小的基因值，并進(jìn)行相應(yīng)的刪減或增加，從而使得染色體滿足條件。

（4）變異。對(duì)每代種群以一定的變異率pm進(jìn)行染色體變異。變異方式采用將兩個(gè)父代染色體上相同基因位置的進(jìn)箱量相互置換，同時(shí)重新獲得箱區(qū)各貝位所對(duì)應(yīng)的進(jìn)口箱分配量。

（5）結(jié)束規(guī)則。當(dāng)算法迭代到設(shè)定的最大次數(shù)時(shí)，結(jié)束并輸出結(jié)果。

3.2 提箱作業(yè)啟發(fā)式算法

貝內(nèi)集裝箱在被提取之前均有確定且唯一的優(yōu)先級(jí)別。視一個(gè)優(yōu)先級(jí)別的集裝箱為一個(gè)提箱階段。當(dāng)提箱過程中出現(xiàn)翻箱操作時(shí)，需尋找其余未滿額定層數(shù)的堆棧，作為翻出障礙箱的可落堆棧集。為障礙箱選擇合適的落箱堆棧是降低二次翻箱率的關(guān)鍵，而可落堆棧集中包含以下6種堆存情形：（1）僅有一個(gè)堆棧內(nèi)現(xiàn)有集裝箱最高優(yōu)先級(jí)低于p0，且無空棧;（2）與情形1的類似，但有多個(gè)滿足條件的合適堆棧;（3）存在某堆棧內(nèi)現(xiàn)有集裝箱最高優(yōu)先級(jí)低于p0的情況，也存在空棧;（4）存在空棧且唯一，而不存在某堆棧內(nèi)現(xiàn)有集裝箱最高優(yōu)先級(jí)低于p0的情況;（5）存在多個(gè)空棧，而不存在某堆棧內(nèi)現(xiàn)有集裝箱最高優(yōu)先級(jí)低于p0的情況;（6）除以上5種情形之外的其他堆存情況。

確定貝內(nèi)每個(gè)提箱階段目標(biāo)箱所存箱位，將可落堆棧集中包含的6種堆存情形設(shè)置成相應(yīng)的落箱規(guī)則，并將其嵌套至提箱作業(yè)過程中;外集卡依次到港，提取已經(jīng)獲得優(yōu)先級(jí)的集裝箱。嵌套落箱規(guī)則的提箱作業(yè)啟發(fā)式算法流程見圖3。

4 算 例

4.1 第一階段箱區(qū)貝位分配模型求解

以自動(dòng)化碼頭單箱區(qū)為研究對(duì)象，參數(shù)設(shè)置見表1。箱區(qū)每個(gè)貝位初始集裝箱堆存量已知，見表2。遺傳算法中：最大迭代次數(shù)設(shè)為1 000;種群規(guī)模設(shè)為100（100條染色體）;交叉系數(shù)pc設(shè)為0.85;變異系數(shù)pm設(shè)為0.15;交叉線性組合系數(shù)k設(shè)為0.8。應(yīng)用MATLAB R2018a編程環(huán)境，使用遺傳算法計(jì)算10次，取最優(yōu)的收斂效果，見圖4。當(dāng)算法迭代至第550次時(shí)，目標(biāo)函數(shù)值趨向收斂，用時(shí)約54 638 s，場(chǎng)橋期望提箱完工時(shí)間達(dá)到最短。當(dāng)目標(biāo)函數(shù)值最小時(shí)，對(duì)應(yīng)決策變量的最優(yōu)結(jié)果見表2，例如，1號(hào)貝位所分配的進(jìn)口箱最優(yōu)箱量為16 TEU，總堆存量最優(yōu)結(jié)果為20 TEU。

4.2 第二階段提箱優(yōu)化模型求解

在提箱優(yōu)化問題上，為驗(yàn)證所設(shè)計(jì)的啟發(fā)式算法的有效性，采用MATLAB R2018a編程將本文提出的啟發(fā)式算法與KH算法[16]、IH算法[10]和OH算法[14]的運(yùn)行結(jié)果進(jìn)行對(duì)比。選取10種不同的貝位規(guī)模，并在每種規(guī)模下隨機(jī)生成100個(gè)算例實(shí)驗(yàn)，最后以平均翻箱量和平均運(yùn)行時(shí)間這兩個(gè)指標(biāo)評(píng)比算法的性能，見表3。第二階段進(jìn)一步落實(shí)到在提箱過程中對(duì)貝內(nèi)集裝箱的翻箱問題，基于第一階段箱區(qū)貝位分配模型的求解結(jié)果，取其中優(yōu)化的進(jìn)箱量為28 TEU的貝位進(jìn)行研究，應(yīng)用翻箱策略提取集裝箱的示例見圖5。集裝箱上的編號(hào)表示提箱優(yōu)先級(jí)，0表示空箱位。當(dāng)場(chǎng)橋提取1號(hào)集裝箱時(shí)，存在壓箱現(xiàn)象，需依次將24號(hào)和22號(hào)箱翻倒至同貝的其他堆棧。根據(jù)提箱作業(yè)啟發(fā)式算法，結(jié)合當(dāng)前階段貝內(nèi)集裝箱的堆存狀態(tài)，直接將24號(hào)和22號(hào)集裝箱分別移至堆棧5和3的頂部。

當(dāng)目標(biāo)箱上的障礙箱全部翻倒后，場(chǎng)橋可提取1號(hào)箱離場(chǎng)。同理，應(yīng)用設(shè)計(jì)的啟發(fā)式算法，依次提取每個(gè)提箱階段的目標(biāo)箱，最終將貝內(nèi)集裝箱提取完畢。

4.3 結(jié)果分析

第一階段優(yōu)化以單箱區(qū)為研究對(duì)象，目標(biāo)是降低箱區(qū)首次翻箱率。當(dāng)場(chǎng)橋的期望提箱完工時(shí)間最短時(shí)，從船舶卸載的集裝箱分配到箱區(qū)的箱量在每個(gè)貝位內(nèi)大致呈現(xiàn)均衡性。而由于場(chǎng)橋的單位翻箱作業(yè)時(shí)間遠(yuǎn)大于其在單位貝位內(nèi)的移動(dòng)時(shí)間，故箱區(qū)內(nèi)鄰近陸側(cè)貝位的進(jìn)箱量比海側(cè)端口的多。

第一階段僅僅從研究堆存策略的角度，初步使堆場(chǎng)箱區(qū)首次翻箱率得到降低。在進(jìn)行第二階段提箱優(yōu)化后，通過合理設(shè)計(jì)障礙箱倒箱落位的啟發(fā)式算法，并結(jié)合算例表明，即使擴(kuò)大貝位規(guī)模，本文設(shè)計(jì)的啟發(fā)式算法在翻箱量和程序運(yùn)行時(shí)間上也優(yōu)于已有算法，驗(yàn)證了本文算法的有效性。因此，可將設(shè)計(jì)的算法應(yīng)用于提箱過程中，進(jìn)一步降低堆場(chǎng)二次翻箱率。

5 結(jié) 論

本文為制定合理的堆存與翻箱策略，對(duì)進(jìn)口箱堆存與提取過程進(jìn)行了優(yōu)化研究，并考慮兩個(gè)過程的關(guān)聯(lián)性，構(gòu)建了兩階段規(guī)劃模型。第一階段以場(chǎng)橋期望提箱完工時(shí)間最短為目標(biāo)，在初步降低箱區(qū)期望首次翻箱率的情況下，建立箱區(qū)貝位分配模型進(jìn)行堆存優(yōu)化。將第一階段求解的貝位進(jìn)箱分配結(jié)果作為輸入，并在集裝箱堆存狀態(tài)和提箱順序已知的條件下，優(yōu)化障礙箱的落箱位，建立第二階段提箱優(yōu)化模型，并設(shè)計(jì)啟發(fā)式算法進(jìn)行求解。該算法與已有算法的對(duì)比結(jié)果表明，應(yīng)用該算法降低貝內(nèi)二次翻箱率的效果更顯著，從而表明提出的兩階段模型能夠降低堆場(chǎng)翻箱率，使進(jìn)口箱提箱作業(yè)效率得到提升。

參考文獻(xiàn)：

[1]BAZZAZI M， SAFAEI N， JAVADIAN N， et al. A genetic algorithm to solve the storage space allocation problem in a container terminal[J]. Computers & Industrial Engineering， 2009， 56： 44-52. DOI： 10.1016/j.cie.2008.03.012.

[2]YU Mingzhu， QI Xiangtong. Storage space allocation models for inbound containers in an automatic container terminal[J]. European Journal of Operational Research， 2013， 226： 32-45. DOI： 10.1016/j.ejor.2012.10.045.

[3]周鵬飛， 李丕安. 集裝箱堆場(chǎng)不確定提箱次序與卸船箱位分配[J]. 哈爾濱工程大學(xué)學(xué)報(bào)， 2013， 34（9）： 1119-1123. DOI： 10.3969 /j.issn.1006-7043.201301022.

[4]嚴(yán)偉， 朱夷詩(shī)， 黃有方， 等. 基于聚類分析的集裝箱碼頭堆場(chǎng)策略[J]. 上海海事大學(xué)學(xué)報(bào)， 2014， 35（1）： 35-40， 59. DOI： 10.13340 /j.jsmu.2014.01.008.

[5]周思方， 張慶年. 基于提箱同步的進(jìn)口箱堆存策略研究[J]. 交通運(yùn)輸系統(tǒng)工程與信息， 2018， 18（5）： 151-157. DOI： 10.16097/j.cnki.1009-6744.2018.05.022.

[6]梁承姬， 賈茹， 盛揚(yáng). 基于網(wǎng)絡(luò)流的自動(dòng)化集裝箱碼頭堆場(chǎng)空間分配[J]. 計(jì)算機(jī)應(yīng)用與軟件， 2018， 35（1）： 77-84. DOI： 10.3969 /j.issn.1000-386x.2018.01.013.

[7]武慧榮， 朱曉寧， 鄧紅星. 集裝箱海鐵聯(lián)運(yùn)港口混堆堆場(chǎng)箱區(qū)均衡分配模型[J]. 重慶交通大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2018， 37（4）： 109-115. DOI： 10.3969 /j.issn.1674-0696.2018.04.17.

[8]KIM K H. Evaluation of the number of rehandles in container yards[J]. Computers & Industrial Engineering， 1997， 32（4）： 701-711.

[9]LEE Yusin， CHAO Shih-Liang. A neighborhood search heuristic for pre-marshalling export containers[J]. European Journal of Operational Research， 2009， 196： 468-475. DOI： 10.1016/j.ejor.2008.03.011.

[10]徐亞， 陳秋雙， 龍磊， 等. 集裝箱倒箱問題的啟發(fā)式算法研究[J]. 系統(tǒng)仿真學(xué)報(bào)， 2008， 20（14）： 3666-3669， 3674. DOI： 10.16182/j.cnki.joss.2008.14.006.

[11]PETERING M E H， HUSSEIN M I. A new mixed integer program and extended look-ahead heuristic algorithm for the block relocation problem[J]. European Journal of Operational Research， 2013， 231： 120-130. DOI： 10.1016/j.ejor.2013.05.037.

[12]鄭斯斯， 王愛虎. 路徑優(yōu)化算法求解集裝箱碼頭堆場(chǎng)翻箱問題[J]. 工業(yè)工程與管理， 2017， 22（3）： 31-40. DOI： 10.19495/j.cnki.1007-5429.2017.03.005.

[13]郭瑞智， 史瑪君， 林昊堃. 集裝箱倒箱問題的模型與啟發(fā)式算法研究[J]. 數(shù)學(xué)雜志， 2017， 37（4）： 805-810. DOI： 10.13548/j.sxzz.2017.04.005.

[14]KIM K H， HONG G-P. A heuristic rule for relocating blocks[J]. Computers & Operations Research， 2006， 33（4）： 940-954. DOI： 10.1016/j.cor.2004.08.005.

[15]朱明華， 程奐翀， 范秀敏. 基于定向搜索算法的集裝箱堆場(chǎng)翻箱問題[J]. 計(jì)算機(jī)集成制造系統(tǒng)， 2012， 18（3）： 639-644. DOI： 10.13196/j.cims.2012.03.193.zhumh.013.

[16]KIM K H， KIM H B. Segregating space allocation models for container inventories in port container terminals[J]. International Journal of Production Economics， 1999， 59： 415-423.

（編輯 賈裙平）