CN113205676A - 一种基于元胞传输模型的建模方法及参数标定方法 - Google Patents

一种基于元胞传输模型的建模方法及参数标定方法 Download PDF

Info

Publication number
CN113205676A
CN113205676A CN202110337985.7A CN202110337985A CN113205676A CN 113205676 A CN113205676 A CN 113205676A CN 202110337985 A CN202110337985 A CN 202110337985A CN 113205676 A CN113205676 A CN 113205676A
Authority
CN
China
Prior art keywords
ramp
traffic
cell
flow
cells
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
CN202110337985.7A
Other languages
English (en)
Other versions
CN113205676B (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.)
Zhejiang Gaoxin Technology Co Ltd
Original Assignee
Zhejiang Highway Information Engineering Technology Co ltd
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 Zhejiang Highway Information Engineering Technology Co ltd filed Critical Zhejiang Highway Information Engineering Technology Co ltd
Priority to CN202110337985.7A priority Critical patent/CN113205676B/zh
Publication of CN113205676A publication Critical patent/CN113205676A/zh
Application granted granted Critical
Publication of CN113205676B publication Critical patent/CN113205676B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/01Detecting movement of traffic to be counted or controlled
    • G08G1/0104Measuring and analyzing of parameters relative to traffic conditions
    • G08G1/0125Traffic data processing
    • G08G1/0129Traffic data processing for creating historical data or processing based on historical data
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/01Detecting movement of traffic to be counted or controlled
    • G08G1/0104Measuring and analyzing of parameters relative to traffic conditions
    • G08G1/0137Measuring and analyzing of parameters relative to traffic conditions for specific applications
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/065Traffic control systems for road vehicles by counting the vehicles in a section of the road or in a parking area, i.e. comparing incoming count with outgoing count

Abstract

本发明公开一种基于元胞传输模型的建模方法及参数标定方法,建模方法首先将元胞划分为三类,然后根据流量守恒理论建立高速路网元胞传输理论模型,仿真预测方法然后再用历史数据对元胞传输理论模型的基本参数进行标定,并估算因事件引起的道路剩余通行能力,然后再对匝道分流区、合流区的流量比例进行估算,对匝道管控方案实施后交通需求进行预测;最后基于标定的所有参数,运行S2的高速路网元胞传输理论模型,对路网交通状态的演变进行仿真,输出仿真结果。本发明的方法可作为匝道管控方案的制定与选择的基础依据。

Description

