李夢(mèng)娟 任靜

集裝箱在使用過(guò)程中難免遭受不同程度的碰撞、摩擦和污染，需要及時(shí)維修。合理布局集裝箱維修點(diǎn)能夠提高集裝箱周轉(zhuǎn)率，緩解箱源緊張問(wèn)題。目前，我國(guó)尚未形成科學(xué)合理的集裝箱維修點(diǎn)規(guī)劃布局體系，設(shè)立集裝箱維修點(diǎn)時(shí)往往僅考慮貨運(yùn)量因素，無(wú)法有效滿足集裝箱維修需求。業(yè)界對(duì)集裝箱維修點(diǎn)選址的研究較少：張陽(yáng)[1]考慮需求和距離兩大因素，以鐵路集裝箱辦理站與鐵路集裝箱維修點(diǎn)距離最短為目標(biāo)構(gòu)建模型，優(yōu)化鐵路集裝箱維修點(diǎn)布局。從服務(wù)設(shè)施選址領(lǐng)域的研究成果來(lái)看：以往的研究以傳統(tǒng)選址模型為主，考慮的因素和目標(biāo)較少，一般僅限于服務(wù)響應(yīng)速度提升或成本優(yōu)化等[2-4];近年來(lái)的研究以擴(kuò)展選址模型為主，在綜合考慮各種影響因素的基礎(chǔ)上，著力解決服務(wù)半徑與成本優(yōu)化的矛盾，實(shí)現(xiàn)單一優(yōu)化目標(biāo)向整體優(yōu)化目標(biāo)的轉(zhuǎn)變[5-10]。本文在考慮建站成本、運(yùn)輸成本、維修需求、服務(wù)半徑等因素的基礎(chǔ)上，建立以運(yùn)營(yíng)成本最小化為目標(biāo)的集裝箱維修點(diǎn)選址模型，并采用改進(jìn)遺傳算法對(duì)實(shí)例求解，再?gòu)牡乩砦恢?、交通條件、基礎(chǔ)設(shè)施、政策環(huán)境等角度對(duì)求解結(jié)果進(jìn)行評(píng)價(jià)分析，最終得出合理的集裝箱維修點(diǎn)布局方案。

1 集裝箱維修點(diǎn)選址模型

某區(qū)域內(nèi)存在若干修箱需求點(diǎn)和備選集裝箱維修點(diǎn)。破損集裝箱從修箱需求點(diǎn)運(yùn)至集裝箱維修點(diǎn)，經(jīng)維修后恢復(fù)使用功能。為了在提高集裝箱周轉(zhuǎn)率和客戶滿意度的同時(shí)降低集裝箱維修點(diǎn)的運(yùn)營(yíng)成本，本文構(gòu)建模型求解集裝箱維修點(diǎn)的最優(yōu)數(shù)量和位置，以實(shí)現(xiàn)集裝箱維修點(diǎn)運(yùn)營(yíng)成本和服務(wù)半徑的雙重優(yōu)化。

1.1 符號(hào)定義

1.1.1 集合

I為所有修箱需求點(diǎn)的集合，I={1，2，…，i};J為所有備選集裝箱維修點(diǎn)的集合，J={1，2，…，j}。

1.1.2 參數(shù)

f0為集裝箱維修點(diǎn)投入設(shè)備的固定費(fèi)用，元;fj為集裝箱維修點(diǎn)j的土地和建設(shè)費(fèi)用，元;dij為修箱需求點(diǎn)i與集裝箱維修點(diǎn)j之間的距離，km;c0為單個(gè)集裝箱的運(yùn)輸價(jià)格，元/（km 個(gè)）;dmax為集裝箱維修點(diǎn)的最大服務(wù)半徑，km;Qi為修箱需求點(diǎn)i需要維修的集裝箱數(shù)量，個(gè);T為集裝箱維修點(diǎn)的設(shè)備折舊期，年;P為集裝箱維修點(diǎn)為超過(guò)其服務(wù)半徑的修箱需求點(diǎn)提供服務(wù)時(shí)所產(chǎn)生的單箱懲罰成本，元/個(gè);C為集裝箱維修點(diǎn)服務(wù)半徑內(nèi)單個(gè)集裝箱的最高運(yùn)輸費(fèi)用，元/個(gè)。

