数学建模国赛国家二等奖优秀论文

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

承 诺 书

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

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

我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。

我们知道，抄袭别人的成果是违反竞赛章程和参赛规则的，如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。

我们郑重承诺，严格遵守竞赛章程和参赛规则，以保证竞赛的公正、公平性。如有违反竞赛章程和参赛规则的行为，我们将受到严肃处理。

我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。

我们参赛选择的题号是（从A/B/C/D中选择一项填写）： A 我们的参赛报名号为（如果赛区设置报名号的话）： 所属学校（请填写完整的全名）： 参赛队员 (打印并签名) ：1. 刘冲 2. 3. 指导教师或指导教师组负责人 (打印并签名)

日期： 2013 年 9 月 16 日

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

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

编 号 专 用 页

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

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

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

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

车道被占用对城市道路通行能力的研究

摘要

本文就交通事故对通行能力的影响进行分析研究，主要对实际通行能力的变化、排队长度、事故持续时间、交通流量等问题建立相应的数学模型，并运用

SPSS、MATLAB等软件工具对模型求解。

针对问题一，首先对视频一进行数据采集和提取，利用插值法对缺失数据进行补充。然后以基本通行能力、可能通行能力为基础，综合考虑外界动态因素，构建出“合流难度系数”模型，进而得出实际通行能力的函数式，由此详细地描述出事故横断面处实际通行能力的变化过程。

针对问题二，首先应用配对样本t检验法得出所占车道不同对通行能力的确存在显著性差异的结论。然后构建出视频二中的实际通行能力函数，与问题一的函数进行对比分析。再结合综合分析模型，从不同车道的车流量、拥堵车道的车流容量以及拥堵时间比例等角度进行对比，分析出差异产生原因在于：各车道车流辆不同导致合流密度不同，合流密度越大，换道难度越大，通行能力下降越多。

针对问题三，首先构建理想条件下的“到达—离开模型”，构建出车辆排队长度与实际通行能力、事故持续时间、路段上游车流量之间的关系；其次，引入交通波理论，构建出“车流波动理论模型”；最后结合交通信号灯对交通流有周期性影响的实际情况，建立“基于二流理论的动态排队模型”，得到在一个周期内对长的相对增量，再通过累加得出车队长的表达式。

针对问题四，考虑小区进出车辆的影响，以及在更高车流量下合流系数的改变，对上述模型参数做出修正，估算出排队时间大约为7.3分钟。接着应用“基于元胞自动机的交通模型”进行仿真模拟，得出误差在接受的范围之内，证明了模型的合理性。

最后，针对模型的优缺点进行评价，并提出了进一步改进的优化方向及在其它方面的推广应用。

关键词：实际通行能力、到达—离开模型、交通波、二流理论动态排队模型

1

一、问题重述

车道被占用是指因交通事故、路边停车、车道施工等因素，导致道路横断面通行能力在单位时间内降低的现象。由于城市道路具有交通流密度大、连续性强等特点，一条车道被占用，可能降低路段所有车道的通行能力，即使时间短，也可能引起车辆排队，出现交通堵塞。

车道被占用的种类繁多、复杂，正确估算车道被占用对城市道路通行能力的影响程度，将为交通管理部门正确引导车辆行驶、审批占道施工、设计道路渠化方案、设置路边停车位等提供理论依据。 根据题意，本文需要解决的问题有：

（1）描述视频一中从交通事故发生到撤离期间内，横断面内实际通行能力的变化过程。

（2）根据问题一的结论并结合视频二，说明同一横截面交通事故所占车道不同对该横截面实际通行能力的影响差异。

（3）构建恰当的数学模型，分析车辆排队长度与实际通行能力、事故持续时间、路段上游车流量的关系。

（4）先给定交通事故所处横断面距离上游路口为140米，路段上游车流量

1500pcu/h；并且假定事故发生时车辆初始排队长度为零，事故持续不撤离。计

算车辆排队长度到达上游路口时所用的时间。

二、问题分析

在车道被占用对通行能力的影响中，车道被占用的情况是多种多样的，本题只考虑交通事故对其横断面实际通行能力的影响。

