数学建模论文车道被占用对城市道路通行能力的影响.docx
《数学建模论文车道被占用对城市道路通行能力的影响.docx》由会员分享,可在线阅读,更多相关《数学建模论文车道被占用对城市道路通行能力的影响.docx(36页珍藏版)》请在冰豆网上搜索。
数学建模论文车道被占用对城市道路通行能力的影响
高教社杯全国大学生数学建模竞赛
承诺书
我们仔细阅读了《全国大学生数学建模竞赛章程》和《全国大学生数学建模竞赛参赛规则》(以下简称为“竞赛章程和参赛规则”,可从全国大学生数学建模竞赛网站下载)。
我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题。
我们知道,抄袭别人的成果是违反竞赛章程和参赛规则的,如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺,严格遵守竞赛章程和参赛规则,以保证竞赛的公正、公平性。
如有违反竞赛章程和参赛规则的行为,我们将受到严肃处理。
我们授权全国大学生数学建模竞赛组委会,可将我们的论文以任何形式进行公开展示(包括进行网上公示,在书籍、期刊和其他媒体进行正式或非正式发表等)。
我们参赛选择的题号是(从A/B/C/D中选择一项填写):
A
我们的参赛报名号为(如果赛区设置报名号的话):
所属学校(请填写完整的全名):
HPU
参赛队员(打印并签名):
1.
2.
3.
指导教师或指导教师组负责人(打印并签名):
数模指导组
(论文纸质版与电子版中的以上信息必须一致,只是电子版中无需签名。
以上内容请仔细核对,提交后将不再允许做任何修改。
如填写错误,论文可能被取消评奖资格。
)
日期:
2013年9月16日
赛区评阅编号(由赛区组委会评阅前进行编号):
高教社杯全国大学生数学建模竞赛
编号专用页
赛区评阅编号(由赛区组委会评阅前进行编号):
赛区评阅记录(可供赛区评阅时使用):
评
阅
人
评
分
备
注
全国统一编号(由赛区组委会送交全国前编号):
全国评阅编号(由全国组委会评阅前进行编号):
车道被占用对城市道路通行能力的影响
摘要
针对问题一,为了描述视频1中交通事故发生至撤离期间事故所处横断面实际通行能力的变化过程,我们从两个方面进行研究:
从定性方面来看,首先通过视频1数出每30s的时间间隔内通过的各种车的数量,其次根据车辆折算系数转换成标准车当量数,最后通过matlab编程进行插值分析,得出道路通行能力随时间变化的曲线图。
由图可判断出事故处道路横断面的通行能力随时间的变化呈不规则的锯齿状;从定量方面来看,我们定义一个评价事故处横断面实际通行能力值
其中
为事故处横断面实际通行能力、
为基本通行能力、
为车道宽度及侧向净空修正系数、
为交通信号灯的影响系数、
为驾驶员适应性修正系数,最终得出
其中
为比例系数,进而可断定
与车辆通过事故横断面的速度
呈正相关。
针对问题二,我们采用图像对比法和单因素方差分析法进行分析,首先运用matlab软件在同一图中画出两视频中通行能力随时间的变化的对比图,由对比图像直观地看出了两者确实存在明显差异,再采用spss软件进行单因素方差分析,得出两视频中的事故横断面的通行能力
,说明同一横断面交通事故所占车道不同对该横断面实际通行能力影响确实存在明显的差异。
针对问题三,我们以30秒为时间间隔,通过统计最大队长作为对应时刻的路段车辆排队长度,以相同间隔统计出对应时段的事故路段上游车流量,并利用问题一统计出的事故横断面实际通行能力,运用spss软件进行多元线性回归分析,得出交通事故所影响的路段车辆排队长度
的回归方程为
其中
为横断面实际通行能力,
为事故持续时间,
为事故路段上游车流量。
为了更加精确的找出所选时刻的队长,我们利用了流体力学理论对本模型进行优化并作出了流量-密度图及波形时距图,通过分析得出了排队长度指标的一些变化规律。
针对问题四,我们利用模型三的结论,将车辆排队长度
的回归方程应用于本问题,其中
,
,带入回归方程可得车辆排队长度到达事故路段上游路口所需时间为
。
关键词:
插值法单因素方差分析流体力学理论多元回归分析
一、问题重述
车道被占用是指因交通事故、路边停车、占道施工等因素,导致车道或道路横断面通行能力在单位时间内降低的现象。
由于城市道路具有交通流密度大、连续性强等特点,一条车道被占用,也可能降低路段所有车道的通行能力,即使时间短,也可能引起车辆排队,出现交通阻塞。
如处理不当,甚至出现区域性拥堵。
车道被占用的情况种类繁多、复杂,正确估算车道被占用对城市道路通行能力的影响程度,将为交通管理部门正确引导车辆行驶、审批占道施工、设计道路渠化方案、设置路边停车位和设置非港湾式公交车站等提供理论依据。
视频1(附件1)和视频2(附件2)中的两个交通事故处于同一路段的同一横断面,且完全占用两条车道。
请研究以下问题:
1.根据视频1(附件1),描述视频中交通事故发生至撤离期间,事故所处横断面实际通行能力的变化过程。
2.根据问题1所得结论,结合视频2(附件2),分析说明同一横断面交通事故所占车道不同对该横断面实际通行能力影响的差异。
3.构建数学模型,分析视频1(附件1)中交通事故所影响的路段车辆排队长度与事故横断面实际通行能力、事故持续时间、路段上游车流量间的关系。
4.假如视频1(附件1)中的交通事故所处横断面距离上游路口变为140米,路段下游方向需求不变,路段上游车流量为1500pcu/h,事故发生时车辆初始排队长度为零,且事故持续不撤离。
请估算,从事故发生开始,经过多长时间,车辆排队长度将到达上游路口。
附件1:
视频1
附件2:
视频2
附件3:
视频1中交通事故位置示意图
附件4:
上游路口交通组织方案图
附件5:
上游路口信号配时方案图
注:
只考虑四轮及以上机动车、电瓶车的交通流量,且换算成标准车当量数
二、问题分析
2.1问题一的分析
为了描述出视频1中交通事故从发生至撤离期间事故所处横断面实际通行能力的变化过程,我们分别采用了定性和定量的方法进行分析。
首先我们根据视频1搜集原始数据,以30s为时间间隔,查找出了每个时间间隔内的小车数、大车数,接着根据车辆折算系数(见附录二)将所有车辆化为标准车当数量。
为了能够直观的看出实际通行能力的变化过程,我们以时间点为横坐标、道路的通行能力为纵坐标通过matlab编程进行插值运算,作出道路通行能力随时间变化的曲线图。
为了进一步定量的分析该变化的过程,我们定义了评价事故所处横断面实际通行能力的值
其中
为基本通行能力、
为车道宽度及侧向净空修正系数、
为交通信号灯的影响系数、
为驾驶员适应性修正系数,通过查找数据分析得出了实际通行能力与通过事故处车辆的速度呈正相关的结论,表达式如下
其中
为比例系数通过事故横断面的速度、
为车辆通过事故横断面的速度。
2.2问题二的分析
为了说明同一横断面交通事故所占车道不同对该横断面实际通行能力影响的差异,我们采用了图像分析法和单因素方差分析法。
对于图像分析,我们在同一个图像中做出了两视频中事故发生时横断面的实际通行能力的变化趋势,得到了视频1中的通行能力明显低于视频2,即视频1中的事故发生处所占车道对横断面实际通行能力影响较大。
对于单因素方差分析法,由表六可得方差分析结果
,显著性水平
此结果表明同一横断面交通事故所占车道不同对该横断面实际通行能力的影响存在显著性差异。
综合两种方法可以说明同一横断面交通事故所占车道不同对该横断面实际通行能力影响的差异还是很明显的,考虑到实际问题,由于附件3给出的不同车道的车辆分布有差异,故所占车道不同对该横断面实际通行能力的影响存在差异是合理的。
2.3问题三的分析
为了确定出事故发生路段车辆的排队长度与事故横断面实际通行能力、事故持续时间、路段上游车流量间的关系,我们首先统计出每隔相同时间间隔内事故发生路段车辆的排队长度、事故横断面实际通行能力、事故持续时间及路段上游车流量,然后利用多元回归分析确定出事故发生路段车辆的排队长度与事故横断面实际通行能力、事故持续时间及路段上游车流量的函数表达式。
2.4问题四的分析
为了求得车辆排队长度到达上游路口所需要的时间,我们将
和
,带入问题三所建立的回归方程,进而可求得车辆排队长度到达事故路段上游路口所需时间为
。
三、符号说明
——道路横断面的实际通行能力;
——理想条件下道路横断面的基本通行能力;
——车道宽度及侧向净空修正系数;
——交通信号灯的影响系数;
——驾驶员适应性修正系数;
——视频1中事故处横断面的实际通行能力;
——视频2中事故处横断面的实际通行能力;
——车辆排队的长度;
四、模型假设
1、在一定的时间内,不同种车型的速度及性能是一致的。
2、只考虑四轮及以上机动车、电瓶车的交通流量,且换算成标准车当量数。
3、视频1、2中的两个交通事故处于同一路段的同一横断面,且完全占用两条车道。
4、在同一段时间内,两视频中所示路段发生的交通事故均为一次。
5、车流中各单个车辆的行驶状态与它前面的车辆完全一致, 而且即使车道下游没有车辆, 某个特定车流状态的领头车辆也不提高它的速度。
五、模型的建立和求解
5.1问题一的求解
道路通行能力是指道路上某一点某一车道或某一断面处,单位时间内可能通过的最大交通实体(车辆)数,亦称道路通行能量,用辆/h或用辆/昼夜或辆/秒表示,车辆多指小汽车,当有其它车辆时,均采用等效通行能力的当量标准车辆为单位(pcu)。
当道路上实际交通量小于其通行能力时,道路上行驶车辆处于自由行驶状态,车速较高,交通密度较小,车头时距分布规律符合负指数分布,车辆能实行超车;当道路上实际交通量接近或等于其通行能力时,道路上行驶的车辆用接近匀速的车速跟踪行驶,出现车队行驶现象,车头时距分布接近均数值;当道路上实际交通量超过其通行能力时,道路上行驶车辆密度增大,车速降低,出现交通拥挤和阻塞现象。
下面我们对视频1当中提到的道路堵塞问题进行如下研究。
5.1.1原始数据的采集
为了得到事故处横断面实际通行能力的变化过程,根据视频1(附件1),我们选取事故发生到事故撤离过程为研究对象,以30秒为时间间隔,数出每隔30秒通过事故横断面处的各种车辆数,并利用不同车辆的转换系数(见附录一)将所有车辆转换成标准当量车数,统计结果如表一
序号
起始时间
终止时间
时间点(min)
小车数(辆)
大车数(辆)
标准车当量数(pcu)
通行能力(pcu/h)
1
16:
42:
32
16:
43:
02
0.5
7
2
10
1200
2
16:
43:
02
16:
43:
32
1
6
2
9
1080
3
16:
43:
32
16:
44:
02
1.5
11
0
11
1320
4
16:
44:
02
16:
44:
32
2
7
1
8.5
1020
5
16:
44:
32
16:
45:
02
2.5
8
0
8
960
6
16:
45:
02
16:
45:
32
3
7
0
7
840
7
16:
45:
32
16:
46:
02
3.5
8
0
8
960
8
16:
46:
02
16:
46:
32
4
8
1
9.5
1140
9
16:
46:
32
16:
47:
02
4.5
9
0
9
1080
10
16:
47:
02
16:
47:
32
5
6
0
6
720
11
16:
47:
32
16:
48:
02
5.5
8
1
9.5
1140
12
16:
48:
02
16:
48:
32
6
10
0
10
1200
13
16:
48:
32
16:
49:
02
6.5
10
0
10
1200
14
16:
49:
02
16:
49:
32
7
9
0
10
1200
15
16:
49:
32
16:
50:
02
7.5
7
1
8.5
1020
16:
50:
02
16:
50:
26
16
16:
50:
26
16:
50:
56
8
10
0
10
1200
17
16:
50:
56
16:
51:
26
8.5
9
0
9
1080
18
16:
51:
26
16:
51:
56
9
9
0
9
1080
19
16:
51:
56
16:
52:
26
9.5
8
1
9.5
1140
20
16:
52:
26
16:
52:
56
10
9
0
9
1080
21
16:
52:
56
16:
53:
26
10.5
6
2
9
1080
22
16:
53:
26
16:
53:
52
11
8
1
9.5
1140
16:
54:
03
16:
54:
07
23
16:
54:
07
16:
54:
37
11.5
10
1
11.5
1380
24
16:
54:
37
16:
55:
07
12
6
1
7.5
900
25
16:
55:
07
16:
55:
37
12.5
10
0
10
1200
26
16:
55:
37
16:
56:
05
13
7
1
8.5
1020
16:
57:
54
16:
57:
56
27
16:
57:
56
16:
58:
18
13.5
8
1
9.5
1140
16:
59:
07
16:
59:
15
28
16:
59:
15
16:
59:
31
14
11
1
12.5
1500
16:
59:
43
16:
59:
57
表一视频1中每隔30秒通过事故处横断面的各种车辆数
5.1.2通行能力变化的趋势
以时间点为横坐标、道路的通行能力为纵坐标通过matlab编程(程序见附录一)进行曲线拟合,作出道路通行能力随时间变化的曲线图,如图一
图一视频1事故发生处横断面通行能力的变化曲线
由图一可以看出,视频1中事故发生处横断面的通行能力整体上呈周期性变化,在一个周期内,相邻两个极小值之间的时间间隔约为60s,又由于上游的十字路口的信号灯的周期也为60s,这就说明了事故处横断面的实际通行能力确实与信号灯有一定的关系,并且当事故发生14分钟以后,此时通行能力明显呈上升趋势。
5.1.3道路实际通行能力评价指标体系的构建
为了定量地分析出道路的实际通行能力,我们首先依据设计速度确定基本通行能力,然后结合实际道路、交通条件对基本通行能力进行车道宽度及侧向净空修正、交通信号灯的影响系数和驾驶员适应性修正,其基本公式定义为
其中
为事故处横截面实际通行能力、
为基本通行能力、
为车道宽度及侧向净空修正系数、
为交通信号灯的影响系数、
为驾驶员适应性修正系数。
(1)基本通行能力的确定
基本通行能力是指在理想的道路、交通条件以及某一控制速度下,一条车道的一个断面在一小时内能通过的最大车辆数,是一个经过大量理想实验所确定的数值(查资料数据如表二),首先我们测出车辆通过事故横断面的速度
,接着通过表二中数据运用插值法来确定基本通行能力
,显然
与
呈正相关,不妨设
车辆通过事故横断面的速度
与基本通行能力
的对应关系如图表二
设计速度/(km/h)
基本通行能力/[pcu/(h*ln)]
120
2200
100
2100
80
2000
60
1800
表二理想条件下的基本通行能力
(2)对于车道的宽度及侧向的净空修正系数
的确定
查阅相关资料可知车道的宽度及侧向的净空修正系数
为
其中
为一条车道的宽度。
对于本题有
,因此
(3)交通信号灯影响系数
的确定
其中
为绿信比,即绿灯的有效时间和周期的比;
为事故发生处到上游交叉路口间的距离。
对于本问题
,因此
(4)驾驶员适应性修正系数的确定
驾驶员条件对通行能力的修正系数,一般在0.9—1之间,对于事故现场处,由于侧向净空受限,加之车辆系数的调整,此时结合各因素确定
。
综上可得事故处横断面实际通行能力的综合评价指标为
由于
(
),故事故横断面实际通行能力与车辆通过事故横断面的速度
呈正相关,由于交通信号灯的影响,导致从上游过来的车辆数呈现周期性的变化,到达事故处的速度也发生相应的变化,进而导致事故横断面的实际通行能力呈周期性变化。
5.2问题二模型的建立和求解
5.2.1原始数据的采集
为了说明同一横断面交通事故所占车道不同对该横断面实际通行能力影响的差异,如同5.1.1,根据视频2,我们选取事故发生到事故撤离过程为研究对象,以30秒为时间间隔,数出每隔30秒通过事故处横断面的各种车辆数,并利用不同车辆的转换系数将所有车辆转换成标准当量车数,如表三
序号
起始时间
终止时间
时间点(min)
小车数(辆)
大车数(辆)
标准车当量数(pcu)
通行能力(pcu/h)
1
17:
34:
16
17:
34:
46
0.5
6
0
6
720
2
17:
34:
46
17:
35:
16
1
13
2
16
1920
3
17:
35:
16
17:
35:
46
1.5
12
1
13.5
1620
4
17:
35:
46
17:
36:
16
2
8
1
9.5
1140
5
17:
36:
16
17:
36:
46
2.5
7
0
7
840
6
17:
36:
46
17:
37:
16
3
12
1
13.5
1620
7
17:
37:
16
17:
37:
46
3.5
10
1
11.5
1380
8
17:
37:
46
17:
38:
16
4
7
2
10
1200
9
17:
38:
16
17:
38:
46
4.5
12
0
12
1440
10
17:
38:
46
17:
39:
16
5
10
1
11.5
1380
11
17:
39:
16
17:
39:
46
5.5
4
4
10
1200
12
17:
39:
46
17:
40:
16
6
10
0
10
1200
13
17:
40:
16
17:
40:
46
6.5
9
1
10.5
1260
14
17:
40:
46
17:
41:
16
7
8
1
9.5
1140
15
17:
41:
16
17:
41:
46
7.5
11
0
11
1320
16
17:
41:
46
17:
42:
16
8
9
1
10.5
1260
17
17:
42:
16
17:
42:
46
8.5
14
0
14
1680
18
17:
42:
46
17:
43:
16
9
9
2
12
1440
19
17:
43:
16
17:
43:
46
9.5
11
0
11
1320
20
17:
43:
46
17:
44:
16
10
6
1
7.5
900
21
17:
44:
16
17:
44:
46
10.5
8
0
8
960
22
17:
44:
46
17:
45:
16
11
9
1
10.5
1260
23
17:
45:
16
17:
45:
46
11.5
10
0
10
1200
24
17:
45:
46
17:
46:
16
12
9
0
9
1080
25
17:
46:
16
17:
46:
46
12.5
8
2
11
1320
26
17:
46:
46
17:
47:
16
13
6
0
6
720
27
17:
47:
16
17:
47:
46
13.5
5
2
8
960
28
17:
47:
46
17:
48:
16
14
8
2
11
1320
29
17:
48:
16
17:
48:
46
14.5
10
0
10
1200
30
17:
48:
46
17:
49:
16
15
8
1
9.5
1140
31
17:
49:
16
17:
49:
46
15.5
11
0
11
1320
32
17:
49:
46
17:
50:
16
16
10
1
11.5
1380
33
17:
50:
16
17:
50:
46
16.5
11
0
11
1320
34
17:
50:
46
17:
51:
16
17
11
0
11
1320
35
17:
51:
16
17:
51:
46
17.5
7
1
8.5
1020
36
17:
51:
46
17:
52:
16
18
10
0
10
1200
37
17:
52:
16
17:
52:
46
18.5
9
1
10.5
1260
38
17:
52:
46
17:
53:
16
19
9
0
9
1080
39
17:
53:
16
17:
53:
46
19.5
6
1
7.5
900
40
17:
53:
46
17:
54:
16
20
12
0
12
1440
41
17:
54:
16
17:
54:
46
20.5
10
1
11.5
1380
42
17:
54:
46
17:
55:
16
21
11
0
11
1320
43
17:
55:
16
17:
55:
46
21.5
6
2
9
1080
44
17:
55:
46
17:
56:
16
22
10
1
11.5
1380
45
17:
56:
16
17:
56:
46
22.5
8
2
11
1320
46
17:
56:
46
17:
57:
16
23
8
1
9.5
1140
47
17:
57:
16
17:
57:
46
23.5
10
0
10
1200
48
17:
57:
46
17:
58:
16
24
7
2
10
1200
49
17:
58:
16
17:
58:
46
24.5
7
2
10
1200
50
17:
58:
46
17:
59:
16
25
4
1
5.5
660
51
17:
59:
16
17:
59:
46
25.5
6
2
9
1080
52
17:
59:
46
06:
00:
16
26
12
0
12
1440
53
18:
00:
16
18:
00:
46
26.5
9
1
10.5
1260
54
18:
00:
46
18:
01:
16
27
8
1
9.5
1140
55
18:
01:
16
18:
01:
46
27.5
9
0
9
1080
56
18:
01:
46
18:
02:
16
28
9
0
9
1080
57
18:
02:
16
18:
02:
46
28.5
9
1
10.5
1260
58
18:
02:
46
18:
03:
16
29
10
1
11.5
1380
表三视频2中每隔30秒通过事故处横断面的各种车辆数
5.2.2两次事故处道路横断面实际通行能力对比图