CN115032997B - 一种基于蚁群算法的第四方物流运输路径规划方法 - Google Patents

Abstract

本发明公开了一种基于蚁群算法的第四方物流运输路径规划方法，涉及第四方物流运输优化领域，该方法采用有向图的邻接表方式保存第四方物流运输网络的信息，并提出一种情景式编码方法，基于问题特征设计了启发式信息和信息素的计算方式，并结合节点禁忌表和路径禁忌表来确定可访问路径，利用蚁群算法可以在运输成本约束和运载量约束下得到运输路径，以实现物流作业整体的总运输时间最短的目标，求解难度较低且解的质量较高。

Description

一种基于蚁群算法的第四方物流运输路径规划方法

技术领域

本发明涉及第四方物流运输优化领域，尤其是一种基于蚁群算法的第四方物流运输路径规划方法。

背景技术

随着信息技术的发展，物流体系逐渐完善，人们享受着足不出户便利购物体验。在这种大背景下，各类商家为了进一步降低成本提高效率，将物流业务转交给第三方物流(third party logistics，简称3PL)公司。但由于物流市场的不断拓展下电商平台订单膨胀式的增加，3PL逐渐展现出在服务种类、资源协调、信息化程度等方面的短缺。第四方物流(fourth party logistics，简称4PL)为了弥补这些不足应运而生，它开创一种崭新的物流模式，通过整合3PL的物流资源，使用新的信息技术来构建一套完整的供应链解决方案。4PL在商家和第三方物流公司之间起到了至关重要的作用，对于商家可以专注与自身的核心业务，第三方物流公司则可以进一步提升物流效率。

4PL的优势是显而易见的，得到了国内外物流领域的广泛关注。4PL的业务流程可分为三个不同的阶段：客户委托、运输、交付。其中运输阶段是重中之重，4PL需要在接到委托任务后，设计运输路线，选择合适的3PL代理商，在一定的成本约束预算下和客户期望时间内找到解决方案。所以路径与物流时间优化问题对于4PL而言是极为重要的，相比于传统的3PL需要额外综合考量供应商的运输能力、运输成本等因素，这使得4PL路径问题求解空间更大，求解难度更高。所以对第四方物流优化的研究具有很强的学术价值和现实意义。

进化计算和运筹学方法是求解4PL相关优化问题的常用解决方法。卢福强和任亮等人分别在第十三届自动化科学与工程国际会议和第二十六届中国控制与决策会议上采用遗传算法和蚁群算法求解4PL优化问题。但是，前者在面对较为复杂的模型下算法效果欠佳，后者算法时间复杂度较高，从而会消耗较多的时间求解问题。

发明内容

本发明人针对上述问题及技术需求，提出了一种基于蚁群算法的第四方物流运输路径规划方法，本发明的技术方案如下：

一种基于蚁群算法的第四方物流运输路径规划方法，该方法包括：

构建第四方物流运输网络的有向邻接表gList，有向邻接表gList中包括N个城市节点，每两个有向邻接的城市节点之间通过G条运输路径相连，每条运输路径表示使用对应的一个运输代理商在两个城市节点之间进行物流运输；

初始化所有运输路径上的信息素，在每一轮迭代中，对于任意蚂蚁个体k，初始化蚂蚁个体k位于有向邻接表gList中的起点城市节点处，根据有向邻接表gList在信息素τ以及启发式信息η的引导下、从当前所在的城市节点的可访问路径中确定目标运输路径，并沿着目标运输路径移动至下一个有向邻接的城市节点，直至蚂蚁个体k移动到达终点城市节点时完成解的构建，其中信息素τ和启发式信息η均与运输路径所需的运输花费相关且启发式信息η与运输路径所需的运输花费成反比；

在所有m个蚂蚁个体都完成解的构建而完成一轮迭代后，全局更新所有运输路径上的信息素，并执行下一轮迭代直至达到迭代终止条件，得到从起点城市节点至终点城市节点的总运输时间最短的运输路径；

