基于廣義費用函數(shù)的多式聯(lián)運路徑優(yōu)化研究

戶佐安 賈葉子 孫 燕 薛 鋒* 趙 寶

(西南交通大學交通運輸與物流學院1) 成都 611756) (西南交通大學綜合交通大數(shù)據(jù)應(yīng)用技術(shù)國家工程實驗室2) 成都 611756) (大秦鐵路股份有限公司科學技術(shù)研究所3) 太原 030013)

0 引 言

多式聯(lián)運參與的主體眾多，托運人、承運人及運輸管理者等決策主體均出于自身利益考慮對運輸組織方案有不同的需求，進而影響多式聯(lián)運的運輸路徑和運輸方式選擇.因此，在多式聯(lián)運路徑優(yōu)化問題中綜合考慮不同決策主體的目標，有利于整合社會資源、提高多式聯(lián)運效率，實現(xiàn)整個社會綜合效益最大化.

國內(nèi)外學者在多式聯(lián)運路徑優(yōu)化方面已有諸多研究.Qu等[1]構(gòu)建了考慮中轉(zhuǎn)成本和碳排放成本的多商品流多式聯(lián)運整數(shù)非線性規(guī)劃模型；Athanasios等[2]構(gòu)建了以運輸時間最短為目標的模型，并設(shè)計相關(guān)算法得到最優(yōu)運輸方案；Min等[3]在研究多式聯(lián)運運輸方式選擇問題時，以成本、時間和風險為目標，構(gòu)建了不考慮運輸能力限制的優(yōu)化模型；劉璘等[4]在研究海鐵聯(lián)運路徑優(yōu)化問題時，將總成本劃分為制冷成本、運輸成本和轉(zhuǎn)運成本，并構(gòu)建了以總成本最小為目標的冷藏集裝箱路徑優(yōu)化模型；李珺等[5]在研究不確定條件下的多式聯(lián)運路徑優(yōu)化問題時，考慮了運輸成本、碳排放成本和時間懲罰成本等因素對運輸方案的影響，并構(gòu)建綠色多式聯(lián)運路徑優(yōu)化模型；彭勇等[6]在構(gòu)建模型時，除了考慮運輸成本、運輸時間等常見因素外，還增加了對運輸風險和班期限制的研究；魏航等[7]考慮了危險物品在運輸途中可能對運輸路線周邊人口產(chǎn)生危險，將危險物品途經(jīng)處的人口覆蓋率以及成本、風險共同作為模型目標；甄遠迪等[8]在構(gòu)建不確定情況下的多式聯(lián)運路徑優(yōu)化模型時，選擇的優(yōu)化目標包括成本、時間、碳排放量及轉(zhuǎn)運節(jié)點能力；雷定猷等[9]在研究長大貨物多式聯(lián)運路徑優(yōu)化問題時，在模型構(gòu)建中增加了適用于長大貨物運輸?shù)募s束條件.

綜上，現(xiàn)有關(guān)于多式聯(lián)運路徑優(yōu)化方面的研究較多，但針對于多式聯(lián)運全過程構(gòu)建廣義費用函數(shù)的相關(guān)研究較少，且多是從運輸客戶需求角度，而從交通管理部門出發(fā)的研究較少，所考慮的因素不夠全面.文中針對當前研究存在的問題，在已有研究成果的基礎(chǔ)上，從托運人、承運人和運輸管理者等不同決策主體的角度出發(fā)，采用廣義費用函數(shù)體現(xiàn)多式聯(lián)運對成本、時效性、安全性和環(huán)保等多方面的現(xiàn)實要求，構(gòu)建基于廣義費用函數(shù)的多式聯(lián)運路徑優(yōu)化模型，設(shè)計算例進行優(yōu)化，為不同決策主體提供參考意見.

1 多式聯(lián)運廣義費用函數(shù)構(gòu)建

1.1 托運人目標分析

1) 運輸成本

(1)