对于问题一，要想描述通行能力的变化过程，必须先定义出“实际通行能力”的概念：单位时间内通过道路上某一均匀段或者某一横断面上的最大车辆数。再结合基本通行能力(C0)、可能通行能力(C1)，外加外界动态因素f1(t)的影响给出实际通行能力随时间的变化函数。然后从视频一中采集与函数相关的数据，画出函数的变化图象。通过观察图像的变化情况，外加文字叙述描述该变化过程。

2

对于问题二，在问题一函数的基础上，再结合视频二中的有用信息，写出实际通行能力的表达式C22?C1?f2(t)。然后利用SPSS软件对两组实际通行能力值进行配对样本t检验，得出两者间确实存在显著性差异。接下来为体现两者间的明显差异，把两个函数图象画到同一坐标系中进行比较。最后，根据统计得到的数据分析此差异形成的原因。

对于问题三，以问题一、二为基础，首先说明车辆排队长度与实际通行能力、事故持续时间、路段上游车流量确实存在某种关系。然后先构建理想化的模型体现四个量之间的关系；其次，把模型的影响因素考虑进去，对模型进逐渐修正，得到更为精确的模型。

对于问题四，观察附件五，当排队距离变为140米时，交通事故正好发生在小区门口，所以此时应考虑小区进出车辆的影响。接着对问题三中的模型参数做微型修正。再结合已知的排队长度、上游车流量估算出相应时间。

三、模型假设

1. 在研究路段内，路边停车现象对通行能力的影响忽略不计。 2. 假定具有良好的气候条件和路面条件下的通行能力。 3. 该段时间内分析的道路及道路上的基础设施不发生变化。 4. 假设司机反应时间相同且遵守交通规则。 5. 假设上游汽车当量服从离散型随机分布。

四、符号说明

C?? 通行能力；

Qa??规定时间周期内驶出交通事故横断面的车辆数；

Qb??规定时间周期内驶入标准路段长度内的车辆总数； Ki??各汽车代表车型与车辆折算系数，(i?1,2)；

f?(t)??交通事故中车辆改道对通行能力的合流难度系数；

g(LN)?? 不同车道编号机器组合；

3

6.2.2 模型建立

首先从视频二中采集与问题一中相类似的数据信息，结合模型一中建立的“合流难度系数”，构建视频二中的实际通行能力模型，模型如下：

C22(t)?C1?f2(t) （5）

Qa??KnPnQb??KnPni?1i?122f2(t)? （6）

首先，结合上述函数画出在“交通事故发生到事故撤离期间内”交通事故横断面上实际通行能力变化曲线，图形如下图4所示：

图4 视频二通行能力变化曲线

其次，再根据 C21?g(LN)?g(2,3)与C22?g(LN)?g(1,2)，将两个视频中“合流难度”系数的变化趋势画在一起进行比较，并分析两者之间的显著性差异，如下图5所示：

图5 合流难度系数变化比较

之所以存在显著性差异，是因为每一车道车的数量不同，导致在车辆换道时的难易程度不同，下图为在交通事故期间不同车道的平均流量图6：

9

图6 不同车道平均流量图

说明：由上图6可以看出，中间车道和右车道的车辆总数相对于左车道多，且公交车的数量右车道相对较多。所以当中间车道和右车道被堵时，两车道上原来行驶的车辆需要换道。考虑到公交车的长度对换到有很大影响，所以此时车道

Q2、Q3对车道Q1的合流密度会增大，使得通行能力下降比较明显；而Q1、Q2相对于Q3需要换道的车辆及公交车数比较少，所以对应的合流密度比较小，对实际通行能力的抑制作用比较小。

从而，造成了比较明显的差异。

6.3问题三分析与建模

根据问题一、二中的视频数据及函数分析，车辆排队长度与事故横断面实际通行能力、事故持续时间、路段上游车流量确实存在着某种函数关系，于是建立以下模型进行求解。 6.3.1累计到达——离开模型

累计到达一离去模型[3]是在排队论的基础上演化而来，应用车辆到达和离去之间的关系，主要通过累计到达的车辆数和累计离去车辆数之间的差值来估算排队长度。

