摘要
为了提高恶劣天气下高速公路交通安全水平,实施更加科学有效的高速公路车速管理策略,提出了基于安全风险的恶劣天气下高速公路建议车速确定方法。综合考虑交通安全风险影响因素,基于贝叶斯Logistic回归建立交通风险评估模型,利用模糊C⁃均值聚类算法进行交通风险等级划分并得到了风险阈值。在此基础上提出了不同天气条件下安全车速的两步计算方法,即先根据交通风险分布,在风险小于阈值的交通流状态中选取车辆运行速度的第85 %位车速作为初定安全行驶车速,后利用停车视距原理进行修正,最终得到恶劣天气下高速公路建议车速。并以G15沈海高速(上海段)为例,展示了计算方法,给出了雾天不同能见度条件下高速公路建议车速取值。
高速公路在带来巨大的社会经济效益的同时,由恶劣天气引发的重、特大高速公路交通事故也日渐频发。恶劣天气会对交通流运行特征产生影响,进而也会带来交通安全隐患。在高速公路交通安全管理中,速度管理是重要的管理手段。
绝大多数国家对车辆在高速公路上行驶都有最高车速限制,如意大利最高限速为140 km·
近年来,随着主动式交通管理与控制的理念愈加成熟,可变限速(也叫速度协调控制或动态限速)的相关研究逐渐增多,部分学者分析了可变限速策略对于降低路段交通事故风险的效果,如Sah
本文首先给出了恶劣天气下高速公路建议车速确定方法的思路和流程,随后明确了恶劣天气下高速公路交通安全风险的评估与等级划分方法,在此基础上提出了恶劣天气下高速公路建议车速计算方法,并通过实例分析展示了方法的应用,给出了G15沈海高速(上海段)雾天不同能见度条件下高速公路建议车速取值。
高速公路建议车速是指在特定的道路、环境、交通组成的综合情境下,驾驶员受到道路条件、车辆稳定性、驾驶员处理能力等因素的影响,能够保持安全行驶所建议采取的车速。通常来讲,驾驶员会在行驶过程中依据道路条件、车流状况、所驾驶车辆性能等因素经综合考虑后形成心理上自认为的安全行驶车速,即主观期望车
交通安全风险是指,在某些随机因素可能引起的实际情况与预定目标产生偏离的交通安全状况下,发生交通安全事故的可能性或概

图 1 恶劣天气下高速公路建议车速制定流程
Fig. 1 Flowchart of determination of freeway recommended speed in adverse weather
恶劣天气下高速公路建议车速制定流程中,关键是通过建立高速公路交通安全风险模型来进行安全风险评估,并利用交通风险值计算安全行驶车速。
传统集计交通事故模型是一种静态、被动的“事后”分析方法,难以反映交通流动态特征对事故风险的影响。近年来,交通安全学界开始利用高精度交通流数据研究实时事故风险模
交通安全风险评估模型中,响应变量事故发生与否是一个二分变量,Logistic回归模型在处理该问题时被广泛应用。不同于传统基于最大似然估计求解的Logistic回归模型,贝叶斯Logistic回归给出模型参数的后验概率分布,提供了灵活的框架,通过设定先验概率分布,可以实现模型的空间移植。另外,贝叶斯Logistic回归能够有效避免小样本、多变量数据引起的过拟合问
设y为一个二分类因变量,事故发生时y=1,无事故发生时y=0。贝叶斯Logistic回归的公式如下:
(1) |
(2) |
式(1)~(2)中:为交通事故发生的概率;为效用函数;为样本i中变量j的值;为回归截距;为解释变量j的回归系数。
模型的似然函数表达式如下:
(3) |
式中:N为样本数。
模型中的参数均采用无信息先验概率分布:
式中:为期望;为方差。
通常具有大方差的先验分布可以代表无信息先验概率分布,令=0,=10 000。
根据贝叶斯定理,参数的后验联合概率密度分布正比例与似然函数和先验概率分布的乘积,即:
(4) |
将高速公路交通风险分为安全状态、可允许风险、中度风险、重大风险、不可接受风险5个等级,具体风险等级划分的标准如
将恶劣天气下高速公路交通风险评估模型计算得到的交通风险值进行聚类,对交通风险进行等级划分。由于交通风险不同等级之间的界限有一定的模糊性,本文使用模糊C-均值聚类算
模糊C-均值聚类算法基于样本与c个聚类中心间的加权相似性测度,对目标函数进行迭代最小化,以确定其最佳的类别。目标函数定义如下:
(5) |
且满足条件:
(6) |
式中:X 为聚类样本集合;n为聚类空间的样本个数;V 为c个聚类中心;c为聚类的类别数;为与之间的归一化距离;U 为维的矩阵;为第k个样本对i类的隶属度值;m为初始幂指数参数。
模糊C-均值聚类法计算步骤如下:
第1步:根据样本划分类的个数c,幂指数和初始隶属度矩阵,取上的均匀分布随机数来确定。令表示第1步迭代。
第2步:计算第步的聚类中心:
(7) |
第3步:修正隶属度矩阵,计算目标函数值。
(8) |
(9) |
式中:。
第4步:对给定的隶属度终止容限,或目标函数终止容限,或对于最大迭代步长,当,或当,有或时,迭代停止,否则,然后重复第2步,第3步。
经过上述的循环迭代之后,当目标函数达到最小值时,根据最终的隶属度矩阵U中元素的取值确定所有样本的归属,当 时,可将样本归为第类。
以G15沈海高速(上海段)数据为样本,进行案例分析。采集了2013年1月1日至2013年12月31日时段沈海高速(上海段)部分路段的交通流数据、天气数据及道路特征数据。
采用非配对病例-对照研究方法提取正常情况下的交通流数据、天气数据和道路特征数据,并进行数据匹配。其中,选取天气状况为雾天和雨天的事故数据,共86个;对于样本中的每一个事故数据,按照事故和非事故数据1:4的比
基于样本数据集,利用随机森林算法,筛选模型变量,最终选入模型中的变量如
利用R软件的rstanarm包,实现贝叶斯Logistic模型的建立,并通过马尔可夫链蒙特卡罗方法(Markov Chain Monte Carlo,MCMC)方法计算各个回归系数的后验概率分布。模型标定结果如
根据模型标定结果,计算交通风险值。通过模糊C-均值聚类算法,根据高速公路交通风险等级划分标准确定聚类数为5,将交通风险值进行聚类分析,结果如

