林 勇，張 立，汪貽生，肖 驊

(1.陸軍勤務(wù)學(xué)院 軍事物流系，重慶 401311；2.重慶財(cái)經(jīng)學(xué)院 科研處，重慶 401320)

0 引言

應(yīng)急物資政企協(xié)同預(yù)置儲備因能有效破解政府單方面預(yù)置儲備存在的供應(yīng)風(fēng)險高、財(cái)政負(fù)擔(dān)重、輪換更新難等問題[1-3]，被認(rèn)為是高效應(yīng)對災(zāi)害的必然出路。有關(guān)應(yīng)急物資政企協(xié)同預(yù)置儲備優(yōu)化的研究已經(jīng)引起學(xué)者們的關(guān)注，但研究成果主要集中在災(zāi)前預(yù)置量、災(zāi)后采購量和災(zāi)后運(yùn)輸量的聯(lián)合決策方面且研究成果較少[4-13]。同時，現(xiàn)有研究均未考慮預(yù)置儲備過程中的輪換更新環(huán)節(jié)，即默認(rèn)預(yù)置儲備期間不會進(jìn)行應(yīng)急物資輪換更新。但實(shí)際輪換更新是應(yīng)急物資預(yù)置儲備實(shí)踐中必不可少的環(huán)節(jié)[14-15]。鑒于此，本文在分析輪換更新環(huán)節(jié)基礎(chǔ)上，構(gòu)建計(jì)及輪換更新費(fèi)用的應(yīng)急物資政企協(xié)同預(yù)置儲備優(yōu)化模型，并通過數(shù)值算例驗(yàn)證模型的有效性。

1 問題描述

1)應(yīng)急準(zhǔn)備階段。政府救災(zāi)部門與應(yīng)急物資生產(chǎn)企業(yè)(以下簡稱“企業(yè)”)達(dá)成協(xié)議，承諾采購不低于約定數(shù)量的應(yīng)急物資，并儲備至政府下屬的儲備庫。協(xié)議同時約定，企業(yè)必須預(yù)置一定規(guī)模的庫存物資以備災(zāi)后應(yīng)急所需。為確保應(yīng)急物資長儲常新，儲備庫和企業(yè)都需要定期對應(yīng)急物資進(jìn)行輪換更新。

2)應(yīng)急響應(yīng)階段。政府救災(zāi)部門可通過應(yīng)急調(diào)撥儲備庫物資和應(yīng)急采購企業(yè)庫存物資2種方式進(jìn)行應(yīng)急救援。應(yīng)急物資輪換更新涉及折舊出售和補(bǔ)充更新2步操作：儲備庫與企業(yè)將擬輪換的應(yīng)急物資折舊出售給零售商(如大型超市)；儲備庫與企業(yè)通過協(xié)議采購以及組織生產(chǎn)的方式補(bǔ)充更新因輪換而產(chǎn)生的應(yīng)急物資空缺。在輪換更新周期內(nèi)，受外部環(huán)境等因素影響，應(yīng)急物資價值通常隨儲備時間推移不斷下降，且不同類型應(yīng)急物資的價值損失變化率存在差異[16]。應(yīng)急物資價值與預(yù)置儲備時間函數(shù)關(guān)系如式(1)所示：

(1)

式中：C0為應(yīng)急物資初生產(chǎn)出來時的價值，當(dāng)應(yīng)急物資預(yù)置儲備達(dá)到保質(zhì)期T0時，其價值為0；T0為保質(zhì)期，a；t為預(yù)置儲備時間，a；n為形狀參數(shù)。

2 模型構(gòu)建

2.1 模型假設(shè)

為明確研究邊界，作如下假設(shè)：

1)應(yīng)急物資的保質(zhì)期T0以及輪換更新周期Tg已知，且Tg

2)應(yīng)急物資輪換產(chǎn)生的勞務(wù)費(fèi)以及運(yùn)輸費(fèi)由零售商承擔(dān)。

3)儲備庫與企業(yè)每次輪換更新過程中補(bǔ)充更新的應(yīng)急物資為剛生產(chǎn)出來的全新物資。

4)協(xié)議期T內(nèi)，無災(zāi)害事件發(fā)生。

2.2 數(shù)學(xué)模型

構(gòu)建考慮輪換更新的應(yīng)急物資政企協(xié)同預(yù)置儲備優(yōu)化模型。模型目標(biāo)函數(shù)如式(2)所示：

(2)

式中：f表示總成本，美元；fc為儲備庫設(shè)立成本，美元；pc為初次儲備采購成本，美元；tc為采購以及輪換更新產(chǎn)生的運(yùn)輸成本，美元；hc為儲備庫以及企業(yè)預(yù)置應(yīng)急物資產(chǎn)生的庫存維持成本，美元；dc為儲備庫和企業(yè)進(jìn)行R次輪換更新產(chǎn)生的差價(折舊)補(bǔ)貼，美元；pcs為災(zāi)后應(yīng)急采購成本，美元；tcs為災(zāi)后運(yùn)輸成本，美元；S為災(zāi)害情景集合；ωs為情景s的發(fā)生概率。各項(xiàng)成本如式(3)～(9)所示：

