劉虹，林楚玥

（福州大學(xué)經(jīng)濟(jì)與管理學(xué)院，福建福州350100）

一、引言

近年來(lái)，國(guó)內(nèi)物流行業(yè)發(fā)展迅猛，第九屆中國(guó)物流信息化大會(huì)以“創(chuàng)新驅(qū)動(dòng)物流產(chǎn)業(yè)升級(jí)”為主題，力求推進(jìn)智慧物流的行業(yè)應(yīng)用。而越庫(kù)（cross-docking）作為一種以供給鏈間信息化為基礎(chǔ)的高效物流新技術(shù)，因其能大幅降低庫(kù)存水平，加快貨物流動(dòng)，在國(guó)外已被廣泛應(yīng)用于配送網(wǎng)絡(luò)中[1]。

越庫(kù)問(wèn)題的研究主要有內(nèi)部運(yùn)作調(diào)度和外部車(chē)輛路徑兩方面，其中車(chē)輛路徑研究考慮了協(xié)調(diào)運(yùn)送、運(yùn)輸策略等方面。Lee 等學(xué)者首次提出越庫(kù)的車(chē)輛路徑問(wèn)題，取貨車(chē)輛應(yīng)同時(shí)到達(dá)越庫(kù)中心，以便于合并處理[2]；Liao 等學(xué)者對(duì)此方法進(jìn)行了改進(jìn)[3]。Wen 等學(xué)者放寬了車(chē)輛在越庫(kù)中同步到達(dá)的假設(shè)，通過(guò)整合來(lái)協(xié)調(diào)運(yùn)送路線[4]。Yu 等學(xué)者研究越庫(kù)的開(kāi)放式車(chē)輛路徑問(wèn)題，協(xié)同取貨車(chē)輛的到庫(kù)時(shí)間[5]。曹克官討論了直送和循環(huán)取貨兩種調(diào)度策略下的越庫(kù)配送網(wǎng)絡(luò)問(wèn)題[6]。葛顯龍等學(xué)者在循環(huán)取送貨策略下，研究供應(yīng)鏈環(huán)境下的越庫(kù)協(xié)同配送網(wǎng)絡(luò)最小成本問(wèn)題[7]。而B(niǎo)aniamerian 等學(xué)者則關(guān)注越庫(kù)配送網(wǎng)絡(luò)的最大利潤(rùn)問(wèn)題[8]。Maknoon 和Laporte 在多越庫(kù)的配送網(wǎng)絡(luò)中考慮車(chē)輛負(fù)載同步性，研究貨物可在越庫(kù)間轉(zhuǎn)移合并策略下的車(chē)輛路徑問(wèn)題，并提出自適應(yīng)大領(lǐng)域搜索算法[9]。

越庫(kù)配送網(wǎng)絡(luò)注重整體的配送效率，而配送效率的高低直接影響著客戶(hù)服務(wù)水平，衡量配送服務(wù)水平的研究，主要集中在建立與時(shí)間相關(guān)的滿(mǎn)意度函數(shù)來(lái)評(píng)價(jià)服務(wù)水平，滿(mǎn)意度越高則服務(wù)水平越好。Fan 認(rèn)為配送服務(wù)水平與客戶(hù)等待時(shí)間呈現(xiàn)反比關(guān)系，由此構(gòu)建滿(mǎn)意度函數(shù)[10]。Chen、樓振凱等學(xué)者用梯形模糊數(shù)表征模糊時(shí)間窗，用模糊度隸屬函數(shù)刻畫(huà)客戶(hù)滿(mǎn)意度，通過(guò)線性分段式函數(shù)體現(xiàn)客戶(hù)對(duì)配送服務(wù)的滿(mǎn)意程度[11-12]。王君同樣是用模糊時(shí)間窗反映客戶(hù)對(duì)時(shí)間的偏好程度，并將客戶(hù)滿(mǎn)意度表示為關(guān)于時(shí)間的凹模糊數(shù)[13]。任亮等學(xué)者考慮客戶(hù)拖期厭惡行為，運(yùn)用前景理論價(jià)值函數(shù)體現(xiàn)客戶(hù)配送服務(wù)水平[14]。

