CN107886196A - 一种用于货物取送的单车调度方法 - Google Patents

Abstract

本发明公开一种用于货物取送的单车调度方法，获取货物取送车辆调度应用场景的系统特征数据，创建取送任务，将取送问题的系统特征数据处理为符合组合结构构建要求的形式；创建用于构建组合结构的约束及其构建组合结构的方法，并根据其对组合结构的影响形式进行归类，随机地构造一个排序结构，利用基于排序结构约束的排序结构构建方法对当前排序结构进行局部调整使其满足约束，在排序结构基础上利用分组结构约束逐一构建每个分组，并利用整体结构约束进行检测和局部的组合结构调整；将其置于元启发式方法，对方案进行迭代进化，根据终止规则停止迭代，输出优化的调度计划。本发明达到实时性和可移植性兼备的目的，满足在实际应用场景中的使用需求。

Description

一种用于货物取送的单车调度方法

技术领域

本发明涉及智能调度技术领域，具体为一种用于货物取送的单车调度方法。

背景技术

货物取送存在于社会经济系统的诸多领域，例如，城市物流配送系统、分拣系统、货柜 自动化装箱系统和铁路车站取送调车系统等。特别地，随着科技的进步和发展，它将面向众 多新的应用领域，如无人机配送调度、无人配送小车调度和无人驾驶出租车调度等。尽管该 问题由来已久，且已有诸多较为成熟的方法，但大多倾向于理论研究，当具体到应用层面时， 由于使用场景的特征各异，使得方法在效率或可移植性上受到诸多制约，进而使得货物取送 的车辆调度方法被持续关注。

车辆调度方法主要分为精确方法和启发式方法。精确方法主要用于寻求问题的解析解， 包括分枝定价法、列生成方法等，其特点在于将问题转换为线性规划问题进行求解。然而， 当具体到应用场景时，往往呈现出大量的非线性特征，尽管通常能将非线性约束转化为线性 约束，但这将不可避免地引入新的决策变量，进而增大搜索空间，使得调度方法的效率受到 制约。启发式方法是一种基于直观经验的方法，它能在可接受的开销基础上得到满意的调度 结果，尽管不能确定它与最优调度之间的偏离程度，但它可较好地解决精确方法中由于线性 转换所产生的搜索空间过大问题，同时，更容易将人的直观经验和知识植入方法中，使得方 法更为合理，此类方法包括蚁群算法、遗传算法、模拟退火算法和神经网络算法等，可见， 启发式方法是一种更倾向于实时性的方法。尽管启发式方法在效率方面往往更具优势，但它 无法避免对具体特征的依赖，这些特征由应用场景所决定，不同场景中系统特征的微小区别 往往使得调度方法截然不同，这使得启发式方法的可扩展性受到制约，进而限制其适用范围。 实践表明，启发式方法效率的高低并非完全取决于元启发式方法的迭代或进化机制，例如蚁 群算法的进化机制为模拟蚂蚁觅食行为，遗传算法的进化机制为模拟基因交叉和变异，而更 多取决于对实际应用场景的编码方式，以及基于这种编码方式和系统特征的特定启发式方法 的设计，换言之，它主要取决于针对可行调度结果构建的启发式方法。

发明内容

针对上述问，本发明的目的在于提供能达到实时性和可移植性兼备的目的，从而满足在 实际应用场景中的使用需求，且更为高效的用于货物取送的单车调度方法。技术方案如下：

一种用于货物取送的单车调度方法，包括以下步骤：

步骤1：获取货物取送车辆调度应用场景的系统特征数据，包括地理信息数据、配送车辆 信息数据、货件信息数据；

步骤2：对货物取送车辆调度应用场景的系统特征数据进行预处理，创建取送任务，以取 送任务为基本单元，将其与系统特征数据进行关联，进而将取送问题的系统特征数据处理为 符合组合结构构建要求的形式；

步骤3：对货物取送车辆调度应用场景的系统特征进行分析，根据系统特征创建用于构建 组合结构的约束及其构建组合结构的方法，并根据其对组合结构的影响形式进行归类，包括 局部组合结构约束、整体组合结构约束和超组合结构约束；

步骤4：以取送任务作为取送方案组合结构的基本构造单位，采用构造可行组合结构的方 式；可行组合结构的构造方法为基于组合结构的约束条件及其分类，随机地构造一个排序结 构，利用基于排序结构约束的排序结构构建方法对当前排序结构进行局部调整使其满足约束， 在排序结构基础上利用分组结构约束逐一构建每个分组，在每个分组创建完成后利用整体结 构约束进行检测和局部的组合结构调整，若组合结构构造过程中存在不满足整体结构约束， 并且不能基于整体结构约束的组合结构构建方法进行结构的局部调整，则重复执行可行组合 结构的构造方法，直至得到一个完整的、可行的组合结构；

步骤5：将步骤4置于元启发式方法，对方案进行迭代进化，最终根据终止规则停止迭代， 输出优化方案。

进一步的，所述地理信息数据包括仓库位置、客户位置、客户位置分区、各位置之间的 连接路径、各路径弧段的平均行驶速度信息；配送车辆信息数据包括可执行任务车辆数、车 辆的可执行任务时间范围信息、车型信息；货件信息数据包括货件的运输量、货件对应的取 送起点和终点位置、货件在其起点和终点位置上的时间窗，以及在客户处进行中间加工的时 间。

更进一步的，所述步骤2具体包括：

步骤21：对货物取送车辆调度应用场景中所涉及的每个位置赋予索引，将仓库的索引赋 值为0，客户的索引从1开始赋值；对取送问题中包含的每个货件赋予索引；对执行取送任 务的每个车辆赋予索引；

步骤22：根据货件所对应的取送起点和终点位置数据为每个货件创建取送任务，包括从 仓库送至客户处的送件类型，从客户处取回仓库的取件类型，从客户处取出并送往另一客户 处的调移类型；其中调移任务被拆分为调移取件任务和调移送件任务；为每项取送任务赋予 唯一的索引；

步骤23：基于货物取送车辆调度应用场景中的特征数据和已创建的取送任务，创建满足 组合结构构建的已知数据变量。

更进一步的，所述步骤23具体包括：

步骤231：创建取送任务及其取送货件的关联变量c_i；其中，i为取送任务的索引，c_i为取 送任务i对应的货件；

步骤232：创建取送任务及其起点和终点位置的关联变量o_i和d_i；其中，o_i为取送任务i的 起点位置索引，d_i为取送任务i的终点位置索引；

步骤233：创建取送任务及其直接关联客户所在位置的关联变量l_i；其中，l_i为取送任务 i的直接关联客户的索引；对于取件类型任务i，有l_i＝o_i；对于取件类型任务i，有l_i＝d_i；

步骤234：创建取送任务在起点和终点位置上的时间窗变量和其中， 和分别为取送任务i在其起点位置o_i上可执行任务的最早和最晚时间，和分别为 取送任务i在其终点位置d_i上可执行任务的最早和最晚时间；对于不涉及最早时间或的 取送任务，有或对于不涉及最晚时间或的取送任务，有或其中M为一个充分大的正数；

步骤235：创建各类型取送任务的集合，包括送件类型任务集合N^-、取件类型任务集合N⁺， 所有货件最后一项取送任务集合N_E，调移送件任务集合和调移取件任务集合其中， 送件类型任务集合N^-为所有送件任务和调移送件任务的集合，取件类型任务集合N⁺为所有 取件任务和调移取件任务的集合；

步骤236：创建取送任务及其运输量的关联变量w_i；其中，w_i为取送任务i所取送货件的 运输量；

步骤237：创建包括仓库与客户，和客户与客户位置之间的距离变量td_ll′，及行驶时间变 量t_ll′；其中，l和l′分别为行驶路径弧段的起点和终点位置的索引，特别地，若两个位置之间 无直接通路，则有td_ll′和t_ll′＝M，M为一个充分大的正数；

步骤238：创建取送任务及其对应货件在客户处的中间加工时间的关联变量T_i，其中，i为 取送任务的索引，T_i为取送任务i所对应货件在客户处的中间加工时间；

步骤239：创建执行取送任务的车辆及其最大装载量的关联变量Q，其中，Q为车辆的最 大装载量。

更进一步的，所述步骤3具体包括：

步骤31：创建局部组合结构约束，包括排序结构约束和分组结构约束，该约束在其他局 部组合结构固定的条件下仅通过改变一种局部组合结构便可满足；

