王亞東， 石全， 尤志鋒， 王芳， 夏偉，2

(1.陸軍工程大學(xué)石家莊校區(qū) 裝備指揮與管理系, 河北 石家莊 050003；2.陸軍步兵學(xué)院石家莊校區(qū) 機(jī)械化步兵系, 河北 石家莊 050083)

0 引言

充足及時(shí)的備件供應(yīng)是維修工作順利開展的前提，也是裝備保障工作的重中之重。備件供應(yīng)優(yōu)化和方案決策需要綜合考慮時(shí)間、風(fēng)險(xiǎn)、成本等多方面因素，優(yōu)化目標(biāo)多元，建模和求解難度較大[1]。因此備件供應(yīng)網(wǎng)絡(luò)優(yōu)化和供應(yīng)方案決策始終是備件保障的重點(diǎn)和難點(diǎn)問題。

目前很多學(xué)者對(duì)備件供應(yīng)優(yōu)化開展了相關(guān)研究。魏國強(qiáng)等對(duì)需求點(diǎn)的優(yōu)先級(jí)聚類排序，研究了資源不足情況下的備件調(diào)度模型，以供應(yīng)開始時(shí)間最早和調(diào)運(yùn)線路最少為目標(biāo)建立模型，并通過賦權(quán)將多目標(biāo)轉(zhuǎn)化為單目標(biāo)優(yōu)化進(jìn)行求解[2]。劉喜春等建立了典型3級(jí)備件供應(yīng)保障結(jié)構(gòu)下的戰(zhàn)時(shí)多階段備件供應(yīng)保障規(guī)劃模型，以期望缺件量為目標(biāo)建立無約束單目標(biāo)優(yōu)化模型，通過迭代方法獲得備件供應(yīng)保障優(yōu)化策略[3]。張斯嘉等建立了基于可信性理論的戰(zhàn)時(shí)多目標(biāo)保障物資供應(yīng)優(yōu)化模型，采用蝙蝠算法求解多目標(biāo)優(yōu)化問題[4]。秦進(jìn)等以最小化系統(tǒng)成本為目標(biāo)，考慮應(yīng)急資源數(shù)量、備用容量和位置數(shù)量等約束，建立了考慮服務(wù)可靠性的應(yīng)急資源布局優(yōu)化模型，提出一種基于矩陣實(shí)數(shù)編碼的遺傳算法求解方法[5]。Fazli-Khalaf等以最小化設(shè)施之間的總供應(yīng)成本和運(yùn)輸時(shí)間，同時(shí)最大化供應(yīng)可靠性為目標(biāo)，建立了應(yīng)急供應(yīng)三目標(biāo)優(yōu)化模型[6]。Zhang等以最大限度地減少運(yùn)輸時(shí)間、運(yùn)輸成本和不滿足度為目標(biāo)，建立了多目標(biāo)三階段隨機(jī)規(guī)劃模型，采用基于隸屬度模糊輔助變量的替代單目標(biāo)模型來處理多目標(biāo)模型[7]。Mohammadi等以最大化總預(yù)期需求覆蓋、最小化總預(yù)期成本并最小化節(jié)點(diǎn)之間滿意率為目標(biāo)，建立了一個(gè)多目標(biāo)隨機(jī)規(guī)劃模型，并采用多目標(biāo)粒子群優(yōu)化算法來解決這個(gè)模型[8]。Su等以完全響應(yīng)時(shí)間和緊急資源成本為目標(biāo)構(gòu)建雙目標(biāo)乘法約束整數(shù)線性規(guī)劃模型，并使用差分進(jìn)化算法搜索模型的最優(yōu)解[9]。綜上所述可以看出，目前相關(guān)研究大多采用多目標(biāo)優(yōu)化模型。這是因?yàn)閱文繕?biāo)優(yōu)化的優(yōu)化目標(biāo)單一，考慮因素較少，不適于備件供應(yīng)優(yōu)化決策。但是，多目標(biāo)優(yōu)化求解難度大且得到的是一系列非支配解集，無法確定唯一的最優(yōu)方案，因此不利于進(jìn)一步?jīng)Q策。

