競爭環(huán)境下的軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)設(shè)計(jì)問題

趙宇哲

(大連海事大學(xué)交通運(yùn)輸管理學(xué)院, 遼寧 大連 116026)

競爭環(huán)境下的軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)設(shè)計(jì)問題

趙宇哲

(大連海事大學(xué)交通運(yùn)輸管理學(xué)院, 遼寧 大連 116026)

針對航運(yùn)企業(yè)的重組與全球擴(kuò)張引起的競爭問題，提出了競爭環(huán)境下的軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)設(shè)計(jì)模型。模型采用基于路徑的變量作為決策變量，利用離散函數(shù)來表示航運(yùn)企業(yè)與航運(yùn)聯(lián)盟的競爭可吸引的流量(或客戶)，目的在于通過設(shè)計(jì)混合軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)，實(shí)現(xiàn)以更低的服務(wù)成本和更短的服務(wù)時間最大化可吸引的流量，建立了樞紐港口數(shù)量約束、航線連接約束、航線中轉(zhuǎn)約束、流量競爭約束等，運(yùn)用多點(diǎn)交叉遺傳算法進(jìn)行求解，最后結(jié)合亞歐航線的集裝箱海運(yùn)市場進(jìn)行實(shí)例分析，對考慮客戶需求多樣性與航運(yùn)聯(lián)盟對策下的軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)進(jìn)行設(shè)計(jì)，并驗(yàn)證了算法的計(jì)算效果。

軸-輻式網(wǎng)絡(luò)；競爭；樞紐選址；遺傳算法；集裝箱海運(yùn)

1 引言

海運(yùn)業(yè)伴隨著經(jīng)濟(jì)全球化的進(jìn)程在全球范圍內(nèi)取得了空前的發(fā)展。目前，海上運(yùn)輸承擔(dān)了國際貿(mào)易總量的89.6%、貿(mào)易總值的70.1%，已成為國際貿(mào)易中最重要的運(yùn)輸方式。貨物運(yùn)輸包括：集裝箱、普通貨物、超尺度或超重貨物，其中集裝箱的比例明顯大于其它種類貨物[1]。過去三十年，海上運(yùn)輸持續(xù)增長、年平均增長率為3.1%，運(yùn)輸需求的增長直接導(dǎo)致了海運(yùn)市場上運(yùn)力供給的快速增長。航運(yùn)企業(yè)作為海運(yùn)市場的主體，也是海運(yùn)資源重組的執(zhí)行者。為爭奪更大的市場份額，航運(yùn)企業(yè)的重組與全球擴(kuò)張，成為對海運(yùn)市場結(jié)構(gòu)和發(fā)展趨勢產(chǎn)生重要影響的新生因素。通過航運(yùn)企業(yè)的合作，建立起來的航運(yùn)聯(lián)盟可采取運(yùn)價協(xié)定、貨載分配協(xié)定等對內(nèi)措施和回扣制度、聯(lián)運(yùn)協(xié)定等對外措施消除競爭，控制海運(yùn)資源并排擠其它航運(yùn)企業(yè)[2]，這使得一些運(yùn)營規(guī)模較小或新進(jìn)入市場的航運(yùn)企業(yè)正逐漸喪失它們的市場份額。集裝箱海運(yùn)網(wǎng)絡(luò)是航運(yùn)企業(yè)提供集裝箱海運(yùn)服務(wù)的運(yùn)營基礎(chǔ)[3]，對于不屬于航運(yùn)聯(lián)盟的中小型航運(yùn)企業(yè)而言，為爭取更多的或者至少恢復(fù)已經(jīng)喪失的市場份額，如何設(shè)計(jì)能與航運(yùn)聯(lián)盟相競爭的集裝箱海運(yùn)網(wǎng)絡(luò)至關(guān)重要。

由于集裝箱海運(yùn)服務(wù)的需求派生于全球經(jīng)濟(jì)增長的需要，軸-輻式網(wǎng)絡(luò)憑借其符合全球貿(mào)易格局的結(jié)構(gòu)特性和規(guī)模經(jīng)濟(jì)效應(yīng)，能夠大幅減少船舶運(yùn)營成本、船舶運(yùn)營數(shù)量，提高船舶負(fù)載率、船舶利用率，已成為全球集裝箱海運(yùn)系統(tǒng)中最重要的網(wǎng)絡(luò)形態(tài)，引起了學(xué)術(shù)界與航運(yùn)界的廣泛關(guān)注[4-5]。軸-輻式網(wǎng)絡(luò)，根據(jù)網(wǎng)絡(luò)中是否存在直接連接可分為純軸-輻式網(wǎng)絡(luò)和混合軸-輻式網(wǎng)絡(luò)[6]。純軸-輻式網(wǎng)絡(luò)是網(wǎng)絡(luò)中所有非樞紐節(jié)點(diǎn)上的流量都要通過樞紐節(jié)點(diǎn)進(jìn)行中轉(zhuǎn)，非樞紐節(jié)點(diǎn)之間不存在直接連接的網(wǎng)絡(luò)；混合軸-輻式網(wǎng)絡(luò)中的連接除了非樞紐節(jié)點(diǎn)與樞紐節(jié)點(diǎn)之間的連接以外，還允許非樞紐節(jié)點(diǎn)之間的直接連接。由于混合軸-輻式網(wǎng)絡(luò)可根據(jù)集裝箱海運(yùn)需求的多樣性、市場競爭的加劇等因素靈活分配樞紐港口中轉(zhuǎn)并進(jìn)行設(shè)計(jì)相應(yīng)的航線連接，因此它比純軸-輻式網(wǎng)絡(luò)具有明顯的競爭優(yōu)勢，但其組織形式也更為復(fù)雜。

