考慮時效與公平性的震后應(yīng)急物資動態(tài)配送優(yōu)化研究

曲沖沖，王　晶，黃　鈞，何明珂

(1.北京工商大學(xué)商學(xué)院，北京　100048;2.中國科學(xué)院大學(xué)工程科學(xué)學(xué)院　100049;3.北京物資學(xué)院物流學(xué)院　101149)

1　概述

大規(guī)模地震嚴重威脅人民生命和財產(chǎn)安全。加強震后應(yīng)急物資保障是應(yīng)對突發(fā)自然災(zāi)害條件下提高救援效率的關(guān)鍵。震后物資保障需要將大量的應(yīng)急物資運輸?shù)绞転?zāi)區(qū)域，應(yīng)急物資運輸?shù)捻憫?yīng)速度和分配策略是救援工作開展的重點任務(wù)。同時，根據(jù)受災(zāi)區(qū)域的破壞程度，合理選擇運輸工具和穩(wěn)定的配送中心是保證救援工作有序開展的前提。

震后應(yīng)急配送中心選址和物資運輸分配問題是國內(nèi)外學(xué)者廣泛研究的課題。震后應(yīng)急救援物流不同于一般的社會性物流，其活動開展具有不確定性、突發(fā)性、時間約束較強、經(jīng)濟性較弱等特點。宮華等[1]從外圍物資集散的中心選址、救災(zāi)物品的分配滿意度等方面進行了研究。沈曉冰等[2]在分析震后應(yīng)急物資的配送特點和應(yīng)急配送中心的選址問題的基礎(chǔ)上，建立了一個綜合考慮應(yīng)急物流運作效率和運輸路徑選擇的多目標(biāo)優(yōu)化模型，對震后救援效率和路徑選擇進行研究。劉長時等[3-4]在對震后應(yīng)急物資運輸?shù)缆芬?guī)劃問題進行研究時，通過建立多種運輸方式、多運輸周期的LRP模型，對應(yīng)急物資開展運輸?shù)臅r間進度進行討論；在對震后災(zāi)民的非理性攀比因素進行研究時，通過構(gòu)建一個靜態(tài)的運輸-分配優(yōu)化模型，對靜態(tài)條件下的震后救援工作的開展進行了討論。王海軍等[5]通過建立一個雙目標(biāo)的OLRP模型，對配送中心的選址和應(yīng)急物資的運輸安排進行了討論。 Sheu[6]對臺灣震后救援工作開展的研究中，通過遞推機制對配送中心和受災(zāi)地區(qū)進行分組分析。Najafi等[7]在對多目標(biāo)、多線程、多種物資開展震后運輸和救助傷員的研究中，通過設(shè)計時間優(yōu)化算法對救援工作有序開展提供建議。Yi Wei等[8]對研究震后響應(yīng)活動中的協(xié)調(diào)后勤支援及疏散作業(yè)的綜合作業(yè)定位分配模式進行了研究。Haghi等[9]通過對震后受傷人員的救治路徑和健康中心布局規(guī)劃研究，提升了救助覆蓋的全面性。Arora等[10]在對公共衛(wèi)生事件下開展物資、藥品分配的研究的基礎(chǔ)上，認為政策因素和距離因素對物資分配的公平性有一定的影響。Altay[11]認為提升救災(zāi)效率和公平分配物資對傷員救治和災(zāi)后區(qū)域恢復(fù)有著積極的促進作用。Liu Cong等[12]認為將救援工作從時間上劃分成若干個時段可以提升震后救援的效率。

