楊慧慧，鄭義彬

0 引言

隨著新冠疫情的沖擊以及“互聯(lián)網(wǎng)+醫(yī)療”模式的不斷推進(jìn)，人們對于藥品需求的及時性、多樣性及個性化需求不斷增加，傳統(tǒng)的醫(yī)藥物流模式已不能滿足客戶的體驗以及速度要求。近年來興起的藥品O2O模式具有將線上和線下特點完美結(jié)合的優(yōu)勢，不僅能夠與消費者近距離接觸，增加消費者信任度，也能以更快的速度響應(yīng)消費者需求。但是由于我國線下零售藥店龐大、藥店的面積有限以及國家藥品的配送具有嚴(yán)格的要求，使得藥品O2O企業(yè)在未能實際滿足消費者需求的情況下，運營成本還一直居高不下。如何才能在保持較低成本下給消費者提供更加便捷、高效的物流服務(wù)，成為醫(yī)藥電商企業(yè)急需解決的問題。

國內(nèi)外關(guān)于藥品物流配送系統(tǒng)優(yōu)化的研究起步較早，目前已有較為成熟的研究成果，且主要集中于對配送系統(tǒng)庫存、選址、配送路徑等方面的研究。劉聰[1]應(yīng)用雙層規(guī)劃模型研究了O2O模式下藥品的配送選址路徑問題，并采用遺傳算法進(jìn)行求解；劉閃[2]針對醫(yī)藥產(chǎn)品所具有的特殊特點，研究了具有時間窗限制的多車型指派和配送路徑集成優(yōu)化問題；Cho，等[3]研究了制藥物流公司在進(jìn)行物流中心選址時應(yīng)當(dāng)考慮的因素，研究結(jié)果可以為藥品分銷公司選擇物流中心和制定相關(guān)戰(zhàn)略提供指導(dǎo)；羅威，等[4]使用改進(jìn)的Dijkstra算法對藥品物流的配送路線進(jìn)行了優(yōu)化研究；胡劍玫[5]根據(jù)普通藥品和冷藏藥品不同的配送特點，分別建立了帶時間窗的醫(yī)藥物流配送模型，研究了其配送路徑的優(yōu)化問題；趙蒙，等[6]采用時空網(wǎng)絡(luò)建模思想，研究了新零售背景下藥品的協(xié)同配送優(yōu)化問題?？偨Y(jié)可得：目前有關(guān)電子商務(wù)模式藥品配送的相關(guān)研究較為匱乏，且大多只集中于選址、庫存、配送等單一方面的研究，多層次集成決策的研究較少。

因此，筆者從藥品O2O企業(yè)在實際運營中存在的問題出發(fā)，研究O2O模式下藥品的選址-庫存-配送問題，將藥品O2O配送系統(tǒng)中涉及到的選址、庫存、配送問題看作一個整體，考慮以集成優(yōu)化的思想對其進(jìn)行優(yōu)化，以達(dá)到藥品O2O配送系統(tǒng)成本最優(yōu)、顧客滿意度最高的目的。

1 問題描述與模型構(gòu)建

1.1 問題描述

藥品O2O配送系統(tǒng)優(yōu)化是指藥品配送系統(tǒng)中“區(qū)域配送中心—藥店—消費者”環(huán)節(jié)的配送優(yōu)化。本文研究醫(yī)藥O2O企業(yè)在某個城市的藥品O2O配送系統(tǒng)優(yōu)化，研究范圍包括藥品從區(qū)域配送中心到藥店再到消費者的環(huán)節(jié)，研究對象包括一個區(qū)域配送中心、多個潛在的藥店以及眾多潛在的消費者點，網(wǎng)絡(luò)結(jié)構(gòu)圖如圖1所示。O2O模式下藥品的選址—路徑—庫存集成優(yōu)化問題可以描述為：隨機(jī)需求下帶軟時間窗的藥品O2O配送系統(tǒng)集成優(yōu)化，優(yōu)化目標(biāo)為在保證配送系統(tǒng)總成本最小、消費者滿意度最高的情況下，確定開放藥店的位置及數(shù)量，確定每個藥店需要服務(wù)的消費者及產(chǎn)生的配送線路，以及確定藥店的最優(yōu)庫存管理策略。

圖1 藥品O2O配送系統(tǒng)優(yōu)化網(wǎng)絡(luò)結(jié)構(gòu)圖

1.2 模型假設(shè)

為了便于研究，本文作出相關(guān)合理的假設(shè)條件如下：

（1）在藥店的運營中，其銷售的藥品復(fù)雜多樣，但大多屬于標(biāo)準(zhǔn)盒裝或者瓶裝，所以本文將配送的藥品都抽象為同類藥品。