軸-輻式網(wǎng)絡(luò)設(shè)計(jì)問題的研究主要集中于不同情景下的樞紐選址模型設(shè)計(jì)：中心問題，中位問題，覆蓋問題等[7]。Campbell等[8]、Alumur等[9]、Campbell等[10]和Faharani等[11]對樞紐選址模型作了完整性的評述。當(dāng)前，樞紐選址模型設(shè)計(jì)的趨勢在于開發(fā)新的構(gòu)想，沿著這個研究方向已經(jīng)有很多研究進(jìn)行了多面性的探討：折扣流量問題[8-9]、樞紐容量問題[12]、網(wǎng)絡(luò)效應(yīng)問題[13]等，研究中仍有兩方面關(guān)注非常少：競爭與軸-輻式網(wǎng)絡(luò)在海運(yùn)上的應(yīng)用。 (1)關(guān)于具有競爭性的樞紐選址模型設(shè)計(jì)研究中，Marianov等[14]第一個針對出現(xiàn)在客運(yùn)和貨運(yùn)中的競爭問題進(jìn)行了樞紐選址設(shè)計(jì)，構(gòu)建的競爭原則是提供運(yùn)輸服務(wù)收入的增加；將市場上現(xiàn)存企業(yè)稱為領(lǐng)導(dǎo)者，同時將新加入市場企業(yè)稱為跟隨者；結(jié)果顯示跟隨者若想進(jìn)入市場需通過設(shè)計(jì)新的樞紐以及降價的方式促使用戶改變選擇，以最大化運(yùn)輸服務(wù)的收入。Sasaki等[15]通過設(shè)計(jì)一個Stackelberg樞紐選址模型用于制定市場中領(lǐng)導(dǎo)者和跟隨者之間的競爭策略，結(jié)果顯示若領(lǐng)導(dǎo)者的決策與跟隨者不一致，則其市場份額將會受到競爭威脅。Lüer-Villagra等[16]針對客戶選擇行為對市場競爭策略的影響，構(gòu)建一個連續(xù)型樞紐選址模型用于評估跟隨者是否該進(jìn)入市場，結(jié)果顯示跟隨者進(jìn)入市場時應(yīng)考慮客戶對價格的敏感度，否則跟隨者無法估計(jì)市場需求和可能的利潤。上述問題的樞紐選址模型設(shè)計(jì)的都是一些傳統(tǒng)的目標(biāo)函數(shù)，只考慮了成本、利潤、費(fèi)用等因素，但在具有競爭性的樞紐選址模型的應(yīng)用研究中，需結(jié)合一些其它的潛在屬性來設(shè)置更客觀的目標(biāo)函數(shù)，Eiselt等[17]拓展了Marianov等的研究，提出基于費(fèi)用/時間的吸引力函數(shù)最大化的非線性模型，用于解決客運(yùn)和貨運(yùn)中的流量競爭問題。另外，Sasaki等[18]提出競爭者之間策略應(yīng)設(shè)立非樞紐節(jié)點(diǎn)之間的直接連接以最大化其利益，且每條連接路徑上允許存在多個樞紐節(jié)點(diǎn)。 (2)關(guān)于軸-輻式網(wǎng)絡(luò)在海運(yùn)上的應(yīng)用研究中， Imai等[19]對常規(guī)船型的多港掛靠方式和大型船舶的軸-輻式網(wǎng)絡(luò)方式進(jìn)行對比分析，設(shè)計(jì)出一個可同時運(yùn)用于多港掛靠和軸-輻式網(wǎng)絡(luò)的優(yōu)化模型，并結(jié)合實(shí)例證明該模型在處理集裝箱空箱調(diào)運(yùn)問題時能獲得顯著效果；計(jì)明軍等[20]綜合考慮樞紐港口的泊位時間限制和集裝箱船舶容量限制，建立了一個以航行時間最小為目標(biāo)函數(shù)的軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)模型，通過一種嵌入式的混合遺傳算法求解出最優(yōu)船舶容量。上述研究在理論上界定了軸-輻式網(wǎng)絡(luò)在集裝箱海運(yùn)應(yīng)用中的重點(diǎn)方向，主要表現(xiàn)在模型設(shè)計(jì)中應(yīng)更多的考慮一些現(xiàn)實(shí)因素，包括關(guān)于集裝箱海運(yùn)中的運(yùn)輸費(fèi)用/時間、中轉(zhuǎn)費(fèi)用/時間、運(yùn)輸流量、基于費(fèi)用/時間的吸引力等競爭因素，關(guān)于集裝箱海運(yùn)中港口的地理位置、航線連接上的樞紐港口數(shù)量、非樞紐港口之間的直達(dá)航線等網(wǎng)絡(luò)因素。

