CN103854072B - 用于路径选择的处理方法及系统 - Google Patents

用于路径选择的处理方法及系统 Download PDF

Info

Publication number
CN103854072B
CN103854072B CN201410097795.2A CN201410097795A CN103854072B CN 103854072 B CN103854072 B CN 103854072B CN 201410097795 A CN201410097795 A CN 201410097795A CN 103854072 B CN103854072 B CN 103854072B
Authority
CN
China
Prior art keywords
value
utility
section
point
transportation
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.)
Expired - Fee Related
Application number
CN201410097795.2A
Other languages
English (en)
Other versions
CN103854072A (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.)
Hunan Institute of Technology
Original Assignee
Hunan Institute of Technology
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 Hunan Institute of Technology filed Critical Hunan Institute of Technology
Priority to CN201410097795.2A priority Critical patent/CN103854072B/zh
Publication of CN103854072A publication Critical patent/CN103854072A/zh
Application granted granted Critical
Publication of CN103854072B publication Critical patent/CN103854072B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/10Internal combustion engine [ICE] based vehicles
    • Y02T10/40Engine management systems

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Traffic Control Systems (AREA)

Abstract

本发明公开了一种用于路径选择的处理方法及系统,该方法包括:构建路径优化模型,路径优化模型包括:各路段行程时间的目标函数、运输成本的目标函数和安全通过概率的目标函数;根据多属性决策理论将各目标函数进行效用转换,分别算出各路段的行程时间、运输成本和安全通过概率的最大值和最小值;经线性变换确定各路段的时间效用值、运输成本效用值和安全通过概率效用值;将时间效用值、运输成本效用值和安全通过概率效用值分别进行无量纲处理转换成相应的无量纲值;通过对各无量纲值加权将路径优化转换成单属性决策;采用最短路径优化算法求解出决策效用值最大的通行路线。本发明简化了路径选择的过程,并能充分满足应急救灾物资的运输要求。

Description