步骤32：创建需改变多种局部组合结构才能满足的整体组合结构约束；建立货件取送时 间窗上界约束，并将其归类于整体组合结构约束，货件取送时间窗上界约束要求每项取送任 务的执行时间都不超过货件取送时间窗上界；基于货件取送时间窗上界约束的组合结构构建 方法具体为：根据外部输入的最大分组范围[k⁺，k^-]，即当前分组中首、末两项取送任务在当 前排序结构下的执行顺序序号，从首项取送任务开始，依排序结构中顺序序号的升序检索 当前分组中的取送任务i_r；

步骤33：创建保证在固定组合结构下取送方案唯一性的超组合结构约束，包括货件取送 时间窗下界约束，货件中间加工完成时间约束和离开客户时间约束。

更进一步的，所述步骤31具体包括：

步骤311：创建排序结构约束，方法为：对同一货件的取送任务之间客观上存在优先序关 系的任务，对取送任务的排序结构加以约束，建立同一货件的取送优先序约束，并将其归类 于排序结构约束；同一货件的取送优先序约束要求在取送任务的排序结构中，同一货件的取 送任务执行顺序满足基于同一货件的取 送优先序约束的排序结构构建方法具体为：当一个取送任务序列中存在某个货件c_i的取送任 务不满足优先序约束时，采用顺序互换的方式进行排序结构的局部调整，使取送任务序列满 足同一货件的取送优先序约束；

步骤312：将一个取送任务分组定义为车辆从仓库出发并回到仓库的过程中所执行的取送 任务集合，创建仅通过改变取送任务的分组结构就能被满足的分组结构约束。

更进一步的，所述步骤312具体包括：

步骤3121：建立相异分区取送任务不同组约束，并将其归类于分组结构约束；相异分区 取送任务不同组约束要求在取送任务的分组中，当顺序相邻的两项取送任务所对应的客户位 置之间无直接通路时，将这两项取送任务依次分配至相邻的两个分组；基于相异分区取送任 务不同组约束的分组结构构建方法为：在一个固定的排序结构下构建一个分组，根据外部输 入的最大分组范围[k⁺，k^-]，即当前分组中首、末两项取送任务在当前排序结构下的执行顺序 序号，从首项取送任务开始，依排序结构中顺序序号的升序检索当前分组中顺序相邻的两 项取送任务i_r和i_r+1所对应的直接客户的位置和之间是否存在直接通路，其中r≥k⁺且r+1≤k^-，如果不存在直接通路，即则将取送任务i_r作为当前分组的最后一项 取送任务，并将当前分组的末项取送任务的执行顺序序号更新为r，返回当前分组中末项取送 任务的执行顺序序号r，其中，i_r为在当前固定排序结构下第r项取送任务的索引；

步骤3122：建立车辆能力约束，并将其归类于分组结构约束；车辆能力约束要求在一个 固定排序结构下执行每个分组的取送任务过程中，车辆负载不超过其最大装载量。

更进一步的，所述步骤3122具体包括：

步骤31221：建立车辆始发能力约束；车辆始发能力约束要求在一个固定排序结构下的取 送任务分组中，起点为仓库和终点为客户的送件类型任务所运送货件的总量不大于任务执行 车辆的最大装载量；基于车辆始发能力约束的分组结构构建方法为：用于在一个固定的排序 结构下构建一个分组，根据外部输入的最大分组范围[k⁺，k^-]，即当前分组中首、末两项取送 任务在当前排序结构下的执行顺序序号，从首项取送任务开始，依排序结构中顺序序号升 序检索当前分组中的取送任务i_r，对起点为仓库和终点为客户位置的送件类型任务所运送货件的运输量进行累加，若出现且则将取送任务i_k作为当前分组的最后一项取送任务，将作为车辆始发负载量，并 将当前分组的末项取送任务i_k的执行顺序序号更新为k，其中，k+1≤k^-，i_r为在当前固定 排序结构下第r项取送任务的索引，σ_i为符号变量，当取送任务i为送件类型任务时，σ_i取值为 -1；为取件类型任务时，σ_i取值为1；

步骤31222、建立车辆在途能力约束；车辆在途能力约束要求在取送任务的分组中，车辆 执行每项取送任务后的负载量不大于车辆的最大装载量；基于车辆在途能力约束的分组结构 构建方法为：用于在一个固定的排序结构下构建一个分组，根据最大分组范围[k⁺，k^-]，即当 前分组中首、末两项取送任务在当前排序结构下的执行顺序序号，从首项取送任务开始， 依排序结构中顺序序号的升序检索当前分组中的取送任务i_r，利用公式计算得到当前取送任务i_k执行完毕后的负载量其中，k≤k^-，i_r为在当前固定排序结构下第r项取送任务的索引，σ_i为符号变量，当取送任务i为送件类型任务 时，σ_i取值为-1；为取件类型任务时，σ_i取值为1；根据的取值情况，更新车辆始发负载量q⁰和当前分组的末项取送任务的执行顺序序号k^-。

更进一步的，所述步骤31222还包括：

步骤312221：如果q_ik>Q，转步骤312223，否则进入下一步；

步骤312222：令k＝k+1，如果k≤k^-，则计算q_ik，转步骤312221；否则，将取送任 务作为当前分组的最后一项取送任务，返回当前分组中末项取送任务的执行顺序序号k^-；

步骤312223：在范围(k，k^-]内按降序检索执行顺序r所对应的取送任务i_r，如果不存在起 点为仓库和终点为客户位置的送件类型任务i_r，则将取送任务i_k-1作为当前分组的最后一项取 送任务，并将当前分组的末项取送任务的执行顺序序号更新为k-1，返回当前分组中末项取 送任务的执行顺序序号k-1；

步骤312224：在范围(k，k^-]内按降序检索执行顺序r所对应的取送任务i_r，将首个出现的 起点为仓库和终点为客户位置的送件类型任务i_r所对应运输量从q⁰中扣除，即有将当前分组的末项取送任务的执行顺序序号k^-更新为r-1；计算q_ik，转步骤312221。

更进一步的，所述步骤32中对取送任务i_r的处理方法包括：

步骤321：如果或转步骤323；

步骤322：令r＝r+1，如果r≤k^-，转步骤321；否则，将取送任务作为当前分组的最后一项取送任务，返回局部调整成功指示器LT＝true和分组序号指示器g＝-1，此处g取值为-1表示不需对已构建的排序结构和分组结构进行局部调整；

步骤323：在范围[0，r)内随机选择一个执行顺序r′，要求取送任务i_r′和i_r所对应的货件不 同当时有当时有如果不存在这 样的取送任务i_r′，则不再基于当前已构建的局部组合结构继续进行后续操作，返回局部调整 成功指示器LT＝false和分组序号指示器g＝-1；

步骤324：交换取送任务i_r′和i_r在当前排序结构中的顺序；

步骤325：若当前排序结构中货件和的取送任务不满足同一货件的取送优先序约束， 则采用基于同一货件的取送优先序约束的排序结构构建方法进行取送任务执行顺序的局部调 整；

步骤326：返回局部调整成功指示器LT＝true和分组序号指示器g＝g_r′，其中，g_r′为执 行顺序为r′的取送任务当前所在分组的分组序号。

更进一步的，所述步骤33具体包括：

步骤331：建立货件取送时间窗下界约束，并将其归类于超组合结构约束；货件取送时间 窗下界约束要求当车辆到达取送任务i的起点时间时车辆需等待至时间之后开始 执行取送任务i，当车辆到达取送任务i的终点时间时，车辆需等待至时间之后开 始执行取送任务i；

步骤332：建立货件中间加工完成时间约束，并将其归类于超组合结构约束；货件中间加 工完成时间约束要求当送件类型任务i将货件c_i送至客户处l_i，经时长为T_i的中间加工后，由取 件类型任务j取出送至另一地点，若车辆到达取件类型任务j起点的时间车辆需等待至时间之后开始执行任务j；

步骤333：建立离开客户时间约束，并将其归类于超组合结构约束；离开客户时间约束要 求，对于送件类型任务i，当货件c_i被送入客户处l_i后，车辆从客户处l_i离开的时间为 对于取件类型任务i，当货件c_i被送入客户处l_i后，车辆从客户处l_i离开的 时间为其中，ta_i和tl_i分别为车辆到达和离开l_i的时间，为车辆可 将货件送入客户处l_i的最早时间，th_i为货件c_i的中间加工完毕时间。