(3)

式中：Fm為m型儲備庫的設(shè)立成本，美元；xim為0～1變量，若在候選點(diǎn)i設(shè)立m型儲備庫，則xim為1；I為儲備庫候選點(diǎn)集合，i∈I；M為儲備庫大小類型，m∈M。

(4)

式中：Cl為l型企業(yè)的應(yīng)急物資單位采購成本，美元；pikl為候選點(diǎn)i的儲備庫向候選點(diǎn)k的l型企業(yè)進(jìn)行災(zāi)前采購的應(yīng)急物資量；K為企業(yè)的候選點(diǎn)集合k∈K；L為企業(yè)大小類型，l∈L。

(5)

式中：R為協(xié)議期內(nèi)應(yīng)急物資的輪換更新次數(shù)，次；TC為單位物資量、單位距離的運(yùn)輸成本，美元；Hik為候選點(diǎn)i的儲備庫與候選點(diǎn)k的企業(yè)之間的距離，km。

(6)

式中：eim為候選點(diǎn)i的m型儲備庫應(yīng)急物資儲備量；HC為單位應(yīng)急物資的年庫存維持成本，美元；T為政府與企業(yè)的合同協(xié)議期限，a；gkl為候選點(diǎn)k的l型企業(yè)需要預(yù)置的庫存量。

(7)

式中：J為需求點(diǎn)集合，j∈J。

(8)

式中：ε為災(zāi)后應(yīng)急采購單價補(bǔ)償系數(shù)；rkjls為情景s下，候選點(diǎn)k的l型企業(yè)發(fā)售給需求點(diǎn)j的應(yīng)急物資量。

(9)

式中：qijs為情景s下，候選點(diǎn)i的儲備庫調(diào)撥給需求點(diǎn)j的應(yīng)急物資量；Hij為候選點(diǎn)i的儲備庫與需求點(diǎn)j的距離，km；Hkj為候選點(diǎn)k的企業(yè)與需求點(diǎn)j的距離，km。

模型約束條件如式(10)～(21)所示。其中，儲備庫和企業(yè)類型數(shù)量限制如式(10)～(11)所示：

(10)

(11)

最低采購量和有效采購量約束如式(12)所示：

(12)

式中：ykl為0～1變量，若選擇候選點(diǎn)k的l型企業(yè)，則為1；Ql為政府向l型企業(yè)進(jìn)行災(zāi)前采購的最低協(xié)議訂購量；B為較大的正數(shù)。儲備庫容量限制如式(13)～(14)所示：

(13)

eim≤xim·Am?i∈I，m∈M

(14)

式中：Am為m型儲備庫的最大儲備容量。災(zāi)后儲備庫的調(diào)撥量不能超過其有效儲備量的限制，如式(15)所示：

(15)

式中：αis為災(zāi)害情景s發(fā)生后，候選點(diǎn)i的儲備庫可用應(yīng)急物資比例。災(zāi)后采購量不能超過企業(yè)的有效預(yù)置庫存量的限制如式(16)所示：

(16)

式中：αks為災(zāi)害情景s發(fā)生后候選點(diǎn)k的企業(yè)可用應(yīng)急物資比例。企業(yè)預(yù)置庫存量不能超過企業(yè)的庫存容量限制，如式(17)所示：

gkl≤ykl·Vl?k∈K，l∈L

(17)

式中：Vl為l型企業(yè)的預(yù)置庫存最大容量。需求滿足限制如式(18)所示：

(18)

式中：Djs為情景s下需求點(diǎn)j的應(yīng)急物資需求量。0～1決策變量約束如式(19)～(20)所示：

xim∈{0,1} ?i∈I，m∈M

(19)

ykl∈{0,1} ?k∈K，l∈L

(20)

變量的非負(fù)整數(shù)約束如式(21)所示：

pikl，gkl，qijs，rkjls，eim∈N，?i,j,k,m,l,s

(21)

式中：N為自然數(shù)集合(非負(fù)整數(shù)集)。

3 算例分析

3.1 算例設(shè)置

以美國東南部墨西哥灣颶風(fēng)災(zāi)害為背景[13]，選取3種代表不同價值損失變化率和輪換更新周期的應(yīng)急物資為研究對象，研究計(jì)及輪換更新費(fèi)用的應(yīng)急物資政企協(xié)同預(yù)置儲備優(yōu)化問題。假定3種代表性的應(yīng)急物資分別為E1，E2，E3。E1代表方便面，設(shè)其n=0.5，Tg=0.8 a，T0=1 a；E2代表桶裝水，設(shè)其n=1，Tg=1.5 a，T0=2 a；E3代表應(yīng)急燈，設(shè)其n=4，Tg=2.6 a，T0=3 a。儲備庫和企業(yè)分為大、中和小3種類型，各類型儲備庫和企業(yè)的參數(shù)見表1。其余參數(shù)設(shè)置如下：ε=0.5，TC=0.003美元，HC=0.65美元，T=3 a。