本文為解決備件供應(yīng)優(yōu)化決策問題，以供應(yīng)時(shí)間、風(fēng)險(xiǎn)、成本為目標(biāo)，構(gòu)建備件供應(yīng)多目標(biāo)優(yōu)化模型。在求解算法上，本文受到了帶精英策略的非支配排序遺傳算法(NSGA-II)的啟發(fā)，在多目標(biāo)進(jìn)化算法(MOEA)的基礎(chǔ)上，采用交叉效率數(shù)據(jù)包絡(luò)分析(DEA)對(duì)存檔中的非支配解進(jìn)行排序。一方面，可以指導(dǎo)種群向最優(yōu)效率個(gè)體進(jìn)化；另一方面，可以計(jì)算得到非支配解的自評(píng)效率和互評(píng)效率值，從而對(duì)眾多非支配解進(jìn)行排序、決策出最優(yōu)解。

1 備件供應(yīng)多目標(biāo)優(yōu)化模型

1.1 問題描述

某備件供應(yīng)網(wǎng)絡(luò)由后方倉庫、野戰(zhàn)倉庫和用裝單位3級(jí)節(jié)點(diǎn)組成。用裝單位的維修力量對(duì)損壞裝備進(jìn)行維修時(shí)產(chǎn)生備件需求。后方倉庫為備件的供應(yīng)點(diǎn)，考慮到安全問題，通常距離任務(wù)前線相對(duì)較遠(yuǎn)。因此在后方倉庫和用裝單位之間設(shè)置野戰(zhàn)倉庫，野戰(zhàn)倉庫通常靠前設(shè)置，接收后方倉庫供應(yīng)的備件并前送至用裝單位。備件供應(yīng)優(yōu)化則是根據(jù)備件需求、節(jié)點(diǎn)以及供應(yīng)網(wǎng)路的狀態(tài)信息確定最佳的備件供應(yīng)方案，從而保證以較短的時(shí)間、較低的總風(fēng)險(xiǎn)和較少的成本滿足備件需求。

本文備件供應(yīng)模型建立在以下假設(shè)之上：1)以某類關(guān)鍵備件為研究對(duì)象。2)節(jié)點(diǎn)之間的備件運(yùn)輸成本和運(yùn)輸時(shí)間已知且為定值，節(jié)點(diǎn)之間運(yùn)力充足。3)野戰(zhàn)倉庫的開放費(fèi)用、庫存費(fèi)用、最大儲(chǔ)備量均為已知且為定值。4)備件未能滿足用裝單位需求而造成的后果用缺件損失表示，缺件損失越大的用裝單位表明其優(yōu)先級(jí)越高；用裝單位的備件需求量和缺件損失已知。5)備件全部到達(dá)野戰(zhàn)倉庫后再統(tǒng)一向用裝單位進(jìn)行分配和供應(yīng)。6)供應(yīng)網(wǎng)路的中斷風(fēng)險(xiǎn)僅發(fā)生在野戰(zhàn)倉庫和用裝單位之間。7)暫不考慮越級(jí)供應(yīng)以及同級(jí)之間橫向供應(yīng)的情況。

1.2 變量及其說明

i：后方倉庫的序號(hào)，i=1，2，…，I，I為后方倉庫的數(shù)量；

j：野戰(zhàn)倉庫的序號(hào)，j=1，2，…，J，J為野戰(zhàn)倉庫的數(shù)量；

k：野戰(zhàn)倉庫的序號(hào)，k=1，2，…，K，K為用裝單位的數(shù)量；

Uj：第j個(gè)野戰(zhàn)倉庫最大容量；

dk：第k個(gè)用裝單位的備件需求量；

Rij：備件從后方倉庫i運(yùn)至野戰(zhàn)倉庫j的中斷風(fēng)險(xiǎn)，假設(shè)中斷發(fā)生在途中；

Rjk：備件從野戰(zhàn)倉庫j運(yùn)至用裝單位k的中斷風(fēng)險(xiǎn)，假設(shè)中斷發(fā)生在途中；