更进一步的，所述步骤4具体包括：

步骤41：随机生成一个以取送任务为基本单元的有序序列seq；

步骤42：使用基于同一货件的取送优先序约束的排序结构构建方法对序列seq中取送作业 进行执行顺序的局部调整，调整后的序列seq满足同一货件的取送优先序约束；

步骤43：依次构造每个分组并基于整体组合结构约束进行局部调整。

更进一步的，所述步骤43具体包括：

步骤431：初始化当前分组序号g＝0，初始化当前分组首项和末项取送任务在序列seq中 的执行顺序序号和其中，n为取送任务总数；

步骤432：根据和确定的最大分组范围，使用基于相异分区取送任务不同组约束的 分组结构构建方法更新序列seq中分组g的末项取送任务在当前序列seq中的顺序序号

步骤433：根据和确定的最大分组范围，使用基于车辆始发能力约束的分组结构构 建方法更新序列seq中分组g的末项取送任务在当前序列seq中的顺序序号使用基于车辆 在途能力约束的分组结构构建方法更新序列seq中分组g的末项取送任务在当前序列seq中的 顺序序号

步骤434：根据和确定的最大分组范围，使用基于货件取送时间窗上界约束的组合 结构构建方法对分组中不满足货件取送时间窗上界约束的局部组合结构进行局部调整，根据 返回结果进行相应操作。

更进一步的，所述步骤434具体包括：

步骤4341：如果局部调整成功指示器LT_g＝false，则转步骤41；

步骤4342：如果基于货件取送时间窗上界约束的组合结构构建方法返回的分组序号指示 器g′的取值为-1，则令g＝g+1、和转步骤432；

步骤4343：对于基于货件取送时间窗上界约束的组合结构构建方法返回的分组序号指示 器的取值g′，如果g′≤1，则令g＝0、和否则，令g＝g′-1、和转步骤432。

本发明的有益效果是：本发明针对货物取送车辆调度方法的发展现状、存在的问题，达 到实时性和可移植性兼备的目的，从而满足在实际应用场景中的使用需求，大大提到了货物 的取送效率，就要很大的实用性。

附图说明

图1为本发明用于货物取送的单车调度方法的流程图。

图2为本发明基于同一货件的取送优先序约束的排序结构构建方法实施示意图。

图3为本发明基于相异分区取送任务不同组约束的分组结构构建方法实施示意图。

图4为本发明基于车辆始发能力约束的分组结构构建方法实施示意图。

图5为本发明基于车辆始在途能力约束的分组结构构建方法实施示意图。

图6基于货件取送时间窗上界约束的组合结构构建方法实施示意图

图7.1为本发明货物取送单车调度计划的组合结构构造方法中随机生成的取送任务的有序序 列示意图。

图7.2本发明货物取送单车调度计划的组合结构构造方法中施以基于同一货件的取送优先序 约束的排序结构构建方法示意图。

图7.3本发明货物取送单车调度计划的组合结构构造方法中取送任务序列seq的分组0示意 图。

图7.4本发明货物取送单车调度计划的组合结构构造方法中取送任务序列seq的分组1示意 图。

图7.5本发明货物取送单车调度计划的组合结构构造方法中取送任务序列seq的分组2示意 图。

图7.6本发明货物取送单车调度计划的组合结构构造方法中最终可行组合结构示意图。

具体实施方式

下面结合附图和具体实施例对本发明做进一步详细说明。

本发明公开了一种用于货物取送的单车调度方法，其技术方案包括如下步骤：

步骤1、获取货物取送车辆调度应用场景的系统特征数据，包括地理信息、配送车辆信息、 货件信息和客户仓容量信息，地理信息数据包括仓库位置、客户位置、客户位置分区、各位 置之间的连接路径、各路径弧段的平均行驶速度信息，配送车辆信息数据包括可执行任务车 辆数、车辆的可执行任务时间范围信息、车型信息，货件信息数据包括货件的运输量、货件 对应的取送起点和终点位置、货件在其起点和终点位置上的时间窗，在客户处进行中间加工 的时间。

步骤2、对货物取送车辆调度应用场景的系统特征数据进行预处理，创建取送任务，以取 送任务为基本单元，将其与系统特征数据进行关联，进而将取送问题的系统特征数据处理为 符合组合结构构建要求的形式。

步骤21、对货物取送车辆调度应用场景中所涉及的每个位置赋予索引，将仓库的索引赋 值为0，客户的索引从1开始赋值；对取送问题中包含的每个货件赋予索引；对执行取送任 务的每个车辆赋予索引；

步骤22、根据货件所对应的取送起点和终点位置数据为每个货件创建取送任务，包括从 仓库送至客户处的送件类型，从客户处取回仓库的取件类型，从客户处取出并送往另一客户 处的调移类型，其中，调移任务被拆分为调移取件任务和调移送件任务，调移取件任务指调 移任务中从起点客户取出货件，调移送件指将从起点客户取出的货件送至终点客户的任务； 为每项取送任务赋予唯一的索引；

步骤23、基于货物取送车辆调度应用场景中的特征数据和已创建的取送任务，创建满足 组合结构构建的已知数据变量，包括：

步骤231、创建取送任务及其取送货件的关联变量c_i，其中，i为取送任务的索引，c_i为取 送任务i对应的货件；

步骤232、创建取送任务及其起点和终点位置的关联变量o_i和d_i，其中，i为取送任务的索 引，o_i为取送任务i的起点位置索引，d_i为取送任务i的终点位置索引；

步骤233、创建取送任务及其直接关联客户所在位置的关联变量l_i，其中，i为取送任务的 索引，l_i为取送任务i的直接关联客户的索引；对于取件类型任务i，有l_i＝o_i；对于取件类型 任务i，有l_i＝d_i；

步骤234、创建取送任务在起点和终点位置上的时间窗变量和其中， i为取送任务的索引，和分别为取送任务i在其起点位置o_i上可执行任务的最早和最晚 时间，和分别为取送任务i在其终点位置d_i上可执行任务的最早和最晚时间；对于不涉 及最早时间或的取送任务，有或对于不涉及最晚时间或的取送任务，有或其中M为一个充分大的正数；

步骤235、创建各类型取送任务的集合，包括送件类型任务集合N^-、取件类型任务集合N⁺， 所有货件最后一项取送任务集合N_E，调移送件任务集合和调移取件任务集合等，其中， 送件类型任务集合N^-为所有送件任务和调移送件任务的集合，取件类型任务集合N⁺为所有 取件任务和调移取件任务的集合；

步骤236、创建取送任务及其运输量的关联变量w_i，其中，i为取送任务的索引，w_i为取 送任务i所取送货件的运输量；

步骤237、创建包括仓库与客户，和客户与客户位置之间的距离变量td_ll′，及行驶时间变 量t_ll′；其中，l和l′分别为行驶路径弧段的起点和终点位置的索引，特别地，若两个位置之间 无直接通路，则有td_ll′＝M和t_ll′＝M，M为一个充分大的正数；；

步骤238、货物取送车辆调度应用场景中会存在货件被送至客户处进行中间加工，随后再 被取出送至另一地点的情况，因此创建取送任务及其对应货件在客户处的中间加工时间的关 联变量T_i，其中，i为取送任务的索引，T_i为取送任务i所对应货件在客户处的中间加工时间；

步骤239、创建执行取送任务的车辆及其最大装载量的关联变量Q，其中，Q为车辆的最 大装载量。

步骤3、对货物取送车辆调度应用场景的系统特征进行分析，根据系统特征创建用于构建 组合结构的约束及其构建组合结构的方法，并根据它们对组合结构的影响形式进行归类，包 括局部组合结构约束、整体组合结构约束和超组合结构约束。

步骤31、创建局部组合结构约束，包括排序结构约束和分组结构约束等在其他局部组合 结构固定的条件下仅通过改变一种局部组合结构便可满足的约束，包括：