其中，任意蚂蚁个体k在任意城市节点i处的可访问路径包括城市节点i与所有可访问城市节点之间除路径禁忌表edgeTabooList中的运输路径之外的所有运输路径，蚂蚁个体k在城市节点i处的可访问城市节点包括城市节点i的所有有向邻接的城市节点中除节点禁忌表nodeTabooList中的城市节点之外的所有城市节点；蚂蚁个体k在城市节点i处的节点禁忌表nodeTabooList包含与城市节点i有向邻接且不满足运输成本约束条件的城市节点，蚂蚁个体k在城市节点i处的路径禁忌表edgeTabooList包含与城市节点i相连且不满足运载量约束条件的运输路径。

其进一步的技术方案为，初始化城市节点i与城市节点j之间对应于运输代理商g的运输路径上的信息素为C^td是由弗洛伊德算法通过运输路径对应的单位运输时间和城市节点i与城市节点j之间的运输距离d_ij构造的路径长度；其中，运输路径对应的单位运输时间表示运输代理商g在城市节点i与城市节点j之间进行物流运输时每单位运输距离以及每单位运载量所需的运输时间。

其进一步的技术方案为，城市节点i与城市节点j之间对应于运输代理商g的运输路径上与运输花费成反比的启发式信息其中，为运输代理商g的运输速度，d_ij为运输路径的运输距离，是运输代理商g在城市节点i与城市节点j之间进行物流运输时的总花费。

其进一步的技术方案为，当蚂蚁个体k位于城市节点i处时，对于城市节点i的任意一个有向邻接的城市节点j，若根据城市节点i与城市节点j之间的运输距离d_ij以及最低运输成本c_{ij_min}确定得到的城市节点i与城市节点j之间的最低运输成本超过额定最大运输成本C，则确定城市节点j不满足运输成本约束条件、且包含在蚂蚁个体k在城市节点i处的节点禁忌表nodeTabooList中；

其中，城市节点i与城市节点j之间的最低运输成本c_{ij_min}是所有运输代理商在城市节点i与城市节点j之间进行物流运输时每单位运输距离以及每单位运载量所需的运输成本的最小值。

其进一步的技术方案为，当蚂蚁个体k位于城市节点i处时，对于城市节点i与任意一个有向邻接的城市节点j之间的任意一条运输路径，若运输路径对应的运输代理商的最大运载量小于待配送总运载量Q，则确定运输路径不满足运载量约束条件、且包含在蚂蚁个体k在城市节点i处的路径禁忌表edgeTabooList中。

其进一步的技术方案为，从当前所在的城市节点的可访问路径中确定目标运输路径，包括：

生成一个随机数q∈[0,1]，若q≤q₀则按照的公式、选择可访问路径中使得最大的运输路径作为目标运输路径，否则通过轮盘赌的方法从可访问路径中选择得到目标运输路径；

其中，表示运输路径上的信息素，表示运输路径上的启发式信息，运输路径表示城市节点i与城市节点j之间对应于运输代理商g的运输路径，α、β、q₀为参数，表示运输路径包含在蚂蚁个体k在城市节点i处的可访问路径构成的集合J_k(i)中。

其进一步的技术方案为，通过轮盘赌的方法从可访问路径中选择得到目标运输路径，包括：

从可访问路径中选择对应概率最大的运输路径作为目标运输路径，任意的运输路径对应的概率为：

其中，t表示蚂蚁个体k在城市节点i处的可访问路径构成的集合J_k(i)中的任意一个运输路径，τ_t是运输路径t上的信息素，η_t表示运输路径t上的启发式信息。

其进一步的技术方案为，在任意一轮迭代中对于任意一个蚂蚁个体k，方法还包括：

若蚂蚁个体k在当前的城市节点处不存在可访问路径，且蚂蚁个体k当前所在的城市节点不是终点城市节点，则控制蚂蚁个体k沿着回溯路径返回至上一个城市节点，回溯路径是蚂蚁个体k在上一个城市节点处确定得到的目标运输路径；将回溯路径添加到蚂蚁个体k在上一个城市节点处的路径禁忌表edgeTabooList中并对应更新蚂蚁个体k在上一个城市节点的可访问路径，重新确定蚂蚁个体k在上一个城市节点处的目标运输路径。

