2013年全国大学生数学建模竞赛A题:车辆排队长度与事故持续时间、道路实际通行能力、路段上流流量间的关系
- 格式:doc
- 大小:113.00 KB
- 文档页数:12
承诺书我们仔细阅读了《全国大学生数学建模竞赛章程》和《全国大学生数学建模竞赛参赛规则》(以下简称为“竞赛章程和参赛规则”,可从全国大学生数学建模竞赛网站下载)。
我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题。
我们知道,抄袭别人的成果是违反竞赛章程和参赛规则的,如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺,严格遵守竞赛章程和参赛规则,以保证竞赛的公正、公平性。
如有违反竞赛章程和参赛规则的行为,我们将受到严肃处理。
我们授权全国大学生数学建模竞赛组委会,可将我们的论文以任何形式进行公开展示(包括进行网上公示,在书籍、期刊和其他媒体进行正式或非正式发表等)。
我们参赛选择的题号是(从A/B/C/D中选择一项填写):A我们的参赛报名号为(如果赛区设置报名号的话):A1106所属学校(请填写完整的全名):龙岩学院参赛队员(打印并签名) :指导教师或指导教师组负责人(打印并签名):指导组(论文纸质版与电子版中的以上信息必须一致,只是电子版中无需签名。
以上内容请仔细核对,提交后将不再允许做任何修改。
如填写错误,论文可能被取消评奖资格。
)日期:2013 年 9 月 16 日赛区评阅编号(由赛区组委会评阅前进行编号):编号专用页赛区评阅编号(由赛区组委会评阅前进行编号):全国统一编号(由赛区组委会送交全国前编号):全国评阅编号(由全国组委会评阅前进行编号):车道占用后通行能力的优化方案摘要本文就所给的视频资料针对交通事故占道对实际通行能力的影响问题,先后建立了无向赋权图、随机排队理论、集结波与消散波模型、车流波动模型。
对各个问题进行分析、求解。
对于问题(1),结合所给的视频1,对其中的数据进行多次、合理采样,并用图表表示事故所处横断面实际通行能力的变化过程。
2013高教社杯全国大学生数学建模竞赛编号专用页赛区评阅编号(由赛区组委会评阅前进行编号):全国统一编号(由赛区组委会送交全国前编号):全国评阅编号(由全国组委会评阅前进行编号):车道被占用对城市道路通行能力的影响摘要本文主要研究车道被占用对城市道路通行能力的影响并建立了相应的数学模型。
针对问题一,考虑到交通信号灯的周期,我们选择1分钟为周期,结合不同车辆的标准车当量的折算系数,求出每个采样点的交通量,通过MATLAB作图,从定性方面对道路通行能力进行分析,然后通过基本通行能力和4个修正系数建立动态通行能力的模型。
图像显示,事故发生后(采样点5附近),实际通行能力下降至一个较低水平,并且横断面处的实际能力变化过程呈先下后上的波形变化,在事故解决(第20个采样点)以后,由图像看出实际通行能力持续上升。
针对问题二,利用问题一建立的模型,结合视频二,比较交通事故所占不同车道时横断面的实际通行能力,可以发现二者实际通行能力变化趋势大致相同,但视频二实际通行能力大于视频一实际通行能力。
可见占用车流量大的车道使道路通行能力降低更多。
针对问题三,首先我们建立单车道排队车辆数目的积分模型,单个车道的滞留车辆为上游车流量和实际通行能力的差值。
我们以30s为一个时间段,对视频一中的车流量进行统计,得到横截面处每个监测段的实际通行能力。
本题要求考虑三车道,总体排队长度不容易通过积分模型确定,所以我们将队列长度问题转化为车辆数目问题,通过视频资料统计120米对应24辆车,据此关系转换,从而得到车辆排队长度与事故横断面实际通行能力、事故持续时间和上游车流量的关系。
针对问题四,在对问题3研究的基础上,根据问题3建立的数学模型,建立起某一段时间间隔车辆排队的长度,然后,通过求得的关系得到当排队长度为140m的时候所对应的时间段,由于每段时间间隔设为30s,因此,可以求得排队长度到达上游时用的时间为347.7273s。
关键词:交通事故车道占用通行能力排队论一、问题的重述车道被占用是指因交通事故、路边停车、占道施工等因素,导致车道或道路横断面通行能力在单位时间内降低的现象。
车道被占用对城市道路通行能力的影响影响道路通行能力的主要因素有道路状况、车辆性能、交通条件、交通管理、环境、驾驶技术和气候等条件。
道路条件是指道路的几何线形组成,如车道宽度、侧向净空、路面性质和状况、平纵线形组成、实际能保证的视距长度、纵坡的大小和坡长等。
车辆性能是指车辆行驶的动力性能,如减速、加速、制动、爬坡能力等。
交通条件是指交通流中车辆组成、车道分布、交通量的变化、超车及转移车道等运行情况的改变。
环境是指街道与道路所处的环境、景观、地貌、自然状况、沿途的街道状况、公共汽车停站布置和数量、单位长度的交叉数量及行人过街道等情况。
气候因素是指气温的高低、风力大小、雨雪状况!公路通行能力的计算方法公路通行能力的计算方法(一)、无平交路段通行能力(1)基本通行能力一般路段是指不受信号、暂停标志、铁公路口等外界因素的中断,保证大体连续的交通流的公路部分。
多车道公路的基本通行能力是以高速公路上观测到的最大交通量为基准确定的。
根据观测结果,城市快速路比城际间高速公路的值来得大一些,在大体接近城市快速路最大交通量处确定了多车道公路的基本通行能力为每车道2200pcu/h。
往返2车道公路的基本通行能力用往返合计值表示。
其理由为往返2车道公路通常不进行往返车道的分离,以供对面车辆超车用,这种方法是比较现实的。
实际上,在往返2车道公路上发生超车时的最大交通量的观测数据非常少,在美国《公路通行能力手册》中写明往返2车道公路的基本通行能力大约为多车道公路中2车道基本通行能力的二分之一,并确定为2500pcu/h。
另外,与多车道公路相同,对单向通行公路,把其基本通行能力定为每车道2200pcu/h。
(2)可能通行能力可能通行能力是用基本通行能力乘以公路的几何结构、交通条件对应的各种补偿系数求出的。
亦即C= CB*γL*γC*γI*……(2.1)式中,C:可能通行能力;CB:基本通行能力;γLγCγI:各种补偿系数。
就多车道公路而言,先用(2.1)式求出每车道的可能通行能力,然后乘以车道数求出公路截面的可能通行能力。
2013高教社杯全国大学生数学建模竞赛B题车道被占用对城市道路通行能力的影响摘要车道被占用是指因交通事故、路边停车、占道施工等因素,导致车道或道路横断面通行能力在单位时间内降低的现象。
由于城市道路具有交通流密度大、连续性强等特点,一条车道被占用,也可能降低路段所有车道的通行能力,即使时间短,也可能引起车辆排队,出现交通阻塞。
如处理不当,甚至出现区域性拥堵。
对于问题一,本文提高结果的精准度,结合两种方法进行研究,且两种方法的结果十分吻合。
由于实际通行能力是建立在基本通行能力和可能通行能力之上的,所以在求解实际通行能力之前,需要算出基本通行能力和可能通行能力,针对问题一创建了一张流程图,并借助软件加以拟合。
对实际通行能力计算,得出实际通行能力的变化过程,根据GREENSHIELD K-V线性算法得出道路越堵,车速越慢,则实际通行能力就越差,反之就会较好。
对于问题二,因为所占的车道不同,并且给的条件中有说明左转车流比例和右转车流比例不同,那只需验证两者是否存在显著性差异,运用配对样本t检验的方法就是要先满足这一方法的两个前提条件,首先必须验证是否满足正态分布,经过SPSS软件的验证可以得出符合正态分布。
然后再进行配对,从配对的结果中可以看出存在显著性差异,再结合左右转的车流量比例,更加可以看出存在显著性差异。
对于问题三,主要是对所推出来的回归方程的判断和分析因变量和各因子之间的关系,在本问中要先求出排队长度,排队长度是根据堵塞密度,进出车辆数之间的差值来求解,再根据最小二乘法来判断所假设的这一模型是否符合多元线性回归关系,本问中得出符合多元线性回归关系。
再在排队长度和最小二乘法的基础之上,运用SPSS软件,在进行结果分析时得出实际通行能力对于排队长度没有影响,所以可以剔除,而事故持续时间和上游车流量对排队长度都有明显的影响,然后得出他们的相关系数,求出最后的相关方程式。
对于问题四,题目中给出了事故发生点到上游路口的距离为140米,并且上游车流量为1500pcu/h,结合视频1中多次出现的120米这一个顶点,推算出120米内大概最大的堵塞车流量,然后按比例分配推算出140米的最大堵塞车流量,视频1中的可以通过加权平均来求出平均的实际通行能力,则事故持续时间就是要靠140米的最大堵塞车流量和平均实际通行能力来计算,最后得出事故持续时间为2.37min。
1 交通事故影响时间分析由于从交通事故发生到检测到事故、接警、事故现场勘测、处理、清理事故现场恢复交通,以及恢复交通后车辆排队不再增加都需要一定的时间。
这部分时间主要由三部分构成: 第一部分是事故发生到警察到达现场的时间1T ; 第二部分是交通事故现场处理时间2T ,由现场测、处理到事故族除、恢复交通;第三部分是交通事故持续影响时间T3,部分时间从恢复事故场交通开始,到事故上游车辆排队不再增加,即排队开始减弱。
在T1内,事故现场保持原状,没有进行处理,这里分两种情况考虑: ( 1) 当交通事故占部分车道时,这时事故点的剩余通行能力Qs ≠0,交通事故越严重,则相应Qs 越小。
若事故点上游的交通需求Q < Qs ,则车辆以较低的速度通过事故点,上游不会形成车辆拥挤排队; 若Q > Qs ,则交通流可按事故点的剩余断面通行能力通过事故点,超过该通行能力的车流在事故点上游排队。
( 2) 当交通事故十分严重时,事故点的剩余通行能力Qs = 0,造成事发路段断流,事故点上游车辆排队,发生交通拥挤堵塞,进而排队一直向上游延伸。
在T2内,确认交通事故发生后,相关部门到现场处理异常事件,在此过程中,事故点交通可能会受到进一步影响,事故断面通行能力也随之发生变化[5],一般会变小,甚至变为0( 全封闭处理) ,视事件处理具体情况而定,事发点上游交通处于严重拥挤状态,车辆排队增加。
由于在交通事故接警时间T1和处理时间T2阶段事故点上游交通车辆产生排队,若没有车辆排队,则T3 = 0; 若有车辆排队,则当事故处理完毕、道路恢复交通时,排队车辆开始消散。
交通事故持续影响时间T3是事故处理完毕、道路恢复交通至车辆排队不再增加这段时间,即交通流消散波从车辆排队队列的头部传到尾部这段时间。
2、事故路段车辆排队长度分析如下图图发生交通事故的路段该事故路段长度为L( m) ,单方向车道数为n,单方向车道宽度为D( m) ,在道路上t = 0 时刻发生了一起交通事故,事故车辆占用道路宽度为b( m) ,长度为a( m) ,事故点上游路段长度为L'。
2013高教社杯全国大学生数学建模竞赛A题2013高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了《全国大学生数学建模竞赛章程》和《全国大学生数学建模竞赛参赛规则》。
我们完全明白,在竞赛开始后参赛队员不能以任何方式与队外的任何人研究、讨论与赛题有关的问题。
我们知道,抄袭别人的成果是违反竞赛章程和参赛规则的,如果引用别人的成果或其他公开的资料,必须按照规定的面的车辆数。
实际通行车流量的采集与处理视频1中出现车辆多种多样,要统计车流量数据,需先统一车流标准,把视频中出现的车辆进行折算,以小轿车做为标准,对各个型号车辆进行折算[2],折算系数如表1所示。
表1 车辆折算系数附件中出现汽车小轿车中型车大客车车辆折算系数在事故发生前,道路的通行能力足以应对上游车流量,当发生事故时,事故点上游共有10辆小轿车与5辆大客车,车流量为20pcu。
之后一分钟(16:42:32-16:43:32),上游又有车流量21pcu,但只通过了21pcu,说明造成了交通拥堵和排队情况。
“附件5”可知,相位时间为30s,红灯时间为30s,即60s为一个周期,进行统计时间周期也为60s,不会造成因交通灯引起的误差。
实际通行流量是指折算后通过事故横断面的车流,上游车流量是指折算后从各个路口驶入事故横断面的车流。
对附件1中事故横断面处的车流量进行统计,得出实际通行车流量情况,并统计横断面上游的车流量,在统计过程中发现视频并不是完全连续的,例如在16:49:40时出现了突变,直接到16:50:04,跳跃间隔为24s,但于堵车情况较重,可以根据车流量守恒原则和车辆追踪,统计出通过横断面处的车流量及上游车流量。
但16:56:04等时间,跳跃时间较长,近2分钟,无法精确统计,如表2处“空缺”所示。
在17:00:07到17:01:20时视频发生跳变,在此期间事故车辆驶离道路,之后为事故恢复时间。
为了描述事故发生开始到车辆离开车道全程的实际通行能力变化情况,将视频中空缺数据通过灰色预测(程序见附录)进行填补,结果如表2所示。
道路上不断增加的交通流经常导致拥挤。
拥挤产生延误、降低流率、带来燃油损耗和负面的环境影响。
为了提高道路系统的效率,国内外许多研究者一直致力于车流运行模型的研究。
Daganzo[1]提出了一种和流体力学LWR 模型相一致的元胞传输模型,这种模型能用来模拟和预测交通流的时空演化,包括暂时的现象,如排队的形成、传播、和消散。
Heydecker 和Addison[2]通过研究车速和密度的因果关系分析和模拟了在变化的车速限制下的交通流。
Jennifer 和Sallissou[3]提出了一种混合宏观模型有效地描述了路网的交通流。
然而,拥挤也会由交通异常事件引起。
交通异常事件定义为影响道路通行能力的意外事件[4],如交通事故、车辆抛锚、落物、短期施工等,从广义角度看,还应包括恶劣天气与特殊勤务等。
异常事件往往造成局部车道阻塞或关闭,形成交通瓶颈,引起偶发性拥挤,这已经逐渐成为高速道路交通拥挤的主要原因[5],越来越多地受到研究者们的重视。
例如M. Baykal-Gursoy[6]等人提出了成批服务受干扰下的稳态M/M/c 排队系统模拟了发生异常事件的道路路段的交通流。
Chung[7]依据韩国高速公路系统监测的准确记录的大型交通事故数据库提出了一种事故持续时间预测模型。
当然,这些研究最终都是为了帮助缓解异常事件引起的交通拥挤。
交通异常事件发生后,事发地段通行能力减小,当交通需求大于事发段剩余通行能力时,车辆排队,产生延误,行程时间增加[8],交通流量发生变化。
本文以高速公路基本路段发生交通事故为例,主要分析了交通事故发生后不同时间段内事故点及其上游下游路段交通流量的变化,用于以后进一步的突发事件下交通流预测工作。
1 交通事故影响时间分析由于从交通事故发生到检测到事故、接警、事故现场勘测、处理、清理事故现场恢复交通,以及恢复交通后车辆排队不再增加都需要一定的时间。
这部分时间主要由三部分构成: 第一部分是事故发生到警察到达现场的时间T1; 第二部分是交通事故现场处理时间T2,由现场勘测、处理到事故族除、恢复交通; 第三部分是交通事故持续影响时间T3,这部分时间从恢复事故现场交通开始,到事故上游车辆排队不再增加,即排队开始减弱[9]。
2018高教社杯全国大学生数学建模竞赛题目<请先阅读“全国大学生数学建模竞赛论文格式规范”)A题车道被占用对城市道路通行能力地影响车道被占用是指因交通事故、路边停车、占道施工等因素,导致车道或道路横断面通行能力在单位时间内降低地现象.因为城市道路具有交通流密度大、连续性强等特点,一条车道被占用,也可能降低路段所有车道地通行能力,即使时间短,也可能引起车辆排队,出现交通阻塞.如处理不当,甚至出现区域性拥堵.车道被占用地情况种类繁多、复杂,正确估算车道被占用对城市道路通行能力地影响程度,将为交通管理部门正确引导车辆行驶、审批占道施工、设计道路渠化方案、设置路边停车位和设置非港湾式公交车站等提供理论依据.视频1<附件1)和视频2<附件2)中地两个交通事故处于同一路段地同一横断面,且完全占用两条车道.请研究以下问题:1.根据视频1<附件1),描述视频中交通事故发生至撤离期间,事故所处横断面实际通行能力地变化过程.根据问题1所得结论,结合视频2<附件2),分析说明同一横断面交通事故所占车道不同对该横断面实际通行能力影响地差异.构建数学模型,分析视频1<附件1)中交通事故所影响地路段车辆排队长度与事故横断面实际通行能力、事故持续时间、路段上游车流量间地关系.假如视频1<附件1)中地交通事故所处横断面距离上游路口变为140M,路段下游方向需求不变,路段上游车流量为1500pcu/h,事故发生时车辆初始排队长度为零,且事故持续不撤离.请估算,从事故发生开始,经过多长时间,车辆排队长度将到达上游路口.附件1:视频1附件2:视频2附件3:视频1中交通事故位置示意图附件4:上游路口交通组织方案图附件5:上游路口信号配时方案图注:只考虑四轮及以上机动车、电瓶车地交通流量,且换算成标准车当量数.附件3视频1中交通事故位置示意图附件4附件5上游路口信号配时方案本题附件1、2地数据量较大,请竞赛开始后从竞赛合作网站“中国大学生在线”网站下载:试卷专题页面:试卷下载地址:2018高教社杯全国大学生数学建模竞赛题目<请先阅读“全国大学生数学建模竞赛论文格式规范”)B题碎纸片地拼接复原破碎文件地拼接在司法物证复原、历史文献修复以及军事情报获取等领域都有着重要地应用.传统上,拼接复原工作需由人工完成,准确率较高,但效率很低.特别是当碎片数量巨大,人工拼接很难在短时间内完成任务.随着计算机技术地发展,人们试图开发碎纸片地自动拼接技术,以提高拼接复原效率.请讨论以下问题:1. 对于给定地来自同一页印刷文字文件地碎纸机破碎纸片<仅纵切),建立碎纸片拼接复原模型和算法,并针对附件1、附件2给出地中、英文各一页文件地碎片数据进行拼接复原.如果复原过程需要人工干预,请写出干预方式及干预地时间节点.复原结果以图片形式及表格形式表达<见【结果表达格式说明】).2. 对于碎纸机既纵切又横切地情形,请设计碎纸片拼接复原模型和算法,并针对附件3、附件4给出地中、英文各一页文件地碎片数据进行拼接复原.如果复原过程需要人工干预,请写出干预方式及干预地时间节点.复原结果表达要求同上.3. 上述所给碎片数据均为单面打印文件,从现实情形出发,还可能有双面打印文件地碎纸片拼接复原问题需要解决.附件5给出地是一页英文印刷文字双面打印文件地碎片数据.请尝试设计相应地碎纸片拼接复原模型与算法,并就附件5地碎片数据给出拼接复原结果,结果表达要求同上.【数据文件说明】(1)每一附件为同一页纸地碎片数据.(2)附件1、附件2为纵切碎片数据,每页纸被切为19条碎片.(3)附件3、附件4为纵横切碎片数据,每页纸被切为11×19个碎片.附件5为纵横切碎片数据,每页纸被切为11×19个碎片,每个碎片有正反两面.该附件中每一碎片对应两个文件,共有2×11×19个文件,例如,第一个碎片地两面分别对应文件000a、000b.【结果表达格式说明】复原图片放入附录中,表格表达格式如下:(1)附件1、附件2地结果:将碎片序号按复原后顺序填入1×19地表格;(2)附件3、附件4地结果:将碎片序号按复原后顺序填入11×19地表格;(3)附件5地结果:将碎片序号按复原后顺序填入两个11×19地表格;(4)不能确定复原位置地碎片,可不填入上述表格,单独列表.。
2013国赛思路解析A题(仅供参考,内部使用,请勿外传)此题为交通运输类问题,可以视作优化类问题,而且本题重点在于目标的选取和目标函数的建立,而最优值的求解反而不是问题的重点(因为哪里会发生交通事故、持续时间、车流量等等都是不可控制的参数,本题几乎没有可决策变量)。
可以用到的知识有排队论,元胞自动机,模拟仿真等等,用这些手段来建立函数关系;关键概念:通行能力,指单位时间内通过断面的最大车辆数TC(traffic capacity)=n/t=vd (n为通过车辆数,t是时间,v为车辆平均速度,d是道路宽度);问题一:求出函数表达式TC=f(t),可以根据视频中的信息,隔一段时间求一次对应的TC值,再通过插值方法求出解f,或者深入研究事故发生时对车辆行进情况的变化机理来求解f,最后用图像或者解析式来表达出结果;问题二:求出泛函数表达式TC=g(LN),LN表示车道编号或其组合,此处TC代表问题一中的f函数,这个处理和问题一是一样的,可以用的方法也可以是直接从视频中读取,可以得到LN=(1,2)或(2,3)时的TC关于t的函数,如果采用机理分析方法,如排队论,元胞自动机来仿真这个过程,则可以求出LN=1,2,3时的情况;比较有两种形式:直观比较:将几个函数图像画在一起相互比较,就可以比较LN不同时,对通行能力的影响;数量化比较:可以将LN不同时的TC关于t的函数作差后积分,求得不同堵车形式对总的通行车辆数的影响;XXXXX第三题。
不让说的。
问题四:用问题三求出的函数表达式计算结果即可。
相关有用的资料我会及时上传群共享,大家加油!A题第一个问这个还是比较简单的,先数出事故未发生的时候,单位时间内通过的车辆,在数出发生事故之后单位时间内通过的车辆,注意,一定要在红灯变成绿灯之后的时候,也就车流量处于饱和的时候第二个问和第一个问差不多,先处理好有事故和没有事故界面的通行的实际能力,再将两个答案相比较即可得出答案第三个问这个应该是本模型最难的一题,首先将该路段的车辆排队长度设为因变量,而其他的都为自变量,此题就为一个权益的问题,通过控制变量的方法,不难理解,排队长度和事故持续时间成正比,和通行能力成反比(在不考虑路段的情况下,假如两根车道堵了一根,则能力变为原来的一半,而在单位时间内长度自然成为2倍,注:理想情况下),如果赋值出权益,则此题可解,一定要得出一个方程,这样才叫建模第四个问这个就相当于是第三个问的方程的实际应用,把所有的已知条件全部带入公式,得出答案。
道路上不断增加的交通流经常导致拥挤。
拥挤产生延误、降低流率、带来燃油损耗和负面的环境影响。
为了提高道路系统的效率,国内外许多研究者一直致力于车流运行模型的研究。
Daganzo[1]提出了一种和流体力学LWR 模型相一致的元胞传输模型,这种模型能用来模拟和预测交通流的时空演化,包括暂时的现象,如排队的形成、传播、和消散。
Heydecker 和Addison[2]通过研究车速和密度的因果关系分析和模拟了在变化的车速限制下的交通流。
Jennifer 和Sallissou[3]提出了一种混合宏观模型有效地描述了路网的交通流。
然而,拥挤也会由交通异常事件引起。
交通异常事件定义为影响道路通行能力的意外事件[4],如交通事故、车辆抛锚、落物、短期施工等,从广义角度看,还应包括恶劣天气与特殊勤务等。
异常事件往往造成局部车道阻塞或关闭,形成交通瓶颈,引起偶发性拥挤,这已经逐渐成为高速道路交通拥挤的主要原因[5],越来越多地受到研究者们的重视。
例如M. Baykal-Gursoy[6]等人提出了成批服务受干扰下的稳态M/M/c 排队系统模拟了发生异常事件的道路路段的交通流。
Chung[7]依据韩国高速公路系统监测的准确记录的大型交通事故数据库提出了一种事故持续时间预测模型。
当然,这些研究最终都是为了帮助缓解异常事件引起的交通拥挤。
交通异常事件发生后,事发地段通行能力减小,当交通需求大于事发段剩余通行能力时,车辆排队,产生延误,行程时间增加[8],交通流量发生变化。
本文以高速公路基本路段发生交通事故为例,主要分析了交通事故发生后不同时间段内事故点及其上游下游路段交通流量的变化,用于以后进一步的突发事件下交通流预测工作。
1 交通事故影响时间分析由于从交通事故发生到检测到事故、接警、事故现场勘测、处理、清理事故现场恢复交通,以及恢复交通后车辆排队不再增加都需要一定的时间。
这部分时间主要由三部分构成: 第一部分是事故发生到警察到达现场的时间T1; 第二部分是交通事故现场处理时间T2,由现场勘测、处理到事故族除、恢复交通; 第三部分是交通事故持续影响时间T3,这部分时间从恢复事故现场交通开始,到事故上游车辆排队不再增加,即排队开始减弱[9]。
在T1内,事故现场保持原状,没有进行处理,这里分两种情况考虑: ( 1) 当交通事故占用部分车道时,这时事故点的剩余通行能力Q s≠0,交通事故越严重,则相应Q s越小。
若事故点上游的交通需求Q <Q s,则车辆以较低的速度通过事故点,上游不会形成车辆拥挤排队; 若Q >Q s,则交通流可按事故点的剩余断面通行能力通过事故点,超过该通行能力的车流在事故点上游排队。
( 2) 当交通事故十分严重时,事故点的剩余通行能力Q s = 0,造成事发路段断流,事故点上游车辆排队,发生交通拥挤堵塞,进而排队一直向上游延伸。
在T2内,确认交通事故发生后,相关部门到现场处理异常事件,在此过程中,事故点交通可能会受到进一步影响,事故断面通行能力也随之发生变化[5],一般会变小,甚至变为0( 全封闭处理) ,视事件处理具体情况而定,事发点上游交通处于严重拥挤状态,车辆排队增加。
由于在交通事故接警时间T1和处理时间T2阶段事故点上游交通车辆产生排队,若没有车辆排队,则T3 = 0; 若有车辆排队,则当事故处理完毕、道路恢复交通时,排队车辆开始消散。
交通事故持续影响时间T3是事故处理完毕、道路恢复交通至车辆排队不再增加这段时间,即交通流消散波从车辆排队队列的头部传到尾部这段时间[9]。
2 事故路段车辆排队长度分析如图1 所示,设某高速公路基本路段长度为L( m) ,单方向车道数为n,单方向车道宽度为D( m) ,在道路上t = 0 时刻发生了一起交通事故,事故车辆占用道路宽度为b( m) ,长度为a( m) ,事故点上游路段长度为L'。
假设车辆的到达率为Q,在同级服务水平上事故发生断面通行能力为Q s,道路在正常条件下的单方向的通行能力为Q i。
图1 发生交通事故的高速公路基本路段本文暂只考虑如图1 所示的基本路段内的车辆排队长度,这里不同于以往文献的“排队长度”,以往文献中的“排队长度”没有区分不同的“阻塞行车道宽度”。
这里的“阻塞行车道宽度”不只是事故车辆实际占用宽度,还包括虚拟占用宽度,比如事故发生位置横跨在两车道之间,导致事故点只能通行一个车道宽度的车流,那么此时“阻塞行车道宽度”为两个车道的宽度。
设Q >Q s,m( T i) 为T i时间内事故点阻塞行车道宽度( 本文把单个车道宽度和车辆宽度看作同宽) ,L m( t) 为t 时刻事故点上游路段L'内车流以阻塞行车道宽度m 的排队长度,且w( T i) = u f 1 -k i1 + k i ( ) 2k ( ) j为T i时间内新产生的交通波[10]的速度,其中u f为该事故路段的自由流速度,即该路段的设计车速,可以通过城市地理信息平台GIS 得到道路基本数据; k i1、k i2分别为T i内事故点上游、事故点瓶颈段的交通密度,可以由交通检测系统监测得到; k j为该路段的交通堵塞密度,由道路的基本数据可以计算得到[9]。
2. 1 0 <t≤T1( 1) 若L'w( T1) <T1,则t <L'w( T1) 时,L m( T1) ( t) = w( T1) t <L';t≥L'w( T1) 时,L m( T1) ( t) = L'。
L m( T2) -m( T1) ( t) = 0。
( 2) 若T1≤L'w( T1),则L m( T1) ( t) = w( T1) t,L m( T2) -m( T1) ( t) = 0,2. 2 T1 <t≤T1 + T2在T2时间内,事发点断面通行能力一般会变化,设变为Q'S,则m( T i) 也会相应发生变化。
这里,还需要考虑一个时间,就是交通波w ( T2) 赶上w( T1) 的时间( 设为T'2) ,赶上之后车流以w( T2) 的速度、m( T i) 的宽度继续排队。
本文由于只考虑L'内的排队长度,所以考虑在T2时间内且在L'段内交通波w( T2) 是否赶上w( T1) ,即T'2若同时满足以下两个条件才需被考虑:L m( T1) ( T1) <L',L m( T1) ( T1) + w( T1) T'2≤L';T'2≤T2。
28 期陈诚,等: 交通事故影响下事发路段交通流量变化分析6905 2. 2. 1 L m( T1) ( T1) <L'若L m( T1) ( T1) + w( T1) T'2 ≤L'且T'2 ≤T2,则:当T2 <( L' -L m( T1) ( T1) -w( T1) T'2) /w( T2) + T'2时,L m( T1) ( t) = L m( T1) ( T1) + | w( T1) | ( t -T1) ,当t -T1≤T'2时L m( T1) ( t) = L m( T1) ( T1) + | w( T1) | T'2 + | w( T2) |( t -T1 -T'2) ,当t -T1 >T'2⎧⎨⎩⎪⎪⎪⎪时。
L m( T2) -m( T1) ( t) = w( T2) ( t -T1) ,当t -T1≤L' / w( T2) 时L m( T2) -m( T1) ( t) = L',当t -T1 >L' / w( T2 ) { 时。
当T2≥( L' -L m( T1) ( T1) -w( T1) T'2) / w( T2) +T'2时,L m( T1) ( t) = L m( T1) ( T1) + | w( T1) | ( t -T1) ,当t -T1≤T'2时L m( T1) ( t) = L m( T1) ( T1) + | w( T1) | T'2 + | w( T2) |( t -T1 -T'2) ,当( L' -L m( T1) ( T1) +| w( T1) | T'2) /| w( T2) | + T'2 ≥t -T1 >t'2时L m( T1) ( t) = L',当( L' -L m( T1) ( T1) + | w( T1) | T'2) /| w( T2) | + T'2 <t -T1≤t2⎧⎨⎩⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪时。
L m( T2) -m( T1) ( t) = w( T2) ( t -T1) ,当t -T1≤L' / w( T2) 时L m( T2) -m( T1) ( t) = L',当t -T1 >L' / w( T2 ) 时{ 。
若L m( T1) ( T1) + w( T1) T'2 >L',则L m( T1) ( t) =L m( T1) ( T1) + w( T1) ( t -T1) ,当t -T1≤L' -L m( T1) ( T1 ( ) ) / w( T1) 时L m( T1) ( t) =L',当t -T1 >L' -L m( T1) ( T1 ( ) ) / w( T1) 时{ 。
L m( T2) -m( T1) ( t) = w( T2) ( t -T1) ,当t -T1≤L' / w( T2) 时L m( T2) -m( T1) ( t) = L',当t -T1 >L' / w( T2) 时{ 。
若T'2 >T2,则L m( T1) ( t) = L m( T1) ( T1) + w( T1) ( t -T1) ,当t -T1≤L' -L m( T1) ( T1 ( ) ) / w( T1) 时L m( T1) ( t) =L',当t -T1 >L' -L m( T1) ( T1 ( ) ) / w( T1) 时{ 。
L m( T2) -m( T1) ( t) = w( T2) ( t -T1) ,当t -T1≤L' / w( T2) 时L m( T2) -m( T1) ( t) = L',当t -T1 >L' / w( T2) 时{ 。
2. 2. 2 L m( T1) ( T1) = L'L m( T1) ( t) = L',L m( T2) -m( T1) ( t) = w( T2) ( t -T1) ,当t -T1≤L' / w( T2) 时L m( T2) -m( T1) ( t) = L',当t -T1 >L' / w( T2 ) 时 。