步骤311、排序结构是取送任务执行顺序的结构化表达，它可由一个取送任务构成的有序 序列显示表达，排序结构约束是仅通过改变取送任务的排序结构就能被满足的约束，创建排 序结构约束，包括：当取送问题中的货件须被多次取送时，由于只有当货件被送至某个地点 后，才能从该地点取出货件送至其他地点，因此同一货件的取送任务之间客观上存在优先序 关系，满足这种优先序关系只须对取送任务的排序结构加以约束，为此，建立同一货件的取 送优先序约束，并将其归类于排序结构约束；同一货件的取送优先序约束要求在取送任务的 排序结构中，同一货件的取送任务执行顺序满足 其中，符号是同一货件紧邻取送任务优先序关系符号，表示前者在执行顺序上优 先于后者；基于同一货件的取送优先序约束的排序结构构建方法为当一个取送任务序列中存 在某个货件c_i的取送任务不满足优先序约束时，采用顺序互换的方式进行排序结构的局部调 整，使取送任务序列满足同一货件的取送优先序约束；

步骤312、将一个取送任务分组定义为车辆从仓库出发并回到仓库的过程中所执行的取送 任务集合，分组结构是分组数量以及各分组包含的取送任务的结构化表达，创建仅通过改变 取送任务的分组结构就能被满足的分组结构约束，包括：

步骤3121、根据取送系统的实际需要常对客户进行分区管理，使得执行取送任务的车辆 不跨区域执行同一组取送任务，即不同分区客户的取送任务不在同一分组中，在取送问题中 表现为不同分区的客户之间不存在直接通路，为此，建立相异分区取送任务不同组约束，并 将其归类于分组结构约束；相异分区取送任务不同组约束要求在取送任务的分组中，当顺序 相邻的两项取送任务所对应的客户位置之间无直接通路时，将这两项取送任务依次分配至相 邻的两个分组；基于相异分区取送任务不同组约束的分组结构构建方法用于在一个固定的排 序结构下构建一个分组，根据外部输入的最大分组范围[k⁺，k^-]，即当前分组中首、末两项取 送任务在当前排序结构下的执行顺序序号，从首项取送任务开始，依排序结构中顺序序号 的升序检索当前分组中顺序相邻的两项取送任务i_r和i_r+1所对应的直接客户的位置和之 间是否存在直接通路，其中r≥k⁺且r+1≤k^-，如果不存在直接通路，即则将 取送任务i_r作为当前分组的最后一项取送任务，并将当前分组的末项取送任务的执行顺序序 号更新为r，返回当前分组中末项取送任务的执行顺序序号r，其中，i_r为在当前固定排序结 构下第r项取送任务的索引；

步骤3122、执行取送任务的车辆存在最大装载量的限制，当待取送货件足够多时，由于 车辆装载能力限制使得它不能一次完成所有取送任务，因此须将取送任务分配到不同分组中， 并逐一完成每个分组的取送任务，为此，建立车辆能力约束，并将其归类于分组结构约束； 车辆能力约束要求在一个固定排序结构下执行每个分组的取送任务过程中，车辆负载不超过 其最大装载量，包括：

步骤31221、车辆的始发负载量须满足车辆能力约束，为此，建立车辆始发能力约束；车 辆始发能力约束要求在一个固定排序结构下的取送任务分组中，起点为仓库和终点为客户的 送件类型任务所运送货件的总量不大于任务执行车辆的最大装载量；基于车辆始发能力约束 的分组结构构建方法用于在一个固定的排序结构下构建一个分组，根据外部输入的最大分组 范围[k⁺，k^-]，即当前分组中首、末两项取送任务在当前排序结构下的执行顺序序号，从首项 取送任务开始，依排序结构中顺序序号升序检索当前分组中的取送任务i_r，对起点为仓库 和终点为客户位置的送件类型任务所运送货件的运输量进行累加，若出现 且则将取送任务i_k作为当前分组的最后一项 取送任务，将作为车辆始发负载量，并将当前分组的末项取送任务i_k的 执行顺序序号更新为k，其中，k+1≤k^-，i_r为在当前固定排序结构下第r项取送任务的索 引，σ_i为符号变量，当取送任务i为送件类型任务时σ_i取值为-1，为取件类型任务时σ_i取值为1；

步骤31222、由于执行送件类型任务后车辆负载量减少，执行取件类型任务后车辆负载量 增加，而一个分组中可同时出现上述两种类型的任务，因此在执行一个分组的取送任务过程 中车辆负载量动态变化，此时要求在执行当前分组中任何一项取送任务后车辆的负载都满足 车辆能力约束，为此，建立车辆在途能力约束；车辆在途能力约束要求在取送任务的分组中， 车辆执行每项取送任务后的负载量不大于车辆的最大装载量；基于车辆在途能力约束的分组 结构构建方法用于在一个固定的排序结构下构建一个分组，根据最大分组范围[k⁺，k^-]，即当 前分组中首、末两项取送任务在当前排序结构下的执行顺序序号，从首项取送任务开始， 依排序结构中顺序序号的升序检索当前分组中的取送任务i_r，利用公式 计算得到当前取送任务i_k执行完毕后的负载量其中，k≤k^-，i_r为 在当前固定排序结构下第r项取送任务的索引，σ_i为符号变量，当取送任务i为送件类型任务 时σ_i取值为-1，为取件类型任务时σ_i取值为1，根据的取值情况，更新车辆始发负载量q⁰和 当前分组的末项取送任务的执行顺序序号k^-，包括：

步骤312221、如果转步骤312223；

步骤312222、令k＝k+1，如果k≤k^-，则计算转步骤312221；否则，将取送任 务作为当前分组的最后一项取送任务，返回当前分组中末项取送任务的执行顺序序号k^-；

步骤312223、在范围(k，k^-]内按降序检索执行顺序r所对应的取送任务i_r，如果不存在起 点为仓库和终点为客户位置的送件类型任务i_r，则将取送任务i_k-1作为当前分组的最后一项取 送任务，并将当前分组的末项取送任务的执行顺序序号更新为k-1，返回当前分组中末项取 送任务的执行顺序序号k-1；

步骤312224、在范围(k，k^-]内按降序检索执行顺序r所对应的取送任务i_r，将首个出现的 起点为仓库和终点为客户位置的送件类型任务i_r所对应运输量从q⁰中扣除，即有将当前分组的末项取送任务的执行顺序序号k^-更新为r-1；计算转步骤312221；

步骤32、创建须改变多种局部组合结构才能满足的整体组合结构约束；取送问题中常对 取送任务i限定任务开始的最晚时间和任务结束的最晚时间当车辆到达取送任务i的 起点时间ts_i大于时，或当车辆到达取送任务i的终点时间te_i大于时，取送任务i在时间 上不满足相应时间窗的上界要求，为此，建立货件取送时间窗上界约束，并将其归类于整体 组合结构约束；货件取送时间窗上界约束要求每项取送任务的执行时间都不超过货件取送时 间窗上界；基于货件取送时间窗上界约束的组合结构构建方法用于对分组中不满足货件取送 时间窗上界约束的局部组合结构进行局部调整，根据外部输入的最大分组范围[k⁺，k^-]，即当 前分组中首、末两项取送任务在当前排序结构下的执行顺序序号，从首项取送任务开始， 依排序结构中顺序序号的升序检索当前分组中的取送任务i_r，对取送任务i_r的处理方法包括：

步骤321、如果或转步骤333；

步骤322、令r＝r+1，如果r≤k^-，转步骤331；否则，将取送任务作为当前分组的最后一项取送任务，返回局部调整成功指示器LT＝true和分组序号指示器g＝-1，此处g取值为-1表示不须对已构建的排序结构和分组结构进行局部调整；

步骤323、在范围[0，r)内随机选择一个执行顺序r′，要求取送任务i_r′和i_r所对应的货件不 同当时有当时有如果不存在这 样的取送任务i_r′，则不再基于当前已构建的局部组合结构继续进行后续操作，返回局部调整 成功指示器LT＝false和分组序号指示器g＝-1；

步骤324、交换取送任务i_r′和i_r在当前排序结构中的顺序；

步骤325、若当前排序结构中货件和的取送任务不满足同一货件的取送优先序约束， 则采用基于同一货件的取送优先序约束的排序结构构建方法进行取送任务执行顺序的局部调 整；

步骤326、返回局部调整成功指示器LT＝true和分组序号指示器g＝g_r′，其中，g_r′为执 行顺序为r′的取送任务当前所在分组的分组序号。

步骤33、创建保证在固定组合结构下取送方案唯一性的超组合结构约束，它们是在组合 结构的基础上对取送任务执行方式上的约定，包括货件取送时间窗下界约束，货件中间加工 完成时间约束和离开客户时间约束等，包括：