在某一时间段内，从上游到达的车流量为?Qbdt，而在同一时间段内流向下游的当量为?Qadt，由于上下游车流量的不同，从而在时间段0~t内就会有车辆积累。随着时间的积累形成车辆排队现象。在理想状态下，车流会依次等距排队，因此可以建立以下车辆排队长度与其它量的函数关系：

???Qadt??Qbdt?Q? （7） ???3L?Q?(lt?lh)

10

Q???在标准时间的滞留车辆； lt??车的平均长度； lh??车与车间的距离； L??车辆排队长度。

图7 到达、离开车辆与时间函数

由于此模型是在上述比较比较理想情况下得出的，假设一条车道上所有的车辆走完之后相邻车道上的车紧跟其后；车在车道上是均匀分布的，且认为向下游的车流流量就是交通事故横断面的实际通行能力。该模型与实际问题相差甚大，接着修正模型模型使得相对误差较小。 6.3.2车流波动理论模型

当车流在交通事故界面处发生拥挤时，后面来的车辆要依次减速行驶或者停车，这种影响会逐渐向后传播，直到最后一辆车明显减速或停车，就像水面上的一个水波不停地向周围传播，因此可以将水波的传播类比到堵塞车辆的依次停车现象中，最终建立车流波动理论来研究车辆排队长度问题。

车流波动理论[4]：当行驶车辆到达事故地点时速度陆续减慢甚至停车而集结成密度较高的队列；事故消除后，由于路段通行能力的恢复，排队车辆又陆续加速而疏散成一列具有适度密度的车队，车流中两种不同密度部分的分界面经过一辆辆车向车队后部传播的现象。波速公式：

u?(Qx?Qy)/(Kx?Ky) （8）

u?? 集散波的波速；Qx、Qy??前后两种车流状态的流量，辆/时；

Kx、Ky??前后两种车流状态的密度，辆/Km.

根据交通流模型可知，交通量Q、行车速度v、车流密度K三着之间的关系为：

11

Q?v?K （9）

速度-密度线性关系模型：

v?vf(1?K/Kj) （10）

由上面（8）（9）（10）推导出波速与密度的函数关系式，即波速与上游车流量、横断面实际通行能力的关系式：

u?vf(1?Kx?KyKj) （11）

其中Kx、Ky与Qx、Qy有关，Kx就是上游车流量，而Ky就是下游的车流量。当集散波向后传播时，车道就可以认为处于最大流量状态，即等于事故横断面处的实际通行能力；从而得出：

Kx?Ky?u?v(1?)?fKj? （12） ?t?L?udt??0?

6.3.3 基于二流理论的动态排队模型

在波动理论中，车流波的传播会受到很多情况的影响而发生停止传播现象。而且实际道路中，每次交通信号灯周期性的变化，都会导致每次流入车流的不均衡现象，因此为了能够更好地描述每个周期内车辆排队长度，结合二流理论的思想，提出了动态排队模型。

图8 实际情况下的变化曲线

根据二流理论思想[5]，将运动车辆形成的交通流称为行驶交通流，停止车辆形成的交通流称为阻塞交通流。这样，交通流实际运行状态中过渡状态的不均匀

12

交通流相当于阻塞交通流和行驶交通流的某种加权和，即任意交通流的实际运行状态可以用二流运行状态来描述。

图9 周期与?t的关系

以上面图9为依据，根据流量守恒原理可知：

N0?NU(t)?ND(t)??N(t) （13）

N0??初始时刻（即事故即将发生时）上游、事故横断面间的车辆数；

NU(t)、ND(t)??t时刻通过上游、横断面的车辆累计数；

?N(t)??t时刻上游与横断面间的车辆数。

结合二流理论：

?N(t)=kjLD(t)?km[L?LD(t)] （14）

联立（13）（14）得到单车道路段当量排队长：

LD(t)?[N0?NU(t)?ND(t)?kmL]/(kj?km) （15）

以此为基础，得到多车道路段平均当量排队长度模型;

N0??NU(i,t)??ND(i,t)?kmLMi?1i?1MM?LD(t)?M(kj?km)?? （16）

