甄遠(yuǎn)迪 楊 斌

(上海海事大學(xué)物流研究中心 上海 201306)

0 引 言

當(dāng)前多式聯(lián)運(yùn)備受重視，它是物流運(yùn)輸?shù)母呒壈l(fā)展階段，然而多式聯(lián)運(yùn)網(wǎng)絡(luò)的復(fù)雜性需要及時(shí)響應(yīng)上游操作不斷變化的動態(tài)和不確定，運(yùn)輸還需要同時(shí)兼顧多個(gè)目標(biāo)，比如成本、時(shí)間、容量約束以及低碳等。但是，運(yùn)輸方式的不同其技術(shù)經(jīng)濟(jì)特點(diǎn)也不一樣，以致這些目標(biāo)之間具有一定的效益悖反規(guī)律。因此，多式聯(lián)運(yùn)問題比單一模式的運(yùn)輸更具有挑戰(zhàn)性，一方面是因?yàn)樵诙嗍铰?lián)運(yùn)中，最優(yōu)路徑不一定是最短路，因?yàn)檫€要考慮一些額外的成本，如時(shí)間、金錢；另一方面，還需要考慮一系列新的約束，這些約束與在節(jié)點(diǎn)轉(zhuǎn)換交通方式有關(guān)，比如將貨物從卡車轉(zhuǎn)移到火車或船上的操作所帶來的一系列問題。

對此，國內(nèi)外學(xué)者近幾十年來也對多式聯(lián)運(yùn)網(wǎng)絡(luò)規(guī)劃做出大量研究，模型構(gòu)建大多數(shù)還是基于最短路或者總成本最低為目標(biāo)。佟璐等[1]用蟻群算法求解以運(yùn)輸時(shí)間和運(yùn)輸路徑最短為優(yōu)化目標(biāo)的多式聯(lián)運(yùn)運(yùn)輸路徑抉擇問題。李雙琳等[2]求解配送總時(shí)間最短和配送物資未滿足的總損失最小的多目標(biāo)選址—多式聯(lián)運(yùn)問題采用的是多目標(biāo)遺傳算法，在其算法中使用了二維編碼的非支配排序。楊秋秋等[3]和何艷梅等[4]則是構(gòu)建了以成本最低和時(shí)間最短為優(yōu)化目標(biāo)的雙目標(biāo)多式聯(lián)運(yùn)最短路徑選擇模型。黃麗霞等[5]以多式聯(lián)運(yùn)過程中總成本最小和總風(fēng)險(xiǎn)最小為目標(biāo)，構(gòu)建了雙目標(biāo)0-1線性規(guī)劃模型。以上都是基于帕累托分析的雙目標(biāo)優(yōu)化模型。一些文章考慮的因素較多，Bhattacharya 等[6]用混合整數(shù)規(guī)劃來優(yōu)化考慮多個(gè)成本和額外的容量約束的多式聯(lián)運(yùn)網(wǎng)絡(luò)時(shí)間表。Modesti等[7]利用一個(gè)特設(shè)的效用函數(shù)來權(quán)衡成本和時(shí)間，選擇最優(yōu)路徑使總成本、時(shí)間和客戶不便最小化。Sun等[8-9]用帕累托最優(yōu)來選擇多式聯(lián)運(yùn)最優(yōu)路線，以達(dá)到總運(yùn)輸成本和時(shí)間最優(yōu)，然后考慮了商品流路徑、鐵路服務(wù)的特定時(shí)刻表、二氧化碳排放量等的總成本最優(yōu)，雖然考慮因素多，但是轉(zhuǎn)化成單一目標(biāo)函數(shù)來求解。Kengpol等[10]為降低成本、提前期、風(fēng)險(xiǎn)和二氧化碳排放，通過層次分析法和數(shù)據(jù)包絡(luò)分析來確定目標(biāo)函數(shù)，最終通過0-1規(guī)劃計(jì)算最優(yōu)路徑。賀竹磬等[11]考慮了時(shí)間和容量約束，優(yōu)先考慮了時(shí)效性物流整體費(fèi)用最小化，以0-1整數(shù)規(guī)劃求解。陳怡霏[12]將碳排放量和時(shí)間都轉(zhuǎn)化為相應(yīng)成本，構(gòu)建集裝箱多式聯(lián)運(yùn)過程總成本的單目標(biāo)模型。李高波[13]將容量約束和服務(wù)時(shí)間作為約束，構(gòu)建了以成本最小化為目標(biāo)函數(shù)的單目標(biāo)模型。付曉鳳等[14]建立了基于成本和時(shí)間一體化的單目標(biāo)多式聯(lián)運(yùn)路徑選擇模型。王正彬等[15]將經(jīng)營收入、轉(zhuǎn)運(yùn)成本、因延誤的保存成本和損失費(fèi)用、時(shí)間價(jià)值這些因素都轉(zhuǎn)換為成本費(fèi)用構(gòu)建了單目標(biāo)規(guī)劃問題。這些多目標(biāo)規(guī)劃研究，雖然考慮的因素較多，但最終研究是將所考慮的因素化作單目標(biāo)函數(shù)來處理。因?yàn)閯討B(tài)不確定性在運(yùn)輸過程中的存在，所以有學(xué)者認(rèn)為客戶需求量和運(yùn)輸時(shí)間可以當(dāng)作運(yùn)輸過程中的動態(tài)不確定性。Wang等[16]將目標(biāo)設(shè)為最小成本，提出了新的帶有時(shí)間約束條件的混合整數(shù)模型，運(yùn)用圖像變換，能夠把問題變成節(jié)點(diǎn)作業(yè)帶有運(yùn)輸時(shí)間約束和隨機(jī)特性的最短路模型。王慧等[17]優(yōu)先考慮運(yùn)輸需求量是模糊的情況下，建立了以總成本最優(yōu)為主要目的的多式聯(lián)運(yùn)運(yùn)輸方式和箱型運(yùn)輸方式的優(yōu)化模型，并以改進(jìn)型粒子蟻群算法求解。張得志等[18]則有限考慮了多式聯(lián)運(yùn)中轉(zhuǎn)運(yùn)時(shí)間和運(yùn)輸?shù)牟淮_定性、時(shí)間窗限制以及轉(zhuǎn)運(yùn)的可能性等情況，建立了具有時(shí)效性的多式聯(lián)運(yùn)協(xié)同優(yōu)化模型。但這些研究中加入的不確定因素較為單一，且均為單目標(biāo)優(yōu)化模型。

