楊曉華 郭健全

摘 要：針對(duì)生鮮品因易腐易損性而產(chǎn)生的高頻次物流配送及不確定需求與退貨量的問(wèn)題，提出了模糊環(huán)境下多周期生鮮閉環(huán)物流網(wǎng)絡(luò)系統(tǒng)，以實(shí)現(xiàn)最小系統(tǒng)成本、最優(yōu)設(shè)施選址與最佳配送路徑的多決策安排。為求解系統(tǒng)對(duì)應(yīng)的模糊混合整數(shù)線性規(guī)劃（FMILP）模型，首先將生鮮需求量和退貨量設(shè)定為三角模糊值，其次運(yùn)用模糊機(jī)會(huì)約束規(guī)劃方法將模糊約束等價(jià)變換為清晰式，最后利用遺傳算法（GA）和粒子群優(yōu)化（PSO）算法搜索案例的最優(yōu)解。實(shí)驗(yàn)結(jié)果表明，多周期閉環(huán)系統(tǒng)比單周期更能兼顧多決策規(guī)劃，同時(shí)三角模糊量的置信水平變化對(duì)企業(yè)最優(yōu)運(yùn)作有著顯著影響，進(jìn)而為相關(guān)決策者提供借鑒。

關(guān)鍵詞：生鮮;多周期配送;閉環(huán)物流網(wǎng)絡(luò);模糊機(jī)會(huì)約束規(guī)劃;混合整數(shù)線性規(guī)劃;遺傳算法;粒子群優(yōu)化算法

Abstract： Concerning the high frequency logistics distribution of fresh products due to the products perishability and vulnerability， as well as the uncertainty of demand and return， a multi-period closed-loop logistics network for fresh products with fuzzy variables was constructed to achieve the multi-decision arrangement of minimum system cost， optimal facility location and optimal delivery route. In order to solve the Fuzzy Mixed Integer Linear Programming （FMILP） model corresponding to the system， firstly， the amounts of demand and return were defined as triangular fuzzy parameters; secondly， the fuzzy constraints were transformed into crisp formula by using fuzzy chance constrained programming method; finally， the optimal solution of case was obtained by using Genetic Algorithm （GA） and Particle Swarm Optimization （PSO） algorithm. The experimental results show that multi-period closed-loop system performs better than single-period system in the aspect of multi-decision programming， meanwhile， the confidence levels of triangular fuzzy parameters have significant influence on the optimal operation of enterprise， thus providing a reference for relevant decision makers.

Key words： fresh product; multi-period delivery; closed-loop logistics network; fuzzy chance constrained programming; mixed integer linear programming; Genetic Algorithm （GA）; Particle Swarm Optimization （PSO） algorithm

0 引言

隨著中國(guó)生鮮電商的迅猛發(fā)展，越來(lái)越多研究聚焦于建設(shè)生鮮物流體系[1]。由于生鮮產(chǎn)品的易腐易損、需求量波動(dòng)大、退貨率高等特征[2]，產(chǎn)品的儲(chǔ)運(yùn)、保鮮成本以及生鮮企業(yè)的設(shè)施運(yùn)營(yíng)成本很高[1，3]。科學(xué)、高效的閉環(huán)物流網(wǎng)絡(luò)設(shè)計(jì)有利于減少配送過(guò)程的物流損失、提高客戶滿意度和實(shí)現(xiàn)資源的優(yōu)化利用[3-4]，然而，現(xiàn)有生鮮產(chǎn)品的國(guó)內(nèi)外研究大多探討其運(yùn)營(yíng)模式、保鮮技術(shù)、配送系統(tǒng)等方面[5-7]，對(duì)于閉環(huán)物流網(wǎng)絡(luò)的系統(tǒng)性分析尚有不足，因此，本文以生鮮產(chǎn)品為研究對(duì)象，構(gòu)建了閉環(huán)物流網(wǎng)絡(luò)系統(tǒng)及其數(shù)學(xué)模型，以期為該類研究提供理論支持，并為相關(guān)生鮮企業(yè)提供實(shí)踐指導(dǎo)。