對應(yīng)急物資配送中心選址、運輸路徑選擇和物資分配公平性開展研究時，眾多文獻是通過建立數(shù)學(xué)模型進行分析的。李孟良等[13]在研究震后多式聯(lián)運問題時，對未滿足災(zāi)區(qū)應(yīng)急物品需求的最小懲罰系數(shù)進行了討論，提出了基于魯棒優(yōu)化的聯(lián)運調(diào)配模型。鄭斌等[14]針對震后物資供不應(yīng)求的狀況，通過設(shè)計一種混合遺傳算法結(jié)合汶川地震案例對物資的運輸時間規(guī)劃、物資分配方式進行了分析。Abounacer等[15]設(shè)計了一種啟發(fā)式算法—ε-約束方法，旨在確定實際所需的分配中心數(shù)量、位置和救援作業(yè)時間規(guī)劃。張偉等[16]在討論考慮運輸時間、運輸距離和道路復(fù)雜程度對路徑選擇影響的問題時設(shè)計了一種多目標(biāo)非線性整數(shù)規(guī)劃模型。Mahmoudabad等[17]提出了一種基于動態(tài)風(fēng)險和動態(tài)損傷烈度網(wǎng)絡(luò)的混沌理論解決突發(fā)情況下信息中斷時救援路徑選擇和物資分配的問題。王旭坪等[18]在借鑒前景理論的基礎(chǔ)上，通過構(gòu)建風(fēng)險感知函數(shù)，從道路運力、物資流動、決策過程等方面建立了物資調(diào)配的系統(tǒng)動力學(xué)模型。Haghani和Oh[19]通過交互式啟發(fā)算法解決救災(zāi)物資分配和帶時間窗問題的多模態(tài)網(wǎng)絡(luò)流問題。Najafi等[6]用成本分析的方法提出了多目標(biāo)、多模式、多商品、多周期的隨機模型為震后物資分配提供指導(dǎo)。

眾多文獻設(shè)計的目標(biāo)函數(shù)通常考慮救援成本和時間效率，對受災(zāi)點物資需求滿足程度與分配公平性考慮較少。本文在解決多時段、多種運輸方式共同參與的應(yīng)急物資保障問題的過程中兼顧時效性與公平性，在綜合物資運輸與配送的兩階段、多時段的動態(tài)救援特點和多種運輸方式共同參與的基礎(chǔ)上，建立了保證救災(zāi)物資運輸時效性和實現(xiàn)災(zāi)區(qū)物資分配公平的最優(yōu)的應(yīng)急物資配送中心選址與運輸配送路徑優(yōu)化模型。針對此模型設(shè)計了一種帶精英策略的非支配排序的遺傳算法，結(jié)合九寨溝地區(qū)震后災(zāi)難情景對模型進行實證研究，通過結(jié)合Pareto前沿面分析，驗證了模型和算法的有效性，為震后應(yīng)急物資保障提供了理論指導(dǎo)與決策支持。

2　模型建立

2.1　問題描述

震后應(yīng)急物資保障工作是依托于物資集散地、應(yīng)急物資配送中心和受災(zāi)點之間的相互聯(lián)系展開的。災(zāi)區(qū)所需的大批量、多批次、多品種的物資需要從災(zāi)區(qū)外圍的救災(zāi)物資集散地通過公路和直升飛機向災(zāi)區(qū)周圍的臨時救災(zāi)物資存放地應(yīng)急物資配送中心進行分配中轉(zhuǎn)，這是應(yīng)急物資運輸?shù)牡谝浑A段；再由臨時配送中心通過公路運輸向災(zāi)區(qū)進行分配，這是應(yīng)急物資運輸?shù)牡诙A段。在兩個階段中根據(jù)不同時段的道路狀況、運輸裝備性能開展的運輸路線調(diào)整和應(yīng)急物資配送中心的備選工作是該階段的重點。兩階段劃分示意圖如圖1所示：

圖1　災(zāi)后應(yīng)急物資運輸系統(tǒng)示意圖

震后救災(zāi)物資主要的供應(yīng)模式是區(qū)域儲備倉庫和外圍物資集散地整合物資-區(qū)域應(yīng)急物資配送中心分配中轉(zhuǎn)-受災(zāi)點，受到運輸車數(shù)量和車輛的負荷裝載能力、路況、環(huán)境、救援人員數(shù)量變化等一系列客觀因素的制約以及隨著震后救援工作不斷深入，對災(zāi)區(qū)實際狀況逐漸了解，需要對較長時間的救援工作進行時間細分，通過多時段的動態(tài)救援滿足災(zāi)區(qū)的物資供應(yīng)。

