2021数学建模竞赛承诺书

2013高教社杯全国大学生数学建模竞赛

承诺书

我们仔细阅读了《全国大学生数学建模竞赛章程》和《全国大学生数学建模竞赛参

赛规则》（以下简称为“竞赛章程和参赛规则”,可从全国大学生数学建模竞赛网站下载）。

我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网

上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。

我们知道，抄袭别人的成果是违反竞赛章程和参赛规则的，如果引用别人的成果或

其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文

引用处和参考文献中明确列出。

我们郑重承诺，严格遵守竞赛章程和参赛规则，以保证竞赛的公正、公平性。如有

违反竞赛章程和参赛规则的行为，我们将受到严肃处理。

我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展

示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。

我们参赛选择的题号是（从A/B/C/D中选择一项填写）：

我们的参赛报名号为（如果赛区设置报名号的话）：

所属学校（请填写完整的全名）：___________________________________________

参赛队员（打印并签名）：1._______________________________________________

2._______________________________________________

3.___________________________________________

指导教师或指导教师组负责人（打印并签名）：

（论文纸质版与电子版中的以上信息必须一致，只是电子版中无需签名。以上内容

请仔细核对，提交后将不再允许做任何修改。如填写错误，论文可能被取消评奖资格。）

日期：2013年9月15日

赛区评阅编号（由赛区组委会评阅前进行编号）：

2013高教社杯全国大学生数学建模竞赛

编号专用页

赛区评阅编号（由赛区组委会评阅前进行编号）：

赛区评阅记录（可供赛区评阅时使用）：

评

阅

人

评

分

备

注

全国统一编号（由赛区组委会送交全国前编号）：

全国评阅编号（由全国组委会评阅前进行编号）：

车道被占用对城市道路通行能力的影响

摘要

本文是研究因交通事故导致车道或道路横断面通行能力下降的问题，对于如何正确估算车道被

占用对城市道路通行能力的影响程度，

在问题一，建立通行能力折减系数模型，公路上发生交通事故，事故点的通行能力

必然降低，实际同行能力与理论通行能力的比值即为折减系数，对视频一中的发生事故

后车流量进行统计，并标准当量化，发现其折减系数较为均匀的下降。从而我们知道事

故发生后，交通能力总体在近似线性地降低

在问题二，利用比例相关法，有附件3可知公路上不同道车流量不同，1和2道流量

和比例占0.45,2和3道流量和比例占0.55,且视频二中的折减系数大于视频一中折减系

数,可知两者呈正相关.所以占不同流量的道对通行能力的影响不同，占用流量比例大的

道则通行能力降低越严重

在问题三，建立交通流量波动模型，并运用AutoCAD对模型部分图像进行绘制，得出

最长排队长度公式和消散时间公式：x=%＞（7一"）=3一勺濡）x（T-Q,将

跖

排队长度与实际通行能力，事故持续时间，上游的需求量关系由方程体现出来。但是该

模型对于排队长的计算与实际偏差较为严重，随后本文对该模型进行改进，运用到Q-K

模型（二次抛物线），并考虑红绿灯信，得到队长公式：

XJS-S2)MQ「R田一R3s2+(R+g+&一T)]Q

并预测了排队最长为146m,

S)+e.)

'KS2{QX-S)+KX(S-2KS(S2-

与实际情况贴近。

问题四，引用流体模型，与问题三中结论，得到每分钟参与排队的车辆数，阻塞密

度得，在事故持续不撤离的条件下，排队到上游路口所用的时间在（10,11）分钟之间,

为交通管理部门正确引导车辆行驶、审批占道施工、设计道路渠化方案、设置路边

停车位和设置非港湾式公交车站等提供理论依据。

关键词：交通拥堵，折减系数，车流波动论，流体模型Q-K模型

一'问题重述

车道可能由于交通事故、路边停车、占道施工等因素而被占用，导致车道横断面通

行能力在单位时间内出现降低的现象。由于城市道路具有交通流密度大、连续性强等特

点，一条车道被占用，也可能降低路段所有车道的通行能力，即使时间短，也可能引起

车辆排队，出现交通阻塞。如处理不当，甚至出现区域性拥堵。