綜上所述，越庫(kù)配送的研究主要集中于車(chē)輛調(diào)度、協(xié)同到庫(kù)策略等，考慮分揀整合策略的越庫(kù)配送相對(duì)較少。另一方面，由于存在交通阻塞等一些不確定因素，可能使產(chǎn)品送達(dá)時(shí)間不確定，進(jìn)而使產(chǎn)品送達(dá)時(shí)間窗偏離，也會(huì)使客戶(hù)的心理感受隨之變化；對(duì)于越庫(kù)配送中心而言，產(chǎn)品配送時(shí)間的不確定對(duì)車(chē)輛和駕駛員的調(diào)度均會(huì)產(chǎn)生影響。本文主要研究的問(wèn)題是：在連鎖零售的越庫(kù)配送網(wǎng)絡(luò)中，針對(duì)越庫(kù)中心內(nèi)部的協(xié)同到庫(kù)、分揀整合和外部車(chē)輛路徑的集成配送，考慮客戶(hù)對(duì)配送時(shí)間的心理感知對(duì)服務(wù)水平評(píng)價(jià)存在影響，引入前景理論，構(gòu)建客戶(hù)時(shí)間窗服務(wù)價(jià)值函數(shù)，尋求使得時(shí)間窗服務(wù)價(jià)值最高和配送成本最低的越庫(kù)車(chē)輛路徑優(yōu)化方案。

二、問(wèn)題描述及建模

（一）問(wèn)題描述與假設(shè)

連鎖零售的越庫(kù)配送網(wǎng)絡(luò)，包括：多個(gè)供應(yīng)商、一個(gè)越庫(kù)配送中心和多個(gè)客戶(hù)零售商。零售商依據(jù)各自的銷(xiāo)售需求向各供應(yīng)商提交貨物訂單，越庫(kù)配送系統(tǒng)由此獲取客戶(hù)訂單，其中客戶(hù)訂單包含三個(gè)屬性信息：取貨地點(diǎn)（即供應(yīng)商節(jié)點(diǎn)）、送貨地點(diǎn)（即零售商節(jié)點(diǎn)）和訂單貨物量。配送網(wǎng)絡(luò)見(jiàn)圖1。

根據(jù)越庫(kù)協(xié)同到庫(kù)和分揀整合策略，調(diào)度取貨車(chē)輛前往取貨點(diǎn)拾取貨物，并協(xié)同抵達(dá)越庫(kù)中心；訂單貨物根據(jù)目的地進(jìn)行分揀重裝；經(jīng)過(guò)整合的貨物由送貨車(chē)輛按一定順序送達(dá)客戶(hù)點(diǎn)，送貨車(chē)輛最后返回越庫(kù)中心?？蛻?hù)根據(jù)貨物送達(dá)時(shí)間會(huì)產(chǎn)生不同的心理感知，影響著時(shí)間窗服務(wù)價(jià)值，越庫(kù)配送系統(tǒng)需要優(yōu)化車(chē)輛配送路徑，以期在降低配送成本的同時(shí)提高客戶(hù)時(shí)間窗服務(wù)價(jià)值。有以下假設(shè)：

假設(shè)1客戶(hù)訂單的三個(gè)屬性信息（取貨點(diǎn)、送貨點(diǎn)、貨物量）以及配送時(shí)間窗已知；

假設(shè)2每個(gè)取/送貨點(diǎn)僅由一輛取/送貨車(chē)輛服務(wù)，并且僅被服務(wù)一次；

假設(shè)3取/送貨車(chē)輛均從越庫(kù)中心出發(fā)，在訪問(wèn)完取/送貨點(diǎn)之后必須返回越庫(kù)中心；

假設(shè)4取貨車(chē)輛采用整車(chē)直運(yùn)方式運(yùn)輸，即每個(gè)取貨點(diǎn)只需一輛車(chē)，且不考慮載重約束；

假設(shè)5送貨車(chē)輛具有載重約束，要求每輛車(chē)所分配的所有訂單貨物量之和不超過(guò)其最大載重量；

假設(shè)6節(jié)點(diǎn)間車(chē)輛運(yùn)輸成本和運(yùn)輸時(shí)間、貨物的單位裝卸時(shí)間和單位處理時(shí)間以及可用車(chē)輛數(shù)量已知。

