首页 > 专利 > 杭州电子科技大学 > 二型模糊环境下的危化品车辆路径规划方法专利详情

二型模糊环境下的危化品车辆路径规划方法   0    0

有效专利 查看PDF
专利申请流程有哪些步骤?
专利申请流程图
申请
申请号:指国家知识产权局受理一件专利申请时给予该专利申请的一个标示号码。唯一性原则。
申请日:提出专利申请之日。
2019-08-22
申请公布
申请公布指发明专利申请经初步审查合格后,自申请日(或优先权日)起18个月期满时的公布或根据申请人的请求提前进行的公布。
申请公布号:专利申请过程中,在尚未取得专利授权之前,国家专利局《专利公报》公开专利时的编号。
申请公布日:申请公开的日期,即在专利公报上予以公开的日期。
2019-12-27
授权
授权指对发明专利申请经实质审查没有发现驳回理由,授予发明专利权;或对实用新型或外观设计专利申请经初步审查没有发现驳回理由,授予实用新型专利权或外观设计专利权。
2022-08-26
预估到期
发明专利权的期限为二十年,实用新型专利权期限为十年,外观设计专利权期限为十五年,均自申请日起计算。专利届满后法律终止保护。
2039-08-22
基本信息
有效性 有效专利 专利类型 发明专利
申请号 CN201910776617.5 申请日 2019-08-22
公开/公告号 CN110533238B 公开/公告日 2022-08-26
授权日 2022-08-26 预估到期日 2039-08-22
申请年 2019年 公开/公告年 2022年
缴费截止日
分类号 G06Q10/04G06Q10/08 主分类号 G06Q10/04
是否联合申请 独立申请 文献类型号 B
独权数量 1 从权数量 0
权利要求数量 1 非专利引证数量 0
引用专利数量 6 被引证专利数量 0
非专利引证
引用专利 CN107451693A、CN109948855A、CN109086914A、CN104933474A、CN103761588A、DE102014215473A1 被引证专利
专利权维持 3 专利申请国编码 CN
专利事件 事务标签 公开、实质审查、授权
申请人信息
申请人 第一申请人
专利权人 杭州电子科技大学 当前专利权人 杭州电子科技大学
发明人 蒋鹏、门金坤、许欢 第一发明人 蒋鹏
地址 浙江省杭州市经济技术开发区白杨街道2号大街 邮编 310018
申请人数量 1 发明人数量 3
申请人所在省 浙江省 申请人所在市 浙江省杭州市
代理人信息
代理机构
专利代理机构是经省专利管理局审核,国家知识产权局批准设立,可以接受委托人的委托,在委托权限范围内以委托人的名义办理专利申请或其他专利事务的服务机构。
浙江千克知识产权代理有限公司 代理人
专利代理师是代理他人进行专利申请和办理其他专利事务,取得一定资格的人。
周希良
摘要
本发明公开了一种二型模糊环境下的危化品车辆路径规划方法。本发明步骤是针对运输风险的不确定性,定义了危化品运输车辆路径规划模型。目标是确定风险最小化的运输路线。由于人员的流动性,本发明在传统运输风险模型的基础上引入了二型模糊变量,根据置信度方法构建了机会约束模型以及其对应的等价确定型。针对模型特性,设计了一种模拟退火算法。提出的SAA以一定的概率来接受一个比当前解要差的解,因此有可能会跳出这个局部的最优解,达到全局的最优解。本发明方法具有开放性、灵活性以及计算复杂度低等特点。
  • 摘要附图
    二型模糊环境下的危化品车辆路径规划方法
  • 说明书附图:图1
    二型模糊环境下的危化品车辆路径规划方法
  • 说明书附图:图2
    二型模糊环境下的危化品车辆路径规划方法
  • 说明书附图:图3
    二型模糊环境下的危化品车辆路径规划方法
  • 说明书附图:图4
    二型模糊环境下的危化品车辆路径规划方法
  • 说明书附图:图5
    二型模糊环境下的危化品车辆路径规划方法