LD(t)??多车道路段时刻t上游、横断面间的平均当量排队长度； NU(i,t)、ND(i,t)??第i条车道时刻t上游、横断面的车辆累计数；

M??车道数。

令时间t?t0得到此刻平均当量排队长度LD(t0)，

13

LD?t0????N0????NU?i,t0??ND?i,t0???kmLMi?1M???M?kj?km????? （17）

同理，当t?t0??t时，平均当量排队长度为

LD(t0??t)?{N0??[NU(i,t0??t)?i?1?M （18）

??ND(i,t0??t)kmLM}\\[M(kj?km)]?在t0时刻，时间增量?t引起的平均排队长度增量为：

?LD(t0)??[NU(i,t0??t)?NU(i,t0)?i?1M （19）

ND(i,t0??t)?ND(i,t0)]\\[M(kj?km)]同理，?t时间内上游、横断面车辆累计数的增量为：

M?M??[NU(i,t0??t)?NU(i,t0)]??QU(i,?t)?i?1i?1 （20） ?MM?[N(i,t0??t)?N(i,t0)]?Q(i,?t)?D?DD?i?1?i?1??事故横断面的车辆数。 QU(i,?t)、QD(i,?t)??第i条车道?t时间内通过上游、将上述(20)式带入（19）式得到：

??LD(t0)??Q(i,?t)??Q(i,?t)abi?1MMM(kj?km)M?M?Q(i,?t)?Q(i,?t)?b??a?i?1i?1? （22） L??????n?1?M(kj?km)????tn?T?i?1? （21）

在这里由于下游车道数只有一个，上游车道数为三个，则最终可以得出

?3?tn??Q(i,?t)?(1,?t)????Q??a?T?t??Q??t???bQb? （23） ????aL???i?1??n?1?Mkj?kmn?1?Mkj?km??????tn?T????n??周期的个数，其中T?60s，与红绿灯的周期相等。

14

6.4问题四分析与建模

（1）针对问题三中的排队时间问题，此时的上游车流视作是稳定的，所以可以根据模型对排队时间进行估算。首先根据排队长度函数关系 ：

???Qadt??Qbdt?Q? （24） ???3L?Q?(lt?lh)把

Qa?1500pcu/h代入到式中，得到模型一的估算时间t1?7.27min。

（2）基于已经建立的车流波动理论模型，首先要得到出上下游车流密度，由于在排队过程中，下游车流一直处于在最大通行能力范围内，并且路段下游方向需求不变，就可以利用视频一中的下游实际通行能力进行估算，上游通行能力将保持为25pcu/min。于是根据车流波动理论模型：

Kx?Ky?u?v(1?)?fKj? （25） ?t?L?udt??0?估算得出车辆排队时间t2?6.29min。

（3）在基于二流理论的动态排队模型中，队长的变化是由每一个周期队长的增量累加得到的。在实际交通中的每个周期内，上下游截面的车流状况是不断变化的，而在本问题中已知上游车流量是固定的，于是就要对本问题的数据进行转化。首先，还以60s为一个周期，认为每个周期的上游车流量为

nQa?25pcu/mi，下游横断面的车流量与问题一中实际通行能力相等，再根据动态排队模型，得到了排队时间的函数表达式：

Kx?Ky?)?u?vf(1?Kj? （26） ?t?L?udt??0?

经过计算，得到了估算时间t3?7min。将三个模型的时间估算结果如下表4：

15

表4 各模型对应的时间

模型 估算时间/(min)

模型一 7.27 模型二 6.29 模型三 7 从表4中可以看出，三个模型对时间的估算是不一样的，主要是由模型算法前提条件不同而造成的。接着利用“基于元胞自动机的交通模型”进行模拟模拟，得出的误差在接受的范围之内，说明了模型的合理性。

七 模型评价

7.1 模型优点

（1）对于问题一，利用题目信息构建的实际通行能力变化函数模型充分考虑了外界因素的作用并合流难度系数，使得相对误差相对较小，能很好的反应通行能力的变化过程；其次，利用插值法补全视频一中缺失的数据，更具有准确性。

（2）对于问题二，用配对成份t检验法有效的验证了两组通行能力之间的显著性差异，算法简便。其次，对两个函数变化进行对比时形象直观，能直接反映出占道不同对通行能力的抑制性不同。

（3）对于问题三，累计到达-离去模型可以用流量来计算，而实际道路网流量相对比较容易获取；而且模型比较灵活，方便根据城市道路特性进行改进。 基于二流理论动态排队长度模型综合考虑了多个动态因素的共同作用，使得结果的可靠性提高。 7.2 模型缺点

（1）在构建实际通行能力函数模型时，有的阻滞系数都采用定值或平均值，对预测结果会产生影响。

（2）在到达一离去理想化模型的运用中，进行了很多假设，将多个外界因素影响排除在外，使得与实际函数之间存在些许误差。

16

八、模型推广

（1）到达—离开模型在实际中还可以用于处理动态的行人流、社会流、服务流等，特别是在候车室、候机楼客流的预测中具有实际的意义。

（2）综合分析是一个最广泛应用的统计方法，在数据处理过程中能发挥有效的作用，并为实际问题的分析提供了有说服力的手段和依据。

（3）二流理论可应用于评价宏观道路网络交通流特点，为路网的合理设计和交通相位的最优分配提供重要的理论基础。在合理利用资源方面二流理论也有潜在的研究价值。

九、参考文献

[1] 韩中庚，数学建模竞赛——获奖论文精选与点评，北京：科学出版社，2007 [2] 贾晓敏，城市道路通行能力影响因素研究，公路，2009.

[3] 孔惠惠，秦超等，交通事故引起的排队长度及消散时间的估算，铁道运输与

经济，第27卷 第5期：65-67，2005.

[4] 姚荣涵，王殿海，拥挤交通流当量排队长度变化率模型，交通运输过程学报，

第9卷 第2期：93-96,2009.

[5] 王彦杰，徐建闽等，城市交叉口交通拥挤最大排队长度估算模型，行业管理，

86-87,2009.

17

十、附录

1.MATLAB程序

X=[1 2 3 4 5 6 7 8 9 10];

Y=[19 16.5 20 18 18 20 16.5 17.5 15.5 16.5]; my=mean(Y) %Y的平均值 sy=std(Y) %Y的均方差(标准差)

M=[1.5 2.5 3.5 4.5 5.5 6.5 7.5 8.5 9.5]; N=interp1(X,Y,M,'cubic') P=polyfit(X,Y,8) polyval(P,[10 11])

plot(X,Y,'r*',M,N,'g+',X,polyval(P,X),'m-'); legend('原始点','插值点','拟合曲线')

2. 为了找到这条道路上的当量车长和平均车距和的大小，经过对视频一中达到最大排队长度时的排队长和当量车流数的数据采集得到了下表1：

表1

队长90 /m 当量8 数 90 60 90 90 60 8 7 10 10 7 3.

表二 视频一中折减系数 时间 小型汽车数 驶入 驶出 中型汽车数 驶入 驶出 大型汽车数 驶入 驶出 驶入 车辆 驶出 车辆 折减 系数 事故前16:42:01 12 16:43:01 13 16:44:01 14 18 17 12 12 12 3 3 2 2 4 4 5 3 3 3 18

0 0 1 0 2 0 1 0 1 1 16.5 17.5 19 12 24 24 25 1.47 1.46 中 16.5 0.86 18.5 1.5 故16:45:01 9

16:46:01 14 18.5 0.79

16:47:01 15 16:48:01 21 16:50:01 16 16:51:01 20 16:52:01 16 16:53:01 16 16:54:01 9 14 14 15 13 12 13 14 11 17 3 5 4 2 4 3 2 4 3 4 3 3 3 3 3 4 3 4 1 0 1 1 1 1 1 0 0 2 0 0 2 1 1 0 1 0 21.5 20 1 28.5 19.5 0.65 24 25 24 19.5 0.83 21.5 0.76 18.5 0.77 22.5 19.5 0.85 14 21 17.5 20 1.44 16:59:01 15 事故后17.5 0.83 21 1.2 15:00：01 13 4.

5.

19

6.

7.

20

本文来源：https://www.bwwdw.com/article/3iud.html