1.1.3 決策變量

（1）Xj為0-1變量：若在j點(diǎn)設(shè)立集裝箱維修點(diǎn)，則Xj=1;否則，Xj=0。

（2）Yij為0-1變量：若集裝箱維修點(diǎn)j為修箱需求點(diǎn)i提供修箱服務(wù)，則Yij=1;否則，Yij=0。

（3）Zij為0-1變量：若集裝箱維修點(diǎn)j為超過(guò)其服務(wù)半徑的修箱需求點(diǎn)i提供修箱服務(wù)，則Zij=1;否則，Zij=0。

1.2 假設(shè)條件

（1）集裝箱維修點(diǎn)的運(yùn)營(yíng)成本僅包含建站成本、運(yùn)輸成本和懲罰成本，不包含其他成本。

（2）集裝箱維修點(diǎn)的容量不受限制，各站點(diǎn)的維修能力和維修效率相同且能完成各項(xiàng)維修任務(wù)。

（3）修箱需求點(diǎn)與集裝箱維修點(diǎn)之間的距離固定不變，運(yùn)輸線路不受天氣等因素的影響。

（4）懲罰成本與集裝箱運(yùn)量呈線性正相關(guān)關(guān)系。

（5）每個(gè)修箱需求點(diǎn)僅被唯一的集裝箱維修點(diǎn)服務(wù)，一個(gè)集裝箱維修點(diǎn)可為多個(gè)修箱需求點(diǎn)服務(wù)。

（6）集裝箱維修點(diǎn)的運(yùn)輸成本與集裝箱運(yùn)量和運(yùn)輸距離成正比。

1.3 模型構(gòu)建

本文建立以運(yùn)營(yíng)成本（包括建站成本、運(yùn)輸成本和懲罰成本）最小化為目標(biāo)的集裝箱維修點(diǎn)選址模型，目標(biāo)函數(shù)為

式中：第一項(xiàng)為建站成本和設(shè)備成本;第二項(xiàng)為運(yùn)輸成本;第三項(xiàng)為懲罰成本。由于集裝箱維修點(diǎn)的服務(wù)半徑覆蓋全部修箱需求點(diǎn)會(huì)導(dǎo)致建站成本過(guò)高，故部分修箱需求點(diǎn)可能在集裝箱維修點(diǎn)的服務(wù)半徑之外，而修箱需求點(diǎn)與集裝箱維修點(diǎn)之間的距離過(guò)長(zhǎng)會(huì)降低集裝箱周轉(zhuǎn)率和客戶滿意度;為此，本文引入懲罰成本。

模型的約束條件為

式（2）表示設(shè)立的集裝箱維修點(diǎn)的數(shù)量最少為1個(gè)，最多為N個(gè);式（3）表示只有在j點(diǎn)設(shè)立集裝箱維修點(diǎn)，才能為修箱需求點(diǎn)i提供服務(wù);式（4）表示修箱需求點(diǎn)i被唯一的集裝箱維修點(diǎn)j服務(wù);式（5）、式（6）和式（7）表示當(dāng)修箱需求點(diǎn)i與集裝箱維修點(diǎn)j之間的集裝箱運(yùn)輸費(fèi)用超過(guò)集裝箱維修點(diǎn)j服務(wù)半徑內(nèi)的最高集裝箱運(yùn)輸費(fèi)用時(shí)，Zij=1，即此時(shí)集裝箱維修點(diǎn)j需要承擔(dān)懲罰成本;式（8）和式（9）為非負(fù)約束，表示修箱需求點(diǎn)i需要維修的集裝箱數(shù)量和修箱需求點(diǎn)i與集裝箱維修點(diǎn)j之間的距離均不小于0;式（10）為取值約束，表示Xj、Yij和Zij的取值只能為0或1。