法律状态
序号 法律状态公告日 法律状态 法律状态信息
1 2022-08-26 授权
2 2019-12-27 实质审查的生效 IPC(主分类): G06Q 10/04 专利申请号: 201910776617.5 申请日: 2019.08.22
3 2019-12-03 公开
权利要求
权利要求书是申请文件最核心的部分,是申请人向国家申请保护他的发明创造及划定保护范围的文件。
1.二型模糊环境下的危化品车辆路径规划方法,其特征在于,该方法具体包括以下步骤:
步骤1:获取基础数据,其中包括运输车辆信息、运输道路信息、人口分布和危化品信息;
步骤2:构建危化品车辆运输路径规划模型;
将危化品运输路径规划模型定义在一个完整的有向图G=(N,L)中;设N={0,1,2,…,n}是有向图中的节点集,节点0是仓储节点,C={1,2,…,n}为客户节点集,qi是节点i对危化品的需求量;L为有向图中的弧集,设arcij∈L是连接节点i之间和节点j的弧,dij为arcij的弧长;K={k1,k2,...,k|K|}是运输车辆集,每辆车k∈K有固定的负载能力限制Qk;
①计算运输风险
根据“概率‑后果”框架,运输风险被定义为事故概率与事故后果的乘积,任意两节点间的运输风险如下:
Rij=Pij×Csij,i,j∈N
式中,Rij是节点i,j∈N间的运输风险,Pij是弧arcij上的事故概率;Csij是弧arcij∈L上的事故后果;
②引入梯形区间二型模糊变量
将人口密度设为一个梯形区间二型模糊变量 如下:
式中, 与 为两个梯形一型模糊变量, 的上、下层隶属度函数分别为:
与 为上
层隶属度函数的参数, 为下层隶属度函数的参数, 和 分别为 和 的
高;
则,弧arcij上的运输风险计算如下:
式中, 为引入梯形区间二型模糊变量后的节点i,j∈N间的运输风险;
③危化品车辆运输路径规划模型
模型决策变量xijk如下:
模型目标函数如下:
模型约束如下:
运输车辆必须从仓储节点出发最终回到仓储节点:
每一个顾客的需求都要被满足,且只能被服务一次:
车辆不能超载:
使用的运输车辆数量不能超过|K|:
∑k∈K∑j∈Cx0jk≤|K|
④机会约束规划模型
根据梯形区间二型模糊变量 的上、下层隶属度函数,采用置信度方法将上述模型转换为两个机会约束规划模型;
对于上层隶属度函数, 机会约束规划模型如下:
对于下层隶属度函数, 机会约束规划模型如下:
U L U L
式中,Z 和Z 为机会约束规划模型的目标函数,α和α是预定义的置信度水平;现记,等价约束如下:
式中, 和 定义如下:
上述机会约束规划模型的等价确定性模型如下:
最终,可将原危化品车辆运输路径规划模型转化为:
步骤3:模型求解
根据模型特性,采用模拟退火算法求解危化品车辆运输路径规划模型的等价确定型;
其中模拟退火算法由一个初始解开始,采用邻域算子在初始解的邻域中搜寻新解;
Metropolis准则用来判断新解是否可以替换当前解;若当前解优于最优解,则采用当前解替换最优解;当达到最大内部迭代数时,当前温度将按照预定义的降温率下降;不断重复上述过程直至满足停止迭代标准;
弧arcij上的事故概率Pij计算如下:
Pij=ARij×Prij×dij
式中,ARij是弧arcij上的事故率;Prij是弧arcij上的危化品泄漏事故概率;
弧arcij∈L上的事故后果Csij计算如下:
式中,pdij是事发地点周边的人口密度; 是事故影响半径。
说明书

技术领域

[0001] 本发明属于危化品风险管理领域,涉及到自动化技术,特别是涉及一种二型模糊环境下危化品车辆路径规划问题的机会约束规划方法。

背景技术