圖1：越庫(kù)配送網(wǎng)絡(luò)

（二）符號(hào)與變量說(shuō)明

1.模型參數(shù)

R{ r |r= 1, … l}：表示訂單集合；

D{i | i= 1, … n}：表示送貨點(diǎn)（零售商）集合；

P{ j |j= 1,…m }：表示取貨點(diǎn)（供貨商）集合；

O：表示越庫(kù)配送中心；

εjr：表示訂單r (r ∈R)與取貨點(diǎn)j (j ∈P)對(duì)應(yīng)關(guān)系，若 εjr=1，表示訂單r 的取貨點(diǎn)為j；

V{k |k= 1, … v}：表示取貨車(chē)輛集合；

V'{k |k= 1, … v'}：表示取貨車(chē)輛集合；

Q：表示送貨車(chē)輛最大載重量（以托盤(pán)為單位）；

Qv：表示可用車(chē)輛數(shù)量；

S =P ∪｛O｝ ∪ D ：表示配送網(wǎng)絡(luò)中所有節(jié)點(diǎn)集合；

Cv：表示車(chē)輛固定運(yùn)營(yíng)成本；

qr：表示訂單r(r∈R)貨物量（以托盤(pán)為單位）；

di：表示送貨點(diǎn)i(iD∈)供貨總量（以托盤(pán)為單位）；

pj：表示取貨點(diǎn)j(j∈P)收貨總量（以托盤(pán)為單位）；

Tab：表示節(jié)點(diǎn)a 到節(jié)點(diǎn)b 的運(yùn)輸時(shí)間(a, b∈S)；

Cab：表示節(jié)點(diǎn)a 到節(jié)點(diǎn)b 的運(yùn)輸成本(a, b∈S)；

γ ：表示在取/送貨點(diǎn)每一單位貨物所需裝/卸貨時(shí)間；

φ：表示在越庫(kù)中心每一單位貨物所需處理時(shí)間；

M ：表示極大的數(shù)；

[T1,T2]：表示送貨點(diǎn)最早最晚工作時(shí)間窗；

EETi：表示送貨點(diǎn)i(i∈D)期望送貨服務(wù)時(shí)間節(jié)點(diǎn)；

LLTi：表示送貨點(diǎn)i(i∈D)可容忍送貨服務(wù)時(shí)間節(jié)點(diǎn)。

2.決策變量

stk：取貨車(chē)輛 k(k∈V)離開(kāi)越庫(kù)中心出發(fā)取貨的時(shí)間；

atk：取貨車(chē)輛k(k∈V)到達(dá)越庫(kù)中心的時(shí)間；