參數(shù)不確定性是閉環(huán)物流網(wǎng)絡(luò)中復(fù)雜卻非常重要的研究因素[8]，其中需求量和退貨量是最重要的兩種不確定參數(shù)[9]。Khatami等[10]在考慮閉環(huán)物流網(wǎng)絡(luò)需求與退貨量不確定下，采用Bender分解法獲取混合整數(shù)模型的最優(yōu)解，并通過(guò)Cholesky分解法分析不同產(chǎn)品需求的相關(guān)性。模型中將回收比率、再生率和廢棄率假定為常數(shù)，然而在實(shí)際規(guī)劃中，此類參數(shù)是難以預(yù)測(cè)的。Hamdouch等[11]在設(shè)定成本函數(shù)的前提下，構(gòu)建了分散式閉環(huán)供應(yīng)鏈網(wǎng)絡(luò)，通過(guò)調(diào)整需求和退貨分布，來(lái)量化閉環(huán)供應(yīng)鏈中出貨量、價(jià)格和預(yù)期利潤(rùn)的均衡關(guān)系，進(jìn)而表明隨機(jī)性對(duì)預(yù)期收益的影響。上述文獻(xiàn)，雖然闡明了不確定需求與退貨對(duì)系統(tǒng)運(yùn)作具有影響，但均缺乏構(gòu)建不確定規(guī)劃模型，以具體分析確切的閉環(huán)物流網(wǎng)絡(luò)運(yùn)作。

隨機(jī)規(guī)劃與模糊規(guī)劃是解決生鮮物流網(wǎng)絡(luò)問(wèn)題常見的兩種不確定規(guī)劃[12-13]。然而隨機(jī)規(guī)劃方法需要生鮮參數(shù)的精確分布情況[14]，由此對(duì)此類模型求解提出了實(shí)質(zhì)性的挑戰(zhàn);越來(lái)越多研究?jī)A向于使用模糊規(guī)劃方法，利用企業(yè)制定的三角模糊數(shù)代表參數(shù)發(fā)生的可能性，并且在去模糊化過(guò)程中，通過(guò)無(wú)差異化組合替代模糊約束條件，更適用于企業(yè)實(shí)際運(yùn)作[15]。

多周期物流配送對(duì)閉環(huán)物流網(wǎng)絡(luò)系統(tǒng)中設(shè)施選址與配送路徑規(guī)劃具有重要影響[15-17]。Zhang等[18]基于變分不等式和互補(bǔ)理論，探討了多周期規(guī)劃下包含碳排放約束的閉環(huán)供應(yīng)鏈網(wǎng)絡(luò)，并采用改進(jìn)的投影壓縮算法對(duì)模型進(jìn)行求解。此模型分析了多周期規(guī)劃的全局優(yōu)勢(shì)，但缺乏進(jìn)一步探討在多周期各個(gè)階段維度下的系統(tǒng)運(yùn)作問(wèn)題。Mohammed等[19]為實(shí)現(xiàn)最優(yōu)閉環(huán)物流網(wǎng)絡(luò)的設(shè)計(jì)和規(guī)劃決策，提出了一個(gè)包含多周期（三個(gè)時(shí)間段）、多產(chǎn)品、兩種不確定因素的閉環(huán)供應(yīng)鏈網(wǎng)絡(luò)模型。多周期基于場(chǎng)景的隨機(jī)方法多適用于大規(guī)模網(wǎng)絡(luò)，而此模型僅考慮較少的場(chǎng)景分析，使得結(jié)論的信度與效度有所不足。此外，當(dāng)前文獻(xiàn)較少考慮到需求與退貨不確定下，如何將多周期規(guī)劃運(yùn)用于生鮮閉環(huán)物流網(wǎng)絡(luò)設(shè)計(jì)。

本文構(gòu)建了多周期配送下的生鮮閉環(huán)物流網(wǎng)絡(luò)系統(tǒng)及模糊混合整數(shù)線性規(guī)劃（Fuzzy Mixed Integer Linear Programming， FMILP）模型，在此基礎(chǔ)上，利用模糊機(jī)會(huì)約束規(guī)劃將模型等價(jià)轉(zhuǎn)化為清晰式。為實(shí)現(xiàn)最小系統(tǒng)成本決策，將物流總成本作為模型的目標(biāo)函數(shù)，包含設(shè)施的固定成本、設(shè)施維持成本與運(yùn)營(yíng)成本、企業(yè)的信息處理成本、產(chǎn)品的運(yùn)輸成本與碳交易成本。以某生鮮企業(yè)在上海市6個(gè)區(qū)配送蔬菜為實(shí)例，利用遺傳算法（Genetic Algorithm， GA）和粒子群優(yōu)化（Particle Swarm Optimization， PSO）算法搜索最優(yōu)解，從而驗(yàn)證了模型的信度與效度。

1 問(wèn)題描述

由于生鮮品的易腐易損、需求量波動(dòng)大、退貨率高等特性，閉環(huán)物流網(wǎng)絡(luò)的單周期配送遠(yuǎn)遠(yuǎn)不能滿足生鮮企業(yè)實(shí)際運(yùn)作，而多周期物流配送系統(tǒng)將生鮮產(chǎn)品按易腐度、需求波動(dòng)度等因素來(lái)劃分多周期的不同階段，這使得企業(yè)更能精確把握消費(fèi)者動(dòng)態(tài)需求、設(shè)施長(zhǎng)期選址、配送車輛路徑優(yōu)化等信息，進(jìn)而實(shí)現(xiàn)生鮮閉環(huán)物流網(wǎng)絡(luò)系統(tǒng)的配置優(yōu)化。