步骤331、货物取送车辆调度应用场景中常对取送任务i限定任务开始的最早时间和任 务结束的最早时间并结合客观实际允许车辆提前到达相应地点进行等待，为此，建立货 件取送时间窗下界约束，并将其归类于超组合结构约束；货件取送时间窗下界约束要求当车 辆到达取送任务i的起点时间时车辆须等待至时间之后开始执行取送任务i，当车 辆到达取送任务i的终点时间时，车辆须等待至时间之后开始执行取送任务i；

步骤332、取送问题中存在货件被送至客户处进行诸如包装之类中间加工作业，中间加工 完毕的货件由车辆择机取出并送至另一地点的情况，结合客观实际允许车辆提前到达相应地 点进行等待，为此，建立货件中间加工完成时间约束，并将其归类于超组合结构约束；货件 中间加工完成时间约束要求当送件类型任务i将货件c_i送至客户处l_i，经时长为T_i的中间加工 后，由取件类型任务j取出送至另一地点，若车辆到达取件类型任务j起点的时间 车辆须等待至时间之后开始执行任务j；

步骤333、为了使得取送任务的执行尽可能地高效，要求车辆尽可能早地离开当前客户以 执行下一项取送任务，为此，建立离开客户时间约束，并将其归类于超组合结构约束；离开 客户时间约束要求，对于送件类型任务i，当货件c_i被送入客户处l_i后，车辆从客户处l_i离开的 时间为对于取件类型任务i，当货件c_i被送入客户处l_i后，车辆从客户处l_i离开的时间为其中，ta_i和tl_i分别为车辆到达和离开l_i的时间，为 车辆可将货件送入客户处l_i的最早时间，th_i为货件c_i的中间加工完毕时间。

步骤4、以取送任务作为取送方案组合结构的基本构造单位，由于超组合结构约束的存在， 使得每个取送方案唯一对应一个组合结构，因此当组合结构确定时取送方案也确定，为此采 用构造可行组合结构的方式替代可行方案的构造；可行组合结构的构造方法为基于组合结构 的约束条件及其分类，随机地构造一个排序结构，利用基于排序结构约束的排序结构构建方 法对当前排序结构进行局部调整使其满足约束，在排序结构基础上利用分组结构约束逐一构 建每个分组，在每个分组创建完成后利用整体结构约束进行检测和局部的组合结构调整，若 组合结构构造过程中存在不满足整体结构约束，并且不能基于整体结构约束的组合结构构建 方法进行结构的局部调整，则重复执行可行组合结构的构造方法，直至得到一个完整的、可 行的组合结构。具体包括：

步骤41、随机生成一个以取送任务为基本单元的有序序列seq；

步骤42、使用基于同一货件的取送优先序约束的排序结构构建方法对序列seq中取送作业 进行执行顺序的局部调整，调整后的序列seq满足同一货件的取送优先序约束；

步骤43、依次构造每个分组并基于整体组合结构约束进行局部调整，包括：

步骤431、初始化当前分组序号g＝0，初始化当前分组首项和末项取送任务在序列seq中 的执行顺序序号和其中，n为取送任务总数；

步骤432、根据和确定的最大分组范围，使用基于相异分区取送任务不同组约束的 分组结构构建方法更新序列seq中分组g的末项取送任务在当前序列seq中的顺序序号

步骤433、根据和确定的最大分组范围，使用基于车辆始发能力约束的分组结构构 建方法更新序列seq中分组g的末项取送任务在当前序列seq中的顺序序号使用基于车辆 在途能力约束的分组结构构建方法更新序列seq中分组g的末项取送任务在当前序列seq中的 顺序序号

步骤434、根据和确定的最大分组范围，使用基于货件取送时间窗上界约束的组合 结构构建方法对分组中不满足货件取送时间窗上界约束的局部组合结构进行局部调整，根据 返回结果进行相应操作，包括：

步骤4341、如果局部调整成功指示器LT_g＝false，则转步骤41；

步骤4342、如果基于货件取送时间窗上界约束的组合结构构建方法返回的分组序号指示 器g′的取值为-1，则令g＝g+1、和转步骤432；

步骤4343、对于基于货件取送时间窗上界约束的组合结构构建方法返回的分组序号指示 器的取值g′，如果g′≤1，则令g＝0、和否则，令g＝g′-1、和转步骤432。

步骤5、将步骤4置于诸如蚁群算法的元启发式方法，对方案进行迭代进化，最终根据终 止规则停止迭代，输出优化方案。

所述的一种用于货物取送的单车调度方法，所述步骤4还包括：

如附图1所示，本发明方法按照以下步骤进行：

步骤1、获取货物取送车辆调度应用场景中系统特征数据，包括地理信息、配送车辆信息、 货件信息和客户仓容量信息等，地理信息数据包括仓库位置、客户位置、客户位置分区、各 位置之间的连接路径、各路径弧段的平均行驶速度信息等，配送车辆信息数据包括可执行任 务车辆数、车辆的可执行任务时间范围信息、车型信息等，货件信息数据包括货件运输量、 货件对应的取送起点和终点位置、货件在其起点和终点位置上的时间窗，在客户处进行中间 加工的时间等；

步骤2、对货物取送车辆调度应用场景中系统特征数据进行预处理，创建取送任务，以取 送任务为基本单元，将其与系统特征数据进行关联，进而将取送问题的系统特征数据处理为 符合组合结构构建要求的形式；

根据仓库位置、客户位置、客户位置分区和各位置之间的连接路径信息，计算得到各位 置之间的路径长度td_od，再结合各路径弧段的平均行驶速度信息计算得到配送车辆在各路径 上的行驶时间t_od，实施例中数据如表1所示；整理得到配送车辆信息，实施例为单车问题， 车辆最大负载量Q假设为15个单位当量，可执行任务时间范围为[0，∞)；进行数据预处理， 创建每个取送任务i，取送任务及其取送货件的关联变量c_i，取送任务及其起点和终点位置的 关联变量o_i和d_i，取送任务及其直接关联客户所在位置的关联变量l_i，取送任务在起点和终点 位置上的时间窗变量和取送任务及其运输量的关联变量w_i，取送任务及 其对应货件在客户处的中间加工时间的关联变量T_i，实施例中数据如表2所示；各类型取送 任务的集合，包括送件类型任务集合N^-、取件类型任务集合N⁺，所有货件最后一项取送任 务集合N_E，调移送件任务集合和调移取件任务集合实施例中数据如表3所示；

表1实施例中各位置之间的走行时间(min)

表2实施例中取送任务及其关联变量数据信息

表3实施例中各类型取送任务集合

步骤3、对货物取送车辆调度应用场景中系统特征进行分析分析，根据系统特征创建用于 构建组合结构的约束及其构建组合结构的方法，并根据它们对组合结构的影响形式进行归类， 包括局部组合结构约束、整体组合结构约束和超组合结构约束；

建立同一货件的取送优先序约束，并将其归类于排序结构约束，它要求同一货件的取送 任务执行顺序满足为此，构建基于同一 货件的取送优先序约束的排序结构构建方法，如图2所示；图2中的符号含义为，△、□和 ○中的数字分别表示c_i取值为0、1和2的取送任务序号，虚线下方数字为执行顺序序号，为分组标识符，两个之间的取送任务为同一分组的任务，为空分组标识符，表示可在此处 插入分组标识符；图2中，seq₀为随机生成的一个取送任务有序序列，其中，取送任务2和 3不满足同一货件的取送优先序约束，因此采用基于同一货件的取送优先序约束的排序结构 构建方法进行取送任务的执行顺序交换，得到取送任务的有序序列seq。

建立相异分区取送任务不同组约束，并将其归类于分组结构约束，它要求在取送任务的 分组中，当顺序相邻的两项取送任务所对应的客户位置之间无直接通路时，将这两项取送任 务依次分配至相邻的两个分组，为此，构建基于相异分区取送任务不同组约束的分组结构构 建方法，如图3所示，在seq₀中根据外部输入的最大分组范围[k⁺，k^-]，若k⁺＝0，k^-＝7， 从首项取送任务i₀＝2开始，依排序结构中顺序序号的升序检索当前分组中顺序相邻的两项取 送任务i_r和i_r+1所对应直接客户的位置和之间是否存在直接通路，其中r≥0且 r+1≤7，在检索过程中，取送任务i₂＝5和i₃＝0所对应的直接客户的位置l₅＝2和l₀＝3之 间不存在直接通路，即有td_2，3＝t_2，3＝M，则将取送任务5作为当前分组的最后一项取送任务， 并将当前分组的末项取送任务的执行顺序序号k^-更新为2，返回当前分组中末项取送任务的 执行顺序序号2，此时，得到如序列seq所示的一个新分组标识符。