總體來講，目前多式聯(lián)運(yùn)路徑選擇所建立的模型基本為單目標(biāo)優(yōu)化模型，即使考慮多個(gè)因素，也是將多個(gè)因素轉(zhuǎn)化為單目標(biāo)規(guī)劃問題進(jìn)行處理求解，少數(shù)文章是建立雙目標(biāo)優(yōu)化模型求解，而對含有不確定信息的多式聯(lián)運(yùn)研究則多以單一不確定因素及單目標(biāo)函數(shù)為主。對于多式聯(lián)運(yùn)的發(fā)展趨勢，SteadieSeifi 等[19]指出在多式聯(lián)運(yùn)中多目標(biāo)交通規(guī)劃值得更多研究，需要考慮整合重置資源及同時(shí)規(guī)劃多個(gè)資源，動態(tài)性和數(shù)據(jù)隨機(jī)性也是研究的重要挑戰(zhàn)。Mathisen等[20]指出一些關(guān)于環(huán)境等問題越來越受關(guān)注，如低碳、可持續(xù)發(fā)展。因此，本文站在多式聯(lián)運(yùn)經(jīng)營人的角度，充分考慮影響集裝箱多式聯(lián)運(yùn)過程路徑和運(yùn)輸方式選擇的多個(gè)因素。同時(shí)，由于交通工具的機(jī)械故障、道路檢修、天氣因素，或者突發(fā)交通事故等種種原因，易造成運(yùn)輸時(shí)間和碳排放量不確定，故本文構(gòu)建了包含不確定的運(yùn)輸時(shí)間、碳排放量的多目標(biāo)優(yōu)化模型，使模型盡可能接近真實(shí)情況。然后引入含有多個(gè)隨機(jī)變量的多目標(biāo)粒子群算法進(jìn)行求解，以避免模型在變換過程中失真，得出整個(gè)運(yùn)輸過程可選擇的路徑和每條路徑上的運(yùn)輸方式，以及整個(gè)運(yùn)輸過程所花費(fèi)的總成本、消耗的總時(shí)間和產(chǎn)生的二氧化碳總排放量，從而便于決策者根據(jù)實(shí)際運(yùn)輸情況做出合理的運(yùn)輸決策。

1 模型構(gòu)建

1.1 問題描述

集裝箱多式聯(lián)運(yùn)過程中，貨物從起始地到最終目的地間會經(jīng)過多個(gè)轉(zhuǎn)運(yùn)節(jié)點(diǎn)，在每兩個(gè)轉(zhuǎn)運(yùn)節(jié)點(diǎn)間可以選擇通過卡車、火車和船舶這三種運(yùn)輸方式來運(yùn)輸大量的集裝箱貨物，在每個(gè)轉(zhuǎn)運(yùn)節(jié)點(diǎn)處可以選擇更換另一種運(yùn)輸方式。

具體的運(yùn)輸網(wǎng)絡(luò)結(jié)構(gòu)如圖1所示。Pi代表運(yùn)輸節(jié)點(diǎn)，i=1,2,…,n，Mg代表兩節(jié)點(diǎn)之間所選擇的運(yùn)輸模式，集裝箱多式聯(lián)運(yùn)過程中有三種運(yùn)輸方式可選，所以通常g=1,2,3，分別代表公路運(yùn)輸、鐵路運(yùn)輸、水路運(yùn)輸三種運(yùn)輸方式。

圖1 多式聯(lián)運(yùn)網(wǎng)絡(luò)規(guī)劃結(jié)構(gòu)圖

本文以集裝箱多式聯(lián)運(yùn)過程中的成本、時(shí)間、碳排放量為三個(gè)優(yōu)化目標(biāo)，每個(gè)目標(biāo)都由兩部分組成：一是產(chǎn)生于每兩個(gè)轉(zhuǎn)運(yùn)節(jié)點(diǎn)間的運(yùn)輸過程；二是產(chǎn)生于在轉(zhuǎn)運(yùn)節(jié)點(diǎn)轉(zhuǎn)換運(yùn)輸方式的轉(zhuǎn)換過程。此外，本文考慮了運(yùn)輸過程中時(shí)間和碳排放量的不確定性。根據(jù)不確定性因素產(chǎn)生機(jī)理，整個(gè)運(yùn)輸過程中的時(shí)間、碳排放量等參數(shù)的不確定性主要是隨機(jī)性，可用概率論、數(shù)理統(tǒng)計(jì)或隨即過程理論等來處理，本文采用隨機(jī)變量來處理運(yùn)輸過程中的時(shí)間和碳排放量的不確定性。

1.2 問題假設(shè)

(1) 集裝箱貨物在多式聯(lián)運(yùn)過程中不拆箱、不拼箱。

(2) 不考慮集裝箱貨物類型以及運(yùn)輸方式的車型、船型。

(3) 在任意兩個(gè)存在路徑的轉(zhuǎn)運(yùn)節(jié)點(diǎn)之間，只能選擇一種運(yùn)輸方式。

(4) 只能在節(jié)點(diǎn)處更換運(yùn)輸方式，且在一個(gè)節(jié)點(diǎn)只能更換一次運(yùn)輸方式。

(5) 在轉(zhuǎn)運(yùn)節(jié)點(diǎn)處不考慮更換運(yùn)輸方式時(shí)的設(shè)備條件限制及容量限制。

(6) 假設(shè)運(yùn)輸總能滿足貨物最終送達(dá)的時(shí)間要求。

(7) 假設(shè)根據(jù)調(diào)查，城市之間相應(yīng)的運(yùn)輸時(shí)間、碳排放量，以及不同運(yùn)輸方式之間的轉(zhuǎn)換時(shí)間和轉(zhuǎn)換的碳排放量符合正態(tài)分布。

1.3 符號說明

1.4 模型建立

minZ=(minc,mint,mine)T

(1)

(2)

(3)

(4)

(5)

(6)

(7)

(8)

i,j,k∈Pa,b∈M

(9)

式(1)為集裝箱多式聯(lián)運(yùn)過程中所需優(yōu)化的多目標(biāo)函數(shù)，包括：運(yùn)輸成本最小，運(yùn)輸時(shí)間最小，以及運(yùn)輸過程中碳排放量最少。式(2)-式(4)分別為整個(gè)運(yùn)輸過程的運(yùn)輸成本、運(yùn)輸時(shí)間、運(yùn)輸過程中碳排放量的表達(dá)式，包括在途運(yùn)輸時(shí)和中轉(zhuǎn)節(jié)點(diǎn)更換運(yùn)輸方式時(shí)所產(chǎn)生的成本、花費(fèi)的時(shí)間及產(chǎn)生的碳排放量。式(5)代表每兩個(gè)節(jié)點(diǎn)間只能選擇一種運(yùn)輸方式。式(6)代表在每個(gè)轉(zhuǎn)運(yùn)節(jié)點(diǎn)只能發(fā)生一次運(yùn)輸方式的轉(zhuǎn)變。式(7)保證了在轉(zhuǎn)運(yùn)節(jié)點(diǎn)j處運(yùn)輸方式轉(zhuǎn)換的連續(xù)性。式(8)為0-1變量約束。式(9)中i、j、k屬于節(jié)點(diǎn)集合，a、b屬于運(yùn)輸方式集合。

2 算法流程

在Coello Coello和Lechuga[21]提出的多目標(biāo)粒子群算法(簡稱CMOPSO)中，外部粒子群用來指導(dǎo)該群體外其他粒子的飛行，將Pareto檔案進(jìn)化策略(PAES[22])中的自適應(yīng)網(wǎng)格法用于外部粒子群的維護(hù)。根據(jù)CMOPSO，算法流程如圖2所示。

圖2 多式聯(lián)運(yùn)多目標(biāo)粒子群算法流程圖

其中，所有粒子的位置和速度的更新公式如下：

(10)

(11)

3 算例分析

假設(shè)一次集裝箱多式聯(lián)運(yùn)過程中會經(jīng)過5個(gè)中轉(zhuǎn)點(diǎn)，用數(shù)字1～5代表這5個(gè)中轉(zhuǎn)點(diǎn)，以0和6分別代表起始地和目的地，則以0開頭和6結(jié)尾的一個(gè)數(shù)列即為一條路徑，如0→2→1→5→3→4→6。那么這段路徑中共有6段路程，每段路程需要選擇運(yùn)輸方式，用1～3分別代表公路、鐵路、水路這三種運(yùn)輸方式，由這3個(gè)數(shù)字隨機(jī)構(gòu)成的一個(gè)6個(gè)元素的數(shù)組即對應(yīng)這條路徑中每段路程所采用的運(yùn)輸方式，如{1,3,1,2,3,1}代表在0→2這段路程中選擇1即公路運(yùn)輸。

如圖3所示，0為起始地，6為目的地，現(xiàn)在要運(yùn)輸8個(gè)單位的貨物。

圖3 交通網(wǎng)絡(luò)拓?fù)鋱D

兩個(gè)城市之間的運(yùn)輸方式所對應(yīng)的單位運(yùn)輸成本如表1所示，兩兩城市之間對應(yīng)的符合正態(tài)分布的隨機(jī)運(yùn)輸時(shí)間和隨機(jī)單位碳排放量表2和表3所示。當(dāng)兩個(gè)城市之間不存在某種運(yùn)輸方式時(shí)，其對應(yīng)的運(yùn)輸成本、運(yùn)輸時(shí)間和碳排放量的值選取無窮大，在表中用“-”表示。

表1 城市之間相應(yīng)運(yùn)輸方式的單位運(yùn)輸成本 元

表2 城市之間相應(yīng)運(yùn)輸方式的隨機(jī)單位運(yùn)輸時(shí)間的均值和方差

表3 城市之間相應(yīng)運(yùn)輸方式的隨機(jī)單位碳排放量的均值和方差

在某一城市進(jìn)行中轉(zhuǎn)時(shí)，不同運(yùn)輸方式間轉(zhuǎn)換所產(chǎn)生的費(fèi)用、花費(fèi)的時(shí)間以及產(chǎn)生的碳排放量如表4所示，其中轉(zhuǎn)換時(shí)間和碳排放量的值也是符合正態(tài)分布的隨機(jī)變量的均值和方差。

表4 不同運(yùn)輸方式之間轉(zhuǎn)換費(fèi)用、隨機(jī)轉(zhuǎn)換時(shí)間及隨機(jī)碳排放量

迭代曲線不代表任何非支配解結(jié)果，而是表示迭代過程，代表粒子群的結(jié)果，曲線不斷下降, 說明算法有效, 曲線下降速度較快，說明算法收斂快，效果好，最后曲線平了說明算法已經(jīng)收斂。圖4-圖6為各種迭代曲線。

圖4 運(yùn)輸成本的迭代曲線

圖5 運(yùn)輸時(shí)間的迭代曲線

圖6 碳排放量的迭代曲線

最終共得到如表5的12個(gè)運(yùn)輸方案，多式聯(lián)運(yùn)經(jīng)營人可根據(jù)不同的運(yùn)輸要求來進(jìn)行選擇。

表5 運(yùn)輸方案

續(xù)表5

4 結(jié) 語

本文探究了以運(yùn)輸成本、不確定的運(yùn)輸時(shí)間、不確定的碳排放量為目標(biāo)的集裝箱多式聯(lián)運(yùn)多目標(biāo)優(yōu)化問題。運(yùn)用含隨機(jī)變量的多目標(biāo)粒子群算法進(jìn)行求解，得出集裝箱多式聯(lián)運(yùn)整個(gè)運(yùn)輸過程可供選擇的運(yùn)輸路徑和每段路程所選擇的運(yùn)輸方式，并計(jì)算每種運(yùn)輸方案所花費(fèi)的運(yùn)輸成本、運(yùn)輸時(shí)間和造成的碳排放量，多式聯(lián)運(yùn)經(jīng)營人可根據(jù)不同的偏好和目標(biāo)來選擇運(yùn)輸方案。下一步研究可考慮時(shí)間窗限制，此外，還需對多式聯(lián)運(yùn)過程中的每段運(yùn)輸路程中每種運(yùn)輸方式的運(yùn)輸時(shí)間、碳排放量做深入調(diào)研，從而確定它們具體符合那種隨機(jī)變量的數(shù)字特征。在更為成熟的研究后，或可將此算法設(shè)計(jì)一個(gè)小程序，輸入各個(gè)節(jié)點(diǎn)的坐標(biāo)及各個(gè)節(jié)點(diǎn)間各種運(yùn)輸方式的運(yùn)輸成本、運(yùn)輸時(shí)間和碳排放量后，可以通過封裝起來的算法直接輸出可選擇的運(yùn)輸方案，以供集裝箱多式聯(lián)運(yùn)經(jīng)營人選擇。

參考文獻(xiàn)

[1] 佟璐, 聶磊, 付慧伶. 多式聯(lián)運(yùn)路徑優(yōu)化模型與方法研究[J]. 物流技術(shù), 2010, 29(5):57-60.

[2] 李雙琳, 馬祖軍, 鄭斌,等. 震后初期應(yīng)急物資配送的模糊多目標(biāo)選址-多式聯(lián)運(yùn)問題[J]. 中國管理科學(xué), 2013, 21(2):144-151.

[3] 楊秋秋, 王輝. 基于遺傳算法的多目標(biāo)集裝箱多式聯(lián)運(yùn)運(yùn)輸優(yōu)化模型[J]. 物流科技, 2006, 29(12):29-31.