一种基于元胞传输模型的建模方法及参数标定方法
技术领域
本发明涉及交通仿真技术领域,尤其涉及一种基于元胞传输模型的建模方法及参数标定方法。
背景技术
高速公路作为公路网的重要组成部分,对交通运输经济的发展起着重要的作用。当高速某路段发生一起或多起交通事件时,高速公路管理部门需采取多种交通治理措施以缓解事件造成的拥堵,其中,对事件发生上游匝道流量进行管控,以减少进入事件发生路段的车流量,被认为是一种行之有效的方法。理论上来说,当管控后的上游流量小于或等于事件路段的通行能力时,可达到抑制拥堵扩散的效果;但当管控流量过小时,可能导致收费站车辆排队溢出至城市道路,进而导致进入收费站的城市道路的拥堵,则会造成进入收费站的城市道路的拥堵。然而,在实际应用中,制定合理的匝道管控方案难度较大,其主要原因是缺乏对匝道管控方案下发后的预测与评估方法。
元胞传输模型是以流体动力学为基础的LWR宏观交通流模型的离散化近似,其基本思想为:将路段划分为一系列首尾相连的元胞,根据流量守恒定律、最大容量约束与元胞的当前状态计算元胞的输入和输出交通流量,由此确定元胞的下一个状态,如此不断迭代,便可获得元胞的动态交通特征。高速公路场景中,可以采用基于元胞传输模型的匝道管控效果仿真预测方法,对交通事件时的事件路段、及上游上下闸道车流进行仿真预测,从而对管理部门开展上游流量疏导提供理论依据。因此,仿真就需要首先构建元胞传输模型的建模方法,并考虑元胞传输模型在实际使用中需考虑匝道分合流区的问题,同时也需要对模型的参数进行标定。
发明内容
本发明的目的在于针对现有技术的不足,提供一种基于元胞传输模型的建模方法及参数标定方法,建模方法对现有元胞传输模型进行改进,设计针对匝道合流区的元胞划分和建模方法,元胞传输模型的参数标定方法提出一种基于最小二乘法的交通宏观基本图参数标定方法,解决了元胞传输模型仿真在实际应用中的参数标定问题。
本发明的目的通过如下的技术方案来实现:
一种基于元胞传输模型的建模方法,其特征在于,该方法具体包括如下步骤:
S1:以事件发生地点为终止位置,起始位置为事件发生地点上游,起始终止位置之间路段包括距离事件发生地点为终止位置最近的上匝道、下匝道,对上述高速路网路段进行元胞划分,根据有无匝道,分为无匝道元胞,有上匝道元胞和有下匝道元胞三类;元胞划分需满足的原则:(1)同一个元胞内最多只有一个上匝道或下匝道;(2)所述的有上匝道元胞的起始位置位于上匝道入口的上游,所述的有下匝道元胞的终止位置位于下匝道出口位置的下游;
S2:基于流量守恒理论,建立S1得到的三类元胞的高速路网元胞传输理论模型。
进一步的,所述的S2具体通过如下的子步骤来实现:
S2.1:基于流量守恒理论,元胞i内车密度随时间变化的更新公式如下:
Figure BDA0002998311210000021
其中:
ρi(k)表示第k个时间步内,元胞i内的车密度;
Figure BDA0002998311210000022
分别表示第k个时间步内高速路主线元胞i的输入和输出车流量;
ΔT表示元胞更新时间步的时长;
Li表示元胞i的长度;
S2.2:根据S1中的三类元胞,
Figure BDA0002998311210000023
的计算也分为三种情况:
(1)对于无匝道元胞,元胞内仅存在主线车流
Figure BDA0002998311210000024
Si-1(k)=min(ρi-1(k)vf,i-1(k),qm,i-1) (3)
Ri(k)=min(qm,i,wiJi(k))) (4)
其中,Si-1(k)为元胞i-1在自由流状态下的最大发送流率,Ri(k)为允许进入元胞i的最大接收流量,vf,i-1表示元胞i-1处于自由流状态下的平均车速;qm,i、qm,i-1分别表示元胞i、i-1的最大通行流量;wi表示元胞i处于拥堵状态下的反向激波速度;ρJ表示车流完全堵塞下的车密度;
(2)针对有上匝道元胞,元胞内存在上匝道和主线两股车流的合流区
Figure BDA0002998311210000025
Figure BDA0002998311210000026
Figure BDA0002998311210000027
其中,ri(k)为元胞i内由匝道i进入主线元胞i的流量,di(k)为与元胞i相连的上匝道i的交通需求;αi(k)为上匝道车流占合流区总车流的比例;
Figure BDA0002998311210000031
表示由主线元胞i-1进入元胞i的车流量;
(3)针对有下匝道元胞,元胞内存在下匝道和主线两股车流的分流区
Figure BDA0002998311210000032
Figure BDA0002998311210000033
其中,βi(k)为元胞i内分流区进入下匝道的车流比例,
Figure BDA0002998311210000034
为输入元胞i+1的流量。
进一步的,通过下式对所述的对匝道管控方案实施后交通需求di(k)进行预测:
di(k)=min(Di(k),ql,c×l) (17)
其中,Di(k)为前1小时内收费站的车流量,l为实施管控后收费站剩余车道数,ql,c为收费站单车道的最大通行能力。
本发明第二个目的在于一种元胞传输模型的参数标定方法,在基于元胞传输模型的建模方法的步骤S1、S2之上,还包括:提供S3:基于高速主线上的ETC门架和匝道的卡口检测器获取的历史交通数据,采用基于最小二乘法的元胞传输模型参数标定方法对高速路网元胞传输理论模型的基本参数进行标定;并估算因事件引起的道路剩余通行能力;
S4:对匝道分流区、合流区的流量比例进行估算,并对匝道管控方案实施后的上匝道交通需求进行预测。
进一步的,所述的S3中,元胞中的基本参数qm、vf、w、ρJ的标定的过程如下:
(1)根据高速路段交通宏观基本图的理论模型,车流量q和车密度ρ的相关关系用下式表示:
q=min{ρvf,qm,w(ρJ-ρ)} (10)
(2)基于ETC门架和匝道的卡口检测器获取的历史交通数据,假设数据离散化后,时间步k=0,1,2…,则元胞i的历史车流量数据q′i、平均车速v′i和车密度ρ′i分别表示为:
q′i=[q′i(0),q′i(1),…,q′i(k)]T (11)
v′i=[v′i(0),v′i(1),…,v′i(k)]T (12)
ρ′i=[ρ′i(0),ρ′i(1),…,ρ′i(k)]T=[q′i(0)/v′i(0),q′i(1)/v′i(1),…,q′i(k)/v′i(k)]T (13)
(3)利用最小二乘法对公式(10)进行标定,标定方法用下式表示:
Figure BDA0002998311210000035
Figure BDA0002998311210000041
vf≤vmax,qm≤qmax,ρJ≤ρmax (16)
其中,
Figure BDA0002998311210000042
为qi(k)基于历史交通数据的的估计值;
(4)运用非线性优化问题求解器求解公式(14)~(16),获得元胞的基本参数qm、vf、w、ρJ
当道路发生拥堵时,用剩余通行能力q′m代替qm,q′m的标定的方法如下:
①对于影响车道数n已知的事件,则按照
Figure BDA0002998311210000043
计算,N为路段总车道数;
②对于影响车道数未知的事件,q′m取事件发生路段下游检测器检测到的车流量值。
进一步的,所述的S4中,
对于有上匝道元胞,合流区上匝道流量比例αi(k)的标定方式为,将上匝道卡口检测器和主线ETC门架的车流数据基于车牌信息进行匹配,并计算合流区的前1小时内上匝道的车流量占比;
对于有下匝道元胞,下匝道分流比例βi(k)的标定方法为:将主线ETC门架和下匝道卡口检测器的车流数据基于车牌信息进行匹配,计算分流区的前1小时内下匝道的车流量占比。
本发明的有益效果如下:
本发明的基于元胞传输模型的建模方法及元胞传输模型的参数标定方法,该方法对现有元胞传输模型进行改进,设计针对匝道合流区的元胞划分和建模方法,提出一种基于最小二乘法的交通宏观基本图参数标定方法,解决了元胞传输模型仿真在实际应用中的参数标定问题,并将上匝道(进口匝道)、收费站车道数作为管控上游流量出口的工具,并作为交通需求进行估算,从而对匝道管控方案流量模拟仿真提供了基础。
附图说明
图1是本发明中基于元胞传输模型的建模方法及参数标定方法工作流程图;
图2是高速路段元胞划分示例图;
图3是宏观交通基本图示例。
具体实施方式
下面根据附图和优选实施例详细描述本发明,通过对实施例的描述,本发明的目的和效果将变得更加明白,应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
首先,本发明的基于元胞传输模型,构建一种以事件为终止位置、包括匝道的建模方法,对现有元胞传输模型进行改进,设计针对匝道合流区的元胞划分和建模方法;
其次,提出的交通宏观基本图理论模型参数标定方法、剩余通行能力估算方法、匝道分合流区流量比例估算方法,解决了元胞传输模型仿真在实际应用中的参数标定问题。
如图1~3所示,本发明的基于元胞传输模型的建模方法及参数标定方法,建模方法为步骤S1、S2,参数标定方法还包括步骤S3、S4,参数标定方法是在建模方法上的进一步延申,下面将两种方法合并描述,包括如下步骤:
S1:以事件发生地点为终止位置,起始位置为事件发生地点上游,起始终止位置之间路段包括距离事件发生地点为终止位置最近的一个上匝道、一个下匝道,对上述高速路网路段进行元胞划分,根据有无匝道,分为无匝道元胞,有上匝道(进口匝道)元胞和有下匝道(出口匝道)元胞三类;元胞划分需满足的原则:(1)同一个元胞内最多只有一个上匝道或下匝道;(2)所述的有上匝道元胞的起始位置位于上匝道入口的上游,所述的有下匝道元胞的终止位置位于下匝道出口位置的下游。如图2所示,是高速路网中的一种典型场景,该场景中有一对邻近的下匝道和上匝道,常见于高速公路服务区、收费站附近。针对该场景,高速主线被划分为5个元胞,即元胞i=0,1,2,3,4。高速路网经元胞划分后,可分为如下三类:
①元胞内无匝道(即图2中的元胞0、2和4);
②元胞内有一个上匝道(即图2中的元胞3);
③元胞内有一个下匝道(即图2中的元胞1)。
S2:基于流量守恒理论,建立S1得到的三类元胞的高速路网元胞传输理论模型;
S2.1:基于流量守恒理论,元胞i内车密度随时间变化的更新公式如下:
Figure BDA0002998311210000051
其中:
ρi(k)表示第k个时间步内,元胞i内的车密度;
Figure BDA0002998311210000052
分别表示第k个时间步内高速路主线元胞i的输入和输出车流量;
ΔT表示元胞更新时间步的时长;
Li表示元胞i的长度;
S2.2:根据S1中的三类元胞,
Figure BDA0002998311210000053
的计算也分为三种情况:
(1)对于无匝道元胞,元胞内仅存在主线车流
Figure BDA0002998311210000054
Si-1(k)=min(ρi-1(k)vf,i-1(k),qm,i-1) (3)
Ri(k)=min(qm,i,wiJi(k))) (4)
其中,Si-1(k)为元胞i–1在自由流状态下的最大发送流率,Ri(k)为允许进入元胞i的拥堵状态下最大接收流量的最小值,最小值是指min函数比较最大流率与反向激波Wi乘以密度两者中的最小值,其依据是图3中的宏观交通基本图的原理;vf,i-1表示元胞i-1处于自由流状态下的平均车速;qm,i、qm,i-1分别表示元胞i、i-1的最大通行流量;wi表示元胞i处于拥堵状态下的反向激波速度;ρJ表示车流完全堵塞下的车密度;
公式2-4,元胞i的流量q受两个因素限制:上游i-1输出流量(send)的最大值和元胞i能够接收(receive)的最大流量值。元胞i实际
Figure BDA0002998311210000061
流量输入为元胞i-1实际自由流量输出,但受元胞i-1最大流量、以及元胞i本身最大流量、拥堵导致的方向激波流量限制。
(2)针对有上匝道元胞,元胞内存在上匝道和主线两股车流的合流区,此处存在两种情况:1)元胞i可同时接收来自上匝道i和主线上游元胞i-1的交通需求;2)上匝道i和主线上游元胞i-1的交通需求大于元胞i可接收的最大流量。对于这两种情况,由主线元胞i-1流入元胞i的流量
Figure BDA0002998311210000062
和由匝道i进入主线元胞i的流量ri(k)可分别用下式表示:
Figure BDA0002998311210000063
Figure BDA0002998311210000064
Figure BDA0002998311210000065
其中,ri(k)为元胞i内由匝道i进入主线元胞i的流量,di(k)为与元胞i相连的上匝道i的交通需求;αi(k)为上匝道车流占合流区总车流的比例;
Figure BDA0002998311210000066
表示由主线元胞i-1进入元胞i的车流量。公式(5)和(6)表示的是有上匝道(进口匝道)情形下元胞i的流量计算方式。其中,公式(5)中表示的是有高速路元胞i-1输入元胞i的流量qi,i-1,公式(6)表示的是由匝道输入元胞i的车流ri。然后,公式(7)对两股流量进行加和。
公式(5)(6)中的判断条件是判断元胞i能否接受分别来自高速上游元胞i-1和匝道ri的两股车流(即Ri是否大于Si-1和进口匝道需求di)。当;Si-1(k)+di(k)≤Ri(k),
Figure BDA0002998311210000067
为元胞i-1在自由流状态下的最大发送流率,ri(k)为元胞i相连的上匝道i的交通需求。当Si-1(k)+di(k)>Ri(k),按照αi(k)为上匝道、主线车道分配车流比例。
(3)针对有下匝道元胞,元胞内存在下匝道和主线两股车流的分流区,如图2中的相邻元胞1和2,元胞1内存在分流区。假设元胞i内分流区进入下匝道的车流比例为βi(k),则元胞i的输出流量
Figure BDA0002998311210000071
和输入元胞i+1的流量
Figure BDA0002998311210000072
可分别由下式计算:
Figure BDA0002998311210000073
Figure BDA0002998311210000074
其中,βi(k)为元胞i内分流区进入下匝道的车流比例,
Figure BDA0002998311210000075
为输入元胞i+1的流量;
公式5、6,主线和匝道流量大于元胞i的Ri(k)情况下,采用ai(k)分配,公式8、9元胞i内分流区进入下匝道的车流比例为βi(k),这两个比例由步骤S3标定。Ri+1表示的是高速路主线下游元胞i+1的接收流量,但对于下匝道(出口匝道)的情形,假设元胞i内分流区进入下匝道的车流比例为βi,进入元胞i+1的比例为1-βi,那么下游元胞i+1和下匝道的总接收流量为Ri+1/(1-βi)。
S3:基于高速主线上的ETC门架和匝道的卡口检测器获取的历史交通数据,即元胞起始和终点路段之间所有的门架、匝道,如果有多个门架,可对多个ETC门架获得的流量数据取平均值,基于均值进行参数标定。采用基于最小二乘法的元胞传输模型参数标定方法对高速路网元胞传输理论模型进行标定;并估算因事件引起的道路剩余通行能力;
其中,元胞中的基本参数qm、vf、w、ρJ的标定的过程如下:
(1)根据如图3所示的高速路段交通宏观基本图以及其理论模型,交通宏观基本图定义了路段交通流的基本性质,即当路段车密度小于饱和密度时,路段车流量与车密度成正相关;当车密度大于饱和密度时,路段开始出现拥堵,车流量随车密度增大而减小,车流量随车密度的变化率被定义为反向激波速度,车流量q和车密度ρ的相关关系用下式表示:
q=min{ρvf,qm,w(ρJ-ρ)} (10)
(2)基于ETC门架和匝道的卡口检测器获取的历史交通数据,假设数据离散化后,时间步k=0,1,2…,则元胞i的历史车流量数据q′i、平均车速v′i和车密度ρ′i分别表示为:
q′i=[q′i(0),q′i(1),…,q′i(k)]T (11)
v′i=[v′i(0),v′i(1),…,v′i(k)]T (12)
ρ′i=[ρ′i(0),ρ′i(1),…,ρ′i(k)]T=[q′i(0)/v′i(0),q′i(1)/v′i(1),…,q′i(k)/v′i(k)]T (13)
(3)利用最小二乘法对公式(10)进行标定,标定方法用下式表示:
Figure BDA0002998311210000076
Figure BDA0002998311210000081
vf≤vmax,qm≤qmax,ρJ≤ρmax (16)
其中,
Figure BDA0002998311210000082
为qi(k)基于历史交通数据的的估计值;
公式(14)-(16)是一个非线性优化问题。其中,公式(14)表示最小化由宏观交通图计算得到的车流量qi(k)的估计值
Figure BDA0002998311210000085
和检测器测得的qi(k)的差异的平方和。公式(15)是约束条件1,表示
Figure BDA0002998311210000083
是公式(13)基于历史检测数据的计算结果;公式(16)是约束条件2,分别表示自由流速度vf小于或等于路段限速(一般为120km/h),路段通行能力qm小于或等于理论最大通行能力qmax,路段堵塞车密度ρJ小于等于道路理论最大车密度
(4)运用非线性优化问题求解器(如LINGO)求解公式(14)~(16),获得元胞的基本参数qm、vf、w、ρJ
当道路发生拥堵时,用剩余通行能力q′m代替qm,q′m的标定的方法如下:
①对于影响车道数n已知的事件,则按照
Figure BDA0002998311210000084
计算,N为路段总车道数;
②对于影响车道数未知的事件,q′m取事件发生路段下游检测器检测到的车流量值。这个车流量为统计时间段内(比如15分钟)的近似值。
S4:对匝道分流区、合流区的流量比例进行估算,并对匝道管控方案实施后交通需求进行预测;
对于有上匝道元胞,合流区上匝道流量比例αi(k)的标定方式为,将上匝道卡口检测器和主线ETC门架的车流数据基于车牌信息进行匹配,并计算合流区的前1小时内上匝道的车流量占比;
对于有下匝道元胞,下匝道分流比例βi(k)的标定方法为:将主线ETC门架和下匝道卡口检测器的车流数据基于车牌信息进行匹配,计算分流区的前1小时内下匝道的车流量占比;
通过下式对所述的对匝道管控方案实施后交通需求di(k)进行预测:
di(k)=min(Di(k),ql,c×l) (17)
其中,Di(k)为前1小时内收费站的车流量,l为实施管控后收费站剩余车道数,ql,c为收费站单车道的最大通行能力。
S5:基于S3和S4标定的所有参数,运行S2的高速路网元胞传输理论模型,对路网交通状态的演变进行仿真,输出仿真结果,具体如下:
(1)初始化车密度ρi(0)、时间步k=0、元胞i=0;
(2)基于S3和S4标定的所有参数,当k≤K,i≤I,根据S2中的公式(1)~(9),遍历计算每个元胞第k+1时间步的车密度;根据公式
Figure BDA0002998311210000091
计算第k步的元胞的平均车速;其中,K为总的预测时间步数,I为路网中最后一个元胞的连续编号;
(3)返回所有元胞的仿真结果[ρi(0),…,ρi(k)],[vi(0),…,vi(k)];
(4)根据(3)得到的所有元胞的车密度和车速数据,评价主线路段的交通状态,即匝道管控方案实施后的效果。
本领域普通技术人员可以理解,以上所述仅为发明的优选实例而已,并不用于限制发明,尽管参照前述实例对发明进行了详细的说明,对于本领域的技术人员来说,其依然可以对前述各实例记载的技术方案进行修改,或者对其中部分技术特征进行等同替换。凡在发明的精神和原则之内,所做的修改、等同替换等均应包含在发明的保护范围之内。