同樣地，生鮮閉環(huán)物流網(wǎng)絡(luò)設(shè)計(jì)中還存在參數(shù)值不確定性的復(fù)雜問(wèn)題。企業(yè)很難獲得生鮮行業(yè)或者特定生鮮品的統(tǒng)計(jì)數(shù)據(jù)及其函數(shù)分布，產(chǎn)品參數(shù)往往由企業(yè)決策者制定，并以一定置信水平下的模糊量形式代入決策分析中[20-21]。通過(guò)對(duì)各區(qū)域消費(fèi)者在每個(gè)周期需求與退貨的統(tǒng)計(jì)分析，企業(yè)評(píng)估各區(qū)域的置信區(qū)間上下界、極可能值及置信度，從而為生鮮閉環(huán)物流網(wǎng)絡(luò)系統(tǒng)作出合理的決策安排。

2 生鮮多周期閉環(huán)物流網(wǎng)絡(luò)模型

2.1 生鮮閉環(huán)物流網(wǎng)絡(luò)結(jié)構(gòu)

生鮮閉環(huán)物流網(wǎng)絡(luò)結(jié)構(gòu)如圖1所示。在一個(gè)周期內(nèi)，中轉(zhuǎn)中心分配倉(cāng)作為生鮮閉環(huán)物流網(wǎng)絡(luò)的分銷商，需對(duì)從生鮮供應(yīng)商采購(gòu)的產(chǎn)品進(jìn)行分類、包裝、貼簽等作業(yè)。配送運(yùn)輸過(guò)程中，乘坐配送車輛的工作人員將產(chǎn)品按企業(yè)預(yù)估的訂單需求配送至正向物流終端（各區(qū)域配送網(wǎng)點(diǎn)），與此同時(shí)，完成提取上一周期退貨品的回收工作。在逆向物流中，中轉(zhuǎn)中心回收倉(cāng)將各區(qū)域配送網(wǎng)點(diǎn)的退貨品整合送至生鮮供應(yīng)商處，分類處理，將具有一定殘值的退貨品二次加工，后續(xù)把生鮮殘?jiān)椭翉U物處理廠焚燒或填埋[4]。

2.2 FMILP模型

2.2.1 模型假設(shè)

根據(jù)本文生鮮閉環(huán)物流網(wǎng)絡(luò)特點(diǎn)，作出如下模型假設(shè)：

1）各區(qū)域第三方配送中心有大、小存儲(chǔ)容量。

2）退貨生鮮品殘值低、保質(zhì)期短及二次銷售難度大，其回收處理方式只考慮第三方廢棄處理工序。

3）供應(yīng)商、中轉(zhuǎn)中心、各區(qū)域配送網(wǎng)點(diǎn)和廢物處理廠的候選位置與數(shù)量已知。

4）生鮮企業(yè)的信息處理成本與消費(fèi)者的需求量和退貨量成線性關(guān)系。

5）中轉(zhuǎn)中心到生鮮供應(yīng)商及廢物處理廠的運(yùn)輸成本與運(yùn)輸量和運(yùn)輸距離成線性關(guān)系;各區(qū)域配送網(wǎng)點(diǎn)之間的運(yùn)輸成本、中轉(zhuǎn)中心到各區(qū)域配送網(wǎng)點(diǎn)的運(yùn)輸成本只與運(yùn)輸距離成線性關(guān)系[16-17]。

6）各個(gè)區(qū)域至少有一個(gè)配送中心網(wǎng)點(diǎn)，且網(wǎng)點(diǎn)僅提供該區(qū)域服務(wù)，不能跨區(qū)域配送。

7）若個(gè)別區(qū)域的候選網(wǎng)點(diǎn)與相鄰區(qū)域網(wǎng)點(diǎn)距離相差太小，則僅開放需求量最大的區(qū)域?qū)?yīng)的配送網(wǎng)點(diǎn)，其他區(qū)域不額外設(shè)立網(wǎng)點(diǎn)。

2.2.2 變量描述

