不確定條件下配送回收中心選址配送問題研究

康 凱，王小宇，馬艷芳

河北工業(yè)大學(xué) 經(jīng)濟管理學(xué)院，天津 300401

1 引言

資源緊缺和環(huán)境惡化已威脅到人類的生存和發(fā)展，我國做出了“發(fā)展循環(huán)經(jīng)濟、建設(shè)節(jié)約型社會”的重大戰(zhàn)略決策，再制造不僅節(jié)約資源還有助于保護環(huán)境，受到政府的重視和社會各界的關(guān)注，閉環(huán)供應(yīng)鏈將正向物流與逆向物流相結(jié)合，解決產(chǎn)品回收再制造的問題。在閉環(huán)供應(yīng)鏈的管理中，選址配送是閉環(huán)供應(yīng)鏈網(wǎng)絡(luò)設(shè)計與運營的重要問題。

不確定閉環(huán)供應(yīng)鏈問題引起了學(xué)者的廣泛關(guān)注，在閉環(huán)供應(yīng)鏈網(wǎng)絡(luò)中存在很多不確定因素，逆向物流中的不確定因素要多于正向物流。現(xiàn)階段對閉環(huán)供應(yīng)鏈中不確定因素的考慮集中在正向物流中的需求不確定及逆向物流中的回收數(shù)量不確定。Keyvanshokooh等[1]用模糊隨機表示需求和回收量不確定，建立混合整數(shù)線性規(guī)劃（MILP）模型解決閉環(huán)供應(yīng)鏈網(wǎng)絡(luò)設(shè)計問題。Khatami等[2]運用基于場景的隨機MILP模型研究了市場需求和產(chǎn)品回收數(shù)量不確定的多階段多產(chǎn)品閉環(huán)供應(yīng)鏈網(wǎng)絡(luò)設(shè)計問題，并運用Benders和Cholesky’s分解求解。范小三等[3]探究了需求和回收隨機條件下分散決策與集中決策中產(chǎn)品定價策略，在產(chǎn)品回收波動性增大時，第三方物流可以提高回收價格轉(zhuǎn)嫁風(fēng)險。仍有部分學(xué)者關(guān)注到回收產(chǎn)品質(zhì)量的不確定，但關(guān)于回收質(zhì)量的不確定多采用對回收產(chǎn)品分類的方法。高陽等[4]將廢舊產(chǎn)品分為可維修、可拆解、可分解三類，研究了回收渠道選擇問題，得出制造商的風(fēng)險偏好對回收渠道選擇有重要影響的結(jié)論。劉枚蓮等[5]通過設(shè)定回收產(chǎn)品質(zhì)量類別，研究了回收質(zhì)量不確定的逆向供應(yīng)鏈定價問題。本文采用回收產(chǎn)品可再利用率表示回收產(chǎn)品質(zhì)量的不確定，更具有現(xiàn)實意義。

由于環(huán)境污染嚴(yán)重，社會的關(guān)注及政府的政策法規(guī)要求，企業(yè)開始關(guān)注運營過程中的環(huán)境問題。注重環(huán)保要求并采取措施改善，有助于提高企業(yè)的社會形象。同時有很多文獻(xiàn)在閉環(huán)供應(yīng)鏈優(yōu)化中考慮減少碳排放和成本優(yōu)化雙重目標(biāo)。高舉紅等[6]研究了市場不確定有碳補貼政策的閉環(huán)供應(yīng)鏈網(wǎng)絡(luò)優(yōu)化設(shè)計問題。戴卓等[7]人以碳排放與成本最低建立多目標(biāo)優(yōu)化模型解決需求不確定閉環(huán)供應(yīng)鏈網(wǎng)絡(luò)優(yōu)化設(shè)計問題。Garg等[8]提出了交互多目標(biāo)優(yōu)化算法解決閉環(huán)供應(yīng)鏈網(wǎng)絡(luò)設(shè)計中的環(huán)境問題，確定了最優(yōu)設(shè)施間流量以及最優(yōu)運輸車輛數(shù)量。Talaei等[9]以碳排放和供應(yīng)鏈總成本最低為目標(biāo)，運用ε-約束法解決可變成本和需求率不確定條件下多產(chǎn)品閉環(huán)供應(yīng)鏈中的選址配送問題。Jindal[10]運用模糊隨機理論處理不確定變量，考慮經(jīng)濟和環(huán)境因素建立了多目標(biāo)閉環(huán)供應(yīng)鏈網(wǎng)絡(luò)設(shè)計模型。