最近，Gelareh等[21]綜合考慮了多港掛靠和軸-輻式網(wǎng)絡(luò)的運(yùn)輸方式、非樞紐港口之間直達(dá)航線的潛在優(yōu)勢、基于費(fèi)用/時間的吸引力等因素，建立了競爭環(huán)境下的班輪運(yùn)輸網(wǎng)絡(luò)設(shè)計(jì)模型，提出了海運(yùn)市場中跟隨者的競爭策略，進(jìn)一步豐富了具有競爭性的樞紐選址模型在海運(yùn)上的應(yīng)用研究，但仍存在以下局限：(1)模型采用的是基于弧的變量，盡管基于弧的變量分解具有較好的收斂性，但它比基于路徑的變量分解需要更多的時間進(jìn)行比較[22]；(2)模型未能考慮海運(yùn)市場中跟隨者對客戶需求(費(fèi)用/時間)多樣性的競爭策略，更未能關(guān)注領(lǐng)導(dǎo)者對于跟隨者進(jìn)入海運(yùn)市場的對策。

鑒于上述分析，本文針對航運(yùn)企業(yè)的重組與全球擴(kuò)張引起的競爭問題，以混合軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)為基本網(wǎng)絡(luò)結(jié)構(gòu)，引入考慮服務(wù)成本和服務(wù)時間的(流量)吸引力函數(shù)；采用基于路徑的變量，建立競爭環(huán)境下的軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)設(shè)計(jì)問題的數(shù)學(xué)模型，設(shè)計(jì)求解這一問題的遺傳算法；最后，結(jié)合亞歐航線的集裝箱海運(yùn)市場進(jìn)行實(shí)例分析，對考慮客戶需求多樣性與航運(yùn)聯(lián)盟對策下的軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)進(jìn)行設(shè)計(jì)，并驗(yàn)證了算法的計(jì)算效果。

2 模型描述

競爭環(huán)境下的軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)設(shè)計(jì)問題可定義：假設(shè)存在一個航運(yùn)聯(lián)盟A(領(lǐng)導(dǎo)者)運(yùn)營的集裝箱海運(yùn)網(wǎng)絡(luò)GA=(NA,VA)，其中NA={1,2,…,n}為港口集合，VA={(viA,vjA)│viA,vjA∈VA,j≠i}為航線集合；對應(yīng)于港口i到港口j的集裝箱OD流i→j(j≠i)上的服務(wù)成本CijA和服務(wù)時間TijA。對于一個同為n個港口服務(wù)，但不屬于航運(yùn)聯(lián)盟A的中小型航運(yùn)企業(yè)B(跟隨者)，在調(diào)查客戶需求及其位置的基礎(chǔ)上，需要選擇位于同一集裝箱海運(yùn)網(wǎng)絡(luò)范圍內(nèi)的p(p≤n)個港口作為樞紐港口，通過設(shè)計(jì)混合軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)GB=(NB,VB)，實(shí)現(xiàn)以更低的服務(wù)成本Cij和更短的服務(wù)時間Tij為客戶提供服務(wù)，爭取集裝箱OD流i→j上的流量(或客戶)最大化其市場份額。

3 模型分析

3.1 吸引力函數(shù)

對于一個已經(jīng)存在的航運(yùn)聯(lián)盟A運(yùn)營的集裝箱海運(yùn)網(wǎng)絡(luò)GA(NA,VA)，若航運(yùn)企業(yè)B由于其設(shè)計(jì)的軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)GB的結(jié)構(gòu)特性，對應(yīng)于港口i到港口j的集裝箱OD流i→j上的服務(wù)水平優(yōu)于A，不失一般性，假設(shè)Cij≤CijA且Tij≤TijA，則航運(yùn)聯(lián)盟A在集裝箱OD流i→j上的流量將被吸引至航運(yùn)企業(yè)B上[14]。但是，是否只要滿足Cij≤CijA，且Tij≤TijA時，航運(yùn)聯(lián)盟A在集裝箱OD流i→j上的流量就將全部吸引至航運(yùn)企業(yè)B上? 現(xiàn)實(shí)情況中，服務(wù)水平的提升并不會導(dǎo)致流量全部轉(zhuǎn)變，因其還會受到地理位置、環(huán)境等因素影響。但一般來說，服務(wù)水平提升的程度越高，可獲取的流量要顯著大于那些服務(wù)水平只有輕微提升的轉(zhuǎn)變。因此，本文采用離散函數(shù)Oijm,m∈{1,2,…,M}來表示航運(yùn)企業(yè)B對于航運(yùn)聯(lián)盟A在集裝箱OD流i→j上流量的吸引力，見圖1。對應(yīng)于集裝箱OD流i→j，考慮服務(wù)成本的吸引力函數(shù)為Oijmc考慮服務(wù)時間的吸引力函數(shù)為Oijmt，于是有：

(1)

其中Oij1c,Oij2c,…,OijMc∈(0,1)表示航運(yùn)企業(yè)B利用更低的服務(wù)成本在集裝箱OD流i→j上爭取的流量比例，0<θ1c<θ2c<…<θMc≤1，且有Oij(M+1)c=0；

(2)

其中Oij1t,Oij2t,…,OijMt∈(0,1)表示航運(yùn)企業(yè)B利用更短的服務(wù)時間在集裝箱OD流i→j上爭取的流量比例，0<δ1t<δ2t<…<δMt≤1，且有Oij(M+1)t=0。