（2）藥品O2O企業(yè)采用的是自建藥房的形式。

（3）由于藥店所服務(wù)的消費者屬于分布較為分散的個體消費者，且對于藥品的需求量較少，所以本文可將某個區(qū)域內(nèi)地理位置較近的所有消費者看作為一個消費者需求點，這樣可以在一定程度上保證每天每個消費者需求點都存在用藥需求。

（4）由于消費者對于藥品的需求量、需求時間具有隨機(jī)不確定性，消費者在一天24h內(nèi)都可能存在用藥需求，為了便于研究，本文在計算時僅選取某個時段進(jìn)行計算。

（5）藥店的配送模式為“自營+第三方物流”混合配送模式，當(dāng)藥店的物流團(tuán)隊能夠滿足藥店的訂單量需求時，僅采用自營物流進(jìn)行配送；當(dāng)藥店的訂單量超過其自營物流所能承擔(dān)的最大能力時，將超出的訂單量交由第三方即時物流進(jìn)行配送。

（6）藥店在為消費者進(jìn)行配送時，每條配送線路上僅有一輛配送車輛為線路上的消費者提供服務(wù)，且配送車輛的配送路徑應(yīng)以藥店為起點，經(jīng)過其所需要進(jìn)行配送服務(wù)的消費者后，最后返回到藥店，即形成一個閉回路的行駛路徑。

（7）由于藥品O2O的配送為“最后一公里”的末端配送，所以假設(shè)藥店到消費者的距離都在藥店的配送范圍內(nèi)。

（8）本文所涉及的配送車輛均為同種類型的車輛，即車輛的運載能力相同；在進(jìn)行配送時每條配送路徑上的消費者需求量之和不能大于車輛的最大運載量。

（9）假設(shè)藥品O2O企業(yè)的區(qū)域配送中心可以滿足其負(fù)責(zé)區(qū)域所有藥店的補貨需求，且為藥店補貨時的訂貨提前期一致。

（10）由于大部分藥店都屬于全年營業(yè)狀態(tài)，且本文研究的O2O模式下的藥店更是屬于24h不間斷營業(yè)的藥店，故本文涉及藥店的工作天數(shù)按365d記。

（11）假設(shè)藥店采用連續(xù)型檢查庫存策略，且所有藥店設(shè)立的庫存服務(wù)水平均相同。

（12）由于消費者對于藥品的需求具有隨機(jī)性，本文假設(shè)消費者每日對于藥品的需求相互獨立且服從正態(tài)分布。

1.3 參數(shù)和符號

將模型中用到的相關(guān)變量及參數(shù)進(jìn)行定義如下：

I：所有確認(rèn)開放藥店的集合，I={i |i=1,2,...,m}；

J：所有潛在消費者需求點的集合，J={ j|j=1,2,...,n}；

V：所有從藥店到消費者需求點的配送車輛的集合，V={v |v=1,2,...,k}；

F：藥品O2O配送系統(tǒng)的總成本；

Fi：開放藥店i的年運營成本，包括藥店的建設(shè)成本和實際運營成本；

f1：藥店的總選址成本；

f2：藥店的總庫存成本；

ri：藥店i的再訂貨點的實際庫存量；

Qi：藥店i的固定訂貨批量；

L：藥店i的訂貨提前期，即藥店向區(qū)域配送中心發(fā)起兩次相鄰訂單之間的間隔期；

Ti：藥店i的補貨周期；

ssi：藥店i的安全庫存量；

RCi：藥店i的年補貨成本；

dj：客戶點j的每日需求期望值；

Di：藥 店i的 每 日 需 求 期 望 值，且

C1：區(qū)域配送中心到藥店的單位配送成本；

Si：藥店i一個訂貨周期內(nèi)的平均缺貨量；

SCi：藥店i的年平均缺貨成本；

C2：藥店的單位缺貨成本；

θ：藥品的變質(zhì)率；

Ii(t)：藥店i在t時刻的庫存水平；

HCi：藥店i的年庫存持有成本；

C3：單位藥品的日庫存持有成本；

α：藥店在補貨提前期內(nèi)的庫存服務(wù)水平；

zα：藥店服務(wù)水平為α?xí)r藥店的安全庫存系數(shù)；

qm：藥店一個訂貨周期內(nèi)的變質(zhì)藥品數(shù)量；

MCi：藥店i的年過期成本；

C4：藥品的單位變質(zhì)成本；

f3：藥品從藥店到達(dá)消費者需求點的配送成本；

C5：配送藥品車輛單位距離的運輸成本；