xij：后方倉庫i向野戰(zhàn)倉庫j供應(yīng)備件的數(shù)量；

xjk：野戰(zhàn)倉庫j向用裝單位k供應(yīng)備件的數(shù)量；

yj：二進(jìn)制變量，用來標(biāo)注野戰(zhàn)倉庫j的開放情況，yj=1表示開放，yj=0表示關(guān)閉。

1.3 模型建立

目標(biāo)函數(shù)1：備件供應(yīng)的總時(shí)間最短。

(1)

目標(biāo)函數(shù)2：備件供應(yīng)的總風(fēng)險(xiǎn)最低。

(2)

目標(biāo)函數(shù)3：備件供應(yīng)總成本最小。

minC=Co+Ct+Ci+Cs,

(3)

式中：Co表示野戰(zhàn)倉庫的開放成本；Ct表示運(yùn)輸成本；Ci表示備件在野戰(zhàn)倉庫的庫存成本；Cs表示缺件損失。各項(xiàng)成本計(jì)算如下：

(4)

(5)

(6)

(7)

滿足以下約束：

(8)

(9)

(10)

(11)

(12)

xij∈N+,xjk∈N+,yj={0,1}.

(13)

式中：N+為正整數(shù)。約束(8)式、(9)式表示未開放的野戰(zhàn)倉庫不參與備件供應(yīng)，同時(shí)還規(guī)定了供入的備件數(shù)量不能超過野戰(zhàn)倉庫的最大容量；約束(10)式規(guī)定了用裝單位的備件需求應(yīng)得到滿足；約束(11)式規(guī)定了野戰(zhàn)倉庫備件的供出量應(yīng)不超過供入量；約束(12)式規(guī)定了對(duì)于每個(gè)用裝單位，備件供應(yīng)的延遲時(shí)間應(yīng)不超過允許的最大截止時(shí)間。這里規(guī)定的延遲時(shí)間與目標(biāo)函數(shù)中的總供應(yīng)時(shí)間有所區(qū)別。供應(yīng)總時(shí)間用來體現(xiàn)調(diào)運(yùn)車輛數(shù)量和總里程等信息，而延遲時(shí)間則表示備件從供應(yīng)開始到全部到達(dá)所消耗的時(shí)間；約束(13)式規(guī)定了決策變量的類型和取值范圍。

2 二次目標(biāo)交叉效率排序多目標(biāo)進(jìn)化算法

為了求解備件供應(yīng)多目標(biāo)優(yōu)化模型的最優(yōu)解以及得到最終方案，本文提出了一種二次目標(biāo)交叉排序多目標(biāo)進(jìn)化算法(SGCES-MOEA)。該算法采用進(jìn)化計(jì)算的基本框架，在每次迭代中將種群中的非支配解放入存檔，并從存檔中選擇精英個(gè)體指導(dǎo)種群的進(jìn)化。與傳統(tǒng)進(jìn)化計(jì)算不同的是，SGCES-MOEA采用了類似于NSGA-Ⅱ的排序機(jī)制來選擇精英個(gè)體。NSGA-Ⅱ先根據(jù)種群中個(gè)體被支配的次數(shù)對(duì)個(gè)體分層，再對(duì)同一層的個(gè)體計(jì)算擁擠度距離，并根據(jù)擁擠度距離對(duì)每一層個(gè)體進(jìn)行排序。第一層為非支配解(被支配次數(shù)為0)，該層擁擠度距離最大的個(gè)體為精英個(gè)體[10]。

SGCES-MOEA從NSGA-Ⅱ得到啟發(fā)，用非支配個(gè)體的效率值代替擁擠度距離進(jìn)行排序，從而選擇出精英個(gè)體。相對(duì)于擁擠度距離，效率可以反映出個(gè)體在優(yōu)化模型中的綜合指標(biāo)，具有實(shí)際意義。個(gè)體的效率值計(jì)算則采用改進(jìn)的DEA法。每個(gè)非支配個(gè)體相當(dāng)于一個(gè)決策單元(DMU)，根據(jù)優(yōu)化模型計(jì)算得到各指標(biāo)值后，再利用DEA評(píng)估模型計(jì)算個(gè)體的效率值。SGCES-MOEA的基本步驟如圖1所示。

