CN102521438B - 交通冲突仿真两阶段参数标定方法 - Google Patents

交通冲突仿真两阶段参数标定方法 Download PDF

Info

Publication number
CN102521438B
CN102521438B CN2011103942487A CN201110394248A CN102521438B CN 102521438 B CN102521438 B CN 102521438B CN 2011103942487 A CN2011103942487 A CN 2011103942487A CN 201110394248 A CN201110394248 A CN 201110394248A CN 102521438 B CN102521438 B CN 102521438B
Authority
CN
China
Prior art keywords
conflict
traffic
emulation
crossing
centerdot
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.)
Active
Application number
CN2011103942487A
Other languages
English (en)
Other versions
CN102521438A (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.)
Southeast University
Original Assignee
Southeast University
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 Southeast University filed Critical Southeast University
Priority to CN2011103942487A priority Critical patent/CN102521438B/zh
Publication of CN102521438A publication Critical patent/CN102521438A/zh
Application granted granted Critical
Publication of CN102521438B publication Critical patent/CN102521438B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

一种交通冲突仿真的两阶段参数标定方法,步骤包括第一阶段对交通仿真软件(VISSIM)中影响交通运行的参数标定和第二阶段对冲突分析软件(SSAM)中仿真冲突判定阈值的标定。本发明在现有交通仿真模型参数标定方法的基础上,提出了面向交通冲突仿真的两阶段参数标定方法,克服了传统参数标定方法模拟仿真冲突精度低的缺陷,在利用微观仿真软件对交通设计方案进行安全研究方面具有实际工程运用价值。

Description