建立车辆能力约束，并将其归类于分组结构约束，它要求要求在一个固定排序结构下执 行每个分组的取送任务过程中，车辆负载不超过其最大装载量，进一步地，将车辆能力约束 分解为车辆始发能力约束和车辆在途能力约束；车辆始发能力约束要求在一个固定排序结构 下的取送任务分组中，起点为仓库和终点为客户的送件类型任务所运送货件的总量不大于任 务执行车辆的最大装载量，为此，构建基于车辆始发能力约束的分组结构构建方法，如图4 所示，在seq₀中根据外部输入的最大分组范围[k⁺，k^-]，若k⁺＝3，k^-＝7，从首项取送任务 i₃＝0开始，依排序结构中顺序序号的升序检索当前分组中的取送任务i_r，对起点为仓库和终 点为客户位置的送件类型任务所运送货件的运输量进行累加，由于 (1+σ₀)w₀＝18＜2Q＝30且(1+σ₀)w₀+0+0+(1+σ₆)w₆＝40＞2Q＝30，则将取送任 务i₅＝1作为当前分组的最后一项取送任务，将作为车辆始发 负载量，并将当前分组的末项取送任务的执行顺序序号k^-更新为5，此时，得到如序列seq所 示的一个新分组标识符；车辆在途能力约束要求在取送任务的分组中，车辆执行每项取送任 务后的负载量不大于车辆的最大装载量，为此，构建基于车辆在途能力约束的分组结构构建 方法，如图5所示，在seq₀中根据最大分组范围[k⁺，k^-]，其中，k⁺＝3，k^-＝5从首项取送 任务i₃开始，依排序结构中顺序序号的升序检索当前分组中的取送任务i_r，利用公式计算得到当前取送任务i_k执行完毕后的负载量其中，k≤k^-，并按如下步骤执行，首先令k＝k⁺计算q₀＝q⁰+σ₀w₀＝0，由于因此令 k＝k+1，又因为k＝4≤k^-＝5，则计算因为因此令k＝k+1， 又因为k＝5≤k^-＝5，则计算因为所以在(k，k^-]内按降序检索执 行顺序r所对应的取送任务i_r，由于此时k＝k^-＝5，即区间(k，k^-]不存在，故而不存在起点 为仓库和终点为客户位置的送件类型任务i_r，则将取送任务i_k-1＝i₄＝3作为当前分组的最后 一项取送任务，并将当前分组的末项取送任务的执行顺序序号k^-更新为4，此时，得到如序 列seq所示的一个新分组标识符。

建立货件取送时间窗上界约束，并将其归类于整体组合结构约束，它要求每项取送任务 的执行时间都不超过货件取送时间窗上界，为此，构建基于货件取送时间窗上界约束的组合 结构构建方法，如图6所示，在seq₀中根据外部输入的最大分组范围[k⁺，k^-]，若k⁺＝0， k^-＝2，从首项取送任务i₀＝2开始，依排序结构中顺序序号的升序检索当前分组中的取送任 务i_r，基于超组合结构约束中货件取送时间窗下界约束、货件中间加工完成时间约束和离开 客户时间约束对车辆在取送任务执行方式上的约定，假设车辆在当前分组中待命起点时间为8:30，货件中间加工完成时间、车辆到达和离开各位置的时间如表4所示，由于 所以在范围[0，2)内随机选择一个执行顺序r′，要求取送任务i_r′和i₂＝3所对 应的货件不同且有此时仅有r′＝1对应的取送任务i₁＝4满足随机 选择条件，因此，交换取送任务i₁＝4和i₂＝3的执行顺序，得到，返回局部调整成功指示器 LT＝true和须重新确定分组结构的起始分组序号g＝g_r′＝0，此时，得到如序列seq所示的 一个新序列，以及更新后的分组标识符。

表4图6中seq₀的部分车辆时间信息

建立货件取送时间窗下界约束，并将其归类于超组合结构约束，它要求当车辆到达取送 任务i的起点时间时车辆须等待至时间之后开始执行取送任务i，当车辆到达取送 任务i的终点时间时，车辆须等待至时间之后开始执行取送任务i；

建立货件中间加工完成时间约束，并将其归类于超组合结构约束，它要求当送件类型任 务i将货件c_i送至客户处l_i，经时长为T_i的中间加工后，由取件类型任务j取出送至另一地点，若车辆到达取件类型任务j起点的时间车辆须等待至时间th_j之后开始执行任务j；

建立离开客户时间约束，并将其归类于超组合结构约束，它要求对于送件类型任务i，当 货件c_i被送入客户处l_i后，车辆从客户处l_i离开的时间为对于取件类型任务i，当货件c_i被送入客户处l_i后，车辆从客户处l_i离开的时间为其中， ta_i和tl_i分别为车辆到达和离开l_i的时间，为车辆可将货件送入客户处l_i的最早时间，th_i为 货件c_i的中间加工完毕时间；

步骤4、以取送任务作为取送方案组合结构的基本构造单位，由于超组合结构约束的存在， 使得每个取送方案唯一对应一个组合结构，因此当组合结构确定时取送方案也确定，为此采 用构造可行组合结构的方式替代可行方案的构造；可行组合结构的构造方法为基于组合结构 的约束条件及其分类，随机地构造一个排序结构，利用基于排序结构约束的排序结构构建方 法对当前排序结构进行局部调整使其满足约束，在排序结构基础上利用分组结构约束逐一构 建每个分组，在每个分组创建完成后利用整体结构约束进行检测和局部的组合结构调整，若 组合结构构造过程中存在不满足整体结构约束，并且不能基于整体结构约束的组合结构构建 方法进行结构的局部调整，则重复执行可行组合结构的构造方法，直至得到一个完整的、可 行的组合结构，如图7所示；

随机生成一个以取送任务为基本单元的有序序列seq，如图7.1所示；

使用基于同一货件的取送优先序约束的排序结构构建方法对序列seq中取送作业进行执 行顺序的局部调整，调整后的序列seq满足同一货件的取送优先序约束，如图7.2所示；

基于已构建的排序结构，从首个分组0开始构造每个分组；构造每个分组的过程为依次 施以基于相异分区取送任务不同组约束的分组结构构建方法、基于车辆始发能力约束的分组 结构构建方法、基于车辆在途能力约束的分组结构构建方法和基于货件取送时间窗上界约束 的组合结构构建方法；假设车辆初始待命时间为8:30，初始化当前分组序号g＝0，初始化当 前分组首项和末项取送任务在序列seq中的执行顺序序号和构建序列seq的分 组0，序列seq新增一个关于分组0的分组标识符如图7.3所示；构建序列seq的分组1， 序列seq新增一个关于分组1的分组标识符如图7.4所示；构建序列seq的分组2，在施以 基于货件取送时间窗上界约束的组合结构构建方法时，将取送任务4和取送任务3在当前任 务序列中的执行顺序互换，得到一个新取送任务序列seq和分组标识符，如图7.5所示，根据 基于货件取送时间窗上界约束的组合结构构建方法的返回值，将当前分组序号g的取值更新为 0，相应地更新当前分组首项和末项取送任务在序列seq中的执行顺序序号和对当前取送任务序列seq进行分组结构构建，最终得到一个可行的取送任务方案，如图7.6所 示，其中各时间信息如表5所示。

表5图7.6中最终方案的时间信息

步骤5、将步骤4置于诸如蚁群算法的元启发式方法，对方案进行迭代进化，最终根据终 止规则停止迭代，输出优化方案。