其进一步的技术方案为，全局更新所有运输路径上的信息素，包括：

将城市节点i与城市节点j之间对应于运输代理商g的运输路径上的信息素更新为ρ是参数，表示至今最优的蚂蚁个体在运输路径上释放的信息素，若运输路径包含在至今最优的运输路径中则若运输路径不包含在至今最优的运输路径中则

其进一步的技术方案为，至今最优的蚂蚁个体在运输路径上释放的信息素为：

其中，speed_best为至今最优的运输路径的总的运输时间，distance_best为至今最优的运输路径的总距离，cost_best为至今最优的运输路径的总花费，R_k为至今最优的运输路径。

本发明的有益技术效果是：

本申请公开了一种基于蚁群算法的第四方物流运输路径规划方法，该方法采用有向图的邻接表方式保存第四方物流运输网络的信息，并提出一种情景式编码方法，基于问题特征设计了启发式信息和信息素计算方式，从不同角度提高解的质量，利用蚁群算法可以在运输成本约束和运载量约束下得到运输路径，以实现物流作业整体的总运输时间最短的目标，求解难度较低且解的质量较高，可以兼顾求解质量和求解速度，路径规划效率较高。

在蚁群算法应用过程中，使用节点禁忌表和路径禁忌表来减少不可行解产生的概率，另外使用了路径回溯的方法代替初始化新解来修复不可行路径，蚁群算法效率较高。

附图说明

图1是一个实施例中的第四方物流运输路径规划方法的流程示意图。

图2是一个实例中的第四方物流运输网络的示意图。

图3是另一个实施例中的第四方物流运输路径规划方法的流程示意图。

图4是一个实验对比例中本申请的方法与其他两种算法的算法收敛图。

具体实施方式

下面结合附图对本发明的具体实施方式做进一步说明。

本申请公开了一种基于蚁群算法的第四方物流运输路径规划方法，请参加图1所示的流程图，该方法包括如下步骤：

步骤1，构建第四方物流运输网络的有向邻接表gList。本申请采用有向图的邻接表方式保存第四方物流运输网络的信息，有向邻接表gList中包括N个城市节点，每两个有向邻接的城市节点之间通过G条运输路径相连，每条运输路径表示使用对应的一个运输代理商在两个城市节点之间进行物流运输。

比如如图2示出了一个实例中的第四方物流运输网络的示意图，则对应的有向邻接表中包括7个城市节点，分别表示图2中的数值1～7的圆圈，每两个有向邻接的城市节点之间通过2条运输路径相连，比如城市节点1与城市节点2之间通过运输路径a_1,2和运输路径b_1,2相连，城市节点1与城市节点4之间通过运输路径a_1,4和运输路径b_1,4相连，以此类推。图2以不同的下标表示不同城市节点之间的运输路径，任意两个城市节点之间的运输路径a对应运输代理商A、运输路径b对应运输代理商B，则城市节点1与城市节点2之间的运输路径a_1,2表示使用运输代理商A在城市节点1与城市节点2之间进行物流运输，其他含义也是类似的。

需要说明的是，本申请以每两个有向邻接的城市节点之间通过相同数量、相同类型的运输代理商进行运输为例，实际应用过程中若数量和类型存在不同也可以沿用本申请的方法。

请参加图2所示的第四方物流运输网络的结构，当需要从一个起点城市节点向一个终点城市节点进行物流运输时，有多种路径选择，如图2中，当需要从城市节点1向城市节点5进行物流运输时，可以从城市节点1途径城市节点2到达城市节点5，或者从城市节点1途径城市节点4到达城市节点5，而在每两个城市节点之间也可以选用不同的运输代理商进行运输，比如可以通过运输路径a_1,2和运输路径a_2,5到达城市节点5，也可以通过运输路径b_1,4和运输路径a_4,5到达城市节点5，其他还有更多种可选的运输路径，当第四方物流运输网络更为复杂时，可选的运输路径更多。