[0002] 随着我国工业的快速发展,危化品已经成为工业生产中必不可少的重要组成部分。由于危化品的特殊性质,任何与其使用相关的活动,都伴随着巨大的风险。在危化品运输过程中,普通交通事故所导致危化品泄漏事故的概率极高,此类危化品运输事故可造成
大规模人员伤亡,环境恶化与财产损失。
[0003] 由于工业发展需要,危化品运输风险是无法避免的,只能通过一系列风险管理措施降低事故概率与事故后果,危化品运输路径规划是主要的运输风险管理措施之一。许多
传统的方法将危化品运输问题视为一个确定性问题,忽略了运输风险的不确定性,因而距
离实际应用还存在很大差距。此外,在采用分支定界法、割平面法以及PILP等精确算法在求解此类大规模优化问题时容易产生维数灾难,因此,危化品运输路径规划非常困难。

发明内容

[0004] 本发明的目标是针对危化品运输路径规划中的一些难题,在满足所有给定的约束条件的背景下,确定运输风险最小的车队运输路线。
[0005] 由于不确定性会导致运输风险的显著差异,本发明的技术方案是在危化品运输路径规划模型中引入二型模糊变量,通过置信度方法,提出一种机会约束规划方法将不确定
性模型转换为其等价确定型。根据模型特性,设计一种模拟退火算法求解问题,最终确立了二型模糊环境下危化品车辆路径规划问题的机会约束规划方法。
[0006] 本发明具体包括以下步骤:
[0007] 步骤1:获取基础数据,其中包括运输车辆信息、运输道路信息、人口分布和危化品信息。
[0008] 步骤2:构建危化品车辆运输路径规划模型。
[0009] 将危化品运输路径规划模型定义在一个完整的有向图G=(N,L)中。N={0,1,2,…,n}是有向图中的节点集,节点0是仓储节点,C={1,2,…,n}为客户节点集,qi是节点i对危化品的需求量;L为有向图中的弧集,设arcij∈L是连接节点i和节点j的弧,dij为arcij的弧长;K={k1,k2,...,k|K|}是运输车辆集,每辆车k∈K有固定的负载能力限制Qk。
[0010] ①计算运输风险
[0011] 根据“概率‑后果”框架,运输风险被定义为事故概率与事故后果的乘积,任意两节点间的运输风险如下:
[0012] Rij=Pij×Csij,i,j∈N
[0013] 式中,Rij是节点i,j∈N间的运输风险,Pij是弧arcij上的事故概率;Csij是弧arcij∈L上的事故后果。
[0014] ②引入梯形区间二型模糊变量
[0015] 将人口密度设为一个梯形区间二型模糊变量 如下:
[0016]
[0017] 式中, 与 为两个梯形一型模糊变量, 的上、下层隶属度函数分别为: 与
为上层隶属度函数的参数, 为下层隶属度函数的参数, 和 分别
为 和 的高。
[0018] 则,弧arcij上的运输风险计算如下:
[0019]
[0020] 式中, 为引入梯形区间二型模糊变量后的节点i,j∈N间的运输风险。
[0021] ③危化品车辆运输路径规划模型
[0022] 模型决策变量如下:
[0023]
[0024] 模型目标函数如下:
[0025]
[0026] 模型约束如下:
[0027] 运输车辆必须从仓储节点出发最终回到仓储节点:
[0028]
[0029] 每一个顾客的需求都要被满足,且只能被服务一次:
[0030]
[0031] 车辆不能超载:
[0032]
[0033] 使用的运输车辆数量不能超过|K|:
[0034] ∑k∈K∑j∈Cx0jk≤|K|
[0035] ④机会约束规划模型
[0036] 根据梯形区间二型模糊变量 的上、下层隶属度函数,采用置信度方法将上述模型转换为两个机会约束规划模型。
[0037] 对于上层隶属度函数, 机会约束规划模型如下:
[0038]
[0039] 对于下层隶属度函数, 机会约束规划模型如下:
[0040]
[0041] 式中,ZU和ZL为机会约束规划模型的目标函数,αU和αL是预定义的置信度水平;现记,等价约束如下:
[0042]
[0043] 式中, 和 定义如下:
[0044]
[0045]
[0046] 上述机会约束规划模型的等价确定性模型如下:
[0047]
[0048]
[0049] 最终,可将原危化品车辆运输路径规划模型转化为:
[0050]
[0051] 步骤3:模型求解
[0052] 根据模型特性,采用模拟退火算法求解危化品车辆运输路径规划模型的等价确定型。其中模拟退火算法由一个初始解开始,采用邻域算子在初始解的邻域中搜寻新解。
Metropolis准则用来判断新解是否可以替换当前解。若当前解优于最优解,则采用当前解
替换最优解。当达到最大内部迭代数时,当前温度将按照预定义的降温率下降。不断重复上述过程直至满足停止迭代标准。
[0053] 本发明的有益效果:本发明结合危化品运输特性在传统的车辆路径规划方法的基础上,考虑运输风险的不确定性,构建了更贴近危化品运输实际情况的二型模糊环境下危
化品车辆路径规划问题的机会约束规划方法,本发明具有开放性、灵活性以及计算复杂度
低等特点。