圖1 SGCES-MOEA的基本流程Fig.1 Basic flow chart of SGCES-MOEA

2.1 進(jìn)化策略

進(jìn)化策略是對(duì)父代種群個(gè)體進(jìn)行的變異、交叉、選擇等操作，從而產(chǎn)生子代種群的過程。通過進(jìn)化，種群不斷得到改善，從而獲得較優(yōu)解。

步驟1變異。令ps=ps1,ps2,…,…,psh,…,psD為第s個(gè)父代個(gè)體向量，psD為第s個(gè)個(gè)體向量ps上的第D個(gè)變量；vs=vs1,vs2,…,…,vsh,…,vsD為根據(jù)變異策略產(chǎn)生的變異個(gè)體向量，vsD為第s個(gè)變異個(gè)體向量vs上的第D個(gè)變量。本文采用DE/current-to-best/1策略進(jìn)行變異：

vs=ps+F×(pe-ps)+F×(pr-ps)，

(14)

式中：pe為當(dāng)前種群中的最優(yōu)個(gè)體；pr為從存檔中隨機(jī)選取的個(gè)體，且ps≠pe≠pr；F為變異率。

(15)

式中：rsh為[0,1]之間均勻分布的隨機(jī)數(shù)；CR為交叉率，CR∈[0,1]，即當(dāng)變異個(gè)體vs上第h個(gè)變量vsh對(duì)應(yīng)的隨機(jī)數(shù)小于交叉率時(shí)，交叉?zhèn)€體向量us上對(duì)應(yīng)變量ush取父代個(gè)體ps對(duì)應(yīng)變量xsh，否則保留變異個(gè)體上的變量值。

步驟3選擇。根據(jù)種群中的父代個(gè)體和交叉?zhèn)€體的適應(yīng)度函數(shù)，一一比較其支配關(guān)系。若某父代個(gè)體被對(duì)應(yīng)的交叉?zhèn)€體支配，則用該交叉?zhèn)€體替換對(duì)應(yīng)父代個(gè)體。最終得到的個(gè)體構(gòu)成選擇個(gè)體種群，作為子代個(gè)體進(jìn)入下一次迭代。

步驟4更新存檔。將子代種群中的非支配解放入存檔中，與存檔個(gè)體混合后，再次比較支配關(guān)系，選出混合種群的非支配解作為新的存檔。

2.2 適應(yīng)度函數(shù)及編碼

由于本文構(gòu)建的多目標(biāo)優(yōu)化模型屬于帶約束的優(yōu)化問題，本文采用動(dòng)態(tài)罰函數(shù)法處理模型約束[11]：

(16)

式中：x為決策變量向量；f(x)為適應(yīng)度函數(shù)；O(x)為目標(biāo)函數(shù)；M0為初始懲罰因子；ni為當(dāng)前算法迭代次數(shù)，maxni為算法最大迭代次數(shù)；MK為懲罰系數(shù)；P(x)為約束違反度，

(17)

m為約束條件的總個(gè)數(shù)，ch(x)為約束違反度。若不等式約束表示為lh(x)≤0的形式，其約束違反度ch(x)=max(0,lh(x))；若等式約束表示為hh(x)=0的形式，其約束違反度ch(x)=max (0,|hh(x)-ε|)，ε為一個(gè)很小的實(shí)數(shù)。

如圖1所示，取直徑為170 mm，直徑偏差范圍為0～0.8 mm的帶輪進(jìn)行實(shí)驗(yàn)。對(duì)帶輪邊緣輪廓進(jìn)行測(cè)量，通過雙目視覺圓弧輪廓特征提取和匹配之后，得出邊緣外輪廓點(diǎn)的空間三維坐標(biāo)。圖2所示為求得的帶輪邊緣外輪廓點(diǎn)云圖及空間圓弧擬合圖，共獲得266個(gè)邊緣外輪廓點(diǎn)的三維坐標(biāo)。利用本文擬合優(yōu)化方法，對(duì)上述點(diǎn)云進(jìn)行8次擬合，平均迭代數(shù)為6 895次。擬合結(jié)果如表1所示。