s代表供應(yīng)商，s∈{1，2，…，S};q代表中轉(zhuǎn)中心，q∈{1，2，…，Q};i代表區(qū)域，i∈{1，2，…，I};j代表配送網(wǎng)點(diǎn)， j∈{1，2，…，J};e代表生鮮企業(yè)，e∈{1，2，…，E};l代表廢物處理廠，l∈{1，2，…，L};t代表配送周期，t∈{1，2，…，T};k代表配送網(wǎng)點(diǎn)存儲(chǔ)容量，k∈{s，l|k∈K};v代表運(yùn)輸車輛，v∈V;R代表運(yùn)輸路線，R=1，2，…，RR∈N*此句感覺不太準(zhǔn)確，R=……R，R最后等于R，表述不規(guī)范，需調(diào)整?；貜?fù)：此處調(diào)整為R∈N*。。

2.2.3 參數(shù)

Fq代表中轉(zhuǎn)中心的固定建設(shè)成本;maqt代表第t周期中轉(zhuǎn)中心的設(shè)施維持成本;makijt代表第t周期區(qū)域i配送網(wǎng)點(diǎn)j（類型k）的設(shè)施維持成本;opqt代表第t周期中轉(zhuǎn)中心的運(yùn)營(yíng)成本;opkijt代表第t周期區(qū)域i配送網(wǎng)點(diǎn)j（類型k）的運(yùn)營(yíng)成本;uipet代表第t周期生鮮企業(yè)的單位信息處理成本;utsq，utqij，utiji′j′，utql分別代表兩個(gè)物流節(jié)點(diǎn)間的單位運(yùn)輸成本;ndsq，ndqij，ndiji′j′，ndql分別代表兩個(gè)物流節(jié)點(diǎn)間的距離;hcq，hckij，hcl分別代表各節(jié)點(diǎn)的最大處理能力;δ代表退貨品的廢棄率;K代表單位車輛最大運(yùn)載量;a代表交通工具碳排放系數(shù);b/c代表交通工具能源消耗系數(shù);Eq代表生鮮企業(yè)的碳排放配額;Up代表單位碳排放交易價(jià)格[4];α為三角模糊量的置信水平;Dit代表第t周期區(qū)域i附近的消費(fèi)者需求量，模糊值;Rit代表第t周期區(qū)域i附近的消費(fèi)者退貨量，模糊值。

2.2.4 決策變量

1）Xqt，Xijt，Xlt為0-1變量，若第t周期開放中轉(zhuǎn)中心、區(qū)域i配送網(wǎng)點(diǎn)j、廢物處理廠，則Xqt=1，Xijt=1，Xlt=1;否則Xqt=0，Xijt=0，Xlt=0。

2）XRvt為0-1變量，若第t周期車輛v在第R路線運(yùn)輸生鮮產(chǎn)品，則XRvt=1;否則XRvt=0。

3）XRqijt，XRiji′j′t，XRijqt為0-1變量，若在第t周期的第R路線上，分別將產(chǎn)品從中轉(zhuǎn)中心分配倉(cāng)運(yùn)輸?shù)脚渌途W(wǎng)點(diǎn)、區(qū)域i配送網(wǎng)點(diǎn)j運(yùn)輸?shù)絽^(qū)域i配送網(wǎng)點(diǎn)j，以及從配送網(wǎng)點(diǎn)運(yùn)輸?shù)街修D(zhuǎn)中心回收倉(cāng)，則XRqijt=1，XRiji′j′t=1，XRijqt=1;否則XRqijt=0，XRiji′j′t=0，XRijqt=0。

4）Ysqt，Yqit，Yiqt，Yqst，Yqlt為整數(shù)變量，表示第t周期節(jié)點(diǎn)間的運(yùn)輸量。

2.2.5 數(shù)學(xué)模型的建立

單周期的固定成本與設(shè)施維持成本Z1（single）：

多周期的固定成本與設(shè)施維持成本Z1（multiple）：

2.3 多周期下模糊機(jī)會(huì)約束清晰化

生鮮閉環(huán)物流網(wǎng)絡(luò)系統(tǒng)中需求量Dit以及退貨量Rit為模糊參數(shù)，因此本文求解的不再是約束清晰的混合整數(shù)線性規(guī)劃（Mixed Integer Linear Programming， MILP）問(wèn)題。為有效規(guī)避模糊不確定參數(shù)對(duì)此類問(wèn)題的影響，宜采用模糊機(jī)會(huì)約束規(guī)劃（Fuzzy Chance Constrained Programming請(qǐng)補(bǔ)充FCCP的英文全稱， FCCP）方法[20]。首先將Dit與Rit設(shè)為三角模糊參數(shù)，記Dit=（Di1t，Di2t，Di3t）、Rit=（Ri1t，Ri2t，Ri3t），其中：Di1t（Ri1t）與Di3t（Ri3t）分別為生鮮企業(yè)在第t周期制定的置信水平α的上下界，Di2t（Ri2t）為α的最可能值。其次應(yīng)保證約束條件成立概率控制在企業(yè)制定的置信水平之上，通過(guò)對(duì)模糊對(duì)應(yīng)式進(jìn)行等價(jià)變換，從而使不確定模糊規(guī)劃轉(zhuǎn)化為清晰的確定規(guī)劃。每個(gè)周期對(duì)應(yīng)的模糊隸屬函數(shù)如下表示：