C6：配送車輛的啟用成本，包括使用該車輛進(jìn)行配送的人員成本以及車輛的使用成本；

disabv：配送節(jié)點a至配送節(jié)點b之間的距離；

f4：藥品在配送過程中，將基于時間窗的顧客滿意度轉(zhuǎn)化成的時間懲罰成本；

C7：采用第三方即時物流配送時，配送車輛的單位運輸成本；

C8：當(dāng)配送車輛到達(dá)消費者需求點的時間早于消費者能夠接受的最早服務(wù)時間時，單位等待時間的損失成本；

C9：當(dāng)配送車輛到達(dá)消費者需求點的時間晚于消費者能接受的最晚服務(wù)時間時，單位遲到時間的懲罰成本；

gj：車輛到達(dá)消費者需求點j的時間；

LTj：消費者j可接受車輛提供配送服務(wù)的最晚時間；

ETj：消費者j可接受車輛提供配送服務(wù)的最早時間；

CQ：藥店庫存的最大容量；

Cq：為藥店提供配送服務(wù)的車輛的最大運載能力；

ni：藥店i啟用的自營配送車輛數(shù)量；

wa：配送車輛在配送節(jié)點的服務(wù)持續(xù)時間；

s：配送車輛單位距離的旅行時間。

1.4 模型構(gòu)建

本文的目標(biāo)是藥品O2O配送系統(tǒng)的總成本最小化，故構(gòu)建如下數(shù)學(xué)模型：

其中：

約束條件：

其中：式（1）為模型的目標(biāo)函數(shù)，表示由藥店的選址成本式（2）、藥店的庫存成本式（3）、藥店到消費者需求點的配送成本式（4）以及藥店為消費者服務(wù)的時間懲罰成本式（5）組成的配送系統(tǒng)總成本最小化；式（6）表示消費者需求點與藥店、配送車輛的關(guān)系為多對一關(guān)系；式（7）表示一輛配送車輛最多只能為一個藥店提供服務(wù)；式（8）表示配送車輛的行駛路徑應(yīng)為一個巡回閉回路，且應(yīng)依次訪問各消費者需求點；式（9）表示在所有的備選藥店開放點中，只有最終確定開放的藥店才能為消費者提供服務(wù)；式（10）表示開放的藥店具有庫存容量限制，藥店所服務(wù)的所有消費者對于藥品的總需求量不能超過藥店的最大庫存容量；式（11）表示配送車輛單次服務(wù)的每條配送線路上的所有消費者需求量之和不能超過車輛的最大運載量；式（12）表示藥店自營物流的服務(wù)能力限制；式（13）表示配送車輛在配送點a的等待服務(wù)時間；式（14）表示配送車輛到達(dá)配送節(jié)點b的時間與到達(dá)其前一個配送節(jié)點a的數(shù)量關(guān)系；式（15）-式（18）為決策變量，分別表示車輛v是否從配送節(jié)點a經(jīng)過配送節(jié)點b，藥店i是否采用第三方即時物流進(jìn)行配送，藥店備選點i是否開放，藥店備選點i是否為消費者需求點j服務(wù)。

2 模型求解

2.1 算法設(shè)計

本文采用遺傳算法對提出的藥品選址-庫存-路徑集成優(yōu)化模型進(jìn)行求解，其主要思想在于對參數(shù)模型使用特殊編碼方法，對隨機(jī)生成的初始解進(jìn)行選擇、變異、交叉等操作后，最終得出優(yōu)化和改進(jìn)后的最優(yōu)解。

（1）染色體編碼。為降低算法搜索復(fù)雜度，避免冗長的計算過程，將求解問題通過編碼轉(zhuǎn)換，以達(dá)到算法可執(zhí)行搜索的搜索空間，這個過程稱為染色體編碼。根據(jù)所求解問題的不同情況，一般選用更適合的編碼方式進(jìn)行編碼，二進(jìn)制編碼的編碼能力強，但交叉變異的種群多樣性較低，浮點數(shù)編碼的編碼能力比二進(jìn)制編碼弱，但交叉變異的種群多樣性較高。