2.3 指標(biāo)計(jì)算

根據(jù)DEA模型指標(biāo)的性質(zhì)，將評(píng)價(jià)指標(biāo)分為輸入型指標(biāo)和輸出型指標(biāo)，輸入型指標(biāo)的值越小、對(duì)應(yīng)DMU的效率越高，輸出型指標(biāo)的值越大、DMU的效率越高。

本文的輸入型指標(biāo)包括：

1)供應(yīng)成本：計(jì)算公式即目標(biāo)函數(shù)(3)式；

2)供應(yīng)時(shí)間：計(jì)算公式即目標(biāo)函數(shù)(1)式；

3)約束違反度：由于模型還涉及容量以及流量平衡等約束，將優(yōu)化結(jié)果的約束違反度作為指標(biāo)之一。模型的總體約束違反度計(jì)算公式為(17)式。

輸出型指標(biāo)包括：

圖2 個(gè)體編碼方式示意Fig.2 Chromosome individual coding pattern

1)可靠性指標(biāo)：用備件供應(yīng)的風(fēng)險(xiǎn)衡量，其值為總風(fēng)險(xiǎn)的倒數(shù)。風(fēng)險(xiǎn)計(jì)算公式為目標(biāo)函數(shù)(2)式；

2)時(shí)效性指標(biāo)：用延遲時(shí)間即供應(yīng)開始到最后一波備件到達(dá)的時(shí)間衡量，取延遲時(shí)間的倒數(shù)。延遲時(shí)間的計(jì)算公式為約束(12)式；

3)滿足度：用各用裝單位的備件滿足度表示，計(jì)算公式即約束(10)式。

2.4 二次目標(biāo)交叉效率排序

算法采用DEA模型計(jì)算個(gè)體效率。傳統(tǒng)DEA由Cooper等提出，其構(gòu)建的固定規(guī)模收益的CCR模型[12]如下：

(18)

(19)

μra≥ε,υba≥ε,r=1,2,…,so,b=1,2,…,mi.

(20)

式中：Edd表示DMUd的自評(píng)效率，DMUd為第d個(gè)DMU；r為輸出指標(biāo)序號(hào)；so為輸出指標(biāo)的總個(gè)數(shù)；Yrd表示DMUd的第r個(gè)輸出指標(biāo)；mi為輸入指標(biāo)的總個(gè)數(shù)；μrd為Yrd的權(quán)重；Xbd表示DMUd的第b個(gè)輸入指標(biāo)；υbd為Xbd的權(quán)重；Ea為第a個(gè)決策單元DMUa的自評(píng)效率；Yra表示第a個(gè)決策單元DMUa的第r個(gè)輸出指標(biāo)；μra為Yra的權(quán)重；n為決策單元的個(gè)數(shù)；Xba為第a個(gè)DMUb的第b個(gè)輸入指標(biāo);νba為Xba的權(quán)重；ε為比任何正數(shù)都小的非阿基米德數(shù)。約束(19)式規(guī)定了所有DMU的效率值應(yīng)介于0～1之間。約束(20)式規(guī)定了所有權(quán)重應(yīng)大于0.

傳統(tǒng)DEA模型還存在兩個(gè)缺點(diǎn)：1)由于其采用的是自評(píng)效率，只能區(qū)別有效單元(效率值等于1)和無效單元(效率值小于1)，無法對(duì)有效單元進(jìn)行排序；2)計(jì)算得到的輸入和輸出權(quán)重不唯一[13]。針對(duì)第1個(gè)問題，SGCES-MOEA通過計(jì)算各DMU的交叉效率，可以對(duì)所有DMU進(jìn)行評(píng)估和排序。針對(duì)第2個(gè)問題，SGCES-MOEA通過構(gòu)建二次目標(biāo)模型，使得權(quán)重向量取值唯一。

交叉效率計(jì)算方法如下。

(21)

根據(jù)CCR模型對(duì)每個(gè)DMU分別求其自評(píng)效率，再根據(jù)(21)式求每個(gè)DMU的相互之間的互評(píng)效率，得到互評(píng)效率矩陣：

則DMUa的交叉效率為

(22)

由于CCR中求得的最優(yōu)權(quán)重組合可能不唯一，在利用(21)式計(jì)算交叉效率時(shí)其值也不唯一。為解決此類問題，Doyle等提出了通過引入二次目標(biāo)的方法，在多重最優(yōu)解中選擇一組最優(yōu)解，以最終確定每個(gè)DMU的交叉效率[14]。常用的二次目標(biāo)方法有激進(jìn)型策略和慈善型策略[15]。本文的二次目標(biāo)交叉效率DEA模型如下：

(23)

(24)

(25)

(26)

(27)

(28)

(29)

(30)

μra≥ε，r=1,2,…,so；

(31)

υba≥ε，b=1,2,…,mi.

(32)

3 算例及其分析

3.1 算例描述

以某次演習(xí)任務(wù)的維修備件的支援保障為例，共有6個(gè)用裝單位在前線執(zhí)行演習(xí)任務(wù)。根據(jù)裝備的損壞情況，通過維修任務(wù)確定需要的備件種類和數(shù)量。為檢驗(yàn)支援保障能力，部分關(guān)重件未隨部隊(duì)攜行，僅由后方供應(yīng)。維修備件的儲(chǔ)備由兩個(gè)后方倉庫承擔(dān)，同時(shí)在后方倉庫和用裝單位之間有4個(gè)備選可開設(shè)野戰(zhàn)倉庫。備件由后方倉庫運(yùn)送至野戰(zhàn)倉庫，再根據(jù)前線需求由野戰(zhàn)倉庫進(jìn)行分配轉(zhuǎn)運(yùn)至各個(gè)用裝單位。后方倉庫、野戰(zhàn)倉庫以及用裝單位的相關(guān)數(shù)據(jù)見表1～表5.

表1 單位備件在各節(jié)點(diǎn)間的運(yùn)輸時(shí)間Tab.1 Transportation times of spare parts between nodes h

表2 單位備件在各節(jié)點(diǎn)間的運(yùn)輸成本Tab.2 Transportation costs of spare parts between nodes 元

表3 野戰(zhàn)倉庫和用裝單位間路線中斷風(fēng)險(xiǎn)Tab.3 Disruption risk of transit routes from fieldwarehouses to demand units

3.2 多目標(biāo)優(yōu)化結(jié)果

利用SGCES-MOEA對(duì)算例模型進(jìn)行求解，本文采用MATLAB 2014b軟件進(jìn)行編程，平臺(tái)為個(gè)人筆記本(Intel Core i5-6300HQ (2.3 GHz/L3) CPU和4.00 GB RAM)。相關(guān)參數(shù)設(shè)置如下：種群規(guī)模為100，外部存檔規(guī)模為100，最大迭代次數(shù)為500，變異因子為0.9，交叉率為0.8，選擇閾值為0.8，初始懲罰因子為10 000 000，懲罰系數(shù)為10 000.

表4 野戰(zhàn)倉庫的庫存容量、單位備件庫存成本、開設(shè)費(fèi)用Tab.4 Inventory capacity, unit spare parts inventorycost, and running cost of field warehouse

通過計(jì)算，算法共求得13個(gè)非支配解，每個(gè)解代表一個(gè)供應(yīng)方案。表6給出了所有方案的輸入指標(biāo)、輸出指標(biāo)以及交叉效率和自評(píng)效率的值。其中滿足度1～6為6個(gè)用裝單位的備件滿足度。由表6可以看出，13個(gè)方案中有10個(gè)方案為有效方案，其自評(píng)效率均為1，僅通過自評(píng)效率無法進(jìn)一步從這10個(gè)方案中選出最優(yōu)方案。通過計(jì)算和比較各方案的交叉效率，可以選出最優(yōu)方案，其交叉效率值最大為0.927 8.

表5 用裝單位的備件需求、缺件損失、最大延遲時(shí)間Tab.5 Spare parts requirements, missing loss, maximumdelay time of demand units