图 2 交通风险等级聚类结果
Fig. 2 Traffic risk level clustering result
依据交通风险值的各聚类类别的最大值、最小值和风险等级划分的标准,确定恶劣天气下高速公路交通不同预警风险等级所对应的交通风险值范围,如
由此得到风险阈值为0.2,即选取风险值小于0.2的交通流状态下的车辆运行速度作为安全行驶车速的取值范围。
对交通风险进行控制需要从有潜在交通发生的时候进行干预,降低高速公路交通运行风险。选择安全风险等级为安全状态和可允许风险两个等级的交通流状态,将该状态下车辆行驶速度作为安全行驶车速的取值范围。统计交通风险值小于阈值的所有车辆运行速度,绘制累积车速分布曲线。
累积车速分布曲线中,累积频率为85 %时对应的车速,即85 %位车速,是最高车速限制的重要依
因此,基于交通风险分布的高速公路安全行驶车速计算方法为,根据不同恶劣天气条件下的交通安全风险分布,选取交通风险值小于风险阈值的交通流状态作为安全车速的确定范围,绘制该范围内的累积车速分布曲线,将车辆运行速度的第85 %位车速作为相应天气条件下的安全行驶车速。
停车视距指的是同一车道上,车辆行驶时遇到前方障碍物而必须采取制动停车时所需要最短行车距离。停车视距可分解为反应距离、制动距离和安全距离三部分,如