dtk：送貨車(chē)輛k (k ∈V')離開(kāi)越庫(kù)中心時(shí)間；

tik：送貨車(chē)輛k (k ∈V')到達(dá)送貨點(diǎn)i (i∈D)時(shí)間。

（三）客戶(hù)時(shí)間窗服務(wù)價(jià)值模型

考慮到客戶(hù)在評(píng)價(jià)配送水平時(shí)，往往會(huì)受行為和心理因素的影響，且對(duì)準(zhǔn)時(shí)性要求也并非完全剛性。前景理論價(jià)值函數(shù)，運(yùn)用到服務(wù)評(píng)價(jià)中能夠表示為“收益”與“損失”，結(jié)合模糊時(shí)間窗，構(gòu)建客戶(hù)時(shí)間窗服務(wù)價(jià)值函數(shù)，以此衡量配送服務(wù)水平。

1.價(jià)值函數(shù)的表示

客戶(hù)對(duì)配送到貨時(shí)間產(chǎn)生主觀感受的滿(mǎn)意度或價(jià)值，前景理論可用來(lái)描述客戶(hù)對(duì)時(shí)間窗服務(wù)的不同感知，客戶(hù)i 的感知價(jià)值函數(shù)可表示為如下：

其中：α ,β,λ為參數(shù)，α 表示收益系數(shù)，β 表示損失系數(shù)，λ 表示客戶(hù)對(duì)損失的敏感程度。任何收益或損失都需要根據(jù)參考點(diǎn)做出判斷，一般選取價(jià)值為0 的點(diǎn)作為參考點(diǎn)，在模糊時(shí)間窗中，若送貨時(shí)間在客戶(hù)可容忍時(shí)間之前，則客戶(hù)的時(shí)間窗服務(wù)價(jià)值為正值，可看作是收益；反之則價(jià)值為負(fù)，即損失。因此，本文將客戶(hù)可容忍服務(wù)時(shí)間（LLTi）作為參考點(diǎn)?？紤]到在可容忍時(shí)間端點(diǎn)上是符合配送要求的，取值為0 不符合實(shí)際情況，故假設(shè)在此端點(diǎn)上取值為0.5。

2.背景假設(shè)

為表達(dá)客戶(hù)時(shí)間窗服務(wù)價(jià)值，作如下背景假設(shè)：

假設(shè)7在越庫(kù)配送網(wǎng)絡(luò)中，所有客戶(hù)（零售商）最早和最晚工作時(shí)間分別為 T1、 T2；

假設(shè)8所有客戶(hù)擁有相同的評(píng)價(jià)系數(shù)，即每個(gè)客戶(hù)的α 和β 取值相同；

假設(shè)9客戶(hù)在評(píng)價(jià)配送服務(wù)水平時(shí)，選擇可容忍服務(wù)時(shí)間端點(diǎn)作為參考點(diǎn)，送貨車(chē)輛在此端點(diǎn)到達(dá)，時(shí)間窗服務(wù)價(jià)值取值為0.5；

假設(shè)10當(dāng)送貨車(chē)輛在客戶(hù)i 期望時(shí)間EETi之前到達(dá)客戶(hù)點(diǎn)，時(shí)間窗服務(wù)價(jià)值最高，取值為1；

假設(shè)11當(dāng)送貨車(chē)輛在客戶(hù)i 期望時(shí)間EETi之后，且可容忍時(shí)間LLTi之前到達(dá)客戶(hù)點(diǎn)，客戶(hù)將車(chē)輛到達(dá)時(shí)間與參考點(diǎn)（可容忍時(shí)間端點(diǎn)）進(jìn)行比較，認(rèn)為符合自身時(shí)間要求，時(shí)間窗服務(wù)價(jià)值為正值，可看作收益，應(yīng)選用價(jià)值函數(shù)xi≥ 0部分；

假設(shè)12當(dāng)送貨車(chē)輛在客戶(hù)i 可容忍時(shí)間LLTi之后，且在最晚工作時(shí)間前到達(dá)客戶(hù)點(diǎn)，客戶(hù)接受服務(wù)但認(rèn)為不符合自身時(shí)間要求，時(shí)間窗服務(wù)價(jià)值為負(fù)值，可看作損失，應(yīng)選用xi＜ 0部分。

假設(shè)13當(dāng)送貨車(chē)輛在最晚工作時(shí)間端點(diǎn)到達(dá)客戶(hù)點(diǎn)時(shí)，時(shí)間窗服務(wù)價(jià)值取值為-1，而在2T 之后，客戶(hù)將不再接受服務(wù)。

3.客戶(hù)時(shí)間窗服務(wù)價(jià)值

定理1在越庫(kù)配送網(wǎng)絡(luò)中，若送貨車(chē)輛k 服務(wù)客戶(hù)點(diǎn)i，到達(dá)時(shí)間為 tik，在可容忍時(shí)間之內(nèi)（ [T1,EETi]和(EETi,LLTi]兩個(gè)區(qū)間），則時(shí)間窗服務(wù)價(jià)值函數(shù)U (tik)表達(dá)如下：

證明：由假設(shè)11 可知，客戶(hù)將送貨車(chē)輛到達(dá)時(shí)間與參考點(diǎn)（可容忍時(shí)間端點(diǎn)）進(jìn)行比較，車(chē)輛到達(dá)時(shí)間若在可容忍時(shí)間之前，則認(rèn)為符合自身時(shí)間要求，即時(shí)間窗服務(wù)價(jià)值為正，可看作收益，應(yīng)選用價(jià)值函數(shù)xi≥ 0部分，即

（1）當(dāng)tik∈[T1,EETi]，由假設(shè)10 可知，在期望時(shí)間內(nèi)，客戶(hù)時(shí)間窗服務(wù)價(jià)值最高，數(shù)值恒定為1，即：

（2）當(dāng)tik∈(EETi,LLTi]，設(shè)時(shí)間窗服務(wù)價(jià)值函數(shù) U (tik)為：

由假設(shè)10 和假設(shè)11 可知送貨車(chē)輛到達(dá)時(shí)間為EETi時(shí)，客戶(hù)時(shí)間窗服務(wù)價(jià)值數(shù)值為1，而LLTi時(shí)，時(shí)間窗服務(wù)價(jià)值數(shù)值為1/2，即：

將式（6）除以（5），得：

將1a 代入式（5），得：

即：

綜上，式（2）得證。

定理2在越庫(kù)配送網(wǎng)絡(luò)中，若送貨車(chē)輛k 服務(wù)客戶(hù)點(diǎn)i，到達(dá)時(shí)間為ikt ，在可容忍時(shí)間之后且在最晚工作時(shí)間內(nèi)（(LLTi, T2]區(qū)間），則時(shí)間窗服務(wù)價(jià)值函數(shù)U (tik)如下：

證明：

由假設(shè)12 和假設(shè)13 可知，客戶(hù)將送貨車(chē)輛到達(dá)時(shí)與參考點(diǎn)（可容忍時(shí)間端點(diǎn)）進(jìn)行比較，若車(chē)輛到達(dá)時(shí)間在可容忍時(shí)間之后，客戶(hù)認(rèn)為不符合自身時(shí)間要求，即時(shí)間窗服務(wù)價(jià)值為負(fù)，可看作損失，應(yīng)

選用價(jià)值函數(shù)中xi＜ 0部分，即V ( xi)=-λ · (-xi)β。

當(dāng)tik∈(LLTi, T2,]，時(shí)間窗服務(wù)價(jià)值函數(shù) U (tik)為：

由假設(shè)12 和假設(shè)13 可知，送貨車(chē)輛到達(dá)時(shí)間為L(zhǎng)LTi時(shí)，時(shí)間窗服務(wù)價(jià)值數(shù)值可看作0，而 T2時(shí)，時(shí)間窗服務(wù)價(jià)值數(shù)值為-1，即：

由式（14）得：

將 a2代入式（15），得：

即：

式中λ 可約去，得：

式（12）得證。

推論在假設(shè)7 至假設(shè)13 背景下的越庫(kù)配送網(wǎng)絡(luò)中，如果送貨車(chē)輛k 在 tik時(shí)間到達(dá)客戶(hù)點(diǎn)i，則客戶(hù)時(shí)間窗服務(wù)價(jià)值函數(shù)模型為：

其圖形見(jiàn)圖2，其中，T1、T2分別為客戶(hù)最早、最晚工作時(shí)間，EETi是客戶(hù)期待服務(wù)時(shí)間節(jié)點(diǎn)、LLTi是客戶(hù)可容忍時(shí)間節(jié)點(diǎn)。

圖2：時(shí)間窗服務(wù)價(jià)值函數(shù)

車(chē)輛在不同時(shí)間段到達(dá)客戶(hù)點(diǎn)，客戶(hù)對(duì)配送服務(wù)的感知會(huì)有所不同，相應(yīng)的時(shí)間窗服務(wù)價(jià)值含義也有所不同：

（1）車(chē)輛在iEET 節(jié)點(diǎn)之前到達(dá)，客戶(hù)對(duì)配送服務(wù)感知最為滿(mǎn)意，時(shí)間窗服務(wù)價(jià)值取值均為1；

（2）車(chē)輛在(EETi,LLTi]內(nèi)到達(dá)，時(shí)間窗服務(wù)價(jià)值則隨時(shí)間減小，特別的，在LLTi節(jié)點(diǎn)上，仍滿(mǎn)足客戶(hù)配送要求，亦可看作收益，因此取值0.5；

（3）車(chē)輛在iLLT 之后到達(dá)，則看作損失，且在2T 節(jié)點(diǎn)上取得最小值-1。

4.數(shù)學(xué)模型

（1）目標(biāo)函數(shù)

（2）約束條件

其中：式（21）和式（22）為目標(biāo)函數(shù)，式(21)是總成本最小化目標(biāo)，包括車(chē)輛運(yùn)輸成本以及車(chē)輛固定成本；式（22）是平均客戶(hù)時(shí)間窗服務(wù)價(jià)值最大化目標(biāo)；式（23）和式（24）分別是各取貨點(diǎn)和送貨點(diǎn)取貨或送貨的貨物量計(jì)算公式；式（25）是送貨車(chē)輛最大載重約束；式（26）表示整個(gè)配送過(guò)程所用車(chē)輛數(shù)目不超過(guò)可用車(chē)輛數(shù)；式（27）表示每個(gè)取/送貨點(diǎn)僅被服務(wù)一次且僅由一輛車(chē)服務(wù)；式（28）-式（31）表示每個(gè)訂單僅被分配一輛車(chē)，且同一個(gè)取/送貨點(diǎn)的訂單分配同一輛車(chē)；式（32）-式（33）表示變量 jrky 、'irky 與abkx 之間的正確關(guān)系；式（34）保證車(chē)輛路線的連續(xù)性；式（35）表示取貨車(chē)輛協(xié)同到達(dá)越庫(kù)中心；式（36）是取貨車(chē)輛到達(dá)越庫(kù)中心時(shí)間的計(jì)算公式；式（37）表示送貨車(chē)輛離開(kāi)越庫(kù)中心的時(shí)間約束；式（38）和式（39）表示送貨車(chē)輛在配送階段的時(shí)間約束；式（40）和式（41）表示決策變量屬性。

三、求解算法

越庫(kù)的路徑優(yōu)化屬于NP-hard 問(wèn)題，且雙目標(biāo)優(yōu)化模型的求解難度遠(yuǎn)大于單目標(biāo)優(yōu)化，宜采用啟發(fā)式算法求解。禁忌搜索算法（Tabu Search，TS）通過(guò)使用禁忌列表，能夠避免算法在局部最優(yōu)點(diǎn)迂回搜索，實(shí)現(xiàn)全局最優(yōu)化[15]。在此基礎(chǔ)上，本文結(jié)合局部搜索算法（Local Search，LS）設(shè)計(jì)了混合多目標(biāo)求解算法。整個(gè)算法主要分為兩層：外層使用TS 算法來(lái)優(yōu)化越庫(kù)車(chē)輛路徑問(wèn)題，形成配送方案，其中在內(nèi)部嵌套了一層LS 算法判斷送貨車(chē)的載量，嵌套程序?qū)ζ浞桨高M(jìn)行判斷，將結(jié)果返回主程序，主程序再根據(jù)其返回值對(duì)路徑方案進(jìn)行調(diào)整，以此遞推，最終得到較優(yōu)的車(chē)輛路徑方案。