通過結(jié)合O2O模式下藥品的選址-庫存-路徑集成優(yōu)化的特點，本文選擇采用浮點數(shù)編碼的方式對藥店備選點、消費者需求點進(jìn)行編碼，即用一維數(shù)組儲存藥店備選點i與消費者需求點j，其中i與j均用一定范圍內(nèi)的浮點數(shù)表示，其中的基因取值即為該位置號對應(yīng)的任務(wù)分配到計算資源上的資源編號。假設(shè)有10個消費者需求點，3個藥店備選點，即染色體的長度為13，將藥店備選點、消費者需求點的取值用1到2之間的整數(shù)表示，若產(chǎn)生的染色體為{1，2，2，1，2，1，2，1，1，2}，通過解碼可以知道1、4、6、8、9的消費者由選擇出開放的藥店一為其服務(wù)，2、3、5、7、10的消費者由選擇出開放的藥店二為其服務(wù)。通過對染色體解碼，得到開放藥店的服務(wù)范圍分布，得到每個藥店開放需要服務(wù)的消費者需求點序列A分別為：A1{1、4、6、8、9}、A2{2、3、5、7、10}。

（2）初始化種群。本文種群規(guī)模設(shè)為Sizepop，藥店備選點的數(shù)量為m，消費者需求點的數(shù)量為n，則染色體長度為n,每個染色體的每個基因位值采用隨機(jī)生成的方式從{1，2，…，m}集合中取值。

（3）適應(yīng)度函數(shù)。適應(yīng)度函數(shù)就是在遺傳算法中對種群進(jìn)行適應(yīng)度評判的函數(shù)，它在某種程度上可以說是算法的核心要素，使得每一次迭代的結(jié)果都向著更優(yōu)的方向發(fā)展。本文中，在迭代過程中解的優(yōu)劣也是通過轉(zhuǎn)化后的適應(yīng)度函數(shù)來評判是否滿足最優(yōu)。

（4）選擇操作。選擇操作的核心在于對個體進(jìn)行選擇，適應(yīng)條件的個體存活，否則淘汰。本文所采用的方法是將選中個體的概率與適應(yīng)度函數(shù)相關(guān)聯(lián)，即輪盤賭選擇法。

（5）交叉、變異操作。遺傳算法中的交叉與變異即生物學(xué)里的基因重組、基因變異，是算法獲取新的優(yōu)良個體最重要的手段，決定了遺傳算法的全局搜索能力。本文采用算數(shù)交叉的方法對染色體進(jìn)行交叉操作。即經(jīng)過交叉操作產(chǎn)生的染色體是由兩個個體的線性組合產(chǎn)生的。變異操作是針對單個染色體，根據(jù)指定的變異率改變單個染色體編碼中某個或某些位置上的基因，用新的基因值替換原來的基因值，生成新的染色體。本文采用較小的步長進(jìn)行變異。變異算子為：M’=M±1，式中M為變異前的變量取值，M’為變異后的變量取值，M和M’要在變量的取值范圍內(nèi)。求解思路圖如圖2所示。

圖2 求解算法的基本思路

2.2 算法流程

綜上，本文設(shè)計的藥品O2O的LIRP遺傳算法流程如下：

（1）設(shè)定種群容量、變異交叉概率等預(yù)置參數(shù)值。

（2）采用浮點數(shù)編碼方法進(jìn)行編碼。

（3）通過在所設(shè)集合中生成偽隨機(jī)數(shù)的方式初始化種群。

（4）首先計算各部分成本的染色體適應(yīng)度函數(shù)值，然后按數(shù)值對其進(jìn)行由小到大的排序，最后采用輪盤賭選擇法對其進(jìn)行評判篩選，保留較好的前S個染色體，其余刪除。

（5）使用前文所述的算術(shù)交叉法對所選元素個體進(jìn)行交叉變異操作，產(chǎn)生新生群體，以此為基礎(chǔ)計算群體中個體的適應(yīng)度函數(shù)值，重復(fù)步驟（4）中排序取舍、輪盤賭優(yōu)化的操作，最終得到當(dāng)前最優(yōu)解。

（6）與預(yù)先設(shè)定好的迭代次數(shù)對比，若當(dāng)前次數(shù)等于迭代次數(shù)，則終止計算，輸出結(jié)果；否則，重新轉(zhuǎn)入步驟（4）。具體流程圖如圖3所示。

圖3 隨機(jī)需求下帶軟時間窗的藥品O2O的LIRP遺傳算法流程

3 算例驗證

3.1 算例描述

某藥品O2O企業(yè)計劃在北京市豐臺區(qū)試點布局零售藥店，為消費者提供配送服務(wù)，基于豐臺區(qū)社區(qū)分布和周邊環(huán)境的調(diào)查分析，本文確定了藥店的10個備選點，將距離備選藥店點較近的社區(qū)看作是一個消費者需求點，得到藥店的備選點及消費者需求點分布情況，如圖4所示。本文使用excel取0~5km隨機(jī)數(shù)得出備選藥店與消費者需求點的最近公路距離。

圖4 藥店備選點與消費者需求點分布情況