通過將應(yīng)急物資保障分為若干時段，可以及時的調(diào)整救援工作，從而保證救援工作的高效率運作。同時要兼顧物資運輸?shù)臅r效性和物資分配的的公平性。

2.2　模型構(gòu)建

2.2.1模型假設(shè)

針對特定條件下救災(zāi)物資的配送特征，作出如下假設(shè)：

(1)備選配送中心都是不易受到余震以及其他次生災(zāi)害影響且距離受災(zāi)區(qū)域距離較短的的交通樞紐區(qū)域點，可以長時間進行救災(zāi)作業(yè)。

(2)應(yīng)急物資從區(qū)域儲備倉庫和外圍集散地還未運輸?shù)絽^(qū)域應(yīng)急配送中心時，受災(zāi)區(qū)域儲備的物資可以滿足自身的最低需求；區(qū)域應(yīng)急配送中心向受災(zāi)點運輸?shù)膽?yīng)急物資數(shù)量可以滿足所有受災(zāi)點的最低需求。

2.2.2符號說明

(1)參數(shù)

W：救援物資種類的集合(w∈W)

I：應(yīng)急物資集散點的集合(i∈I)

J：應(yīng)急物資配送中心的集合(j∈J)

K：受災(zāi)點的集合(k∈K)

Ξ：需要決策的時段的集合(ξ∈Ξ)

M=

T(ξ)：決策時段ξ重新決策的時刻

qwi(ξ)：ξ時段下集散點i供給物資w在的量

dwk(ξ)：受災(zāi)區(qū)域k在ξ時段對物資w的實際需求量

Pijm(ξ)：在ξ時段下從集散點i到配送中心j通過運輸方式m單次最大運輸量

Tijm(ξ)：在ξ時段道路未修復(fù)情況下，從集散點i到配送中心j通過運輸方式m所需要的時間

Pjk(ξ)：ξ時段下從配送中心j到受災(zāi)點k的單次最大運輸量

Tjk(ξ)：ξ時段道路未修復(fù)情況下，從配送中心j到受災(zāi)點k所需要的時間

ITijm：從集散點i到配送中心j通過運輸方式m發(fā)車的時間間隔

θw：受災(zāi)點人均需求救災(zāi)物資w的最少量

LTkw(ξ)：ξ時段下受災(zāi)點k對救濟物資w供應(yīng)滿意度最大時，所能夠接受的最長等待時間

STkw(ξ)：ξ時段下受災(zāi)點k對救濟物資w供應(yīng)不滿意度最大時，所能接受的最短等待時間

N(ξ)：ξ時段下允許配送中心運營的最多數(shù)量

PDkw(ξ)：ξ時段下進行下一時段(ξ+1)決策時，受災(zāi)點k接收救災(zāi)物資w必須等待的時間

Tkw(ξ)：ξ時段下受災(zāi)點k完全接受救災(zāi)物資w的時刻

QKw(ξ)：ξ時段下計劃分給受災(zāi)區(qū)域K的救災(zāi)物資w的總量

Θ：可用車輛總數(shù)

Ω：一個大數(shù)

(2)決策變量

Qijwm(ξ)：ξ時段下計劃采用運輸方式m從物資集散點i運輸?shù)脚渌椭行膉物資w的數(shù)量

?j(ξ)：ξ時段下計劃分配給配送中心j進行運輸物資的車輛數(shù)目

Qjkw(ξ)：ξ時段下計劃從配送中心j運輸?shù)绞転?zāi)點k的救災(zāi)物資w的量

2.2.3數(shù)學(xué)模型構(gòu)建

(1)

(2)

(3)

(4)

(5)

(6)

(7)

(8)

(9)