3 算法設(shè)計(jì)

由于系統(tǒng)的多周期設(shè)計(jì)、模糊機(jī)會(huì)約束規(guī)劃、復(fù)雜約束條件等因素，本文求解的是NP-hard問(wèn)題最優(yōu)解，利用元啟發(fā)式算法求解此類模型被認(rèn)為是有效的解決方式[23-24]。

遺傳算法（GA）和粒子群優(yōu)化（PSO）算法都是搜索近似最優(yōu)解的元啟發(fā)式優(yōu)化技術(shù)[4，25]，均被廣泛應(yīng)用于解決大規(guī)模物流網(wǎng)絡(luò)構(gòu)建、選址庫(kù)存路徑等問(wèn)題[4，17]。二者的迭代過(guò)程具有互補(bǔ)優(yōu)勢(shì)，而單一運(yùn)算存在一定的局限性，因此，本文同時(shí)將GA與PSO運(yùn)用到算例求解中，進(jìn)一步驗(yàn)證模型的信度與效度。

3.1 GA

在模型中，式（12）～（13）、式（18） ～ （20）、式（26）、式（28） ～ （29）的決策變量是根據(jù)模糊機(jī)會(huì)約束規(guī)劃，在設(shè)定置信水平中以特定置信區(qū)間值代入模型計(jì)算，基于此，本文僅需要探討包含0-1變量約束的整數(shù)線性規(guī)劃問(wèn)題。

為了使經(jīng)過(guò)遺傳算子操作后的染色體編碼依舊滿足模型約束，同時(shí)提高收斂速度和結(jié)果精確度，本文選擇合適的GA[26]，有效地設(shè)計(jì)遺傳算子，使約束條件被合理地運(yùn)用到GA各個(gè)環(huán)節(jié)中。

步驟1 染色體編碼與初始化。

本文的染色體為IJ+Q+L+V（R）+SQ（R）+QIJ（R）+IJI′J′（R）+IJQ（R）+QS+SI數(shù)組。一類染色體為二進(jìn)制串（1或0），表示設(shè)施J、P、L是否建立，車輛V是否在路線R上運(yùn)輸以及路線R上節(jié)點(diǎn)是否有被選擇;二類染色體對(duì)應(yīng)的數(shù)值表示實(shí)例中節(jié)點(diǎn)間的生鮮運(yùn)輸量。在GA的過(guò)程中，自動(dòng)確定各區(qū)域配送網(wǎng)點(diǎn)IJ、中轉(zhuǎn)中心Q、及廢物處理廠L是否開放以及規(guī)劃相應(yīng)路線。

步驟2 適應(yīng)度評(píng)估與選擇。

本文引入懲罰函數(shù)將0-1變量約束（式（9）～ （11）、式（17）、式（21）～（25）轉(zhuǎn)化為無(wú)約束規(guī)劃。在GA中，變量A代表模型的0-1決策變量a，Z（A）為目標(biāo)函數(shù)，hi（A）≤0，i=1，2，…，n代表所有不等式約束，lj（A）=0， j=n+1，n+2，…，nn代表所有等式約束。定義函數(shù)U（A）如下所示：

步驟3 交叉與變異。

交叉使子代同時(shí)繼承父母代的基因，保持更高的適應(yīng)性。變異使某個(gè)基因以一定概率隨機(jī)發(fā)生突變。交叉和變異不僅產(chǎn)生新個(gè)體，而且提高種群局部搜索能力。

步驟4 終止條件。

若達(dá)到最大進(jìn)化代數(shù)，則終止算法，并將輸出的最高適應(yīng)度值的個(gè)體作為模型是我最優(yōu)解此句不通順，請(qǐng)調(diào)整將輸出最高適應(yīng)度值的個(gè)體作為模型的最優(yōu)解;否則轉(zhuǎn)向步驟2。

4 算例結(jié)果及分析

4.1 實(shí)驗(yàn)數(shù)據(jù)來(lái)源

本文以上海市某生鮮企業(yè)新鮮直送蔬菜業(yè)務(wù)為例，配送網(wǎng)點(diǎn)選址范圍在寶山區(qū)、嘉定區(qū)、浦東新區(qū)、靜安區(qū)、徐匯區(qū)以及青浦區(qū)。假設(shè)配送車輛運(yùn)輸蔬菜產(chǎn)品單周期長(zhǎng)度為一周，多周期中一周分為三個(gè)周期（精確的周期長(zhǎng)度劃分依照產(chǎn)品易腐程度、市場(chǎng)環(huán)境以及企業(yè)實(shí)際運(yùn)營(yíng)變化時(shí)間為準(zhǔn)）。由企業(yè)調(diào)研可知，蔬菜供應(yīng)商坐標(biāo)為（55.7，3.7）;候選的中轉(zhuǎn)中心坐標(biāo)分別為（36.1，21.1）、（42.2，9.9）、（43.8，8.3）、（36.4，11.4）;候選的廢物處理廠坐標(biāo)為（60.4，2.6）、（59.6，7.7）。各區(qū)域的候選配送網(wǎng)點(diǎn)坐標(biāo)及各區(qū)域每個(gè)周期的需求與退貨三角模糊量如表1所示，其他參數(shù)值如表2。

需求與退貨的三角模糊量受生鮮企業(yè)制定的置信水平（αd、αr）影響。考慮到αd與αr值對(duì)企業(yè)最優(yōu)決策的影響，本文采用均勻分布法對(duì)αd與αr分別賦值70%、80%、90%，并將之兩兩交叉，組成如下9種場(chǎng)景：場(chǎng)景1（αd=70%，αr=70%）、場(chǎng)景2（αd=70%，αr=80%）、場(chǎng)景3（αd=70%，αr=90%）、場(chǎng)景4（αd=80%，αr=70%）、場(chǎng)景5（αd=80%，αr=80%）、場(chǎng)景6（αd=80%，αr=90%）、場(chǎng)景7（αd=90%，αr=70%）、場(chǎng)景8（αd=90%，αr=80%）、場(chǎng)景9（αd=90%，αr=90%）置信水平兩兩交叉組成的9種場(chǎng)景9種場(chǎng)景是如何劃分的？不好理解。需在正文中補(bǔ)充說(shuō)明一下，這樣圖3中的場(chǎng)景序號(hào)才好解讀，如圖3所示。

4.2 結(jié)果分析

本文利用Matlab 2016a編寫GA與PSO代碼，迭代計(jì)算所得的多周期與單周期結(jié)果以及算法之間的差值如表3所示。

在預(yù)設(shè)最大迭代次數(shù)150代情況下，GA和PSO分別在第124代和第106代獲得最優(yōu)適應(yīng)度曲線。迭代結(jié)果表明，二者差距Gap僅為-1.19%～0.21%，由此驗(yàn)證了模型的信度與效度。此外，不論多周期還是單周期情況，通過(guò)GA獲得的目標(biāo)值大體上都比PSO值小。為方便分析，后文的設(shè)施選址方案、不同置信水平系統(tǒng)成本對(duì)比情況，以及配送路徑圖均為GA迭代所得，如圖4～5所示。

從模糊角度分析可知：①在任意周期下，需求與退貨量隨著對(duì)應(yīng)三角模糊量的置信水平αd與αr增加而增加，遞增的在途運(yùn)輸量直接使運(yùn)輸成本及碳交易成本呈上升趨勢(shì)，進(jìn)而導(dǎo)致了系統(tǒng)總成本的不斷上升。②生鮮退貨量是需求量的11%～19%，因此在相似的置信水平變化下，需求置信水平對(duì)物流成本與設(shè)施選址影響更大。③當(dāng)αd=90%，αr=70%時(shí)，各區(qū)域的需求量均大幅增加，但由于部分區(qū)域配送網(wǎng)點(diǎn)的容量存儲(chǔ)有限，配送網(wǎng)點(diǎn)選址進(jìn)行了適當(dāng)調(diào)整，同時(shí)配送路徑也作出相應(yīng)的靈活變動(dòng)，這使得優(yōu)化成本出現(xiàn)驟減情況。綜上，生鮮閉環(huán)物流網(wǎng)絡(luò)的模糊機(jī)會(huì)約束規(guī)劃充分考慮了不確定環(huán)境下的系統(tǒng)設(shè)計(jì)，同時(shí)三角模糊量的置信水平變化對(duì)企業(yè)最優(yōu)運(yùn)作有著顯著影響。