实施方案

[0059] 以下结合附图对本发明作进一步说明。
[0060] 本发明方法具体是:
[0061] 步骤1:获取基础数据,其中包括运输车辆信息、运输道路信息、人口分布和危化品信息。
[0062] 步骤2:构建危化品车辆运输路径规划模型。
[0063] 本发明将危化品运输路径规划模型定义在一个完整的有向图G=(N,L)中。N={0,1,2,…,n}是有向图中的节点集,节点0是仓储节点,C={1,2,…,n}为客户节点集,qi,i∈C是节点i对危化品的需求量;L为有向图中的弧集,arcij∈L是节点i,j之间的运输路径,dij为arcij的弧长;K={k1,k2,...,k|K|)是运输车辆集,每辆车k∈K有固定的负载能力限制Qk。
模型约束如下:
[0064] 运输车辆必须从仓储节点出发最终回到仓储节点;
[0065] 每一个顾客的需求都要被满足;
[0066] 每一个顾客只能被服务一次;
[0067] 车辆不能超载;
[0068] 使用的运输车辆数量不能超过|K|。
[0069] 2‑1、计算运输风险
[0070] 根据“概率‑后果”框架,运输风险被定义为事故概率与事故后果的乘积,任意两节点间的运输风险如下:
[0071] Rij=Pij×Csij,i,j∈N
[0072] 式中,Rij是节点i,j∈N间的运输风险,Pij是弧arcij∈L上的事故概率;Csij是弧arcij∈L上的事故后果。
[0073] 弧arcij∈L上的事故概率Pij计算方法如下:
[0074] Pij=ARij×Prij×dij,i,j∈N
[0075] 式中,ARij是弧arcij∈L上的事故率;Prij是弧arcij∈L上的危化品泄漏事故概率。
[0076] 弧arcij∈L上的事故后果严重程度计算方法如下:
[0077]
[0078] 式中,pdij是事发地点周边的人口密度; 是事故影响半径,通常设为1公里。
[0079] 2‑2、引入二型模糊变量
[0080] 由于人口的流动性,本发明将人口密度设为一个梯形二型模糊变量 如下:
[0081]
[0082] 式中, 与 为两个梯形一型模糊变量, 的上、下层隶属度函数分别为: 与
为上层隶属度函数的参数, 为下层隶属度函数的参数, 和 分别为 和
的高。
[0083] 因此,弧arcij∈L上的运输风险计算如下:
[0084]
[0085] 2‑3、危化品车辆运输路径规划模型
[0086] 模型决策变量如下:
[0087]
[0088] 模型目标函数如下:
[0089]
[0090] 模型约束如下:
[0091] 运输车辆必须从仓储节点出发最终回到仓储节点
[0092]
[0093] 每一个顾客的需求都要被满足,且只能被服务一次
[0094]
[0095] 车辆不能超载;
[0096]
[0097] 使用的运输车辆数量不能超过|K|
[0098] ∑k∈K∑j∈Cx0jk≤|K|
[0099] 2‑4、机会约束规划模型
[0100] 根据区间二型模糊变量 的边界隶属度函数,本工作采用了置信度方法将上述模型转换为2个机会约束规划模型。
[0101] 对于上层隶属度函数, 机会约束规划模型如下:
[0102]
[0103] 对于下层隶属度函数, 机会约束规划模型如下:
[0104]U L
[0105] 式中,α 和α 是预定义的置信度水平;现记,等价约束如下:
[0106]
[0107] 式中, 和 定义如下:
[0108]
[0109]
[0110] 上述机会约束模型的等价确定性模型如下:
[0111]
[0112]
[0113] 最终,可将原危化品车辆运输路径规划模型转化为:
[0114]
[0115] 步骤3:模型求解
[0116] 根据模型特性,本发明提出了一个有效的模拟退火算法求解危化品车辆运输路径规划模型的等价确定型。算法流程图如图1。所提出的模拟退火算法由一个初始解开始,采用邻域算子在初始解的邻域中搜寻新解。Metropolis准则用来判断新解是否可以替换当前
解。若当前解优于最优解,则采用当前解替换最优解。当达到最大内部迭代数时,当前温度将按照预定义的降温率下降。算法之后会重复上述过程直至满足停止迭代标准。
[0117] ①初始解构造
[0118] 初始解由一个随机序列生成。用一个具体的例子说明了初始解的构造。如图2所示,三辆车的运输车队需要服务15个客户,即,K={k1,k2,k3},N={0,1,2…,15}以及C={1,
2,…,15}。一组随机序列为:[12,1,5,11,16,2,10,4,9,3,8,17,6,13,14,7,15]。现以随机序列中大于客户数量的数字为断点构成初始解。
[0119] ②目标函数
[0120] 为了加快算法的收敛速度,本发明在目标函数中引入惩罚因子如下:
[0121]
[0122]
[0123] 式中,δk表示车辆k超载程度;pf为所有车辆超载程度。
[0124] 带惩罚因子的目标函数如下:
[0125] Z*=Zr(1+pc*pf)
[0126] 式中,pc为惩罚系数,pc越大,算法对不可行解的容忍度越小。
[0127] ③邻域算子
[0128] 根据模型特性,本发明采用了三种邻域算子产生新解。
[0129] Swap算子:随机交换解对应的序列内的两个数字的位置,如图3所示,解χ={{0,12,1,5,11,0},{0,2,10,4,9,3,8,0},{0,6,13,14,7,15,0}}中的‘4’和‘7’的位置被Swap算子交换,则新解χnew={{0,12,1,5,11,0},{0,2,10,0},{0,9,3,8,4,6,13,14,7,15,0}}。
[0130] Reversion算子:逆反解对应序列内的随机区域,如图4所示,解χ={{0,12,1,5,11,0},{0,2,10,4,9,3,8,0},{0,6,13,14,7,15,0}}对应序列内的[4,9,3,8,17]区域被Reversion算子逆反,则新解χnew={{0,12,1,5,11,0},{0,2,10,0},{0,8,3,9,4,6,13,14,
7,15,0}}。
[0131] Insertion算子:随机将解对应序列中的一个数字插入到另一个位置,如图5,解χ={{0,12,1,5,11,0},{0,2,10,4,9,3,8,0},{0,6,13,14,7,15,0}}对应序列内的‘4’被重新插入到‘17’的后面,则χnew={{0,12,1,5,11,0},{0,2,10,9,3,8,0},{0,4,6,13,14,7,15,0}}。
[0132] ④Metropolis准则
[0133] Metropolis准则用于判定邻域算子作用到当前解χ所产生的新解χnew是否可替换* *
当前解,记Δ为当前解与新解之间的目标函数增量,Δ=Z (χnew)‑Z (χ)。对于最小化问题,Δ<0表明新解χnew优于当前解χ,则采用新解替换当前解。若Δ>0,当前解被替换掉的概率为exp(‑(Δ/T)),T为算法当前温度。

附图说明

[0054] 图1为算法流程图;
[0055] 图2为初始解构造示意图;
[0056] 图3为Swap算子示意图;
[0057] 图4为Reversion算子示意图;
[0058] 图5为Insertion算子示意图。
版权所有:盲专网 ©2023 zlpt.xyz  蜀ICP备2023003576号