(10)

(11)

(12)

Qijwm(ξ)≥0,i∈I,j∈J,w∈W,m∈M

(13)

uijmw(ξ)=(0,1),i∈I,j∈J,w∈W,m∈M

(14)

(15)

(16)

(17)

Tjk(ξ),k∈K,w∈W,(j,k)∈Ak(ξ)

(18)

(19)

Qjkw(ξ)≤Ω·zjkw(ξ),j∈J,k∈K,w∈W

(20)

(21)

(22)

zjkw(ξ)=(0,1)

(23)

目標(biāo)函數(shù)(1)表示在ξ時段下，各種應(yīng)急物資配送到災(zāi)區(qū)的加權(quán)時刻最短，進而保證應(yīng)急物資的時效性。目標(biāo)函數(shù)(2)表示最大化所有受災(zāi)點各種物資滿意度的最小值，保證救援物資分配的公平性。約束(3)表示在ξ時段下計劃運輸物資到配送中心的總量等于該集散點物資總供應(yīng)量；約束(4)表示受災(zāi)點物資最小滿足率；約束(5)表示配送中心的最大承載量；約束(6)表示在不同運輸方式下，將救災(zāi)物資從集散點運輸?shù)脚渌椭行乃钑r間表達式；約束(7)表示第一階段物資運輸時間函數(shù)表達式；約束(8)表示受災(zāi)點完全接收救災(zāi)物資時刻的函數(shù)表達式；約束(9)表示ξ時段下只有選擇某個配送中心時，才會把物資運輸?shù)皆撆渌椭行摹＜s束(10)表示救災(zāi)物資只能向運轉(zhuǎn)的配送中心運輸救災(zāi)物資；約束(11)表示配送中心開始運轉(zhuǎn)后，任何階段都不會關(guān)閉；約束(12)表示可以運轉(zhuǎn)的配送中心小于等于初始設(shè)定的最大數(shù)目限制；約束(13)表示非負約束；約束(14)、(15)表示0-1約束。約束(16)表示配送中心接受物資的量與輸出給災(zāi)區(qū)的量均等；約束(17)表示各個受災(zāi)點的物資接收量總和等于配送中心運往災(zāi)區(qū)的量；約束(18)表示通過公路運輸從配送中心到受災(zāi)點的時間函數(shù)表達式；約束(19)表示由于物資的供不應(yīng)求，運輸?shù)绞転?zāi)地區(qū)的物資總量小于災(zāi)區(qū)實際需求量；約束(20)表示只有選擇為某個受災(zāi)點服務(wù)時，才會輸送物資到該點；約束(21)表示受災(zāi)點接受的救災(zāi)物資只能來自于已經(jīng)開啟運轉(zhuǎn)的配送中心。約束(22)表示受災(zāi)點接受的物資只能來自于一個配送中心。約束(23)表示0-1約束。震后應(yīng)急響應(yīng)時間滿意度通過以下線性時間滿意度函數(shù)表達式描述：

在第一階段中，公路運輸與直升飛機運輸共同參與到應(yīng)急物資的輸送工作，救援車輛從集散點到配送中心的出發(fā)時刻根據(jù)物資運輸方式的不同可以分為以下兩種情況：

i.當(dāng)運輸方式為公路運輸且當(dāng)前時段沒有救援車輛向災(zāi)害地點運輸救災(zāi)物資。

ii.運輸方式為公路運輸，且在下一時段開始

決策時，當(dāng)前時段的運輸工作已經(jīng)完成。

i.當(dāng)運輸方式為公路運輸時，且下一時段救援工作開始時，當(dāng)前時段的車輛仍在途中。

ii.運輸方式為直升飛機時，由于受到航線的限制，下一時段的車輛發(fā)出時刻是當(dāng)前時段車輛返回的最早時刻。

(24)

(25)