采用離散函數(shù)來表示航運(yùn)企業(yè)的(流量)吸引力函數(shù)，需注意以下兩方面[21]：(1)與空運(yùn)、鐵運(yùn)等公共運(yùn)輸中客戶對服務(wù)成本/時間等高度敏感不一致，集裝箱海運(yùn)中的客戶在更換航運(yùn)企業(yè)為其服務(wù)時，除了服務(wù)成本、服務(wù)時間外還考慮地理位置、環(huán)境等因素，因此航運(yùn)企業(yè)B作為客戶選擇的基礎(chǔ)是具有有限的臨界值，即θmc與δmt的設(shè)定應(yīng)在一定的范圍內(nèi)；(2)對于提供相同服務(wù)水平的兩家企業(yè)，與空運(yùn)、鐵運(yùn)等公共運(yùn)輸中客戶易于選擇品牌更好/規(guī)模更大的服務(wù)企業(yè)不一致，集裝箱海運(yùn)中的客戶則更傾向于同時選擇兩家航運(yùn)企業(yè)為其服務(wù)，這可促進(jìn)航運(yùn)企業(yè)競爭以制止價格壟斷行為，因此航運(yùn)企業(yè)B在服務(wù)成本上更易爭取到不錯的市場份額，即Oijmc與θmc之間，Oijmt與δmt之間關(guān)系的設(shè)定存在一定的差異，見圖1(a)與圖1 (b)。

3.2 混合軸-輻式網(wǎng)絡(luò)

對于一個不屬于航運(yùn)聯(lián)盟A的中小型航運(yùn)企業(yè)B，進(jìn)入已知存在的集裝箱海運(yùn)市場時，所有港口(n個)的位置是已知的，但并不確定它們是樞紐港口還是非樞紐港口，軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)GB設(shè)計(jì)的關(guān)鍵在于確定哪些港口作為樞紐港口 (剩余港口自然作為非樞紐港口)以及樞紐港口與非樞紐港口之間、樞紐港口之間和非樞紐港口之間的航線連接。根據(jù)上述分析，競爭環(huán)境下的軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)GB屬于混合軸-輻式網(wǎng)絡(luò)，非樞紐港口不僅可分配給多個樞紐港口， 還可直接分配給多個非樞紐港口，但中轉(zhuǎn)只能在樞紐港口進(jìn)行，見圖2。這種分配關(guān)系與對應(yīng)于港口i到港口j的集裝箱OD流i→j上的流量、運(yùn)輸費(fèi)用/時間，中轉(zhuǎn)費(fèi)用/時間等因素相關(guān)。

圖1 航運(yùn)企業(yè)B的離散吸引力函數(shù)

圖2 軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)GB結(jié)構(gòu)

軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)GB=(NB,VB)中，港口集合NB=HB∪SB={1,2,…,n}，其中HB表示樞紐港口集合，SB表示非樞紐港口集合；航線集合VB={(viB,vjB)│viB,vjB∈VB,j≠i}，其中viB表示離開起點(diǎn)港口i∈NB的集裝箱OD流向集合，vjB表示到達(dá)訖點(diǎn)港口j∈NB的集裝箱OD流向集合。對應(yīng)于起點(diǎn)港口i到訖點(diǎn)港口j的集裝箱OD流i→j上不允許超過2個樞紐港口，集裝箱OD流i→j上的航線連接只可能是2,3或4個港口，于是有[23]：2個港口的航線連接：(i,j) ∈HB×HB∪HB×SB∪SB×HB∪SB×SB；3個港口的航線連接：(i,k,j) ∈SB×HB×HB∪SB×HB×SB∪HB×HB×SB；4個港口的航線連接：(i,k,l,j) ∈SB×HB×HB×SB。

上述的軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)GB=(NB,VB)中，對應(yīng)于起點(diǎn)港口i到訖點(diǎn)港口j的集裝箱OD流i→j上的航線連接不具有唯一性。因此，需要在幾條航線連接中選擇其一作為集裝箱OD流i→j。

4 模型構(gòu)建

4.1 基本假設(shè)

假設(shè)1 軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)結(jié)構(gòu)相對穩(wěn)定，不考慮災(zāi)害、戰(zhàn)爭等突發(fā)狀況對網(wǎng)絡(luò)結(jié)構(gòu)的影響。

假設(shè)2 軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)中所有集裝箱OD流i→j上的流量已知且固定，不考慮集裝箱船舶運(yùn)力的限制。

假設(shè)3 軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)中所有集裝箱OD流i→j上的服務(wù)成本Cij和服務(wù)時間Tij為已知；考慮到樞紐港口之間運(yùn)輸?shù)囊?guī)模經(jīng)濟(jì)效應(yīng)，折扣系數(shù)α為已知。

假設(shè)4 軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)中各個港口的中轉(zhuǎn)費(fèi)用Cktr和中轉(zhuǎn)時間Tktr為已知。

假設(shè)5 軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)中各個港口的容量滿足所有集裝箱流i→j上的流量之和，不考慮港口容量限制約束。

4.2 模型參數(shù)

i為起點(diǎn)港口；j為訖點(diǎn)港口；k,l為樞紐港口；p為樞紐港口的數(shù)量；n為所有港口的數(shù)量。

Cik、Ckl、Clj分別為起點(diǎn)港口i與樞紐港口k、樞紐港口k與樞紐港口l、樞紐港口l與訖點(diǎn)港口j之間的運(yùn)輸費(fèi)用；Cij為起點(diǎn)港口i與訖點(diǎn)港口j(i,j為非樞紐港口)之間的直達(dá)運(yùn)輸費(fèi)用。

α為樞紐港口k,l之間運(yùn)輸?shù)恼劭巯禂?shù)，0<α<1(規(guī)模經(jīng)濟(jì)效應(yīng))。