用于路径选择的处理方法及系统
技术领域
本发明涉及交通运输领域,特别地,涉及一种用于路径选择的处理方法及系统。
背景技术
城市道路交通系统是城市复杂巨系统的重要组成部分,它直接影响着整个城市的经济发展及居民日常生活活动。城市道路交通网络复杂交织,当地震灾害发生时,不同程度的道路损害以及坍塌建筑物覆盖在路面上,都将严重影响整个道路网络的通行能力,进而降低通行效率。地震灾害发生时,不仅会造成直接的经济损失,倘若因救援工作不及时引起的间接损失也是非常严重的,因此,最优路径选择是应急救灾中的一个关键问题。
而目前常规的路径选择模型是对一系列出发点和需求点,选择适当的行车路径,在满足一定的约束条件的情况下,达到一定的优化目标,优化目标通常是运输的费用、距离或时间等。如公开号为CN103246932A的中国专利文献《一种基于利润最大化的货车最优路径搜索方法与装置》,就公开了一种针对货运路径的处理方法与装置,以使得每次货运利润最大化,提高运输生产计划制定的精准度,充分发挥运输资源能力,减少物流运营成本。但是,应急救灾车辆路径选择模型的目标是要在尽可能短的时间里赶到救灾现场并且尽最大可能的满足所有需求,不惜以运费为代价。且在震后环境下,我们要充分考虑道路交通的通行能力变化对救援车辆行程时间的影响,另外,应急救援最优路径并非是最短路径,而需要综合考虑资源运输的时效性和经济性两因素,从而最大化地满足灾后应急救灾物质的运输要求,以减少救援不及时导致的直接和间接损失。
发明内容
本发明目的在于提供一种针对应急救灾的用于路径选择的处理方法及系统,以解决现有的应急救灾运输路径选择难以满足时效性、经济性及安全性的综合要求的技术问题。
为实现上述目的,本发明采用的技术方案如下:
一种用于路径选择的处理方法,用于从通行网络中根据起始点和目的点寻找最佳的通行路线,包括以下步骤:
针对路径通行网络构建路径优化模型,路径优化模型包括:各路段行程时间的目标函数、各路段运输成本的目标函数和各路段安全通过概率的目标函数;
根据多属性决策理论将各目标函数进行效用转换,分别算出各路段的行程时间、运输成本和安全通过概率的最大值和最小值;
经线性变换确定各路段的时间效用值、运输成本效用值和安全通过概率效用值;
将时间效用值、运输成本效用值和安全通过概率效用值分别进行无量纲处理转换成相应的无量纲值;通过对各无量纲值加权将路径选择的优化转换成单属性决策;
采用最短路径优化算法求解出决策效用值最大的通行路线。
进一步地,行程时间的目标函数为:
上式中,Tst表示通过救灾道路<Vs,Vt>所用的时间,Xst是表示道路<Vs,Vt>是否在救灾车辆选定的路径上的决策变量,其中,Vs代表起点,Vt代表目的点;T为救灾车辆到灾区目的点运输线路的总时间;s为代表起点的整数,t为代表目的点的整数。
进一步地,运输成本的目标函数为:
上式中,Cst表示通过救灾道路<Vs,Vt>所用的运输成本,Xst是表示道路<Vs,Vt>是否在救灾车辆选定的路径上的决策变量,其中,Vs代表起点,Vt代表目的点;C为救灾车辆到灾区目的点运输线路的总成本;s为代表起点的整数,t为代表目的点的整数。
进一步地,安全通过概率的目标函数为:
上式中,Pst表示通过救灾道路<Vs,Vt>的安全通过概率,Xst是表示道路<Vs,Vt>是否在救灾车辆选定的路径上的决策变量,其中,Vs代表起点,Vt代表目的点;P为救灾车辆到灾区目的点运输线路的安全通过概率;s为代表起点的整数,t为代表目的点的整数。
进一步地,时间效用值的计算公式为:其中,Gtst为通过线路Lst的时间效用值,Tmax为从起点到目的点线路上最长的行程时间,tst为过线路Lst的行程时间;
运输成本效用值的计算公式为:其中,Gcst为通过线路Lst的运输成本效用值,Cmax为从起点到目的点线路上最大的运输成本,Cst为过线路Lst的运输成本;
安全通过概率效用值的计算公式为:其中,Gpst为通过线路Lst的安全通过概率效用值,其中,p′st=lgpst,p′max=lg pmax,Pmax为从起点到目的点线路上最大的安全通过概率,Pst为过线路Lst的安全通过概率。
进一步地,最短路径优化算法采用改进型的迪克斯特拉算法,包括:
A、将起点标记为固定节点,与起点直接连接的节点标记为临时节点,网络中它节点为未标记节点;
B、判断出发时刻t0所属时段,并计算出该固定节点至临时节点的行程时间、运输成本和安全通过概率;
C、计算每一个临时节点I的d值,d值是指临时节点的决策效用值;
D、选择d值最大的临时节点,判断该点是否为目标节点,是则停止算法,否则,将该临时点标记为固定节点,并将与该固定节点直接相连的节点标记为临时节点,并计算新标志为临时节点的d值,直至找到目标节点,或者网络中所有节点都标记为固定节点。
进一步地,还包括:
采集道路状况的实时数据,通过处理后获得当前时段各路段行程时间的预测值,如果计算得到的车辆最优路径与初始方案不同,并且新路线的决策效用值与原路线的决策效用值相比超过预定阈值,则按新路线行驶;否则,维持原路线不变。
根据本发明的另一方面,还提供一种用于路径选择的处理系统,包括:
路径优化模型构建模块,用于针对救灾路径通行网络进行路径选择优化,包括:各路段行程时间的目标函数、各路段运输成本的目标函数和各路段安全通过概率的目标函数;
目标函数效用转换模块,用于根据多属性决策理论将各目标函数进行效用转换,分别算出各路段的行程时间、运输成本和安全通过概率的最大值和最小值;
线性变换模块,用于经线性变换确定各路段的时间效用值、运输成本效用值和安全通过概率效用值;
无量纲处理模块,用于将时间效用值、运输成本效用值和安全通过概率效用值分别进行无量纲处理转换成相应的无量纲值;通过对各无量纲值加权将路径选择的优化转换成单属性决策;
优化决策模块,用于采用最短路径优化算法求解出决策效用值最大的通行路线。
进一步地,优化决策模块采用改进型的迪克斯特拉算法,包括:
A、将起点标记为固定节点,与起点直接连接的节点标记为临时节点,网络中它节点为未标记节点;
B、判断出发时刻t0所属时段,并计算出该固定节点至临时节点的行程时间、运输成本和安全通过概率;
C、计算每一个临时节点I的d值,d值是指临时节点的决策效用值;
D、选择d值最大的临时节点,判断该点是否为目标节点,是则停止算法,否则,将该临时点标记为固定节点,并将与该固定节点直接相连的节点标记为临时节点,并计算新标志为临时节点的d值,直至找到目标节点,或者网络中所有节点都标记为固定节点。
进一步地,还包括动态调整模块,用于采集道路状况的实时数据,通过处理后获得当前时段各路段行程时间的预测值,如果计算得到的车辆最优路径与初始方案不同,并且新路线的决策效用值与原路线的决策效用值相比超过预定阈值,则按新路线行驶;否则,维持原路线不变。
本发明具有以下有益效果:
本发明路径选择处理方法及处理系统,在设计路径优化模型时,综合考虑了路段的行程时间、运输成本和安全通过概率三个属性,并通过将多属性问题转换成单属性问题,并采用最短路径优化算法求解优化路径,简化了路径选择的过程,并能充分满足灾后应急救灾物资的运输要求,使得救灾物资在保障安全通过概率的前提下,以最小的成本和最快的速度到达目的点。
除了上面所描述的目的、特征和优点之外,本发明还有其它的目的、特征和优点。下面将参照图,对本发明作进一步详细的说明。
附图说明
构成本申请的一部分的附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1是本发明优选实施例路径选择处理方法的步骤流程示意图;
图2是本发明优选实施例最短路径优化算法的流场示意图;
图3是本发明优选实施例单属性决策的模型求解流程图;
图4是本发明优选实施例中带有动态调整方案的步骤流程示意图;
图5是本发明具体实施例中救灾路径网路的节点示意图;
图6是本发明具体实施例中车辆行驶速度随时间变化的分段函数示意图;以及
图7是本发明优选实施例路径选择处理系统的原理方框示意图。
具体实施方式
以下结合附图对本发明的实施例进行详细说明,但是本发明可以由权利要求限定和覆盖的多种不同方式实施。
本发明的优选实施例提供了一种用于路径选择的处理方法,用于在灾后从通行网络中根据起始点和救灾目的点寻找最佳的通行路线。在突发灾害,如地震、泥石流等发生后,应急救援需要解决的最根本问题就是,如何快速、安全、高效将应急救灾物资运送到受灾点,并将受灾人员撤离受灾地区,送往安全地带,从而达到减少人员伤亡和财产损失的目的。因此,突发灾害下,在进行应急救援车辆路径选取时应重点保证选取的路径能够快速且安全的到达目的地。一般不会考虑选择的路径是否舒适性的因素。震后道路评估主要考虑行程时间和路网的可靠性两个因素,通过分析路段行程时间保证应急救援车辆的快速性,通过分析路网可靠性保证路径选取的安全性,并考虑行程距离及运输成本,以考虑经济性。
参照图1,本实施例中,路径的选择处理方法包括以下步骤:
步骤S10,针对救灾路径通行网络构建路径优化模型,路径优化模型包括:各路段行程时间的目标函数、各路段运输成本的目标函数和各路段安全通过概率的目标函数;
步骤S20,根据多属性决策理论将各目标函数进行效用转换,分别算出各路段的行程时间、运输成本和安全通过概率的最大值和最小值;
步骤S30,经线性变换确定各路段的时间效用值、运输成本效用值和安全通过概率效用值;
步骤S40,将时间效用值、运输成本效用值和安全通过概率效用值分别进行无量纲处理转换成相应的无量纲值;通过对各无量纲值加权将路径选择的优化转换成单属性决策;
步骤S50,采用最短路径优化算法求解出决策效用值最大的通行路线。
本实施例在设计路径优化模型时,综合考虑了路段的行程时间、运输成本和安全通过概率三个属性,并通过将多属性问题转换成单属性问题,并采用最短路径优化算法求解优化路径,简化了路径选择的过程,并能充分满足灾后应急救灾物资的运输要求,使得救灾物资在保障安全通过概率的前提下,以最小的成本和最快的速度到达目的点。
本实施例中,在构建路径优化模型时,假定以下条件:
①所有物资都可以混装,由于灾后应急物资运输量都较大,因此只考虑车辆满载情况;
②货物的需求量远大于供应量,每辆车只服务一个需求点;
③运输车辆规格一致,运载量相同。
根据震后路网条件变化实际情况,通过道路通行能力应急评估方法判定各路段的通行能力,分配震后用户交通需求量得出各路段的交通量,从而得出各路段单元的行程时间,并综合考虑行程时间、经济花费和安全性要求,建立震后应急救灾车辆的路径选择模型。
优选地,本实施例中,行程时间的目标函数为:
上式中,Tst表示通过救灾道路<Vs,Vt>所用的时间,Xst是表示道路<Vs,Vt>是否在救灾车辆选定的路径上的决策变量,其中,Vs代表起点,Vt代表目的点;T为救灾车辆到灾区目的点运输线路的总时间;s为代表起点的整数,t为代表目的点的整数。
路段行程时间是应急救援路径选择时一个非常重要的因素。突发灾害发生后,会导致救灾道路上短时交通量的骤然增加,道路交通延误增大,路段的行程时间相较于常态下发生了很大的变化。因此有必要对灾后道路的路段行程时间的计算进行分析。利用道路交通中车流波动理论,灾后救援道路的行驶时间主要与交通流中的流量、速度、密度有关。由于灾后发生后,有些道路没有受到影响,有些通行能力下降甚至完全被阻塞,因而路段行驶时间可分为四种情况:a、道路不受灾害影响,行程时间为正常下的行驶时间;b、灾害救援道路的通行能力降低,但不阻塞;c、灾害救援道路的通行能力降低,且产生阻塞;d、道路被完全损坏,行程时间无限大。根据路段对应的不同行驶时间,可以确定相应路段的行程时间;由于路段时间服从加法原则,故救灾车辆到灾区目的点运输线路的总时间T可通过行程时间的目标函数求得。
本实施例中,运输成本的目标函数为:
上式中,Cst表示通过救灾道路<Vs,Vt>所用的运输成本,Xst是表示道路<Vs,Vt>是否在救灾车辆选定的路径上的决策变量;C为救灾车辆到灾区目的点运输线路的总成本;s为代表起点的整数,t为代表目的点的整数。由于运输成本亦服从加法原则,故救灾车辆到灾区目的点运输线路的总成本C亦可通过运输成本的目标函数求得。
本实施例中,安全通过概率的目标函数为:
上式中,Pst表示通过救灾道路<Vs,Vt>的安全通过概率,Xst是表示道路<Vs,Vt>是否在救灾车辆选定的路径上的决策变量;P为救灾车辆到灾区目的点运输线路的安全通过概率;s为代表起点的整数,t为代表目的点的整数。此处,路段的安全通过概率是指网络上任意节点之间保持连通的概率,本实施例通过将安全通过概率进行指标转化,使得转化后的指标服从加法运算法则。地震发生后,对道路的破坏主要有:路面破坏、桥梁倒塌、建筑物倒塌等。考虑到地震发生后道路基本上会阻断,要靠机械设备进行抢修。这里只研究建筑物的倒塌所造成的路段阻断风险,针对这种风险,考虑以下几个方面:
a、建筑物的倒塌概率
影响建筑物的倒塌几率的因素主要有地震的等级以及建筑物的耐震设计强度。假设m为建筑物的抗震强度,但随着使用年限的增加,建筑物的耐震强度会随之降低。则建筑物的倒塌几率P1为:
其中:m0、mu分别代表预测地震规模的上、下限,β表示地震的等级。
b、路段的阻断严重性
建筑物倒塌是否会阻断道路,除了受建筑物本身破坏机率影响外,也受建筑物高度与其所处位置道路宽度两者之影响,建筑物高度越高或道路宽度越小的路段,其阻断发生机率将比其他路段高,因此可将街道调和比定义为路段阻断严重性。即:Si=Hi/Di,其中Si是指路段i的阻断调和比;Hi是指路段i的建筑物平均高度,Di是路段i的平均宽度。
c、路段的阻断曝光量
阻断曝光量是指新闻媒体实时跟踪报道的信息量和救援车辆参与救援后行驶途中的实际阻断路报量,因而路段的阻断曝光量是两者的相加,计算机公式可用Ei=Ni+Ji来计算。其中Ei是指路段i的阻断曝光量;Ni是指路段i的阻断新闻曝光量;Ji是指路段i的阻断参与救援车辆实时路况曝光量。
综合上面三个方面,路网可靠性Ri=Pi·Si·Ei,其中Ri是指路段i的可靠性风险机率值,Pi是指路段i的建筑物的倒塌概率,Si是路段i的阻断调和比,Ei是指路段i的阻断曝光量。路段的安全通过概率根据路网的可靠性风险几率值确定,路网的可靠性风险几率值越高,则路段的安全通过概率则越大。
本实施例中,设定运输线路Lst(表示起点到目的点)由K条路段串联而成,设车辆通过第i路段的安全通过概率为Pi,则车辆通过Lst的安全通过概率为:
P(Lst)=P1*P2*Pi......*Pk
对上式两边取对数后,再乘以-1,得到路段的安全通过概率的目标函数,该目标函数为:
上式中,Pst表示通过救灾道路<Vs,Vt>的安全通过概率,Xst是表示道路<Vs,Vt>是否在救灾车辆选定的路径上的决策变量;P为救灾车辆到灾区目的点运输线路的安全通过概率;s为代表起点的整数,t为代表目的点的整数。
参照图3,在步骤S20中,根据本实施例的行程时间的目标函数、运输成本的目标函数和安全通过概率的目标函数,并依多属性决策理论对各目标函数进行效用转换,分别计算出各个路段上的行程时间的最大值和最小值(Tmax和Tmin),运输成本的最大值和最小值(Cmax和Cmin),安全通过概率的最大值和最小值(Pmax和Pmin)。
在步骤S30中,时间效用值的计算公式为:其中,Gtst为通过线路Lst的时间效用值,Tmax为从起点到目的点线路上最长的行程时间,tst为过线路Lst的行程时间;
运输成本效用值的计算公式为:其中,Gcst为通过线路Lst的运输成本效用值,Cmax为从起点到目的点线路上最大的运输成本,Cst为过线路Lst的运输成本;
安全通过概率效用值的计算公式为:其中,Gpst为通过线路Lst的安全通过概率效用值,其中,p′st=lgpst,p′max=lg pmax,Pmax为从起点到目的点线路上最大的安全通过概率,Pst为过线路Lst的安全通过概率。
在步骤S40中,对各属性的效用值进行线性变换,以实现0-1的编码,具体如下:
对行程时间属性,令由于Tmin≤tst≤Tmax,所以Gtst∈[0,1]。说明tst越少,Tmax-tst越大,Gtst越大。即运输时间越小时效性越好的线路,其无量纲效用值Gtst越大;
对运输成本属性,令同上所述运输成本越少即经济成本越好的线路,其无量纲效用值Gcst越大。
对安全通过概率属性,令其无量纲效用值Gpst越大,线路安全通过的概率越大。
本实施例中,将各无量纲效用值加权,得到路径选择的单属性决策函数,该单属性决策函数为:
上式中Wt表示时间的权重因素,取值范围为[0,1];Wp表示安全性的权重因素,取值范围为[0,1];Wc表示成本的权重因素,取值范围为[0,1];且Wt+Wp+Wc=1。
步骤S50中,利用最短路径优化算法,通过该单属性决策函数,选取决策效用值最大的路径为最优的救灾通行路径。
优选地,参照图2,本实施例中,最短路径优化算法采用改进型的迪克斯特拉算法(Dijkstra算法),具体步骤如下:
A、将起点标记为固定节点,与起点直接连接的节点标记为临时节点,网络中它节点为未标记节点;
B、判断出发时刻t0所属时段,并计算出该固定节点至临时节点的行程时间、运输成本和安全通过概率;
C、计算每一个临时节点I的d值,d值是指临时节点的决策效用值;
D、选择d值最大的临时节点,判断该点是否为目标节点,是则停止算法,否则,将该临时点标记为固定节点,并将与该固定节点直接相连的节点标记为临时节点,并计算新标志为临时节点的d值,直至找到目标节点,或者网络中所有节点都标记为固定节点。
原始的Dijkstra算法是一种向后搜索的方法,原始的Dijkstra算法各个弧的权值是固定不变的,由于本实施例考虑了时变性,导致每增加一个后续节点,车辆的行驶速度、安全通过率都可能发生变化,即路网的权矩阵就会更新一次,因此本实施例相应地改变未标号节点的弧权,然后用改进Dijkstra算法求解各属性分量的最短最长路并进行无量纲处理,最后求解决策效用函数。
由于地震灾害发生后各种不确定因素可能对道路情况产生影响,即车辆行驶过程中的路况不可能完全符合历史规律,因此车辆在行驶途中需要根据实时交通信息并随时调整行驶路径,这里的实时交通信息通常包括交通流量信息、平均速度、拥堵状态等,而它们的采集及获取可借助GPS、计算机等信息技术来完成。优选地,本实施例中,还包括对通行路线进行动态调整的步骤,具体为:采集道路状况的实时数据,通过处理后获得当前时段各路段行程时间的预测值,如果计算得到的车辆最优路径与初始方案不同,并且新路线的决策效用值与原路线的决策效用值相比超过预定阈值,则按新路线行驶;否则,维持原路线不变。本实施例中,如果计算得到的车辆最优路径与初始方案不同,并且新路线的决策效用值比原路线的决策效用值增大达到或者超过原决策效用值的30%,那么就接受新的路径方案,按新路线行驶;否则,维持原路线不变。具体的行驶途中车辆路径方案调整步骤如图4所示。
下面结合一个具体的应用来对本实施例路径选择的处理方法进行介绍:
1、问题描述:A地区平时道路行车车速一般在40km/h左右,但是由于在7:30-8:30和14:00-15:00这两个上下班高峰期内车流量太大,车辆行驶速度不超过25km/h,而凌晨0:00-5:00时车流量最小,车辆行驶速度超过60km/h。现假定某县为地震灾害发生地点,需要距灾害发生地最近的未发生灾害区火速赶赴灾区进行现场紧急救援。此问题优化目标是在时变条件下,选择一条最佳路径,使应急救援车辆行车到场的时效性、安全性、经济性总效用最大,为便于实施,图5给出对应的以道路交叉口为顶点,道路为边的有向网络图。要解决的问题即为源点1和终点17之间最佳车辆救援路径选择。
2、问题分析:
(1)属性分量分析
由于地震灾害发生后,常规的通信方式和路网信息已经完成破坏,灾区现场可能是一片混乱。救灾指挥中心对受灾区的信息完全不知情,而时间对于灾区来说尤其宝贵,救灾指挥中心或相邻的未发生灾害区接到灾情求救信息后,必须快速做出快速响应,以极短的响应时间到达受灾地区进行现场救灾,才能减少人员伤亡和财产损失。也就是说对于应急救援车辆路径选择问题其实是有时间窗口限制的。但实际上由于灾害发生后道路情况十分复杂、通往灾区的救援车辆也比较多或者是现场指挥调度不力等问题,很多情况下都是超过时间窗限制,故而在此案例中只要求时间最短即可。根据案例描述,可将总时间段大致分为6段:[00:00,05:00||05:00,07:30||07:30,08:30||08:30,14:00||14:00,15:00||15:00,24:00],图6是车辆行驶速度随时间变化的分段函数。
(2)数据分析处理
表1是从各个大时间段里选取的几个小时间段的相关数据作为路网信息,包括各路段的路长,各时段各路段的行驶速度、安全通过概率等。
表1路网信息
3、模型求解
对模型中时间、安全、经济属性的无量纲化指标求解公式分别为:
其中:为第y条路线行驶时间的无量纲指标,为第y条路线运输成本的无量纲指标,为第y条路线运输安全性的无量纲指标。
先分别计算出最短和最长时间Tmax、Tmin,最大和最小安全通过概率Pmax、Pmin,最大和最小行驶成本Cmax、Cmin。可采用改进Djkstra算法分别进行计算。通过最优路线对比分析得到如表2所示的不同时刻出发的最优路线选择模型。
表2不同时刻出发的最优路线选择模型
根据本发明的另一方面,还提供一种用于路径选择的处理系统,此处理系统与本发明方法实施例中的步骤相对应,参照图7,包括:
路径优化模型构建模块10,用于针对救灾路径通行网络进行路径选择优化,包括:各路段行程时间的目标函数、各路段运输成本的目标函数和各路段安全通过概率的目标函数;
目标函数效用转换模块20,用于根据多属性决策理论将各目标函数进行效用转换,分别算出各路段的行程时间、运输成本和安全通过概率的最大值和最小值;
线性变换模块30,用于经线性变换确定各路段的时间效用值、运输成本效用值和安全通过概率效用值;
无量纲处理模块40,用于将时间效用值、运输成本效用值和安全通过概率效用值分别进行无量纲处理转换成相应的无量纲值;通过对各无量纲值加权将路径选择的优化转换成单属性决策;
优化决策模块50,用于采用最短路径优化算法求解出决策效用值最大的通行路线。
优选地,优化决策模块采用改进型的迪克斯特拉算法,具体步骤如下:
A、将起点标记为固定节点,与起点直接连接的节点标记为临时节点,网络中它节点为未标记节点;
B、判断出发时刻t0所属时段,并计算出该固定节点至临时节点的行程时间、运输成本和安全通过概率;
C、计算每一个临时节点I的d值,d值是指临时节点的决策效用值;
D、选择d值最大的临时节点,判断该点是否为目标节点,是则停止算法,否则,将该临时点标记为固定节点,并将与该固定节点直接相连的节点标记为临时节点,并计算新标志为临时节点的d值,直至找到目标节点,或者网络中所有节点都标记为固定节点。
优选地,本发明路径选择处理系统,还包括动态调整模块60,用于采集道路状况的实时数据,通过处理后获得当前时段各路段行程时间的预测值,如果计算得到的车辆最优路径与初始方案不同,并且新路线的决策效用值与原路线的决策效用值相比超过预定阈值,则按新路线行驶;否则,维持原路线不变。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (4)

