摘要
突发事件下的应急车辆调度,既要优先考虑需求紧迫度高的关键节点,又要避免其他节点配送延时过大产生次生影响。基于此,该文建立了考虑需求紧迫度与配送延时公平的应急车辆调度模型。首先,采用熵权-TOPSIS 法计算节点需求紧迫度,并引入配送延时均衡指标避免节点之间的配送延时差异过大;其次,以配送的总时间、总成本和配送延时均衡度为多目标函数,配送中心库存量、路径连续性等为约束条件,建立优化模型;最后,通过 K-Means 算法和 CW 节约算法划分子路网,结合可变邻域局部搜索算法改进蚁群算法。结果表明:相比于基本模型,所提模型下,需求紧迫度高的节点的平均配送延时和最大配送延时分别降低 25.26%、22.98%,配送延时方差增加 34.30%;相比于仅考虑需求紧迫度的模型,所提模型下,需求紧迫度高节点的平均配送延时、最大配送延时和配送延时方差分别降低 10.79%、20.26% 和 40.57%。这表明所提模型能更好地兼顾关键节点的优先和普通节点的公平。
Abstract
In emergency vehicle dispatching after unexpected events, it is crucial to prioritize critical nodes with high demand urgency while minimizing excessive delivery delays at other nodes to prevent secondary impacts. On this basis, this study developed an emergency vehicle dispatching model that considered demand urgency and delivery delay fairness. Firstly, the entropy weight-TOPSIS method was used to calculate the demand urgency for each node, and a delivery delay balance index was introduced to minimize disparities in delivery delays between nodes. Secondly, a multi-objective optimization model was established with total delivery time, overall cost, and delivery delay balance as objectives, while considering constraints such as inventory levels at distribution centers and route continuity. Finally, the K-Means algorithm and the Clarke-Wright (CW) saving algorithm were employed to partition the sub-networks, and an improved ant colony algorithm incorporating a variable neighborhood local search was applied. The results indicate that compared with the basic model, the proposed model reduces the average and maximum delivery delays at nodes with high demand urgency by 25.26% and 22.98%, respectively, while increasing the delivery delay variance by 34.30%. Compared with a model considering only demand urgency, the proposed model decreases the average delivery delay, maximum delivery delay, and delivery delay variance at nodes with high demand urgency by 10.79%,20.26%, and 40.57%, respectively. These findings demonstrate that the proposed model effectively balances priority for critical nodes with fairness for standard nodes.
0 引言
应急事件突发时,受应急车辆的数量和配送能力的约束,无法及时满足所有受灾区域的需求,亟须合理优化车辆调度方案,避免损失加重[1]。因此,配送能力约束下的应急车辆调度优化问题成为重要研究课题。应急车辆调度方案的优劣取决于全体受灾节点的配送时效性[2],既要保证需求紧迫度高的关键节点能优先获得应急物资,又要减少其他节点的配送延时以避免受灾程度恶化。因此,如何兼顾优先和公平是应急车辆调度的难题。
国内外学者针对节点的应急物资需求紧迫度的车辆调度问题开展了广泛研究[3-4],Silva 等[5]评估节点优先级,应用排队论得到应急车辆救援路径;Yang 等[6] 依据救援政策确定优先级,以救援需求满足率最大、成本最低为目标,构建应急车辆调度优化模型; 赵建有等[7] 提出医疗物资需求紧迫度评价方法,以总配送费用最少为目标,建立考虑需求紧迫度的配送路径优化模型;Zheng 等[8] 提出考虑需求紧迫度的路径优化模型能保证受灾严重的节点优先救援,但会导致配送路径末端节点的配送延时过长。
受灾节点的损失增长呈线性,配送延时过长将导致受灾严重,同时会加剧受灾群众的不满情绪,甚至引发群体不理智行为。一些学者将配送延时等公平性指标引入目标函数中,对多目标下应急路径优化问题展开研究。Balcik 等[9] 以公平分配资源和降低运输成本为目标,优化物资配送时间表,解决末端配送难题;Cook 等[10] 以选择需求满足率作为配送公平性指标,并将其纳入应急车辆调度模型的目标函数; Sun 等[11]提出以节点伤亡加权总和最小为公平性指标,结合系统总成本与未满足救援需求的惩罚成本,建立应急车辆调度模型;王旭坪等[12]考虑受灾群众的心理因素,建立了公众心理风险感知与物资未满足度最小的双目标规划模型;郝正博等[13]考虑基于车道数量、路段长度等静态属性与拥堵延迟指数、平均车速等动态属性对应急车辆救援路径进行优化; Ghasemi 等[14]提出一种人道主义物流救援路径规划方法,以出行成本、需求满足率和路径疏散率共同作为多目标函数。
综上,现有研究主要侧重于应急车辆调度优化中关键节点的优先配送问题。虽然部分学者将节点的物资需求满足率引入优化目标中,构建多目标优化模型,但未考虑节点的配送延时均衡,导致配送路径末端节点配送延时过大的问题。本文将节点需求紧迫度作为影响因子纳入模型中,引入配送延时均衡指标,与总配送时间、总成本一同作为优化目标,构建多目标优化模型,并设计一种改进的蚁群算法提高求解速度和质量。最后,应用实际案例验证模型的有效性。
1 研究问题
1.1 问题描述
在突发事件发生时,需求节点的应急物资储备往往有限,难以满足节点需求。同时,应急车辆数量和配送能力通常受到限制,使得配送中心无法直接对所有需求节点进行同时配送。为解决这一问题,本文将由 H 个配送中心和 D 个需求节点构成的二级配送网络合理划分为多个子路网,每个子路网包含 1 个配送节点和若干需求节点。通过建立考虑配送紧迫度与公平的应急车辆调度模型,以期合理分配 K 辆车的配送次序和物资装载量,实现总配送时间最小、优先满足需求紧迫度高节点的物资需求,并确保所有节点的配送延时尽可能均衡。建立模型前,先做以下假设:
(1)配送车辆为同一型号,装载容量为 A,速度为 V。
(2)每个需求节点的配送次数不超过 1 次。
(3)应急物资配送中心不存在物资转运。
(4)各种应急物资统一计算为重量,为单种物资配送。
1.2 需求紧迫度
应急事件发生时,需求节点的应急物资往往供不应求,亟须评估各节点对物资的需求程度并依此优先配送[15]。本文依据不同应急事件特点,定义目标层、准则层和指标层[16],采用熵权-TOPSIS 法计算节点需求紧迫度 Di。
首先,熵权法计算指标权重:将指标数据采用 0~1 标准化处理,得到标准化矩阵并计算指标比重 pij 和信息熵冗余 dj,得到指标权重 wj,如式(1)所示:
(1)
(2)
其次,通过 TOPSIS 法评价节点需求紧迫度:将标准化矩阵与权重 wj 相乘,得到紧迫度矩阵 B,通过正、负理想解 B+ 与 B-,点与最优、劣解的距离 、 计算需求点紧迫度 Di。以紧迫度最小的节点 D min 为基准,计算相对需求紧迫度 λi,并作为影响因子纳入优化模型中。具体公式如式(3)~(5)所示:
(3)
(4)
(5)
1.3 配送延时均衡
配送延时的增大会导致受灾程度加重,也会使受灾者感到不公平,进而引发焦虑情绪[17]。引入配送均衡延伸 ρ 作为多目标函数之一,避免因节点配送延时过大和差异引发不良影响。
配送延时 Tl 为实际救援时间 和理论最快配送时间 的差值。Tl 越小,受灾程度越小,受灾群体心理越稳定。为确保需求紧迫度较高的节点配送延时最小化,将需求紧迫度作为加权因子 λi 纳入配送延时指标 Tl中,得到加权配送延时 ,如式(6)所示:
(6)
应急车辆救援时,既要确保高紧迫度节点配送延时的最小化,又要尽可能减少低紧迫度节点的配送延时,保证节点之间的公平性。因此,提出配送延时均衡指标 ρ,ρ 越小,节点间的配送延时差异越小,从而配送更公平,ρ 值计算如式(7)所示:
(7)
2 考虑紧迫度与配送公平的应急车辆调度模型
2.1 参数定义
H为配送中心集合;D为需求节点集合;K 为运输车辆集合; 为配送中心 h 的应急物资模糊供应量, , 对应的去模糊量为 Qh;Rh 为配送中心 h 的配送路径集合;dh 为配送中心 h 负责的需求节点集合;d h r 为配送中心 h 的配送路径 r 上的需求节点集合;λi为节点的需求紧迫度;F0为物资运输成本;Fv为车辆使用成本;A为车辆容量限制;V 为车辆行驶速度;sij、tij为需求节点 i、j 间的距离和行驶时间;为配送中心 h 的车辆到达需求节点 i 的时间; 为需求节点 i 的理论最快配送时间;Tu、Td为车辆 k 装载、卸载物资所需要的时间;Su、Sd为车辆 k 装载、卸载物资需要的费用;T k uh、T kdi为车辆 k 在物资配送节点 i 的物资装载、卸载时间;Sk uh、Sk di 为车辆 k 在物资配送节点 i 的物资装载、卸载费用;y h k 为配送中心h 的车辆 k 启用,取值为 1,否则为 0;q h ijk为配送中心h 的车辆 k 从需求点 i 到 j 配送的实际物资数量;xh ijk为配送中心 h 的车辆 k 从需求点 i 到 j 进行服务为 1,否则为 0;Tl'为考虑需求影响因子的节点配送延时;ρ 为节点间配送延时均衡度;ϑ 为置信水平;pi'为节点的实际模糊需求量。
2.2 模型建立
本研究属于多人多目标旅行商问题,将节点需求紧迫度作为影响因子纳入应急车辆调度模型中,以总配送时间、总配送成本和配送延时均衡最小为目标函数,考虑配送中心库存数量、车辆运载能力和路径连续性等约束条件下,建立应急车辆调度模型。
主目标函数为:
(8)
子目标函数为:
(9)
(10)
(11)
优化模型包括 3 个子目标,选择线性加权法将多目标问题转换为单目标问题。式(9)为总配送时间最短,包括运输、装载和卸载物资;式(10)为总配送成本最小,包括总配送费用、车辆固定成本、物资装载和卸载成本;式(11)为配送延时均衡度最小。
模型主要约束为:
(12)
(13)
(14)
(15)
(16)
(17)
(18)
(19)
(20)
(21)
(22)
(23)
(24)
(25)
(26)
(27)
(28)
式(12)~(15)为物资装卸时间和费用;式(16) 为节点 h 的物资实际到达时间;式(17)为确保配送中心间不会出现物资转运;式(18)为确保运出物资少于库存;式(19)为分配给需求点的物资少于需求; 式(20)为不超过车辆的装载能力;式(21)和(22)为保证配送路径的连续性和节点不重复出现在路径中;式(23)与(24)规定配送路径从配送中心开始; 式(25)为车辆返回至出发地配送中心;式(26)~(28) 为变量取值范围。
突发事件发生时节点的应急物资供应和需求具有强不确定性,区间数学和三角模糊数学能够灵活、准确地表达不确定性[18],本文引入其来描述不确定性,对模糊参数进行确定性转换。此方法能够灵活且直接表达不确定性,减少对数据的依赖,计算效率高。在物资需求水平 γ 与物资供应置信水平 ϑ 下,约束条件(18)与(19)转换得到式(29)和(30):
(29)
(30)
式中:为节点 i 的模糊需求量, 通常表示上下限值;pi为物资需求量。
3 求解算法
应急车辆调度模型属于非确定性多项式问题,求解十分复杂[19-20]。蚁群算法具有收敛快、求解简单等优点,但存在求解质量一般、容易陷入局部最优等缺陷[21-22],本文应用 CW 算法划分初始节点方案,引入可变邻域局部搜索算法以增加可行解空间对蚁群算法改进。算法流程如图1所示。
3.1 配送中心-节点分组
图1 算法流程
Figure1 Algorithm flow
(1)综合考虑容量限制与距离的 K-Means 聚类
随机选取需求节点为初始聚类中心并确定聚类簇数量,将节点分配到距离最近且容量满足需求的簇中,并持续更新聚类中心位置。
(2) CW 算法构造初始解
将需求节点与对应的配送中心相连形成环路。将任意两条路径融合为一条路径时,计算距离减少量。以路径的配送总量是否超过载重限制为约束,距离减少量最大为目标,对路径进行迭代更新,直至所有需求节点分配完毕。
3.2 局部搜索
可变邻域搜索算法通过交换算子、交叉算子、动态插入算子等策略,增加可行解空间,避免陷入局部最优。其中,交换算子是指对路径内的访问顺序进行改变,若路径的适应度函数增大,则保存改进后的最优路径,否则保持原有路径;交叉算子是指在可行解 g 中选择某个区段,在该区段中将节点 j 随机插入区段中除首、末位外的任意位置,得到新可行解 g';动态插入算子是指随机从可行解中选择 2 个区段 g1、 g2,分别在各自区段选择两个需求节点交换。
3.3 算法效果
选用 Solomon VRPTW 算例中 C 类数据集的前 25、50 和 100 个节点共 27 个数据集进行试验,以验证算法的通用性和普适性。本文的改进蚁群算法与蚁群算法、遗传算法的对比结果表明,在 25 节点的情况下,改进蚁群算法求解时间为 46.67 s,偏差值为 0;蚁群算法和遗传算法分别需要 24.14 s 和 95.82 s,且偏差值为 6.47 和 0.25。在 50 节点和 100 节点情况下,也表现出相同的优势。
4 案例分析
4.1 案例概况
以某市的配送医疗物资的应急车辆调度情况为案例,包括 3 个配送中心和 33 个需求节点。配送中心的基本信息包括经纬度坐标、车辆数、物资供应量;需求节点的基本信息包括所属区域人口统计、住院患者、医护人员、危重病人及物资需求量。考虑疫情感染、医疗基础设施、应急医疗物资需求 3 个因素,评估各节点需求紧迫度,结果如图2所示,可看出节点 1、27、19、25、2、17、12、18 的需求紧迫度最高,需要优先配送。
图2 节点需求紧迫度
Figure2 Demand urgency for nodes
4.2 试验结果
在 R7-5800H CPU 主频 3.2 GHz、内存 16 GB 环境下,应用 Matlab 求解。相关参数设置:蚁群规模 m=50,最大迭代次数 F max=100,物资供应置信水平 ϑ=0.95。
为检验节点需求紧迫度、配送延时均衡对配送结果的影响,设计 3 组方案:方案 1:不考虑需求紧迫度和配送延时均衡,λi = 0,w1 = 0.5,w2 = 0.5,w3 = 0;方案 2:考虑需求紧迫度,不考虑配送延时均衡, λi = 1,w1 = 0.5,w2 = 0.5,w3 = 0;方案 3:考虑需求紧迫度和配送延时均衡,λi= 1,w1 = 0.3,w2 = 0.3, w3 = 0.4。
求解结果如表1和图3所示,方案 2 和 3 优先配送需求紧迫度高的节点,与方案 1 相比,配送总时间、总成本和总配送延时均有小幅增加;方案 3 考虑了配送延时均衡,与方案 1 相比,平均配送延时减小,节点间公平性更高。
表1 相关指标
Table1 Related indicators
图3 配送路径
Figure3 Delivery route
为进一步检验模型的效果,分别针对需求紧迫度高和低的节点的配送延时进行分析,结果如表2和图4所示。方案 2 优先配送紧迫度高的节点,与方案 1 相比,紧迫度排名前 8 位节点的平均配送延时 、最大配送延时 T M ρ 和配送延时方差 σρ 显著下降,分别减少 48.55%、49.36% 和 36.82%;但紧迫度排名后 25 位节点的指标分别增加 43.50%、20.83% 和 42.59%。结果表明:仅考虑关键节点的优先配送方案将导致普通节点的配送延时大幅增加。
表2 各方案配送延时指标对比
Table2 Comparison of delivery delay indicators of various solutions
图4 节点配送延时分布
Figure4 Distribution delay distribution of nodes
方案 3 同时考虑了关键节点优先配送和普通节点配送延时均衡,与方案 1相比,紧迫度排名前 8位节点的 、T M ρ 分别减少 25.26%、22.98%,σρ 增加 34.30%;紧迫度排名后 25位节点配送效率明显优于方案2,、T M ρ、σρ分别减少 10.79%、20.26% 和 40.57%。这表明方案 3 能更好地兼顾关键节点的优先和普通节点的公平。
5 结语
在面对突发事件下的应急车辆调度问题时,既要优先考虑需求紧迫度高的关键节点,又要兼顾其他普通节点利益来确保公平。基于此,本文提出考虑需求紧迫度与配送延时公平的应急车辆调度模型,并设计改进蚁群算法,应用 CW 算法划分初始节点方案,引入可变邻域局部搜索算法,增加可行解空间,提高求解质量与速度。试验结果表明:该模型能提高需求紧迫度高节点的配送效率,并大幅降低整体配送延时,改善配送时间的稳定性。但未考虑车型和物资的差异,后期可以继续开展相关研究。