1. 引言
在当前农村电商蓬勃发展的背景下,农村物流的高成本问题日益凸显。相关数据显示,2024年全国规模以上农产品加工企业营业收入已约18万亿元[1],这表明农村电商的发展潜力巨大,但落后的物流体系严重制约了其进一步发展。为加快发展农村物流,有效促进“农产品进城”和“工业品下乡”,助力乡村振兴,亟须对当前存在的问题予以探究解决[2]。目前农村地区的末端物流成本居高不下,相关数据显示,末端物流成本占农村地区配送总支出的30%至40%,且为城市配送成本的2至5倍[3]。农村地区经济社会发展还比较滞后,农村物流基础设施不完善、运力资源不足、信息化水平偏低、城乡一体化程度不高等问题仍较为突出,严重制约了农村电商的进一步发展[4]。交通基础设施提质是农村物流效率提高的基础[5]。而对于通过修缮路网、拓宽道路等改善交通基础设施的方式来提升物流效率的方法来说,其投资规模大、建设周期长、维护成本高,且在地形复杂、人口稀疏的农村地区边际效益有限。近年来,兴起的无人机技术为解决农村物流配送难题提供了新思路。首先,卡车–无人机协同配送模式无需大规模改造现有路网,即可通过空中路径直接覆盖偏远地区。其次,无人机与货运车辆的联合配送模式,相较于单独使用无人机或货车既能充分发挥货运车辆大容量、长距离的优势,又能利用无人机高效解决近距离、零散的末端配送需求。通过无人机分担短途任务,有效减少了货运车辆的往返次数和行驶里程,从而在提升配送覆盖率和时效性的同时,显著节约了配送时间和运营成本[6]。
本文构建了考虑客户配送方式分型、货物重量约束与无人机续航约束的混合整数规划模型,明确区分卡车与无人机的服务范围,强化时空协同约束,并提出一种两阶段启发式求解框架:第一阶段通过K-Means聚类实现客户分群与任务预分配,压缩搜索空间;第二阶段设计混合遗传算法(HGA)嵌入变邻域搜索(VNS),在优化卡车路径的同时动态调整无人机任务序列,提升求解质量与算法鲁棒性。
2. 国内外研究综述
对于农村电子商务发展的研究来说Yuxuan Hou (2024) [7]指出数字鸿沟和技能不足是阻碍农村电商发展的主要障碍,并呼吁通过政府投资基础设施和提供技术培训来解决这些问题。Lingxuan Zhao (2024) [8]通过实证分析发现城市数字化转型能通过电子商务等中介作用促进乡村振兴,且该效应在东部地区更为显著。潘悦(2024) [9]研究了基于江苏数据,通过实证分析指出农村物流与电商已实现高质量协同,并针对发展挑战提出建议。
就农村物流配送而言,林晓伟(2024) [10]等人探讨了促使共同配送相关主体间达成全面合作的机制和建议,为共配主体间达成的合作提供理论依据。王莉(2025) [11]等人基于农村客货邮融合发展背景,研究了客运班车–电动车–无人机协同配送模式下山区农村寄递物流路径规划问题。Xiaojing Zhu (2023) [12]探讨了大数据在农村物流配送中的应用,提出了基于物联网和云计算的智能农村物流配送模式。Luping Song (2025) [13]构建了农村物流末端配送到家自提混合模式的LRP模型,并将蚁群算法与人工势场算法相结合,提出了一种改进的势场蚁群算法,实现了人工智能辅助下的物流末端配送效率优化。
在卡车–无人机协同配送的模型层面,蒋丽(2024) [14]等人在考虑无人机载重及飞行范围的条件下,建立了以卡车与无人机综合配送成本最小为目标建立混合整数规划模型。付朝晖(2025) [15]等人综合考虑山区道路通行状况、无人机能耗、容量、飞行速度、生鲜农产品新鲜度、卡车容量与速度等因素,以总集货成本最小为目标,构建卡车–无人机联合集货的路径规划模型。Xianlong Ge (2025) [16]等人针对如何合理分配有限的物流资源提高客户满意度,提出了一种利用客户价值分类和模糊时间窗口管理的卡车–无人机协同路由多目标优化模型。
在卡车–无人机协同配送研究的算法层面,Zhang Zaifang (2025) [17]等人开发了一种改进的平衡优化器模拟退火(IEOSA)算法来最大限度地降低运输成本。Wenhao Peng (2025) [18]等人提出了一种基于深度强化学习的端到端算法,能够高效求解考虑动态路况、时间窗与无人机协同的复杂车辆路径问题,并在求解质量、效率和鲁棒性上均明显优于传统方法。Lingrui Hong (2025) [19]等人针对结合卡车与无人机的旅行商问题,提出了以自适应大邻域搜索(ALNS)为主框架,并嵌入融合变邻域下降与数学规划的两阶段局部搜索(2P-LS)的混合算法,通过与Gurobi对比验证了其有效性。Shakoor Barzanjeh (2025) [20]等人开发了一种基于精确逻辑的Benders分解算法通过最早发布日期规则和动态规划算法来解决最小化总完工时间为目标的卡车–无人机配送问题。
综上,现有关于农村电商物流配送模式的研究日益受到关注,学者们已充分认识到农村物流的瓶颈问题,并在协同配送模式探索、复杂约束建模以及先进算法应用方面取得了显著进展。为提升农村物流效率提供了重要的理论依据和实践参考。但多数研究模型较为通用,未能紧密结合农村电商物流特有的场景,如客户极度分散、道路网络稀疏、货物类型差异大等。其次,在算法层面,虽有研究引入了高性能算法,但在解决农村大规模、多约束的协同路径问题时,往往面临求解效率与稳定性不足的挑战,缺乏对问题结构的有效利用以压缩搜索空间。
3. 模型构建
3.1. 问题描述与假设
本研究的问题为:针对农村电商物流网络,本研究聚焦于从县域中心仓到分散末端节点的配送优化问题。为应对单纯卡车配送效率低下的挑战,本文探索引入无人机进行协同配送。无人机能够突破道路限制、快速投递小件电商包裹,弥补卡车短板,但其作业范围受电池续航约束。在此背景下,研究卡车与无人机协同的路径规划与调度优化,旨在综合提升配送效率、降低运营成本,是解决农村电商末端配送难题的有效技术路径。
3.2. 价值驱动场景描述
3.2.1. 主要假设
1) 配送网络有一个配送中心
、多个客户点
;
2) 卡车必须经过配送中心,可以多次停靠不同配送点;
3) 无人机续航能力不受任何干扰;
4) 一架无人机一次可访问多个客户点,配送完成后无人机必须返回至发射点;
5) 客户需求分为卡车可配送和仅无人机可配送;
6) 无人机到达客户点即为配送完成,不考虑客户实际取货时间;
7) 无人机的飞行距离受电池容量限制,每次执行任务必须在航程范围内返回卡车。
3.2.2. 变量定义
1) 基础变量见表1。
2) 决策变量见表2。
Table 1. Basic variable symbols and meanings
表1. 基础变量符号及其含义
变量符号 |
变量含义 |
|
顾客需求点集合,
;包括配送中心(编号0)以及n个客户点 |
|
客户集合点,
|
|
无人机的编号集合,
|
|
无人机服务的客户集合 |
|
卡车的最大载重 |
|
无人机的最大载重 |
|
无人机的最大飞行时间 |
|
无人机的最大飞行半径 |
|
节点i到j的欧式距离 |
|
卡车从节点i行驶到j所需的时间 |
|
卡车从节点i到k所需的飞行时间 |
|
卡车单位距离运输成本 |
|
无人机单位距离飞行成本 |
|
客户点的i货物重量 |
|
客户i的最迟服务时间 |
|
时间惩罚系数 |
Table 2. Decision variables and their meanings
表2. 决策变量及其含义
变量符号 |
变量含义 |
|
若卡车从i到j行驶,则
,否则为0; |
|
若无人机k为顾客点i服务,则
,否则为0 |
|
客户i实际被服务时间 |
|
客户k由卡车服务,则
,否则为0 |
|
客户i的访问顺序 |
|
若无人机k的起始服务时间 |
|
若无人机k的任务持续时间 |
3.3. 目标函数和约束条件
(1)
(2)
(3)
(4)
(5)
(6)
(7)
(8)
(9)
(10)
(11)
(12)
(13)
式中:式(1)表示目标函数为卡车运输成本、无人机飞行成本以及时间惩罚成本的和的最小值;式(2)表示每个卡车停靠点最多只能访问一次;式(3)表示卡车路径为封闭路径,确保卡车路径连续,卡车从配送中心出发并返回配送中心,形成闭环;式(4)表示卡车载重约束,卡车运载货物总重量不得超过卡车最大载重;式(5)和式(6)表示防止卡车子回路;式(7)表示无人机载重约束;式(8)表示无人机飞行距离约束,无人机飞行距离不得超过无人机的最大飞行半径;式(9)表示无人机飞行时间限制;式(10)表示无人机任务闭合约束,无人机从卡车停靠点释放后,必须返回同一停靠点或指定集结点;式(11)表示每个货物只能由一种方式配送;式(12)表示客户服务的时间窗限制;式(13)表示卡车与无人机的协同时间同步。
4. 算法设计
未解决具有混合路径约束和协同调度关系的“卡车–无人机协同配送问题(Truck-Drone Routing Problem, TDRP)”本文采用两阶段启发式算法框架,结合聚类预分配、遗传算法和贪心调度策略,有效提升计算效率和质量。
4.1. 算法框架
本研究算法分为两个阶段。第一阶段通过K-Means将客户根据地理位置和货物属性进行划分,将目标区域划分为“卡车服务集合”和“无人机优先集合”并以卡车服务点作为聚类中心对无人机客户点集合进行聚类,降低整体搜索空间。第二阶段使用遗传算法优化卡车路径,并基于卡车路径构造贪心式无人机任务分配方案,实现两者协同工作。算法框架如图1所示。
Figure 1. Algorithm process
图1. 算法流程
4.2. 协同路径构造
4.2.1. 第一阶段K-Means聚类算法设计
农村电商物流配送的最大特点在于客户地理位置分布不均,部分客户集中于主干公路周边,部分位于交通不便、道路狭窄或海拔较高的偏远地区。同时,货物类型呈现出显著差异:工业品下行以小件、轻量为主,农产品上行常为重货或大宗。在这样的背景下,统一优化所有客户路径会导致问题规模指数膨胀,并难以在考虑续航、载重和可达性等多重约束下得到高质量解。因此,有必要在算法初始阶段对客户进行聚类,本阶段算法对所有客户点进行聚类,并确定聚类中心,即卡车的配送点,然后对同一聚类中的客户,实现无人机到客户点的最大距离小于无人机最大飞行半径,同时客户点需求不大于无人机最大载重。具体步骤为:
步骤1:初始化聚类中心。首先对所有配送客户点进行分类,根据货物的重量及无人机配送的适应性标准将客户分为卡车服务集合和无
人机服务集合,当货物重量大于
无人机最大载重时强制划分到卡车服务集合。从所有卡车可服务的客户点集合随机选择一个客户点作为初始聚类中心,即设聚类数K=1,初始聚类中心为
。
步骤2:客户点划分。对于每一个无人机客户点
,计算其与当前所有聚类中心点
的欧式距离
将无人机客户点
分配给其到聚类中心集合中距离最近的中心点
。
步骤3:更新聚类中心。对当前每一个子簇计算该簇内所有客户点的几何中心位置,并在该类中重新选择距离该形心最近的卡车客户点作为新的聚类中心。如果新中心与旧中心相同(或者所有簇的中心都不变),则进行有效性验证;否则,用新中心替换旧中心,并回到步骤2。
步骤4:有效性验证。先对每个无人机客户点i,计算其与簇中心
的欧式距离
,若有距离
大于无人机最大服务半径
的则转至步骤5,若无则转至步骤6。
步骤5:重新划分簇。将聚类中心数加一,即K = K + 1;从未被选择的聚类中心的卡车客户点中,随机选择新的卡车点加入聚类中心集合;回到步骤2。
步骤6:输出聚类结果。
4.2.2. 第二阶段遗传算法算法设计
本阶段在第一阶段聚类结果基础上,采用改进的混合遗传算法(HGA)求解卡车–无人机协同配送路径问题,算法融合遗传算法的全局搜索能力与变邻域搜索(VNS)的局部优化能力,通过双层编码机制实现资源协同优化。
1) 编码与初始化
1、双层编码方案:卡车路径层采用排列编码表示路径顺序,染色体
,其中
(卡车簇客户),首尾0表示仓库;无人机任务层采用二进制矩阵
表示无人机j服务客户i (
),每列对应一架无人机的任务序列。
2、初始种群生成:卡车路径采用节约算法(Clarke-Wright Savings)构建初始路径,通过计算合并路径的节约值(
),优先合并节约值高的节点对,在合并路径时实时计算载重,若超限则跳过该合并。无人机分配对无人机簇客户,随机选择满足
且
的分配方案。
2) 适应度函数:综合考虑运输成本与时间惩罚,定义适应度值F:
其中
,
为时间惩罚系数,
为实际服务时间。
3) 遗传操作设计
选择算子:在遗传操作设计中,选择算子采用锦标赛选择(Tournament Size = 3)保留优质个体,同时引入精英保留策略(Elitism Rate = 10%),交叉算子针对卡车层使用顺序交叉(OX)保留父代路径片段并补全剩余客户,无人机层则采用均匀交叉以概率
交换父母代的列向量,变异算子对卡车层执行2-opt局部优化随机选择两个节点反转路径,无人机层则对变异个体(概率
)随机交换两行赋值。这一系列操作通过锦标赛选择确保种群多样性,精英保留策略维持优秀个体,顺序交叉和均匀交叉分别优化卡车路径和无人机任务分配,2-opt变异提升卡车路径局部最优性,无人机层变异则增强任务分配的灵活性,共同保证算法在搜索空间中的全局探索与局部开发能力。
4) 变邻域搜索增强(VNS)
在每代遗传操作后对最优个体执行变邻域搜索(VNS)以增强局部优化能力:首先构建三种邻域结构——
对卡车路径相邻客户进行交换,
将当前无人机任务重新分配至续航更短的无人机,
合并两个相邻卡车停靠点以释放无人机服务能力;采用循环搜索策略依次在
中寻找改进解,若任一邻域操作获得更优解则立即返回
重启搜索过程,通过动态切换邻域空间避免陷入局部最优。
5) 约束处理机制
在解码阶段,针对卡车载重约束采用贪婪插入法进行客户点分配,若当前路径载重超过卡车容量上限,则优先移除边际收益(单位重量配送成本收益)最低的客户点,并将其重新分配至无人机配送队列;对于无人机续航约束,通过动态调整服务序列中的客户点访问顺序,确保单架无人机的任务总时间满足
的续航限制,若无法满足则将该无人机的超时任务拆分为多个子任务并由不同无人机协同完成;同时,为满足时空协同约束,强制要求卡车停靠点的停留时间
。
6) 终止条件
遗传算法应设置终止条件避免算法无休止运行下去,本算法中满足以下任一条件即终止:当最大迭代次数
时;当最优适应度值
连续20代未提升时,判定算法陷入局部最优或收敛,提前终止;适应度值的标准差
小于阈值
时提前终止。
5. 算例分析
5.1. 实验参数设置
为验证本文提出的模型和两阶段启发式算法的有效性和可行性,应用python软件,借助Solomon的数据集生成随机算例进行仿真测试,实验的运行环境为Intel(R) Core(TM) i7-9750H CPU @ 2.60GHz,2592 Mhz,实验参数设置见表3,此外本文设置单架无人机单位距离成本为0.5元,一架无人机固定费用为20元,无人机最大载重10 kg,最大飞行半径15 km,最大飞行时间45 min,派遣准备时间3 min,回收需1 min,无人机配取货需1 min,无人机飞行速度为80 km/h。卡车单位距离成本2元/km,一辆卡车的固定费用(折旧费)为80元,卡车行驶速度为50 km/h,卡车最大载重750 kg,卡车1次配取货时间为3 min。时间惩罚系数为0.1。
Table 3. Parameter settings
表3. 参数设置
参数 |
数值 |
种群规模 |
200 |
最大迭代次数 |
200 |
交叉概率 |
0.9 |
变异概率 |
0.1 |
Table 4. Algorithm results
表.4. 算法结果
算例信息 |
无人机与车辆协同配送 |
卡车单独配送 |
算例编号 |
客户个数 |
卡车成本 |
无人机成本 |
时间惩罚成本 |
综合成本 |
卡车成本 |
时间惩罚成本 |
综合成本 |
RC101 |
100 |
1173.99 |
252.73 |
3038.76 |
4465.49 |
1464.29 |
9150.59 |
10614.88 |
Rc102 |
100 |
1173.99 |
252.73 |
2942.18 |
4368.90 |
1464.29 |
8948.89 |
10413.18 |
RC103 |
100 |
1173.99 |
252.73 |
2883.28 |
4310.00 |
1464.29 |
8775.51 |
102239.80 |
R101 |
100 |
1065.80 |
352.42 |
2060.85 |
3479.06 |
1616.93 |
9773.59 |
11390.52 |
R102 |
100 |
1065.80 |
352.42 |
1819.63 |
3237.85 |
1616.93 |
9514.19 |
11131.12 |
R103 |
100 |
1065.80 |
352.42 |
1747.68 |
3165.89 |
1616.93 |
9302.19 |
10919.12 |
C201 |
100 |
946.74 |
290.71 |
10238.18 |
11475.63 |
1407.94 |
38919.61 |
40327.55 |
C202 |
100 |
946.74 |
290.71 |
9270.51 |
10507.96 |
1407.94 |
37241.49 |
38649.42 |
C203 |
100 |
946.74 |
290.71 |
7626.75 |
8864.20 |
1407.94 |
33886.27 |
35294.16 |
5.2. 算法求解及结果分析
Solomon数据集中R类数据集的节点坐标较为分散,C类数据集中,节点坐标有明显聚簇,RC类数据集则混合了随机和聚簇的特点。R1、C1、RC1类的节点时间窗间隔显著小于R2、C2、RC2类。在C201、C202、C203算例中,协同模式的优势最为巨大。因为客户点聚集,卡车可以停靠在少数几个中心点,无人机就能高效覆盖周边的大量客户,协同效应发挥到极致。
综合成本是衡量配送方案整体经济效益的核心指标。如表4所示在所有测试算例中,卡车–无人机协同配送模式的综合成本均远低于传统卡车单独配送模式。协同模式相比单独配送,综合成本平均降低了约60.23%。在某些算例中,优势尤为惊人,例如在C201算例中,协同配送综合成本为11475.63,而单独配送为40327.55,成本降低了71.54%;在R101算例中,成本降低了69.45%。时间惩罚成本是两种模式差异最显著的指标,直接反映了配送方案的时效性。如表所示,协同模式的时间惩罚成本极低。如,在RC101中,时间惩罚从9150.59降至3038.76降低66.8%;在C201中,从38919.61降至10238.18降低73.7%。
根据算法求解结果分析显示,卡车–无人机协同配送模式能够显著降低时间惩罚成本,从而提高客户满意度,提高物流效率,同时在卡车成本方面,协同配送模式利用了无人机来协助降低总运输成本。
6. 结论
本研究针对“最后一公里”配送成本高、效率低的难题,构建了一种卡车与无人机协同配送的路径优化模型,并通过对一系列标准算例的仿真计算,验证了所提出模型与算法的有效性。核心研究结论如下:
1) 协同配送模式在综合成本上具有绝对优势:如计算结果所示,在所有测试算例中,卡车–无人机协同配送的综合成本均显著低于传统卡车单独配送。成本降低幅度巨大,普遍节省了50%以上的综合成本。例如,在C201算例中,协同配送的综合成本为11475.63,而单独配送成本高达40327.55,成本降低了约71.5%;在R101算例中,成本也从11390.52降至3479.06,降低约69.5%。这充分证明了协同模式在经济上的巨大潜力和可行性。
2) 协同配送能极大提升时效性,减少时间惩罚成本:数据中最显著的差异体现在“时间惩罚成本”上。协同配送通过无人机进行并行、快速的末端送达,极大地缩短了客户等待时间,从而将时间惩罚成本降低了65%至85%,RC101从9150.59降至3038.76;C201从38919.61降至10238.18。这表明该模式不仅能降低成本,更能大幅提升客户服务水平与满意度,为解决时效性要求高的配送场景提供了有效方案。
但是本文模型基于静态环境客户需求、位置、路况、天气在规划期内已知且不变,并假设无人机性能不受任何干扰,这种简化假设可能导致实验结果比实际结果乐观。基于以上分析,未来的研究可以发展动态路径规划和随机规划模型,考虑随机出现的客户订单、动态变化的交通状况和天气信息;构建更真实的无人机能耗模型,将风速、风向、气温和负重作为变量,建立更精确的无人机能耗与飞行速度模型。