從多周期角度分析可得：①?gòu)谋?～4、圖3～5可知，考慮到生鮮產(chǎn)品的多周期配送，細(xì)分的單位周期的需求量與退貨量對(duì)于設(shè)施的選址更加具有針對(duì)性，從而減少了開放更多的設(shè)施而產(chǎn)生的固定建設(shè)成本，因此多周期系統(tǒng)成本比單周期少。②隨著置信水平不斷增大，因開放設(shè)施而產(chǎn)生的固定建設(shè)成本比重不斷減少，多周期與單周期的差值呈遞減趨勢(shì)。③觀測(cè)上述最優(yōu)置信水平下不同周期的路徑圖可知，多周期各階段的設(shè)施量與配送路線均低于單周期，同時(shí)，多周期三個(gè)階段的中轉(zhuǎn)中心與廢物處理廠選址雖一致，但各自的配送路線根據(jù)各區(qū)域?qū)嶋H需求而作出相應(yīng)調(diào)整，更符合企業(yè)實(shí)際決策。綜上，考慮生鮮閉環(huán)物流網(wǎng)絡(luò)的多周期配送使得企業(yè)在保證成本最低情況下，實(shí)現(xiàn)最優(yōu)設(shè)施選址及路徑規(guī)劃。

5 結(jié)語(yǔ)

本文以生鮮品為研究對(duì)象，將解決多周期配送、需求與退貨模糊不確定性問(wèn)題，以及實(shí)現(xiàn)最小系統(tǒng)成本、最優(yōu)設(shè)施選址與最佳配送路徑的最優(yōu)決策作為出發(fā)點(diǎn)，設(shè)計(jì)了生鮮品閉環(huán)物流網(wǎng)絡(luò)系統(tǒng)。GA與PSO的實(shí)驗(yàn)分析表明，多周期物流配送系統(tǒng)比單周期更能均衡多決策安排，同時(shí)模糊機(jī)會(huì)約束規(guī)劃對(duì)企業(yè)最優(yōu)運(yùn)作有著合理的借鑒意義，從而驗(yàn)證了本文模型的有效性與可行性。

本文通過(guò)模糊機(jī)會(huì)約束規(guī)劃來(lái)分析需求與退貨不確定性，然而求解閉環(huán)物流網(wǎng)絡(luò)的不確定規(guī)劃方法是多樣的，如何將隨機(jī)規(guī)劃、魯棒優(yōu)化等方法運(yùn)用到模型，還待進(jìn)一步分析。

參考文獻(xiàn) （References）

[1] YANG B， ZHANG D. Research on coordination fresh product supply chain under new retailing model[J]. Computational Intelligence and Intelligent Systems， 2018， 874： 439-445.

[2] 王晶，劉昊天，趙然.基于食品安全的生鮮食品冷鏈運(yùn)營(yíng)優(yōu)化研究[J].系統(tǒng)工程理論與實(shí)踐，2018，38（1）：122-134.（WANG J， LIU H T， ZHAO R. The optimization of cold chain operation based on fresh food safety[J]. Systems Engineering-Theory and Practice， 2018， 38（1）： 122-134.）

[3] WU Q， MU Y， FENG Y. Coordinating contracts for fresh product outsourcing logistics channels with power structures [J]. International Journal of Production Economics， 2015， 160： 94-105.

[4] 郭健全，王心月.碳交易下生鮮電商跨區(qū)域閉環(huán)物流網(wǎng)絡(luò)及路徑[J].計(jì)算機(jī)集成制造系統(tǒng)，2017，23（4）：874-882.（GUO J Q， WANG X Y. Network and route planning of cross-regional closed-loop logistics for fresh food e-commerce under environment of carbon trading[J]. Computer Integrated Manufacturing Systems， 2017， 23（4）： 874-882.）

[5] FERREIRA G， ARRUDA E， MARUJO L. Inventory management of perishable items in long-term humanitarian operations using Markov decision processes [J]. International Journal of Disaster Risk Reduction， 2018， 31： 460-469.

[6] SUN X， BALDWIN E， BAI J. Applications of gaseous chlorine dioxide on postharvest handling and storage of fruits and vegetables-a review [J]. Food Control， 2019， 95： 18-26.

[7] 張文峰，梁凱豪.生鮮農(nóng)產(chǎn)品冷鏈物流網(wǎng)絡(luò)節(jié)點(diǎn)和配送的優(yōu)化[J].系統(tǒng)工程，2017，35（1）：119-123.（ZHANG W F， LIANG K H. Optimization of cold-chain network nodes and delivery for fresh agricultural products[J]. Systems Engineering， 2017， 35（1）： 119-123.）

[8] CUI Y， GUAN Z， SAIF U， et al. Close loop supply chain network problem with uncertainty in demand and returned products： genetic artificial bee colony algorithm approach[J]. Journal of Cleaner Production， 2017， 162： 717-742.

[9] GOVINDAN K， SOLEIMANI H， KANNAN D. Reverse logistics and closed-loop supply chain： a comprehensive review to explore the future[J]. European Journal of Operational Research， 2015， 240（3）： 603-626.