1.一种用于路径选择的处理方法,用于在地震灾后从通行网络中根据起始点和目的点寻找最佳的通行路线,其特征在于,包括以下步骤:
针对路径通行网络构建路径优化模型,所述路径优化模型包括:各路段行程时间的目标函数、各路段运输成本的目标函数和各路段安全通过概率的目标函数;在构建路径优化模型时,假定以下条件:①所有物资都能混装,由于灾后应急物资运输量都较大,仅考虑车辆满载情况;②货物的需求量远大于供应量,每辆车只服务一个需求点;③运输车辆规格一致,运载量相同;根据震后路网条件变化实际情况,通过道路通行能力应急评估方法判定各路段的通行能力,分配震后用户交通需求量得出各路段的交通量,从而得出各路段单元的行程时间,并综合考虑行程时间、经济花费和安全性要求,建立震后应急救灾车辆的路径选择模型;
根据多属性决策理论将各目标函数进行效用转换,分别算出各路段的行程时间、运输成本以及安全通过概率的最大值和最小值;
经线性变换确定各路段的时间效用值、运输成本效用值和安全通过概率效用值;
将所述时间效用值、运输成本效用值和安全通过概率效用值分别进行无量纲处理转换成相应的无量纲值;通过对各无量纲值加权将路径选择的优化转换成单属性决策;
采用最短路径优化算法求解出决策效用值最大的通行路线;
所述最短路径优化算法采用改进型的迪克斯特拉算法,包括:
A、将起点标记为固定节点,与起点直接连接的节点标记为临时节点,网络中它节点为未标记节点;
B、判断出发时刻t0所属时段,并计算出该固定节点至临时节点的行程时间、运输成本和安全通过概率;
C、计算每一个临时节点I的d值,所述d值是指所述临时节点的决策效用值;
D、选择d值最大的临时节点,判断该点是否为目标节点,是则停止算法,否则,将该临时点标记为固定节点,并将与该固定节点直接相连的节点标记为临时节点,并计算新标志为临时节点的d值,直至找到目标节点,或者网络中所有节点都标记为固定节点;
还包括:
采集道路状况的实时数据,通过处理后获得当前时段各路段行程时间的预测值,如果计算得到的车辆最优路径与初始方案不同,并且新路线的决策效用值与原路线的决策效用值相比超过预定阈值,则按新路线行驶;否则,维持原路线不变;
所述时间效用值的计算公式为:其中,Gtst为通过线路Lst的时间效用值,Tmax为从起点到目的点线路上最长的行程时间,tst为过线路Lst的行程时间;
所述运输成本效用值的计算公式为:其中,Gcst为通过线路Lst的运输成本效用值,Cmax为从起点到目的点线路上最大的运输成本,Cst为过线路Lst的运输成本;
所述安全通过概率效用值的计算公式为:其中,Gpst为通过线路Lst的安全通过概率效用值,其中,pst=lg pst,pmax=lg pmax,Pmax为从起点到目的点线路上最大的安全通过概率,Pst为过线路Lst的安全通过概率;
所述安全通过概率的目标函数为:
上式中,Pst表示通过救灾道路<Vs,Vt>的安全通过概率,Xst是表示道路<Vs,Vt>是否在救灾车辆选定的路径上的决策变量,其中,Vs代表起点,Vt代表目的点;P为救灾车辆到灾区目的点运输线路的安全通过概率;s为代表起点的整数,t为代表目的点的整数;路段选择需考虑路网可靠性,所述路网可靠性Ri=Pi·Si·Ei,其中Ri是指路段i的可靠性风险机率值,Pi是指路段i的建筑物的倒塌概率,Si是路段i的阻断调和比,Ei是指路段i的阻断曝光量;路段的安全通过概率根据路网的可靠性风险几率值确定,路网的可靠性风险几率值越高,则路段的安全通过概率则越大。
2.根据权利要求1所述的用于路径选择的处理方法,其特征在于,
所述行程时间的目标函数为:
上式中,Tst表示通过救灾道路<Vs,Vt>所用的时间,Xst是表示道路<Vs,Vt>是否在救灾车辆选定的路径上的决策变量,其中,Vs代表起点,Vt代表目的点;T为救灾车辆到灾区目的点运输线路的总时间;s为代表起点的整数,t为代表目的点的整数。
3.根据权利要求1所述的用于路径选择的处理方法,其特征在于,
所述运输成本的目标函数为:
上式中,Cst表示通过救灾道路<Vs,Vt>所用的运输成本,Xst是表示道路<Vs,Vt>是否在救灾车辆选定的路径上的决策变量,其中,Vs代表起点,Vt代表目的点;C为救灾车辆到灾区目的点运输线路的总成本;s为代表起点的整数,t为代表目的点的整数。
4.一种用于路径选择的处理系统,用于在地震灾后从通行网络中根据起始点和目的点寻找最佳的通行路线,其特征在于,包括:
路径优化模型构建模块,用于针对救灾路径通行网络进行路径选择优化,包括:各路段行程时间的目标函数、各路段运输成本的目标函数和各路段安全通过概率的目标函数;在构建路径优化模型时,假定以下条件:①所有物资都能混装,由于灾后应急物资运输量都较大,仅考虑车辆满载情况;②货物的需求量远大于供应量,每辆车只服务一个需求点;③运输车辆规格一致,运载量相同;根据震后路网条件变化实际情况,通过道路通行能力应急评估方法判定各路段的通行能力,分配震后用户交通需求量得出各路段的交通量,从而得出各路段单元的行程时间,并综合考虑行程时间、经济花费和安全性要求,建立震后应急救灾车辆的路径选择模型;
目标函数效用转换模块,用于根据多属性决策理论将各目标函数进行效用转换,分别算出各路段的行程时间、运输成本和安全通过概率的最大值和最小值;
线性变换模块,用于经线性变换确定各路段的时间效用值、运输成本效用值和安全通过概率效用值;
无量纲处理模块,用于将所述时间效用值、运输成本效用值和安全通过概率效用值分别进行无量纲处理转换成相应的无量纲值;通过对各无量纲值加权将路径选择的优化转换成单属性决策;
优化决策模块,用于采用最短路径优化算法求解出决策效用值最大的通行路线;
所述优化决策模块采用改进型的迪克斯特拉算法,包括:
A、将起点标记为固定节点,与起点直接连接的节点标记为临时节点,网络中它节点为未标记节点;
B、判断出发时刻t0所属时段,并计算出该固定节点至临时节点的行程时间、运输成本和安全通过概率;
C、计算每一个临时节点I的d值,所述d值是指所述临时节点的决策效用值;
D、选择d值最大的临时节点,判断该点是否为目标节点,是则停止算法,否则,将该临时点标记为固定节点,并将与该固定节点直接相连的节点标记为临时节点,并计算新标志为临时节点的d值,直至找到目标节点,或者网络中所有节点都标记为固定节点;
还包括动态调整模块,用于采集道路状况的实时数据,通过处理后获得当前时段各路段行程时间的预测值,如果计算得到的车辆最优路径与初始方案不同,并且新路线的决策效用值与原路线的决策效用值相比超过预定阈值,则按新路线行驶;否则,维持原路线不变;
所述时间效用值的计算公式为:其中,Gtst为通过线路Lst的时间效用值,Tmax为从起点到目的点线路上最长的行程时间,tst为过线路Lst的行程时间;
所述运输成本效用值的计算公式为:其中,Gcst为通过线路Lst的运输成本效用值,Cmax为从起点到目的点线路上最大的运输成本,Cst为过线路Lst的运输成本;
所述安全通过概率效用值的计算公式为:其中,Gpst为通过线路Lst的安全通过概率效用值,其中,pst=lg pst,pmax=lg pmax,Pmax为从起点到目的点线路上最大的安全通过概率,Pst为过线路Lst的安全通过概率;
所述安全通过概率的目标函数为:
上式中,Pst表示通过救灾道路<Vs,Vt>的安全通过概率,Xst是表示道路<Vs,Vt>是否在救灾车辆选定的路径上的决策变量,其中,Vs代表起点,Vt代表目的点;P为救灾车辆到灾区目的点运输线路的安全通过概率;s为代表起点的整数,t为代表目的点的整数;路段选择需考虑路网可靠性,所述路网可靠性Ri=Pi·Si·Ei,其中Ri是指路段i的可靠性风险机率值,Pi是指路段i的建筑物的倒塌概率,Si是路段i的阻断调和比,Ei是指路段i的阻断曝光量;路段的安全通过概率根据路网的可靠性风险几率值确定,路网的可靠性风险几率值越高,则路段的安全通过概率则越大。
CN201410097795.2A 2014-03-17 2014-03-17 用于路径选择的处理方法及系统 Expired - Fee Related CN103854072B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410097795.2A CN103854072B (zh) 2014-03-17 2014-03-17 用于路径选择的处理方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410097795.2A CN103854072B (zh) 2014-03-17 2014-03-17 用于路径选择的处理方法及系统