算法具體步驟如下：

Step1：初始化算法參數(shù)。包括最大迭代次數(shù)max_iter，禁忌長(zhǎng)度tb_length，鄰域解個(gè)數(shù)neighbor_coun，禁忌表tabulist=?，全局非劣解集局部非劣解集

Step2：生成初始解。采用隨機(jī)方式產(chǎn)生初始送貨路徑序列Initial_TSP，由此得到初始解Initial_Sol，并置當(dāng)前解Current_Sol=Initial_Sol，

Step3：計(jì)算目標(biāo)函數(shù)值 f1、 f2，并判斷是否滿(mǎn)足時(shí)間窗約束條件，若不滿(mǎn)足，則令 f2= f2- M，M是一個(gè)很大的數(shù)。

Step4：調(diào)用局部搜索算法，判斷裝車(chē)是否成功，若成功直接返回 f1、 f2，否則返回 f1= f1+M，f2= f2- M。

Step5：對(duì)當(dāng)前解進(jìn)行鄰域操作，得到候選解集Candidate_Sol。

Step6：重復(fù)Step3 和Step4，得候選解Candidate_Soli 的目標(biāo)函數(shù)值 f1(i )、f2(i )。

Step7：更新局部非劣解、全局非劣解，在候選解集中，應(yīng)用非支配排序算法篩選局部非劣解集與全局非劣解

