CN114611806A - 一种解决多配送中心两级城市协同配送的大邻域搜索方法 - Google Patents

一种解决多配送中心两级城市协同配送的大邻域搜索方法 Download PDF

Info

Publication number
CN114611806A
CN114611806A CN202210260534.2A CN202210260534A CN114611806A CN 114611806 A CN114611806 A CN 114611806A CN 202210260534 A CN202210260534 A CN 202210260534A CN 114611806 A CN114611806 A CN 114611806A
Authority
CN
China
Prior art keywords
network
distribution
latest
operator
solution
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202210260534.2A
Other languages
English (en)
Other versions
CN114611806B (zh
Inventor
邓蕾
贾圣阳
赵泉午
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing University
Original Assignee
Chongqing University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Chongqing University filed Critical Chongqing University
Priority to CN202210260534.2A priority Critical patent/CN114611806B/zh
Publication of CN114611806A publication Critical patent/CN114611806A/zh
Application granted granted Critical
Publication of CN114611806B publication Critical patent/CN114611806B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G06Q10/047—Optimisation of routes or paths, e.g. travelling salesman problem
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06311—Scheduling, planning or task assignment for a person or group
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/08—Logistics, e.g. warehousing, loading or distribution; Inventory or stock management
    • G06Q10/083—Shipping

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Marketing (AREA)
  • Development Economics (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Game Theory and Decision Science (AREA)
  • Educational Administration (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种解决多配送中心两级城市协同配送的大邻域搜索方法，包括以下步骤：S1：根据客户点与网点之间的距离，构造初始可行解；S2：对初始可行解进行大邻域搜索，确定最新第二级配送网络；S3：对最新第二级配送网络进行优化；S4：基于优化后的最新第二级配送网络，对市趟运输线路进行重构，得到可行市趟运输线路；S5：对可行市趟运输线路进行优化，并对由配送中心与网点之间的车辆运输线路构成的一级路径进行优化，得到最新可行解；S6：根据最新可行解确定当前解，并根据模拟退火准则判断是否达到算法终止条件，若是则结束大邻域搜索，否则返回步骤S2。本发明通过整体优化两级配送网络中的市趟运输线路、最后一公里线路以及车辆调度，降低了整体网络成本而不仅是某单一配送中心配送网络成本。

Description

一种解决多配送中心两级城市协同配送的大邻域搜索方法
技术领域
本发明属于配送网络技术领域，具体涉及一种解决多配送中心两级城市协同配送的大邻域搜索方法。
背景技术
在多配送中心的两级城市物流配送网络中(两级网络是指由配送中心到网点、网点到客户点组成的网络)，协调各配送中心的市趟运输(指将货物从配送中心发往城市内各网点)线路与车辆调度可以降低网络运营成本、提高网点运作效率、提升客户时效感知度。
传统规划市趟运输线路和车辆调度的方法是使用人力借助Excel完成市趟运输线路的排布与车辆调度。这种传统的方法在大规模规划任务中难以获得高质量的市趟运输线路和车辆调度方案，也不能根据实际运作情况快速、灵活地调整方案以使得网络高效地运作。目前，某些物流配送公司使用人工智能、运筹优化算法等规划市趟运输线路和车辆调度，这些算法将多配送中心网络中划分成多个单配送中心网络进行规划，得到单配送中心网络的较优市趟运输和车辆调度方案，依然未能解决多配送中心网络的协调运输问题，从而使得整体网络效率偏低，成本偏高。
发明内容
本发明的目的是为了解决满足多配送中心网络协调与同步性要求下的市趟运输线路和车辆调度规划问题，提出了一种解决多配送中心两级城市协同配送的大邻域搜索方法。
本发明的技术方案是：一种解决多配送中心两级城市协同配送的大邻域搜索方法包括以下步骤：
S1：根据客户点与网点之间的距离，构造初始可行解；
S2：对初始可行解进行大邻域搜索，确定最新第二级配送网络；
S3：对最新第二级配送网络进行优化；
S4：基于优化后的最新第二级配送网络，对市趟运输线路进行重构，得到可行市趟运输线路；
S5：对可行市趟运输线路进行优化，并对由配送中心与网点之间的车辆运输线路构成的一级路径进行优化，得到最新可行解；
S6：根据最新可行解确定当前解，并根据模拟退火准则判断是否达到算法终止条件，若是则结束大邻域搜索，否则返回步骤S2。
进一步地，步骤S1包括以下子步骤：
S11：根据客户点与网点之间的距离，利用轮盘赌选择法将客户点分配至各网点，确定各网点对各配送中心的货物需求量；
S12：根据各网点对各配送中心的货物需求量，利用扫描算法确定各配送中心与各网点的市趟运输线路和车辆调度；
S13：根据市趟运输线路和车辆调度中车辆到达网点的时间，利用扫描算法确定各网点到达各客户点的最后一公里运输线路，确定初始可行解。
进一步地，步骤S2中，确定最新二级配送网络的具体方法为：利用移除算子和插入算子对网点与客户点组成的二级配送网络进行大邻域搜索，确定最新第二级网络；其中，移除算子包括随机移除算子、最坏移除算子、相关移除算子和线路移除算子；插入算子包括贪婪插入算子、贪婪扰动插入算子、贪婪禁忌插入算子和遗憾插入算子。
进一步地，随机移除算子进行大邻域搜索的具体方法为：从当前解中随机移除若干个客户点，得到当前解的邻域；
最坏移除算子进行大邻域搜索的具体方法为：从当前解中移除若干个使配送成本降低最多的客户点，得到当前解的邻域；
相关移除算子进行大邻域搜索的具体方法为：迭代运行相关移除算子若干次，随机选择一个客户点作为种子客户，移除与该种子客户相关度量最大的客户点，并从当前解中移除若干个客户点，得到当前解的邻域；
线路移除算子进行大邻域搜索的具体方法为：迭代随机移除各网点到达各客户点的最后一公里运输线路以及最后一公里运输线路中的客户点，直至被移除的客户点达到设定阈值，得到当前解的邻域；
贪婪插入算子进行大邻域搜索的具体方法为：将被移除的若干个客户点重新插入二级配送网络中；
贪婪扰动插入算子进行大邻域搜索的具体方法为：将被移除的若干个客户点重新插入二级配送网络中，并对插入产生的配送成本进行扰动；
贪婪禁忌插入算子进行大邻域搜索的具体方法为：将被移除的若干个客户点重新插入二级配送网络中，且插入位置不与上一次迭代时客户点的运输线路位置重叠；
遗憾插入算子进行大邻域搜索的具体方法为：将遗憾值最大的客户点重新迭代插入二级配送网络中；
其中，将初始可行解作为大邻域搜索的初始当前解，并将大邻域搜索中迭代产生的最新解作为当前解。
进一步地，相关度量R(i,j)的计算公式为：
Figure BDA0003549973410000031
其中，di,j表示客户点i和客户点j之间的距离，
Figure BDA0003549973410000032
表示客户点i对一个配送中心的货物需求量，
Figure BDA0003549973410000033
表示客户点i对另一个配送中心的货物需求量，
Figure BDA0003549973410000034
表示客户j对一个配送中心的货物需求量，
Figure BDA0003549973410000035
表示客户j对另一个配送中心的货物需求量，
Figure BDA0003549973410000036
表示第一权重参数，λ表示第二权重参数，θ表示第三权重参数；
遗憾值最大的客户点i:的计算公式：
Figure BDA0003549973410000037
其中，D表示客户集合，k表示插入方式总数，fi h表示客户i的第h低的插入成本，fi 1表示客户i的最低成本。
进一步地，步骤S3中，对最新二级配送网络进行优化的具体方法为：以市趟运输线路的片段为单位，利用2-opt*算子在最新二级配送网络中进行组合搜索，利用Swap算子交换最新二级配送网络中的所有客户点，得到最优运输线路，完成最新二级配送网络优化。
进一步地，步骤S4中，对市趟运输线路进行重构的具体方法为：将距离配送中心最远的网点插入至空闲市趟运输线路中，迭代将所有未规划的网点插入至空闲市趟运输线路中运输成本最少的线路位置处，并根据运输车辆的容量限制、网点的作业时间窗限制和客户点时效限制，将最新的市趟运输线路分割得到若干条可行市趟运输线路。
进一步地，步骤S5中，对可行市趟运输线路进行优化的具体方法为：利用网点移除算子将可行市趟运输线路中各运输车辆不同步到达的网点移除，利用网点插入算子将被移除的网点重新插入至可行市趟运输线路中，以使优化后的最新二级配送网络的运输成本最小，完成对可行市趟运输线路的优化；
步骤S5中，利用2-opt局部搜索算子、relocate局部搜索算子和swap局部搜索算子对一级路径进行优化。
进一步地，步骤S6中，若最新可行解目标函数值低于设定的当前最优目标函数值，则将最新可行解作为当前最优解；
若最新可行解目标函数值低于当前解目标函数值但高于当前最优目标函数值，则将最新可行解作为当前解；
若最新可行解目标函数值劣于当前解目标函数值，则根据最新可行解的接受概率
Figure BDA0003549973410000038
判断是否将最新可行解目标函数值作为当前解，其中，Tcur表示当前温度，f(snew)和f(scur)分别表示新解snew和当前解scur的目标函数值；
步骤S6中，若模拟退火准则中设定的当前温度Tcur低于结束温度Tend，则输出当前最优解，否则返回步骤S2进入下一次迭代，直至结束大邻域搜索。
本发明的有益效果是：
(1)本发明通过整体优化两级配送网络中的市趟运输线路、最后一公里线路以及车辆调度，降低了整体网络成本而不仅是某单一配送中心配送网络成本；
(2)本发明使用精心设计的网点移除和插入算子调度市趟运输车辆到达各网点的时间，消除了市趟运输车辆到达网点的时间超出网点作业班次时间窗所带来的网点运作痛点；
(3)本发明采用更优的线路排布和车辆调度使得车辆使用数量大幅减少，降低了企业的资金占用。
附图说明
图1为大邻域搜索方法的流程图。
具体实施方式
下面结合附图对本发明的实施例作进一步的说明。
在描述本发明的具体实施例之前，为使本发明的方案更加清楚完整，首先对本发明中出现的缩略语和关键术语定义进行说明：
轮盘赌选择法：是最简单也是最常用的选择方法，在该方法中，各个个体的选择概率和其适应度值成比例，适应度越大，选中概率也越大。但实际在进行轮盘赌选择时个体的选择往往不是依据个体的选择概率，而是根据“累积概率”来进行选择。
扫描算法：采用极坐标来表示各需求点的区位，然后任取一需求点为起始点，定其角度为零度，以顺时钟或逆时钟方向，以车容量为限制条件进行服务区域之分割，再建构车辆排程路线。
2-opt*算子：属于局部搜索算法，局部搜索算法是解决组合优化问题的有效工具。
2-opt局部搜索算子：基本思想为在一条路径中随机选取两个点，将这两点之间的路径翻转，从而得到一条新路径，不断迭代运行以期通过搜索得到更好的解。
relocate局部搜索算子：基本思想为将某条路径中点重新分配到其他路径中，不断迭代运行以期通过搜索得到更好的解
swap局部搜索算子：基本思想为交换相同或不同路径中的两个点的位置以得到更好的线路排布。
如图1所示，本发明提供了一种解决多配送中心两级城市协同配送的大邻域搜索方法，包括以下步骤：
S1：根据客户点与网点之间的距离，构造初始可行解；
S2：对初始可行解进行大邻域搜索，确定最新第二级配送网络；
S3：对最新第二级配送网络进行优化；
S4：基于优化后的最新第二级配送网络，对市趟运输线路进行重构，得到可行市趟运输线路；
S5：对可行市趟运输线路进行优化，并对由配送中心与网点之间的车辆运输线路构成的一级路径进行优化，得到最新可行解；
S6：根据最新可行解确定当前解，并根据模拟退火准则判断是否达到算法终止条件，若是则结束大邻域搜索，否则返回步骤S2。
在本发明实施例中，步骤S1包括以下子步骤：
S11：根据客户点与网点之间的距离，利用轮盘赌选择法将客户点分配至各网点，确定各网点对各配送中心的货物需求量；
S12：根据各网点对各配送中心的货物需求量，利用扫描算法确定各配送中心与各网点的市趟运输线路和车辆调度；
S13：根据市趟运输线路和车辆调度中车辆到达网点的时间，利用扫描算法确定各网点到达各客户点的最后一公里运输线路，确定初始可行解。
在满足网点作业班次时间窗(即当某一班次网点作业的时间间隔)以及客户时效要求下，综合使用轮盘赌选择法和扫描算法构造一个初始可行解。
对于初始解的构造，首先基于两点之间的导航距离使用轮盘赌选择法将客户点分配给各网点，从而确定了各网点对各配送中心的货物需求量，随即使用扫描算法在满足客户时效的要求下规划配送中心到网点的市趟运输线路与车辆调度，再根据市趟运输车辆实际到达网点的时间使用扫描算法规划网点到客户的最后一公里运输线路，从而获得一个初始可行解。
在本发明实施例中，步骤S2中，确定最新二级配送网络的具体方法为：利用移除算子和插入算子对网点与客户点组成的二级配送网络进行大邻域搜索，确定最新第二级网络；其中，移除算子包括随机移除算子、最坏移除算子、相关移除算子和线路移除算子；插入算子包括贪婪插入算子、贪婪扰动插入算子、贪婪禁忌插入算子和遗憾插入算子。
从设计的移除和插入算子池中，根据各算子在之前算法迭代中的效果选择并使用一对移除和插入算子对第二级网络(即网点到客户组成的网络)进行大邻域搜索获得一个新的第二级网络。
在本发明实施例中，随机移除算子进行大邻域搜索的具体方法为：从当前解中随机移除若干个客户点，得到当前解的邻域；
最坏移除算子进行大邻域搜索的具体方法为：从当前解中移除若干个使配送成本降低最多的客户点，得到当前解的邻域；
相关移除算子进行大邻域搜索的具体方法为：迭代运行相关移除算子若干次，随机选择一个客户点作为种子客户，移除与该种子客户相关度量最大的客户点，并从当前解中移除若干个客户点，得到当前解的邻域；
线路移除算子进行大邻域搜索的具体方法为：迭代随机移除各网点到达各客户点的最后一公里运输线路以及最后一公里运输线路中的客户点，直至被移除的客户点达到设定阈值，得到当前解的邻域；
贪婪插入算子进行大邻域搜索的具体方法为：将被移除的若干个客户点重新插入二级配送网络中；
贪婪扰动插入算子进行大邻域搜索的具体方法为：将被移除的若干个客户点重新插入二级配送网络中，并对插入产生的配送成本进行扰动；
贪婪禁忌插入算子进行大邻域搜索的具体方法为：将被移除的若干个客户点重新插入二级配送网络中，且插入位置不与上一次迭代时客户点的运输线路位置重叠；
遗憾插入算子进行大邻域搜索的具体方法为：将遗憾值最大的客户点重新迭代插入二级配送网络中；
其中，将初始可行解作为大邻域搜索的初始当前解，并将大邻域搜索中迭代产生的最新解作为当前解。
在算法首次运行时，当前解就是初始解；算法每次迭代都是从一个当前可行解开始，随着算法不断地迭代运行，当前解也不断更新，此处描述移除算子基于算法的一般迭代过程，因此以当前解作为描述对象。
在本发明实施例中，贪婪插入算子依据当前插入成本最小原则，在满足客户时效要求前提下迭代地将被移除的客户重新插入第二级配送网络中形成完整、可行的第二级配送网络。
贪婪扰动插入算子与贪婪插入算子相似，区别在于将客户重新插入到当前解中时对产生的成本进行扰动，从而使得搜索得到的解具有多样性。
贪婪禁忌插入算子与贪婪插入算子相似，区别在于将客户重新插入到当前解中时禁止将该客户重新插入到上一次迭代时该客户所处的线路位置中，以增加解的多样性。
遗憾插入算子针对贪婪插入算子的短视性进行改进，迭代地将遗憾值最大的客户重新插入到第二级配送网络中。遗憾值的含义是将某客户插入到各个线路前后所产生的成本差值进行升序排列，计算第2到第k种插入方式与第1种插入方式的差值之和，该结果表示若不将该客户插入到使成本增加最少的位置所造成的遗憾。每次迭代所选择的待插入客户计算公式为：
Figure BDA0003549973410000071
然后将该客户插入使成本增加最少的第二级配送网络位置中，并满足客户时效要求，从而形成完整、可行的第二级配送网络。
在本发明实施例中，相关度量R(i,j)的计算公式为：
Figure BDA0003549973410000072
其中，di,j表示客户点i和客户点j之间的距离，
Figure BDA0003549973410000073
表示客户点i对一个配送中心的货物需求量，
Figure BDA0003549973410000074
表示客户点i对另一个配送中心的货物需求量，
Figure BDA0003549973410000075
表示客户j对一个配送中心的货物需求量，
Figure BDA0003549973410000076
表示客户j对另一个配送中心的货物需求量，
Figure BDA0003549973410000077
表示第一权重参数，λ表示第二权重参数，θ表示第三权重参数；
遗憾值最大的客户点i:的计算公式：
Figure BDA0003549973410000078
其中，D表示客户集合，k表示插入方式总数，fi h表示客户i的第h低的插入成本，fi 1表示客户i的最低成本。
在本发明实施例中，步骤S3中，对最新二级配送网络进行优化的具体方法为：以市趟运输线路的片段为单位，利用2-opt*算子在最新二级配送网络中进行组合搜索，利用Swap算子交换最新二级配送网络中的所有客户点，得到最优运输线路，完成最新二级配送网络优化。
2-opt*算子以线路的片段为单位，在所有第二级线路之间进行组合搜索以获得更优的线路排布。Swap算子尝试交换所有第二级线路的所有客户，以获得更优的线路排布。
在使用2-opt*和Swap算子时，使用一种预留时间Gap(间隔)的方法快速检查线路的时间可行性以满足客户时效要求。即我们计算在当前所有线路中的客户距离达到时效要求的Gap值，然后计算各配送中心直发各网点的市趟运输线路所需的时间，通过比较Gap值和直发市趟运输时间可以快速判定时间可行性，提高算法效率。
在本发明实施例中，步骤S4中，对市趟运输线路进行重构的具体方法为：将距离配送中心最远的网点插入至空闲市趟运输线路中，迭代将所有未规划的网点插入至空闲市趟运输线路中运输成本最少的线路位置处，并根据运输车辆的容量限制、网点的作业时间窗限制和客户点时效限制，将最新的市趟运输线路分割得到若干条可行市趟运输线路。
在本发明实施例中，步骤S5中，对可行市趟运输线路进行优化的具体方法为：利用网点移除算子将可行市趟运输线路中各运输车辆不同步到达的网点移除，利用网点插入算子将被移除的网点重新插入至可行市趟运输线路中，以使优化后的最新二级配送网络的运输成本最小，完成对可行市趟运输线路的优化；
使用网点移除算子和网点插入算子对市趟运输线路及车辆调度进行优化。网点移除算子将那些出现了各市趟运输车辆不同步到达的网点移除。网点插入算子以整体二级配送网络成本最小为原则将被移除的网点重新插入到市趟运输线路中，从而获得新的、完整、可行、更优的两级配送网络，并在整体网络成本最低的要求下尽可能的消除了各市趟运输车辆到达各网点的不同步性。
步骤S5中，利用2-opt局部搜索算子、relocate局部搜索算子和swap局部搜索算子对一级路径进行优化。
在本发明实施例中，步骤S6中，若最新可行解目标函数值低于设定的当前最优目标函数值，则将最新可行解作为当前最优解；整个算法运行过程中主要存在2种解：当前解、最优解，二者都随算法运行不断变化，也都分别对应一个目标函数值；最优目标函数值对应最优解的目标函数值，那么优于最优解不断变化，最优目标函数值也随之不断变化，因此当前目标函数值就是算法运行过程中当前得到的最优解对应的目标函数值。
若最新可行解目标函数值低于当前解目标函数值但高于当前最优目标函数值，则将最新可行解作为当前解；
若最新可行解目标函数值劣于当前解目标函数值，则根据最新可行解的接受概率
Figure BDA0003549973410000081
判断是否将最新可行解目标函数值作为当前解，其中，Tcur表示当前温度，f(snew)和f(scur)分别表示新解snew和当前解scur的目标函数值；程序会生成随机数以判断是否接受该新解；
步骤S6中，若模拟退火准则中设定的当前温度Tcur低于结束温度Tend，则输出当前最优解，否则返回步骤S2进入下一次迭代，直至结束大邻域搜索。
本发明的工作原理及过程为：本发明的基本思想是从一个初始可行解出发，使用一对契合多配送中心两级城市物流配送网络的市趟运输线路与车辆调度规划问题的算子得到当前解的邻域，再使用若干局部搜索算子提高该邻域解的质量，在得到第二级的线路与车辆调度方案后，使用重构方法对第一级的市趟运输线路与车辆调度进行重构使其成为可行解，再使用我们设计的网点移除和插入算子对各配送中心同步性要求进行优化并得到新的可行解，再使用若干局部搜索算子对该新可行解进行优化，根据模拟退火准则判定解的接受与否再进入下一次迭代，直到算法达到退出条件输出最优解。
本发明的有益效果为：
(1)本发明通过整体优化两级配送网络中的市趟运输线路、最后一公里线路以及车辆调度，降低了整体网络成本而不仅是某单一配送中心配送网络成本；
(2)本发明使用精心设计的网点移除和插入算子调度市趟运输车辆到达各网点的时间，消除了市趟运输车辆到达网点的时间超出网点作业班次时间窗所带来的网点运作痛点；
(3)本发明采用更优的线路排布和车辆调度使得车辆使用数量大幅减少，降低了企业的资金占用。
本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合，这些变形和组合仍然在本发明的保护范围内。

Claims (9)

1.一种解决多配送中心两级城市协同配送的大邻域搜索方法，其特征在于，包括以下步骤：
S1：根据客户点与网点之间的距离，构造初始可行解；
S2：对初始可行解进行大邻域搜索，确定最新第二级配送网络；
S3：对最新第二级配送网络进行优化；
S4：基于优化后的最新第二级配送网络，对市趟运输线路进行重构，得到可行市趟运输线路；
S5：对可行市趟运输线路进行优化，并对由配送中心与网点之间的车辆运输线路构成的一级路径进行优化，得到最新可行解；
S6：根据最新可行解确定当前解，并根据模拟退火准则判断是否达到算法终止条件，若是则结束大邻域搜索，否则返回步骤S2。
2.根据权利要求1所述的解决多配送中心两级城市协同配送的大邻域搜索方法，其特征在于，所述步骤S1包括以下子步骤：
S11：根据客户点与网点之间的距离，利用轮盘赌选择法将客户点分配至各网点，确定各网点对各配送中心的货物需求量；
S12：根据各网点对各配送中心的货物需求量，利用扫描算法确定各配送中心与各网点的市趟运输线路和车辆调度；
S13：根据市趟运输线路和车辆调度中车辆到达网点的时间，利用扫描算法确定各网点到达各客户点的最后一公里运输线路，确定初始可行解。
3.根据权利要求1所述的解决多配送中心两级城市协同配送的大邻域搜索方法，其特征在于，所述步骤S2中，确定最新二级配送网络的具体方法为：利用移除算子和插入算子对网点与客户点组成的二级配送网络进行大邻域搜索，确定最新第二级网络；其中，移除算子包括随机移除算子、最坏移除算子、相关移除算子和线路移除算子；插入算子包括贪婪插入算子、贪婪扰动插入算子、贪婪禁忌插入算子和遗憾插入算子。
4.根据权利要求3所述的解决多配送中心两级城市协同配送的大邻域搜索方法，其特征在于，所述随机移除算子进行大邻域搜索的具体方法为：从当前解中随机移除若干个客户点，得到当前解的邻域；
所述最坏移除算子进行大邻域搜索的具体方法为：从当前解中移除若干个使配送成本降低最多的客户点，得到当前解的邻域；
所述相关移除算子进行大邻域搜索的具体方法为：迭代运行相关移除算子若干次，随机选择一个客户点作为种子客户，移除与该种子客户相关度量最大的客户点，并从当前解中移除若干个客户点，得到当前解的邻域；
所述线路移除算子进行大邻域搜索的具体方法为：迭代随机移除各网点到达各客户点的最后一公里运输线路以及最后一公里运输线路中的客户点，直至被移除的客户点达到设定阈值，得到当前解的邻域；
所述贪婪插入算子进行大邻域搜索的具体方法为：将被移除的若干个客户点重新插入二级配送网络中；
所述贪婪扰动插入算子进行大邻域搜索的具体方法为：将被移除的若干个客户点重新插入二级配送网络中，并对插入产生的配送成本进行扰动；
所述贪婪禁忌插入算子进行大邻域搜索的具体方法为：将被移除的若干个客户点重新插入二级配送网络中，且插入位置不与上一次迭代时客户点的运输线路位置重叠；
所述遗憾插入算子进行大邻域搜索的具体方法为：将遗憾值最大的客户点重新迭代插入二级配送网络中；
其中，将初始可行解作为大邻域搜索的初始当前解，并将大邻域搜索中迭代产生的最新解作为当前解。
5.根据权利要求4所述的解决多配送中心两级城市协同配送的大邻域搜索方法，其特征在于，所述相关度量R(i,j)的计算公式为：
Figure FDA0003549973400000021
其中，di,j表示客户点i和客户点j之间的距离，
Figure FDA0003549973400000022
表示客户点i对一个配送中心的货物需求量，
Figure FDA0003549973400000023
表示客户点i对另一个配送中心的货物需求量，
Figure FDA0003549973400000024
表示客户j对一个配送中心的货物需求量，
Figure FDA0003549973400000025
表示客户j对另一个配送中心的货物需求量，
Figure FDA0003549973400000026
表示第一权重参数，λ表示第二权重参数，θ表示第三权重参数；
所述遗憾值最大的客户点i:的计算公式：
Figure FDA0003549973400000027
其中，D表示客户集合，k表示插入方式总数，fi h表示客户i的第h低的插入成本，fi 1表示客户i的最低成本。
6.根据权利要求1所述的解决多配送中心两级城市协同配送的大邻域搜索方法，其特征在于，所述步骤S3中，对最新二级配送网络进行优化的具体方法为：以市趟运输线路的片段为单位，利用2-opt*算子在最新二级配送网络中进行组合搜索，利用swap算子交换最新二级配送网络中的所有客户点，得到最优运输线路，完成最新二级配送网络优化。
7.根据权利要求1所述的解决多配送中心两级城市协同配送的大邻域搜索方法，其特征在于，所述步骤S4中，对市趟运输线路进行重构的具体方法为：将距离配送中心最远的网点插入至空闲市趟运输线路中，迭代将所有未规划的网点插入至空闲市趟运输线路中运输成本最少的线路位置处，并根据运输车辆的容量限制、网点的作业时间窗限制和客户点时效限制，将最新的市趟运输线路分割得到若干条可行市趟运输线路。
8.根据权利要求1所述的解决多配送中心两级城市协同配送的大邻域搜索方法，其特征在于，所述步骤S5中，对可行市趟运输线路进行优化的具体方法为：利用网点移除算子将可行市趟运输线路中各运输车辆不同步到达的网点移除，利用网点插入算子将被移除的网点重新插入至可行市趟运输线路中，以使优化后的最新二级配送网络的运输成本最小，完成对可行市趟运输线路的优化；
所述步骤S5中，利用2-opt局部搜索算子、relocate局部搜索算子和swap局部搜索算子对一级路径进行优化。
9.根据权利要求1所述的解决多配送中心两级城市协同配送的大邻域搜索方法，其特征在于，所述步骤S6中，若最新可行解目标函数值低于设定的当前最优目标函数值，则将最新可行解作为当前最优解；
若最新可行解目标函数值低于当前解目标函数值但高于当前最优目标函数值，则将最新可行解作为当前解；
若最新可行解目标函数值劣于当前解目标函数值，则根据最新可行解的接受概率
Figure FDA0003549973400000031
判断是否将最新可行解目标函数值作为当前解，其中，Tcur表示当前温度，f(snew)和f(scur)分别表示新解snew和当前解scur的目标函数值；
所述步骤S6中，若模拟退火准则中设定的当前温度Tcur低于结束温度Tend，则输出当前最优解，否则返回步骤S2进入下一次迭代，直至结束大邻域搜索。
CN202210260534.2A 2022-03-16 2022-03-16 一种解决多配送中心两级城市协同配送的大邻域搜索方法 Active CN114611806B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210260534.2A CN114611806B (zh) 2022-03-16 2022-03-16 一种解决多配送中心两级城市协同配送的大邻域搜索方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210260534.2A CN114611806B (zh) 2022-03-16 2022-03-16 一种解决多配送中心两级城市协同配送的大邻域搜索方法

Publications (2)

Publication Number Publication Date
CN114611806A true CN114611806A (zh) 2022-06-10
CN114611806B CN114611806B (zh) 2024-05-28

Family

ID=81863077

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210260534.2A Active CN114611806B (zh) 2022-03-16 2022-03-16 一种解决多配送中心两级城市协同配送的大邻域搜索方法

Country Status (1)

Country Link
CN (1) CN114611806B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116629586A (zh) * 2023-07-24 2023-08-22 青岛民航凯亚系统集成有限公司 一种基于alns机场保障车辆调度方法及系统

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130159206A1 (en) * 2011-12-14 2013-06-20 International Business Machines Corporation Dynamic vehicle routing in multi-stage distribution networks
CN109583650A (zh) * 2018-11-30 2019-04-05 浙江工商大学 一种电动车换电站选址和物流配送联合调度的方法
CN111798067A (zh) * 2020-07-17 2020-10-20 大连理工大学 基于自适应大邻域搜索算法的自动驾驶汽车配送路径规划方法
CN112700850A (zh) * 2020-12-30 2021-04-23 杭州电子科技大学 一种基于三阶段混合启发式算法的家庭医疗护理调度随机优化方法
CN112734188A (zh) * 2020-12-30 2021-04-30 杭州电子科技大学 一种基于两阶段混合元启发式算法的家庭医疗护理调度优化方法
CN112884409A (zh) * 2021-02-26 2021-06-01 上海东普信息科技有限公司 配送路线推荐方法、装置、设备及存储介质
CN113469473A (zh) * 2021-09-06 2021-10-01 华南理工大学 一种考虑逆向物流的同城配送路线规划方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130159206A1 (en) * 2011-12-14 2013-06-20 International Business Machines Corporation Dynamic vehicle routing in multi-stage distribution networks
CN109583650A (zh) * 2018-11-30 2019-04-05 浙江工商大学 一种电动车换电站选址和物流配送联合调度的方法
CN111798067A (zh) * 2020-07-17 2020-10-20 大连理工大学 基于自适应大邻域搜索算法的自动驾驶汽车配送路径规划方法
CN112700850A (zh) * 2020-12-30 2021-04-23 杭州电子科技大学 一种基于三阶段混合启发式算法的家庭医疗护理调度随机优化方法
CN112734188A (zh) * 2020-12-30 2021-04-30 杭州电子科技大学 一种基于两阶段混合元启发式算法的家庭医疗护理调度优化方法
CN112884409A (zh) * 2021-02-26 2021-06-01 上海东普信息科技有限公司 配送路线推荐方法、装置、设备及存储介质
CN113469473A (zh) * 2021-09-06 2021-10-01 华南理工大学 一种考虑逆向物流的同城配送路线规划方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
DORIAN DUMEZ 等: "A large neighborhood search approach to the vehicle routing problem with delivery options", 《TRANSPORTATION RESEARCH PART B》, vol. 144, 8 January 2021 (2021-01-08), pages 103 - 132, XP086462534, DOI: 10.1016/j.trb.2020.11.012 *
蒋丽: "循环取货物流配送路径优化及车辆数预测系统实现", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》, no. 1, 15 January 2021 (2021-01-15), pages 034 - 1071 *
贾圣阳: "A快递企业多产品城市快递两级车辆路径问题研究", 《万方学位论文》, 1 November 2023 (2023-11-01), pages 1 - 71 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116629586A (zh) * 2023-07-24 2023-08-22 青岛民航凯亚系统集成有限公司 一种基于alns机场保障车辆调度方法及系统

Also Published As

Publication number Publication date
CN114611806B (zh) 2024-05-28

Similar Documents

Publication Publication Date Title
CN110657816B (zh) 一种基于烟花算法的带硬时间窗的车辆路径问题规划方法
CN105513400B (zh) 动态规划出行路径的方法
Vareias et al. Assessing customer service reliability in route planning with self-imposed time windows and stochastic travel times
Yu et al. Transit route network design-maximizing direct and transfer demand density
Li-ying et al. Multiple charging station location-routing problem with time window of electric vehicle.
Mauttone et al. A multi-objective metaheuristic approach for the transit network design problem
CN107909228B (zh) 基于模因计算的动态车辆收发货路径规划方法及装置
Ghoseiri et al. A hybrid genetic algorithm for multi-depot homogenous locomotive assignment with time windows
CN110097218B (zh) 一种时变环境下无人商品配送方法及系统
CN108256969A (zh) 一种公共自行车租赁点调度区域划分方法
Bakas et al. A flexible transportation service for the optimization of a fixed-route public transport network
Brands et al. Performance of a genetic algorithm for solving the multi-objective, multimodel transportation network design problem
CN113887782A (zh) 一种面向维修资源配送调度的遗传-烟花混合方法及系统
CN114707693A (zh) 一种车辆配送路径规划方法及系统
CN101782986A (zh) 一种基于免疫算法的物流配送分区均衡优化方法
CN111709561A (zh) 一种解决动态车辆路径问题的实时添加和优化方法
CN114611806A (zh) 一种解决多配送中心两级城市协同配送的大邻域搜索方法
CN115879657A (zh) 考虑多站点容量设计的电动车换电站选址路径优化方法
CN113344267A (zh) 一种基于合作的物流网络资源配置优化方法
CN112036623A (zh) 一种横向物流联盟的利益协调方法
Aghalari et al. A two-stage stochastic location–routing problem for electric vehicles fast charging
CN115713207A (zh) 一种混合公交服务组合优化方法
CN118228901A (zh) 基于两阶段启发式算法的两级生鲜车辆路径规划方法
Shen et al. A MultiObjective optimization approach for integrated timetabling and vehicle scheduling with uncertainty
Migliavacca et al. The REALISEGRID cost-benefit methodology to rank pan-European infrastructure investments

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant