基于隨機p-魯棒優(yōu)化的國家血液戰(zhàn)略儲備庫選址-庫存模型

周愉峰，李 志，劉思峰

(1.重慶工商大學重慶市發(fā)展信息管理工程技術(shù)研究中心，重慶 400067; 2.南京航空航天大學管理科學與工程博士后流動站，江蘇 南京 210016)

1 引言

大地震、特大恐怖襲擊等嚴重突發(fā)事件可能導(dǎo)致大量人員傷亡，短時間內(nèi)造成血液需求驟增。血液屬于易腐物品，其采集、檢測、制備、儲存、使用等要求嚴格、無法即采即用，且難以被其他物品替代。再加上血液供應(yīng)關(guān)系到傷員的生命健康安全。因而完善血液儲備體系對于災(zāi)后搶救、降低社會恐慌，提高血液保障水平等方面具有重要意義[1]。

目前，我國尚無統(tǒng)一的國家血液戰(zhàn)略儲備庫體系。而以往的救災(zāi)實踐又對國家層面的血液戰(zhàn)略儲備庫建設(shè)提出了迫切要求。因此，國家血液戰(zhàn)略儲備庫的建設(shè)問題在理論與實踐中亟待解決。當前，一些研究者初步設(shè)想在若干省級血液中心建立國家血液戰(zhàn)略儲備庫[2-3]。例如，吳衛(wèi)星等[4]建議在人口密度較大的區(qū)域設(shè)立1到3個由國家衛(wèi)生部控制的國家血液戰(zhàn)略儲備單元。鄭忠偉等[1]提出在我國6大區(qū)域建立若干個國家層級的應(yīng)急血液儲備庫，并執(zhí)行血液定期輪換策略。黎城等[5]提出以省為單位儲存紅細胞，構(gòu)成國家血液儲備單元。沈紅艷[3]提出以全國七大軍區(qū)駐地血液中心為基礎(chǔ)，建設(shè)國家血液儲備庫。

上述文獻均從定性角度，探討國家血液戰(zhàn)略儲備庫建設(shè)的可行性和必要性等問題，缺乏定量研究。鑒于此，孟超[2]考慮調(diào)配時效和可靠性目標，研究了國家血液戰(zhàn)略儲備設(shè)施選址-分配問題。但孟超[2]沒有考慮庫存因素。由于大規(guī)模涉血災(zāi)害的發(fā)生概率非常小，若一味追求災(zāi)時血液保障時效最優(yōu)而完全忽略成本因素，將產(chǎn)生極其高昂的日常運營成本。因此，國家血液戰(zhàn)略儲備庫選址應(yīng)該考慮日常庫存運營成本，建立預(yù)算約束條件下的應(yīng)急調(diào)配時效最優(yōu)化的選址-庫存問題模型(Location-inventory problem，LIP)。此外，孟超[2]假設(shè)國家血液戰(zhàn)略儲備庫在現(xiàn)有省級血液中心基礎(chǔ)上設(shè)置，因而可充分利用既有資源而無需考慮設(shè)施建設(shè)成本。該假設(shè)顯然與現(xiàn)實不符。原因在于各省血液中心的庫存容量是基于日常需求設(shè)計的。設(shè)為國家層面的戰(zhàn)略儲備庫后，災(zāi)時爆炸性增長的血液需求會導(dǎo)致血液中心儲備能力受到限制。因此，須對選定的血液中心進行改擴建處理，并在模型中考慮改擴建所涉及的固定投入成本及庫存運營成本增量。因此，本文集成庫存決策，考慮預(yù)算約束與系統(tǒng)及時度目標，在孟超[2]的基礎(chǔ)上進一步研究國家血液戰(zhàn)略儲備庫LIP模型。

與本研究相關(guān)的問題主要有兩類。第一類為應(yīng)急設(shè)施選址問題(Facility location problem，F(xiàn)LP)及采供血機構(gòu)選址問題。有關(guān)應(yīng)急設(shè)施選址問題的研究較多，例如，張敏和張玲[6]研究了基于失效情景的應(yīng)急設(shè)施選址模型。近年來，以采供血機構(gòu)為背景的FLP研究也引起了研究者們的關(guān)注。Zahiri等[7-8]分別從靜態(tài)與多階段優(yōu)化的角度研究了固定采血點與流動獻血車的選址-指派問題。Sha Yue和Huang Jun[9]建立了震后應(yīng)急血液FLP模型。周愉峰等[10]以血站為例研究了應(yīng)急物資儲備庫選址的可靠性p-中值模型。Chaiwuttisak等[11]建立了兩類采血服務(wù)設(shè)施的選址以及與醫(yī)院、血站之間的指派優(yōu)化模型。Zahiri和Pishvae[12]考慮血型兼容與需求不確定性，建立了血液供應(yīng)鏈網(wǎng)絡(luò)設(shè)計模型。但這些研究一般強調(diào)受災(zāi)點滿意度或?qū)π枨蟮捻憫?yīng)速度等因素，很少同時考慮設(shè)施的日常運營成本與災(zāi)時響應(yīng)速度，且沒有集成庫存問題。第二類為選址-庫存問題(LIP)。以往的LIP文獻大都研究的是一般物資在常規(guī)條件下的集成決策問題。例如，Sadjadi等[13]研究了三級供應(yīng)鏈網(wǎng)絡(luò)設(shè)計中的隨機需求LIP。Zhang Zhihai和Unnikrishnan[14]研究了閉環(huán)供應(yīng)鏈中的隨機需求LIP。Li Yanhui等[15]研究了閉環(huán)物流系統(tǒng)中帶第三方物流商的LIP。Vahdani等[16]研究了三級供應(yīng)鏈網(wǎng)絡(luò)設(shè)計中的LIP。Diabat等[17]研究了需求與提前期均不確定的LIP。這些研究沒有考慮血液制品特性與血液保障的應(yīng)急特性。有關(guān)應(yīng)急設(shè)施LIP的研究較少[2]，但一些文獻基于隨機規(guī)劃[18-19]、魯棒優(yōu)化[20]、可靠性優(yōu)化[21]等方法研究了LIP，可為應(yīng)急條件下的LIP研究提供借鑒。

由于血液產(chǎn)品特性以及應(yīng)急血液保障特性的存在，使得國家血液戰(zhàn)略儲備庫選址-庫存問題需要考慮眾多因素，例如：儲備庫日常運行成本與災(zāi)時響應(yīng)速度的均衡、血液制品的多品種性、災(zāi)后應(yīng)急血液需求的多階段性、日常需求的隨機性、儲備庫建立后的設(shè)施容量限制、改擴建后的庫存狀態(tài)變更、應(yīng)急條件下的協(xié)同定位等因素。上述因素在以往的應(yīng)急FLP與LIP研究中沒有同時考慮，因而需在國家血液戰(zhàn)略儲備庫LIP問題中予以擴展研究。同時，考慮到多情景涉血災(zāi)害發(fā)生的不確定性，本文引入p-魯棒,構(gòu)建了一個考慮多情景、多血型、災(zāi)后需求多階段、日常需求隨機、設(shè)施有容量限制及協(xié)同定位、有預(yù)算約束的國家血液戰(zhàn)略儲備庫LIP模型。針對模型性質(zhì)，設(shè)計了求解問題的兩類亞啟發(fā)式算法：遺傳算法(Genetic algorithm，GA)與免疫優(yōu)化算法(Immune optimization algorithm，IOA)，并通過仿真算例驗證了所提模型及算法的有效性。

2 國家血液戰(zhàn)略儲備庫選址-庫存問題的隨機p-魯棒優(yōu)化模型

2.1 問題描述

設(shè)置國家血液戰(zhàn)略儲備庫的目的是在災(zāi)時為傷員提供快速有效的血液保障。因而以災(zāi)時的系統(tǒng)響應(yīng)時效最快來體現(xiàn)這一目標。由于涉血災(zāi)害發(fā)生的突發(fā)性和不可預(yù)見性，應(yīng)急血液需求難以事先確定。為此，采用情景分析方法來描述應(yīng)急需求。假設(shè)所有情景均只有一個需求點，即任何情景下只有一個省份發(fā)生大規(guī)模涉血災(zāi)害。引發(fā)大量應(yīng)急血液需求的突發(fā)事件主要有大規(guī)模地震與特大恐怖襲擊。在這些突發(fā)事件中，應(yīng)急血液需求具有明顯的階段性，可分為緊急救援期和救援中后期兩個階段。災(zāi)害發(fā)生后的72小時為緊急救援期，即兩個階段的間隔時間為3天。災(zāi)時的應(yīng)急血液需求與救援期所處階段密切相關(guān)。應(yīng)急條件下由于需求量驟增，儲備庫容量受到限制，每種血液制品在出庫后需要進行庫存補充。由于災(zāi)時采集儲備能力受限與需求增長，為了保障受災(zāi)點的血液供應(yīng)，允許協(xié)同定位，即允許一個需求點同時由多個儲備庫提供應(yīng)急血液需求。

由于大規(guī)模涉血災(zāi)害發(fā)生的概率很小，因而不能無限追求應(yīng)急系統(tǒng)響應(yīng)速度，必須考慮設(shè)施在常規(guī)條件下的日常運營成本，即考慮預(yù)算限制。為了最大限度的利用既有資源、降低成本，選擇若干省級血液中心進行改擴建，而非在新址上全新建設(shè)。因此，須考慮改擴建的固定投入成本，以及庫存狀態(tài)變化引起的日常運營成本增量。省級血液中心被選建為國家血液戰(zhàn)略儲備庫后，需大量增加各種血液制品的儲備量以應(yīng)對將來可能發(fā)生的突發(fā)事件。由于血液只能通過人體采集，得來不易，且保管周期有限。為了避免血液制品的大量過期報廢，在平時需對庫存進行輪換更新，因而省級血液中心被選建為戰(zhàn)略儲備庫后，其庫存系統(tǒng)會發(fā)生變化。常規(guī)條件下，省級血液中心日常需求量隨機，可根據(jù)各地的歷史數(shù)據(jù)擬合成一定的概率分布。

綜上，問題可以被描述為一個考慮多情景、多階段不確定應(yīng)急需求、隨機日常需求、多血型、有容量限制、有預(yù)算約束及協(xié)同定位的國家血液戰(zhàn)略儲備LIP。要解決的問題是：應(yīng)該建立多少個儲備庫？建在哪里？如何對應(yīng)急血液需求點與儲備庫進行指派？如何分配應(yīng)急血液需求？如何進行日常庫存決策？

2.2 符號說明

2.2.1 參數(shù)

I為需求點集合，即每個省級行政區(qū)；i∈I；

J為候選地集合，即現(xiàn)有的省級血液中心j∈J；

N為應(yīng)急血液需求階段集合；n∈N，n=1,2；

K為血型集合；k∈K；

tsn,n+1為相鄰應(yīng)急需求階段的時間間隔；

Capjk為候選地j對k型血的儲備能力上限；

β為庫存補充系數(shù)；

L為提前期；

FB為建立國家血液儲備庫的預(yù)算；

hjk為儲備庫j中k型血的單位庫存持有成本；

c為單位運輸成本；

dij為需求點i與候選地j之間的距離；

Ti為需求點i的應(yīng)急需求調(diào)配時限；

tij為需求點i與候選點j之間的應(yīng)急血液調(diào)配時間，tij由運輸時間，需求點的接收耗時tp和候選點調(diào)劑準備耗時tb構(gòu)成，即tij=dij/Vij+tp+tb；

χ為規(guī)劃期長度；

θ為與庫存成本相關(guān)的權(quán)因子；