Step 8：判斷候選解集是否滿(mǎn)足藐視準(zhǔn)則。若是，則將對(duì)應(yīng)的禁忌對(duì)象替換最早進(jìn)入禁忌表的禁忌對(duì)象，更新全局非劣解集；否則，在非禁忌對(duì)象的候選非劣解中選擇一個(gè)作為當(dāng)前解，并將相應(yīng)的禁忌對(duì)象替換最早進(jìn)入禁忌表的禁忌對(duì)象，更新全局非劣解集。

Step9：終止條件判斷。若達(dá)到最大循環(huán)迭代次數(shù)，輸出全局非劣解集結(jié)束算法；否則繼續(xù)Step5-Step8。

四、算例與仿真

（一）算例與實(shí)驗(yàn)設(shè)計(jì)

越庫(kù)車(chē)輛路徑優(yōu)化模型尚無(wú)標(biāo)準(zhǔn)測(cè)試算例，本文選用已有算例數(shù)據(jù)，由10、50 個(gè)節(jié)點(diǎn)組成的兩種規(guī)模問(wèn)題，并根據(jù)研究?jī)?nèi)容對(duì)數(shù)據(jù)進(jìn)行擴(kuò)展[5]，算例參數(shù)值見(jiàn)表1。

兩組實(shí)驗(yàn)設(shè)計(jì)如下：