表1 與儲備庫和企業(yè)相關(guān)的模型參數(shù)設(shè)置Table 1 Setting of model parameters related to reserve warehouses and enterprises

3.2 結(jié)果分析

考慮模型為線性整數(shù)規(guī)劃模型，本文借助Anaconda平臺的Spyder集成開發(fā)環(huán)境，采用Python語言進(jìn)行編碼，并調(diào)用Gurobi求解器進(jìn)行求解。經(jīng)計(jì)算得到以上3種應(yīng)急物資對應(yīng)政企協(xié)同預(yù)置儲備優(yōu)化結(jié)果，見表2。

表2 3種應(yīng)急物資對應(yīng)的政企協(xié)同預(yù)置儲備優(yōu)化結(jié)果Table 2 Optimization results of government-enterprise collaborative preset reserve corresponding to three kinds of emergency supplies

優(yōu)化結(jié)果驗(yàn)證本文模型的有效性，能夠?qū)哂胁煌啌Q更新周期的應(yīng)急物資政企協(xié)同預(yù)置儲備問題進(jìn)行優(yōu)化求解。由表2可知，3種應(yīng)急物資E1，E2，E3的輪換更新次數(shù)R分別為3，2，1次，但3種應(yīng)急物資在儲備過程中需要花費(fèi)的差價(折舊)補(bǔ)貼普遍較高，3者差價(折舊)補(bǔ)貼與總成本的占比分別為71.48%，69.43%，60.34%。由此可見，在應(yīng)急物資政企協(xié)同預(yù)置儲備過程中，由于輪換更新環(huán)節(jié)的存在，差價(折舊)補(bǔ)貼是不可忽視的成本支出。另一方面，從庫存維持成本來看，盡管3種應(yīng)急物資的輪換更新次數(shù)存在差異，但3者的庫存維持成本是一樣的。這是因?yàn)閹齑婢S持成本由預(yù)置儲備量、單位應(yīng)急物資的年庫存維持成本以及合同期限共同決定，與輪換更新次數(shù)無關(guān)。在本文算例中，3種應(yīng)急物資的儲備期限相同，即T=3 a。

此外，由式(5)可知，在應(yīng)急物資政企協(xié)同預(yù)置儲備過程中，輪換更新次數(shù)除對差價(折舊)補(bǔ)貼有直接影響外，還對災(zāi)前運(yùn)輸成本有著直接的影響。在差價(折舊)補(bǔ)貼成本較高的情況下，為實(shí)現(xiàn)總成本最小，模型優(yōu)化結(jié)果勢必會對其它成本進(jìn)行折中，這也是本例中3種應(yīng)急物資的災(zāi)前運(yùn)輸成本普遍較低的主要原因。災(zāi)前運(yùn)輸成本較低，意味著應(yīng)急物資的災(zāi)前儲備采購以及補(bǔ)充更新主要通過與就近企業(yè)對接完成，例如對于應(yīng)急物資E1而言，其災(zāi)前運(yùn)輸成本為0，表示各儲備庫全部通過儲備庫所在地的企業(yè)進(jìn)行災(zāi)前儲備采購和輪換更新過程中的補(bǔ)充更新。

由此可見，輪換更新次數(shù)通過影響災(zāi)前運(yùn)輸成本，對優(yōu)化結(jié)果中儲備庫與企業(yè)的協(xié)同關(guān)系產(chǎn)生影響。在合同期內(nèi)，對于輪換更新較為頻繁的應(yīng)急物資，選擇與儲備庫較近的企業(yè)完成災(zāi)前儲備采購以及輪換更新是政企雙方開展協(xié)同預(yù)置儲備的最佳方式。相反，對于輪換更新頻次不高的應(yīng)急物資，可以選擇離儲備庫有一定距離的企業(yè)完成災(zāi)前儲備采購與輪換更新。

4 結(jié)論

1)合同期內(nèi)，無論輪換更新次數(shù)多少，應(yīng)急物資政企協(xié)同預(yù)置儲備需要支出的差價(折舊)補(bǔ)貼普遍較高。算例中，3種應(yīng)急物資因輪換更新產(chǎn)生的差價(折舊)補(bǔ)貼大于總成本的60%。

2)輪換更新次數(shù)不會影響庫存維持成本，但會影響儲備庫與企業(yè)的協(xié)同關(guān)系。合同期內(nèi)，對于輪換更新較頻繁的應(yīng)急物資，政府最好選擇離儲備庫較近的企業(yè)完成災(zāi)前儲備采購以及輪換更新；對于輪換更新次數(shù)不多的應(yīng)急物資，政府可以選擇離儲備庫有一定距離的企業(yè)完成災(zāi)前儲備采購與輪換更新。