閉環(huán)供應(yīng)鏈選址與配送問題是NP難問題。模糊隨機環(huán)境下閉環(huán)供應(yīng)鏈中配送回收中心選址配送問題更為復(fù)雜。粒子群算法由于實現(xiàn)較為簡單，自適應(yīng)能力強的特點已成為一種重要的智能算法并廣泛應(yīng)用于各個學(xué)科及領(lǐng)域。粒子群算法在解決NP難問題中有明顯優(yōu)勢，改進(jìn)粒子群算法在解決供應(yīng)鏈管理問題應(yīng)用更加廣泛。Ai和Kachitvichyanukul[11]提出的全局-局部-鄰域粒子群算法，Xu等[12]提出的雙層模糊隨機仿真全局-局部-鄰域粒子群算法，為本文模型求解提供了方法基礎(chǔ)。本文提出基于優(yōu)先級的全局-局部-鄰域粒子群算法設(shè)計求解方案并求解。

2 問題描述和模型建立

2.1 問題描述

本文以含制造商、零售商、一個廢棄處理中心和多個配送回收中心構(gòu)成的閉環(huán)供應(yīng)鏈為研究對象，如圖1所示。產(chǎn)品在工廠完成生產(chǎn)，運往配送回收中心進(jìn)行調(diào)度配送，零售商在接收新產(chǎn)品的同時把回收的廢舊產(chǎn)品交給運輸車輛運回配送回收中心，在配送回收中心檢測、分類，可再利用的運回工廠進(jìn)行再制造，不可再利用的運往廢棄處理中心按國家環(huán)保要求進(jìn)行廢棄處理。配送回收中心選址的目標(biāo)是考慮現(xiàn)有約束情況下，確定配送回收中心的數(shù)量和位置，為后期網(wǎng)絡(luò)運營奠定良好的基礎(chǔ)，節(jié)約成本。

圖1 閉環(huán)供應(yīng)鏈網(wǎng)絡(luò)結(jié)構(gòu)示意圖

制造商采用按訂單生產(chǎn)模式生產(chǎn)。由于顧客可能會將產(chǎn)品丟棄或作為他用，廢舊產(chǎn)品回收率考慮為模糊隨機變量；在運輸以及搬運過程中可能會導(dǎo)致回收產(chǎn)品受損甚至破碎，或者已回收產(chǎn)品中存在有損產(chǎn)品不可再修復(fù)，產(chǎn)品可再利用率考慮為模糊隨機變量。運輸成本和運輸碳排放與運輸距離以及運輸產(chǎn)品數(shù)量有關(guān)。本文提出如下假設(shè)：（1）只考慮單一產(chǎn)品單一階段[13]，從制造商開始的正向物流作為起點，再到產(chǎn)品回到制造商進(jìn)行再制造結(jié)束；（2）配送回收中心的備選點已知，工廠、零售商和廢棄處理中心的位置已知；（3）配送回收中心和制造再制造工廠有能力上限[14]，考慮有能力上限的設(shè)施規(guī)劃更加符合實際情況；（4）回收產(chǎn)品再制造的成本低于購買新的原材料進(jìn)行再生產(chǎn)的成本[15]；從經(jīng)濟利益的角度出發(fā)使得企業(yè)愿意采取回收行動，推動閉環(huán)供應(yīng)鏈的建立與運營；（5）企業(yè)采取需求拉動的生產(chǎn)模式，零售商需求根據(jù)訂單已知[16]；（6）交通運輸導(dǎo)致的碳排放量在全球碳排放中占有重要比例，是供應(yīng)鏈減排中最為關(guān)鍵的環(huán)節(jié)，故僅考慮運輸過程中的碳排放量，且與運輸距離以及運輸量有關(guān)[17]。

2.2 符號說明

Ω：配送回收中心備選點集，Ω={1,2,…,I}；

