考慮貨損和碳排放的生鮮產品配送路徑優化

沈麗　李成玉　甘彥　趙剛

摘要：

為減少生鮮產品運輸過程中的價值損耗和碳排放量，考慮生鮮產品易腐性和冷鏈物流高排放的特點，對貨損和碳排放來源進行細化分析，發現貨損主要包括物理損傷、呼吸作用引起的腐敗損失和冷藏車開門導致生鮮產品加速腐敗的損失，碳排放主要由車輛行駛和制冷機組在運輸和卸貨過程中運行所消耗的燃油產生。以固定成本、貨損成本、燃油成本、碳排放成本和時間懲罰成本之和最小為目標，構建生鮮產品配送路徑優化模型，并用遺傳算法進行求解。算例分析表明，當總成本在可接受范圍內變動時，物流企業可以為客戶提供更優質的服務。

關鍵詞：

冷鏈物流; 車輛路徑問題（VRP）; 貨損; 碳排放; 遺傳算法

中圖分類號：? F252.14; U116.2

文獻標志碼：? A

Distribution route optimization of fresh products considering

cargo damage and carbon emission

SHEN Li， LI Chengyu， GAN Yan， ZHAO Gang

（School of Logistics， Beijing Wuzi University， Beijing 101149， China）

Abstract：

In order to reduce the value loss and carbon emission during the transportation of fresh products，considering the perishable nature of fresh products and the high emission characteristic of cold chain logistics， the sources of cargo damage and carbon emission are analyzed in detail. The cargo damage mainly includes the physical damage， the corruption loss caused by respiration and the accelerated corruption loss of fresh products? caused by the opening of refrigerated trucks. Carbon emission is mainly generated by fuel consumption of vehicle driving and refrigeration unit operation during transportation and unloading. The lowest total cost is taken as the objective， where the total cost includes the fixed cost， the cargo damage cost， the fuel cost， the carbon emission cost and the time penalty cost. The optimization model of fresh product distribution route is constructed. The genetic algorithm is used to solve the model. Example analysis shows that logistics enterprises can provide customers? with better service under the conditions that the total cost changes within acceptable range.

Key words：

cold chain logistics; vehicle routing problem （VRP）; cargo damage; carbon emission;genetic algorithm

收稿日期： 2020-04-09

修回日期： 2020-06-15

基金項目：

北京市自然科學基金（2182027）

作者簡介：

沈麗（1978—），女，天津人，副教授，博士，研究方向為物流與供應鏈管理，（E-mail）lishen97@163.com

0 引 言

生鮮產品易腐敗的特性以及人們對生鮮產品品質要求的提升，對生鮮產品的物流配送問題提出了挑戰，冷鏈配送應運而生。

合理有效規劃生鮮產品的冷鏈物流配送路徑，不僅能提高配送時效性，減少生鮮產品的貨損，還能縮短配送里程，節約配送成本，減少能耗和碳排放，在最大限度保證生鮮產品新鮮度的前提下，提高配送服務水平。

國內外學者對生鮮產品配送車輛路徑問題（vehicle routing problem， VRP）進行了深入研究，已經取得了諸多成果。HSU等[1]將能源消耗成本和違反客戶時間窗的懲罰成本作為優化目標，建立了帶時間窗、隨機的易腐品配送路徑優化模型;OSVALD等[2]提出將易腐性作為關鍵考慮因素對新鮮蔬菜進行配送，建立了基于客戶時間窗和配送時間窗的車輛路徑優化模型;SONG等[3]分析了車型混用的多種易腐品的配送路徑問題，將客戶滿意度與生鮮產品新鮮度結合起來，建立非線性數學模型，在保證易腐品新鮮度的同時，減少配送損耗;孫國華[4]研究帶軟時間窗的開放式滿載車輛的路徑優化問題，認為軟時間窗與現實狀況更加契合，因此構建了懲罰成本隨配送時間變動的線性函數關系;葛顯龍等[5]針對生鮮物流如何將最新鮮的產品以最短時間配送到消費者手中的問題，建立了考慮時間窗和生鮮貨損的配送模型，并設計自適應遺傳算法進行求解;崔巖等[6]以累積前景理論為基礎，將代理點需求、車輛載質量和代理點要求服務的時間窗作為約束條件，以總成本最小為優化目標，建立生鮮電商配送路徑優化模型;姚源果等[7]認為冷藏車制冷能耗成本與冷藏車的運行時間、開門時長和開門次數有關，以總成本最小為優化目標，提出在冷鏈配送中合理設置接駁點可以有效降低配送成本并能提高客戶滿意度。黃星星等[8]、康凱等[9]研究了碳規則下生鮮農產品冷鏈配送路徑優化問題，以減少冷鏈物流的碳排放量。VRP模型求解方法多樣，常用的有大鄰域搜索算法[10]、蟻群算法[11-12]、混合自學習粒子群優化算法[13]、自適應禁忌搜索算法[14]以及遺傳算法[15-16]等。