同時由于災(zāi)后救援物資需求具有持需時間長、多時段等特點，本文將救援過程細化成若干時段，目標(biāo)函數(shù)考慮各個時段配送方案的設(shè)計與評價，保證在任意時段下救援物資配送的時效性與公平性。模型可轉(zhuǎn)化為：

s.t.(1)-(25)

3　算法設(shè)計

為了實現(xiàn)災(zāi)后應(yīng)急物資運輸與分配的時效性與公平性，構(gòu)建了包括多時段動態(tài)的配送中心選址和應(yīng)急物資配送路徑優(yōu)化的多目標(biāo)規(guī)劃模型，并設(shè)計了一種基于Pareto最優(yōu)解基礎(chǔ)上改進的帶精英策略非支配排序的多目標(biāo)遺傳算法，將精英策略引入到模型中，提升了算法了運算效率和準(zhǔn)確率。通過基因交叉、多點變異將父代種群與子代種群組合，保證優(yōu)良的種群個體在進化過程中不會被丟棄，進一步提高優(yōu)化結(jié)果的精度。通過對所有個體的分層存放，迅速提高種群水平。

3.1　初始化參數(shù)

設(shè)該遺傳算法的種群規(guī)模popsize，最大迭代次數(shù)為max gen，交叉概率Pc，變異概率Pm。

3.2　染色體編碼

通過矩陣編碼策略描述每個時段配送中心到受災(zāi)點的指派關(guān)系zjkw(ξ)和分配量Qjkw(ξ)。每個受災(zāi)點在每個時段內(nèi)的需求滿足情況是通過列向量基因位表現(xiàn)在染色體上。在每個列向量中第一個基因位是通過在1到|J|的整數(shù)內(nèi)隨機產(chǎn)生，表示配送中心到受災(zāi)點之間的指派關(guān)系。第二基因位到第(w+1)位是在[1，dwk(ξ)]的整數(shù)內(nèi)隨機產(chǎn)生，表示在每個時段下配送中心向受災(zāi)點運送的不同物資的具體的數(shù)量。因此，染色體編碼的總長度為(1+w)×(1+ξ)。

假設(shè)受災(zāi)點的數(shù)量是16個，即K={k|1,2,3,4,5,6,7,8,…14,15,16}，受災(zāi)點對應(yīng)急物資的需求種類為2種：食品和生活用品。救援分為3個時段，即Ξ={ξ|1,2,3}。配送中心候選地點有8個，即|J|=8,選從中選出4個點作為開放的配送中心。因此，可以得出以下染色體編碼，見圖3。

圖3　染色體矩陣編碼示意圖

圖3矩陣含義：選擇2，5，6，8物資配送中心進行運營為16個受災(zāi)地點運輸生活用品和食品。當(dāng)ξ=1時，物資配送中心5為災(zāi)區(qū)1和災(zāi)區(qū)15服務(wù)，向災(zāi)區(qū)1運送24單位和40單位的應(yīng)急物資。同時向災(zāi)區(qū)15運輸27單位和33單位的應(yīng)急物資。

配送中心所屬車輛的數(shù)量是其向受災(zāi)點運輸物資的數(shù)量來決定的，具體的運算公式為：

3.3　種群初始化及適應(yīng)度函數(shù)計算

為了使初始種群中染色體數(shù)量盡可能豐富，在滿足約束條件的基礎(chǔ)上隨機均勻產(chǎn)生各基因位上的值，使得每個受災(zāi)點K都能夠得到救援。在遺傳算法程序設(shè)計上，對各個受災(zāi)點救援次數(shù)進行跟蹤，對于觀測到救援次數(shù)偏少的受災(zāi)點將隨機選擇種群中沒有對其進行救助的個體，隨機將若干個基因位設(shè)置為其對應(yīng)的受災(zāi)點序號。

3.4　帶精英策略的非支配排序解集構(gòu)造

非支配解排序的流程如下：

對于種群中每一個個體i，設(shè)都有兩個參數(shù)ni和Si，ni是指種群中支配個體i的解的個體數(shù)量，Si是指被個體i支配的解的個體集合。