在众多可选的运输路径中，需要考虑在一定的运输成本下选择最佳路径和适合的运输代理商，进而实现物流作业整体的总运输时间最短的目标，此模型优化目标定义为：

其中，表示运输代理商g在城市节点i与城市节点j之间进行物流运输时每单位运输距离以及每单位运载量所需的运输时间。表示利用运输代理商g在城市节点i与城市节点j之间进行物流运输的运输路径，当选用运输代理商g进行运输时，取1，否则取0。

此模型的约束条件包括运输成本约束条件和运载量约束条件，写为：

Q≤Q_g,g∈{1,2,…,G}

其中，d_ij表示城市节点i与城市节点j之间的运输距离，Q为需要从起点城市节点向终点城市节点运输的货物的待配送总运载量，表示利用运输代理商g在城市节点i与城市节点j之间进行物流运输的每单位运输距离以及每单位运载量所需的运输成本。C是用户能够承担的额定最大运输成本。Q_g为运输代理商g的最大运载量。

本申请基于蚁群算法，结合有向邻接表gList所记载的第四方物流运输网络的信息，可以在上述约束条件下求解得到目标模型，从而选择出运输时间最优的运输路径，具体做法如下：

步骤2，初始化所有运输路径上的信息素。

在一个实施例中，对于任意一个城市节点i与城市节点j之间对应于运输代理商g的运输路径上的信息素，将其初始化为

其中，m是蚁群算法中所有蚂蚁个体的总个数。C^td是由弗洛伊德算法通过运输路径对应的单位运输时间和城市节点i与城市节点j之间的运输距离d_ij构造的路径长度。运输路径对应的单位运输时间表示运输代理商g在城市节点i与城市节点j之间进行物流运输时每单位运输距离以及每单位运载量所需的运输时间。

步骤3，在每一轮迭代中，对于任意蚂蚁个体k，初始化蚂蚁个体k位于有向邻接表gList中的起点城市节点处，起点城市节点和终点城市节点都是有向邻接表gList中的一个城市节点，是已知的城市节点。

根据有向邻接表gList在信息素τ以及启发式信息η的引导下、从蚂蚁个体k当前所在的城市节点的可访问路径中确定目标运输路径，并沿着目标运输路径移动至下一个有向邻接的城市节点，直至蚂蚁个体k移动到达终点城市节点时完成解的构建。

其中，信息素τ和启发式信息η均与运输路径所需的运输距离、运输速度和相应花费相关且启发式信息η与运输路径所需的运输花费成反比。具体的，运输路径上与运输花费成反比的启发式信息为运输代理商g的运输速度，d_ij为运输路径的运输距离，是运输代理商g在城市节点i与城市节点j之间进行物流运输时的总花费。

任意蚂蚁个体k在任意城市节点i处的可访问路径包括城市节点i与所有可访问城市节点之间除路径禁忌表edgeTabooList中的运输路径之外的所有运输路径。而蚂蚁个体k在城市节点i处的可访问城市节点包括城市节点i的所有有向邻接的城市节点中除节点禁忌表nodeTabooList中的城市节点之外的所有城市节点。也即，蚂蚁个体k在城市节点i处的可访问路径是城市节点i连接的所有运输路径中排除掉nodeTabooList中的城市节点以及排除掉edgeTabooList中的运输路径后剩余的运输路径。

蚂蚁个体k在城市节点i处的节点禁忌表nodeTabooList包含与城市节点i有向邻接且不满足运输成本约束条件的城市节点。具体的：当蚂蚁个体k位于城市节点i处时，对于城市节点i的任意一个有向邻接的城市节点j，若根据城市节点i与城市节点j之间的运输距离d_ij以及最低运输成本c_{ij_min}确定得到的城市节点i与城市节点j之间的最低运输成本超过额定最大运输成本C。则确定城市节点j不满足运输成本约束条件、且包含在蚂蚁个体k在城市节点i处的节点禁忌表nodeTabooList中。城市节点i与城市节点j之间的最低运输成本c_{ij_min}是所有运输代理商在城市节点i与城市节点j之间进行物流运输时每单位运输距离以及每单位运载量所需的运输成本的最小值，则当d_ij×c_{ij_min}×Q>C时，确定城市节点j不满足运输成本约束条件。