正确估算车道被占用对城市道路通行能力的影响程度，将为交通管理部门正确引导

车辆行驶、审批占道施工、设计道路渠化方案、设置路边停车位和设置非港湾式公交车

站等提供理论依据。

视频1（附件1）和视频2（附件2）中的两个交通事故处于同一路段的同一横断面,

且完全占用两条车道。研究以下问题：

1.根据视频1（附件1）,描述视频中交通事故发生至撤离期间，事故所处横断面实

际通行能力的变化过程。

2.根据问题1所得结论，结合视频2（附件2）,分析说明同一横断面交通事故所占

车道不同对该横断面实际通行能力影响的差异。

3.构建数学模型，分析视频1（附件1）中交通事故所影响的路段车辆排队长度与事

故横断面实际通行能力、事故持续时间、路段上游车流量间的关系。

4.假如视频1（附件1）中的交通事故所处横断面距离上游路口变为140米，路段下

游方向需求不变，路段上游车流量为1500pcu/h,事故发生时车辆初始排队长度为零，且

事故持续不撤离。估算从事故发生开始，经过多长时间，车辆排队长度将到达上游路口。

二'问题的分析

交通车流具有流体性，当由于交通事故会占用部分道路资源时，车流行驶至事故区

域，车流会重新分布，因为部分车道封闭，处于封闭车道上的车辆不得不变换车道，行

驶至开放车道，由此产生交通车流融合的合流。

（1）对于问题一，交通事故的发生改变了公路路段的道路条件，由于事故车道受阻,

从而形成道路瓶颈，造成主线车速降低，并以波的形式向车辆行进的逆方向传递，因此,

公路上发生交通事故，事故点的通行能力必然降低，通常都先计算基本交通能力，再根

据实际同行能力（数据从视频一收集得到）来确定折减系数。

问题一通过对折减系数的描述来反应事故发生至撤离期间，事故所处横断面实际通

行能力的变化过程。

（2）对于问题二，根据问题一所得结论，结合从视频一、视频二中收集的数据计算

得到相应的折减系数，说明同一横断面交通事故所占车道不同对该横断面实际通行能力

影响的差异并分析其原因。

（3）对于问题三，部分车道被出事车辆所占用，通行能力下降，上游交通需求量大

于路段现行交通能力时，就会形成排队，利用车流波动理论模型对排队长度进行定量的

分析研究并进行相应的改进。

（4）对于问题四，交通事故发生后，事故影响不同时间段内车辆排队长度的不同变

化，通过对问题三模型的解答我们将会解得，排队长为140米时所对应的时间。

根据视频可以得到许多关于车流量、通行能力、车速等方面的数据，但这些数据通

常具有不完全性和随机性，所以本文在信息提取和模型建立上提出了适当的简化处理方

法。

三、基本假设

1）视频中所给出的信息能正确地反应集结和消散的车流量；

2）事故地点上下游车流量稳定；

3）进入事故路段时，车辆根据下游方向的选择而选择相应的车道;

四、符号说明

符号说明单位

C事故条件下的通行能力辆/h

G正常情况下的通行能力辆/h

P折减系数\

5.事故解除前通行能力辆/h

事故解除前车流速度km/h

事故解除前车流密度辆/km

事故解除后通行能力辆/h

52

匕2事故解除后车流速度km/h

心事故解除后车流密度辆/km

峪阻塞密度辆/km

2,事故时路段的车流需求量辆/h

0事故后路段的车流需求量辆/h

/事故时路段的车速km/h

事故后路段的车速km/h

V2

事故处理时间min

To

z队列持续增加至最长时间min

T两波相遇时间min

五'模型的分析'建立与求解

5.1问题一的解答:

根据附件视频一收集的数据，可以得到单位时间的不同类型的车流量。本文将车辆

类型分为三类，即公交车、小汽车和电瓶车（用1,2,3代替）。

交通事故的发生改变了公路路段的道路条件，由于事故车道受阻，从而形成道路瓶

颈，造成主线车速降低，并以波的形式向车辆行进的逆方向传递，因此，公路上发生交