（2）運輸時制冷機組所產生的燃油消耗。冷藏車產生的熱負荷主要來自通過太陽輻射傳入車廂內部的熱量和車廂漏氣造成的熱量差值，冷藏車在行駛過程中產生的熱負荷可以表示為

式中：δ為常數，由車廂整體劣化程度決定;R是熱導率，單位是1 000cal/（h·m2·℃）（1 cal=4.186 J）;S為車廂受太陽輻射的面積，通常S=SoutSin，其中Sout和Sin分別表示車廂體外和體內表面積，m2;Tout為車廂環境溫度，℃。因此，此項燃油消耗為

式中：δ1為降低單位熱負荷所消耗的燃油量。

（3）冷藏車停車時制冷機組所產生的燃油油耗。

式中：δ2為降低單位熱負荷所消耗的燃油量。由于冷藏車車廂門打開，制冷機組需消耗更多燃油，因此δ2>δ1。

綜上，生鮮產品在整個配送過程中所消耗的燃油成本是冷藏車行駛和制冷機組運行所消耗的總的燃油成本，其表達式為

2.2.4 碳排放成本

碳排放主要指消耗燃油產生的CO2排放量，本文采用碳稅的形式計算碳排放成本，即碳排放成本=碳稅×碳排放量。記λ′為燃油的碳排放系數，則配送過程中產生的碳排放成本C4為

2.2.5 時間懲罰成本

采用軟時間窗約束。在Te，i，Tl，i時間段外將生鮮產品送達，客戶由于某些原因不能接收時，懲罰成本為無窮大M。懲罰成本函數如下：

則總的時間懲罰成本為

2.3 模型建立

綜上所述，考慮貨損和碳排放的生鮮產品配送路徑優化模型為

約束條件除式（18）的時間窗約束外，還有

式（21）表示每個客戶僅有一輛冷藏車提供配送服務;式（22）表示任何車輛經過某客戶點后也必須從該客戶點離開;式（23）表示每輛車從配送中心出發，也必須返回配送中心;式（24）為車輛載貨質量限制;式（25）為決策變量，取值只能是0或1;式（26）和（27）為變量之間的關系。

3 算法設計

車輛路徑優化問題一直是優化組合領域研究的重難點，本文采用相對成熟的遺傳算法求解模型，通過對生成的初始種群進行改善，提高初始種群的適應度，以免陷入局部最優。首先隨機生成滿足實際問題約束的初始解集，再對初始解集反復進行交叉、變異等遺傳操作，將每條染色體上的優秀基因一代代地傳下去，形成具有高適應度的種群，從而獲得最優解，其基本流程見圖2。

具體操作過程如下：

（1）編碼。采用自然數編碼方法。自然數編碼能清楚地將配送車輛、客戶點與自然數列進行對應，能很好地貼合路徑優化問題。染色體按照從左到右的編碼順序，依次將客戶點的需求量和到達客戶點的行駛路程進行累加，當累計需求量和行駛路程大于第1輛冷藏車的最大載質量和最大行駛距離時，記

錄累計次數i，記錄斷點位置i-1，插入0元素;接著從客戶i+1開始，對每個客戶點的需求量和行駛路程進行累加，當累計需求量超過第2輛冷藏車的最大載質量時，記錄此時的累計次數j，記錄斷點位置i+j-1，插入0元素，以此類推。