2 改進(jìn)遺傳算法設(shè)計(jì)

2.1 編碼

遺傳算法廣泛應(yīng)用于各種問(wèn)題的求解，其主要編碼方式有二進(jìn)制編碼方式、格雷碼編碼方式和浮點(diǎn)數(shù)編碼方式等，其中：二進(jìn)制編碼方式在編碼、解碼方面易于操作，是一種常用的編碼方式。二進(jìn)制編碼符號(hào)集由0和1組成，其構(gòu)成的個(gè)體基因型是一個(gè)二進(jìn)制編碼符號(hào)串，該二進(jìn)制編碼符號(hào)串的長(zhǎng)度與問(wèn)題求解精度有關(guān)。

2.2 產(chǎn)生初始種群

種群規(guī)模會(huì)影響迭代時(shí)間。一般來(lái)說(shuō)，個(gè)體數(shù)量較多的初始種群能夠同時(shí)處理更多解，從而使全局最優(yōu)解更容易被找到;但實(shí)踐表明，種群規(guī)模龐大并不能提高算法性能。本文設(shè)定初始種群的個(gè)體數(shù)量為50個(gè)。

2.3 適應(yīng)度函數(shù)

適應(yīng)度函數(shù)是度量個(gè)體適應(yīng)度值的函數(shù)，適應(yīng)度函數(shù)值可用于檢驗(yàn)個(gè)體優(yōu)劣：適應(yīng)度函數(shù)值越大，個(gè)體遺傳到下一代的概率就越高。適應(yīng)度函數(shù)值必須為非負(fù)值，對(duì)于求解最大值的目標(biāo)函數(shù)而言，如果所有函數(shù)值均為非負(fù)值，則可將目標(biāo)函數(shù)作為適應(yīng)度函數(shù)。本文構(gòu)建的目標(biāo)函數(shù)求解最小值，可通過(guò)變換目標(biāo)函數(shù)得出適應(yīng)度函數(shù)。

2.4 選擇算子

選擇算子又稱“復(fù)制算子”，指在父代群體中選擇具有強(qiáng)大生命力的個(gè)體來(lái)產(chǎn)生子代群體的過(guò)程。選擇算子方法主要有輪盤賭選擇法、隨機(jī)競(jìng)爭(zhēng)選擇法和最佳保留選擇法等，本文采用輪盤賭選擇法選擇算子。設(shè)種群規(guī)模為M，個(gè)體適應(yīng)度函數(shù)為f （xi），選擇算子的步驟如下。

（1）計(jì)算父代群體中個(gè)體的適應(yīng)度函數(shù)值f （xi）。

（2）計(jì)算父代群體中的個(gè)體被遺傳到子代群體中的概率，即

p （xi）=

（3）計(jì)算個(gè)體的累計(jì)概率，即

di=（xj），j=1，2，…，M

（4）在[0，1]區(qū)間內(nèi)產(chǎn)生一個(gè)偽隨機(jī)數(shù)r。

（5）若r

（6）重復(fù)步驟（4）和（5）共M次。

2.5 交叉算子和變異算子

交叉算子又稱“交配算子”或“重組算子”，指兩個(gè)個(gè)體間隨機(jī)交換部分基因信息，從而形成新個(gè)體的過(guò)程。交叉算子決定遺傳算法的全局搜索能力，是產(chǎn)生新個(gè)體的主要方法。變異算子作為遺傳算法的重要步驟，能夠改善遺傳算法的局部搜索能力，并維持群體的多樣性，防止出現(xiàn)早熟現(xiàn)象，是產(chǎn)生新個(gè)體的輔助方法。鑒于交叉算子和變異算子的重要性，本文通過(guò)以下公式對(duì)其進(jìn)行調(diào)整，即

Pc=

Pm=