Tik、Tkl、Tlj分別為起點(diǎn)港口i與樞紐港口k、樞紐港口k與樞紐港口l、樞紐港口l與訖點(diǎn)港口j之間的運(yùn)輸時間；Tij為起點(diǎn)港口i與訖點(diǎn)港口j(i,j為非樞紐港口)之間的直達(dá)運(yùn)輸時間。

Cktr為樞紐港口k的中轉(zhuǎn)費(fèi)用；Tktr為樞紐港口k的中轉(zhuǎn)時間。

4.3 決策變量

hk為樞紐港口選擇變量。若港口k被選擇為樞紐港口，則hk=1；否則hk=0。

yij為起訖港口i,j之間的集裝箱OD流i→j上不存在樞紐港口k(k,l)的選擇變量。若存在2個港口的航線連接(i,j)，其中k≠i,j,l≠i,j，則yij=1；否則yij=0。

4.4 數(shù)學(xué)模型

競爭環(huán)境下的軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)設(shè)計(jì)問題相應(yīng)的數(shù)學(xué)模型為：

(3)

(4)

(5)

(6)

(7)

yij≤1-hi, ?i,j≠i

(8)

yij≤1-hj, ?i,j≠i

(9)

(10)

(11)

(12)

(13)

(14)

(15)

(16)

5 模型算法

(2)當(dāng)m=1時,G ～=S3或G為2n(2n-1)階Frobenius群,其中Sylow 2-子群正規(guī),2n-1為素?cái)?shù).

由于，模型(3)-(16)是NP-難問題，應(yīng)用標(biāo)準(zhǔn)遺傳算法(GA)求解模型(3)-(16)時，基于樞紐港口的航線連接表示的編碼方式是最自然、最常用的個體表示方法，但為保證航運(yùn)企業(yè)B決定設(shè)計(jì)樞紐港口的數(shù)量限制，需要在交叉算子和變異算子上對標(biāo)準(zhǔn)遺傳算法作出適當(dāng)調(diào)整。

5.1 染色體編碼

軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)GB的關(guān)鍵在于樞紐港口集合HB?NB的確定，因此采用實(shí)數(shù)編碼的方式對樞紐港口選擇變量hk進(jìn)行染色體編碼。每一元素的值表示港口是否選擇作為樞紐港口，1表示樞紐港口，0表示非樞紐港口。

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

由于采用輪盤賭算法進(jìn)行選擇，符合遺傳算法的適應(yīng)度函數(shù)值為最大時效果最好[24]。因此適應(yīng)度函數(shù)取目標(biāo)函數(shù)，則有：

FitV=Q

(17)

為保證子代包含的樞紐港口數(shù)量與父代相同，采用多點(diǎn)交叉的方式[25]。對于一組父代個體，以交叉概率Pc選擇需要交叉的個體，從右至左尋找第一次出現(xiàn)父代個體1編碼為“1”而父代個體2編碼為“0”的編碼位置i，同時從左至右尋找第1次出現(xiàn)父代個體1編碼為“0”而父代個體2編碼為“1”的編碼位置j，相互調(diào)換位置j與i處兩個父代個體的編碼；重復(fù)此操作直至j≥i。此時，子代的樞紐港口數(shù)量與父代個體相同。

5.4 變異

為保證變異后的子代包含的樞紐港口數(shù)量不變，采用雙點(diǎn)變異的方式。對于一組子代個體，以變異概率Pm選擇需要變異的個體，尋找出可能的組合，互換其相應(yīng)的位置。將變異生成的個體替代原有的子代個體，可保證遺傳算法的有效性，同時使遺傳算法保持群體的多樣性。

6 實(shí)例分析

6.1 數(shù)據(jù)選取

為驗(yàn)證競爭環(huán)境下的軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)設(shè)計(jì)問題，以亞歐航線的集裝箱海運(yùn)市場為例，數(shù)據(jù)主要來源于CI-Online國際集裝箱在線、環(huán)球運(yùn)費(fèi)網(wǎng)、海歌Seagle以及基于BLM-Shipping軟件的測算。對于某個已知存在的航運(yùn)聯(lián)盟A，由于數(shù)據(jù)的獲取難度較大，選取其運(yùn)營的集裝箱海運(yùn)網(wǎng)絡(luò)GA覆蓋的港口數(shù)量NA={1,2,…,10}，對應(yīng)于起訖港口之間的服務(wù)成本CijA和服務(wù)時間TijA，見表1。對于一個不屬于航運(yùn)聯(lián)盟A的中小型航運(yùn)企業(yè)B，對應(yīng)于起訖港口之間的服務(wù)成本Cij和服務(wù)時間Tij，見表2。

6.2 計(jì)算結(jié)果與算法比較

選取p=3，μ=0.7；設(shè)置模型算法基本參數(shù)：終止進(jìn)化代數(shù)=300，種群規(guī)模=20，交叉概率Pc=0.7，變異概率Pm=0.07，利用Matlab2013a編程運(yùn)行，計(jì)算結(jié)果見圖3。

航運(yùn)企業(yè)B設(shè)計(jì)的軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)GB上，其中香港、新加坡和安特衛(wèi)普選擇作為樞紐港口；此時，航運(yùn)企業(yè)B同航運(yùn)聯(lián)盟A競爭可獲取的總流量為56.68。不難發(fā)現(xiàn)，亞歐航線上大部分港口之間的航線連接都經(jīng)過香港和安特衛(wèi)普中轉(zhuǎn)， 以借助樞紐港口之間干線運(yùn)輸?shù)囊?guī)模經(jīng)濟(jì)效應(yīng)降低服務(wù)成本；另一方面，受限于地理位置的因素，亞洲與歐洲兩個區(qū)域內(nèi)的港口之間更多采用的是直達(dá)運(yùn)輸，包括支線(樞紐港口與非樞紐港口之間)與支線直達(dá)航線(非樞紐港口之間)，以避免繞道產(chǎn)生過多的服務(wù)成本和服務(wù)時間。