[4] 何艷梅, 何俊生. 多目標(biāo)多式聯(lián)運(yùn)配送路徑研究[J]. 物流科技, 2013, 36(5):112-114.

[5] 黃麗霞, 帥斌. 危險(xiǎn)貨物多式聯(lián)運(yùn)路徑優(yōu)化問題的多目標(biāo)優(yōu)化算法[J]. 中國安全生產(chǎn)科學(xué)技術(shù), 2014, 10(9):10-16.

[6] Bhattacharya A, Kumar S A, Tiwari M K, et al. An intermodal freight transport system for optimal supply chain logistics[J]. Transportation Research Part C Emerging Technologies, 2014, 38(1):73-84.

[7] Modesti P, Sciomachen A. A utility measure for finding multiobjective shortest paths in urban multimodal transportation networks[J]. European Journal of Operational Research, 1998, 111(3):495-508.

[8] Sun Y, Lang M, Wang D. Optimization Models and Solution Algorithms for Freight Routing Planning Problem in the Multi-Modal Transportation Networks: A Review of the State-of-the-Art[J]. Open Civil Engineering Journal, 2015, 9(1):714-723.

[9] Sun Y, Lang M. Bi-objective optimization for multi-modal transportation routing planning problem based on Pareto optimality[J]. Journal of Industrial Engineering & Management, 2015, 8(4): 1195-1217.

[10] Kengpol A, Tuammee S, Tuominen M. The development of a framework for route selection in multimodal transportation[J]. International Journal of Logistics Management, 2014, 25(3):581-610.

[11] 賀竹磬, 孫林巖, 李曉宏. 時(shí)效性物流聯(lián)運(yùn)方式選擇模型及其算法[J]. 管理科學(xué), 2007, 20(1):8-12.

[12] 陳怡霏. 碳約束政策下的集裝箱多式聯(lián)運(yùn)路徑優(yōu)化研究[D]. 北京交通大學(xué), 2014.

[13] 李高波. 具有容量約束的混合軸輻式多式聯(lián)運(yùn)樞紐網(wǎng)絡(luò)設(shè)計(jì)[D]. 長安大學(xué), 2014.

[14] 付曉鳳, 馬彬, 張娟,等. 多目標(biāo)一體化的聯(lián)運(yùn)路徑優(yōu)化方法研究[J]. 鐵道運(yùn)輸與經(jīng)濟(jì), 2009, 31(9):83-85.

[15] 王正彬, 鄭康立. 集裝箱多式聯(lián)運(yùn)組合優(yōu)化研究[J]. 大連交通大學(xué)學(xué)報(bào), 2012, 33(5):63-66.

[16] Wang Q B,Han Z,Ji M,et al.Path Optimization of Container Multimodal Transportation Based on Node Operation Randomness [J].Journal of Transportation Systems Engineering and Information Technology,2011,11(6):137-144.

[17] 王慧, 汪傳旭. 模糊需求環(huán)境下集裝箱多式聯(lián)運(yùn)箱型和運(yùn)輸方式的選擇[J]. 公路交通科技, 2012, 29(4):153-158.

[18] 張得志, 李雙艷. 不確定環(huán)境下協(xié)同運(yùn)輸優(yōu)化模型及其求解算法[J]. 鐵道科學(xué)與工程學(xué)報(bào), 2010, 7(4):116-120.

[19] Steadieseifi M, Dellaert N P, Nuijten W, et al. Multimodal freight transportation planning: A literature review[J]. European Journal of Operational Research, 2014, 233(1):1-15.

[20] Mathisen T A, Hanssen T E S. The Academic Literature on Intermodal Freight Transport [J]. Transportation Research Procedia, 2014, 3:611-620.

[21] Coello C C A, Lechuga M S. MOPSO: a proposal for multiple objective particle swarm optimization[C]// Evolutionary Computation, 2002. CEC ’02. Proceedings of the 2002 Congress on. IEEE Xplore, 2002:1051-1056.

[22] Knowles J D, Corne D W. Approximating the Nondominated Front Using the Pareto Archived Evolution Strategy[J]. Evolutionary Computation, 1993, 8(2):149.