Claims (9)

1.一种基于元胞传输模型的建模方法,其特征在于,该方法具体包括如下步骤:
S1:以事件发生地点为终止位置,起始位置为事件发生地点上游,起始终止位置之间路段包括距离事件发生地点为终止位置最近的上匝道、下匝道,对上述高速路网路段进行元胞划分,根据有无匝道,分为无匝道元胞,有上匝道元胞和有下匝道元胞三类;元胞划分需满足的原则:(1)同一个元胞内最多只有一个上匝道或下匝道;(2)所述的有上匝道元胞的起始位置位于上匝道入口的上游,所述的有下匝道元胞的终止位置位于下匝道出口位置的下游;
S2:基于流量守恒理论,建立S1得到的三类元胞的高速路网元胞传输理论模型。
2.根据权利要求1所述的基于元胞传输模型的建模方法,其特征在于,
所述的S2具体通过如下的子步骤来实现:
基于流量守恒理论,元胞i内车密度随时间变化的更新公式如下:
Figure FDA0002998311200000011
其中:
ρi(k)表示第k个时间步内,元胞i内的车密度;
Figure FDA0002998311200000012
分别表示第k个时间步内高速路主线元胞i的输入和输出车流量;
ΔT表示元胞更新时间步的时长;
Li表示元胞i的长度。
3.根据权利要求1所述的基于元胞传输模型的建模方法,其特征在于,
根据S1中的三类元胞中无匝道元胞,元胞内仅存在主线车流
Figure FDA0002998311200000013
Si-1(k)=min(ρi-1(k)vf,i-1(k),qm,i-1) (3)
Ri(k)=min(qm,i,wiJi(k))) (4)
其中,Si-1(k)为元胞i-1在自由流状态下的最大发送流率,Ri(k)为允许进入元胞i的最大接收流量,vf,i-1表示元胞i-1处于自由流状态下的平均车速;qm,i、qm,i-1分别表示元胞i、i-1的最大通行流量;wi表示元胞i处于拥堵状态下的反向激波速度;ρJ表示车流完全堵塞下的车密度。
4.根据权利要求1所述的基于元胞传输模型的建模方法,其特征在于,
根据S1中的三类元胞中的上匝道元胞,元胞内存在上匝道和主线两股车流的合流区
Figure FDA0002998311200000021
Figure FDA0002998311200000022
Figure FDA0002998311200000023
其中,ri(k)为元胞i内由匝道i进入主线元胞i的流量,di(k)为与元胞i相连的上匝道i的交通需求;αi(k)为上匝道车流占合流区总车流的比例;
Figure FDA0002998311200000024
表示由主线元胞i-1进入元胞i的车流量。
5.根据权利要求4所述的基于元胞传输模型的建模方法,其特征在于,
通过下式对所述的交通需求di(k)进行预测:
di(k)=min(Di(k),ql,c×l) (17)
其中,Di(k)为前1小时内收费站的车流量,l为实施管控后收费站剩余车道数,ql,c为收费站单车道的最大通行能力。
6.根据权利要求1所述的基于元胞传输模型的建模方法,其特征在于,
根据S1中的三类元胞中的下匝道元胞,元胞内存在下匝道和主线两股车流的分流区
Figure FDA0002998311200000025
Figure FDA0002998311200000026
其中,βi(k)为元胞i内分流区进入下匝道的车流比例,
Figure FDA0002998311200000027
为输入元胞i+1的流量。
7.一种元胞传输模型的参数标定方法,其特征在于,包括权利要求1-6中任一项所述的建模方法,还包括:
S3:基于高速主线上起始终止路段中的ETC门架和匝道的卡口检测器获取的历史交通数据,采用基于最小二乘法的元胞传输模型参数标定方法对高速路网元胞传输理论模型的基本参数进行标定;并估算因事件引起的道路剩余通行能力;
S4:对匝道分流区、合流区的流量比例进行估算,并对匝道管控方案实施后的上匝道交通需求进行预测。
8.根据权利要求7所述的元胞传输模型的参数标定方法,其特征在于,所述的S3中,元胞中的基本参数qm、vf、w、ρJ的标定的过程如下:
(1)根据高速路段交通宏观基本图的理论模型,车流量q和车密度ρ的相关关系用下式表示:
q=min{ρvf,qm,w(ρJ-ρ)} (10)
(2)基于ETC门架和匝道的卡口检测器获取的历史交通数据,假设数据离散化后,时间步k=0,1,2…,则元胞i的历史车流量数据q′i、平均车速v′i和车密度ρ′i分别表示为:
q′i=[q′i(0),q′i(1),…,q′i(k)]T (11)
v′i=[v′i(0),v′i(1),…,v′i(k)]T (12)
ρ′i=[ρ′i(0),ρ′i(1),…,ρ′i(k)]T=[q′i(0)/v′i(0),q′i(1)/v′i(1),…,q′i(k)/v′i(k)]T (13)
(3)利用最小二乘法对公式(10)进行标定,标定方法用下式表示:
Figure FDA0002998311200000031
Figure FDA0002998311200000032
vf≤vmax,qm≤qmax,ρJ≤ρmax (16)
其中,
Figure FDA0002998311200000033
为qi(k)基于历史交通数据的的估计值;
(4)运用非线性优化问题求解器求解公式(14)~(16),获得元胞的基本参数qm、vf、w、ρJ
当道路发生拥堵时,用剩余通行能力q′m代替qm,q′m的标定的方法如下:
①对于影响车道数n已知的事件,则按照
Figure FDA0002998311200000034
计算,N为路段总车道数;
②对于影响车道数未知的事件,q′m取事件发生路段下游检测器检测到的车流量值。
9.根据权利要求7所述的元胞传输模型的参数标定方法,其特征在于,所述的S4中,
对于有上匝道元胞,合流区上匝道流量比例αi(k)的标定方式为,将上匝道卡口检测器和主线ETC门架的车流数据基于车牌信息进行匹配,并计算合流区的前1小时内上匝道的车流量占比;
对于有下匝道元胞,下匝道分流比例βi(k)的标定方法为:将主线ETC门架和下匝道卡口检测器的车流数据基于车牌信息进行匹配,计算分流区的前1小时内下匝道的车流量占比。
CN202110337985.7A 2021-03-30 2021-03-30 一种基于元胞传输模型的建模方法及参数标定方法 Active CN113205676B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110337985.7A CN113205676B (zh) 2021-03-30 2021-03-30 一种基于元胞传输模型的建模方法及参数标定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110337985.7A CN113205676B (zh) 2021-03-30 2021-03-30 一种基于元胞传输模型的建模方法及参数标定方法