2) 運輸時間價值成本 貨物運輸總時間見式(2),運輸時間價值成本優(yōu)化目標見式(3).

(2)

式中：CT為貨物運輸時間價值成本，元；Vmax為單位集裝箱滿載時的貨物價值，元；η為貨物入箱率，%；ρ為年貨物持有成本系數(shù)，%；Vr為單位集裝箱單位時間租賃費用，元/(箱·d).

1.2 承運人目標分析

1) 運輸安全性

(4)

2) 運輸準時性

minf(tD)=pE·max(ET-tD,0)+

pL·max(tD-LT,0)

(5)

式中：f(tD)為時間偏離懲罰函數(shù)；tD為貨物到達目的地的時間；pE為貨物因提前到達需等待的單位時間滯留成本，元；pL為貨物因延遲到達需支付的單位時間懲罰成本，元；ET為客戶要求時間窗的最早時間；LT為客戶要求時間窗的最晚時間.

3) 運輸可靠性

(6)

1.3 運輸管理者目標分析

1) 大氣污染成本

(7)

2) 噪聲污染成本

enk=φ·enuk·ψg·ψp

(8)

式中：enk為我國第k種運輸方式的單位噪聲污染成本，元/(t·km)；enuk為歐洲第k種運輸方式的單位噪聲污染成本，歐元/(t·km)；φ為歐元兌人民幣平均匯率；ψg為人均GDP調(diào)整系數(shù)；ψp為人口密度調(diào)整系數(shù).

3) 交通事故成本

(9)

根據(jù)該地區(qū)該種運輸方式的貨物周轉(zhuǎn)量，可計算得單位交通事故外部成本，為

(10)

式中：uck為第k種運輸方式的單位交通事故外部成本；TCk為第k種運輸方式的交通事故總外部成本；TQk為第k種運輸方式的總貨運周轉(zhuǎn)量.

由上述分析可知，多式聯(lián)運廣義費用函數(shù)由運輸成本、運輸時間成本、運輸安全性、運輸準時性、運輸可靠性和運輸外部成本六部分組成，即

s.t.λ1+λ2+λ3+λ4=1,λ1,λ2,λ3,λ4>0

In order to translate Zhuang drama into English and maintain the faithfulness of cultural transmission,understanding the cultural characteristicsof Zhuangdramawill beutmost priority.

(11)

式中：GC為廣義費用函數(shù)；λi為權(quán)重系數(shù)，反映各影響因素在廣義費用函數(shù)中所占比重.

2 模型構(gòu)建

2.1 問題描述與模型假設(shè)

某貨主有一批貨物需要從起始地O運送到目的地D，貨物從起始地到目的地會經(jīng)過若干個轉(zhuǎn)運節(jié)點，每兩個相互連通的節(jié)點之間可通過鐵路、公路、水路三種運輸方式進行運輸，每種運輸方式的運輸成本和運輸時間不同.每個轉(zhuǎn)運節(jié)點處可選擇一種運輸方式轉(zhuǎn)換，由此產(chǎn)生相應(yīng)的轉(zhuǎn)運時間和轉(zhuǎn)運成本.在以上條件下，以廣義費用函數(shù)最小為目標，確定合理的運輸方案.

為貼近實際情況并方便對模型求解，提出以下假設(shè)：①同一批貨物在運輸及中轉(zhuǎn)過程中不可拆分；②不考慮不同類別貨物間的拼箱運輸；③集裝箱在運輸及中轉(zhuǎn)過程中的箱型均為TEU標準箱；④任意兩個相鄰節(jié)點之間最多只能選擇一條路徑，每條路徑只能選擇一種運輸方式；⑤不同運輸方式的轉(zhuǎn)換只能發(fā)生在節(jié)點處且每個節(jié)點最多只能完成一次中轉(zhuǎn)；⑥不同運輸方式的運輸能力均能滿足托運人貨運量的要求.

2.2 模型建立

對相關(guān)集合、參數(shù)的說明見表1.

表1 模型相關(guān)集合、參數(shù)說明