λ為與輪換更新成本相關(guān)的權(quán)因子；

fj為血液中心j改擴建為國家血液戰(zhàn)略儲備庫的固定投入成本；

sjk為儲備庫j對k型血的安全庫存保有量；

M為一個大的正數(shù)。

2.2.2 變量

Xj表示在j地建立儲備庫時為1，否則為0

Sjk為儲備庫j中k型血的庫存定至點

Qjk1為候選地j未建儲備庫前對k型血的采血批量

Qjk2為候選地j被建為儲備庫后對k型血的采血批量

2.3 庫存成本分析

為了應(yīng)對未來可能發(fā)生的應(yīng)急需求，選擇對現(xiàn)有的若干省級血液中心進行改擴建，并固化增持一部分安全庫存sjk。此時，原有的一些庫存參數(shù)，包括儲庫存容量、安全庫存，血液采集批量等發(fā)生了變化。因而需考慮建設(shè)國家儲備庫所支付的變動成本，包括：改擴建的固定投入成本，由改擴建引發(fā)的新增庫存持有成本，新增安全庫存持有成本及新增血液輪換更新成本。

(1)選作戰(zhàn)略儲備庫之前的運營成本分析：

戰(zhàn)略儲備庫的日常庫存控制參數(shù)可表示如下：

安全庫存：

(1)

式(1)中，zjk為儲備庫j中k型血的血液保障水平系數(shù)。

訂貨點：

(2)

平均庫存：

(3)

j地未設(shè)為戰(zhàn)略儲備庫時的日常周轉(zhuǎn)庫存成本可以表示為：

(4)

式(4)，TPjk為j地兩次連續(xù)訂購的時間間隔：

TPjk=Qjk1/μjk,?j∈J,k∈K

(5)

假設(shè)血液本地輪換成本為v(Qjk1)，令：

v(Qjk1)=gj+ajQjk1,?j∈J,k∈K

(6)

式(6)中g(shù)j表示血液固定輪換更新成本；aj表示血液可變輪換更新成本。

對式(4)除以TPjk可以得到成本率的表達式：

(7)

對式(7)中Qjk1求導(dǎo)，等于0，得到：

(8)

則j地未建戰(zhàn)備庫前的運營成本可以表示為：

(9)

(2)選作戰(zhàn)略儲備庫之后的運營成本分析:

通過設(shè)置安全庫存sjk以保障災(zāi)時的應(yīng)急需求，為此，相應(yīng)的庫存控制參數(shù)需重新定義：

安全庫存：

(10)

訂貨點：

(11)

庫存定至點：

(12)

平均庫存：

(13)

j地設(shè)為戰(zhàn)略儲備庫后的周轉(zhuǎn)庫存成本為：

(14)

式(14)除以TPjk2得到成本率的表達式：

+λμjkv(Qjk2)/Qjk2,?j∈J,k∈K

(15)

建為儲備庫后，災(zāi)時應(yīng)急需求劇增，血庫容量受限。此時，不能再用求導(dǎo)的方法計算最優(yōu)血液采集量。為此，設(shè)置兩個約束來處理應(yīng)急血液采集量Qjk2。

約束(16)用來處理最大采集量。 假設(shè)儲備庫j對k型血的最大可采集量用Qjk2max表示，則

Qjk2≤Qjk2max,?j∈J,k∈K

(16)

約束(17)為候選儲備庫的服務(wù)水平約束。假設(shè)δ為違背庫存容量約束的最大概率，SDjk(LTj)為提前期內(nèi)的隨機需求量。采用機會約束規(guī)劃處理服務(wù)水平的滿足率：

(17)

在此隨機約束中，允許達到100%的服務(wù)水平。

(18)

同時，為了保證安全庫存保有量sjk也滿足儲備庫已開放的前提，將成本率的表達式(15)改寫為：

+λμjkv(Qjk2)/Qjk2,?j∈J,k∈K

(19)

因此，在j地由于建設(shè)戰(zhàn)略儲備庫所引起的新增成本為：

(20)

2.4 模型構(gòu)建

至此，可建立LIP模型MP1：

(21)

(22)

(23)

(24)

(25)

(26)

(27)

(28)

(29)

(30)

Qjk2≤Qjk2max,?j∈J,k∈K

(31)

(32)

Xj∈{0,1},?j∈J

(33)

目標(21)表示災(zāi)時的加權(quán)需求時間最短，用來表示設(shè)施網(wǎng)絡(luò)的系統(tǒng)及時度。式(22)為預(yù)算約束。約束(23)表示任意儲備庫的任意血液制品在任意階段的供應(yīng)量不能超過可用的最大庫存量。約束(24)表示候選地在任意需求階段開始時的最大可用庫存量不超過安全庫存保有量與儲備能力上限。式(25)為每階段開始時最大可用庫存量的表達式。式(26)為災(zāi)害情景發(fā)生時最大可用庫存量表達式；式(27)表示應(yīng)急血液需求不出現(xiàn)短缺。約束(28)為應(yīng)急調(diào)配時間限制。約束(29) 表示血液調(diào)配只能在已經(jīng)設(shè)立的儲備庫進行。式(30)為非負約束。式(31)為血液采集量上限約束。約束(32)表示應(yīng)急血液儲備必須在開放儲備庫進行，同時必須滿足血液保障水平。約束(33)為0-1整數(shù)變量約束。

確定情景下的隨機優(yōu)化可產(chǎn)生目標值最優(yōu)的解, 但是在特定的情景中可能表現(xiàn)很差[22]。為了規(guī)避這種情況，建立國家血液戰(zhàn)略儲備庫LIP的隨機p-魯棒優(yōu)化模型。

由模型MP1可知，確定情景ω下，國家血液戰(zhàn)略儲備庫選址-庫存問題的模型MP2為

MP2：

(34)

(35)

(36)

(37)

(38)

(39)

(40)

(41)

(42)

MP3：

(43)

(44)

(45)

式(44)為p-魯棒優(yōu)化的約束條件。式(45)表示任意需求點所有需求情景的概率之和為1。約束(22) 、約束(31) -約束(33)與模型MP1、MP2相同，約束(35)-(42)與模型MP2相同。

故性質(zhì)1得證。

性質(zhì)2：令I(lǐng)i={j|?ij=1,?j}表示在時限約束內(nèi)可對需求點i提供服務(wù)的儲備庫集合。按應(yīng)急血液調(diào)配時間tij,j∈Ii升序排列Ii內(nèi)的所有儲備庫，得到{j1,j2,…jr|j∈Ii}。若儲備庫已選定，則最優(yōu)的分配方案為按照j1,j2,…jr的順序依次對需求點i提供血液保障。

證明：由于任意情景均不會有兩個以上省份同時發(fā)生大規(guī)模涉血災(zāi)害，所以若開放儲備庫給定，則系統(tǒng)總成本也固定。此時，將需求點i指派給最近的持有救災(zāi)庫存的儲備庫能提高供應(yīng)速度。故性質(zhì)2得證。

3 求解算法

MP3也為NP-hard問題，難以精確求解。為此，設(shè)計一種0-1編碼的遺傳算法對模型進行求解。遺傳算法具有并行計算能力強，收斂速度快，不易陷入局部最優(yōu)等特點，對于解決復(fù)雜組合優(yōu)化問題具有突出優(yōu)勢。同時，為了驗證所提遺傳算法的性能，另外設(shè)計了一種免疫優(yōu)化算法，并在后文對這兩種經(jīng)典進化算法進行比較分析。

所設(shè)計的遺傳算法流程描述如下：

①編碼/解碼

(46)

否則：

(47)

②初始化種群

③適應(yīng)度評價

由于種群中可能存在非可行解，需對不滿足約束的染色體進行懲罰。

設(shè)lost為某條非可行解染色體r的應(yīng)急血液需求短缺量，則：

(48)

然后對需求短缺、預(yù)算超支、p-魯棒性三個約束進行懲罰，目標值Zr可以表示為：

(49)

式(49)中，第2項表示對血液短缺量進行懲罰，第3項表示對預(yù)算超支的懲罰,第4項表示對p-魯棒約束未滿足的懲罰。κ，π，γ分別為染色體需求未滿足、預(yù)算超支、p-魯棒約束未滿足的懲罰因子，應(yīng)為較大的正數(shù)。

④選擇、交叉和變異

采用輪盤賭選擇，并采用精英保留策略將每代的最優(yōu)個體保留進入下一代種群，以加快收斂。采用單點交叉和互換變異操作。

⑤終止條件

達到最大迭代代數(shù)時，停止計算并輸出結(jié)果。

4 算例分析

以中國大陸地區(qū)31個省級血液中心作為儲備庫的候選建設(shè)點，31個省級行政區(qū)為需求點。省級血液中心的日常需求分布見孟超[2]和沈紅艷[3]。在破壞性地震、大規(guī)?？植酪u擊等典型涉血災(zāi)害下，災(zāi)后應(yīng)急血液需求具有明顯的階段性，可分為緊急救援期和救援中后期兩個階段。災(zāi)害發(fā)生后的72小時為緊急救援期，即兩個階段的間隔時間為3天。臨床輸血搶救大多發(fā)生在緊急救援階段，后期需求趨緩，因而假設(shè)中后期血液需求量為緊急救援階段的70%[2]。考慮3種情景，情景1為最壞情景，情景2較嚴重，情景3程度最輕。根據(jù)中國人口危險性等級分布圖及對照表得到31個省級行政區(qū)(需求點)的相對災(zāi)害發(fā)生概率[24-25]。假設(shè)“5.12”地震為四川省的最壞情景?！?.12”期間成都血液中心關(guān)鍵救援期內(nèi)緊急調(diào)運的3類血液制品紅細胞、全血、血漿分別為4235.5U、504U、1780.75U。以此數(shù)據(jù)為基準設(shè)置情景1的總需求量。假設(shè)情景2和3的需求量分別為情景1的50%和20%(災(zāi)害越嚴重，臨床血液需求量越大)。一般情況下，某個需求點傷員的血型比例分布大致與當?shù)乜側(cè)丝诘难捅壤嗤?，因?種血型血液制品的需求量也按照當?shù)乜側(cè)丝诘难捅壤齺矸峙洚a(chǎn)生。各省總?cè)丝谘捅壤植?、省級血液中心編號與經(jīng)緯度分布見周愉峰[24]。候選地的儲備能力上限系數(shù)見孟超[2]。庫存補充系數(shù)為0.3/天[2]。需求點的應(yīng)急血液接收耗時，各候選地調(diào)劑準備耗時均為3小時，應(yīng)急調(diào)配時限為8小時。運輸速度為800km/h，運輸時間通過運輸距離除以運輸速度計算。運輸距離根據(jù)31個城市血液中心的經(jīng)緯度進行計算[3，24]。單位運輸成本為10元/噸公里，庫存持有成本為10元/單位。任意血液制品的庫存服務(wù)水平zjk為1.96。魯棒概率參數(shù)p=0.03。