Publications (2)

Publication Number Publication Date
CN113205676A true CN113205676A (zh) 2021-08-03
CN113205676B CN113205676B (zh) 2022-03-11

Family

ID=77025841

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110337985.7A Active CN113205676B (zh) 2021-03-30 2021-03-30 一种基于元胞传输模型的建模方法及参数标定方法

Country Status (1)

Country Link
CN (1) CN113205676B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114387781A (zh) * 2021-12-30 2022-04-22 北京建筑大学 车辆引导控制方法
CN114627658A (zh) * 2022-04-22 2022-06-14 河北上元智能科技股份有限公司 一种重大特殊车队通行高速公路的交通控制方法
CN114973671A (zh) * 2022-05-23 2022-08-30 浙江数智交院科技股份有限公司 公路网od数据处理方法、装置、设备及存储介质

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105930614A (zh) * 2016-05-13 2016-09-07 东南大学 一种针对可变限速控制的元胞传输模型参数标定与验证方法
CN112885088A (zh) * 2021-01-25 2021-06-01 浙江工业大学 一种基于动态交通流的多匝道协调控制方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105930614A (zh) * 2016-05-13 2016-09-07 东南大学 一种针对可变限速控制的元胞传输模型参数标定与验证方法
CN112885088A (zh) * 2021-01-25 2021-06-01 浙江工业大学 一种基于动态交通流的多匝道协调控制方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
黄超,陈日强: "《基于元胞传输模型的高速公路交通流仿真模型》", 《中国交通信息化》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114387781A (zh) * 2021-12-30 2022-04-22 北京建筑大学 车辆引导控制方法
CN114387781B (zh) * 2021-12-30 2024-04-09 北京建筑大学 车辆引导控制方法
CN114627658A (zh) * 2022-04-22 2022-06-14 河北上元智能科技股份有限公司 一种重大特殊车队通行高速公路的交通控制方法
CN114627658B (zh) * 2022-04-22 2023-10-27 河北上元智能科技股份有限公司 一种重大特殊车队通行高速公路的交通控制方法
CN114973671A (zh) * 2022-05-23 2022-08-30 浙江数智交院科技股份有限公司 公路网od数据处理方法、装置、设备及存储介质
CN114973671B (zh) * 2022-05-23 2024-03-12 浙江数智交院科技股份有限公司 公路网od数据处理方法、装置、设备及存储介质