交通冲突仿真两阶段参数标定方法
技术领域
本发明属于交通仿真领域,为一种利用微观仿真软件VISSIM和冲突分析软件SSAM进行交通冲突模拟的参数标定方法。
背景技术
微观仿真技术的不断成熟,使其在交通工程领域中的应用日益广泛,特别是在通行能力研究方面,逐渐形成了利用微观仿真技术模拟实际交通运行状况的参数标定流程。然而在交通冲突仿真方面的工程实践表明,传统面向通行能力研究的参数标定流程在模拟实际交通冲突方面存在很大的缺陷,并不能正确反映实际交通冲突状况。
发明内容
如何改善传统微观交通仿真模型的标定流程,使其能够应用于交通冲突仿真是本发明主要解决的问题。
为了弥补传统微观仿真模型参数标定过程不能正确反映实际交通冲突状况的不足,本发明提出了交通冲突仿真两阶段参数标定方法,包括以下步骤:
第1步,现场数据采集,建立交叉口VISSIM仿真模型。现场数据采集主要包括交叉口几何尺寸、交通流量、车速、大小车比例、交通冲突以及信号配时等数据的采集。然后根据交叉口航拍图或者CAD图建立交叉口仿真模型,并定义仿真模型中交通流量、大小车比例、车辆速度以及交通信号配时等信息。
第2步,参数校准。对VISSIM中车辆跟驰模型、变道模型以及优先规则等所包含的参数进行校准。初次参数校准可采用软件默认值,若默认值无法满足研究需求,可根据经验和国内外研究成果设定参数的变化区间和变化幅度,重复试验进行参数校准。
第3步,模型检验。首先选择能够反映交叉口交通运行状况的检验指标,如车头时距、排队长度等,然后进行皮尔森卡方检验计算该指标实测值与仿真输出值的之间差异,计算公式为公式(1)。
χ 2 = Σ i = 1 k ( s i - r i ) 2 r i - - - ( 1 )
其中k表示对检验指标数值划分成k个区间,si代表检验指标的仿真输出值在第i个区间的个数,ri代表检验指标的实际观测值在第i个区间的个数。在给定显著性水平α下,若
Figure GDA00002806836700021
则认为两者之间的差异过大,返回第2步,继续进行参数校准;反之则认为仿真模型能够反映实际交通运行状况,进行第4步。
第4步,运行仿真模型,输出轨迹文件。通过第3步模型检验之后,运行仿真模型,输出车辆的轨迹文件。
第5步,仿真冲突阈值校准。冲突分析软件SSAM通过计算车辆之间的距冲突时间(TTC)来判定是否有交通冲突发生。SSAM中默认的冲突阈值为1.5s,当计算得到的TTC值小于或等于1.5s时,认为有冲突产生;当TTC值大于1.5s时,认为没有冲突产生。而实际冲突观测中,观测者通过观察车辆是否有明显的避险行为,如突然减速、变向等来判断是否有冲突产生。因此为了使仿真冲突与实际冲突具有一致性,必须对仿真冲突阈值进行校准。同样,仿真冲突阈值的初始值可以定为软件默认值,若初始值不可接受,则可设定其变化区间和变化幅度,重复试验进行参数校准。
第6步,SSAM冲突分析。运行冲突分析软件SSAM,对交叉口仿真模型输出的车辆轨迹文件进行冲突分析,输出仿真冲突数据。
第7步,冲突检验。利用公式(2)计算交叉口实际观测到的冲突与分析得出的仿真冲突的差异:
MAPE = 1 n Σ i = 1 n | a i - b i a i | * 100 % - - - ( 2 )
其中MAPE代表实际冲突数与仿真冲突数之间的平均绝对百分比误差,n代表交叉口个数,ai代表第i个交叉口的实际冲突个数,bi代表第i个交叉口的仿真冲突个数。若计算得到实际冲突与仿真冲突之间的平均绝对百分比误差大于10%,则认为误差过大,返回第5步重新进行仿真冲突阈值校准;反之,若两者之间的平均绝对百分比误差不大于10%,认为误差可接受,则完成整个参数标定过程。
附图说明
图1为本发明实施例的现场数据采集示意图。
图2为本发明进行参数标定的流程图。
具体实施方式
本发明的交通冲突仿真两阶段参数标定方法步骤包括:
1、建立交叉口VISSIM仿真模型:
第1步,构建交叉口模型。微观交通仿真软件VISSIM是在微观交通仿真领域得到广泛应用的仿真软件,构建交叉口仿真模型首先需要对交叉口进行现场数据采集,包括交叉口几何尺寸,交通流量,大小车比例,交通信号配时等信息。然后根据交叉口航拍图或者CAD图建立仿真交叉口,并根据实际观测数据定义仿真模型中交通流量、大小车比例、车辆速度以及交通信号配时等信息。
第2步,参数校准。VISSIM通过车辆跟驰模型、变道模型以及优先规则等来定义仿真模型中车辆的运行动作,研究者可以通过对这些模型的参数进行修正来模拟交叉口车辆的实际运行状况。初次参数校准可采用软件默认值,若默认值无法满足研究需求,可根据经验和国内外研究成果设定参数的变化区间和变化幅度,并重复试验进行参数校准。
第3步,模型检验。首先选择能够反映交叉口交通运行状况的检验指标,如车头时距、排队长度等,然后进行皮尔森卡方检验,计算该指标实测值与仿真输出值的之间差异,计算公式为公式(1)。
χ 2 = Σ i = 1 k ( s i - r i ) 2 r i - - - ( 1 )
其中k代表对检验指标数值划分成k个区间,si代表检验指标的仿真输出值在第i个区间的个数,ri代表检验指标的实际观测值在第i个区间的个数。在给定显著性水平α下,若
Figure GDA00002806836700032
则认为两者之间的差异过大,返回第2步,继续进行参数校准;反之则认为仿真模型能够反映实际交通运行状况,完成交叉口仿真建模。
2、面向交通冲突的参数标定:
第1步,运行权利要求1所建立的交叉口VISSIM仿真模型,输出车辆轨迹文件。
第2步,仿真冲突阈值校准。冲突分析软件SSAM通过计算车辆之间的距冲突时间TTC来判定是否有交通冲突发生。SSAM中默认的冲突阈值为1.5s,当计算得到车辆之间的TTC值小于或等于1.5s时,认为有冲突产生;当TTC值大于1.5s时,认为没有冲突产生。而实际冲突观测中,观测者通过观察车辆是否有明显的避险行为,如突然减速、变向等来判断是否有冲突产生。因此为了使仿真冲突与实际冲突具有一致性,必须对仿真冲突阈值进行校准。同样,仿真冲突阈值的初始值可以定为软件默认值,若初始值不可接受,则可通过设定其变化区间和变化幅度,并重复试验进行参数校准。
第3步,SSAM冲突分析。运行冲突分析软件SSAM,对交叉口仿真模型输出的车辆轨迹文件进行冲突分析,输出仿真冲突数据。
第4步,冲突检验。利用公式(2)计算交叉口实际观测到的冲突与分析得出的仿真冲突的差异:
MAPE = 1 n Σ i = 1 n | a i - b i a i | * 100 % - - - ( 2 )
其中MAPE代表实际冲突数与仿真冲突数之间的平均绝对百分比误差,n代表对交叉口进行冲突仿真的小时数,ai代表交叉口第i小时内的实际冲突个数,bi代表交叉口第i小时内的仿真冲突个数。若计算得到实际冲突与仿真冲突之间的平均绝对百分比误差大于10%,则认为误差不可接受,重新进行仿真冲突阈值校准;反之,若两者之间的平均绝对百分比误差不大于10%,认为误差可接受,则完成整个参数标定过程。
下面结合南京某交叉口对本标定方法进行演示。
1)现场数据采集,建立交叉口VISSIM仿真模型。调查人员利用米尺、测距滚轮测量交叉口几何尺寸,并在交叉口各进口道分别架设一台摄像机拍摄进口道和交叉口内部的交通运行状况,如图1所示。然后,通过对录像进行分析,获取交叉口的交通流量、大小车比例以及交通冲突数据,共获得该交叉口5个小时的交通信息。根据交叉口几何尺寸和分析得到的交通流量,大小车比例等信息,建立交叉口仿真模型。
2)参数校准。交叉口仿真模型中车辆跟驰模型、变道模型以及优先规则等包含的参数均采用软件默认值。
3)模型检验。选取车头时距分布作为检验模型有效性的评价指标,运行仿真模型,输出车头时距,此外通过对交叉口视频进行分析,获取实际车头时距。根据样本数据分类原则,车头时距共划分成7组,进行皮尔森卡方检验,由公式1可计算得到,χ2值为15.72,而在显著性水平α下,
Figure GDA00002806836700042
此时
Figure GDA00002806836700043
认为仿真车头时距分布与实际车头时距分布相差过大,需要对默认参数进行校准。
根据以往研究,影响仿真模型中车头时距分布的主要影响因素为期望安全距离的加项BXadd和乘项BXmulti,两者的默认值分别为2.0和3.0,设定两者的变化范围都为(1.0,5.0),变化步长为0.5。反复试验得到:当BXadd=2.5,BXmulti=6.0时,仿真车头时距与实际车头时距之间的皮尔森卡方值χ2=2.98,该值小于
Figure GDA00002806836700044
因此通过模型检验。
4)运行交叉口仿真模型,输出在5个不同小时交通流量下的车辆轨迹文件。
5)运行SSAM进行交通冲突分析。在默认的仿真交通冲突阈值(TTC=1.5s)下,运行冲突分析软件SSAM,对第4步得到的5个车辆轨迹文件分别进行冲突分析。
6)冲突检验。根据该交叉口5个小时的实际交通冲突数和仿真冲突数(如表1所示),由公式2计算实际交通冲突与仿真交通冲突之间的差异:
Figure GDA00002806836700051
计算得到的MAPE值大于10%,说明在默认仿真冲突阈值情况下,仿真冲突与实际冲突之间差异较大,需要对仿真冲突阈值进行校准。
根据国内外研究成果,设定仿真冲突阈值TTC的变化范围为1.0s—2.0s,变化幅度为0.1s。经过反复试验发现,在TTC=1.6s时,仿真冲突与实际冲突之间的差异最小,具体数据见表1。计算得到此时的MAPE值为8.3%<10%,通过冲突检验。至此,完成交通冲突仿真两阶段参数标定过程。
表1实际交通冲突数与仿真冲突数
编号 实际冲突数 仿真冲突数1 仿真冲突数2
1 20 7 19
2 21 10 22
3 17 7 18
4 16 6 13
5 14 12 15
注:仿真冲突数1是指在默认仿真冲突阈值下,即TTC=1.5s,SSAM分析得到的仿真冲突数;仿真冲突数2是指在TTC=1.6s时,SSAM分析得到的仿真冲突数。
由于该交叉口事故发生较为频繁,拟对其进行交通安全改善,提出以下两种改善方案:方案一,交叉口道路限速由60km/h降低到50km/h;方案二,增设左转专用信号。传统的安全改善效果评价需要采集方案实施之后几年内的交通事故数据,或在实施之后进行现场人工冲突观测,这两种方法不但需要耗费大量的人力物力,而且都是事后评价方法,无法预先比较不同方案的优劣。然而,利用本发明提出的交通冲突仿真两阶段参数标定方法,可在仿真模型中实现不同的交通安全改善方案,模拟改善后的交通运行状况,并输出仿真冲突数据,从而进行方案比选。以本交叉口为例,模拟实施方案一后,交叉口1小时内平均产生仿真冲突16起;模拟实施方案二后,交叉口1小时内平均产生仿真冲突9起,因此就实施效果而言,方案二优于方案一。