通事故，事故点的通行能力必然降低，通常都先计算基本交通能力，再根据因素来确定

折减系数，不妨设折减数为“，用〃值表示其通行能力的大小。其通行能力计算公式如

下：

c=品•P

式中：

C。为事故条件下的通行能力；

C为正常情况下的通行能力；

P为交通事故下通行能力的折减系数。

根据《城市道路设计规范》川（CJJ37-90）一条道路的理论通行能力值，见下表：

表1：城市道路与车流量规范表

V(km/h)6050403020

Co(pcu/17301690164015501380

h)

查阅资料后，得知在城市内的道路行驶时一般不得超过40km/h,可以得到:

Co=1640pcu/h

在视频中可以得到每分钟通过的各种车辆数，由于不同种类的车的尺寸不同，这样

就对通过能力的计算产生影响。本文先对各种车辆进行当量标准化，根据《公路工程技

术标准》⑵，将小汽车作为标准车型给出车辆换算系数，见表2：

表2：不同车辆的换算系数

车型公交车小汽车电瓶车

换算系数210.5

根据从视频中收集的数据和换算系数，可以得到自视频开始时每分钟标准当量后的

车流量（部分，全部数据间附录1）。（单位：辆/分）数据见表3：

表3：通过事故截面车流量变化表

时刻(min)123・・・15161718192021

1111.・・1211577

车辆类型

281316・・・13121916183014

.・・

」35694005252

标准当量12.51822.5・・・17162120.52916.529

将车流量换算为标准单位制后与正常情况下的通行能力相比，得到各个时刻的折减

系数。根据题目要求，本文取事故后一分钟至事故结束前一分钟的数据并作图。

表4：各时刻的折减系数

事故后

时间1234567891011

通行能

力1470114012301050129012001260111012301020990

折减系

数0.900.700.750.640.790.730.770.680.750.620.60

其中，折减系数的平均值为0.72,方差是0.0065

折减系数

图1：视频一中事故发生后各时刻的折减系数

折减系数表示事故后的实际通行量与正常情况下通行量的比，其值越小表示拥堵情

况越严重。根据图1可知，在事故发生后折减系数呈波动下降趋势，但其方差是0.0065,

可以认为折减系数比较均匀的下降，即拥堵情况越来越严重。

5.2问题二的解答：

对视频二处理和视频一基本一致，从视频中收集每分钟的车流量，通过计算得到其

车流量和折减系数（部分，全部见附录2）,见表5：

表5：附件2中事故发生后各时刻的折减系数

时刻1234•••232425262728

通行

1710135014701560・・・189012301140120014101260

能力

折减

1.040.820.900.95・・・1.150.750.700.730.860.77

系数

其中，折减系数的平均值为0.84,方差为0.03

根据表中的关于折减系数的数值可以得到折减系数在事故发生后各个时刻的值作

图。

折减系数

事故后时间

图2：事故后各时刻的折减系数

如图所示，折减系数波动较大，方差为0.03,是附件一中数据的6倍，说明其拥堵

程度并不均匀。平均值为0.84,说明其拥堵情况总体来说并不是很严重。

通过表4和表5、图1和图2的对比，可知视频一中拥挤情况比二中拥挤情况严重

且一中的拥堵程度随时间变化比较均匀。

原因分析：当车流进入合流后，封闭车道上的车辆需要进行换道行驶，这就需要与

开放车道上的车辆争夺道路使用权，但是二者的机会并不均等，开方车道上的车辆将享

有优先使用权，而封闭车道上的车辆等到开发车道的车辆走后才有权使用车道。示意图

如下（方向：从左到右；从下至上一次为1,2,3车道）：

图3：附件二中车辆通过事故截面示意图

由附件3可知，车辆在下游路径选择上有不同。

表6：车辆在下游换道方式的比例

下游换道方式左转直行右转

所占比例35%44%21%

直行的车辆在事故截面无论向左还是向右对拥堵的贡献量假设是相同的。由表可以

看出，选择左转人数的比例(35%)比选择右转(21%)的人多了近67虬根据交通规则，

选择左转的在3车道，选择右转的在1车道。当交通事故占领2,3车道时，选择右转的