式中：fmax為群體中最大適應(yīng)度值;favg為群體中平均適應(yīng)度函數(shù);為群體中的兩個(gè)交叉?zhèn)€體中的較大適應(yīng)度值;f為變異個(gè)體的適應(yīng)度值。

2.6 算法終止條件

選取最大迭代數(shù)作為終止條件，一般為100～ 500，本文選取200。

3 集裝箱維修點(diǎn)選址模型應(yīng)用

某公司主營(yíng)集裝箱租賃和承運(yùn)業(yè)務(wù)。隨著該公司擁有的集裝箱數(shù)量逐年增多，受損集裝箱數(shù)量也不斷增多。為了提高修箱效率和集裝箱周轉(zhuǎn)率，緩解箱源緊張等問(wèn)題，該公司計(jì)劃在華北區(qū)域現(xiàn)有的42個(gè)貨運(yùn)站點(diǎn)（見(jiàn)圖1）附近設(shè)立不多于5個(gè)集裝箱維修點(diǎn)。上述42個(gè)貨運(yùn)站點(diǎn)既是修箱需求點(diǎn)，也是備選集裝箱維修點(diǎn)。采用前述集裝箱維修點(diǎn)選址模型和改進(jìn)遺傳算法計(jì)算集裝箱維修點(diǎn)選址結(jié)果，從中選出4組較優(yōu)結(jié)果（見(jiàn)表1），并對(duì)選址結(jié)果進(jìn)行評(píng)價(jià)分析。

（1）安棚站點(diǎn) 安棚站點(diǎn)位于河南省南陽(yáng)市。南陽(yáng)市地理位置優(yōu)越，交通便捷，焦枝鐵路與寧西鐵路在城區(qū)交會(huì)，高速公路通車?yán)锍逃?00 km，是承東啟西、連南貫北的重要交通樞紐;境內(nèi)水資源豐富，河流眾多，分屬長(zhǎng)江、黃河和淮河三大水系;通信和電力條件優(yōu)越，并且地質(zhì)結(jié)構(gòu)極為穩(wěn)定，符合集裝箱維修點(diǎn)選址要求;招商引資和建設(shè)項(xiàng)目的政策環(huán)境較好，有利于集裝箱維修點(diǎn)建設(shè)和發(fā)展。安棚站點(diǎn)主要輻射洛陽(yáng)、漯河、運(yùn)城等修箱需求點(diǎn)，雖然其輻射的修箱需求點(diǎn)數(shù)量較少，但考慮到其自身修箱需求較大，在該站點(diǎn)設(shè)立集裝箱維修點(diǎn)較為合理。

（2）京唐港站點(diǎn) 京唐港站點(diǎn)位于河北省唐山市樂(lè)亭縣。樂(lè)亭縣位于環(huán)渤海地區(qū)中心地帶，被北京、天津、秦皇島和唐山四市環(huán)繞，地理位置優(yōu)越;唐港鐵路和遷曹鐵路貫穿南北，高速公路和省道覆蓋全境，境內(nèi)還有京唐港和曹妃甸港兩大港區(qū)，已形成現(xiàn)代化的交通運(yùn)輸體系;地勢(shì)低平，水電條件優(yōu)越，國(guó)電華北電力有限公司、唐山冀東氯堿有限公司等企業(yè)在此落戶。憑借優(yōu)越的自然、交通和政策等條件，樂(lè)亭縣入選2018年度全國(guó)投資潛力百?gòu)?qiáng)縣市?？紤]到京唐港站點(diǎn)的修箱需求將持續(xù)增長(zhǎng)，修箱業(yè)務(wù)發(fā)展前景廣闊，在該站點(diǎn)設(shè)立集裝箱維修點(diǎn)較為合理。