步驟1：找出種群中所有ni=0的個體，表示個體i支配其他個體，將他們存入當(dāng)前非支配集合Z1中；

步驟2：非支配集合Z1中的每一個個體j，遍歷其所支配的個體集合Sj中每一條染色體t，令nt=nt-1。若nt=0，將個體t存入集合H。

步驟3：Z1作為第一級非支配個體的集合，在Z1中的解得個體都是最優(yōu)的。令其非支配序rank(i)=1，在此基礎(chǔ)對集合H作上述分級操作，直到所有的個體都被分級，同時賦予相應(yīng)的非支配排序。

3.5　擁擠度與比較算子

擁擠度的計算是確保種群多樣性的一個重要因素，計算步驟如下：

(1)每個點的擁擠度id置為0；

(2)根據(jù)每個優(yōu)化目標(biāo)，令邊界上的兩個個體擁擠度為無限大，即od=ld=∞；

(3)對種群中的其它個體的擁擠度進行運算：

其中，id表示i點的擁擠度，fji+1表示i+1點的第j個目標(biāo)函數(shù)值，fji-1表示i-1點的第j個目標(biāo)函數(shù)值。

個體i與另一個體j進行比較，只要有以下一個條件滿足，則個體i勝利。

(1)個體i所處的支配層優(yōu)于個體j所處的支配層，即ranki

(2)若種群中有兩個個體有相同等級(處在相同的非支配層)，且個體i比個體j有一個更大的擁擠距離，即irank=jrank且id>jd。

3.6　精英策略與遺傳操作

將父代種群與其產(chǎn)生的子代種群混合后進行非支配排序，可以有效的避免父代種群中的精英個體流失。通過輪盤選舉對個體進行選擇操作，采用單點、多點交叉和基因互換進行變異操作。設(shè)計算法的流程圖如圖2。

圖2　帶精英策略的非支配模型遺傳算法流程

4　算例分析

4.1　算例介紹與參數(shù)設(shè)置

根據(jù)四川九寨溝地震區(qū)域的災(zāi)難情景，選取8個應(yīng)急物資配送中心(編號為1～8)，以及16個受災(zāi)點(編號1～16)同時選取2中主要的應(yīng)急物資(選取食品和生活用品)。震后救援初始階段選擇4個應(yīng)急配送中心開放運營。九寨溝縣各應(yīng)急物資配送中心的倉庫容量如表1，主要受災(zāi)點的物資需求量見表2。

表1　應(yīng)急物資配送中心庫存容量

表2　主要受災(zāi)點的物資需求如下

圖3　Pareto前沿面解的分布圖

帶精英策略的遺傳算法的參數(shù)設(shè)置如下：種群規(guī)模popsize=400,最大迭代次數(shù)maxgen400，交叉概率pc=0.92，變異概率為pm=0.01。

4.2　運算結(jié)果分析

通過Matlab2017a編程，在2.9 GHz Intel Core i5，8 GB內(nèi)存，macOS操作系統(tǒng)的計算機上進行運算。算法終止時，所有解收斂到Pareto前沿面上，得到140個Pareto解，Pareto最優(yōu)解分布見圖3。表3是從中抽取的5組(A-E)典型Pareto解。

4.3　結(jié)果對比與分析

為了保證算法運算的穩(wěn)定性與可操作性，文章根據(jù)不同種群數(shù)量和算法迭代次數(shù)設(shè)計了4組參照組。表4是參照組1、2、3、4各主要指標(biāo)的具體數(shù)值以及不同種群數(shù)量和迭代次數(shù)下的解數(shù)量和運算時間，各參照組Pareto前沿面結(jié)果如圖4。

表3　典型Pareto解

表4　各參照組主要指標(biāo)及運行時長

通過對結(jié)果的分析我們可以得出以下論：

(1)通過帶精英策略的非支配排序的遺傳算法的運算，我們得出不同參數(shù)下的Pareto前沿面解集，證明了模型和算法的有效性。