以廣義費用函數(shù)最小為目標，建立模型為

λ3·(pE·max(ET-tD,0)+pL·max(tD-LT,0))+

(12)

s.t.

(13)

(14)

(15)

(16)

(17)

(18)

(19)

2.3 模型求解

本文建立的模型為多目標混合整數(shù)規(guī)劃模型，由于各單目標量綱不同，無法直接采用線性加權(quán)法進行求解.理想點法的基本思想是先構(gòu)造滿足所有目標的理想點，之后在約束條件下尋找與該理想點最接近的可行解并以最小偏差量作為目標函數(shù)進行求解[10]，在多目標決策和多目標優(yōu)化中已有廣泛應(yīng)用.本文選擇該方法進行求解，具體步驟如下.

本文采用二項分布賦權(quán)法對含有n個目標的多目標模型設(shè)定權(quán)重，可用概率公式對二項分布描述為

k=0,1,2,…,mp∈(0,1)

設(shè)權(quán)重組合為r組，則各組權(quán)重對應(yīng)p=(c+1)/(r+1)，其中c可取0,1,2,…,r-1.

通過上述步驟，多目標規(guī)劃問題轉(zhuǎn)化為以加權(quán)歐氏距離最小為目標的單目標規(guī)劃問題，借助優(yōu)化求解器LINGO可進行求解.

3 算例分析

3.1 算例描述

某貨主現(xiàn)有一批360 t的集裝箱貨物，采用TEU標準集裝箱從起始節(jié)點運輸?shù)侥康墓?jié)點，假定集裝箱滿載量為24 t，則所需要標箱數(shù)量為15 TEU.運輸途中有8個集裝箱中轉(zhuǎn)站可進行運輸方式轉(zhuǎn)換.

表2 集裝箱各運輸方式運價核算表

表3 不同運輸方式之間單位中轉(zhuǎn)費用 單位：元/TEU

表4 各節(jié)點運輸距離 單位：km

表5 不同運輸方式之間單位節(jié)點換裝時間 單位：h/TEU

表6 各節(jié)點運輸方式班期時刻表

鐵路運輸速度為55 km/h，公路運輸速度為65 km/h，水路運輸速度為30 km/h.單位集裝箱滿載運輸價值為60萬元，貨物入箱率為100%，年貨物持有成本系數(shù)為25%，單位集裝箱單位時間租賃費用為100元/(TEU·d).每百公里條件下，公路運輸貨損率為0.01%，鐵路運輸為0.015%，水路運輸為0.02%，轉(zhuǎn)運節(jié)點處貨損率為0.04%/次.貨物在凌晨零點出發(fā)，客戶要求的軟時間窗為[15 h,25 h]，因提前到達需等待的單位時間滯留成本為40元/h，因延遲到達需支付的單位時間懲罰成本為60元/h，運輸弧上不同運輸方式的各因素可靠度見表7.

大氣污染成本計算中，鐵路運輸以HXD3D型電力機車為例，公路運輸以一汽解放J6P重卡為例，水路運輸以內(nèi)河全集裝箱船舶相關(guān)參數(shù)為例，見表8～14.鐵路運輸中，假定電力機車能源結(jié)構(gòu)中77.7%來源于火力發(fā)電，平均傳輸損耗率為93%，充電損耗率為20%，燃燒1 t標準煤生產(chǎn)的可用電量為8 131 kW·h.

表7 運輸弧上不同運輸方式的各因素可靠度

表8 HXD3D型電力機車相關(guān)參數(shù)表

表9 標準煤排放率 單位：kg/t

表10 一汽解放J6P重卡相關(guān)參數(shù)表

表11 公路運輸污染物功率排放表 單位：g/(kW·h)

表12 內(nèi)河全集裝箱船舶相關(guān)參數(shù)表

表13 水路運輸污染物功率排放表 單位：g/(kW·h)

表14 污染物排放單價表 單位：元/kg