實(shí)驗(yàn)1：對(duì)小規(guī)模算例進(jìn)行仿真求解，研究協(xié)同到庫(kù)和分揀整合策略在越庫(kù)車(chē)輛路徑中的作用。

實(shí)驗(yàn)2：對(duì)大規(guī)模算例進(jìn)行仿真，對(duì)比研究存在負(fù)值的時(shí)間窗服務(wù)價(jià)值（場(chǎng)景1）與其他文獻(xiàn)中服務(wù)水平最低值僅為0（場(chǎng)景2）。

表1：參數(shù)信息

表2：第1 組實(shí)驗(yàn)取貨調(diào)度方案

表3：第1 組實(shí)驗(yàn)Pareto 解集及送貨方案

（二）實(shí)驗(yàn)結(jié)果與分析

本文應(yīng)用MATLAB R2016a 進(jìn)行仿真實(shí)驗(yàn)，系統(tǒng)是Windows 7，硬件環(huán)境是Intel（R）Core（TM）i3-7100 CPU@3.90GHz，內(nèi)存3.87G。經(jīng)過(guò)反復(fù)測(cè)試，參數(shù)設(shè)置如下：迭代次數(shù)max_iter=500，鄰域解個(gè)數(shù)neighbor_count=160，禁忌長(zhǎng)度tb_length=20。根據(jù)Kahneman 的結(jié)論[16]，α 、β 取值為0.88。通過(guò)算法求解，實(shí)驗(yàn)1 仿真結(jié)果見(jiàn)表2 和表3。

實(shí)驗(yàn)1 結(jié)果分析：

1.由表2 可知，取貨階段采用直運(yùn)方式，共需4 輛車(chē)，路徑1 中，行車(chē)時(shí)間為186.3 個(gè)單位時(shí)間，在供應(yīng)商處裝貨時(shí)間為19 個(gè)單位時(shí)間，路徑總時(shí)間為205.3 個(gè)單位時(shí)間。4 輛車(chē)的路徑總時(shí)間在[169.9,225.5]區(qū)間，根據(jù)協(xié)同到庫(kù)策略要求，以225.5 作為取貨車(chē)輛協(xié)同到庫(kù)基準(zhǔn)，據(jù)此調(diào)度取貨車(chē)輛的發(fā)車(chē)順序，為后續(xù)貨物進(jìn)行分揀整合提供有效保障。

2.由表3 可知，24 份訂單貨物在越庫(kù)中心根據(jù)不同目的地被分揀重裝為6 份，并由3 輛送貨車(chē)進(jìn)行后續(xù)配送。

表4：第2 組實(shí)驗(yàn)仿真結(jié)果

