1、第 卷 第 期 年 月公 路 交 通 科 技 .收稿日期:基金项目:北京市社会科学基金规划项目()。作者简介:赵银婷(),女,河南舞阳人,硕士研究生.(.):.疫情暴发初期应急医疗物资配送优化研究赵银婷,徐晓敏,周 行(.北京信息科技大学 信息管理学院,北京;.中国科学院大学 电子与电气工程学院,北京)摘要:疫情暴发初期,由于病毒传播性和变异性强、传播速度快、涉及范围广等特点给疫情防控带来了巨大困难和挑战,而应急医疗物资能否及时保障供应直接关乎病患救治率和死亡率,是保障人民生命安全的重要因素。因此,充分考虑疫情暴发初期特点,有针对性地研究应急医疗物资的快速精准供给,具有十分重要的现实意义和实际
2、应用价值。以疫情暴发初期医疗物资需求特点为出发点,对疫情发生地(疫灾点)应急医疗物资的精准匹配、及时运送的优化进行了探索。首先,以各疫灾点总的等待时间最短为目标函数,综合考虑配送车型的多样性、医疗物资多样性、供需不平衡等约束,构建了疫情地储备中心到多个疫灾点间的车辆路径集成优化模型。其次,通过改进遗传算法予以模型实现,通过增加信息素浓度提升系数对信息素浓度的更新规则进行了升级,降低了蚁群的搜索空间且提高了算法的收敛速度和稳定性。最后,通过仿真算例对上述模型及算法进行了验证,将其优化结果与改进遗传算法求解结果进行了对比分析,验证了改进蚁群算法求解的适应性和有效性。研究较好地解决了疫情暴发初期应急
3、医疗物资“最后一公里”的优化配送问题,可有效提高应急医疗物资的配送效率,为应对重大公共卫生事件提供一定的决策参考。关键词:物流工程;车辆路径综合优化;改进蚁群算法;医疗物资;应急管理中图分类号:.,文献标识码:文章编号:(),(.,;.,):,.,.,.,().,公路交通科技第 卷.,.,“”,.:;();引言 年末暴发的新冠肺炎疫情对国民经济和人民的安全造成了极大冲击。重大疫情事件往往具有在短时间内传染速度快、传播不确定性强等特点。当疫情发生时,患者救治、医护人员安全、疫灾点卫生防疫等问题需要在第一时间得到处理,会对应急医疗物资产生很大的需求。因此医疗物资的及时配送就成为应对疫情重大公共卫生
4、事件的关键问题。近年来,国内外学者对应急医疗物资配送问题的探索主要集中于以下几个方面。在单目标建模优化进行应急物资配送的研究中主要有卢安文,罗平等学者,这些学者基于成本因素视角考虑运输费用及多车型低耗车路径问题进行了应急物资配送的研究。,徐卫,等侧重考虑时间的影响,基于最短运达时间建立模型对应急物资的配送展开了探究。还有一些学者则考虑公平原则,基于公平视角构建了应急医疗物资模型。在多目标、多模型建模优化进行应急物资配送的研究中,一些学者结合效率与公平及未满足需求量和总物资延误时间因素等建立了多目标模型,且考虑物资的分类分批不同周期等来构建多模型进行研究。综合以上研究成果可以看出,目前专门针对疫
5、情应急医疗物资配送优化问题的研究相对较少,多数研究着眼于通用应急物资的配送优化,而没有考虑突发公共卫生事件所具有的特殊场景。仅仅局限在单一应急物资或单一车型的配送研究,没有集成性考虑多种物资、多类车型的综合性优化。疫情应急物流具有比较鲜明的特点,如疫情初期对医疗物资供不应求、配送时间紧急性、物资需求点的变动性等不确定性因素更多,情况十分复杂,因此需要针对疫情应急事件物资配送的特点进行专门研究。鉴于此,本研究以我国应对新冠肺炎重大疫情事件为背景,深入分析疫情中应急医疗物资配送中存在的问题,结合我国国情,考虑疫情发生时疫灾地应急物资储备中心与当地医疗物资需求点之间,以需求点等待时间最短为目标、构建
6、车辆路径的综合优化模型,并应用改进蚁群算法实现对模型的求解,得出最佳配送方案。模型构建.研究背景本研究主要针对重大疫情应急医疗物资配送优化问题专门进行。应急医疗物资是指突发性自然灾害、公共卫生事件等发生后,用于救治伤员、患者及保护医护人员生命安全所需的一类专用物资,该类物资一般包括:各类医疗设备、医疗药品及医用防护用品等。应急医疗物资具有不确定性、不可替代性、时效性及滞后性等特点。应急医疗物资的及时保障对人民生命安全极为重要,因此应急医疗物资的配送问题成为突发事件应急管理的重中之重。应急医疗物资保障一般包括物资筹措、物资调度、物资配送等环节,不同的突发事件其物资保障方式有所不同,其中尤以疫情医
7、疗物资具有较大的特殊性,如:疫情发生时现有交通网络、道路设置完好,但疫灾点位置、数量及疫灾点对医疗物资需求量和需求类型随着时间呈现动态变化。基于此,本研究针对疫区应急医疗物资配送优化问题进行深入探究。具体地,通过收集各疫灾点所需各类医疗物资信息,构建从疫情地应急物资储备中心到各疫灾点的总等待时间最短的数学模型,对模型进行求解得出优化配送方案,将应急储备中心内有限的医疗物资在最短的时间内快速、精准送达至各个疫灾点。第 期赵银婷,等:疫情暴发初期应急医疗物资配送优化研究.模型假设假设:疫情发生地建有 个应急医疗物资储备中心,记为 点,有 个疫灾医疗物资需求点。假设:应急医疗物资储备中心的位置及医疗
8、物资类型、供应量、疫灾点位置、医疗物资类型和需求量已知,储备中心与各疫灾点及疫灾点与疫灾点之间的线路连通且距离已知。假设:疫情发生时,当地应急管理部门政府临时征用各类社会车辆(如公交车、出租车、私家车等)进行统一调度。车辆类型各异,运输能力差别较大。配送车辆有其不同的核定载重及不同的运载体积限制,所有车辆从应急物资储备中心装载医疗物资出发按照优化方案将指定物资运送至指定的疫灾点。疫情应急医疗物资配送网络模型见图。图 应急医疗物资配送网络模型.模型建立目标函数:以各疫灾点总的等待时间最短来构建数学模型,寻求出最优的配送路径方案,同时找到各疫灾点的总的最短等待时间,即:,()式中,为疫灾点节点集合
9、,;为从应急储备中心出发到疫灾点 点间的距离,(,);(,)为 点到 点的距离;为向量 中的元素;为车辆行驶速度。约束条件 应急储备中心医疗物资供应量和各疫灾点需求量之间的约束。考虑疫情初期医疗物资供不应求的实际情况,则各疫灾点的总需求量应大于应急医疗物资的总供应量,即:,()式中,为医疗物资的类型;为医疗物资的类型总数;为第 种医疗物资的供应量;为第 种医疗物资的需求量。约束条件 各车型的可载医疗物资总体积不能超过各车型的核定载容量,即:,()式中,为配送车辆的类型;为车辆的类型总数;为第 种医疗物资的单位体积;为第 种类型车的载容量;为第 种类型车的数量;为车辆在实际配送中,第 类车辆所装
10、第 类医疗物资的量。约束条件 各车型所装载各种物资的总重不能超过其核定载重,即:,()式中,为第 种医疗物资的单位质量;为第 种类型车的载重量。约束条件 可供应的所有应急医疗物资总体积不得超过所有运载车辆总的载容量,即:。()约束条件 可供应的所有应急医疗物资总质量不得超过所有运载车辆总的载重,即:,()优化实现上述应急医疗物资配送优化问题因其包括车辆路径的集成优化,模型变量及约束条件较多,求解比较复杂。针对该模型的这些特点,选择具有较强鲁棒性、并行性、稳定性及能够快速求解的蚁群算法比较合适。蚁群算法(,)是由意大利学者 等在 年创造的一种随机搜索算法。该算法的理念来自于自然界中蚁群寻找食物的
11、现象。许多研究表明蚁群算法公路交通科技第 卷在解决物资分配、路径优化等问题时表现出良好性能。假设蚁群中共有 只蚂蚁,蚂蚁代表配送的车辆,共有 个疫灾点,为所有疫灾点的点集,设()为 时刻在疫灾点 处蚂蚁的数量,则 ()为蚂蚁总量,()为在 时刻路径 上积累的信息素量,为信息素挥发因子。但在应用蚁群算法优化求解过程中发现,基本蚁群算法求解的模型多为任一只蚂蚁需走完全局路径的情况,而本问题的优化求解要求蚂蚁所代表的车辆只需负责个别疫灾点的配送,用基本蚁群算法求解该模型时易局部最优,使得模型求解结果不收敛。因此,需要根据问题的特点对基本蚁群算法进行相应改进,具体改进如下。.信息素更新规则在每次寻找路
12、线的迭代中,每只蚂蚁代表的车辆都从应急物资储备中心出发开始访问,在访问任意 个疫灾点后都有一定的概率决定是否继续配送后面的疫灾点,而一旦这只蚂蚁决定配送下一个疫灾点,则当前疫灾点完成配送,并将之放到禁忌表的点集中。经过 次迭代,可以得到各车辆所配送分路径,而蚂蚁代表的车辆在每次运动过程中都会释放信息素到所走的路径上。到疫灾点各路线上信息素的更新规则为:()()()|,()式中,信息素浓度 为当前配送车辆在对疫灾点 配送完成后,下一个疫灾点 对该车辆产生的吸引程度;为信息素强度(此处);为本次迭代中各疫灾点总的等待时间。疫灾点总的等待时间对信息素浓度的增量有很大的影响,总等待时间越短,则信息素浓
13、度的增量则会越大。在算法的设计中,一方面维护储备中心到每个疫灾点的最短等待时间矩阵,每次迭代观察是否找到的到 点的支路上等待时间比原来的短,若不是,则令信息素浓度为,限制后面迭代中蚂蚁再走这条支路到 点。若找到的当前到 点的支路比原来短,则将新的信息素浓度转移到新找到的支路上,并且更新到每个疫灾点的最短等待时间矩阵。另一方面,增加了信息素浓度提升系数,|是自适应变化的,随着蚂蚁本次行走支路径的时间而进行动态调整。在某次迭代中,若到疫灾点 有更短的配送路径时(即,为初始迭代的到 点最短的路径的时间;为本次迭代中从 点到 点所在路径的时间)有信息素浓度自适应提升系数 与原结合对新配送路径的信息素浓
14、度进行更新,即增加 信息素浓度。信息素浓度提升系数对路径上信息素浓度 进行动态提升,主要依据新支路径上的行驶时间与原最短路径行驶时间进行比较,若该路径行驶时间小于历史路径行驶的最短时间,信息素浓度的自适应提升系数加强当前支路径上信息素的增量(即正反馈),且每条路径都会根据自身值得优劣来获得相应得信息素增量。.状态转移规则车辆在选择配送疫灾点的过程中,各分路径线上的信息素浓度同时会影响车辆下一步的转移方向。随着路径上不断积累新的信息素,车辆在配送的过程中会根据各路径上信息素的强度及各路径的启发式信息决定下一步行动方向,即状态转移概率。上面对信息素改进中增加了提升信息素浓度系数,因此在蚂蚁由疫灾点
15、向另一个疫灾点的转移概率也会随之变化。其中,蚂蚁 由当前疫灾点 到可访问的疫灾点集合 中的状态转移概率为:()()()()若 其他|,(),()式中(),为蚂蚁 当前可选择配送的所有疫灾点集合;为疫灾点 到 上的信息素浓度;为信息素重要程度因子;为启发函数重要程度因子;为待配送疫灾点 到 的距离;为蚂蚁从疫灾点 到 的期望程度。当蚂蚁在某疫灾点完成配送,则该疫灾点以后将不能被访问。通过上述改进,克服了标准蚁群算法在搜索空间和时间性能上的不协调,较好地实现了每只蚂蚁(所代表相应的车辆)找到所负责的 个或多个疫灾点的 条路径。将所有蚂蚁找到的路径组合在一起得到完整疫灾点的配送路线,继而得到各个疫灾
16、点总的最短等待时间。基于改进蚁群算法的疫情应急物资配送优化求解步骤为:()初始化各参数,将 设定为迭代次数,开始时间,设定最大的迭代次数。将 只蚂蚁放在应急物资储备中心,初始时刻()。第 期赵银婷,等:疫情暴发初期应急医疗物资配送优化研究()建立禁忌表,根据状态转移式(),确定下一个配送的疫灾点(通过蚁群算法确定的概率由轮盘算法选择 个疫灾点),并将这个疫灾点列入禁忌表。直到蚂蚁走完全部的疫灾点,终结蚂蚁的循环活动。()记录本次迭代各疫灾点总的最短等待时间,当前迭代次数中的最短等待时间即是所求的各车辆运行最佳路线上对应各疫灾点最短等待时间的和,按信息素更新式()进行信息素全局更新。()记录各疫灾点总的最短等待时间,并多次执行()(),使之达到最大迭代次数。()结果输出。模型应用以下通过实际算例对疫情应急医疗物资配送模型及优化算法的有效性进行验证。.应用数据假定某地有 个疫灾点,个应急医疗物资储备中心 点,待配送应急医疗物资的种类总数,车辆的类型总数,用 表示,且共有车 。等待装配车辆基本信息如表 所示。表 车辆基本信息.车辆类型载重量 载容量.车辆数 待配送应急医疗物资的单位重量和单位