Publications (2)

Publication Number Publication Date
CN103854072A CN103854072A (zh) 2014-06-11
CN103854072B true CN103854072B (zh) 2017-12-12

Family

ID=50861702

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410097795.2A Expired - Fee Related CN103854072B (zh) 2014-03-17 2014-03-17 用于路径选择的处理方法及系统

Country Status (1)

Country Link
CN (1) CN103854072B (zh)

Families Citing this family (30)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105447995B (zh) * 2014-08-27 2018-05-15 日立金融设备系统(深圳)有限公司 装填现金自动计划装置及方法
CN104537439B (zh) * 2014-12-29 2017-10-10 中南大学 一种用于缓解轨道交通拥挤的最小成本路径和混合路径分配方法
CN104537836B (zh) * 2014-12-30 2016-11-30 北京通博科技有限公司 路段行驶时间分布预测方法
CN106209635B (zh) * 2015-04-30 2019-06-11 深圳市中兴微电子技术有限公司 一种加权多路径wcmp的路由控制方法、装置和系统
US9646493B2 (en) 2015-06-19 2017-05-09 International Business Machines Corporation Management of moving objects
US10749734B2 (en) 2015-07-07 2020-08-18 International Business Machines Corporation Management of events and moving objects
CN106570584A (zh) * 2016-11-02 2017-04-19 北京工商大学 基于改进帝国主义竞争算法的城市轨道交通乘客路径选择方法
US10535266B2 (en) 2017-06-21 2020-01-14 International Business Machines Corporation Management of mobile objects
US10600322B2 (en) 2017-06-21 2020-03-24 International Business Machines Corporation Management of mobile objects
US10585180B2 (en) 2017-06-21 2020-03-10 International Business Machines Corporation Management of mobile objects
US10504368B2 (en) 2017-06-21 2019-12-10 International Business Machines Corporation Management of mobile objects
US10546488B2 (en) 2017-06-21 2020-01-28 International Business Machines Corporation Management of mobile objects
US10540895B2 (en) * 2017-06-21 2020-01-21 International Business Machines Corporation Management of mobile objects
CN107330608B (zh) * 2017-06-27 2021-04-23 宁波创元信息科技有限公司 一种基于神经元网络技术的生产排程方法
CN109934411B (zh) * 2018-11-09 2022-11-04 安庆师范大学 考虑紧急度的应急救灾智能路径规划方法及装置
CN109801491B (zh) * 2019-01-18 2022-09-13 深圳壹账通智能科技有限公司 基于风险评估的智能导航方法、装置、设备及存储介质
CN109900291A (zh) * 2019-04-03 2019-06-18 南京林业大学 一种综合安全性和出行时间的机动车导航方法
CN110175799B (zh) * 2019-05-05 2022-06-07 武汉理工大学 一种应急物资配送路径支路限流方案生成方法
CN110608738B (zh) * 2019-08-14 2021-08-03 青岛科技大学 一种无人船全局气象航线动态规划方法和系统
CN110853349A (zh) * 2019-10-24 2020-02-28 杭州飞步科技有限公司 车辆调度方法、装置及设备
CN111784260A (zh) * 2020-07-14 2020-10-16 国网北京市电力公司 运输规划方法、装置、存储介质及处理器
CN111968396A (zh) * 2020-08-22 2020-11-20 刘秀萍 智能车联网驱动的应急救援车行驶优化方法
CN113178092B (zh) * 2021-03-16 2022-05-27 淮阴工学院 一种基于效用函数的限高线路选取方法
CN113420920B (zh) * 2021-06-22 2023-01-03 哈尔滨工业大学 一种应急资源运送路径与交通管制措施同步决策方法及系统
TWI789876B (zh) * 2021-08-19 2023-01-11 奇點無限有限公司 最低風險路線的產生方法
CN114240001B (zh) * 2021-11-09 2024-06-18 北京京东振世信息技术有限公司 物流路由网络确定方法和装置
CN115994635B (zh) * 2023-03-23 2023-06-16 广东鉴面智能科技有限公司 一种皮带最优卸料运输路径检测方法、系统和介质
CN117114262B (zh) * 2023-04-07 2024-07-12 天津大学 一种碳利用点供给源的选择与分配方法
CN117876196A (zh) * 2024-03-11 2024-04-12 西安国信物联技术有限公司 一种地质灾害信息管理系统
CN117952498A (zh) * 2024-03-26 2024-04-30 云南电投绿能科技有限公司 一种光伏运输车的配置方法、装置、设备及存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101315731A (zh) * 2008-07-17 2008-12-03 上海途锐信息技术有限公司 一种利用gps智能车载终端实现车辆监控调度的系统
CN101346000A (zh) * 2008-07-29 2009-01-14 方圆信通科技(北京)有限公司 基于gps和移动终端的提供本地移动地理信息服务的系统
CN101782982A (zh) * 2009-07-07 2010-07-21 上海海事大学 集装箱多式联运路径选择的多目标整数线性规划方法
CN101840533A (zh) * 2009-03-20 2010-09-22 中华电信股份有限公司 电子化物流车队管理系统及其方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101315731A (zh) * 2008-07-17 2008-12-03 上海途锐信息技术有限公司 一种利用gps智能车载终端实现车辆监控调度的系统
CN101346000A (zh) * 2008-07-29 2009-01-14 方圆信通科技(北京)有限公司 基于gps和移动终端的提供本地移动地理信息服务的系统
CN101840533A (zh) * 2009-03-20 2010-09-22 中华电信股份有限公司 电子化物流车队管理系统及其方法
CN101782982A (zh) * 2009-07-07 2010-07-21 上海海事大学 集装箱多式联运路径选择的多目标整数线性规划方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
交通网络最优路径分析研究;陈虓;《中国优秀硕士学位论文全文数据库 信息科技辑》;20130615;58-70页 *