图 3 跟车状态下的停车过程
Fig. 3 Stopping process when cars follow
一般情况下,驾驶员发现前面车辆时,前车处于制动状态,此时后车停车所需的安全距离应满足:
(10) |
式中:为后方车辆驾驶员反应时间内的车辆行驶距离,m;为后方车辆制动时间内的行驶距,m;为安全距离,一般取值为5~10 m,为保障恶劣天气下的行车安全,取值为20 m;为路段的可视距离,m。
在恶劣天气下,驾驶员有效视距和路面附着系数会发生变化,为了保障车辆的安全行驶,应该考虑最不利的情况,即由于车辆故障、轮胎损坏、抛锚、货物洒落及事故等原因,前方物体的速度为零,车流中出现严重的速度差,后车必须进行紧急制动。一般情况下驾驶员的反应时间为0.5~1.7
将恶劣天气下的能见度作为可视距离,计算得到基于停车视距的高速公路安全车速,如式(11)所示:
(11) |
式中:为安全车速,km·
对比两种方法等到的不同雾天等级条件下高速公路安全行驶车速,综合考虑不同恶劣天气等级下的交通风险分布和停车视距理论,给出不同恶劣天气等级下高速公路建议车速。为保障恶劣天气下高速公路行车安全,选择更低的安全车速作为该天气等级下的建议车速。
以G15沈海高速(上海段)为研究对象,采集了2013年1月1日至2013年12月31日的交通流数据及天气数据,并进行数据匹配。为减少数据噪声,将交通流数据以5 min进行集计。依据恶劣天气下高速公路交通风险评估模型计算雾天条件下的交通风险分布,基于交通风险分布计算高速公路安全行驶车速,并通过停车视距理论进行修正,最终得到雾天条件下高速公路建议车速。
参照《高速公路交通气象等级》(QX/T 111—2010

图 4 不同能见度等级下交通风险分布
Fig. 4 Traffic risk distribution in different visibilities
在不同雾天等级条件下,选择高速公路交通风险值小于等于0.2的交通流状态作为安全车速的确定范围,绘制该状态下车速累积频率分布曲线并计算车辆运行速度的第85 %位车速,如

图 5 不同能见度等级下速度累积频率分布曲线
Fig. 5 Accumulated frequency distribution of speed in different visibilities
依照速度累积频率分布曲线图,计算得到不同能见度等级下车辆运行速度的第85 %位车速,为方便高速公路驾驶人员接受和管理部门发布信息,将第85 %位车速取整后给出安全行驶车速,如
雾天天气时,由于雾水落于路面,使得路面潮湿,附着系数降低,取潮湿路面的附着系数。依据《公路工程技术标准》(JTG B01—2014
根据基于停车视距的高速公路安全行驶车速模型,通过对能见度及道路纵坡取不同数值,计算得到不同能见度等级下安全车速。为方便高速公路驾驶人员接受和管理部门发布信息,将车速取整后给出不同等级下安全行驶车速,如
当能见度在200~500 m范围内,基于停车视距理论,无需对高速公路上行驶车辆进行速度控制,基于该天气条件下的交通风险分布,高速公路建议安全车速为85 km·
从交通安全风险角度出发,提出了恶劣天气下高速公路建议车速的确定方法。通过建立交通安全风险评估模型,划分交通风险预警等级,根据不同天气等级下的交通风险分布,按照安全交通风险预警等级确定不同天气等级下的初定安全行驶车速,并利用停车视距原理对基于安全风险计算得到的安全行驶车速进行修正,给出恶劣天气下高速公路建议车速。以G15沈海高速(上海段)为例,展示了计算方法的应用,并给出了雾天不同能见度条件下高速公路建议车速取值。
论文综合考虑了影响车速的多种因素,尤其考虑了交通流条件和天气条件等可变因素对速度的影响,从而可以实现动态地确定建议车速。引入了交通安全风险作为建议车速制定的主要依据,从影响交通安全风险的因素出发,分析了交通安全风险和建议车速的关系。利用基于停车视距原理对安全行驶车速进行修正,保证了恶劣天气下高速公路建议车速满足停车视距的安全要求。
由于数据、时间等限制,基于安全风险的恶劣天气下高速公路建议车速的实施效果以及具体的实施细节还需要进一步研究,后续可进一步探索基于交通风险评估的恶劣天气下高速公路交通安全主动控制策略。
参考文献
中华人民共和国全国人民代表大会常务委员会. 中华人民共和国道路交通安全法[J]. 中华人民共和国全国人民代表大会常务委员会公报, 2003(6): 536. [百度学术]
Standing Committee of the National People’s Congress of the People’s Republic of China. Law of the People’s Republic of China on road traffic safety[J]. Communiqué of the Standing Committee of the National People’s Congress of the People’s Republic of China, 2003(6): 536. [百度学术]
Committee for Guidance on Setting and Enforcing Speed Limits National Research Council. Managing speed: Review of current practice for setting and enforcing speed limits[M]. Washington D C: National Academy Press, 1998. [百度学术]
中华人民共和国交通运输部. 公路路线设计规范:JTG D20—2017[S]. 北京:人民交通出版社,2017. [百度学术]
Ministry of Transport of the People’s Republic of China. Design specification for highway alignment:JTG D20—2017[S]. Beijing:China Communications Press,2017. [百度学术]
顾叶华, 王伟. 高速公路限速值调整交通工程论证指标研究[J]. 智能城市, 2019(9): 85. [百度学术]
GU Yehua, WANG Wei. Study on the demonstration index of expressway limited speed adjustment traffic engineering[J]. Intelligent City, 2019(9): 85. [百度学术]
SAHA P, YOUNG R K. Weather-based safety analysis for the effectiveness of rural vsl corridors[C]//93rd Annual Meeting of the Transportation Research Board. Washington D C: [s.n.], 2014. [百度学术]
KWON E, BRANNAN D, SHOUMAN K, et al. Development and field evaluation of variable advisory speed limit system for work zones[J]. Transportation Research Record, 2007, 2015(1): 12. [百度学术]
WALLER S T, NG M W, Ferguson E, et al. Speed harmonization and peak-period shoulder use to manage urban freeway congestion[R]. Austin: University of Texas at Austin. Center for Transportation Research, 2009. [百度学术]
LI D, RANJITKAR P, CEDER A. A logic tree based algorithm for variable speed limit controllers to manage recurrently congested bottlenecks[C]//93rd Annual Meeting of the Transportation Research Board. Washington D C: [s.n.], 2014:1-14. [百度学术]
王艳丽, 李晓庆, 王忠宇, 等.面向出口匝道拥挤的快速路速度协调控制模型[J]. 同济大学学报(自然科学版), 2018, 46(7): 49. [百度学术]
WANG Yanli, LI Xiaoqing, WANG Zhongyu, et al. Speed harmonization model for off-ramp bottlenecks on urban expressway[J]. Journal of Tongji University (Natural Science), 2018, 46(7): 49. [百度学术]
LUOMA J, RÄMÄ P, PENTTINEN M, et al. Effects of variable message signs for slippery road conditions on reported driver behaviour[J]. Transportation Research Part F: Traffic Psychology and Behaviour, 2000, 3(2): 75. [百度学术]
郑安文. 期望车速的意义及其影响因素分析[J]. 武汉科技大学学报(自然科学版), 2005,28(1):61. [百度学术]
ZHENG Anwen. Expectation speed: Significance and influential factors[J]. Journal of Wuhan University of Science and Technology(Natural Science Edition), 2005, 28(1): 61. [百度学术]
涂辉招, 牟涛, 毕玉峰, 等. 基于安全车速分析的高速公路安全评价方法[J]. 内蒙古公路与运输, 2015(6): 1. [百度学术]
TU Huizhao, MOU Tao , BI Yufeng, et al. Highway safety evaluation method based on safe vehicle speed analysis[J]. Highways & Transportation in Inner Mongolia, 2015(6): 1. [百度学术]
赵学刚. 区域路网交通安全风险动态预警关键技术研究[D]. 西安: 长安大学, 2010. [百度学术]
ZHAO Xuegang. Study on key technology of dynamic early warning for traffic safety risk of regional road network[D]. Xi’an: Chang’an University, 2010. [百度学术]
ABDEL-ATY M, UDDIN N, PANDE A. Split models for predicting multivehicle crashes during high-speed and low-speed operating conditions on freeways[J]. Transportation Research Record, 2005, 1908(1): 51. [百度学术]
杨奎, 余荣杰, 王雪松. 基于车道集计交通流数据的事故风险评估分析[J]. 同济大学学报(自然科学版), 2016, 44(10): 1567. [百度学术]
YANG Kui, YU Rongjie, WANG Xuesong. Application of aggregated lane traffic data from dual-loop detector to crash risk evaluation[J]. Journal of Tongji University (Natural Science), 2016, 44(10): 1567. [百度学术]
XU C, WANG W, LIU P, et al. Using the Bayesian updating approach to improve the spatial and temporal transferability of real-time crash risk prediction models[J]. Transportation Research Part C: Emerging Technologies, 2014, 38: 167. [百度学术]
BEZDEK J C, EHRLICH R, FULL W. FCM: The fuzzy c -means clustering algorithm[J]. Computers & Geosciences, 1984, 10(2): 191. [百度学术]
YU R, ABDEL-ATY M. Utilizing support vector machine in real-time crash risk evaluation[J]. Accident Analysis & Prevention, 2013, 51: 252. [百度学术]
中国公路学会. 交通工程手册[M]. 北京:人民交通出版社. 2001. [百度学术]
China Highway & Transportation Society. Traffic engineering manual [M]. Beijing: People’s Communications Press. 2001. [百度学术]
AGENT K, PIGMAN J, WEBER J. Evaluation of speed limits in Kentucky[J]. Transportation Research Record, 1998, 1640(1): 57. [百度学术]
FITZPATRICK K, BLASCHKE J D, SHAMBURGER C B, et al. Compatibility of design speed, operating speed , and posted speed[R]. Austin: Texas Department of Transportation, College Station, 1995. [百度学术]
吴斌, 朱西产, 沈剑平. 基于自然驾驶数据的驾驶员紧急制动行为特征[J].同济大学学报(自然科学版), 2018, 46(11), 1514. [百度学术]
WU Bin, ZHU Xichan, SHEN Jianping. Driver emergency braking behavior based on naturalistic driving data[J]. Journal of Tongji University (Natural Science), 2018, 46(11), 1514. [百度学术]