Also Published As

Publication number Publication date
CN113205676B (zh) 2022-03-11

Similar Documents

Publication Publication Date Title
CN112800550B (zh) 一种基于元胞传输模型的匝道流量管控效果仿真预测方法
CN113205676B (zh) 一种基于元胞传输模型的建模方法及参数标定方法
Srivastava et al. Empirical observations of capacity drop in freeway merges with ramp control and integration in a first-order model
Papamichail et al. Coordinated ramp metering for freeway networks–a model-predictive hierarchical control approach
CN102005122B (zh) 多匝道调节方法及其系统
Li et al. Adaptive coordinated traffic control for stochastic demand
Kurzhanskiy et al. Guaranteed prediction and estimation of the state of a road network
Mahut et al. Calibration and application of a simulation-based dynamic traffic assignment model
CN104484514A (zh) 一种规避拥堵路段的弹性公交线路设计方法
Papamichail et al. Balancing of queues or waiting times on metered dual-branch on-ramps
Lee et al. Group-based approach to predictive delay model based on incremental queue accumulations for adaptive traffic control systems
Li et al. Optimal control of motorways by ramp metering, variable speed limits, and hard-shoulder running
Chitturi et al. Work zone queue length and delay methodology
Zhao et al. CTM based real-time queue length estimation at signalized intersection
Hou et al. Development of travel time functions for disrupted urban arterials with microscopic traffic simulation
Rarità et al. Sensitivity analysis of permeability parameters for flows on Barcelona networks
Chow et al. TOPL: Tools for operational planning of transportation networks
Akgungor et al. Analytical delay models for signalized intersections
Mehran et al. Considering travel time reliability and safety for evaluation of congestion relief schemes on expressway segments
Ghanbarikarekani et al. Optimization of signal timing of intersections by internal metering of queue time ratio of vehicles in network scale
Afandizadeh Zargari et al. Improving SCATS operation during congestion periods using internal/external traffic metering strategy
Feng et al. Improving Minnesota's stratified ramp control strategy
Liu et al. A two-stage stochastic model considering randomness of demand in variable speed limit and boundary flow control
Koutsopoulos et al. Effect of intersection delay modeling on the performance of traffic equilibrium models
Tian Development and evaluation of operational strategies for providing an integrated diamond interchange ramp-metering control system

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
CP01 Change in the name or title of a patent holder

Address after: 303 Wenhui Road, Hangzhou, Zhejiang 310000

Patentee after: Zhejiang Gaoxin Technology Co.,Ltd.

Address before: 303 Wenhui Road, Hangzhou, Zhejiang 310000

Patentee before: ZHEJIANG HIGHWAY INFORMATION ENGINEERING TECHNOLOGY CO.,LTD.

CP01 Change in the name or title of a patent holder