蚂蚁个体k在城市节点i处的路径禁忌表edgeTabooList包含与城市节点i相连且不满足运载量约束条件的运输路径。具体的：当蚂蚁个体k位于城市节点i处时，对于城市节点i与任意一个有向邻接的城市节点j之间的任意一条运输路径，若运输路径对应的运输代理商的最大运载量小于待配送总运载量Q，则确定运输路径不满足运载量约束条件、且包含在蚂蚁个体k在城市节点i处的路径禁忌表edgeTabooList中。

比如，基于图2所示的第四方物流运输网络，当蚂蚁个体k在任意城市节点1处时，假设其节点禁忌表nodeTabooList中包括城市节点3，路径禁忌表edgeTabooList中包括运输路径b_1,4。则此时，城市节点i的所有有向邻接的城市节点包括城市节点2、城市节点3和城市节点4，排除掉节点禁忌表nodeTabooList中的城市节点3后，确定此时所有可访问城市节点包括城市节点2和城市节点4。城市节点1与城市节点2和城市节点4之间的所有运输路径包括a_1,2、b_1,2、a_1,4、b_1,4，排除掉edgeTabooList中包含的b_1,4，就能确定此时蚂蚁个体k在城市节点1处的可访问路径包括a_1,2、b_1,2、a_1,4。

在蚂蚁个体k位于每个城市节点i处时，都按照上述方法确定此时的节点禁忌表nodeTabooList和路径禁忌表edgeTabooList，继而确定相应的可访问路径，请参考图3所示的流程图，然后从中选取目标运输路径以移动到下一个城市节点，确定目标运输路径的方法包括：

生成一个随机数q∈[0,1]，若q≤q₀则按照的公式、选择可访问路径中使得最大的运输路径作为目标运输路径，表示运输路径上的信息素，表示运输路径上的启发式信息，运输路径表示城市节点i与城市节点j之间对应于运输代理商g的运输路径，α、β、q₀为参数，表示运输路径包含在蚂蚁个体k在城市节点i处的可访问路径构成的集合J_k(i)中。

若q>q₀，则通过轮盘赌的方法从可访问路径中选择得到目标运输路径。在通过轮盘赌的方法从可访问路径中选择得到目标运输路径时，可访问路径中选择对应概率最大的运输路径作为目标运输路径，任意的运输路径对应的概率为：

其中，t表示蚂蚁个体k在城市节点i处的可访问路径构成的集合J_k(i)中的任意一个运输路径，τ_t是运输路径t上的信息素，η_t表示运输路径t上的启发式信息。

请参考图3，在蚂蚁个体k逐步移动的过程中，若蚂蚁个体k在当前的城市节点处不存在可访问路径，且蚂蚁个体k当前所在的城市节点是终点城市节点，则表示已经完成了解的构建。若蚂蚁个体k在当前的城市节点处不存在可访问路径，且蚂蚁个体k当前所在的城市节点不是终点城市节点，则控制蚂蚁个体k沿着回溯路径返回至上一个城市节点，回溯路径是蚂蚁个体k在上一个城市节点处确定得到的目标运输路径。将回溯路径添加到蚂蚁个体k在上一个城市节点处的路径禁忌表edgeTabooList中并对应更新蚂蚁个体k在上一个城市节点的可访问路径，重新确定蚂蚁个体k在上一个城市节点处的目标运输路径，继而重新从上一个城市节点处移动至重新确定的城市节点处。edgeTabooList的一部分作用是保证蚂蚁个体能够避开不符合运输商能够承担的最大运输量约束的路径，提高算法的效率，另一部分作用是负责记录失败的路径信息，当蚂蚁个体所在的最后的一个城市节点不是期望的终点城市节点时，通过路径的回溯可以回到上一城市节点来修复不可行路径。