遺傳算法參數(shù)設(shè)置如下：最大迭代次數(shù)為Maxgen=500，種群規(guī)模popsize=100，交叉概率pc=0.9，變異概率pm=0.1。3個懲罰因子κ=π=γ=100000?；贛atlab?2010a平臺，在Intel Pentium(R) G3260T 2.90GHz CPU，64位操作系統(tǒng)，4G RAM的PC機上運行程序。根據(jù)不同預(yù)算值FB，設(shè)計6個算例，計算結(jié)果見表1(取5次運算的最優(yōu)值作為最終結(jié)果)。圖1為算例5的算法收斂曲線。表2為開放儲備庫的采血批量(算例5)。圖2、圖3分別為算例5(情景1)和算例6(情景2)下的選址-分派方案。計算結(jié)果表明：①系統(tǒng)及時度與預(yù)算投入正相關(guān)，但及時度的改善速度隨著預(yù)算的增加而下降：在預(yù)算投入較少時，少量增加成本可以大幅度提升系統(tǒng)及時度；當預(yù)算增加到一定程度時，再增加投入引起的系統(tǒng)及時度目標改善非常有限(例如，算例5與算例6相比，成本增加14.29%，系統(tǒng)及時度目標改善1.53%；算例1與算例2相比，成本增加66.67%，系統(tǒng)及時度目標僅改善1.73%)。②少數(shù)情景下的魯棒模型對應(yīng)的目標值小于確定性模型的最優(yōu)值，但兩者的相對偏差較小(最大偏差≤1.08%)。因而可得出結(jié)論：雖然p-魯棒解不能在所有情景下實現(xiàn)最優(yōu)目標值, 但適用于所有情景，與確定性模型最優(yōu)值的差距很小。因而p-魯棒優(yōu)化可降低不確定性風險。③血液戰(zhàn)略儲備庫布局呈現(xiàn)出西部較稀疏、中東部較稠密的特征。在庫存決策上，西部地區(qū)建立的儲備庫采血批量較小，中、東部地區(qū)儲備庫采血批量較大。主要原因可能在于西部地區(qū)地廣人稀，血液供需規(guī)模較小，導(dǎo)致西部地區(qū)設(shè)施點數(shù)量、指派規(guī)模、采血批量等低于中、東部地區(qū)。

為了測試遺傳算法的性能，與免疫優(yōu)化算法進行比較。在此，設(shè)計兩組算例。第1組算例基于我國大陸地區(qū)31個省級血液中心與省級行政區(qū)的數(shù)據(jù)，根據(jù)不同預(yù)算值FB，設(shè)計6個算例。第2組算例為6個不同規(guī)模的模擬算例。最大規(guī)模為60×60節(jié)點算例，節(jié)點生成規(guī)則為：在當前31個節(jié)點的基礎(chǔ)上再生成29個隨機節(jié)點，緯度在[20, 45.77]之間隨機均勻產(chǎn)生，經(jīng)度在[87.57,126.68]之間隨機均勻產(chǎn)生(其中，20，45.77分別是31個省級血液中心緯度的最小值與最大值；87.57，126.68分別是31個省級血液中心經(jīng)度的最小值與最大值)。其他參數(shù)在31節(jié)點算例參數(shù)的最小值與最大值之間隨機均勻產(chǎn)生。其他規(guī)模算例的節(jié)點則在60節(jié)點中隨機抽取。6個不同規(guī)模算例的預(yù)算參數(shù)FB分別設(shè)為：1000000，2000000，3000000，4000000，45000000，5000000；庫存補充系數(shù)設(shè)為0.1/天。其他參數(shù)設(shè)置同上文。由表3與表4的12個算例可知，在計算精度上，遺傳算法在7個算例上的表現(xiàn)較好，免疫優(yōu)化算法在3個算例上的表現(xiàn)更優(yōu)，另外2個算例結(jié)果相同。而在計算效率上，遺傳算法明顯優(yōu)于免疫優(yōu)化算法。因此本文提出的遺傳算法具有較好的性能。

表1 6個算例的計算結(jié)果

圖1 遺傳算法收斂曲線(算例5)

表2 已選儲備庫各型血的采血批量(算例5)

圖2 算例5下的選址-分派方案(情景1)

圖3 算例6下的選址-分派方案(情景2)

此外，表3、表4的結(jié)果同樣說明了p-魯棒優(yōu)化可降低不確定性風險(平均偏差≤1.12%)。

表3 遺傳算法與免疫優(yōu)化算法性能對比(31節(jié)點算例)

表4 遺傳算法與免疫優(yōu)化算法性能對比(規(guī)模分析)

