“雙碳”背景下的冷鏈共同配送路徑優(yōu)化研究

郭丹蓉 陳碧玲

（福州理工學(xué)院,福建 福州 350000）

“雙碳”是碳達(dá)峰與碳中和的簡稱。2020年9月中國明確提出2030年“碳達(dá)峰”與2060年“碳中和”目標(biāo)[1]。作為僅次于工業(yè)、建筑業(yè)的第三大碳排放行業(yè)，交通運輸業(yè)的減排有利于實現(xiàn)國家的“雙碳”目標(biāo)[2]。物流服務(wù)離不開交通運輸，隨著物流需求的不斷增加，碳排放總量控制難度很大[3]。目前，發(fā)展低碳環(huán)保高效的物流已成為必然趨勢[4]，通過優(yōu)化配送路徑來實現(xiàn)不失為一種方法。車輛路徑問題由Dantzig等[5]在1959年提出，是指配送中心以最小配送成本的方式向客戶送貨。Bae等[6]提出了將多配送中心問題轉(zhuǎn)化為單配送中心再進行求解的思路。Solomon[7]考慮了時間窗對車輛路徑問題的限制。

當(dāng)前，社會經(jīng)濟的發(fā)展提高了對冷鏈物流的要求。共同配送是物流企業(yè)開展運輸時廣泛使用的一種模式，多家物流企業(yè)共享訂單與車輛，通過統(tǒng)一調(diào)度的方式實施配送活動[8]。丁艷[9]構(gòu)建多溫共配冷鏈物流車輛配送路徑優(yōu)化模型，結(jié)合蟻群算法求解問題。趙志學(xué)等[10]將道路擁堵因素融入冷鏈配送路徑優(yōu)化模型中。呂成瑤等[11]在考慮冷鏈配送時，用指數(shù)函數(shù)表示生鮮食品的腐敗規(guī)律。任騰等[12]基于現(xiàn)實配送情況，將顧客滿意度和道路擁堵狀況考慮在內(nèi)。

在“雙碳”背景下，冷鏈配送領(lǐng)域面臨新挑戰(zhàn)。本文將同時考慮冷鏈企業(yè)共同配送過程中的碳排放成本及冷鏈配送成本，在低碳與冷鏈之間尋找平衡點，以配送總成本最小化為目標(biāo)建立模型，通過粒子群算法得到優(yōu)化后的路徑，降低企業(yè)配送成本，并減少碳排放。

1 條件假設(shè)

本文基于“雙碳”背景，以減少碳排放為需求，以降低配送成本為目標(biāo)，在平衡碳排放和冷鏈共同配送的前提下，使用同一型號的若干冷鏈配送車從一個固定的冷鏈配送中心向多個客戶交付生鮮產(chǎn)品，通過規(guī)劃路徑，使配送過程總成本降至最低。

為了便于分析以下文本中的關(guān)鍵問題，首先做出了一些條件假設(shè)：（1）車輛從冷鏈配送中心出發(fā)，完成路線上所有客戶的交付，然后返回配送中心，車輛可以隨意部署；（2）送貨車輛數(shù)量沒有限制，車型相同，不允許超載；（3）節(jié)點（包括配送中心和客戶）的位置是已知的，客戶需求是一個固定的和已知的常數(shù)；（4）配送中心的貨物不短缺，每個客戶的需求都必須得到滿足，并且只能由一輛車提供一次服務(wù)。

2 模型構(gòu)建

2.1 車輛啟動成本

2.1.1 車輛的固定成本

固定成本是固定的和已知的常數(shù)，僅和車輛折舊等有關(guān)。記車輛固定成本為c11，則

2.2 貨損成本

2.3 制冷成本

2.4 碳排放成本

2.5 時間懲罰成本

2.6 總模型建立

目標(biāo)函數(shù)要求配送總成本最小，條件（1）表示配送車輛不超載；條件（2）表示服務(wù)完所有客戶；條件（3）和（4）表示每個客戶只被服務(wù)一次；條件（5）表示車輛起止點均為配送中心；條件（6）表示客戶的時間窗要求。

3 算法求解

3.1 粒子群算法原理

在問題的解空間中隨機獲得一組粒子作為初始解，每個粒子代表一個可行解；每個粒子隨機獲得包括大小和方向的初始速度，并開始局部搜索最優(yōu)解。對粒子進行連續(xù)迭代，每次迭代操作都導(dǎo)致粒子更新其位置pbest，并且總?cè)簩⑾鄳?yīng)地更新最佳值gbest，最終，粒子群將聚集在某區(qū)域作為最優(yōu)解。

