CN113204734A - 基于排队论的过饱和状态下交通系统多尺度供需关系的系统建模方法 - Google Patents
基于排队论的过饱和状态下交通系统多尺度供需关系的系统建模方法 Download PDFInfo
- Publication number
- CN113204734A CN113204734A CN202110350726.8A CN202110350726A CN113204734A CN 113204734 A CN113204734 A CN 113204734A CN 202110350726 A CN202110350726 A CN 202110350726A CN 113204734 A CN113204734 A CN 113204734A
- Authority
- CN
- China
- Prior art keywords
- time
- function
- delay
- queue length
- virtual queue
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000005094 computer simulation Methods 0.000 title claims abstract description 10
- 230000008569 process Effects 0.000 claims abstract description 13
- 238000011156 evaluation Methods 0.000 claims abstract description 6
- 230000001186 cumulative effect Effects 0.000 claims description 21
- 238000005457 optimization Methods 0.000 claims description 6
- 238000009440 infrastructure construction Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 3
- 238000007726 management method Methods 0.000 description 3
- 230000003068 static effect Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011144 upstream manufacturing Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- 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"
-
- 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/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- 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
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/40—Business processes related to the transportation industry
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- General Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Mathematical Physics (AREA)
- Tourism & Hospitality (AREA)
- Development Economics (AREA)
- General Business, Economics & Management (AREA)
- Game Theory and Decision Science (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Educational Administration (AREA)
- Pure & Applied Mathematics (AREA)
- Quality & Reliability (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Traffic Control Systems (AREA)
Abstract
本发明公开了一种基于排队论的过饱和状态下交通系统多尺度供需关系的系统建模方法,通过引入动力系统方程来描述具有虚拟队列演化过程的确定性队列模型,并基于多项式函数近似逼近的到达率,解析地推导出了交通系统性能的各种评价指标,比如,虚拟队列长度、平均延迟、物理队列长度和时变路段通行时间等,讨论了过饱和因子的不同取值范围的适用性情况,通过采集多源数据对该系统模型中的关键参数进行了校准,验证了该系统建模方法的有效性。本发明可以被用于联合优化需求管理政策和基础设施建设工作,针对不同规模的复杂且过饱和的动态排队系统,决策者可以应用本发明系统地制定需求侧和供给侧的拥堵缓解策略。
Description
技术领域
本发明涉及城市管理与控制领域,尤其是一种基于排队论的过饱和状态下交通系统多尺度供需关系的系统建模方法。
背景技术
在很多网络系统中,当需求时间上或空间上超过其容量时,都会出现拥堵现象。在世界各地,由于共享方式逐步趋向于机动性和自动化,使得交通系统正在经受一个重大的变革,而很多的区域规划组织和交通管理部门仍然面临着巨大的挑战,因为需要去减缓严重的交通拥堵以提高居民服务水平。原则上来说,有两种方式可以增强流动性以降低拥堵情况,一是主动地管理需求;二是通过基础设施建设来增加供给。因此,在需求管制和基础设施建设之间实现良好的平衡,而不是孤立地采取面向地方的减少拥堵战略,可以更好地缓解交通拥堵和提高城市流动性。
作为交通需求与供给的联系功能之一,宏观的流量-延迟函数已经得到了广泛的应用,比如,BPR(即美国公路局)函数自20世纪60年代以来就被广泛地用于城市规划,以此来量化需求和供给之间的非线性函数形式。而且,有着多项式形式的BPR函数容易刻画交通流和延迟(阻抗)的关系,同时,该函数计算效率高,易于标定,也在交通规划实践中易于实现;然而,交通规划部门早就意识到静态的BPR函数无法刻画动态交通流特征以及队列演化过程,尤其是队列的形成、传播和消散。此外,BPR函数也难以使用平均通行时间来描述一个高密度但低流量的过饱和瓶颈,而与静态交通分配模型相比,动态交通分配模型的目的是通过引入排队模型或其他类型的动态交通流模型来刻画交通拥堵的演化过程。但是,由于时间和空间维度的离散性(比如,元胞传输模型和路段传输模型),动态流量分配模型必须解决许多计算难题,而针对动态交通分配问题以及系统性能评估提供一个合理的时变路段通行时间函数是许多规划应用中一个新兴的研究需求。
发明内容
本发明所要解决的技术问题是:通过引入排队模型来刻画动态交通流特征以及队列的演化过程,从而弥补传统的静态BPR函数的不足,并提供了一种基于排队论的过饱和状态下交通系统多尺度供需关系的系统建模方法,能够被用于联合优化需求管理政策和基础设施能力建设工作,而且决策者也可以应用本发明系统地制定需求侧和供给侧的拥堵缓解策略。
技术方案:为解决上述技术问题,本发明提供了一种基于排队论的过饱和状态下交通系统多尺度供需关系的系统建模方法,包括如下步骤:
(1)通过引入动力系统方程来描述具有虚拟队列演化过程的确定性队列模型;
(2)根据确定性队列模型,并基于多项式函数近似逼近的到达率,获取交通系统性能评价指标,包括虚拟队列长度、时变延迟、总延迟、平均延迟、物理队列长度和时变路段通行时间,构建交通系统性能模型;
(3)确定过饱和因子的取值范围,根据过饱和因子不同的取值范围确定基于多项式函数近似逼近的到达率所适用的不同饱和程度的排队系统;
(4)采集道路上的传感器数据,以离开率、虚拟队列长度、时变延迟和时变达到率非负为约束,同时考虑过饱和因子的取值范围,对步骤(2)中构建的交通系统性能模型的参数进行校准,并与实际的观测数据进行对比。
进一步的,所述步骤(1)中,通过如下一系列的动态系统方程来描述具有虚拟队列演化过程的确定性队列模型:
约束条件:
λ(t)-μ(t)>0,t0<t<t2,λ(t)-μ(t)<0,t2<t<t3
其中λ(t),μ(t)和Q(t)分别表示时变的到达率,时变的离开率以及在任意时间t时的虚拟队列长度,λ(t)-μ(t)表示在时间t时的净流率,A(t),D(t)和W(t)分别表示累计到达数、累计离开数以及从t0到t时的总延迟,t0,t1,t2和t3分别表示到达率第一次超过离开率的时间即排队开始的时间,到达率为最大时的时间,队列长度达到最大时的时间以及队列完全消散时的时间。
进一步的,所述步骤(2)中,基于多项式函数近似逼近的到达率,确定虚拟队列长度函数Q(t),时变延迟函数w(t),总延迟函数W(t3),平均延迟函数w,物理队列长度函数Qp(t)以及时变路段通行时间函数tt;具体方法如下:
(2-1)拥堵时期的到达率函数λ(t)表示为一个三次多项式函数,表达式为:
其中γi表示第i阶变量的系数;在动态系统方程中存在边界条件λ(t0)=λ(t2)=μ,则到达率用净流率函数的因式分解形式为:其中是除了t0和t2之外的三次净流率函数的一个根,μ表示离开率,为常数,γ为三次多项式函数的形状参数;
(2-2)虚拟队列长度函数为:
其中m表示过饱和因子,通过队列开始形成到队列长度达到最大时的所用时间与整个拥堵期的总时间的比值来定义,其定义式如下:
(2-3)通过虚拟队列长度函数,得到时变延迟函数为:
(2-4)对虚拟队列长度函数进行积分,得到整个拥堵期t0到t3时的总延迟函数为:
W(t3)=γ·g(m)·(t3-t0)5
(2-5)将平均延迟定义为w=W/D,其中W表示总延迟,D表示在整个高峰期间的总出行需求;交通拥堵的持续时间为t3-t0=D/μ,则得到平均延迟函数为:
(2-6)基于中观层次下车辆在时空平面上的行驶轨迹,计算物理队列长度Qp(t):
其中Q(t)为虚拟队列长度,vf为自由流速度,vμ为实际的行驶速度;
(2-7)在整个高峰期内车辆在道路上的时变路段通行时间tt为:
其中tf表示自由流时间。
进一步的,所述步骤(3)中,当过饱和因子m和三次多项式函数的形状参数γ的取值范围不同时,提出的基于三次多项式函数近似逼近的到达率将适用于不同饱和程度的排队系统,具体的分类如下:
过饱和因子m的取值范围为:
当γ<0和m∈(2/3,3/4]时,基于三次多项式函数近似逼近的到达率适用于刻画轻度饱和的排队系统;当γ>0和m∈[1/2,2/3)时,该到达率适用于刻画轻度饱和以及过饱和的动态排队系统。
进一步的,所述步骤(4)中,采用约束优化模型和非线性最小二乘法对系统性能模型中的参数进行校准,具体方法如下:
使用排队系统中车辆数的累计数来校准离开率μ,并通过时变的虚拟队列长度和出行延迟来校准参数γ和过饱和因子m;首先将排队系统中车辆数的累计数、虚拟队列长度和出行延迟归一化到0-1的范围内,接着通过将平方和误差最小化来校准三个参数;得到如下优化模型:
约束条件:
μ>0,Q(t)≥0,w(t)≥0
上式中,|P|表示高峰期内的时间间隔;N(t),Q(t)和w(t)分别表示在t时,车辆的累计数量、虚拟队列长度以及延迟时间的理论值;N(t)=t·μ;和分别表示在t时,车辆的累计数量、虚拟队列长度以及延迟时间的观察值;Nmax,Qmax和wmax分别表示在t时,车辆的累计数量、虚拟队列长度以及延迟时间的理论的最大值,而 和分别表示在t时,车辆的累计数量、虚拟队列长度以及延迟时间的观察的最大值;Nmin,Qmin和wmin分别表示在t时,车辆的累计数量、虚拟队列长度以及延迟时间的理论的最小值,而和分别表示在t时,车辆的累计数量、虚拟队列长度以及延迟时间的观察的最小值。
有益效果:本发明针对过饱和的动态排队系统,提出了一个介于宏观和中观尺度之间的需求和供给的交叉解析模型,并基于多项式函数来近似逼近瓶颈处的车辆到达率,从中观的时变排队系统出发,建立了宏观的平均需求-延迟函数;本发明在中观尺度下,可以将拥堵期间的车辆时空轨迹映射到具有多项式函数近似逼近的时变到达率的动态排队系统方程中;本发明分析并推导了系统性能评价指标,包括时变情形下的队列长度和出行延迟等中观层次上的性能,以及平均延迟和时变路段通行时间等宏观层次上的性能;本发明针对不同饱和程度的排队系统,明确定义了过饱和率,并基于三次多项式函数近似逼近的到达率讨论了过饱和因子的不同取值范围的适用性情况;本发明提出的系统模型只需选择少量的参数,就可以很容易地与真实数据进行校准。因此,针对不同规模的复杂且过饱和的动态排队系统,决策者可以应用本发明系统地制定需求侧和供给侧的拥堵缓解策略。
附图说明
图1为本发明的方法流程示意图;
图2为本发明的确定性排队论模型解释图;
图3为本发明的区域DS1的时空速度图;
图4为本发明的区域DS1的速度和占有率变化过程图;
图5为本发明的区域DS1中不同探测器下的流量-道路占有率图;
图6为本发明的区域DS1的参数校准图。
具体实施方式
下面将结合采集的多源数据和说明书附图对本发明作进一步的说明。
如图1所示,一种基于排队论的过饱和状态下交通系统多尺度供需关系的系统建模方法,包括如下步骤:
(1)通过引入动力系统方程来描述具有虚拟队列演化过程的确定性队列模型;图2描述了虚拟队列演化过程的确定性排队论模型,其中,图2(a)展示了中观层次下车辆在时空平面上的行车轨迹,灰色双曲线描述了物理队列长度范围,黑色虚线描述了车辆的行车轨迹,在进入队列之前车辆以自由流速度vf行驶,而在进入队列之后以速度vμ行驶;图2(b)描述了可以被用来校准系统模型的多源数据的采集过程,包括有回路探测器数据、探测车辆数据等;图2(c)描述了针对过饱和交通系统的时变到达率曲线λ(t),以及常数离开率μ,很明显可以看到λ(t0)=λ(t2)=μ,而且时刻t2前后的阴影区域面积是相等的;图2(d)刻画了虚拟队列长度的演化过程,可以看到在时刻t2时的虚拟队列长度最大,而在时刻t0和t3的虚拟队列长度为0;图2(e)描述了系统中车辆的累计数,斜率为μ的灰色直线D(t)表示高峰期间的车辆累计离开曲线,而黑色曲线A(t)表示车辆的累计到达曲线,在t0之前没有发生拥堵,所以此时车辆的累计到达数和累计离开数相等,而在t0到t3的高峰期内,由于拥堵效应,会使得此时车辆的累计到达数超过累计离开数,在任意时刻t时的累计到达曲线和累计离开曲线的垂直差即为虚拟队列长度Q(t),而且由于累计离开曲线的斜率为常数,所以可以很容易就得到相应的出行延迟为w(t)。
图2所展示的虚拟队列演化过程的确定性排队论模型可以通过如下一系列的动态系统方程来规划:
约束条件:
λ(t0)=μ(t0) (5)
λ(t2)=μ(t2) (6)
λ(t)-μ(t)>0,t0<t<t2 (7)
λ(t)-μ(t)<0,t2<t<t3 (8)
Q(t0)=0 (10)
Q(t3)=0 (11)
其中λ(t),μ(t)和Q(t)分别表示时变的到达率,时变的离开率以及在任意时间t时的虚拟队列长度,λ(t)-μ(t)表示在时间t时的净流率,A(t),D(t)和W(t)分别表示累计到达数、累计离开数以及从t0到t时的总延迟,t0,t1,t2和t3分别表示到达率第一次超过离开率的时间即排队开始的时间,到达率为最大时的时间,队列长度达到最大时的时间以及队列完全消散时的时间。
(2)根据确定性队列模型,并基于多项式函数近似逼近的到达率,获取交通系统性能评价指标,包括虚拟队列长度、时变延迟、总延迟、平均延迟、物理队列长度和时变路段通行时间,构建交通系统性能模型。基于多项式函数近似逼近的到达率,解析地推导出虚拟队列长度函数Q(t),时变延迟函数w(t),总延迟函数W(t3),平均延迟函数w,物理队列长度函数Qp(t)以及时变路段通行时间函数tt;具体方法如下:
(2-1)拥堵时期的到达率函数λ(t)表示为一个三次多项式函数,表达式为:
其中γi表示第i阶变量的系数;
考虑到在动态系统方程中存在边界条件λ(t0)=λ(t2)=μ,则到达率用净流率函数的因式分解形式改写为:
(2-2)将式(12)代入式(3),然后对结果进行积分,并通过变换积分上下限就可以推导出虚拟队列长度函数为:
(2-4)通过将式(14)代入式(13),则可以进一步得到虚拟队列长度函数为:
其中m表示过饱和因子,通过队列开始形成到队列长度达到最大时的所用时间与整个拥堵期的总时间的比值来定义,其定义式如下:
(2-5)通过虚拟队列长度函数式(15),得到时变延迟函数为:
(2-6)对虚拟队列长度函数式(15)进行积分,得到整个拥堵期t0到t3时的总延迟函数为:
W(t3)=γ·g(m)·(t3-t0)5 (18)
(2-7)将平均延迟定义为w=W/D,其中W表示总延迟,D表示在整个高峰期间的总出行需求;交通拥堵的持续时间为t3-t0=D/μ,则得到平均延迟函数为:
(2-8)基于中观层次下车辆在时空平面上的行驶轨迹,通过虚拟队列长度Q(t),自由流速度vf和实际的行驶速度vμ来计算物理队列长度Qp(t),其公式为:
(2-9)根据平均延迟函数式(20),得到在整个高峰期内车辆在道路上的时变路段通行时间tt为:
其中tf表示自由流时间。
(3)确定过饱和因子的取值范围,根据过饱和因子不同的取值范围确定基于多项式函数近似逼近的到达率所适用的不同饱和程度的排队系统;
当过饱和因子m和三次多项式函数的形状参数γ的取值范围不同时,提出的基于三次多项式函数近似逼近的到达率将适用于不同饱和程度的排队系统,具体的分类如下:
根据步骤(2)中提出的基于三次多项式函数近似逼近的到达率,当γ<0时,通过式(13)可以推导出条件也成立,又因为m=(t2-t0)/(t3-t0),则基于上述条件对式(14)进行化简,就可以得到m的取值范围如下:
因此,过饱和因子m的取值范围为:
上述解析式并不是都适合于过饱和的动态排队系统,比如,当γ<0和m∈(2/3,3/4]时,基于三次多项式函数近似逼近的到达率更适用于刻画轻度饱和的排队系统;其原因是针对过度拥堵的排队系统,在拥堵期结束时的达到率λ(t)会明显地小于离开率μ,这是因为在接近t3时到达率λ(t)会急剧下降,此时λ(t3)的估计值可能为负值,这就违反了车流量为正的假设;当γ>0和m∈[1/2,2/3)时,基于三次多项式函数近似逼近的到达率适用于刻画轻度饱和以及过饱和的动态排队系统。
(4)采集道路上的传感器数据,以离开率、虚拟队列长度、时变延迟和时变达到率非负为约束,同时考虑过饱和因子的取值范围,对步骤(2)中构建的交通系统性能模型的参数进行校准,并与实际的观测数据进行对比。
考虑到车辆的到达率、虚拟队列长度和出行延迟均是非负的,不能运用传统的最小二乘法来进行校准,因此,本发明采用约束优化模型和非线性最小二乘法对系统性能模型中的参数进行校准,具体方法如下:
使用排队系统中车辆数的累计数来校准离开率μ,并通过时变的虚拟队列长度和出行延迟来校准参数γ和过饱和因子m;考虑到这三种校准所使用数据的单位和尺度均不同,因此,首先将排队系统中车辆数的累计数、虚拟队列长度和出行延迟归一化到0-1的范围内,接着通过将平方和误差最小化来校准三个参数;得到如下优化模型:
约束条件包括式(25)以及式(27)~(30):
μ>0 (27)
Q(t)≥0 (28)
w(t)≥0 (29)
上式中,|P|表示高峰期内的时间间隔;N(t),Q(t)和w(t)分别表示在t时,车辆的累计数量、虚拟队列长度以及延迟时间的理论值;N(t)=t·μ;Q(t)和w(t)分别通过式(15)和式(17)计算;和分别表示在t时,车辆的累计数量、虚拟队列长度以及延迟时间的观察值;Nmax,Qmax和wmax分别表示在t时,车辆的累计数量、虚拟队列长度以及延迟时间的理论的最大值,而和分别表示在t时,车辆的累计数量、虚拟队列长度以及延迟时间的观察的最大值;Nmin,Qmin和wmin分别表示在t时,车辆的累计数量、虚拟队列长度以及延迟时间的理论的最小值,而和分别表示在t时,车辆的累计数量、虚拟队列长度以及延迟时间的观察的最小值。
为了对步骤(2)中提出的交通系统性能模型中的关键参数进行校准,首先,本实例通过相关网站以及实地观测等方法,在洛杉矶的区域DS1内采集了相关数据,有关实验数据的详细描述如表1所示:
表1洛杉矶某区域的实验数据集的描述
针对区域DS1,首先可以应用采集的数据来绘制时空速度图以获得瓶颈的位置,如图3所示,可以看到区域DS1的瓶颈位置在Abs=13.51英里处,而且高峰期是从t0=13:10到t3=19:45;然后,通过绘制速度和占有率变化过程图(参见图4(a)~4(c))可以确定t0,t3和自由流速度vf,可以看到,在高峰期间瓶颈下游的车辆速度稳定在45英里/小时,而瓶颈处和瓶颈上游的车辆速度会急剧下降,同时,瓶颈下游的道路占有率低于0.11,而瓶颈处和瓶颈上游的道路占有率到达了0.25;接着,将瓶颈位置邻近下游的车辆累计数视为观察的车辆离开数的测量,并通过总通行时间减去自由流时间来计算观察的出行延迟,而针对虚拟队列长度的测量,可以先基于流量-道路占有率散点图(参见图5)来得到占有率的临界值,然后将占有率的临界值转换成密度,并最终用密度数据来计算观察的虚拟队列长度,如图5所示,可以看到瓶颈位于第二行第二列,此时排队系统的道路占有率的临界值接近于0.13,这将被用于去计算虚拟队列长度,并进一步计算物理队列长度,从而校准参数。
图6展示了本发明提出的基于排队论的过饱和状态下交通系统多尺度供需关系的系统建模方法有着较好的拟合效果,其中,图6(a)描述了区域DS1中参数μ的校准结果,图6(b)描述了区域DS1中物理队列长度的校准结果,而图6(c)描述了区域DS1中到达率以及μ的校准结果。校准的结果分别如下:针对区域DS1,μDS1=3860车辆/小时,γDS1=13车辆/小时4以及mDS1=0.537。可以发现,区域DS1的离开率变化较小,其原因是在区域DS1中大多是有着较高限速的高速公路,因此车辆行驶速度较稳定,不会出现走走停停现象;而针对特定的分析时段,区域DS1在交通严重拥堵时期的有效离开率是较高的,但相较于理论最大值是明显较低的,其主要原因是受到瓶颈处的交通密度、复杂的道路地形以及驾驶员行为的影响。在本发明中提出的交通系统性能模型的主要特征之一是,对于一个过饱和的交通系统,将出行高峰时刻t1时的到达率λ(t1)(或者说“需求”)和离开率μ(或者说“供给”)进行了量化分析,这可以进一步被定义为系统利用率ρ=λ(t1)/μ,而针对洛杉矶该区域的系统利用率为ρDS1=1.059,通过分析可以发现,针对洛杉矶的数据集DS1,其拥堵的主要原因是和长期累计的过剩需求有关。
上述实施例仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和等同替换,这些对本发明权利要求进行改进和等同替换后的技术方案,均落入本发明的保护范围。
Claims (5)
1.基于排队论的过饱和状态下交通系统多尺度供需关系的系统建模方法,其特征在于,包括如下步骤:
(1)通过引入动力系统方程来描述具有虚拟队列演化过程的确定性队列模型;
(2)根据确定性队列模型,并基于多项式函数近似逼近的到达率,获取交通系统性能评价指标,包括虚拟队列长度、时变延迟、总延迟、平均延迟、物理队列长度和时变路段通行时间,构建交通系统性能模型;
(3)确定过饱和因子的取值范围,根据过饱和因子不同的取值范围确定基于多项式函数近似逼近的到达率所适用的不同饱和程度的排队系统;
(4)采集道路上的传感器数据,以离开率、虚拟队列长度、时变延迟和时变达到率非负为约束,同时考虑过饱和因子的取值范围,对步骤(2)中构建的交通系统性能模型的参数进行校准,并与实际的观测数据进行对比。
2.根据权利要求1所述的方法,其特征在于,所述步骤(1)中,通过如下一系列的动态系统方程来描述具有虚拟队列演化过程的确定性队列模型:
约束条件:
λ(t)-μ(t)>0,t0<t<t2,λ(t)-μ(t)<0,t2<t<t3
其中λ(t),μ(t)和Q(t)分别表示时变的到达率,时变的离开率以及在任意时间t时的虚拟队列长度,λ(t)-μ(t)表示在时间t时的净流率,A(t),D(t)和W(t)分别表示累计到达数、累计离开数以及从t0到t时的总延迟,t0,t1,t2和t3分别表示到达率第一次超过离开率的时间即排队开始的时间,到达率为最大时的时间,队列长度达到最大时的时间以及队列完全消散时的时间。
3.根据权利要求2所述的方法,其特征在于,所述步骤(2)中,基于多项式函数近似逼近的到达率,确定虚拟队列长度函数Q(t),时变延迟函数w(t),总延迟函数W(t3),平均延迟函数w,物理队列长度函数Qp(t)以及时变路段通行时间函数tt;具体方法如下:
(2-1)拥堵时期的到达率函数λ(t)表示为一个三次多项式函数,表达式为:
其中γi表示第i阶变量的系数;在动态系统方程中存在边界条件λ(t0)=λ(t2)=μ,则到达率用净流率函数的因式分解形式为:其中是除了t0和t2之外的三次净流率函数的一个根,μ表示离开率,为常数,γ为三次多项式函数的形状参数;
(2-2)虚拟队列长度函数为:
其中m表示过饱和因子,通过队列开始形成到队列长度达到最大时的所用时间与整个拥堵期的总时间的比值来定义,其定义式如下:
(2-3)通过虚拟队列长度函数,得到时变延迟函数为:
(2-4)对虚拟队列长度函数进行积分,得到整个拥堵期t0到t3时的总延迟函数为:
W(t3)=γ·g(m)·(t3-t0)5
(2-5)将平均延迟定义为w=W/D,其中W表示总延迟,D表示在整个高峰期间的总出行需求;交通拥堵的持续时间为t3-t0=D/μ,则得到平均延迟函数为:
(2-6)基于中观层次下车辆在时空平面上的行驶轨迹,计算物理队列长度Qp(t):
其中Q(t)为虚拟队列长度,vf为自由流速度,vμ为实际的行驶速度;
(2-7)在整个高峰期内车辆在道路上的时变路段通行时间tt为:
其中tf表示自由流时间。
5.根据权利要求3所述的方法,其特征在于,所述步骤(4)中,采用约束优化模型和非线性最小二乘法对系统性能模型中的参数进行校准,具体方法如下:
使用排队系统中车辆数的累计数来校准离开率μ,并通过时变的虚拟队列长度和出行延迟来校准参数γ和过饱和因子m;首先将排队系统中车辆数的累计数、虚拟队列长度和出行延迟归一化到0-1的范围内,接着通过将平方和误差最小化来校准三个参数;得到如下优化模型:
约束条件:
μ>0,Q(t)≥0,w(t)≥0
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110350726.8A CN113204734B (zh) | 2021-03-31 | 2021-03-31 | 基于排队论的过饱和状态下交通系统多尺度供需关系的系统建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110350726.8A CN113204734B (zh) | 2021-03-31 | 2021-03-31 | 基于排队论的过饱和状态下交通系统多尺度供需关系的系统建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113204734A true CN113204734A (zh) | 2021-08-03 |
CN113204734B CN113204734B (zh) | 2024-04-09 |
Family
ID=77026013
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110350726.8A Active CN113204734B (zh) | 2021-03-31 | 2021-03-31 | 基于排队论的过饱和状态下交通系统多尺度供需关系的系统建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113204734B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE10108611A1 (de) * | 2001-02-22 | 2002-09-05 | Daimler Chrysler Ag | Verfahren zur Simulation und Prognose der Bewegung von Einzelfahrzeugen auf einem Verkehrswegenetz |
US20050122987A1 (en) * | 2003-12-09 | 2005-06-09 | Michael Ignatowski | Apparatus and method for modeling queueing systems with highly variable traffic arrival rates |
CN106251649A (zh) * | 2016-08-09 | 2016-12-21 | 南京航空航天大学 | 基于缓解过饱和状态下道路交叉口拥堵状况的控制策略 |
CN106652458A (zh) * | 2017-02-20 | 2017-05-10 | 东南大学 | 基于虚拟车辆轨迹重构的在线城市道路路径行程时间估计方法 |
CN109360432A (zh) * | 2018-11-27 | 2019-02-19 | 南京航空航天大学 | 一种基于延误最小和饱和度均衡的多交叉口的控制方法 |
-
2021
- 2021-03-31 CN CN202110350726.8A patent/CN113204734B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE10108611A1 (de) * | 2001-02-22 | 2002-09-05 | Daimler Chrysler Ag | Verfahren zur Simulation und Prognose der Bewegung von Einzelfahrzeugen auf einem Verkehrswegenetz |
US20050122987A1 (en) * | 2003-12-09 | 2005-06-09 | Michael Ignatowski | Apparatus and method for modeling queueing systems with highly variable traffic arrival rates |
CN106251649A (zh) * | 2016-08-09 | 2016-12-21 | 南京航空航天大学 | 基于缓解过饱和状态下道路交叉口拥堵状况的控制策略 |
CN106652458A (zh) * | 2017-02-20 | 2017-05-10 | 东南大学 | 基于虚拟车辆轨迹重构的在线城市道路路径行程时间估计方法 |
CN109360432A (zh) * | 2018-11-27 | 2019-02-19 | 南京航空航天大学 | 一种基于延误最小和饱和度均衡的多交叉口的控制方法 |
Non-Patent Citations (1)
Title |
---|
马书红;孙朝旭;: "基于系统动力学的城市交通系统供需分析模型", 长安大学学报(社会科学版), no. 03 * |
Also Published As
Publication number | Publication date |
---|---|
CN113204734B (zh) | 2024-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Srivastava et al. | Empirical observations of capacity drop in freeway merges with ramp control and integration in a first-order model | |
CN111161537B (zh) | 考虑拥堵叠加效应的公路拥堵态势预测方法 | |
CN114783183B (zh) | 一种基于交通态势算法的监控方法与系统 | |
CN106781499B (zh) | 一种交通网络效率评价系统 | |
WO2023216793A1 (zh) | 一种混合交通流环境高速公路瓶颈路段动态限速控制方法 | |
Rakha et al. | Modeling inclement weather impacts on traffic stream behavior | |
CN109118770A (zh) | 一种基于交通监测数据的路段通行能力挖掘方法 | |
Zu et al. | Real-time energy-efficient traffic control via convex optimization | |
CN105006149A (zh) | 交通路况估计动态迭代方法 | |
CN104182633B (zh) | 分层次的交通运行评价方法 | |
CN113159374A (zh) | 一种数据驱动的城市交通流速模式识别与实时预测预警方法 | |
Xie et al. | Modeling safety impacts of off-hour delivery programs in urban areas | |
Aron et al. | Estimating travel time distribution for reliability analysis | |
CN116052435A (zh) | 一种城市道路拥堵影响范围界定与道路影响力计算方法 | |
CN116186994A (zh) | 基于微观交通仿真的高速公路拥堵溯源分析及管控方法 | |
Chitturi et al. | Work zone queue length and delay methodology | |
CN111429717A (zh) | 一种城市快速路道路运行能力评估方法 | |
Ivanjko et al. | Simulational analysis of two controllers for variable speed limit control | |
CN106991817B (zh) | 多级路网路段交通容量确定方法 | |
CN113204734A (zh) | 基于排队论的过饱和状态下交通系统多尺度供需关系的系统建模方法 | |
CN113849976B (zh) | 规划用地的开发强度评估方法、装置及设备 | |
Geroliminis et al. | Development of the Next Generation Stratified Ramp Metering Algorithm Based on Freeway Density | |
Liang et al. | Identifying expressway accident black spots based on the secondary division of road units | |
Li et al. | An approach to modeling the impact of traffic incident on urban expressway | |
Moses et al. | Evaluation of the 2010 highway capacity manual urban street free-flow speed prediction model |
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 |