车辆较多，导致排队情况比较严重；当交通事故占领1,2车道时，选择左转的车辆较

少，导致排队情况相对轻松。

5.3问题三的解答：

5.3.1模型的理论基础

从交通事故发生到事故消除，这期间由于部分车道被出事车辆所占用，因此

该路段的通行能力下降，上游交通需求量大于路段现行通行能力时，就会形成排

队。当事故解除以后，路段通行能力有所回升，此时排队仍然存在，所以根据流

量密度关系，时的通行能力还达不到该路段原有通行能力。当排队彻底消散以后，

通行能力恢复到原有水平，该路段恢复正常行车。

假设上游交通需求量大于事发路段现有通行能力，到达车流在事故地点陆续

减慢速度甚至停车而集结成密度较高的队列，事故解除后，由于路段通行能力的

恢复，排队车辆又陆续加速而疏散成一列具有适当密度的车队，流中两种不同密

度部分的分界面经过一辆辆车向车队后部传播的现象，称为车流波动⑶。

上游车流由高速低密的畅通状态转变为低速高密的拥挤状态，从而形成集结

波，波面以一定的速度向车队的后方传播；事故解除后，除了集结波继续向车队

后方传播外，在车队的前方又形成了消散波，波面同样向车队后方传播。当消散

波的速度大于集结波的速度时排队消散终能完成。

5.3.2模型建立

由车流波动理论可知，波速

公式为：

Wp=(Qx-Qy)/(K「K,)(1)

Q=vxK(2)

Um

B

Or/min

A

图4：车队运行变化图

1993年，格林希尔茨（Green-shields）提出了速度一密度线性关系模型:

u=O（l-£）（3）

Kj

由以上⑴、（2）、⑶式可以推导

出波速与密度的关系

忆，一（1一一Q一）⑷

车

辆

数

、

辆

o

T*时间/分

-到达车辆数

通过车辆数

图5：车辆累计消散过程图

图5为事故发生后累计车辆一时间图，实线表示交通需求流量，点划线表示

通过能力。叙述简便，对所用符号说明如下：事故发生时堵塞了部分车道，该

路段通行能力下降M；相应密度上升Si；交通事故处理所需时间7；；事故解除

后到车队消散前通行能力回升为§2；车流密度相应地下降为舄2。其中路段的通

行能力由图4中点划线的斜率表示。路段上游交通需求流量为。、由图5

中实线斜率表示；持续时间为刀、心，相应的车流密度K|，K

在图5中可以看出当两条折线疝交时表示车队消散，所需2O时间T。但无法计

算出排队长，可用车流波动理论进行求解。

t/min

图中每条曲线表示每一辆车运行的时间-空间轨迹。横轴表示时间，纵轴表示与事故点

的相对位置，原点0表示事故发生点，纵轴的负半轴表示事故点的上游，正半轴表示事

故点的下游，虚线OA,0B表示集结波，CB表示消散波，其斜率的绝对值表示波速，斜

率的正负号表示波传播的方向。两波相遇的时间为T,当集结波与消散波在T>0的范围

内有交点时，表示车队可以在有限时间内消散，否则不能消散。

首先假设两波相遇之前该路段需求流量始终为9,0A与CB相交处表示排队向上游延伸

到的最远处，设两波相遇时的时间为T,集结波波速为叱.2,消散波波速为购.3,则根

据两波相遇时波传播的距离相等这一关系可知：

W12xT=W2ix（T-T0）（5）

其中：

Sj—S2KS]+KS2、

"2,3=*=勺（1

KKK

则：T='~^-~^XT0

K「KS2

若T>T],则说明在车队消散之前该路段上游需求量发生了变化，需求流量变为Q2,相

应的密度变化为K2。所以（5）式改写为：

叱.2X（+也.2、（7-7；）=%.3*（7一册）

其中：=0（1—

则，丁_（马-"Ks2）xT0+（七-&）xE

一一

K2-KS2

根据公式：x=叱,3X（T—4）=V/（1_KSL：K.）x（T_4）

Kj

从而解出此次事故的排队长，即集结波的最大范围。