Also Published As

Publication number Publication date
CN103854072A (zh) 2014-06-11

Similar Documents

Publication Publication Date Title
CN103854072B (zh) 用于路径选择的处理方法及系统
Kim et al. Integrating travel demand modeling and flood hazard risk analysis for evacuation and sheltering
WO2016034142A1 (zh) 一种电力应急处置物资的智能化调拨方法
Ouallane et al. Overview of road traffic management solutions based on IoT and AI
Wey Model formulation and solution algorithm of traffic signal control in an urban network
Yuan et al. Proposed framework for simultaneous optimization of evacuation traffic destination and route assignment
Madireddy et al. An agent based model for evacuation traffic management
Dubey et al. Internet of Things based adaptive traffic management system as a part of Intelligent Transportation System (ITS)
CN108520353A (zh) 一种综合交通运输网络安全风险点辨识方法
Khalid et al. EvacSys: A cloud-based service for emergency evacuation
Zhang et al. Analysis of dynamic road risk for pedestrian evacuation
CN205050341U (zh) 一种视频排队检测装置及城市交通告警控制系统
Lu et al. Bi-level programming model for dynamic reversible lane assignment
Hamza-Lup et al. Leveraging e-transportation in real-time traffic evacuation management
Kumar et al. Optimizing infrastructure enhancements for evacuation planning
CN104091460A (zh) 一种用于树状结构网的疏散控制方法及系统
Patel et al. Assessment of network traffic congestion through Traffic Congestability Value (TCV): a new index
Elmitiny et al. Emergency evacuation planning and preparedness of transit facilities: traffic simulation modeling
Yuan et al. Evacuation strategy optimization study based on system theory
CN106448134A (zh) 一种应急期间地面临时交通管制方案的优化方法及系统
Fujihara et al. On the use of congestion information for rerouting in the disaster evacuation guidance using opportunistic communication
CN106803148B (zh) 一种应用于智慧城市的数据安全管理系统
Zhao et al. Vehicle route assignment optimization for emergency evacuation in a complex network
Ren et al. Bus contraflow lane: improved contraflow approach in freeway evacuation
CN104658271B (zh) 一种交通状态监控方法及系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20171212

Termination date: 20200317