（2）種群初始化。考慮到遺傳算法對初始種群的質量具有較強的依賴性，對隨機生成的初始種群進行處理：首先對隨機生成的初始種群進行適應度計算，將適應度高的染色體留在種群內，淘汰適應度低的染色體，然后隨機生成部分染色體，直到染色體數量達到種群規模。這樣得到的初始種群具有較高的適應度，提升了種群質量，使算法能較好地找到全局最優解。

（3）交叉操作。采用單點交叉，先確定一個交叉位置，再將雙親中一方染色體的前一部分與另一方染色體的后一部分進行組合，構成一個新的子代染色體。

（4）變異操作。自然數編碼下基因代表客戶點，變異操作就是指該客戶點被替換成其他客戶點，本文采用交換變異，即在染色體上任意選擇兩個基因的位置，然后進行互換操作，互換后產生新的個體。

（5）選擇操作。選用輪盤賭方法進行選擇操作，父代染色體的適應度高低決定該染色體被傳到下一代的機會的多少，適應度越高，被選中的概率越高，這個個體的遺傳信息在種群中擴散的概率就越高。

（6）終止條件。在遺傳代數達到預先設定好的進化代數后，在產生的最后一代種群中，挑選出適應度最高的染色體，將該染色體所對應的配送路徑作為該模型的最優解。

4 算例分析

采用隨機生成數值試驗的方法驗證模型和算法的有效性。假設某城市里有一家冷鏈配送企業，擁有一個大型配送中心（編號為0）和10輛江鈴冷藏車，同時為20個生鮮超市提供配送服務;冷藏車型號相同;冷藏車車廂內的溫度保持在4 ℃;冷藏車自身質量為3 200 kg，最大載質量為1 100 kg;車輛從配送中心出發，完成配送服務后返回配送中心;每輛車的固定成本為300元/輛，平均行駛速度為50 km/h，燃油價格為6.5元/L，車輛提前到達和延遲到達超市的懲罰系數分別為20元/h和15元/h;生鮮產品單位價值為10 000元/t，在裝卸和配送過程中產生的破損率為0.05，生鮮產品變質率與溫度的系數a為0.001 25。各超市的位置坐標、生鮮產品需求量、到達的時間窗約束以及冷藏車服務停留時間均已給出，見表1。

4.1 模型求解

考慮或不考慮貨損和碳排放的模型求解結果見表2。不考慮貨損和碳排放的車輛配送路徑見圖3，考慮貨損和碳排放的車輛配送路徑見圖4;通過對比可知，考慮貨損和碳排放的車輛配送路徑更優。

4.2 數據結果分析

（1）在總成本方面，與不考慮貨損和碳排放的情況相比，考慮貨損和碳排放的配送車輛數量由原來的3輛增加為4輛，總成本提高了12.69%，但是車輛數量的增加使得配送能更好地遵守客戶的時間窗約束，懲罰成本降低了34.80%。綜合來看，在成本可以接受的情形下，考慮貨損和碳排放制定的配送方案能更好地提高客戶滿意度。

（2）在燃油成本方面，考慮貨損和碳排放比不考慮貨損和碳排放的燃油成本降低了29.09%，這是因為考慮貨損和碳排放的配送路徑更優，降低了配送總里程。

（3）從客戶角度分析，配送里程越短，生鮮產品的貨損量就越少，生鮮產品的新鮮度就越高，從而能夠提高客戶滿意度;從配送企業角度分析，考慮貨損

和碳排放后，雖然總成本小幅增加，但能為客戶提供較好的服務體驗，進而提高產品銷量。

（4）貨損成本和燃油成本在總成本中的占比較大，分別占總成本的18.10%、40.53%。配送企業應加強冷鏈配送管理來降低燃油成本和貨損成本。例如：培養司機良好的開車習慣以便減少油耗，裝卸貨時減少冷藏車箱的開門時間和次數以便穩定車廂溫度，規劃良好的行車路徑縮短行駛的總里程等，都可以降低燃油成本和貨損成本。

5 結 論