步骤4，在所有m个蚂蚁个体都完成解的构建而完成一轮迭代后，全局更新所有运输路径上的信息素，并执行下一轮迭代直至达到迭代终止条件，得到从起点城市节点至终点城市节点的总运输时间最短的运输路径。

具体的迭代过程为，在开始一轮迭代时，初始化k＝0，通过上述方法控制蚂蚁个体k移动到达终点城市节点时完成解的构建，并令k＝k+1，若满足k<m，则利用更新后的k对下一个蚂蚁个体进行解的构建。若不满足k<m，则确定完成一轮迭代，并全局更新所有运输路径上的信息素执行下一轮迭代直至达到迭代终止条件。迭代终止条件可以预先设定，比如设定迭代的总轮数等等。

在全局更新所有运输路径上的信息素时，将城市节点i与城市节点j之间对应于运输代理商g的运输路径上的信息素更新为ρ是参数，表示至今最优的蚂蚁个体在运输路径上释放的信息素，若运输路径包含在至今最优的运输路径中则若运输路径不包含在至今最优的运输路径中则也即全局信息素的更新仅在至今最优的路径的边上进行更新。其中，至今最优的蚂蚁个体在运输路径上释放的信息素为：

其中，speed_best为至今最优的运输路径的总的运输时间，distance_best为至今最优的运输路径的总距离，cost_best为至今最优的运输路径的总花费，R_k为至今最优的运输路径。

为了说明本申请方法的有效性，本申请设置如下实例进行实验数据对比：以算法一为本申请提供的方法(IACO)。算法二为2021年发表在智能系统学报中的《模糊粒子群优化算法的第四方物流运输时间优化》一文中的引入收敛因子和隶属度函数的模糊粒子群优化算法(CFPSO)。算法三为经典的遗传算法(GA)。

以算法二的论文提供的三个不同规模的算例为实验数据，分别利用算法一、二、三进行求解，这三个算例包含了不同的数据规模与成本约束，能够很好的验证求解算法的性能，三个算例的具体内容如下表所示：

(1)、算例一的城市节点间距离如下表所示，且每两个有向邻接的城市节点之间包括运输代理商分别为代理商1和代理商2，代理商1的运输成本为0.16RMB/吨、运输速度为80km/h。代理商2的运输成本为0.08RMB/吨、运输速度为40km/h。

城市	杭州	南京	上海	南通	泰州	淮安
							杭州	0	330	195	—	—	—
南京	330	0	248	306	—	355
							上海	195	248	0	103	217	—
南通	—	306	103	0	160	354
							泰州	—	—	217	160	0	193
淮安	—	355	—	354	193	0

(2)算例二的城市节点间距离如下表所示，且每两个有向邻接的城市节点之间包括运输代理商分别为代理商1、代理商2和代理商3，代理商1的运输成本为0.16RMB/吨、运输速度为80km/h。代理商2的运输成本为0.08RMB/吨、运输速度为40km/h。代理商3的运输成本为1.5RMB/吨、运输速度为750km/h。

(3)算例三的城市节点间距离如下表所示，且每两个有向邻接的城市节点之间包括运输代理商分别为代理商1、代理商2和代理商3，代理商1的运输成本为0.16RMB/吨、运输速度为80km/h。代理商2的运输成本为0.08RMB/吨、运输速度为40km/h。代理商3的运输成本为1.5RMB/吨、运输速度为750km/h。

以上三个算例的运输网络分别有6、12和18个城市节点，每个算例中的起点城市节点和终点城市节点分别是算例中第一个城市节点和最后一个城市节点。算法一、二、三的参数通过控制变量方法进行设置，具体的参数如下表。

CFPSO与GA适应度函数的构建都采用罚函数法，可表示为：