（3）農(nóng)中站點(diǎn)和湖田站點(diǎn) 農(nóng)中站點(diǎn)和湖田站點(diǎn)分別位于山東省淄博市恒臺(tái)縣和張店區(qū)。恒臺(tái)縣交通條件優(yōu)越，境內(nèi)鐵路、高速公路、國(guó)道、省道縱橫交錯(cuò);地勢(shì)南高北低，平原居多，水資源豐富。張店區(qū)則是濟(jì)南都市圈的節(jié)點(diǎn)之一。農(nóng)中站點(diǎn)與湖田站點(diǎn)直線距離11 km，交通和基礎(chǔ)設(shè)施條件相差不大。根據(jù)淄博市人民政府2018年發(fā)布的《淄博市新舊動(dòng)能轉(zhuǎn)換重大工程實(shí)施規(guī)劃》：恒臺(tái)縣將積極培育現(xiàn)代物流等服務(wù)業(yè);張店區(qū)將著力發(fā)展電子信息、會(huì)展經(jīng)濟(jì)等服務(wù)業(yè)。從發(fā)展規(guī)劃來(lái)看，農(nóng)中站點(diǎn)比湖田站點(diǎn)更適合設(shè)立集裝箱維修點(diǎn)。

（4）滄州站點(diǎn) 滄州站點(diǎn)位于河北省滄州市。滄州市東臨渤海，西鄰保定市和衡水市，北接天津市和廊坊市，南與山東省濱州市和德州市相望，地理位置優(yōu)越，是環(huán)渤海地區(qū)重要的現(xiàn)代化工業(yè)城市;境內(nèi)不僅有四通八達(dá)的高速公路和京滬、京九、朔黃等鐵路，還有集礦石碼頭、雜貨碼頭、集裝箱碼頭、原油碼頭、煤炭碼頭等于一體的綜合性大港黃驊港，已基本形成現(xiàn)代化的交通運(yùn)輸體系;地勢(shì)低平，以平原為主，地處海河流域，境內(nèi)有20多條河流，水資源豐富。從地理位置和交通條件來(lái)看，滄州站點(diǎn)適合設(shè)立集裝箱維修點(diǎn)。

（5）濮陽(yáng)東站點(diǎn) 濮陽(yáng)東站點(diǎn)位于河南省濮陽(yáng)市濮陽(yáng)縣。濮陽(yáng)縣位于河南省與山東省交界處，銜接中原經(jīng)濟(jì)區(qū)和環(huán)渤海經(jīng)濟(jì)圈;緊鄰大廣高速公路和范輝高速公路，境內(nèi)晉魯豫鐵路與京廣鐵路、京九鐵路相連，國(guó)道貫穿南北，交通四通八達(dá);地處平原，黃河、金堤河流經(jīng)境內(nèi)，水資源豐富;產(chǎn)業(yè)基礎(chǔ)較好，已形成電子電氣裝備制造、新興化工等產(chǎn)業(yè)集群，且發(fā)展前景良好。2018年，濮陽(yáng)縣入選中國(guó)工業(yè)百?gòu)?qiáng)縣（市）。從地理位置、基礎(chǔ)設(shè)施和發(fā)展前景來(lái)看，濮陽(yáng)東站點(diǎn)適合設(shè)立集裝箱維修點(diǎn)。

（6）煙臺(tái)站點(diǎn) 煙臺(tái)站點(diǎn)位于山東省煙臺(tái)市，為一等客貨運(yùn)站。煙臺(tái)市是我國(guó)近代工業(yè)的發(fā)祥地，也是環(huán)渤海地區(qū)的重要港口城市和海上絲綢之路的重要節(jié)點(diǎn);經(jīng)公路至青島和威海兩地僅需1～，境內(nèi)有藍(lán)煙鐵路等3條鐵路和煙臺(tái)港，交通網(wǎng)絡(luò)發(fā)達(dá);地形為低山丘陵，地質(zhì)結(jié)構(gòu)穩(wěn)定，境內(nèi)河流眾多，長(zhǎng)度超過(guò)5 km的河流共有百余條，水資源豐富。雖然煙臺(tái)站點(diǎn)的地理位置和交通條件較好，但從煙臺(tái)市的地形和周邊修箱需求來(lái)看，目前煙臺(tái)站點(diǎn)不適合設(shè)立集裝箱維修點(diǎn)，未來(lái)可以重點(diǎn)考慮。