Ψ ：工廠地點集，Ψ={1,2,…,J}；

Φ：零售商集合，Φ={1,2,…,K}；

Y：廢棄處理中心集合，Y={1,2,…,N}；

i：配送回收中心備選點，i∈Ω={1,2,…,I}；

j：已知工廠位置，j∈Ψ={1,2,…,J}；

k：已知零售商位置，k∈Φ={1,2,…,K}；

n：已知廢氣處理中心位置，n∈Y={1,2,…,N}；

U：配送回收中心數(shù)量上限；

Dk：零售商需求；

αi：配送回收中心i的處理能力上限；

γj：工廠 j的生產(chǎn)能力；

Pji：從工廠 j運往配送回收中心i的產(chǎn)品數(shù)量；

Qik：從配送回收中心i運往零售商k的產(chǎn)品數(shù)量；

C0：碳排放目標(biāo)系數(shù)；

βij：配送回收中心i到工廠 j之間的單位產(chǎn)品碳排放量；

βik：配送回收中心i到零售商k之間的單位產(chǎn)品碳排放量；

βin：配送回收中心i到廢棄處理中心n之間的單位產(chǎn)品碳排放量；

xi：為0-1變量，備選點i被選中建立配送回收中心為1，否則為0；

yik：為0-1變量，零售商k由配送回收中心i服務(wù)為1，否則為0。

2.3 模型構(gòu)建

2.3.1 目標(biāo)函數(shù)

本文以建設(shè)配送回收中心的經(jīng)濟成本和產(chǎn)品配送回收運輸過程碳排放最低為目標(biāo)建立如下模型。

經(jīng)濟目標(biāo)：

碳排放目標(biāo)：

2.3.2 約束條件

約束條件主要包括流量均衡約束，配送回收中心數(shù)量約束，處理能力上限約束。

（1）由工廠運往配送回收中心的新產(chǎn)品與由零售商處回收的廢舊產(chǎn)品之和不超過配送回收中心的處理能力上限。

（2）新產(chǎn)品的需求量和廢舊產(chǎn)品的可再利用量不超過工廠的生產(chǎn)能力上限。

（3）產(chǎn)品回收總量不超過產(chǎn)品生產(chǎn)總量。

（4）零售商k的需求被滿足。

（5）可再利用產(chǎn)品數(shù)量不多于回收產(chǎn)品數(shù)量。

（6）配送回收中心的數(shù)量不少于1個。

（7）配送回收中心的數(shù)量不超過設(shè)定的上限。

（8）每個零售商都有且只有一個配送回收中心向其提供服務(wù)。

（9）xi是0-1變量，備選點i被選中建立配送回收中心為1，否則為0。

（10）yik是0-1變量，零售商k由配送回收中心i服務(wù)為1，否則為0。

模型中所涉及的配送回收中心經(jīng)濟成本最小和碳排放最小，由于其屬于不同范疇，量綱不同，故進(jìn)行加權(quán)平均求解[20]。本文采用改進(jìn)加權(quán)平均求解法，給碳排放目標(biāo)賦予系數(shù)C0，將其理解為碳排放產(chǎn)生的經(jīng)濟成本，構(gòu)成目標(biāo)函數(shù)min F=Z1+C0Z2，將多目標(biāo)規(guī)劃轉(zhuǎn)化為單目標(biāo)規(guī)劃問題。

3 算法改進(jìn)方案

粒子群算法是對鳥群覓食行為模擬仿真的智能算法，由Kennedy和Eberhart[21]在1995年首次提出。粒子群算法經(jīng)常被用來解決NP-難問題。近些年的研究發(fā)現(xiàn)標(biāo)準(zhǔn)粒子群算法存在一定缺陷。在標(biāo)準(zhǔn)粒子群算法當(dāng)中，每個粒子在每次迭代過程中都要進(jìn)行兩個方面的學(xué)習(xí)，個體最優(yōu)以及全局最優(yōu)，并通過學(xué)習(xí)到的數(shù)據(jù)調(diào)整自己即將飛行的方向和速度。但是一個種群中的粒子易向全局最優(yōu)粒子附近聚集，導(dǎo)致該種群頻繁陷入局部最優(yōu)解且不再更新。為了處理這種過早收斂于局部最優(yōu)解的缺陷，Ai和Kachitvichyanukul[11]提出了全局-局部-鄰域粒子群算法，有效地解決了這一問題。本文在此基礎(chǔ)上，提出基于優(yōu)先級的全局-局部-鄰域粒子群算法來解決多目標(biāo)配送回收中心的選址配送問題。