歐元兌人民幣平均匯率為7.783 4(2019年)，人均GDP調(diào)整系數(shù)為0.26，人口密度調(diào)整系數(shù)為1.19，歐洲噪聲污染情況見表15.

表15 歐洲噪聲總成本和平均單位噪聲成本預(yù)測值

根據(jù)我國實際情況，事故死亡1人平均賠付30萬元，受傷1人平均賠付5萬元，公路事故管理救助成本為1 000元/起，鐵路事故為10萬元/起，水路事故為15萬元/起.通過查閱相關(guān)統(tǒng)計年鑒和行業(yè)數(shù)據(jù)，得到不同運輸方式2018年的交通事故統(tǒng)計情況和貨運周轉(zhuǎn)量，見表16.

表16 各運輸方式交通事故相關(guān)數(shù)據(jù)

3.2 算例求解

采用二項分布賦權(quán)法對六個目標設(shè)定權(quán)重，設(shè)權(quán)重組合為10組，得到各組權(quán)重對應(yīng)p=(1/11,2/11,…,9/11,10/11).借助數(shù)學規(guī)劃軟件LINGO12.0對不同權(quán)重組合下的單目標規(guī)劃模型進行求解，得到多組最優(yōu)解下的多式聯(lián)運運輸方案，見表17.

表17 不同權(quán)重組合下的運輸方案

由表17可知，當p=1/11時，運輸外部成本目標權(quán)重取最大，運輸方案為鐵水聯(lián)運，可見公路運輸所產(chǎn)生的外部成本遠大于鐵路和水路運輸，從運輸管理者角度制定運輸方案時，為降低運輸外部成本，應(yīng)盡量選擇鐵水聯(lián)運模式；當p從2/11增加至6/11時，運輸外部成本目標權(quán)重逐步減小，而承運人所關(guān)注的安全性、準時性和可靠性的目標權(quán)重逐步增大，此階段運輸方案為公鐵聯(lián)運且鐵路運輸所占比例遠大于公路運輸，說明公鐵聯(lián)運在時效性方面強于鐵水聯(lián)運，而鐵路運輸?shù)陌踩?、可靠性等指標?yōu)于公路運輸和水路運輸，因此從承運人角度制定運輸方案時，應(yīng)優(yōu)先選擇鐵路運輸和其他運輸方式組合的運輸模式；當p從7/11增加至10/11時，托運人所關(guān)注的運輸成本和運輸時間價值成本的目標權(quán)重不斷增加，而承運人和運輸管理者的目標權(quán)重都有所減小，此階段運輸方案經(jīng)歷了由公鐵聯(lián)運轉(zhuǎn)為全程公路運輸后又轉(zhuǎn)為公鐵聯(lián)運的過程，且公路運輸與鐵路運輸所占比例相近，可看出公鐵聯(lián)運模式可同時滿足時效性和成本兩方面要求，因此從托運人角度制定運輸方案時，從時效性和成本兩方面考慮時應(yīng)盡量選擇公鐵聯(lián)運.

4 結(jié) 束 語

本文考慮了多式聯(lián)運不同決策主體的需求，分別對其目標進行分析和量化處理，建立了包含運輸費用、運輸時間價值成本、時間偏離懲罰成本、運輸外部成本、運輸安全性和運輸可靠性在內(nèi)的多式聯(lián)運廣義費用函數(shù)，并以廣義費用函數(shù)最優(yōu)為目標構(gòu)建了多式聯(lián)運路徑優(yōu)化模型.設(shè)計理想點法并應(yīng)用LINGO12.0對算例進行求解，獲得了不同權(quán)重組合下的多組全局最優(yōu)解，算例求解結(jié)果表明本文所建模型能夠為不同決策主體提供參考依據(jù)，同時能夠較好平衡不同決策主體目標之間的關(guān)系，具有一定的有效性和實用性.為簡化模型，本文對多式聯(lián)運外部成本的核算方法進行了部分簡化，這也是后續(xù)將進行的工作.