初始化元启发式方法中相关参数，如蚁群算法中，以取送任务为节点设置信息素矩阵并 将各路径弧段的信息素初始化为相同的正实数；施以步骤4，其中，随机生成取送任务的有 序序列可根据采用的具体元启发式方法执行，如蚁群算法中每只人工蚂蚁按照状态转移公式 选择下一访问客户节点j，直至得到一个完整的取 送任务有序序列，式中τ_ij为(i，j)弧段之间的信息素测度，η_k(j)是节点j所提供的启发式信息(可 视化信息)，k(j)为节点j对应的车组序号，α为信息素强度因子，β为启发式信息因子；根据 元启发式方法的参数更新机制对其进行参数更新，如蚁群算法中利用 τ_ij(t+1)＝(1-ρ_t)τ_ij(t)+Δτ_ij(t)对信息素矩阵进行更新，式中ρ_t∈(0，1]为当前迭代周期t 的信息素衰减系数，Δτ_ij(t)为当前迭代周期t中路径弧段(i，j)上应累积的信息量；元启发式方 法保持方案的迭代进化，直至满足终止条件，输出货物取送的单车调度计划。

Claims (14)

1.一种用于货物取送的单车调度方法，其特征在于，包括以下步骤：

步骤1：获取货物取送车辆调度应用场景的系统特征数据，包括地理信息数据、配送车辆信息数据、货件信息数据；

步骤2：对货物取送车辆调度应用场景的系统特征数据进行预处理，创建取送任务，以取送任务为基本单元，将其与系统特征数据进行关联，进而将取送问题的系统特征数据处理为符合组合结构构建要求的形式；

步骤3：对货物取送车辆调度应用场景的系统特征进行分析，根据系统特征创建用于构建组合结构的约束及其构建组合结构的方法，并根据其对组合结构的影响形式进行归类，包括局部组合结构约束、整体组合结构约束和超组合结构约束；

步骤4：以取送任务作为取送方案组合结构的基本构造单位，采用构造可行组合结构的方式；可行组合结构的构造方法为基于组合结构的约束条件及其分类，随机地构造一个排序结构，利用基于排序结构约束的排序结构构建方法对当前排序结构进行局部调整使其满足约束，在排序结构基础上利用分组结构约束逐一构建每个分组，在每个分组创建完成后利用整体结构约束进行检测和局部的组合结构调整，若组合结构构造过程中存在不满足整体结构约束，并且不能基于整体结构约束的组合结构构建方法进行结构的局部调整，则重复执行可行组合结构的构造方法，直至得到一个完整的、可行的组合结构；

步骤5：将步骤4置于元启发式方法，对方案进行迭代进化，最终根据终止规则停止迭代，输出优化方案。

2.根据权利要求1所述的用于货物取送的单车调度方法，其特征在于，所述地理信息数据包括仓库位置、客户位置、客户位置分区、各位置之间的连接路径、各路径弧段的平均行驶速度信息；配送车辆信息数据包括可执行任务车辆数、车辆的可执行任务时间范围信息、车型信息；货件信息数据包括货件的运输量、货件对应的取送起点和终点位置、货件在其起点和终点位置上的时间窗，以及在客户处进行中间加工的时间。

3.根据权利要求1所述的用于货物取送的单车调度方法，其特征在于，所述步骤2具体包括：

步骤21：对货物取送车辆调度应用场景中所涉及的每个位置赋予索引，将仓库的索引赋值为0，客户的索引从1开始赋值；对取送问题中包含的每个货件赋予索引；对执行取送任务的每个车辆赋予索引；

步骤22：根据货件所对应的取送起点和终点位置数据为每个货件创建取送任务，包括从仓库送至客户处的送件类型，从客户处取回仓库的取件类型，从客户处取出并送往另一客户处的调移类型；其中调移任务被拆分为调移取件任务和调移送件任务；为每项取送任务赋予唯一的索引；

步骤23：基于货物取送车辆调度应用场景中的特征数据和已创建的取送任务，创建满足组合结构构建的已知数据变量。

4.根据权利要求3所述的用于货物取送的单车调度方法，其特征在于，所述步骤23具体包括：

步骤231：创建取送任务及其取送货件的关联变量c_i；其中，i为取送任务的索引，c_i为取送任务i对应的货件；

步骤232：创建取送任务及其起点和终点位置的关联变量o_i和d_i；其中，o_i为取送任务i的起点位置索引，d_i为取送任务i的终点位置索引；

步骤233：创建取送任务及其直接关联客户所在位置的关联变量l_i；其中，l_i为取送任务i的直接关联客户的索引；对于取件类型任务i，有l_i＝o_i；对于送件类型任务i，有l_i＝d_i；

步骤234：创建取送任务在起点和终点位置上的时间窗变量和其中，和分别为取送任务i在其起点位置o_i上可执行任务的最早和最晚时间，和分别为取送任务i在其终点位置d_i上可执行任务的最早和最晚时间；对于不涉及最早时间或的取送任务，有或对于不涉及最晚时间或的取送任务，有或其中M为一个充分大的正数；

步骤235：创建各类型取送任务的集合，包括送件类型任务集合N^-、取件类型任务集合N⁺，所有货件最后一项取送任务集合N_E，调移送件任务集合和调移取件任务集合其中，送件类型任务集合N^-为所有送件任务和调移送件任务的集合，取件类型任务集合N⁺为所有取件任务和调移取件任务的集合；

步骤236：创建取送任务及其运输量的关联变量w_i；其中，w_i为取送任务i所取送货件的运输量；

步骤237：创建包括仓库与客户，和客户与客户位置之间的距离变量td_ll′，及行驶时间变量t_ll′；其中，l和l′分别为行驶路径弧段的起点和终点位置的索引，特别地，若两个位置之间无直接通路，则有td_ll′＝M和t_ll′＝M，M为一个充分大的正数；

步骤238：创建取送任务及其对应货件在客户处的中间加工时间的关联变量T_i，其中，i为取送任务的索引，T_i为取送任务i所对应货件在客户处的中间加工时间；

步骤239：创建执行取送任务的车辆及其最大装载量的关联变量Q，其中，Q为车辆的最大装载量。

5.根据权利要求1所述的用于货物取送的单车调度方法，其特征在于，所述步骤3具体包括：

步骤31：创建局部组合结构约束，包括排序结构约束和分组结构约束，该约束在其他局部组合结构固定的条件下仅通过改变一种局部组合结构便可满足；

步骤32：创建需改变多种局部组合结构才能满足的整体组合结构约束；建立货件取送时间窗上界约束，并将其归类于整体组合结构约束，货件取送时间窗上界约束要求每项取送任务的执行时间都不超过货件取送时间窗上界；基于货件取送时间窗上界约束的组合结构构建方法具体为：根据外部输入的最大分组范围[k⁺，k^-]，即当前分组中首、末两项取送任务在当前排序结构下的执行顺序序号，从首项取送任务开始，依排序结构中顺序序号的升序检索当前分组中的取送任务i_r；

步骤33：创建保证在固定组合结构下取送方案唯一性的超组合结构约束，包括货件取送时间窗下界约束，货件中间加工完成时间约束和离开客户时间约束。

6.根据权利要求5所述的用于货物取送的单车调度方法，其特征在于，所述步骤31具体包括：

步骤311：创建排序结构约束，方法为：对同一货件的取送任务之间客观上存在优先序关系的任务，对取送任务的排序结构加以约束，建立同一货件的取送优先序约束，并将其归类于排序结构约束；同一货件的取送优先序约束要求在取送任务的排序结构中，同一货件的取送任务执行顺序满足送件任务＞调移取件任务＞调移送件任务＞取件任务；基于同一货件的取送优先序约束的排序结构构建方法具体为：当一个取送任务序列中存在某个货件c_i的取送任务不满足优先序约束时，采用顺序互换的方式进行排序结构的局部调整，使取送任务序列满足同一货件的取送优先序约束；

步骤312：将一个取送任务分组定义为车辆从仓库出发并回到仓库的过程中所执行的取送任务集合，创建仅通过改变取送任务的分组结构就能被满足的分组结构约束。

7.根据权利要求6所述的用于货物取送的单车调度方法，其特征在于，所述步骤312具体包括：

步骤3121：建立相异分区取送任务不同组约束，并将其归类于分组结构约束；相异分区取送任务不同组约束要求在取送任务的分组中，当顺序相邻的两项取送任务所对应的客户位置之间无直接通路时，将这两项取送任务依次分配至相邻的两个分组；基于相异分区取送任务不同组约束的分组结构构建方法为：在一个固定的排序结构下构建一个分组，根据外部输入的最大分组范围[k⁺，k^-]，即当前分组中首、末两项取送任务在当前排序结构下的执行顺序序号，从首项取送任务开始，依排序结构中顺序序号的升序检索当前分组中顺序相邻的两项取送任务i_r和i_r+1所对应的直接客户的位置和之间是否存在直接通路，其中r≥k⁺且r+1≤k^-，如果不存在直接通路，即则将取送任务i_r作为当前分组的最后一项取送任务，并将当前分组的末项取送任务的执行顺序序号更新为r，返回当前分组中末项取送任务的执行顺序序号r，其中，i_r为在当前固定排序结构下第r项取送任务的索引；