Claims (4)

1.一种交通冲突仿真两阶段参数标定方法,其特征是包括步骤:
1)现场数据采集,建立交叉口VISSIM仿真模型;
2)参数校准:对VISSIM中车辆跟驰模型、变道模型以及优先规则所包含的参数进行校准;
3)模型检验:
首先,选择能够反映交叉口交通运行状况的检验指标;
然后,进行皮尔森卡方检验计算该检验指标实测值与该检验指标仿真输出值的之间差异,计算公式为公式(1): &chi; 2 = &Sigma; i = 1 k ( s i - r i ) 2 r i &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; ( 1 ) ,
其中,k对检验指标数值划分成k个区间,si代表检验指标的仿真输出值在第i个区间的个数,ri代表检验指标的实测值在第i个区间的个数;
在给定显著性水平α下,若
Figure FDA00003492199900012
则返回步骤2),继续进行参数校准;若
Figure FDA00003492199900013
则认为仿真模型能够反映实际交通运行状况,进行步骤4);
4)运行仿真模型,输出轨迹文件:通过步骤3)的模型检验之后,运行仿真模型,输出车辆的轨迹文件;
5)仿真冲突阈值校准:冲突分析软件SSAM进行冲突分析是通过计算车辆之间的距冲突时间TTC,并与仿真冲突阈值进行比较来判定是否有交通冲突发生;仿真冲突阈值的初始值定为软件默认值,若初始值不可接受,则设定仿真冲突阈值变化区间和变化幅度,重复试验进行参数校准;
6)SSAM冲突分析:在默认的仿真交通冲突阈值1.5s下,运行冲突分析软件SSAM,对步骤4)得到的车辆轨迹文件分别进行冲突分析;
7)冲突检验:利用公式(2)计算交叉口实际观测到的冲突与分析得出的仿真冲突的差异: MAPE = 1 n &Sigma; i = 1 n | a i - b i a i | * 100 % &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; ( 2 )
其中,MAPE代表实际冲突数与仿真冲突数之间的平均绝对百分比误差,n代表交叉口个数,ai代表第i个交叉口的实际冲突个数,bi代表第i个交叉口的仿真冲突个数;
若计算得到实际冲突与仿真冲突二者之间的平均绝对百分比误差大于10%,则认为误差过大,返回步骤5)重新进行仿真冲突阈值校准;
若计算得到二者之间的平均绝对百分比误差不大于10%,认为误差可接受,则完成整个参数标定过程;
所述步骤1)~4)为第一阶段参数标定,所述步骤5)~7)为第二阶段参数标定。
2.根据权利要求1所述的交通冲突仿真两阶段参数标定方法,其特征是所述步骤1)中,
首先,采集现场数据,现场数据采集主要包括交叉口几何尺寸、交通流量、车速、大小车比例、交通冲突以及信号配时数据;
然后,根据交叉口航拍图或者CAD图建立交叉口仿真模型,并定义仿真模型中交通流量、大小车比例、车辆速度以及交通信号配时信息。
3.根据权利要求1所述的交通冲突仿真两阶段参数标定方法,其特征是所述步骤2)中,初次参数校准采用软件默认值,若默认值无法满足研究需求,根据经验和现有参数来设定参数的变化区间和变化幅度,重复试验进行参数校准。
4.根据权利要求1所述的交通冲突仿真两阶段参数标定方法,其特征是所述步骤3)中,能够反映交叉口交通运行状况的检验指标,包括车头时距和排队长度。
CN2011103942487A 2011-12-02 2011-12-02 交通冲突仿真两阶段参数标定方法 Active CN102521438B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2011103942487A CN102521438B (zh) 2011-12-02 2011-12-02 交通冲突仿真两阶段参数标定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2011103942487A CN102521438B (zh) 2011-12-02 2011-12-02 交通冲突仿真两阶段参数标定方法