實(shí)驗(yàn)2 仿真結(jié)果如表4 所示，場(chǎng)景1 所得解集中，總成本最低為7305.6307，時(shí)間窗服務(wù)價(jià)值最高為0.8234，當(dāng)時(shí)間窗服務(wù)價(jià)值上升時(shí)，總成本會(huì)亦有所增加。時(shí)間窗服務(wù)價(jià)值的提高往往是以增加配送成本為代價(jià)。

對(duì)比場(chǎng)景1 和場(chǎng)景2 可知，通過(guò)改變超過(guò)可容忍時(shí)間的價(jià)值函數(shù)，函數(shù)值取零和取負(fù)的結(jié)果是不同的。在相近成本下，場(chǎng)景2 所衡量的服務(wù)水平會(huì)高于場(chǎng)景1。為考慮具體影響，選取場(chǎng)景1 總成本最低的實(shí)驗(yàn)結(jié)果進(jìn)行下一步分析，根據(jù)相應(yīng)的配送方案計(jì)算每一個(gè)客戶(hù)在場(chǎng)景1 和場(chǎng)景2 下的時(shí)間窗服務(wù)價(jià)值，結(jié)果對(duì)比見(jiàn)圖3。其中，場(chǎng)景1 中時(shí)間窗服務(wù)價(jià)值平均值為0.7277，場(chǎng)景2 平均值為0.7446。

圖3：時(shí)間窗服務(wù)價(jià)值對(duì)比分析圖

實(shí)驗(yàn)2 結(jié)果分析：

（1）時(shí)間窗服務(wù)價(jià)值和配送總成本之間存在相關(guān)性，兩者互相沖突，決策者可根據(jù)實(shí)際情況選擇合適的配送方案。

（2）由圖3 可知，該方案有6 個(gè)客戶(hù)點(diǎn)出現(xiàn)了配送延遲，超過(guò)了客戶(hù)可容忍時(shí)間，兩種方法計(jì)算出的客戶(hù)時(shí)間窗服務(wù)價(jià)值平均值分別為0.7446 和0.7277，由于存在負(fù)數(shù)，降低了整個(gè)配送的服務(wù)水平。

（3）場(chǎng)景1 中的客戶(hù)時(shí)間窗服務(wù)價(jià)值存在負(fù)值，直觀地表示出客戶(hù)負(fù)向的心理認(rèn)知，大大影響了其總體配送服務(wù)水平，應(yīng)引起決策者的重視。

五、結(jié)論

在越庫(kù)配送系統(tǒng)中，通過(guò)協(xié)同到庫(kù)策略，相同目的地的貨物得以在越庫(kù)中心得到分揀整合，合理規(guī)劃車(chē)輛路徑，有助于進(jìn)一步降低配送成本，同時(shí)提高配送服務(wù)水平。本文同時(shí)考慮了協(xié)同到庫(kù)、分揀和客戶(hù)時(shí)間窗服務(wù)，以配送總成本最小、客戶(hù)時(shí)間窗服務(wù)價(jià)值最大化為優(yōu)化目標(biāo)，建立越庫(kù)配送路徑多目標(biāo)優(yōu)化模型；針對(duì)模型特點(diǎn)，設(shè)計(jì)了禁忌搜索算法和局部搜索算法相結(jié)合的混合算法尋求最優(yōu)配送路徑。得出以下結(jié)論：

（1）越庫(kù)中心的協(xié)同到庫(kù)及分揀策略、配送車(chē)輛運(yùn)輸路徑是越庫(kù)問(wèn)題中較重要的影響因素，將兩階段集成研究能有效提高越庫(kù)的整體配送效率，更符合實(shí)際。

（2）采用客戶(hù)時(shí)間窗服務(wù)價(jià)值，符合客戶(hù)的心理變化，反映出客戶(hù)對(duì)越庫(kù)配送服務(wù)的感知水平。

（3）禁忌搜索算法和局部搜索算法相結(jié)合的混合多目標(biāo)求解算法，能降低無(wú)效路徑生成的概率，提高算法求解效率。

在現(xiàn)實(shí)中，越庫(kù)配送網(wǎng)絡(luò)的應(yīng)用更為復(fù)雜，未來(lái)將進(jìn)一步研究多個(gè)越庫(kù)中心、多個(gè)周期的動(dòng)態(tài)越庫(kù)系統(tǒng)調(diào)度。