步骤3122：建立车辆能力约束，并将其归类于分组结构约束；车辆能力约束要求在一个固定排序结构下执行每个分组的取送任务过程中，车辆负载不超过其最大装载量。

8.根据权利要求7所述的用于货物取送的单车调度方法，其特征在于，所述步骤3122具体包括：

步骤31221：建立车辆始发能力约束；车辆始发能力约束要求在一个固定排序结构下的取送任务分组中，起点为仓库和终点为客户的送件类型任务所运送货件的总量不大于任务执行车辆的最大装载量；基于车辆始发能力约束的分组结构构建方法为：用于在一个固定的排序结构下构建一个分组，根据外部输入的最大分组范围[k⁺，k^-]，即当前分组中首、末两项取送任务在当前排序结构下的执行顺序序号，从首项取送任务开始，依排序结构中顺序序号升序检索当前分组中的取送任务i_r，对起点为仓库和终点为客户位置的送件类型任务所运送货件的运输量进行累加，若出现且则将取送任务i_k作为当前分组的最后一项取送任务，将作为车辆始发负载量，并将当前分组的末项取送任务i_k的执行顺序序号更新为k，其中，k+1≤k^-，i_r为在当前固定排序结构下第r项取送任务的索引，σ_i为符号变量，当取送任务i为送件类型任务时，σ_i取值为-1；为取件类型任务时，σ_i取值为1；

步骤31222、建立车辆在途能力约束；车辆在途能力约束要求在取送任务的分组中，车辆执行每项取送任务后的负载量不大于车辆的最大装载量；基于车辆在途能力约束的分组结构构建方法为：用于在一个固定的排序结构下构建一个分组，根据最大分组范围[k⁺，k^-]，即当前分组中首、末两项取送任务在当前排序结构下的执行顺序序号，从首项取送任务开始，依排序结构中顺序序号的升序检索当前分组中的取送任务i_r，利用公式计算得到当前取送任务i_k执行完毕后的负载量其中，k≤k^-，i_r为在当前固定排序结构下第r项取送任务的索引，σ_i为符号变量，当取送任务i为送件类型任务时，σ_i取值为-1；为取件类型任务时，σ_i取值为1；根据的取值情况，更新车辆始发负载量q⁰和当前分组的末项取送任务的执行顺序序号k^-。

9.根据权利要求8所述的用于货物取送的单车调度方法，其特征在于，所述步骤31222还包括：

步骤312221：如果q_ik>Q，转步骤312223，否则进入下一步；

步骤312222：令k＝k+1，如果k≤k^-，则计算q_ik，转步骤312221；否则，将取送任务作为当前分组的最后一项取送任务，返回当前分组中末项取送任务的执行顺序序号k^-；

步骤312223：在范围(k，k^-]内按降序检索执行顺序r所对应的取送任务i_r，如果不存在起点为仓库和终点为客户位置的送件类型任务i_r，则将取送任务i_k-1作为当前分组的最后一项取送任务，并将当前分组的末项取送任务的执行顺序序号更新为k-1，返回当前分组中末项取送任务的执行顺序序号k-1；

步骤312224：在范围(k，k^-]内按降序检索执行顺序r所对应的取送任务i_r，将首个出现的起点为仓库和终点为客户位置的送件类型任务i_r所对应运输量从q⁰中扣除，即有将当前分组的末项取送任务的执行顺序序号k^-更新为r-1；计算q_ik，转步骤312221。

10.根据权利要求6所述的用于货物取送的单车调度方法，其特征在于，所述步骤32中对取送任务i_r的处理方法包括：

步骤321：如果或转步骤323；

步骤322：令r＝r+1，如果r≤k^-，转步骤321；否则，将取送任务作为当前分组的最后一项取送任务，返回局部调整成功指示器LT＝true和分组序号指示器g＝-1，此处g取值为-1表示不需对已构建的排序结构和分组结构进行局部调整；

步骤323：在范围[0，r)内随机选择一个执行顺序r′，要求取送任务i_r′和i_r所对应的货件不同当时有当时有如果不存在这样的取送任务i_r′，则不再基于当前已构建的局部组合结构继续进行后续操作，返回局部调整成功指示器LT＝false和分组序号指示器g＝-1；

步骤324：交换取送任务i_r′和i_r在当前排序结构中的顺序；

步骤325：若当前排序结构中货件和的取送任务不满足同一货件的取送优先序约束，则采用基于同一货件的取送优先序约束的排序结构构建方法进行取送任务执行顺序的局部调整；

步骤326：返回局部调整成功指示器LT＝true和分组序号指示器g＝g_r′，其中，g_r′为执行顺序为r′的取送任务当前所在分组的分组序号。

11.根据权利要求5所述的用于货物取送的单车调度方法，其特征在于，所述步骤33具体包括：

步骤331：建立货件取送时间窗下界约束，并将其归类于超组合结构约束；货件取送时间窗下界约束要求当车辆到达取送任务i的起点时间时车辆需等待至时间之后开始执行取送任务i，当车辆到达取送任务i的终点时间时，车辆需等待至时间之后开始执行取送任务i；

步骤332：建立货件中间加工完成时间约束，并将其归类于超组合结构约束；货件中间加工完成时间约束要求当送件类型任务i将货件c_i送至客户处l_i，经时长为T_i的中间加工后，由取件类型任务j取出送至另一地点，若车辆到达取件类型任务j起点的时间车辆需等待至时间之后开始执行任务j；

步骤333：建立离开客户时间约束，并将其归类于超组合结构约束；离开客户时间约束要求，对于送件类型任务i，当货件c_i被送入客户处l_i后，车辆从客户处l_i离开的时间为对于取件类型任务i，当货件c_i被送入客户处l_i后，车辆从客户处l_i离开的时间为其中，ta_i和tl_i分别为车辆到达和离开l_i的时间，为车辆可将货件送入客户处l_i的最早时间，th_i为货件c_i的中间加工完毕时间。

12.根据权利要求8所述的用于货物取送的单车调度方法，其特征在于，所述步骤4具体包括：

步骤41：随机生成一个以取送任务为基本单元的有序序列seq；

步骤42：使用基于同一货件的取送优先序约束的排序结构构建方法对序列seq中取送作业进行执行顺序的局部调整，调整后的序列seq满足同一货件的取送优先序约束；

步骤43：依次构造每个分组并基于整体组合结构约束进行局部调整。

13.根据权利要求12所述的用于货物取送的单车调度方法，其特征在于，所述步骤43具体包括：

步骤431：初始化当前分组序号g＝0，初始化当前分组首项和末项取送任务在序列seq中的执行顺序序号和其中，n为取送任务总数；

步骤432：根据和确定的最大分组范围，使用基于相异分区取送任务不同组约束的分组结构构建方法更新序列seq中分组g的末项取送任务在当前序列seq中的顺序序号

步骤433：根据和确定的最大分组范围，使用基于车辆始发能力约束的分组结构构建方法更新序列seq中分组g的末项取送任务在当前序列seq中的顺序序号使用基于车辆在途能力约束的分组结构构建方法更新序列seq中分组g的末项取送任务在当前序列seq中的顺序序号

步骤434：根据和确定的最大分组范围，使用基于货件取送时间窗上界约束的组合结构构建方法对分组中不满足货件取送时间窗上界约束的局部组合结构进行局部调整，根据返回结果进行相应操作。

14.根据权利要求13所述的用于货物取送的单车调度方法，其特征在于，所述步骤434具体包括：

步骤4341：如果局部调整成功指示器LT_g＝false，则转步骤41；

步骤4342：如果基于货件取送时间窗上界约束的组合结构构建方法返回的分组序号指示器g′的取值为-1，则令g＝g+1、和转步骤432；

步骤4343：对于基于货件取送时间窗上界约束的组合结构构建方法返回的分组序号指示器的取值g′，如果g′≤1，则令g＝0、和否则，令g＝g′-1、和转步骤432。