表1 航運(yùn)聯(lián)盟A的服務(wù)成本CijA與服務(wù)時間TijA

表2 航運(yùn)企業(yè)B的服務(wù)成本Cij與服務(wù)時間Tij

圖3 權(quán)重μ=0.7時軸-輻式集裝箱海運(yùn)網(wǎng)路GB的設(shè)計(jì)結(jié)果

為驗(yàn)證多點(diǎn)交叉遺傳算法的有效性，分別選取p=2,3,4,5,6，將其與允許不可行解的遺傳算法進(jìn)行對比，見表3。允許不可行解的遺傳算法在進(jìn)行交叉操作時，不考慮子代的樞紐港口數(shù)量限制，直接對所有子代進(jìn)行適應(yīng)度函數(shù)計(jì)算。由表3可知，利用允許不可行解的遺傳算法求出的(近似)最優(yōu)解中，樞紐港口數(shù)量存在波動，且均與預(yù)期確定的樞紐港口數(shù)量p不符，為非可行解。而多點(diǎn)交叉遺傳算法的(近似)最優(yōu)解中，樞紐港口數(shù)量分別為2,3,4,5,6，均為可行解。顯然，相比于允許不可行解的遺傳算法，多點(diǎn)交叉遺傳算法可避免非可行解的產(chǎn)生，對競爭環(huán)境下的軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)設(shè)計(jì)問題具有較好的求解效果。

表3 允許不可行解的遺傳算法與多點(diǎn)交叉遺傳算法的對比分析

6.3 不同權(quán)重系數(shù)下的設(shè)計(jì)分析

航運(yùn)企業(yè)B設(shè)計(jì)其運(yùn)營的軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)同時受到服務(wù)成本和服務(wù)時間兩個因素的影響，為爭取到更多的流量，航運(yùn)企業(yè)B可根據(jù)集裝箱海運(yùn)市場中客戶需求來設(shè)定不同的重視程度(μ, 1-μ)。選取p=3；μ=0, 0.1, 0.2, …, 1。設(shè)置上述模型算法基本參數(shù)，計(jì)算結(jié)果見表4。

由表4可知，隨著μ的增加，航運(yùn)企業(yè)B逐漸重視其提供的服務(wù)成本，樞紐港口的選擇發(fā)生了一定的變化，由(2, 4, 5)經(jīng)(5, 6, 7)變?yōu)?5, 6, 10)，同時樞紐港口的中轉(zhuǎn)次數(shù)由40增至81?？梢姡瑯屑~港口在軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)中具有重要的地位。當(dāng)μ=0時，航運(yùn)企業(yè)B僅以減少服務(wù)時間入手，試圖憑借較少的服務(wù)時間吸引集裝箱海運(yùn)市場中的客戶；此時，樞紐港口的中轉(zhuǎn)次數(shù)只有40，支線直達(dá)航線數(shù)量最大，為24。這說明，對于一些對服務(wù)時間要求較高的客戶，航運(yùn)企業(yè)可通過降低在樞紐港口的中轉(zhuǎn)次數(shù)以減少服務(wù)時間。當(dāng)μ=0.1時，樞紐港口的中轉(zhuǎn)次數(shù)增至64，增長幅度達(dá)到60%，支線直達(dá)航線數(shù)量由24減至18。這表明，航運(yùn)企業(yè)B在設(shè)計(jì)軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)時開始考慮服務(wù)成本，通過提高在樞紐港口的中轉(zhuǎn)次數(shù)，實(shí)現(xiàn)樞紐港口之間的干線運(yùn)輸帶來的規(guī)模經(jīng)濟(jì)效應(yīng)，減少服務(wù)成本以爭取獲得更多的流量。

表4 不同權(quán)重系數(shù)的軸-輻式集裝箱海運(yùn)

除此之外，在μ由0.1增至0.9的過程中，樞紐港口的選擇一直未發(fā)生變化，中轉(zhuǎn)次數(shù)的變化也不明顯，支線直達(dá)航線數(shù)量也始終未變，但航運(yùn)企業(yè)B同航運(yùn)聯(lián)盟A競爭可獲取的總流量由34.38增至63.16；這在一定程度上說明了服務(wù)成本在集裝箱海運(yùn)市場中的重要性?；蛘?，由于相同航線在服務(wù)時間上差異不大，航運(yùn)企業(yè)B在剛剛進(jìn)入市場時為更好地同航運(yùn)聯(lián)盟A展開競爭，應(yīng)考慮將其提供的服務(wù)成本盡量壓低；另一方面，受限于地理位置的因素，部分起訖港口之間的支線直達(dá)航線比較固定，可見，支線直達(dá)航線數(shù)量受μ的影響較小。

6.4 不同降價系數(shù)下的設(shè)計(jì)分析

