1. 引言
在全球电子商务的发展和数字支付技术的普及下,线上外卖服务需求迎来显著增长,据统计,全球外卖配送服务的市场规模从2019年的1074.4亿美元增长至2023年的1543.4亿美元,年增长率为11.51% [1]。随着以食品为主的外卖配送市场趋于饱和,近年来外卖平台开始提供药品、日用百货等配送服务以提高竞争力[2]。现阶段外卖配送主要采用零工经济模式服务客户,外卖平台将订单快速匹配给附近的配送员配送,将商品从商家送至客户指定地点。在零工经济模式下,配送员需要配送的订单数量和行使路线具有较强的自主性,并对商品的完好性全程负责。这种模式存在较大的安全风险,配送员为了避免超时或提高收入,会采取一些不安全行为,如超速、逆向行驶等,根据上海交警部门的统计,2024年全市快递外卖行业日均交通违法记录达到2462起。当遇到订单激增的情况,会出现运力不足,等待时间大幅延长。在COVID-19大流行期间,政府强制关闭或限制线下餐饮服务的同时,外卖平台推出“非接触式配送”功能使消费者接触病毒的风险降低,居民对线上外卖配送服务需求激增,依赖零工经济的配送网络无法按时满足客户的需求[3] [4]。
网络设计是提升物流配送效率的重要环节[5] [6],目前许多学者从不同的角度对如何提高城市物流配送效率进行了研究,李锐[7]等研究带硬时间窗的4PL多到多网络设计问题,设计迭代局部搜索算法求解,最小化物流总成本。Dayarian [8]等提出了一种基于循环的启发式方法用于设计城市内和部分城市间包裹配送的发车频率和商品路径,通过确定成本最低的周期集以确定最优时刻表,保证了较高的准点率。Crainic [9]等建立了考虑设施资源限制的二级城市物流服务网络设计模型,可以用于确定用户选择合适的车队,路线与时间表。杨从平[10]等改变以往学者基于轴辐式网络理论的快递网络优化方法,基于复杂网络理论对快递网络进行优化。Zhao [11]等采用两层轮辐结构来构建城市内地铁物流网络,建立熵加权TOPSIS模型选择枢纽节点和辐条节点,缓解城市交通拥堵和环境污染问题的可持续解决方案。
在上述城市物流网络设计研究中,研究问题大多聚焦在同城配送领域,其配送范围覆盖整座城市,配送时效性较差,对配送资源的统筹规划策略不适用于时效性要求较强的外卖配送服务网络。Nguyen [12]等为了提高线上外卖订购平台满足客户需求的能力,提出一种新的算法用于计算配送员在满足时间与空间约束时每个配送员的可行订单分配,并提出一种新的网络可靠性指标用于评估外卖配送网络的性能。徐蓉蓉[13]研究了派单模式下专送运力的餐饮外卖即时配送运力调度问题,设计了鲸鱼群算法求解订单分配和配送路径的联合优化。He [14]等提出餐饮外卖市场的众多独立餐厅正在通过外卖订购平台竞争顾客,而餐厅的食品质量和地理位置是两个主要考虑因素,并探究了不同配送策略对餐厅食品质量和选址决策的影响。Correa [15]等通过收集供应商与客户之间全部可能的交通路线,评估交通状况对在线送餐服务关键绩效指标的影响,发现交通状况对交易量和交货时间的履行没有实际影响。
上述针对外卖配送网络的研究仍然建立在以零工经济为基础的网络框架内,研究内容包括控制成本,提高客户满意度,设施选址等。目前较少有文献研究集散型网络在外卖配送中的应用,以提高外卖订单的配送效率。本文通过研究集散性配送网络在城市外卖配送系统中的使用场景,在不同订单规模背景下与传统零工经济模式进行比较,分析不同配送模式下的成本优势,达到整合配送资源控制配送成本,提高配送效率的目的。
2. 问题描述
集散型网络的配送区域可以区分两种配送活动,取货活动和送货活动,外卖配送网络需要满足的特殊性要求是配送员必须先在商家取到商品,才能对顾客进行送货服务,且配送员要为多个商家进行联合配送。集散型配送网络的商品配送时间等于取货时间和送货时间之和,且不超过规定的配送时间窗。
针对所研究的配送区域商家与客户的地理分布特征,可以将地理位置比较接近的商家划分成一个商家群体,同理也可以划分出若干客户群体。本文设
表示研究区域的配送网络,其中
是商家群集合,
是客户群集合,配送站点是取货员与送货员交接商品的地点,
是连接商家群与配送站点之间的链路集合,
是连接客户群与配送站点之间的链路集合。如图1所示说明了集散型外卖配送服务网络的设计结构。
Figure 1. Design structure of distributed takeout delivery service network
图1. 集散型外卖配送服务网络设计结构
3. 数学模型构建
在城市外卖配送活动中,配送平台会根据人口密度,城市布局,交通便利程度等因素建设配送站点,配送区域内的服务对象(包括商家聚集区,市民住宅区,商业办公区)的地理位置分布相对固定,配送员数量可以根据地区常驻人口数量合理分配。通过对区域内商家节点与客户节点的分布状况提前聚类,将配送员限制在某一商家群内或客户群内活动,可以避免配送员选择较远的订单配送,最大限度缩短配送员行使路程,提高整个网络配送效率。
3.1. 节点群体划分模型
3.1.1. 样本空间定义
由于集散型配送网络的配送员不能同时服务商家和客户,所以网络节点群体划分需要将商家节点和客户节点分离出来,分别做聚类分析。设商家样本集合用X表示,假设有n个节点样本,每个节点有m个属性,则商家样本集合可以用矩阵X表示为:
(1)
集散型外卖配送网络中,商家节点包含的属性有商家坐标值、商家取货时间窗、商家出货时间、商家订单数量、商家订单对应的客户坐标等。聚类的核心是根据节点的相似度划分类,节点的相似度是根据节点具有的属性定义的,选择不同的节点属性作为相似度的定义得到的聚类结果是不同的[16]。本文假设所有商家的外卖出货时间,订单数量一致,因此本文中配送员的行使距离越短,配送时间就越短,配送效率就越高,选择商家坐标值作为聚类属性。同理客户节点也用客户坐标值进行聚类分析。
3.1.2. K均值聚类数学模型
K均值聚类的策略是通过误差平方和的最小化对样本空间选取最优的划分[17],式(2)采用欧式距离的平方作为节点之间的距离
,式(3)定义节点与所属簇中心之间的距离总和为误差平方和WCSS。本节设K均值聚类的目标是将样本分为k个类,第μ个类的簇中心用
表示。
是指示函数,若节点i属于类μ取值为1,否则取值0。
(2)
(3)
K均值聚类的数学模型即求解最优化问题:
(4)
3.2. 配送路径优化模型
3.2.1. 相关假设
为了保证问题的可行性和有效性,对集散型外卖配送服务网络的路径规划问题作如下假设和说明:
(1) 外卖配送的服务区域已知,配送速度恒定,最大载重量已知不变。
(2) 假设每件外卖的重量和体积一致,因此用外卖取送的数量表示车辆目前的载重量。
(3) 配送站点,客户,商家节点的位置坐标已知,假设所有节点之间的行使时间与距离成正比,本文取各节点间的距离为欧式距离。
(4) 在取货环节,每件外卖订单只能由一位配送员取货;同理在送货环节每件外卖只能有一位配送员送货。
(5) 客户时间窗为软时间窗,允许提前送达的订单不影响客户满意度,但是延迟送达订单会降低客户满意度,需要考虑惩罚函数。
(6) 假设所有订单均送到外卖柜,并设置有服务等待时间。
3.2.2. 符号说明
本文符号定义见表1。
Table 1. Symbol definition
表1. 符号定义表
集合 |
|
取货配送员集合 |
|
送货配送员集合 |
|
商家节点集合 |
|
客户节点集合 |
|
订单集合 |
参数 |
|
节点i与节点j之间的距离 |
|
配送员从节点i到节点j的旅行时间 |
|
商家节点i的时间窗 |
|
客户节点j的时间窗 |
|
配送中心的时间窗 |
|
节点i的服务等待时间 |
|
取货员配送订单
到达商家节点i的时刻 |
|
送货员配送订单
到达客户节点j的时刻 |
、
|
最大载重量 |
|
配送中心节点 |
|
第k个配送员的配送路径 |
|
第k个配送员拜访商家客户节点的总数量 |
|
商家i的第
个订单 |
|
客户j的第
个订单 |
决策变量 |
|
订单g由取货员k配送取1,否则取0 |
|
订单g由送货员k配送取1,否则取0 |
3.3. 数学模型
基于上述分析,集散型外卖配送服务网络设计及路径优化模型可表示如下:
(5)
(6)
(7)
(8)
(9)
(10)
(11)
(12)
(13)
(14)
(15)
(16)
(17)
(18)
(19)
(20)
(21)
式(5)是目标函数,表示最小化全局配送员行使总路程,包括取货行程和送货行程;式(6)表示配送员行使时间计算公式;式(7)和式(8)表示配送员k的取货路径集合和送货路径集合,其中的元素
和
分别表示第k条路径的第
个取货顺位和第
个送货顺位,
表示配送中心且
;式(9)和式(10)是时间流平衡约束,其中式(9)表示取货环节中,第k个取货员到达商家节点i − 1的时刻加上在i − 1节点的等候服务时间再加上前往商家节点i的旅行时间等于到达商家节点i的时刻;式(10)同理为送货环节送货员到达客户节点j的时刻;式(11)和式(12)表示配送员k从配送中心出发的时间和返回配送中心的时间必须在配送中心的左右时间窗之间;式(13)表示若取货员早于出货时间
到达商家节点i,则总等待时间为等待出货时间
加上0.5分钟的固定服务时间(如排队,整理包裹等);式(14)表示送货服务时间为0.5分钟,这里假设外卖允许提前送达并放置在外卖柜;式(15)表示送货员不得晚于客户最晚要求的送达时间;式(16)和式(17)表示每件外卖订单必须得到服务并且要求取货服务一次,送货服务一次;式(18)和式(19)表示配送员的装载容量约束;式(20)和式(21)为决策变量。
4. 基于K-Means聚类与遗传算法的模型求解
集散型外卖配送网络问题涉及到商家群体与客户群体的划分以及订单配送过程中配送员的订单分配与路径规划问题。因此,该问题属于时间与容量约束下的路径–分配组合优化问题[18]-[20]。本文使用K-means聚类算法和带有局部搜索策略的遗传算法分两阶段求解集散型外卖配送网络。
4.1. 聚类算法求解步骤
K-means聚类需要分两个阶段的算法完成:第一阶段利用手肘法确定聚类类别数k。在手肘法中,聚类结果的质量可以用类中每个点到簇中心的距离平方和(WCSS)来衡量。在算法第二阶段选择k个数据点作为初始簇中心,然后将离初始簇中心较近的数据节点划分到此簇中心所属的类别,再针对每个已经划分好的类别重新计算簇中心,使节点到簇中心的误差平方和最小。本文采用如下流程,使用Python求解器求得最优k值和聚类结果:
Step 1 数据初始化 导入sklearn聚类数据库,绘图数据库;将商家与客户节点经纬度坐标转化为平面直角坐标数据导入。
Step 2 计算WCSS数据 根据式(3)计算k取值1到11时的类中每个点到簇中心的距离平方和(WCSS)并记录。
Step 3 绘制肘部图 绘制k取不同值时WCSS数值的变化趋势图,取变化率最大的拐点作为最优k值。
Step 4 初始化数据中心点 随机选择k个数据点作为初始簇中心。
Step 5 划分初始群体 将每个数据节点分配到与其最近的中心点所属的类别。
Step 6 迭代优化 对于每个类别,根据式(4)重新计算中心点位置,使WCSS最小。
Step 7 输出聚类结果 分别得到配送区域内商家节点划分结果与客户节点划分结果。
4.2. 遗传算法求解步骤
4.2.1. 染色体编码
本文使用了一种带有局部搜索策略的遗传算法求解配送员的订单分配与路径规划问题。在遗传算法中每条染色体代表一个配送员,其标记了订单的配送顺序。如图2所示,与传统零工经济模式下的外卖配送染色体编码方式不同,集散型配送网络的染色体编码只需要区分取货人员和送货人员,取货人员只经过商家节点,送货人员只经过客户节点。
Figure 2. Coding diagram
图2. 编码示意图
4.2.2. 构建适应度函数
集散型外卖配送服务网络路径优化模型的适应度函数(见公式23)包括两个部分:配送员行使总距离F1和违反时间窗约束的惩罚函数F2。其中式(5)表示了配送员的行使总距离,式(22)表示配送员只有晚于右时间窗到达客户节点才算违反约束。式(23)中,a1和a2分别为F1和F2的权重系数,可以根据各自在总的目标函数中所占的重要性来设定。本文根据所研究的外卖配送网络的实际情况,并采用实际数据模拟之后将权重设置为
,
较为合理。
(22)
(23)
4.3. 局部搜索策略
在局部搜索操作中,可以破坏当前迭代子种群中的部分染色体并进行修复,以进一步提高种群质量[21]。算法流程如下:
Step 1 输入当前子代个体,计算子代染色体适应度值CF。
Step 2 移除当前染色体中某一节点元素。
Step 3 将此元素重新插入当前染色体其它位置。
Step 4 计算重新插入元素后当前染色体的适应度值RCF。
Step 5 如果RCF > CF,则重新将此元素插入其它位置,并设置最大重新插入次数为
或
,即不超过配送员的最大载重数量。
Step 6 如果RCF < CF,则替换原子代个体形成新种群。
5. 仿真案例分析
5.1. 案例数据生成
本文选取某平台2019年12月7日的实际运营数据作为预设场景进行案例分析,订单信息包括订单产生时间,订单出货时间,预期送达时间,商家坐标和客户坐标。
因为外卖订单具有连续产生的特点,集散型外卖配送网络需要取货环节和送货环节配合进行。因此本案例仿真通过规范配送员的取货时间窗和送货时间窗,以保证每个订单在规定时间窗内取货和送货。配送时间窗划分如图3所示,以11:00~11:10分产生的订单为例,外卖取货人员需要在时间段11:00~11:25内完成取货任务并返回配送中心,外卖送货人员需要在时间段11:25~11:50内完成订单的配送任务并返回配送中心。案例通过模拟10分钟产生的订单分配与路径优化结果推算90分钟的高峰时段所需配送员数量和成本,减少了模型计算复杂度和运算时间。
Figure 3. Order delivery time window
图3. 订单配送时间窗
仿真案例的实验参数设计如下:根据已有文献的调查,外卖快递配送人员的每小时使用成本为5元;零工经济模式下使用二轮电动摩托车配送每公里运营成本是0.2元;外卖配送站点的覆盖半径在5公里以内,城市内配送速度为30公里每小时;商家和客户节点的服务时间均设置为0.5分钟。为保证实验结果的有效性,零工经济模式和集散型配送网络均使用相同的参数。
仿真案例中使用的不同订单规模信息见表2。其中订单规模指的是每10分钟产生订单数量。
Table 2. Scale of the case
表2. 案例数据规模
订单规模 |
商家节点数 |
客户节点数 |
20 |
15 |
20 |
30 |
20 |
30 |
40 |
23 |
40 |
55 |
27 |
55 |
70 |
31 |
70 |
80 |
35 |
80 |
90 |
36 |
90 |
100 |
38 |
100 |
110 |
40 |
110 |
120 |
41 |
120 |
130 |
43 |
130 |
5.2. 案例仿真实验结果
5.2.1. 案例节点聚类结果
首先对43个商家节点的聚类结果进行分析,商家群体1和商家群体3节点分布较为密集,商家间距较短,仿真过程中使用较少的配送员便可完成配送任务,有利于体现集散型网络的规模效应。商家群体3中商家间距较长,配送员配送耗时较长,但是由于节点数量较少,订单规模不大,仿真过程中使用较少的配送员也可以完成配送。如图4所示为仿真案例中商家节点的聚类结果。
Figure 4. Clustering results of merchant nodes
图4. 商家节点聚类结果
5.2.2. 案例订单分配与路径优化结果
本文列举订单规模为20单时分别采用零工经济配送网络和集散型配送网络求得的算法结果,见表3所示。其中客户节点用数字1~20表示,商家节点用数字21~40表示。数组(21, 1)表示一个订单从商家21送到客户1。从仿真结果可以看出,在满足配送时间窗的约束下,订单均能在40分钟以内送达,同时集散型网络需要更多的配送资源,但是每个配送员的行使里程大幅缩短。
Table 3. Comparison of case simulation data under two distribution modes
表3. 两种配送模式下的案例仿真数据对比
配送模式 |
配送员 |
配送路线 |
配送耗时min |
行使距离km |
集散型配送网络 |
取货1 |
0-24-29-37-30-28-36-32-21-0 |
8.17 |
1.39 |
取货2 |
0-39-25-31-22-33-40-0 |
14.01 |
3.67 |
取货3 |
0-35-34-27-26-23-38-0 |
21.45 |
6.15 |
送货1 |
0-1-17-16-12-2-0 |
13.68 |
7.45 |
送货2 |
0-5-15-4-8-18-19-0 |
12.60 |
6.40 |
送货3 |
0-20-10-13-9-0 |
16.81 |
9.87 |
送货4 |
0-3-14-6-7-11-0 |
20.29 |
11.86 |
零工经济配送网络 |
配送1 |
0-31-24-38-18-27-4-33-7-11-13-0 |
39.92 |
23.28 |
配送2 |
0-35-29-33-9-40-34-15-14-20-12-0 |
37.94 |
21.96 |
配送3 |
0-39-22-26-28-6-23-3-2-19-8-0 |
38.83 |
22.55 |
配送4 |
0-25-21-36-30-37-10-17-16-1-5 |
25.49 |
13.66 |
Figure 5. Trend of the number of deliverymen under different order sizes
图5. 不同订单规模下配送员人数变化趋势
图5是不同订单规模下两种外卖配送模式所需配送员人数趋势变化。集散型外卖配送网络在订单规模小于90件/10min时,所需的配送员人数略高于传统零工经济模式,这是由于集散型网络在划分客户群体之后,需要保证每个群体内都有配送员可以提供服务,并且为了保证较高的配送效率,集散型网络不允许配送员跨群体服务商家客户,即使某些配送员仍然有多余的运力可以提供。在订单规模大于90件/10min时,集散型网络所需配送资源明显小于传统零工经济网络,在高峰时期配送资源不足的情景下集散型网络具有明显优势。
5.3. 每公里运营成本对网络设计的影响
外卖配送成本包括两部分配送人员使用成本和人员运营成本。其中人员使用成本 = 单位人员使用成本 * 所需配送员数量;人员运营成本 = 每公里运营成本 * 配送员行使总里程。图6是两种配送模式下配送成本的对比。
当两种配送模式的单位运营成本均为0.2元/公里时,随着市场需求的不断增加,集散性配送网络与零工经济模式下的配送成本差距不断增大,因此集散型网络具有明显的成本优势。原因在于集散型网络限制了配送员的行动路线和活动区域,比零工经济模式更有效的缩短了行驶里程,减少了运营成本,同时订单规模超过一定的阈值,可以使用较少的配送员,人员使用成本也有所下降。
为了更深入分析单位运营成本的变化对配送总成本的影响,本节将集散型配送网络的运营成本提高至0.4元/公里,0.6元/公里和0.8元/公里,并同现有零工经济模型下的配送成本对比研究发现,在运营成本提高至原来的3倍时,集散型网络就不再具有明显的成本优势。因此,外卖配送平台采用集散型配送网络,可以通过提高配送员载货量,采用更先进,速度更快的外卖配送车辆以提高配送效率,由此造成的运营成本增加控制在0.6元/公里以内是较为经济的。
Figure 6. Comparison of distribution cost change trends of the two distribution modes
图6. 两种配送模式的配送成本变化趋势对比
6. 总结
集散型网络并不是一个全新的概念,在物流配送领域尤其是同城配送领域具有广泛的研究和使用,然而在小范围内的外卖配送领域较少有研究和使用。此类网络设计的核心是对提前对配送区域内的商家和客户进行群体划分,让配送员尽可能在所划分好的群体内部进行配送活动,并通过路径优化尽可能减少配送员行走重复路线,达到缩短配送行使里程,控制成本的目的。在客观上限制配送员的服务区域,便于管理配送员的行驶路线,从网络设计层面提高了配送活动的安全性。
本文通过案例仿真验证了集散型网络在外卖配送服务网络中使用的可行性,但是本文案例中使用的客户节点分布较为均匀,不适合聚类分析,因此仅对商家群体进行了划分,为进一步评估集散型网络的性能特点,未来可以结合实际情况,分析仅对客户群体进行聚类划分,同时对商家群体和客户群体进行聚类划分等不同情景下集散型网络的适用条件。
基金项目
国家社会基金一般项目:新冠疫情中“信息疫情”的形成机理及应对策略研究(项目编号:22BGL240)。