Publications (2)

Publication Number Publication Date
CN102521438A CN102521438A (zh) 2012-06-27
CN102521438B true CN102521438B (zh) 2013-12-04

Family

ID=46292351

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2011103942487A Active CN102521438B (zh) 2011-12-02 2011-12-02 交通冲突仿真两阶段参数标定方法

Country Status (1)

Country Link
CN (1) CN102521438B (zh)

Families Citing this family (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102915636A (zh) * 2012-10-19 2013-02-06 北京工业大学 交通信号机控制策略综合性能有效性测试系统及方法
CN103077275B (zh) * 2013-01-06 2016-04-13 东南大学 高速公路匝道仿真模型的参数标定方法
CN103077279B (zh) * 2013-01-10 2016-06-29 东南大学 一种城市交通环境下平面交叉口安全水平的评价方法
CN103106811A (zh) * 2013-01-15 2013-05-15 东南大学 一种基于两车碰撞时间的机动车有效交通冲突识别方法
CN103761138B (zh) * 2014-01-16 2017-01-18 昆明理工大学 一种交通仿真软件的参数校正方法
CN103942398B (zh) * 2014-05-05 2016-10-19 江苏物联网研究发展中心 基于遗传算法及广义回归神经网络的交通仿真校正方法
CN104820763B (zh) * 2015-05-25 2017-12-15 西华大学 一种基于vissim的交通事故三维仿真方法
CN104900064B (zh) * 2015-06-19 2018-01-09 南京全司达交通科技有限公司 一种获取高速公路优选交通组织方式的方法及系统
CN105243876B (zh) * 2015-11-06 2017-05-24 东南大学 一种用于互通立交的交通冲突严重性分析方法
CN105513357B (zh) * 2015-12-30 2017-11-03 银江股份有限公司 一种基于微波数据的交通流基本图参数标定方法
CN107085382A (zh) * 2016-02-15 2017-08-22 北京福田康明斯发动机有限公司 一种车辆参数标定设备和方法
CN105930565B (zh) * 2016-04-13 2019-11-19 中山大学 基于线性策略的交叉熵算法的交通仿真模型参数标定方法
CN106529076B (zh) * 2016-11-28 2019-07-30 东南大学 高速公路交通安全仿真分析的两阶段参数标定方法
CN107248290B (zh) * 2017-07-14 2019-07-12 东南大学 基于车辆视频识别的混合式主线收费站交通冲突评价方法
CN107766636A (zh) * 2017-10-12 2018-03-06 东南大学 一种基于极值理论与微观仿真的城市交叉口安全评价方法
CN109087534A (zh) * 2018-10-09 2018-12-25 王业宝 一种基于车辆行驶轨迹的交通冲突检测方法
CN109544934A (zh) * 2018-12-19 2019-03-29 同济大学 基于城市交叉口混合交通流三维轨迹的效率安全监控系统
CN110992675B (zh) * 2019-05-10 2021-06-01 北京五一视界数字孪生科技股份有限公司 交叉口交通仿真方法、装置、介质及电子设备
CN112150046A (zh) * 2020-11-06 2020-12-29 上海交通大学 一种道路交叉口安全风险指数计算方法
CN112487617B (zh) * 2020-11-19 2024-03-08 平安科技(深圳)有限公司 基于碰撞模型的风险预防方法、装置、设备及存储介质
CN112580954B (zh) * 2020-12-11 2021-09-07 东南大学 一种基于交通冲突极值分布的仿真模型标定方法
CN112598236B (zh) * 2020-12-14 2021-09-07 东南大学 一种基于交通冲突的事前事后安全评价方法
CN113094875B (zh) * 2021-03-16 2024-04-09 东南大学 一种城市快速路交织区微观交通仿真系统标定方法及装置
CN114611932B (zh) * 2022-03-10 2023-09-05 哈尔滨工业大学 基于仿真和变道冲突的入口加速车道服务水平分级方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1844580A (zh) * 2006-03-21 2006-10-11 东南大学 确定城市单向交通配对道路最佳间距的方法

Also Published As

Publication number Publication date
CN102521438A (zh) 2012-06-27

Similar Documents

Publication Publication Date Title
CN102521438B (zh) 交通冲突仿真两阶段参数标定方法
CN108573601B (zh) 一种基于wim数据的交通安全风险场构建方法
CN103439070B (zh) 一种桥梁长期挠度效应的分离方法
CN101364344B (zh) 一种基于压力测试的路网极限容量确定方法
CN105741551A (zh) 一种无信号交叉口交通安全评价方法
CN103077275B (zh) 高速公路匝道仿真模型的参数标定方法
CN103761138A (zh) 一种交通仿真软件的参数校正方法
CN103093623B (zh) 一种平面信号交叉口直左冲突数的预测方法
CN105930614A (zh) 一种针对可变限速控制的元胞传输模型参数标定与验证方法
CN111695241A (zh) 一种基于vissim仿真的左入匝道合流区加速车道长度确定方法
CN104282149A (zh) 一种基于交通状态精度指标评价的路网浮动车配置方法
Ge et al. Construction and simulation of rear-end conflicts recognition model based on improved TTC algorithm
CN105824987A (zh) 一种基于遗传算法的风场特征统计分布模型建立方法
Luo et al. Research of vehicle rear-end collision model considering multiple factors
CN107092794A (zh) 用于台风混合气候地区极值风速估算的分析方法及装置
CN108898857A (zh) 一种考虑安全可靠性的交叉口机动车绿灯间隔设置方法
CN105046958A (zh) 一种高速公路交通信息采集点非等距优化布局方法
Liu et al. Digital twin of highway entrances and exits: A traffic risk identification method
CN105912775A (zh) 桥梁动态称重系统车辆轴重数据的多模态建模方法
Chang et al. A multi-stage linear prediction model for the irregularity of the longitudinal level over unit railway sections
CN105095649B (zh) 建筑物倾斜角度的数据分析方法和系统
CN109783893B (zh) 一种考虑路段上下游线形设计参数的连续建模方法
Li et al. VISSIM parameter calibration based on traffic characteristics distribution at signalized intersections
CN107945515A (zh) 基于车辆运行速度的道路交通能见度分级方法
CN105043786A (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
C14 Grant of patent or utility model
GR01 Patent grant