其中w₁、w₂、w₃都为惩罚项系数。由于IACO采用了修复策略，所以种群中的解都为可行解。三个算法都部署在同一实验环境下，且不同问题规模和成本约束都作为一个独立的问题运行二十次，三种算法的实验结果对比如下：

(1)三种算法在算例一上的实验结果对比如下

(2)三种算法在算例二上的实验结果对比如下

(3)三种算法在算例三上的实验结果对比如下

实验结果表明，本申请的方法在不同的数据规模和成本约束下均能取得较对比算法更好的结果。随着数据规模越大时，求解难度也显著增加，但本申请方法取得的结果较小规模相比却领先其他算法越多，说明了本申请的方法有较强的处理大规模问题的能力。观察实验结果可以发现，本申请方法在各个问题规模和成本约束下的最小和平均运输时间差异很小说明本发明方法具有较强的稳定性。

为了进一步的分析本申请方法与其他对比算法的性能，基于算例3且成本约束为300000元问题下，三种算法的迭代次数都设置为500代的算法收敛图如图4所示，由图4可以看出三个算法随着迭代都可以很好的收敛。而本申请方法较其它对比算法收敛快且曲线收敛平稳，说明具有较强的算法收敛能力和稳定性。且本申请方法的曲线在其它算法曲线下方，说明其解的质量稳定高于其他算法。

以上所述的仅是本申请的优选实施方式，本发明不限于以上实施例。可以理解，本领域技术人员在不脱离本发明的精神和构思的前提下直接导出或联想到的其他改进和变化，均应认为包含在本发明的保护范围之内。

Claims (9)

1.一种基于蚁群算法的第四方物流运输路径规划方法，其特征在于，所述方法包括：

构建第四方物流运输网络的有向邻接表gList，所述有向邻接表gList中包括N个城市节点，每两个有向邻接的城市节点之间通过G条运输路径相连，每条运输路径表示使用对应的一个运输代理商在两个城市节点之间进行物流运输；

初始化所有运输路径上的信息素，包括初始化城市节点i与城市节点j之间对应于运输代理商g的运输路径

上的信息素为

C^td是由贪婪算法通过运输路径

对应的单位运输时间

和城市节点i与城市节点j之间的运输距离d_ij构造的路径长度；其中，运输路径

对应的单位运输时间

表示运输代理商g在城市节点i与城市节点j之间进行物流运输时每单位运输距离以及每单位运载量所需的运输时间；在每一轮迭代中，对于任意蚂蚁个体k，初始化蚂蚁个体k位于所述有向邻接表gList中的起点城市节点处，根据所述有向邻接表gList在信息素τ以及启发式信息η的引导下、从当前所在的城市节点的可访问路径中确定目标运输路径，并沿着所述目标运输路径移动至下一个有向邻接的城市节点，直至蚂蚁个体k移动到达终点城市节点时完成解的构建，其中，信息素τ和启发式信息η均与运输路径所需的运输花费相关且启发式信息η与运输路径所需的运输花费成反比；

在所有m个蚂蚁个体都完成解的构建而完成一轮迭代后，全局更新所有运输路径上的信息素，并执行下一轮迭代直至达到迭代终止条件，得到从所述起点城市节点至所述终点城市节点的总运输时间最短的运输路径；

其中，任意蚂蚁个体k在任意城市节点i处的可访问路径包括城市节点i与所有可访问城市节点之间除路径禁忌表edgeTabooList中的运输路径之外的所有运输路径，蚂蚁个体k在城市节点i处的可访问城市节点包括城市节点i的所有有向邻接的城市节点中除节点禁忌表nodeTabooList中的城市节点之外的所有城市节点；蚂蚁个体k在城市节点i处的节点禁忌表nodeTabooList包含与城市节点i有向邻接且不满足运输成本约束条件的城市节点，蚂蚁个体k在城市节点i处的路径禁忌表edgeTabooList包含与城市节点i相连且不满足运载量约束条件的运输路径。

2.根据权利要求1所述的方法，其特征在于，城市节点i与城市节点j之间对应于运输代理商g的运输路径

上与运输花费成反比的启发式信息

其中，

为运输代理商g的运输速度，d_ij为运输路径

的运输距离，

是运输代理商g在城市节点i与城市节点j之间进行物流运输时的总花费。

3.根据权利要求1所述的方法，其特征在于，当蚂蚁个体k位于城市节点i处时，对于城市节点i的任意一个有向邻接的城市节点j，若根据城市节点i与城市节点j之间的运输距离d_ij以及最低运输成本c_{ij_min}确定得到的城市节点i与城市节点j之间的最低运输成本超过额定最大运输成本C，则确定城市节点j不满足运输成本约束条件、且包含在蚂蚁个体k在城市节点i处的节点禁忌表nodeTabooList中；

其中，城市节点i与城市节点j之间的最低运输成本c_{ij_min}是所有运输代理商在城市节点i与城市节点j之间进行物流运输时每单位运输距离以及每单位运载量所需的运输成本的最小值。

4.根据权利要求1所述的方法，其特征在于，当蚂蚁个体k位于城市节点i处时，对于城市节点i与任意一个有向邻接的城市节点j之间的任意一条运输路径，若所述运输路径对应的运输代理商的最大运载量小于待配送总运载量Q，则确定所述运输路径不满足运载量约束条件、且包含在蚂蚁个体k在城市节点i处的路径禁忌表edgeTabooList中。

5.根据权利要求1所述的方法，其特征在于，所述从当前所在的城市节点的可访问路径中确定目标运输路径，包括：

生成一个随机数q∈[0,1]，若q≤q₀则按照

的公式、选择可访问路径中使得

最大的运输路径作为所述目标运输路径，否则通过轮盘赌的方法从可访问路径中选择得到所述目标运输路径；

其中，

表示运输路径

上的信息素，

表示运输路径

上的启发式信息，运输路径

表示城市节点i与城市节点j之间对应于运输代理商g的运输路径，α、β、q₀为参数，

表示运输路径

包含在蚂蚁个体k在城市节点i处的可访问路径构成的集合J_k(i)中。

6.根据权利要求5所述的方法，其特征在于，所述通过轮盘赌的方法从可访问路径中选择得到所述目标运输路径，包括：

从可访问路径中选择对应概率最大的运输路径作为所述目标运输路径，任意的运输路径

对应的概率

为：

其中，t表示蚂蚁个体k在城市节点i处的可访问路径构成的集合J_k(i)中的任意一个运输路径，τ_t是运输路径t上的信息素，η_t表示运输路径t上的启发式信息。

7.根据权利要求1所述的方法，其特征在于，在任意一轮迭代中对于任意一个蚂蚁个体k，所述方法还包括：

若蚂蚁个体k在当前的城市节点处不存在可访问路径，且蚂蚁个体k当前所在的城市节点不是所述终点城市节点，则控制蚂蚁个体k沿着回溯路径返回至上一个城市节点，所述回溯路径是蚂蚁个体k在上一个城市节点处确定得到的目标运输路径；将所述回溯路径添加到蚂蚁个体k在上一个城市节点处的路径禁忌表edgeTabooList中并对应更新蚂蚁个体k在上一个城市节点的可访问路径，重新确定蚂蚁个体k在上一个城市节点处的目标运输路径。

8.根据权利要求1所述的方法，其特征在于，所述全局更新所有运输路径上的信息素，包括：

将城市节点i与城市节点j之间对应于运输代理商g的运输路径

上的信息素更新为

ρ是参数，

表示至今最优的蚂蚁个体在运输路径

上释放的信息素，若运输路径

包含在至今最优的运输路径中则

若运输路径

不包含在至今最优的运输路径中则

9.根据权利要求8所述的方法，其特征在于，至今最优的蚂蚁个体在运输路径

上释放的信息素为：

其中，s_peed_best为至今最优的运输路径的总的运输时间，distance_best为至今最优的运输路径的总距离，cost_best为至今最优的运输路径的总花费，R_k为至今最优的运输路径。

Images