3.1 粒子群算法中的符號

基于優(yōu)先級的全局局部鄰域粒子群算法中的符號如下所示：

τ：迭代代數(shù)，τ=1,2,…,T ；

d：維度，d=1,2,…,D；

l：粒子，l=1,2,…,L；

ωτ：第τ代粒子的慣性權(quán)重；

vld(τ)：第l個粒子在第τ代在d維度上的速度分量；

cp：個體加速度；

cl：局部加速度；

cg：全局加速度；

cn：鄰域加速度；

Pmax：粒子的最大位置分量；

Pmin：粒子的最小位置分量；

Vl：粒子的速度向量；

Pl：粒子的位置向量；

Gbest(τ)：第τ代粒子群的歷史最優(yōu)位置向量；

r1,r2,r3,r4：[0，1]上的相互獨立的隨機數(shù)；

Fitness(Pl)：第l個粒子所在位置對應(yīng)的適應(yīng)值函數(shù)值。

3.2 初始化及適應(yīng)值計算

根據(jù)Gen和Altiparmak[22]提出的優(yōu)先級編碼解碼方法進(jìn)行編碼設(shè)計。分兩階段設(shè)置優(yōu)先級，第一階段是配送回收中心與零售商共同設(shè)置優(yōu)先級，優(yōu)先級為粒子，對應(yīng)配送回收中心和零售商，優(yōu)先級最高的配送回收中心選擇距離最近的零售商進(jìn)行配送（優(yōu)先級最高的零售商選擇距離最近的配送回收中心服務(wù)），具體過程如下：

步驟1隨機生成I+K個數(shù)字表示對應(yīng)零售商與配送回收中心備選點。

步驟2判斷最高優(yōu)先級對應(yīng)的是配送回收中心備選點或零售商的位置。

步驟3選擇距離最近得零售商（或備選點）進(jìn)行配送，如果是配送回收中心選擇距離最近的零售商進(jìn)行配送后，仍有剩余則選擇下一個距離近的進(jìn)行配送，直到不夠配送新的零售商；完成配送后零售商需求為0，對應(yīng)位置的距離變?yōu)闊o窮大。

步驟4重復(fù)前兩步，直到所有零售商的需求被滿足。

第二階段過程與第一階段相同，優(yōu)先級設(shè)定是選中的備選點和工廠。

3.3 約束處理

根據(jù)公式計算適應(yīng)值之后，對于非可行解引入懲罰因子，使其適應(yīng)度值表現(xiàn)較可行解相差很多，這樣可以減少在后續(xù)優(yōu)化過程中非可行解的出現(xiàn)概率[23]。針對本文中的問題，針對流量均衡約束以及處理能力約束采用拒絕法將其取值限制在給定數(shù)值范圍內(nèi)[24]，針對配送回收中心數(shù)量約束設(shè)立懲罰因子μ，當(dāng)選擇的配送回收中心備選點多余上限值時，則在適應(yīng)值計算結(jié)果上加一個極大值作為懲罰因子，將非可行解淘汰。如果配送回收中心數(shù)量滿足要求，則適應(yīng)值計算函數(shù)為：

若配送回收中心數(shù)量不滿足約束條件，則適應(yīng)值計算函數(shù)為：

3.4 更新公式

根據(jù)Ai和Kachitvichyanukul[11]提出的全局局部鄰域粒子群算法中的符號，慣性權(quán)重和粒子速度以及位置根據(jù)以下公式進(jìn)行更新。

3.5 算法流程

本文中，上述全局-局部-鄰域粒子群算法被用于解決設(shè)施選址及配送問題。由于不確定性的存在和問題的獨特性，提出基于優(yōu)先級的全局-局部-鄰域粒子群算法來解決本文模型。如圖2所示。

圖2 算法流程圖