航運(yùn)聯(lián)盟A若發(fā)現(xiàn)集裝箱海運(yùn)市場中有競爭者出現(xiàn)，且感知到來自航運(yùn)企業(yè)B的威脅?？紤]到服務(wù)時間受限于航海技術(shù)，短時間內(nèi)很難通過減少服務(wù)時間奪回流量，航運(yùn)聯(lián)盟A將通過降價策略 (除以降價系數(shù)ψ)來制止航運(yùn)企業(yè)B的競爭。選取ψ=1.3,1.4,1.5,1.6,1.7,1.8, 1.9, 2.0；μ=0.7時；設(shè)置上述模型算法基本參數(shù)，計(jì)算結(jié)果見表5。

表5 不同降價倍數(shù)下的軸-輻式集裝箱海運(yùn)

由表5可知，隨著ψ的增加，航運(yùn)企業(yè)B競爭所獲取的總流量由39.2降至21.88，同時樞紐港口的中轉(zhuǎn)次數(shù)由74降到46；這說明，原本依靠規(guī)模經(jīng)濟(jì)效應(yīng)產(chǎn)生的成本優(yōu)勢已經(jīng)喪失，航運(yùn)企業(yè)B不得不保留某些具有時間優(yōu)勢的航線，支線直達(dá)航線數(shù)量由18增至24。相比于表3中μ=0的情況，當(dāng)ψ=1.7后，樞紐港口的中轉(zhuǎn)次數(shù)穩(wěn)定于46，稍大于僅依靠時間優(yōu)勢時樞紐港口的中轉(zhuǎn)次數(shù)40，這主要是由航運(yùn)聯(lián)盟A運(yùn)營的集裝箱海運(yùn)網(wǎng)絡(luò)GA中某些起訖港口之間航線的服務(wù)成本過高造成的。即使航運(yùn)聯(lián)盟A的降價系數(shù)達(dá)到2.0，航運(yùn)企業(yè)B仍可在這些航線上獲取流量。這在一定程度上表明，對于某些在航運(yùn)聯(lián)盟A運(yùn)營下但不受重視的航線，航運(yùn)企業(yè)B可通過合理的樞紐港口選擇策略來減少服務(wù)成本。以獲取該航線上的流量。

值得注意的是，由表4-5還可發(fā)現(xiàn)，在不同權(quán)重系數(shù)μ與不同降價系數(shù)ψ下的軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)GB設(shè)計(jì)中，港口5香港一直被航運(yùn)企業(yè)B選擇作為樞紐港口。這說明，具有地理優(yōu)勢(到其它港口的服務(wù)成本更低或服務(wù)時間更短)的港口在軸-輻式集裝箱海運(yùn)網(wǎng)絡(luò)設(shè)計(jì)中更容易選擇作為樞紐港口。

7 結(jié)語

[1]NasrinA,RezaZF,MarkG.Networkdesignapproachforhubports-shippingcompaniescompetitionandcooperation[J].TransportationResearchPartA, 2013, 48(2): 1-18.

[2] 王成金. 世界航運(yùn)企業(yè)重組及其對航運(yùn)網(wǎng)絡(luò)結(jié)構(gòu)的影響——兼論對中國港口體系的影響[J]. 世界地理研究, 2008, 17(1): 94-104.

[3]DucruetC,NotteboomT.Theworldwidemaritimenetworkofcontainershipping:spatialstructureandregionaldynamics[J].Globalnetworks, 2012, 12(3): 395-423.

[4]WangChengjin,WangJiaoe.Spatialpatternoftheglobalshippingnetworkanditshub-and-spokesystem[J].ResearchinTransportationEconomics, 2011, 32(1): 54-63.

[5]MengQiang,WangXinchang.Intermodalhub-and-spokenetworkdesign:Incorporatingmultiplestakeholdersandmulti-typecontainers[J].TransportationResearchPartB, 2011, 45(4): 724-742.

[6] 李陽. 軸-輻式網(wǎng)絡(luò)理論及應(yīng)用研究[D]. 上海: 復(fù)旦大學(xué), 2006.

[7] 傅少川, 胡夢飛, 唐方成. 禁忌搜索算法在單分配多樞紐軸輻式物流網(wǎng)絡(luò)中的應(yīng)用[J]. 中國管理科學(xué), 2012, 20(3): 145-151.

[8]CampbellJF.Hublocationfortimedefinitetransportation[J].ComputersandOperationsResearch, 2009, 36(12): 3107-3116.

[9]AlumurS,KaraBY.Networkhublocationproblems:thestateoftheart[J].EuropeanJournalofOperationalResearch, 2008, 190(1): 1-21.

[10]CampbellJF,O'KellyME.Twenty-fiveyearsofhublocationresearch[J].TransportationScience, 2012, 46(2): 153-169.

[11]FaharaniR,HetmakfarM,BolooriA,etal.Hublocationproblems:areviewofmodels,classification,solutiontechniques,andapplications[J].ComputersandIndustrialEngineering, 2013, 64(4): 1096-1109.

[12]GabrielaM,BerndW.Hublocator:Anexactsolutionmethodforthemultipleallocationhublocationproblem[J].ComputersandOperationsResearch, 2002, 29(6): 715-739.

[13]EconomidesN.Theeconomicsofnetworks[J].InternationalJournalofIndustrialOrganization, 1996, 1(14): 673-700.

[14]MarianovV,SerraD,ReVelleC.Locationofhubsinacompetitiveenvironment[J].EuropeanJournalofOperationalResearch1999, 114(2): 363-371.

[15]SasakiM,SuzukiA,DreznerZ.Ontheselectionofhubairportsfortheairlinehub-and-spokesystem[J].ComputersandOperationsResearch, 1999, 26(14):1411-1422.