（7）小結(jié) 綜上所述，雖然第四組選址結(jié)果的運(yùn)營(yíng)成本最低，但與第二組選址結(jié)果的運(yùn)營(yíng)成本相差不大，加之農(nóng)中站點(diǎn)比湖田站點(diǎn)更適合設(shè)立集裝箱維修點(diǎn)，建議將第二組的安棚、滄州、農(nóng)中、京唐港、濮陽(yáng)東等站點(diǎn)作為集裝箱維修點(diǎn)的最終選址結(jié)果（見(jiàn)圖2）。

4 結(jié)束語(yǔ)

本文首先在考慮建站成本、運(yùn)輸成本、維修需求、服務(wù)半徑等因素的基礎(chǔ)上，建立以運(yùn)營(yíng)成本最小化為目標(biāo)的集裝箱維修點(diǎn)選址模型;然后采用改進(jìn)遺傳算法對(duì)實(shí)例求解，得出4組較優(yōu)的選址結(jié)果;最后從地理位置、交通條件、基礎(chǔ)設(shè)施、政策環(huán)境等角度對(duì)選址結(jié)果進(jìn)行評(píng)價(jià)分析，并最終確定選址結(jié)果，從而較好地解決集裝箱維修點(diǎn)選址問(wèn)題。需要注意的是：本文僅研究集裝箱維修點(diǎn)選址問(wèn)題;未來(lái)有必要在考慮各種影響因素的基礎(chǔ)上，對(duì)集裝箱維修點(diǎn)調(diào)配問(wèn)題開(kāi)展研究。

參考文獻(xiàn)：

[1] 張陽(yáng). 沈陽(yáng)鐵路局鐵路通用集裝箱維修網(wǎng)點(diǎn)選址研究[D]. 北京： 北京交通大學(xué)，2016.

[2] C罭OVAS L， GARC虯 S， MAR蚇 A. Solving the uncapaci-tated multiple allocation hub location problem by means of a dual-ascent technique[J]. European Journal of Operational Research， 2007， 179（3）： 990-1007.

[3] MELO M T， NICKEL S， SALDANHA-DA-GAMA F. Faci-lity location and supply chain management―A review[J]. European Journal of Operational Research， 2009， 196（2）： 401-412.

[4] TOHYAMA H， IDA K， MATSUEDA J. A genetic algorithm for the uncapacitated facility location problem[J]. Electronics and Communications in Japan， 2011， 94（5）： 47-54.

[5] 聶鑫路，魏慶朝. 基于 PSO 的城市軌道交通應(yīng)急救援站選址研究[J]. 鐵道工程學(xué)報(bào)，2015，32（7）：100-105.

[6] RAZI N， KARATAS M. A multi-objective model for locating search and rescue boats[J]. European Journal of Operational Research， 2016， 254（1）： 279-293.

[7] 鄭茹苑. 南海水域半潛維修船最佳配備數(shù)量及選址問(wèn)題研究[D]. 大連：大連海事大學(xué)，2017.

[8] 支曉宇. 城市軌道交通應(yīng)急公交駐車點(diǎn)選址分配問(wèn)題研究[D]. 北京：北京交通大學(xué)，2018.

[9] 王俊鵬. 基于情景的鐵路應(yīng)急設(shè)施損毀不確定性選址規(guī)劃研究[D]. 南昌：華東交通大學(xué)，2018.

[10] GLAVID， MILENKOVIM， NIKOLIM， et al. Deter-mining the number and location of winter road maintenance depots―a case study of the district road network in Serbia[J]. Transportation Planning and Technology， 2018， 41（2）： 138- 153.

（編輯：張敏 收稿日期：2020-03-31）