步驟1初始化粒子（即優(yōu)先級）。

步驟2約束檢查，是否滿足約束要求，滿足則進(jìn)入步驟3，不滿足，返回步驟1。

步驟3計算目標(biāo)函數(shù)值。

步驟4根據(jù)更新公式更新粒子速度及位置。

計算每個粒子的適應(yīng)值，并記錄個人最優(yōu)位置，并在所有最優(yōu)位置中選擇最好的最為全局最優(yōu)位置；更新粒子的個體最優(yōu)位置，種群最優(yōu)位置，局部最優(yōu)位置，生成鄰域最優(yōu)位置，更新每個粒子的速度與位置，檢查粒子是否超出邊界。

步驟5檢查是否滿足終止條件，若滿足，則進(jìn)入下一步，若不滿足返回步驟4。

步驟6根據(jù)得出的優(yōu)先級解碼得出選中的備選點。

4 案例分析

4.1 數(shù)據(jù)收集

本文以D制造商的某種產(chǎn)品為例，根據(jù)T地區(qū)的相關(guān)情況構(gòu)建一個閉環(huán)供應(yīng)鏈網(wǎng)絡(luò)。顧客將結(jié)束使用的產(chǎn)品送回零售商，然后在配送回收中心檢測、分類、簡單修護，可進(jìn)行再制造的送往制造商，不可再利用的送往配送回收中心進(jìn)行廢棄處理。本閉環(huán)供應(yīng)鏈網(wǎng)絡(luò)中有3個制造再制造工廠，6個備選點作為建設(shè)配送回收中心的可選項，為10個零售商服務(wù)，一個廢棄處理中心負(fù)責(zé)處理廢棄物。根據(jù)《2006年IPCC國家溫室氣體清單指南》提供的方法，二氧化碳排放量的計算公式為：單位化石燃料二氧化碳的估算值=低位發(fā)熱量×碳排放因子×碳氧化率×碳轉(zhuǎn)換系數(shù)，建議開征碳稅成本為每噸20元[25-26]，故本文碳排放目標(biāo)函數(shù)系數(shù)C0取值為20。配送回收中心備選點的相關(guān)信息如表1所示；各零售商的坐標(biāo)及需求見表2；工廠和廢棄處理中心的坐標(biāo)以及能力上限如表3所示。本文中產(chǎn)品回收率及可再利用率為模糊隨機變量，兩個變量用三角模糊數(shù)表示，如表4、5所示。

表1 配送回收中心備選點相關(guān)參數(shù)

表2 零售商相關(guān)參數(shù)

表3 工廠及廢棄處理中心相關(guān)參數(shù)

4.2 算法參數(shù)實驗

對模型進(jìn)行求解前，首先要設(shè)置算法參數(shù)值。對算法運行效果的驗證，探索不同的種群大小及迭代次數(shù)對計算結(jié)果和運行時間的影響。在不同的種群大小和迭代次數(shù)組合下，分別運行10次基于優(yōu)先級的全局-局部-鄰域粒子群算法。表6給出了不同種群大小和迭代次數(shù)組合下算法的最優(yōu)解、最差解、標(biāo)準(zhǔn)差和平均計算時間。從表6中可以看出，種群大小為10時，在給出的迭代次數(shù)范圍內(nèi)，雖然計算時間短，但是結(jié)果較差且不穩(wěn)定。由表可知，種群大小一定，迭代次數(shù)從100增加到300，計算時間變長，但結(jié)果不一定更優(yōu)。在粒子數(shù)為20，迭代次數(shù)為200代時取得最優(yōu)值。當(dāng)粒子數(shù)增加到30時，隨著迭代次數(shù)的增加，計算時間增加，但計算結(jié)果及標(biāo)準(zhǔn)差無明顯規(guī)律，且沒有更優(yōu)于20個粒子，200代的運行結(jié)果。通過對比，發(fā)現(xiàn)當(dāng)?shù)螖?shù)一定時，粒子數(shù)增加，計算時間增加，但是計算結(jié)果不一定更優(yōu)。

表4 回收率

4.3 結(jié)果分析