[16]Lüer-VillagraA,MarianovV.Acompetitivehublocationandpricingproblem[J].EuropeanJournalofOperationalResearch, 2013, 231(3): 734-744.

[17]EiseltHA,MarianovV.Aconditionalp-hublocationproblemwithattractionfunctions[J].ComputersandOperationsResearch, 2009, 36(12): 3128-3135.

[18]SasakiM,CampbellJF,ErnstAT,etal.Hubarclocationwithcompetition[R].WorkingPaper,TechnicalReportoftheNanzanAcademicSocietyInformationSciencesandEngineering, 2009.

[19]ImaiA,ShintaniK,PapadimitriouS.Multi-portvs.hub-and-spokeportcallsbycontainerships[J].TransportationResearchPartE, 2009, 45(5): 740-757.

[20] 計(jì)明軍, 陳哲, 王清斌. 集裝箱船舶支線運(yùn)輸航線優(yōu)化算法[J]. 交通運(yùn)輸工程學(xué)報(bào), 2011, 11(4): 68-75.

[21]GelarehS,NickelS,PisingerD.Linershippinghubnetworkdesigninacompetitiveenvironment[J].TransportationResearchPartE:LogisticsandTransportationReview, 2010, 46(6): 991-1004.

[22]BektasT,ChoumanM,CrainicTG.Lagrangean-baseddecompositionalgorithmsformulticommoditynetworkdesignproblemswithpenalizedconstraints[J].Networks, 2010, 55(3): 171-180.

[23] 倪玲霖, 史峰. 多分配快遞軸輻網(wǎng)絡(luò)的樞紐選址與分配優(yōu)化方法[J]. 系統(tǒng)工程理論與實(shí)踐, 2012, 32(2): 441-448.

[24] 陸婧, 楊忠振, 劉瑞菊. 考慮發(fā)車間隔與乘客人數(shù)關(guān)系的機(jī)場長途巴士時刻表優(yōu)化設(shè)計(jì)[J]. 系統(tǒng)工程理論與實(shí)踐, 2013, 33(8): 2097-2104.

[25] 鄧亞娟, 陳小鴻, 楊超. 需求不確定的樞紐輻射式航線網(wǎng)絡(luò)設(shè)計(jì)[J]. 交通運(yùn)輸工程學(xué)報(bào), 2009, 9(6): 68-74.

Hub-and-spokeContainerShippingNetworkDesigninaCompetitiveEnvironment

ZHAO Yu-zhe

(College of Transportation Management, Dalian Maritime University, Dalian 116026, China)

A novel mixed integer linear programming mathematical model for the hub-and-spoke container shipping network design in a competitive environment is proposed for the competition problem caused by reorganization and global expansion of container shipping liners. Path-based formulations are adopt as decision variables in the model. A discretized attraction function is introduced to indicate the number of network traffic (customers) of the container shipping alliance that can be attracted to the newcomer container shipping liners through competition. The aim of the model is to maximize the number of network traffic attracted which depends on lower transportation cost and shorter service time - by designing its mixed hub-and-spoke container shipping network. There are 13 constraints including the number of hub ports, routes connection, route transition and network traffic competition to ensure competition factors and network design factors the model should consider. As it is a NP-hard problem, a multi-point crossover genetic algorithm is given to solve the problem. By comparing the genetic algorithm allowed infeasible solutions, the results show that the multi-point crossover genetic algorithm which can avoid infeasible solutions has a better solution. Finally, the example of a container shipping market of Asia-Europe routes is analyzed, and the major data sources are obtained from CI (Containerization International)-Online, 100allin, Seagle and BLM-Shipping. Hub-and-spoke container shipping networks are designed for the cases of the various demands of customers and the countermeasure of container shipping alliance. The calculation results show that (1) Small and medium enterprises should consider that the costs reduction offered by the scale economy for transferring origin-destination (O-D) flows among the reasonable set of hub ports. (2) Since direct routes between the origin-destination (O-D) ports are relatively fixed and can avoid excessive costs and time generated by the detour, small and medium enterprises should consider that cutting costs when entering the container shipping market. (3) Small and medium enterprises should focus on the customer demands from the container shipping market not covered by container shipping alliance, in order to design the hub-and-spoke container shipping network by reasonable hub port selection strategy.

hub-and-spoke network; competition; hub location; genetic algorithm; container shipping

2013-11-21;

2014-2-27

長江學(xué)者和創(chuàng)新團(tuán)隊(duì)發(fā)展計(jì)劃資助項(xiàng)目(IRT13048); 國家自然科學(xué)基金資助項(xiàng)目(71273037,71403035); 遼寧省高校創(chuàng)新團(tuán)隊(duì)支持計(jì)劃資助項(xiàng)目(LT2013011); 交通運(yùn)輸部交通軟科學(xué)項(xiàng)目(2013-322-225-240); 中央高校基本科研業(yè)務(wù)費(fèi)專項(xiàng)資金資助項(xiàng)目(3132013080, 3132015218)

趙宇哲(1983-)，男(漢族)，黑龍江大慶人，大連海事大學(xué)交通運(yùn)輸管理學(xué)院，講師，博士，研究方向：港口規(guī)劃與管理、發(fā)展決策與管理.

1003-207(2015)07-0103-10

10.16381/j.cnki.issn1003-207x.2015.07.013

U6-9; O221.4

A