本文在細化貨損和碳排放來源的基礎上，研究帶時間窗的生鮮產品配送路徑優化問題。綜合考慮固定成本、貨損成本、燃油成本、碳排放成本和時間懲罰成本，構建該問題的基本模型，采用遺傳算法進行求解。通過對比考慮貨損和碳排放與不考慮貨損和碳排放兩種情形下的路徑優化結果，發現考慮貨損和碳排放的總成本略高于不考慮貨損和碳排放的總成本，但燃油成本和時間懲罰成本都所有降低。本文可為生鮮配送企業提供一些借鑒：配送企業可適度增加總成本，在減少冷鏈配送碳排放量的同時，為客戶提供更優質的服務。

參考文獻：

[1]HSU C I，HUNGS S F， LI H C. Vehicle routing problem with time-windows for perishable food delivery[J]. Journal of Food Engineering， 2007， 80（2）： 465-475.

[2]OSVALD A，STIRN L Z. A vehicle routing algorithm for the distribution of fresh vegetables and similar perishable food[J]. Journal of Food Engineering， 2008， 85（2）： 285-295.

[3]SONG B D，KO Y D. A vehicle routing problem of both refrigerated- and general-type vehicles for perishable food products delivery[J]. Journal of Food Engineering， 2016， 169： 61-71.

[4]孫國華. 帶軟時間窗的開放式滿載車輛路徑問題研究[J]. 計算機工程與應用， 2011， 47（17）： 13-17.

[5]葛顯龍， 孔陽. 帶有時間窗的生鮮物流配送路徑優化研究[J]. 數學的實踐與認識， 2016， 46（12）： 78-87.

[6]崔巖， 張子祥， 時新， 等. 考慮顧客時間緊迫度的生鮮電商配送路徑優化問題[J]. 鄭州大學學報（工學版）， 2017， 38（6）： 59-63.

[7]姚源果， 賀盛瑜. 基于交通大數據的農產品冷鏈物流配送路徑優化研究[J]. 管理評論， 2019， 31（4）： 240-253.

[8]黃星星， 胡堅堃， 黃有方. 碳稅和碳限規則下生鮮農產品冷鏈配送路徑優化[J]. 上海海事大學學報， 2018， 39（1）： 74-79， 110. DOI： 10.13340/j.jsmu.2018.01.013.

[9]康凱， 韓杰， 普瑋， 等. 生鮮農產品冷鏈物流低碳配送路徑優化研究[J]. 計算機工程與應用， 2019， 55（2）： 259-265.

[10]AMORIM P，PARRAGH S N， SPERANDIO F，et al. A rich vehicle routing problem dealing with perishable food： a case study[J]. Top， 2014， 22（2）： 489-508. DOI： 10.1007/s11750-012-0266-4.

[11]李末芝. 市域連鎖零售業冷鏈物流配送路徑的優化研究[D]. 太原： 太原理工大學， 2011.

[12]張立毅， 王迎， 費騰， 等. 混沌擾動模擬退火蟻群算法低碳物流路徑優化[J]. 計算機工程與應用， 2017， 53（1）： 63-68， 102.

[13]KUMAR R S，KONDAPANENI K， DIXIT V，et al. Multi-objective modeling of production and pollution routing problem with time window： a self-learning particle swarm optimization approach[J]. Computers & Industrial Engineering， 2016， 99： 29-40. DOI： 10.1016/j.cie.2015.07.003.

[14]LI Jin，WANG Danping， ZHANG Jianghua. Heterogeneous fixed fleet vehicle routing problem based on fuel and carbon emissions[J]. Journal of Cleaner Production， 2018， 201（10）： 896-908.

[15]張如云， 劉清. 考慮低碳的城市配送車輛路徑優化模型研究[J]. 工業工程與管理， 2015， 20（4）： 29-34.

[16]張群， 顏瑞. 基于改進模糊遺傳算法的混合車輛路徑問題[J]. 中國管理科學， 2012， 20（2）： 121-128.

[17]呂寧， 謝如鶴， 劉廣海. 冷藏車開門時車內溫濕度變化實驗研究[J]. 制冷學報， 2013， 34（2）： 85-89. DOI： 10.3969/j.issn.0253-4339.2013.02.085.

[18]MUKHOPADHYAY S，MUKHERJEE R N， CHAUDHURI K S. Joint pricing and ordering policy for deteriorating inventory[J]. Computers & Industrial Engineering， 2004， 47： 339-349. DOI： 10.1016/j.cie.2004.06.007.

（編輯 趙勉）