基于上述數(shù)據(jù)，通過MATLAB編程，基于優(yōu)先級的全局-局部-鄰域粒子群算法有效求解了模型。對于算法中的參數(shù)設(shè)置如下：粒子種群大?。篘=20，最大迭代次數(shù)：T=200，慣性權(quán)重：ω(1)=1，ω(T)=0.1，加速常數(shù)cp=cg=cl=cn=2。將程序運行10次后，獲得最優(yōu)目標(biāo)函數(shù)值，即最低成本為2 975.71萬元，為與之對應(yīng)的計算結(jié)果是選擇備選點2，4，5建設(shè)配送回收中心。由D2負(fù)責(zé)配送回收R1，R2，R5；D4負(fù)責(zé)配送回收R3，R4，R6；D5負(fù)責(zé)配送回收R7，R8，R9，R10。如圖3所示。

圖3 最優(yōu)選址配送結(jié)果

4.4 算法對比

為驗證算法的有效性和實用性，對基于優(yōu)先級的全局-局部-鄰域粒子群算法，遺傳算法，蟻群算法以及免疫優(yōu)化算法之間進(jìn)行對比。同樣使用MATLAB對算法進(jìn)行編譯，每種算法在粒子數(shù)為20，迭代次數(shù)為200的條件下運行10次。在全局-局部-鄰域粒子群算法中，慣性權(quán)重ω(1)=1，ω(T)=0.1，加速度cp=cg=cl=cn=2；在遺傳算法中交叉概率為1，變異概率為0.1；在蟻群算法中信息素重要程度因子為1，啟發(fā)函數(shù)重要程度因子為5，信息素?fù)]發(fā)因子為0.1，常系數(shù)1；在免疫優(yōu)化算法中交叉概率為1，變異概率為0.1。

表5 零售商到配送回收中心的產(chǎn)品可再利用率

表6 計算結(jié)果

由圖4可知，四種算法都有同樣的移動趨勢，隨著迭代次數(shù)的增加，算法結(jié)果越來越好。四種算法在迭代初期結(jié)果都比較差，可能是探索空間超出可行域并造成了對函數(shù)適應(yīng)值的懲罰。隨著程序繼續(xù)運行和迭代次數(shù)的增加，函數(shù)適應(yīng)值減小，在迭代結(jié)束時，算法結(jié)果趨于穩(wěn)定且變得更好。

圖4 算法對比圖

基于優(yōu)先級的全局-局部-鄰域粒子群算法優(yōu)于遺傳算法、蟻群算法以及免疫優(yōu)化算法?；趦?yōu)先級的全局-局部-鄰域粒子群算法相較其他算法結(jié)果更優(yōu)收斂更快?；趦?yōu)先級的全局-局部-鄰域粒子群算法和全局-局部-鄰域粒子群算法在70代左右開始收斂，遺傳算法與蟻群算法在80代左右開始收斂，免疫算法則在100代后趨于穩(wěn)定。將四種算法的計算結(jié)果放在表7中通過對比發(fā)現(xiàn)，基于優(yōu)先級的全局-局部-鄰域粒子群算法相較于其他算法的優(yōu)越性。

表7 算法對比結(jié)果

5 結(jié)論

本文以成本最低和碳排放最小為目標(biāo)，考慮能力約束，流量約束以及顧客需求約束，運用模糊隨機理論，構(gòu)建了閉環(huán)供應(yīng)鏈配送回收中心選址配送模型；提出了基于優(yōu)先級的全局-局部-鄰域粒子群算法，并通過案例驗證了算法的先進(jìn)性以及模型的有效性。本文提出的結(jié)合分配中心和回收中心功能的配送回收中心，不僅可以降低企業(yè)建設(shè)成本還可降低運輸車輛空載運輸成本，同時降低運輸碳排放。同時考慮產(chǎn)品回收率及可再利用率的模糊隨機性，建立的以成本最低和碳排放最小為目標(biāo)的閉環(huán)供應(yīng)鏈選址配送模型，為模糊隨機環(huán)境條件下考慮碳排放的閉環(huán)供應(yīng)鏈選址配送決策提供理論依據(jù)與方法。多階段、多產(chǎn)品考慮零售商需求不確定性的閉環(huán)供應(yīng)鏈選址配送是進(jìn)一步研究的方向。