(2)Pareto解的運算時間與種群數(shù)量和迭代次數(shù)有關(guān)。伴隨著種群數(shù)量和迭代次數(shù)的不斷增加，Pareto解的數(shù)量呈現(xiàn)著遞增趨勢，通過參照組的對比分析，種群規(guī)模越大，迭代次數(shù)越多，得到的結(jié)果就更加豐富，分布也更加均勻。

(3)通過表4發(fā)現(xiàn)，在迭代次數(shù)相同的情況下，種群數(shù)量越多，Pareto解的集合也就越多，但消耗更多的運算時間；在種群規(guī)模相同的情況下，迭代次數(shù)越多，Pareto解的數(shù)量更加豐富，同時消耗更多的運算時間。Pareto解的數(shù)量和種群規(guī)模以及迭代次數(shù)在一方面條件確定的情況下存在著正相關(guān)的關(guān)系。

圖4　各參照組Pareto前沿面解的分布圖

通過圖4三個參照組的對比可以發(fā)現(xiàn)，公平性目標(biāo)與時效性目標(biāo)存在著悖反現(xiàn)象。從Pareto前沿面的參照圖來看，救災(zāi)物資的時效性目標(biāo)波動較小，其分配物資的公平性目標(biāo)波動較大。從圖3可以看出，當(dāng)公平性目標(biāo)越大，時效性目標(biāo)略微增加時，公平性目標(biāo)則大幅度降低。同時災(zāi)區(qū)需求的增加會使公平性目標(biāo)受到損失。

在救援的初級階段，決策者更多的關(guān)注時效性目標(biāo)，在救援的后期階段，決策者更多的關(guān)注公平性目標(biāo)。通過Pareto前沿面的分析，為救援決策者根據(jù)在不同階段對決策目標(biāo)的偏好提供了更多的選擇。

5　結(jié)語

震后應(yīng)急物資的分配與運輸是震后救援保障工作的重要環(huán)節(jié)，配送中心的選址的穩(wěn)定性與物資保障效率的提升對于減少財產(chǎn)損失、保證生命安全有著重要意義。震后物資運輸與配送過程包括將應(yīng)急物資從區(qū)域儲備倉庫和外圍物資集散地運輸?shù)絽^(qū)域應(yīng)急配送中心、區(qū)域配送中心配送到各個受災(zāi)點兩個階段。但是應(yīng)急物資保障具有持續(xù)時間長、需求緊迫性強和物資相對短缺的特點，因此將物資分配時效性與公平性，多時段、多種運輸方式共同參與的應(yīng)急物資保障策略有機的結(jié)合起來，綜合考慮震后應(yīng)急物資的運輸與分配特性，建立了一個考慮時效性-公平性的震后配送中心選址和物資動態(tài)分配的多目標(biāo)規(guī)劃模型。基于模型設(shè)計了一種帶精英策略的非支配排序的遺傳算法，通過非線性混合整數(shù)規(guī)劃結(jié)合Pareto前沿面對時效性-公平性雙目標(biāo)動態(tài)救援模型進行求解。通過算法對實際案例的求解分析，發(fā)現(xiàn)模型可以產(chǎn)生有效的Pareto最優(yōu)解和前沿面，準(zhǔn)確反映了時效性目標(biāo)和公平性目標(biāo)的悖反現(xiàn)象，通過數(shù)學(xué)模型的計算將時效性與公平性的關(guān)系通過數(shù)學(xué)結(jié)果對模型和算法的有效性進行了驗證，為今后進一步研究該問題奠定了基礎(chǔ)。

在構(gòu)建的模型中，救援過程涉及的相關(guān)系數(shù)的是確定的數(shù)值，但是在實際的運輸過程中是隨機變化的。如道路安全系數(shù)、車輛的安全性能、運輸時間受道路清障的影響等因素，這些帶有隨機參數(shù)的目標(biāo)函數(shù)模型是未來研究的方向。