由图3可知，车队消散时间为

r^T+—

其中：匕“一路段通行能力，S2时的行车速度,

5.3.3模型检验

5.3.3.1数据准备

根据表3,可得在事故发生时的车流量为20辆/分，BPS,=1200pcu/h;

根据视频此时车速匕i=12km/h；

此时车流密度K“=S1/.=100辆/km；

事故持续时间游视频可以测得，7；=18.8min；

事故解除后通行能力邑=1500pcu/h；

根据视频测得此时车速匕2=20km/h；

此时车流密度&2=§2/匕2=75辆/km；

阻塞密度，通过网络查阅城市道路的KJ=110辆/km；

到达车辆数可以看做通过车辆数（表3）与等待车辆数之和。等待车辆数数据亦

可从视频中收集，最终得到车辆到达数。见图7

0=2485pcu/h；。2=1740pcu/h；

匕=36km/h；匕=50.4km/h；

可得K1=2/W=69辆/km；K2=Q2/V2=35辆/km

5.3.3.2数据计算

假设两波相遇时车流需求量保持不变0=2485pcu/h,则

K,—K-K”13—,、个

T=I\

%KS2O

由此可知，假设不成立。两波相遇时的车流需求量为02=1740pcu/h；则

(K-K,.-K„)xT+(K-K.)xT.

T=」i————以n~~~^=38.1min

a-KS2

Ks+

x=W^x(T-T0)=vf(1-'^-)x(T-7^)=7.55km

KJ

T*=T+土=38.4min

%

5.3.3.3结果分析

由数据计算结果可以看出，次模型对两波相遇的时间（排队时间）预测比较接近实

际，对排队长度的空间预测与实际有偏差。

可能的原因是：车辆的集结与消散并不是均匀的；

本文的数据在收集上具有一定的随机性；

路段不是无限长，模型中没有考虑红绿灯对车流量的影响。

由以上几个原因需要对模型进行改进。

5.3.4模型改进⑷

图（7）中上面一条折线的斜率表示需求流量。和其中。持续时间为工，对应的密度

记为K2,随后流量下降到口，对应密度记为K2,因交通事故堵塞了部分车道，使通行能

力下降为S,，密度相应地上升为Ksl,持续时间为r,在随后的时间均为排除故障的完

全封路期,通行能力变为零,密度达到最大值K）,随着故障被部分排除,通行能力恢复到

S2,对应密度为Ks?,持续时间为冬，故障完全排除后,通行能力达到最大值S,对应密度

记为Ks,S的持续时间记为T,,如果7；2a+R?+R?+Ts

从图中不难推出：

丁_（N+氏2+氏3）。1一R|S]—R3s2（]）

lS-Q,

如果7J</?|+/??+&+（

二T◎-峪—R3s2+（片+%+%—工）。2，八

'S-Q2

图（8）中下边的折线0BCHM表示流量的供给，即通行能力。值的指出的是，上下折线之

间的垂直距离并不是排队车数，由上图是不可能求出排队占用道路的长度。

2用集散波理论表示排队的真实长度

图8之下图是与上图相对应的，时间轴的刻度完全一样，下图能把车流的阻塞一一消

散过程详细的表示出来，它包含了车流的所有信息一一流量、速度和密度的变化，

排队车辆数及其所占路长总延误及消散时间等等。下图的纵轴表示道路中心线上的不

同位置X。图中由左下方折向右上方的带编码a、b、c、d、e、f的每一根折线表示

一辆车在时间空间二维平面上的运动轨迹。折线的斜率表示该车的速度，其两折线之

间的水平距离表示相应两车的车头时距中间可能还有其它车辆)，纵向距离为相应两

车车头的空间距离。车流中密度不相同的两部分的分界称为集散波，众所周知的波速

公式为：

Ws=(Qx-Qy)/(K,-K,)

有如下公式；：

K,=0.5储(1-1-2)

s

降=0.5储(1+1-立)

s

式中：勺如前文为车流停车时的最大密度；K；为畅流流量。，所对应的密度；K.为拥

挤流量S,所对应的密度。

3系列公式

(1)图1下图是以集散波理论为基础，对照图1上图画出来的，现以tp和xp分

别表示图1中P点的时间坐标和位置坐标。在波速公式WQ,⑨中，可省略Q,只留

下下标i,简记为叱应，如果。尸0,简记为W，，若S,等于零，简记为叱，密度

为K,简记为K,。

(2)可求出最远距离：

XD=[-](3)以公式(1)

%2,S%,S

以及波速公式代入式(3),可推导出

X_(5一§2)[(凡+&+&)Qi-R4-&S2](4)

D

~KS2(Qi-S)+Kl(S-S2)-^Ks(S2-Qi)

Y

j=a+&+&+——(5)

…-年

t,=tF-xF/vie

如果T^t,,且Q2<Q{<S2,则F是最远点，通过图8可得

x-KN+&)。1-RISJS2(6)

F-KJQfl+KM-KszQ

V

tF=R}+R2+777-(7)

%2

若%<7；4勺,且Q|>S2>Q2，最远点为

x-。冏(or)-夫用+国+/?2电]⑻

K£-KS2Qt

若7；%，且2>Q2>立则最远点为

X=(S-S2)%。—R5-R3s2+(K+&+R「7)1。2

Q)

一KS2(Q1-S)+Ki(S-S2)+KS(S2-

若工W4,。>S?>Q2,则最远点为

X_SzIZQi-R\S\+(/?]+/?2-7])02电(]0)

Kj(Q2—S2)+K2S2—KS2Q2

信号交叉口的信号周期C=60s,相位时间均为30s,黄灯时间为3s。饱和流量S=1700,

辆/h,到达流量在红灯前段9s内为Ql=2300辆/h,在同期内其余时段为Q2=

1080辆/h,假设Q-K模型曲线为二次抛物线，可算出

4=旦=42.405辆/碗

K,=2=35!^/km

v2

Ks=50辆球机

于是由式（4）及式（5）得

XD=3.211928WS2J=一36左m/〃，ro=0.17578Atj=0.116562

由于工"，。＞$2,符合公式（9）的条件，所以

代入数据得X=-0.16432km

由实际情况可知，排队长队为164米，与实际符合较好。

5.4问题四的解答:

由表3数据可知，通过车流量为1200pcu/h。问题四假设路段下游方向需求不变，路段

上游车流量为1500pcu/h且保持不变，事故发生时车辆初始排队长度为零，且事故持续

不撤离。由此可以定性的看出排队是不可避免且队伍越来越长。对于车流波动理论，T

表示集结与消散两波相遇的时刻，x表示最大排队长，由车流波动理论不能直接求得队

伍排到140所需的时间。可以假设队伍在排到140米之前不采取措施，保证事故不撤离。

前段队长的增加是比较均匀的，则每秒队伍增加长度

Ax=(7.55x1000)/(38.4x60)=3.3m/s

则排到140米处需要的时间t=140/3.3=42.4s。

此结果视频实际情况偏差较大。考虑另一个模型，流体模型。

-5

其中，排队车量

。“集结车量

(2口消散车量

Qw>0时排队，2。40时，不排队

由题目要求可知上游车流量为1500pcu/h,则。”=25辆/分，通过车流量为1200

pcu/h,则。2=20辆/分。由于红绿灯的影响，消散过程是连续的，集结过程是不连

续的。

根据问题三可知，阻塞密度K,=110辆/km,阻塞即车辆无法通行时车辆密度。

由此可以得到车辆排队140m时对应的车辆数n=16辆。假设车辆不计较均匀的分布

在三条车道上，根据图3的事故截面示意图，参与排队的车辆数为

N=l+3x(16-1)=46辆。对于队尾，只要有一条道路达到140米即算队长为140米,

此时是最少车辆数为46-2=44辆。故达到140米的排队车辆数属于(44,46)

在第一分钟，假设前30秒绿灯，后30秒红灯。则在第一分钟内能够全部消散

全部车辆，即不发生排队。在第二分钟之后，每分钟有5辆车参加排队。则在第10

分钟就有45辆车参与排队，在第11分钟有50辆车参与排队。故在(10,11)分钟

之间排队长度达到140米。

参