表6 SGCES-MOEA求解結(jié)果Tab.6 Solved results of SGCES-MOEA

3.3 靈敏度分析

通過算例的計(jì)算結(jié)果可以看出，通過SGCES-MOEA可以求得模型的非支配解并根據(jù)效率進(jìn)行排序，從而選出最優(yōu)解。為了進(jìn)一步驗(yàn)證本文的排序策略對(duì)算法和結(jié)果的影響。分別用未進(jìn)行排序的MOEA和采用自評(píng)效率排序的多目標(biāo)進(jìn)化算法(SES-MOEA)對(duì)模型進(jìn)行優(yōu)化，將其結(jié)果與SGCES-MOEA的結(jié)果進(jìn)行比較。

表7為MOEA求得的16個(gè)非支配解的相關(guān)指標(biāo)，由于未采用排序機(jī)制而沒有計(jì)算自評(píng)效率和交叉效率值。16個(gè)DMU僅存在非支配關(guān)系，無法進(jìn)一步?jīng)Q策出最優(yōu)的供應(yīng)方案。表8為SES-MOEA求得的7個(gè)非支配解的相關(guān)指標(biāo)及其自評(píng)效率值。由表8可以看出SES-MOEA結(jié)果中有4個(gè)DMU的效率值等于1，為有效DMU，剩余3個(gè)為非有效DMU. 此時(shí)僅通過自評(píng)效率同樣無法進(jìn)一步選擇最優(yōu)供應(yīng)方案。

表7 MOEA求解結(jié)果Tab.7 Solved results of MOEA

表8 SES-MOEA求解結(jié)果Tab.8 Solved results of SES-MOEA

為了驗(yàn)證3種算法求得方案的優(yōu)劣，將SGCES-MOEA、MOEA和SES-MOEA求得的所有DMU混合后作為整體，進(jìn)一步使用本文的二次目標(biāo)交叉效率DEA法計(jì)算所有DMU的自評(píng)效率和交叉效率，并根據(jù)交叉效率進(jìn)行排序。所有DMU的效率值如圖3所示。從圖3中可以看出，SGCES-MOEA對(duì)應(yīng)的交叉效率值普遍較大，而MOEA對(duì)應(yīng)的交叉效率值普遍較小，SES-MOEA介于二者之間。通過計(jì)算發(fā)現(xiàn)：SGCES-MOEA對(duì)應(yīng)DMU的平均交叉效率值為0.878 5；MOEA對(duì)應(yīng)DMU的平均交叉效率值為0.723 8；SES-MOEA對(duì)應(yīng)DMU的平均交叉效率值為0.752 7. 由此可見，SGCES-MOEA的效果要優(yōu)于其他兩種算法。

圖3 3種算法結(jié)果的效率值對(duì)比Fig.3 Comparison of efficiency values of solved results of three algorithms

4 結(jié)論

本文主要解決備件供應(yīng)優(yōu)化決策問題，同時(shí)考慮時(shí)間、風(fēng)險(xiǎn)和成本3個(gè)目標(biāo)，建立了備件供應(yīng)網(wǎng)絡(luò)優(yōu)化的多目標(biāo)優(yōu)化模型。使用基于交叉效率排序的MOEA求得模型的非支配解集，并決策出最優(yōu)供應(yīng)方案。研究?jī)?nèi)容的主要貢獻(xiàn)有：1)提出的多目標(biāo)優(yōu)化算法和約束處理方法可以解決備件供應(yīng)的多目標(biāo)約束優(yōu)化問題；2)提出的基于效率排序策略可以進(jìn)一步對(duì)非支配解排序，同時(shí)可以指導(dǎo)算法向效率最大的個(gè)體收斂；3)通過計(jì)算交叉效率克服了僅憑自評(píng)效率無法對(duì)有效單元進(jìn)行排序的問題；4)提出的二次目標(biāo)DEA模型解決了評(píng)估權(quán)重不唯一的問題。綜上，本文的模型和算法為解決備件供應(yīng)優(yōu)化以及最優(yōu)備件供應(yīng)方案決策提供了框架和求解思路。