[10] KHATAMI M， MAHOOTCHI M， ZANJIRANI FARAHANI R. Benders decomposition for concurrent redesign of forward and closed-loop supply chain network with demand and return uncertainties [J]. Transportation Research Part E： Logistics and Transportation Review， 2015， 79： 1-21.

[11] HAMDOUCH Y， QIANG P Q， GHOUDI K. A closed-loop supply chain equilibrium model with random and price-sensitive demand and return [J]. Networks and Spatial Economics， 2017， 17（2）：459-503.

[12] DAI Z， AQLAN F， ZHENG X， et al. A location-inventory supply chain network model using two heuristic algorithms for perishable products with fuzzy constraints [J]. Computers and Industrial Engineering， 2018， 119： 338-352.

[13] AHUMADA O， VILLALOBOS J R， MASON A N. Tactical planning of the production and distribution of fresh agricultural products under uncertainty[J]. Agricultural Systems， 2012， 112（13）： 17-26.

[14] YU Y， CHEN X， GAO W， et al. Stochastic leaching analysis on cementitious materials considering the influence of material uncertainty [J]. Construction and Building Materials， 2018， 184： 186-202.

[15] SOHEILA K. Simulation based evolutionary algorithms for fuzzy chance-constrained biogas supply chain design[J]. Applied Energy， 2019， 236： 183-195.

[16] 曹鋒，郭健全，劉欣欣.考慮碳排放的多周期醫(yī)藥逆向物流網(wǎng)絡(luò)聯(lián)建研究[J].華東師范大學(xué)學(xué)報(bào)（自然科學(xué)版），2017（2）：52-60.（CAO F， GUO J Q， LIU X X. A joint construction research for multi-period medicine reverse logistics network with consideration of carbon emissions[J]. Journal of East China Normal University （Natural Science）， 2017（2）： 52-60.）

[17] GUO J， JIANG F， MITSUON G. Dynamic joint construction and optimal strategy of multi-objective multi-period multi-stage reverse logistics network： a case study of lead battery in Shanghai[J]. Procedia Manufacturing， 2018， 17： 1171-1178.

[18] ZHANG G， ZHONG Y， SUN H， et al. Multi-period closed-loop supply chain network equilibrium with carbon emission constraints[J]. Resources Conservation and Recycling， 2015， 104（B）： 354-365.

[19] MOHAMMED F， SELIM S Z， HASSAN A， et al. Multi-period planning of closed-loop supply chain with carbon policies under uncertainty[J]. Transportation Research Part D： Transport and Environment， 2017， 51： 146-172.

[20] 狄衛(wèi)民，馬祖軍，代穎.制造/再制造集成物流網(wǎng)絡(luò)模糊優(yōu)化設(shè)計(jì)方法[J].計(jì)算機(jī)集成制造系統(tǒng)，2008，14（8）：1472-1480.（DI W M， MA Z J， DAI Y. Fuzzy optimal design approach for integrated manufacturing/remanufacturing logistics network[J]. Computer Integrated Manufacturing Systems， 2008， 14（8）： 1472-1480.）

[21] GUO J， WANG X， FAN S， et al. Forward and reverse logistics network and route planning under the environment of low-carbon emissions[J]. Computers and Industrial Engineering， 2017， 106（C）： 351-360.

[22] 劉寶碇.隨機(jī)規(guī)劃與模糊規(guī)劃[M].北京：清華大學(xué)出版社，1998：164-183.（LIU B D. Stochastic Programming and Fuzzy Programming[M]. Beijing： Tsinghua University Press， 1998： 164-183.）

[23] DIVEEV A I， BOBR O V. Variational genetic algorithm for NP-hard scheduling problem solution [J]. Procedia Computer Science， 2017， 103： 52-58.

[24] SAHEBJAMNIA N， FARD A M F， HAJIAGHAEI-KESHTELI M. Sustainable tire closed-loop supply chain network design： Hybrid metaheuristic algorithms for large-scale networks [J]. Journal of Cleaner Production， 2018， 196： 273-296.

[25] LEE C K H. A review of applications of genetic algorithms in operations management [J]. Engineering Applications of Artificial Intelligence， 2018， 76： 1-12.

[26] 嚴(yán)余松，唐莉，嚴(yán)余偉，等.樞紐小運(yùn)轉(zhuǎn)列車0-1規(guī)劃模型及其遺傳算法[J].系統(tǒng)工程，2000，18（6）：67-70.（YAN Y S， TANG L， YAN Y W， et al. An 0-1 integer programming model on the organization of transship trains in railway terminal and its genetic algorithms[J]. Systems Engineering， 2000， 18（6）： 67-70.）