CN106981194B - 一种高速路网关键路段的识别方法 - Google Patents
一种高速路网关键路段的识别方法 Download PDFInfo
- Publication number
- CN106981194B CN106981194B CN201710301222.0A CN201710301222A CN106981194B CN 106981194 B CN106981194 B CN 106981194B CN 201710301222 A CN201710301222 A CN 201710301222A CN 106981194 B CN106981194 B CN 106981194B
- Authority
- CN
- China
- Prior art keywords
- community
- section
- key
- road
- network
- 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
Links
Classifications
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
Landscapes
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Traffic Control Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公布了一种高速路网中关键路段的识别方法,首先利用社群划分算法将路网划分成多个子社群网络;再通过概率假设方法对路网建模,在每个子社群网络中,各求解得到候选关键路段;再采用动态规划方法从候选关键路段中找到最终的关键路段,作为高速路网中的关键路段,使得投资所述关键路段的组合能够实现路网通行效率最高。本发明提供方法根据流量和路网特性挖掘路网中的关键路段,能够准确高效的挖掘出高速路网中的关键路段,解决关键路段挖掘问题,可以成为实际应用中挖掘高速路网中关键路段的有效手段,能够为交通部门的整体规划提供参考依据,为高速公路管理人员高效率管理高速路网提供支撑手段。
Description
技术领域
本发明涉及智能交通和复杂网络分析领域,尤其涉及一种对高速路网运营起关键作用的关键路段的识别方法。
背景技术
近年来,我国高速公路迅猛发展,大大缩短了省际之间、重要城市之间的距离,有效降低了生产运输成本,改变了人们的生活方式。伴随着人们对高速公路出行的愈发依赖,突发事故和道路自身因素导致的高速拥堵问题也引发了社会各界人士的关注,成为一种广泛讨论的社会现象。高速拥堵带来了经济损失,影响了人们出行安全性和舒适度,降低了出行体验。
高速路网中存在着某些极为重要的交通枢纽,这些路段(其中节点被认为是路段的一部分)对路网通行起到非常关键的作用。若在这些关键路段发生事故,会对周边大面积路网通行产生级联影响。当高速路网发生拥堵时,其影响可能会通过排队车辆依次传导至周边路段,形成复杂的传播过程。识别出这些起关键作用的路段,能够为高速公路管理人员高效率管理高速路网提供支撑手段。现有技术中,直观地可以通过计算经过某路段的车流量、以及某路段所对应收费站节点的度,从而对路段的重要性进行排序;或使用贪心思想,每次从路网中选取一条断路后对路网通行效率影响最大的路段,选出的前B条路段即为贪心方法挖掘到的关键路段。但是,利用度和流量计算排名的方法过于简单,特征单一,难以反映复杂路网的真实情况;贪心方法时间复杂度太高,应对一省路网都会花费较长的时间,处理效率低,针对大规模复杂路网的处理能力不足。
发明内容
为了克服上述现有技术的不足,本发明提供一种高速路网关键路段的识别方法,其中,高速路网关键路段(其中节点被认为是路段的一部分)是指高速路网中的极为重要的交通枢纽,即一旦发生交通拥堵或道路损坏,会对路网通行效率造成严重影响的路段;本发明采用概率方法建模,并通过社群划分算法,有效地降低计算复杂度,再将路网关键路段的识别问题转化为动态规划的背包问题,从而又快又准地挖掘出B条高速路网中影响路网通行效率的关键路段。本发明提供方法能够准确高效的挖掘出高速路网中的关键路段,解决关键路段挖掘问题,满足实际应用需求。
本发明的原理是:高速路网结构复杂,贪心等方法挖掘关键路段耗费大量时间,且会陷入局部最优解。因此,本发明首先对路网进行社群划分,优化时间效率,并进一步将路网信息构建成概率化模型,再用背包思路准确求解,从而给出一套高速路网中的关键路段识别方法。依据真实高速路网的特点,每条路段均有一定的概率发生交通拥堵或路面损坏事件,而当交管部门向某些路段充分部署人力和物力时,正常情况下该路段发生交通拥堵或路面损坏的概率将会有效降低;依据这一特点,我们对路网的通行能力采用概率的方式建模,帮助挖掘一批更加具有道路通行影响力的关键路段。本发明首先将包含N个收费站、M条路段的路网拓扑结构划分为l个社群;在每个社群内部分别挖掘B条候选关键路段,使得当交管部门的人力物力部署在这B条路段上时,路网通行时间最短;再将从K*B条候选关键路段中选择至多B条从而使全路网通行时间最短的问题转化为背包问题;最终我们将动态规划方法求解出的B条关键路段,作为最终挖掘到的关键路段输出。本发明所提出的方法能够准确高效的挖掘出有路网通行影响力的关键路段,解决高速路网中的关键路段识别问题,满足实际应用需求。
本发明提供的技术方案是:
一种高速路网中关键路段的识别方法,所述关键路段指当高速路网发生交通拥堵或道路损坏时,会对路网通行效率造成严重影响的路段,其中路网节点被认为是路段的一部分;所述方法首先利用社群划分算法将路网划分成K个社群;再通过概率假设对路网进行建模,在每个社群中,分别求解得到候选关键路段;再从这些候选关键路段中,采用动态规划中的背包问题求解方法,找到B条最终的关键路段,由此识别高速路网中的关键路段。
一种高速路网中关键路段的识别方法,为简化模型,本发明具体实施中设定每个社群中都选取相同数量的候选关键路段,最终得到相同数量的关键路段。
上述方法以当前高速路网拓扑结构(即收费站和路段的地理位置以及连接情况)、通行车辆历史数据(即每辆车的OD信息和每段途径路段的通行时间,OD表示一段行程的起点origin和终点destination情况)作为输入数据,根据投资哪些路段的组合能够充分利用有限的资源,实现路网通行效率最高,确定至多B条关键路段;具体包括如下步骤:
1)首先生成路网拓扑结构G=(V,E),计算边(即路段)上包含的属性;
2)对路网拓扑结构G进行社群划分,得到l个社群的集合C={C1,…,Cl};
3)在每个社群内部识别候选关键路段每个社群内选择B条,
4)利用动态规划的方法,挖掘出至多B条能够使得全路网道路通行效率达到最高的关键路段;
5)将这B条路段输出,作为最终挖掘出的高速路网中的关键路段。
针对上述高速路网中的关键路段识别方法,进一步地,
步骤1)所述生成路网拓扑结构及边上属性的方法,具体包含如下步骤:
11)提取被研究省份高速路数据中的所有N个收费站{v1,…,vN},形成路网拓扑结构G中的节点集合V;
12)提取被研究省份高速路数据中的所有M条路段{e1,…,eM},形成路网拓扑结构G中的边集合E;
13)提取和计算每条边ei上的属性信息包括:
131)相应路段的长度,记为
132)历史数据中车辆经过该条路段的平均通行时间,记为
133)历史数据中该条路段经过的车流量,记为
步骤2)所述对路网拓扑结构G进行社群划分的方法,具体包括如下步骤:
21)每个节点初始化为一个社群;
22)计算每个社群S的分区性能f(S),其中cs表示一个顶点在S中、一个顶点不在S中的边的个数;ms表示完全包含在S中的边的个数:
23)对于每个节点i,计算每个与i相邻的社群Ji;
24)对于i的每个相邻社群Ji,计算假设i从原本的社群中删除、加入Ji中后,i原本所在社群及Ji的分区性能的增益;
25)若i的所有相邻社群中,存在某些相邻社群在接纳i之后,使得步骤24)所述的社群分区性能增益大于0,则将i从原本的社群中删除、加入能使增益最大的相邻社群;若不存在,则本轮对节点i的调整结束,跳至网络中下一个节点,直至对网络中所有节点进行一轮调整;
26)重复步骤22)至25),直至所有节点所属的社群不再发生变化,从而得到l个社群的集合C={C1,…,Cl},每个社群包含若干节点,社群之间无交集;
步骤3)所述识别每个社群内部的候选关键路段的方法,具体包括如下步骤:
31)对于社群Ci,计算Ci的含有B条边的所有子图表示所有可能的候选投资路段的组合;
32)计算每种组合中,度大于1的节点的个数
33)计算每种组合发生的可能性由于真实世界中,当某条路段发生拥堵时,与它相连的路段有较大可能通过排队车辆受到拥堵的级联影响,因此倾向于选择有一定联系的B条边的组合:
其中,|Ei|表示社群Ci中的边数;表示从|Ei|条边中选出B条边的排列数;
34)假定每条路段都有一定概率发生拥堵或路面损坏的情况,称为固有事故率;为了简化模型,假设所有路段的固有事故率均为Ppr;
35)对于没有被投资的路段,其事故率就是固有事故率,即:
Pe=Ppr
36)对于被投资的路段,其事故率为固有事故率Ppr与投资收益率Pinv之差,即:
Pe=Ppr-Pinv
投资收益率一般为给定值,在本发明实施例取值为0.1;
37)按事故率对Ci中的每条路段是否发生拥堵或路面损坏等事故进行抽样,将抽样结果中发生事故的路段从Ci中去除,得到1次抽样后的仿真社群拓扑结构为了消除抽样的随机性,进行S次抽样,即得
38)设历史数据中有H辆车经过路网,对于其中每辆车ho,提取其OD信息及经过的所有路段序列zo;
39)计算每辆车ho经过当前仿真社群拓扑结构时的通行效率具体包括步骤391)~392):
391)若路段序列zo中的所有边均存在于第p次仿真社群拓扑结构则车辆ho的通行时间即等于历史数据中的实际通行时间;通行效率为车辆ho通过路径长度与通行时间之比;
392)若路段序列zo中,存在某些边不属于仿真社群拓扑结构则利用迪杰斯特拉算法(Dijkstra's Algorithm)重新在拓扑结构上计算车辆ho从起点到终点的仿真路段序列z′o;则车辆ho的通行效率为:
其中,eq表示仿真路段序列z′o中的第q条路段;表示第q条路段的长度;表示历史数据中车辆经过第q条路段的平均通行时间;
310)社群Ci内部的B条候选关键路段集合为:
社群Ci内部的B条候选关键路段,即为上式所求结果中所包含的B条边;
311)对步骤2)中划分出来的每个社群,进行如上步骤31)~310)的操作,得到每个社群内部的B条候选关键路段;
312)记录当步骤310)式子取最大值时,所对应B条候选关键路段上,每条路段上的通行效率,记为v(i)。
步骤4)所述在全路网中挖掘出至多B条最终的关键路段的动态规划方法,具体包括如下步骤:
41)从l个社群中的l*B条候选关键路段中,挖掘至多B条最终的关键路段;记这l*B条候选关键路段的投资花费为w(i),为简化模型可使w(i)均为1;
42)使用动态规划中的背包思想,首先初始化一个l*B行W列的二维数组Y,将元素初始化为-1;
43)本问题的递推公式为:
Y[i][j]=max(f[i+1][j],f[i+1][j-w[i]]+v[i])
44)以二维数组Y的两个维度分别作为递归函数的两个参数;参数i、j的含义为:从第i条候选路段开始,选择能使总通行效率增益最大的、且总花费不超过j的候选路段;
45)首先从Y[0][W]开始递归;
46)进入递归函数后,若当前二维数组中对应的元素Y[i][j]不等于-1,表示之前已经计算过;直接返回计算过的值;
47)若没有计算过,判断第i条候选路段的花费是否超过当前剩余花费,若不超过,则按照43)所述公式计算即可;若超过,则跳过第i条候选路段,计算Y[i+1][j];
48)直到递归函数完全结束,被选中的B条候选关键路段即为最终所求的全路网关键路段。
与现有技术相比,本发明的有益效果是:
本发明提供一种高速路网中的关键路段识别方法,先利用社群划分算法将路网划分成K个社群;再通过概率假设对路网进行建模,在每个社群中,各求解B条候选关键路段;再从这K*B条候选关键路段中,采用动态规划中的背包问题解决思路,找到B条最终的关键路段,由此识别高速路网中的关键路段。本发明提供方法能够准确高效的挖掘出高速路网中的关键路段,解决关键路段挖掘问题,可以成为实际应用中挖掘高速路网中关键路段的有效手段。
本发明根据流量和路网特性挖掘路网中的关键路段,能够为交通部门的整体规划提供参考依据。采用本发明识别出路网中的关键路段,可为高速公路管理人员高效率管理高速路网提供支撑手段。管理人员可据此将应急资源部署在关键路段附近,一旦发生交通事故或拥堵,可快速舒缓交通压力,减少经济和时间损失;此外,车主可根据本发明结果制定出行规划,有效规避潜在的拥堵因素;本发明能够为高速公路的建设提供一定的指导信息,对于当前易发生拥堵的关键路段,政府部门可考虑采取增设车道、维护路面等措施,长效解决拥堵问题;同时,本发明亦能够为高速公路应急管提供策略支持,增加配流效率、降低损耗。
附图说明
图1是本发明提供的高速路网中的关键路段识别方法的整体流程框图。
图2是本发明方法对路网拓扑结构进行社群划分的流程框图。
图3是本发明方法识别每个社群内部的候选关键路段的流程框图。
图4是从由各社群组成的全路网中挖掘出至多B条最终的关键路段的流程框图。
图5是本发明中递归函数的流程框图。
具体实施方式
下面结合附图,通过实施例进一步描述本发明,但不以任何方式限制本发明的范围。
本发明中,关键路段实质当发生交通拥堵或道路损坏时,会对路网通行效率造成严重影响的路段。当关键路段发生拥堵或道路损坏时,交管部门可以通过在路段配置人力物力等资源,有效降低该条路段发生拥堵或损坏的概率,我们称之为“投资”;但交管部门可用的资源有限,假设最多只能投资B条路段。本发明提供了一种高速路网中的关键路段识别方法,能够达到利用流量和路网拓扑结构挖掘到B条最关键的路段,当资源投资在这B条路段上时,路网通行效率最高的目的。
图1是本发明提供方法的整体流程图,首先利用社群划分算法将路网划分成K个社群;再通过概率假设对路网进行建模,在每个社群中,各求解B条候选关键路段;再从这K*B条候选关键路段中,采用动态规划中的背包问题解决思路,找到B条最终的关键路段。具体包括如下步骤:
1)生成路网拓扑结构,计算边(即路段)上包含的属性
因为本方法将对路网进行社群划分、背包问题等操作,故将网络拓扑结构及边的性质作为输入准备。
2)对路网拓扑结构G进行社群划分,得到一组社群的集合
本步骤为后续挖掘单个社群内部的候选关键路段提供依据,是优化时间效率的关键步骤。
图2是本发明的对路网拓扑结构进行社群划分的流程图。基本流程包括如下过程:
21)将每个节点初始化为一个社群;
22)计算每个社群S的分区性能;
23)对于每个节点i,计算每个与i相邻的社群Ji;
24)对于i的每个相邻社群Ji,计算假设i从原本的社群中删除、加入Ji中后,i原本所在社群及Ji的分区性能的增益;
25)若i的所有相邻社群中,存在某些相邻社群在接纳i之后,使得步骤24)所述的社群分区性能增益大于0,则将i从原本的社群中删除、加入能使增益最大的相邻社群;若不存在,则本轮对节点i的调整结束,跳至网络中下一个节点,直至对网络中所有节点进行一轮调整;
26)重复步骤22)至25),直至所有节点的社群不再发生变化,得到社群划分的结果;
3)识别每个社群内部的候选关键路段
在全路网中直接利用本步骤模型、挖掘B条关键路段需要遍历极多的路段组合,会带来指数级别的时间复杂度。因此,本发明首先将全路网划分为l个社群,在社群内部利用本步骤模型先挖掘出B条候选关键路段;由于社群只是全路网的一个子图,因此能够有效的降低时间复杂度。
图3是本发明的识别每个社群内部的候选关键路段的流程图。基本流程包括如下过程:
31)计算社群Ci的含有B条边的所有子图;
32)计算每种组合中,度大于1的节点的个数;
33)计算每种组合发生的可能性
34)计算固有事故率Ppr;
35)计算没有被投资的路段的事故率;
36)计算被投资的路段的事故率;
37)按事故率对Ci中的每条路段是否发生拥堵或路面损坏等事故进行抽样,将抽样结果中发生事故的路段从Ci中去除,得到一次抽样后的仿真社群拓扑结构为了消除抽样的随机性,进行S次抽样,得到
38)对于经过路网的每辆车ho,提取其OD信息及经过的所有路段序列zo;
39)计算每辆车ho经过当前仿真社群拓扑结构时的通行效率具体包括:
391)若路段序列zo中的所有边均存在于通行效率为车辆ho通过路径长度与历史数据中的实际通行时间之比;
392)若路段序列zo中,存在某些边不属于仿真社群拓扑结构则利用迪杰斯特拉算法(Dijkstra's Algorithm)重新在拓扑结构上计算车辆ho从起点到终点的仿真路段序列z′o;通行效率为路径长度与重新计算的路段的通行时间之比:
310)计算能够使社群Ci内部路网通行效率最大的B条候选关键路段集合
311)对每个社群,进行如上步骤31)~310)的操作;
312)记录步骤310)中取到最大值时,所对应的B条候选关键路段各自的通行效率;
4)全路网中挖掘出至多B条最终的关键路段
由于高速公路管理部门的投资有限,需要我们仅从这l*B条候选关键路段中选择出B条作为最终结果。这一问题与背包的思想相似,本步骤给出了用背包思路求解的递推公式、以及使用递归函数实现背包问题时的主要步骤;
图4是本发明的从全路网中挖掘出至多B条最终的关键路段的流程图。包括如下过程:
41)记这l*B条候选关键路段的投资花费为w(i),为简化模型使w(i)均为1;
42)申请一个l*B行W列的二维数组Y,将元素初始化为-1;
43)Y[i][j]=max(f[i+1][j],f[i+1][j-w[i]]+v[i])为本问题所对应的递推公式;
44)以二维数组Y的两个维度分别作为递归函数的两个参数;
45)首先从Y[0][W]开始递归;
46)若Y[i][j]不等于-1,表示之前已经计算过;直接返回计算过的值;
47)若没有计算过,判断第i+1条候选路段的花费是否超过当前剩余花费,若不超过,则按照43)所述公式计算Y[i][j];若超过,则跳过第i+1条候选路段,计算Y[i+1][j];
48)直到递归函数完全结束,被选中的B条候选关键路段即为最终所求;
通过上述方法,可以筛选出B条路段,即得到挖掘出的全路网关键路段。
本发明以下实施例针对中国安徽省高速路网,输入数据为路网中的所有路段、每条路段的位置信息、经过路网的所有车辆信息(包括车辆id、途经路段、通行时间)。通过本发明提供的高速路网中的关键路段识别方法,提取网络拓扑结构信息、划分社区、每个社群内部识别候选关键路段、并利用动态规划的方法,挖掘出至多B条能够使得全路网道路通行效率达到最高的关键路段。中国安徽省路网分布广泛,拥有海量数据资料。快速挖掘高速路网中的关键路段,能够帮助工作人员更好的进行路网交通管理。
首先,利用数据库工具,提取输入数据,包括路网中的所有路段、每条路段的位置信息、经过路网的所有车辆id、途经路段及通行时间。按照如下步骤,挖掘全路网关键路段:
步骤一:存储网络拓扑结构及边的信息;
步骤二:将每个收费站节点初始化为一个社群;
步骤三:计算每个社群S的分区性能;
步骤四:计算节点i从原本的社群中删除、加入相邻社群后,i原本所在社群及J相邻社群的分区性能的增益;若增益大于0,则将i从原本的社群中删除、加入能使增益最大的相邻社群;若不存在,则跳至网络中下一个节点,直至对网络中所有节点所属社群进行一轮调整;
步骤五:重复步骤三、四,直至所有节点的社群不再发生变化;
步骤六:对于某个社群Ci,计算Ci的含有B条边的所有子图即为所有可能的候选投资路段的组合;例如,中国安徽省的全省路网中,在合肥市附近区域形成一个社群,其中合肥绕城高速的金寨路、包河大道,合六叶路段以及芜合路段构成该社区的一个子图;
步骤七:计算每种组合中,度大于1的节点的个数
步骤八:根据下式,计算每种组合发生的可能性
步骤九:假设所有路段的固有事故率均为Ppr(一般取0.15);
步骤十:对于没有被投资的路段,其事故率即为:Pe=Ppr;
步骤十一:对于被投资的路段,其事故率为固有事故率Ppr与投资收益率Pinv之差,即Pe=Ppr-Pinv(Pinv一般取0.1);
步骤十二:按事故率对Ci中的每条路段是否发生拥堵或路面损坏等事故进行抽样,将抽样结果中发生事故的路段从Ci中去除,得到一次抽样后的仿真社群拓扑结构例如,合肥市附近所形成的社群进行一次抽样后,合六叶路段、蚌合路段被抽样为发生事故,则从社群中去掉这两条边;
步骤十三:为了消除抽样的随机性,重复S次步骤十二,即得
步骤十四:对于经过路网的每辆车ho,提取其OD信息及经过的所有路段序列zo;
计算每辆车ho经过当前仿真社群拓扑结构时的通行效率
步骤十五:若路段序列zo中的所有边均存在于通行效率为车辆ho通过路径长度与历史数据中的实际通行时间之比;
步骤十六:若路段序列zo中,存在某些边不属于仿真社群拓扑结构则利用迪杰斯特拉算法(Dijkstra's Algorithm)重新在拓扑结构上计算车辆ho从起点到终点的仿真路段序列z′o;而用重新计算后的z′o来计算车辆ho的通行效率:
步骤十七:依据下式计算社群Ci内部的B条候选关键路段集合
社群Ci内部的B条候选关键路段,即为上式所求结果中所包含的B条边;
步骤十八:对于划分出来的每个社群,进行如上步骤六到步骤十七的操作,得到每个社群内部的B条候选关键路段;例如,合肥市附近区域所对应的社群中,候选关键路段为:金寨路、包河大道、合六叶及蚌合路段;芜湖市附近区域所对应的社群中,候选关键路段为:芜宣路段、芜合路段、马巢路段、马芜路段;中国安徽省其他社群同样对应着4条候选关键路段;
步骤十九:记录当步骤十七中的式子取最大值时,所对应B条候选关键路段各自的通行效率,记为v(i);
步骤二十:记l*B条候选关键路段的投资花费为w(i),为简化模型可使w(i)均为1;
步骤二十一:申请一个l*B行W列的二维数组Y,将元素初始化为-1;
步骤二十二:依据递推公式Y[i][j]=max(f[i+1][j],f[i+1][j-w[i]]+v[i]),递归的计算Y[0][W]的值,并且记录每次被背包选中的候选关键路段;当递归到某一层时,当前的Y[i][j]不等于-1,表示之前已经计算过,直接返回计算过的值;否则,先判断第i+1条候选路段的花费是否超过当前剩余花费,若不超过,则按照递推公式计算当前的Y[i][j];若超过,则跳过第i+1条候选路段,计算Y[i+1][j];
步骤二十三:直到递归函数完全结束,所记录的B条被背包选中的候选关键路段即为所求的全路网关键路段;例如,最终挖掘出的合肥绕城高速金寨路段、马芜路段、蚌合路段以及滁宁段,为安徽全省路网结构中的关键路段。
上述具体过程运用到中国安徽省高速路网数据,使用一种高效准确的高速路网中的关键路段识别方法找出影响道路通行效率的关键路段。通过本发明提供的技术方案,运维人员可将主要人力部署在关键路段上,一旦发生交通事故或拥堵,可快速纾解交通压力,减少经济和时间损失;并以此作为路网维护管理的合理依据。
需要注意的是,公布实施例的目的在于帮助进一步理解本发明,但是本领域的技术人员可以理解:在不脱离本发明及所附权利要求的精神和范围内,各种替换和修改都是可能的。因此,本发明不应局限于实施例所公开的内容,本发明要求保护的范围以权利要求书界定的范围为准。
Claims (7)
1.一种高速路网中关键路段的识别方法,所述关键路段指当高速路网发生交通拥堵或道路损坏时,会对路网通行效率造成严重影响的路段,其中路网节点是关键路段的一部分;以高速路网中收费站和路段的地理位置及其连接关系构成当前高速路网拓扑结构G;以每辆车的OD信息和每段途经路段的通行时间为通行车辆历史数据,OD表示一段行程的起点origin和终点destination;
所述方法首先利用社群划分算法将路网划分成K个社群网络;再通过概率假设方法对路网建模,在每个社群网络中,各求解得到B条候选关键路段;再从所述候选关键路段中,采用动态规划方法找到最终的关键路段,作为高速路网中的关键路段,使得投资所述关键路段的组合能够实现路网通行效率最高;包括如下步骤:
1)首先提取高速路网中所有N个收费站{v1,...,vN},形成路网拓扑结构G中的节点集合V;提取高速路网中所有M条路段{e1,...,eM},生成路网拓扑结构G=(V,E),计算得到每条边ei上的属性信息
2)对路网拓扑结构G进行社群划分:将每个节点初始化为一个社群,再利用社群分区性能增益进行调整,由此得到一组社群的集合C={C1,...,Cl};l为社群的个数;
3)在每个社群内部识别候选关键路段得到每个社群内部的B条候选关键路段和所述B条候选关键路段上每条路段上的通行效率,
4)利用动态规划的方法,挖掘出至多B条能够使得全路网道路通行效率达到最高的关键路段;
5)将这B条路段输出,作为最终的高速路网中的关键路段。
2.如权利要求1所述识别方法,其特征是,步骤1)计算得到边上的属性信息具体为:
所述边相应路段的长度为所述通行车辆历史数据中车辆经过该条路段的平均通行时间为所述通行车辆历史数据中该条路段经过的车流量为由此得到每条边ei上的属性信息
3.如权利要求1所述识别方法,其特征是,步骤2)对路网拓扑结构G进行社群划分,具体包括如下步骤:
21)每个节点初始化为一个社群;
22)通过式1计算得到每个社群S的分区性能f(S):
其中,cs表示一个顶点在S中、一个顶点不在S中的边的个数;ms表示完全包含在S中的边的个数;
23)对于每个节点i,计算每个与i相邻的社群Ji;
24)对于i的每个相邻社群Ji,计算得到假设i从原本的社群中删除、加入Ji中后,i原本所在社群及Ji的分区性能的增益;
25)根据步骤24)得到的社群分区性能增益对每个节点所述社群进行调整:若i的所有相邻社群中,存在某些相邻社群在接纳i之后,使得步骤24)所述的社群分区性能增益大于0,则将i从原本的社群中删除、加入能使增益最大的相邻社群;若不存在,则本轮对节点i的调整结束,跳至网络中下一个节点,直至对网络中所有节点进行一轮调整;
26)重复步骤22)至25),直至所有节点的社群不再发生变化,从而得到一组社群的集合C={C1,...,Cl}。
4.如权利要求1所述识别方法,其特征是,步骤3)所述在每个社群内部识别候选关键路段的方法,具体包括如下步骤:
31)对于社群Ci,计算得到Ci含有B条边的所有子图表示所有可能的候选投资路段的组合;
32)计算得到每种组合中,度大于1的节点的个数
33)通过式2计算得到每种组合发生的可能性
其中,|Ei|表示社群Ci中的边数;表示从|Ei|条边中选出B条边的排列数;
34)设定所有路段的固有事故率相同,分别得到没有被投资的路段和被投资的路段的事故率Pe;
35)按事故率对Ci中的每条路段进行抽样,将抽样结果中发生事故的路段从Ci中去除,得到抽样后的仿真社群拓扑结构进行S次抽样,即得
36)设历史数据中有H辆车经过路网,对于其中每辆车ho,提取其OD信息及经过的所有路段序列zo;
37)计算得到每辆车ho经过当前仿真社群拓扑结构时的通行效率
38)社群Ci内部的B条候选关键路段集合为式3:
社群Ci内部的B条候选关键路段,即为上式所求结果中所包含的B条边;
39)对每个社群执行步骤31)~38)的操作,得到每个社群内部的B条候选关键路段。
5.如权利要求4所述识别方法,其特征是,步骤34)具体设定所有路段的固有事故率均为Ppr,对于没有被投资的路段,其事故率就是固有事故率;对于被投资的路段,其事故率为固有事故率Ppr与投资收益率Pinv之差。
6.如权利要求4所述识别方法,其特征是,步骤37)计算得到每辆车ho经过当前仿真社群拓扑结构时的通行效率具体包括步骤371)~372):
371)若路段序列zo中的所有边均存在于则车辆ho的通行时间即等于历史数据中的实际通行时间;通行效率为车辆ho通过路径长度与通行时间之比;
372)若路段序列zo中,存在某些边不属于仿真社群拓扑结构则利用迪杰斯特拉算法重新在拓扑结构上计算车辆ho从起点到终点的仿真路段序列z′o;则车辆ho的通行效率通过式4计算得到:
其中,eq表示仿真路段序列z′o中的第q条路段;表示第q条路段的长度;表示历史数据中车辆经过第q条路段的平均通行时间。
7.如权利要求1所述识别方法,其特征是,通过利用所述动态规划方法,在全路网中挖掘得到至多B条最终的关键路段,具体包括如下步骤:
41)从1个社群中的l*B条候选关键路段中挖掘得到至多B条最终的关键路段;记所述l*B条候选关键路段的投资花费为w(i);每条关键路段上的通行效率记为v(i);
42)使用动态规划中的背包问题求解方法,首先初始化一个l*B行W列的二维数组Y,将元素初始化为一固定值;
43)则解决所述在全路网中得到至多B条最终关键路段的问题相应的递推公式为式5:
Y[i][j]=max(f[i+1][j],f[i+1][j-w[i]]+v[i]) (式5)
44)以二维数组Y的两个维度分别作为递归函数的两个参数i、j;参数i、j的含义为:从第i条候选路段开始,选择能使总通行效率增益最大的、且总花费不超过j的候选路段;
45)首先从Y[0][W]开始递归;
46)进入递归函数后,若当前二维数组中对应的元素Y[i][j]不等于初始化值,表示之前已经计算过;直接返回计算过的值;
47)若当前二维数组中对应的元素Y[i][j]等于初始化值,即没有计算过,此时判断第i+1条候选路段的投资花费是否超过当前剩余投资花费,若不超过,则按照式5计算;若超过,则跳过第i+1条候选路段,计算Y[i+1][j];
48)直到递归函数完全结束,被选中的B条候选关键路段即为最终所求的全路网关键路段。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710301222.0A CN106981194B (zh) | 2017-05-02 | 2017-05-02 | 一种高速路网关键路段的识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710301222.0A CN106981194B (zh) | 2017-05-02 | 2017-05-02 | 一种高速路网关键路段的识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106981194A CN106981194A (zh) | 2017-07-25 |
CN106981194B true CN106981194B (zh) | 2019-06-07 |
Family
ID=59341839
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710301222.0A Expired - Fee Related CN106981194B (zh) | 2017-05-02 | 2017-05-02 | 一种高速路网关键路段的识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106981194B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108765958A (zh) * | 2018-06-25 | 2018-11-06 | 交通运输部科学研究院 | 高速公路交通流量分配方法及装置 |
CN109559511B (zh) * | 2018-12-05 | 2019-10-29 | 江苏省城市规划设计研究院 | 一种城市交通拥堵信息定向投放方法 |
CN111654824B (zh) * | 2020-05-29 | 2023-07-25 | Oppo广东移动通信有限公司 | 网络状况提示方法、装置、电子设备及存储介质 |
CN111798663B (zh) * | 2020-08-07 | 2021-12-14 | 交通运输部公路科学研究所 | 一种路网关键路段辨识方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104217579A (zh) * | 2013-07-23 | 2014-12-17 | 同济大学 | 基于路段子网络冗余度的交通网络关键路段搜索方法 |
CN105702031A (zh) * | 2016-03-08 | 2016-06-22 | 北京航空航天大学 | 基于宏观基本图的路网关键路段识别方法 |
CN105894802A (zh) * | 2015-04-10 | 2016-08-24 | 杭州远眺科技有限公司 | 一种基于gps数据的交通拥堵传播路径计算方法 |
CN106157083A (zh) * | 2015-04-14 | 2016-11-23 | 阿里巴巴集团控股有限公司 | 挖掘潜在客户的方法和装置 |
-
2017
- 2017-05-02 CN CN201710301222.0A patent/CN106981194B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104217579A (zh) * | 2013-07-23 | 2014-12-17 | 同济大学 | 基于路段子网络冗余度的交通网络关键路段搜索方法 |
CN105894802A (zh) * | 2015-04-10 | 2016-08-24 | 杭州远眺科技有限公司 | 一种基于gps数据的交通拥堵传播路径计算方法 |
CN106157083A (zh) * | 2015-04-14 | 2016-11-23 | 阿里巴巴集团控股有限公司 | 挖掘潜在客户的方法和装置 |
CN105702031A (zh) * | 2016-03-08 | 2016-06-22 | 北京航空航天大学 | 基于宏观基本图的路网关键路段识别方法 |
Non-Patent Citations (2)
Title |
---|
基于路网可靠性的关键点段辨识;左志等;《交通运输研究》;20151031;第1卷(第5期);第18-24页 |
面向大规模社会网络TOP-K影响力节点挖掘的并行算法;周夏冰等;《计算机研究与发展》;20141231;第216-217页 |
Also Published As
Publication number | Publication date |
---|---|
CN106981194A (zh) | 2017-07-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106981194B (zh) | 一种高速路网关键路段的识别方法 | |
CN101650191B (zh) | 道路网络拓扑抽象的方法及装置 | |
CN105258704A (zh) | 基于快速路网建模的多尺度时空热点路径探测方法 | |
CN107103754A (zh) | 一种道路交通状况预测方法及*** | |
CN103593430A (zh) | 一种基于移动对象时空信息轨迹分段聚类的方法 | |
CN104462190A (zh) | 一种基于海量空间轨迹挖掘的在线的位置预测方法 | |
CN108549952A (zh) | 车辆搭载无人机双层路径的优化方法及装置 | |
CN109359350A (zh) | 一种优化精细施工成本的风电场道路智能设计方法 | |
CN108280463A (zh) | 车辆搭载无人机双层路径的优化方法及装置 | |
CN111507879A (zh) | 一种城市空间结构识别方法 | |
CN109859505A (zh) | 高速站点的预警处理方法、装置、服务器和介质 | |
CN115586557A (zh) | 一种基于路网数据的车辆行驶轨迹纠偏方法及装置 | |
CN105205247A (zh) | 一种基于树结构的仿真路网数据管理方法 | |
CN107071721A (zh) | 基于手机定位数据的停留点提取方法 | |
CN109520499A (zh) | 基于车辆gps轨迹数据实现区域实时等时线的方法 | |
CN105260551B (zh) | 车联网中节点分布特性分析方法 | |
CN109214656A (zh) | 一种用于城市公交网络的节点重要性评估方法 | |
CN116307325A (zh) | 一种配电网的线路规划方法、装置、电子设备及存储介质 | |
CN116772846A (zh) | 一种无人机航迹规划方法、装置、设备及介质 | |
CN107016048A (zh) | 基于追赶法求解圣维南方程组的树状河网智能搜索方法 | |
CN110530391A (zh) | 一种基于分层路网的路径规划方法及装置 | |
CN114637305B (zh) | 一种无人机最短路径规划方法及装置 | |
CN115686070A (zh) | 配电/输电线路无人机自主巡线路径规划方法 | |
CN105427581B (zh) | 基于浮动车数据的交通仿真方法及*** | |
CN102542617A (zh) | 一种复杂路网的路径识别标识站设置方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190607 |
|
CF01 | Termination of patent right due to non-payment of annual fee |