表5 魯棒概率參數(shù)p的敏感性分析

為了進一步分析遺傳算法效率，給出算法的計算復(fù)雜度函數(shù)。設(shè)|I|、|L|、|K|、|Jo|分別表示需求點數(shù)量、情景數(shù)量、血型數(shù)量、開放儲備庫數(shù)量。則一次迭代中：初始化的算法復(fù)雜度為O(popsize)；適應(yīng)度評價的算法復(fù)雜度為O(popsize×|I|×|L|×|K|×(2+|Jo|)；選擇、交叉和變異的復(fù)雜度為：O(6×1/2×popsize)；其他操作為單步，可不計。因此，一代計算過程的算法復(fù)雜度Com為：

Com=[O(popsize)+O(6×1/2×popsize)]+O(popsize×|I|×|L|×|K|×(2+|Jo|)≈O(popsize×|I|×|L|×|K|×|Jo|)

(54)

因為|Jo|的理論最大值為|I|，則復(fù)雜度可進一步簡化為：

(55)

最后對魯棒概率參數(shù)p進行敏感性分析(表5)。為了便于分析結(jié)論，取10*10節(jié)點算例進行試驗(小規(guī)模算例更易得到全局最優(yōu)解)。結(jié)果表明：預(yù)算參數(shù)FB足夠大時，參數(shù)p的取值對解的影響小。反之，F(xiàn)B較小時，p的取值對解的影響較大。原因在于FB越小，模型限制越強，可行解空間越小，對魯棒概率參數(shù)p的限制也越強。

5 結(jié)語

國家血液戰(zhàn)略儲備庫的建設(shè)問題亟待解決。以應(yīng)急血液保障及時度最高為目標，考慮多情景不確定需求、多血型、多階段、設(shè)施容量限制及協(xié)同定位等災(zāi)時應(yīng)急特性，同時考慮日常運營成本與預(yù)算約束，建立了一個國家血液戰(zhàn)略儲備庫LIP的隨機p-魯棒優(yōu)化模型。采用非線性混合整數(shù)規(guī)劃方法描述該模型。并基于模型性質(zhì)，設(shè)計了一種遺傳算法對問題進行求解。通過多組算例驗證并與免疫優(yōu)化算法進行對比, 證明了所提模型和算法的有效性。本文給出了若干組不同預(yù)算值下的選址-庫存決策方案。實踐中，決策者可根據(jù)預(yù)算需求及偏好，選擇合適的選址-庫存決策方案。

研究結(jié)論及啟發(fā)是：①系統(tǒng)及時度與預(yù)算投入正相關(guān)，但及時度的改善速度隨著預(yù)算的增加而大幅下降。因而需要根據(jù)實際情況，權(quán)衡預(yù)算投入與系統(tǒng)及時度目標。在投入較小時，可以少量增加預(yù)算，以大幅度提高應(yīng)急血液保障的效率；反之，若預(yù)算過大，可以考慮適當縮減投資，在應(yīng)急血液保障效率損害很小的前提下大量降低成本，②雖然p-魯棒解不能在所有情景下取得最優(yōu)值, 但與確定性情景模型最優(yōu)值相差不大。由于國家血液儲備網(wǎng)絡(luò)設(shè)計問題屬于戰(zhàn)略性決策問題，因而必須考慮網(wǎng)絡(luò)構(gòu)建的魯棒性，以降低不確定性風險。

本文以國家血液戰(zhàn)略儲備庫為研究對象，實際工作中，可對模型稍作改進，應(yīng)用于同樣存在多情景多階段不確定應(yīng)急需求與日常隨機需求、多產(chǎn)品、有容量限制及協(xié)同定位、改擴建與預(yù)算約束等特性的易腐品(藥品、糧食等)應(yīng)急物資儲備庫LIP中。進一步研究可以考慮血液的新鮮度因素，擴展本文提出的LIP；也可考慮設(shè)施中斷風險等因素，對選址-庫存模型進行可靠性優(yōu)化。