其他相關(guān)參數(shù)設(shè)置如下：

假設(shè)藥品O2O企業(yè)華北地區(qū)的區(qū)域配送中心距離北京市零售藥店的距離為60km，藥店的最大容量為2 000盒，區(qū)域配送中心到藥店配送車輛最大載貨量為6 000盒藥品，整車運輸運費一般按50元/km計費，假設(shè)每次運輸車輛都是滿載，從區(qū)域配送中心到北京市的運費為0.5元/盒。藥店藥品的變質(zhì)率為3%，藥店單位產(chǎn)品每天庫存持有成本2元/盒，單位變質(zhì)成本為8元/盒，單位缺貨成本為3元/盒，庫存服務(wù)水平為90%，訂貨提前期為1天。

藥店的配送人員對于用戶日需求量，每天需要均分為四次送完，算例僅計算一天中18點到21點時段內(nèi)的配送，且在向消費者進(jìn)行送藥時，遲到懲罰成本為10元/h，等待成本為1元/h，假設(shè)將一小時分成12個單位，每個單位時間的懲罰成本為5/12元，每個單位時間的等待成本為1/12元，每個藥店配備的配送車輛數(shù)最多為2輛，若訂單量超過2輛車所能配送的限額，則剩下的訂單將交由第三方即時物流進(jìn)行配送；配送車輛的平均配送速度為20km/h，最大運輸量為150盒。采用自營物流配送時，每輛車的啟用成本為50 200元/年。采用第三方即時物流進(jìn)行配送時，車輛單位距離配送成本為10元。藥店的成本見表1。其他有關(guān)消費者需求點的參數(shù)取值見表2。

表1 藥店備選點成本參數(shù)（單位：元）

表2 消費者需求點參數(shù)

3.2 算例求解

使用遺傳算法的Matlab程序?qū)哂须S機(jī)需求和軟時間窗的藥品O2O配送系統(tǒng)集成優(yōu)化問題進(jìn)行求解。本文遺傳算法的運行參數(shù)設(shè)置為：染色體基因位數(shù)為30，種群規(guī)模為10，迭代次數(shù)為1 000，交叉概率為0.5，變異概率為0.2。算法的最優(yōu)解隨迭代次數(shù)的變化趨勢如圖5所示。由圖5可知，其目標(biāo)函數(shù)隨迭代次數(shù)的變化在前200代表現(xiàn)得非常劇烈，在200代后搜索過程中變化趨于平穩(wěn)。根據(jù)結(jié)果可知，優(yōu)化模型和求解算法都是可行的，能夠達(dá)到配送系統(tǒng)優(yōu)化的目的。

圖5 最優(yōu)解隨遺傳代數(shù)變化的趨勢

考慮隨機(jī)需求下帶客戶服務(wù)時間窗的配送系統(tǒng)集成優(yōu)化模型的優(yōu)化方案中應(yīng)包含最優(yōu)總成本、優(yōu)化的藥店選址-分配策略和最優(yōu)配送路徑以及最優(yōu)訂貨量和再訂貨點，集成優(yōu)化的最優(yōu)方案見表3，最優(yōu)方案對應(yīng)的各項成本見表4。

表3 藥品O2O配送系統(tǒng)模型的最優(yōu)解

本案例的最優(yōu)方案為：從10個藥店備選點中選擇出2個藥店進(jìn)行開放，即9號倉、10號倉，由表4可以看出其選址成本相對較低，且9號倉和10號倉同樣都只啟用了2輛自營的配送車輛，各自生成2條巡回的配送路線；1號倉的再訂貨點為988盒，訂貨批量為1 012盒；2號倉的再訂貨點為1 002盒，訂貨批量為998盒。

由表4可知，在藥品O2O配送系統(tǒng)的各項成本中，庫存成本所占比例最大，約占總成本的49%，選址成本占總成本比例為30%，配送成本占總成本比例為21%。

表4 最優(yōu)方案對應(yīng)的各項成本（單位：元）

4 結(jié)語

本文從O2O模式下藥品配送系統(tǒng)現(xiàn)狀及問題出發(fā)，考慮以多層次的思想對藥品O2O的選址-庫存-配送進(jìn)行優(yōu)化，構(gòu)建了隨機(jī)需求下帶軟時間窗的藥品O2O配送系統(tǒng)集成優(yōu)化模型，并設(shè)計遺傳算法對模型進(jìn)行求解，最后通過算例分析驗證模型的有效性。該研究為藥品O2O模式配送在我國藥品物流配送中的推廣應(yīng)用提供了切實可行的優(yōu)化方案，具有一定的理論價值和實踐意義。