該算法必須滿足以下要求：①所有尋優(yōu)的問題解都是在D維空間中搜索；②需要通過自適應(yīng)函數(shù)來評估所有粒子的質(zhì)量，以確定當(dāng)前位置的質(zhì)量；③每個粒子都可以記住其最佳位置；④ 每個粒子都有一個名為“速度”的屬性，可以清楚地定義飛行的距離和方向，并根據(jù)經(jīng)驗進行動態(tài)調(diào)整；⑤隨著進化的迭代，每個粒子的歷史最佳位置pbest,i和整個種群的歷史最佳位置gbest將不斷更新，從而導(dǎo)致粒子位置和速度的更新。因此，經(jīng)過t次迭代后，粒子i獲得的速度和位置可以表示為

其中，i=1,2,…，N，表示粒子群中的粒子i，總量為N；群體大小N是一個整數(shù)，規(guī)模N通常在20到40之間；k表示迭代次數(shù)；d=1,2,…，D，表示解空間的第d維；z1和z2是可以用來調(diào)整學(xué)習(xí)的最大步長的學(xué)習(xí)因子；r1和r2為區(qū)間[0,1]之間的隨機數(shù)，以增加隨機搜索的可能性；ω是一個非負(fù)慣性權(quán)重，用于調(diào)整解空間的搜索范圍，本文使用的粒子群算法中，將ω值設(shè)置如下：

3.2 粒子群算法步驟

步驟1：初始化粒子群。包括種群規(guī)模大小N，空間維度D，最大迭代次數(shù)T，隨機初始化粒子的速度和位置，以及設(shè)置初始個體的極值以及全局極值。

步驟2：計算適應(yīng)度值。基于適應(yīng)度函數(shù)計算每個粒子i的適應(yīng)度值。

步驟3：更新粒子個體和群體粒子的歷史最優(yōu)位置。將當(dāng)前適應(yīng)度值與在先前迭代期間獲得的個體和組的最優(yōu)適應(yīng)度值進行比較，并更新歷史最優(yōu)位置。

步驟4：更新粒子的速度和位置。計算群體中所有粒子的速度和位置，并確定其是否越界。

步驟5：確定終止條件。如果當(dāng)前迭代代數(shù)達(dá)到預(yù)設(shè)的最大迭代次數(shù)T，或者如果當(dāng)前粒子的適應(yīng)度值達(dá)到一定精度，則停止迭代并輸出結(jié)果；如果條件不滿足，則返回步驟2并繼續(xù)迭代操作。

4 案例分析

4.1 背景介紹

考慮單個配送中心使用同一型號的多輛冷鏈車輛向16個客戶運送貨物的問題，客戶要求送達(dá)時間均為6:00～ 9:00，可以接受的時間為5:30～ 9:30。配送中心O坐標(biāo)為（0，0），相關(guān)參數(shù)如表1和表2所示。

表1 各客戶相關(guān)參數(shù)

表2 其他相關(guān)參數(shù)

其中，各節(jié)點之間的行駛里程按兩點間的直線距離來計，即點A(xa,ya)與點B(xb,yb)之間的距離為

4.2 案例求解

該配送中心車輛調(diào)度人員根據(jù)自己豐富的工作經(jīng)驗，大致確定出初始的配送路徑方案。使用7輛車進行配送：O—P4—P16—P3—O；O—P9—P5—P10—O；O—P2—P7—P11—O；O—P12—P13—O；O—P1—P14—O；O—P6—P8—O；O—P15—O。該配送方案總成本為34 114.0元，其中，啟動成本2 871.7元，貨損成本9 399.5元，制冷成本14 742.1元，碳排放成本7 100.7元，時間懲罰成本0.0元。

使用粒子群算法求解案例，結(jié)果如下。使用7輛車：O—P4—P16—P3—O；O—P10—P9—P5—O；O—P2—P7—P11—O；O—P13—P12—O；O—P14—P1—O；O—P6—P8—O；O—P15—O。總配送成本33 469.0元，其中啟動成本2 879.2元，貨損成本9 146.7元，制冷成本14 385.5元，碳排放成本7 057.6元，時間懲罰成本0.0元。

優(yōu)化方案中，配送總成本降低了645.0元。從長期來看，該方法可以降低配送成本和碳排放成本，滿足客戶需求，為“雙碳”目標(biāo)做出一定的貢獻。

5 結(jié)語

在“雙碳”背景下，減少碳排放是企業(yè)發(fā)展的重中之重，而更好地發(fā)展冷鏈配送，對于冷鏈物流企業(yè)而言既是挑戰(zhàn)也是機遇。本文基于“雙碳”政策，建立了低碳冷鏈共同配送路徑規(guī)劃模型，考慮了冷鏈配送中涉及的啟動成本、制冷成本、貨損成本、碳排放成本和時間懲罰成本等因素，并以總成本最小化為目標(biāo)。利用粒子群算法解決此問題，并通過實例驗證了該算法的可行性。