桂水榮 藍天飛 陳水生 葛世祺

摘要： 為研究交通事故對高速公路交通流的影響，以雙車道安全距離元胞自動機模型為基礎(chǔ)引入強行變道規(guī)則，建立了考慮車輛異質(zhì)性的高速公路車輛堵塞模型。利用MATLAB數(shù)值模擬高速公路異質(zhì)交通流，分析不同交通密度情況下的交通事故持續(xù)時間以及重載貨車占比對高速公路堵塞交通流的影響。研究表明：交通事故持續(xù)時間及重載貨車占比均對高速公路堵塞交通流具有顯著影響；當車輛密度大于4veh/km小于80veh/km時，交通事故持續(xù)時間對交通流有著顯著影響，且持續(xù)時間越長對交通流的影響越大；當車輛密度小于30veh/km時，重載貨車占比對交通流有著顯著的影響，并且占比越大對交通流影響越大。

關(guān)鍵詞： 異質(zhì)交通流；元胞自動機；模型仿真；堵塞

中圖分類號： U491.4；U491.1文獻標識碼： A

收稿日期： 2021-09-17；修回日期：2022-01-03

基金項目： 國家自然科學基金（52268026）；江西省自然科學基金（20181BAB206041，20181BAB206043）

第一作者： 桂水榮（1979-），女，江西九江人，博士，副教授，主要研究方向為公路橋梁車橋耦合振動。

On Congestion Traffic Flow Model Considering Vehicle Heterogeneity

GUI Shuirong，LAN Tianfei，CHEN Shuisheng，GE Shiqi

（School of Civil Engineering and Architecture，East China Jiaotong University，Nanchan 330013， China）

Abstract：In order to study the impact of traffic accidents on expressway traffic flow， a forced lane change rule was introduced based on the two-lane safe distance cellular automaton model， and a highway vehicle congestion model considering the heterogeneity of vehicles was established.The heterogeneous traffic flow on expressways was used by MATLAB numerical simulation， the impact of the duration of traffic accidents and the proportion of heavy trucks under different traffic densities on expressway congestion traffic flow was analyzed. Research shows that the duration of traffic accidents and the proportion of heavy trucks have a significant impact on the traffic flow of highway congestion; when the vehicle density is greater than 4veh/km and less than 80veh/km， the duration of the traffic accident has a significant impact on traffic flow. And the longer the duration， the greater the impact influences on the traffic flow; when the vehicle density is less than 30veh/km， the proportion of heavy trucks has a significant impact on the traffic flow， and the larger the proportion， the greater the impact influences on the traffic flow.

Key words： heterogeneous traffic flow; cellular automata; simulation model; congestion

0 引言

高速公路作為一種現(xiàn)代化的交通運輸方式，對社會經(jīng)濟發(fā)展有著巨大的促進作用，高速公路上的交通流在絕大部分時間里均為混合流，行駛在高速公路上的汽車一般可以分為小汽車及重載貨車，盡管重載貨車數(shù)量遠小于小汽車，但研究表明重載貨車對交通流具有顯著的影響［15］，隨著社會發(fā)展，交通擁堵問題也愈發(fā)明顯。引起高速公路交通產(chǎn)生擁堵的原因多種多樣，其中交通意外事件造成的擁堵居首位。交通意外事件會降低道路通行能力，造成交通混亂，嚴重時會引起新的意外事件。因此分析擁堵狀態(tài)下高速公路交通流特點十分必要。王永明［6］在現(xiàn)有的元胞自動機模型中引入強制換道規(guī)則，模擬了交通擁堵時車流的傳播過程，探索了主要交通流參數(shù)對堵塞傳播的影響規(guī)律。錢勇生等［7］在雙車道交通流元胞自動機中引入意外事件，分析了車道管制下高速公路發(fā)生意外事件對交通流的影響。劉霞等［8］分析了不同事故區(qū)域的換道特點，建立了考慮事故發(fā)生地點的雙車道元胞自動機模型。姬浩等［9］為了探究事故車輛對城市三車道道路交通的影響機制，建立了考慮主動搶道與被動搶道的元胞自動機模型。石京等［10］考慮了多種事故因素聯(lián)合作用，將多種事故因素加入到元胞自動機模型中，計算了多種因素疊加產(chǎn)生的交通事故風險。然而，上述研究均未考慮堵塞交通流異質(zhì)性的特點。

本文考慮交通流異質(zhì)性特點，建立交通事故影響下雙車道元胞自動機模型，分析交通事故堵塞持續(xù)時間及重載貨車占比對車速及車流量的影響，以期為管理部門交通事故管控提供依據(jù)。

1 雙車道安全距離模型

雙車道安全距離模型是基于單車道安全距離模型［1113］考慮了車輛間距、車輛性能以及相鄰車輛車速的雙車道元胞自動機模型。在該元胞自動機模型中，把車輛大致分為小汽車及重載貨車，道路分為2×N個元胞，其中N為單車道元胞個數(shù)。小汽車的長度為Lc，重載貨車的長度為Lt。設小汽車的最大速度為Vmc，重載貨車的最大速度為Vmt。小汽車一般加速度為ac，重載貨車一般加速度為at。小汽車的緊急制動加速度為amc，重載貨車的緊急制動加速度為amt。小汽車的慢化概率為Rc，重載貨車的慢化概率為Rt。

1.1 跟弛規(guī)則

1.1.1 安全距離計算

汽車加速、減速及保持當前速度前進的前提需有足夠的安全距離。為了使汽車不發(fā)生碰撞，在更新汽車速度之前需進行安全距離判斷。安全距離分為3種：汽車加速前進的安全距離daccn；汽車減速前進的安全距離ddecn；汽車勻速前進的安全距離dkeepn。通常情況下，元胞自動機交通流模型安全距離計算方法為

dsafe=sn+s′n-sn-1=v2n2amn+vnτ-v2n-12amn-1（1）

其中，sn為后車的制動距離，s′n為后車反應時間τ內(nèi)行駛的距離，sn-1為前車制動距離。但是，式（1）并不完全正確。若前后兩車加速度不同，兩車的最短距離有可能出現(xiàn)在并非兩車都停止的狀態(tài)。因此，異質(zhì)交通流需考慮車輛行駛速度減速至零時，每一個時間段內(nèi)兩車的距離變化。以daccn為例：

1）加速安全距離daccn。假設t時刻時，車距為dn（t），后車以速度（vn（t）+an）進行位置更新；同一時刻內(nèi)，前車速度為vn-1（t），因突發(fā)情況以加速度amn-1進行緊急制動，前車以（vn-1（t）-amn-1）進行位置更新，并持續(xù)進行緊急制動直至車輛停止。在t+1時刻，后車為避免追尾也進行緊急制動直至車輛停止；假設后車在（t+Δtn）時刻停止運動，前車在（t+Δtn-1）時刻停止運動，在（t+Δt）時刻兩車都停止運動；兩車在緊急制動直至停止過程中，恰好沒有發(fā)生碰撞；則t時刻的車距dn（t）則為后車在t時刻加速運動的臨界安全距離daccn。根據(jù)以上分析過程，可以計算daccn。首先，Δtn與Δtn-1計算表達式為

2 系統(tǒng)仿真

仿真中車輛先正常行駛，然后在某個時間步出現(xiàn)交通事故產(chǎn)生擁堵，隨后交通事故消除，研究交通堵塞情況下的堵塞時間及重載貨車占比對交通流量、速度以及車輛換道率的影響。

2.1 參數(shù)設置

元胞參數(shù)：每個元胞長度為0.5 m，每個車道元胞數(shù)為2 000個，即車道長為1 000 m，共兩個車道。其中小汽車長度為15個元胞（7.5 m），重載貨車長度為30個元胞（15 m）；小汽車最大速度為67個元胞（約為120 km/h），重載貨車最大速度為45個元胞（約為81 km/h）；小汽車加速度為5個元胞（2.5 m/s2），重載貨車加速度為3個元胞（1.5 m/s2）；小汽車緊急制動速度為8個元胞（4 m/s2），重載貨車緊急制動速度為5個元胞（2.5 m/s2）；小汽車隨機慢化概率取0.2，重載貨車隨機慢化概率取0.1。

仿真參數(shù)：每次仿真運行13 600個時間步，選取最后3 600個時間步的車輛數(shù)據(jù)。為了消除隨機性影響，每種情況運行10次，取10次數(shù)據(jù)的平均值。

事故參數(shù)：在時間步為10 000時，1車道的第1 000個元胞前方的最近車輛發(fā)生了交通事故，導致1車道無法通行；并且事故時間持續(xù)了h個時間步。

2.2 堵塞持續(xù)時間對堵塞交通流的影響

分別令堵塞持續(xù)時間h=0 s，h=600 s，h=1 200 s，h=1 800 s，h=2 400 s，h=3 000 s，h=3 600 s，重載貨車數(shù)量占總車輛數(shù)量比例k=0.1，研究不同堵塞時間下，交通流量、速度、換道率隨交通密度的變化規(guī)律。

2.2.1 堵塞持續(xù)時間對交通流量的影響

圖3為隨著堵塞持續(xù)時間增加，堵塞持續(xù)時間與交通流量的流量—密度圖。從圖3可以看出，當ρ≤4veh/km時，不同堵塞持續(xù)時間對應車流量無區(qū)別。在低密度情況下，堵塞持續(xù)時間對混合車流沒有影響，車流可以順暢行駛。當4veh/km<ρ<80veh/km時，不同堵塞時間對混合車流影響較大，在相同密度時，車流量隨著堵塞時間h的增加而減小。在該密度范圍內(nèi)，堵塞持續(xù)時間對混合交通流產(chǎn)生了較大的影響。當ρ≥80veh/km時，不同堵塞持續(xù)時間對應車流量無明顯區(qū)別，此時車輛密度較大，車輛間的相互作用增強，從而交通事故引發(fā)的堵塞對混合交通流影響不大。在高密度情況下，堵塞持續(xù)時間對混合車流影響較小。

2.2.2 堵塞持續(xù)時間對車速的影響

圖4為隨著堵塞持續(xù)時間逐漸增大，堵塞持續(xù)時間與平均速度的速度—密度圖。從圖4可以看出，當ρ≤4 veh/km時，不同堵塞持續(xù)時間對應速度無區(qū)別。在低密度情況下，堵塞持續(xù)時間對混合車流沒有影響，車流可以暢通行駛。當4 veh/km<ρ<80 veh/km時，堵塞時間h=0時，車道上無事故發(fā)生，車速隨著車輛密度的增大而逐漸減小，說明無事故時，車速會隨著車輛密度的增大而減?。辉谙嗤芏葧r，車速隨著堵塞時間的增加而減少，說明在該密度區(qū)間下，堵塞持續(xù)時間對車速影響明顯；盡管堵塞持續(xù)時間不同，但7種情況的車速隨著密度增加都趨于同一個車速，說明堵塞時間對車速影響，隨著車輛密度的增大而減小。當ρ≥80 veh/km時，車速并沒有隨著堵塞持續(xù)時間的不同而發(fā)生變化。在高密度情況下，堵塞持續(xù)時間對混合車流影響已不是很顯著。

2.2.3 堵塞持續(xù)時間對換道率的影響

圖5為分別隨著堵塞持續(xù)時間逐漸增大，重載貨車及小汽車的換道率—密度圖。從圖5可以看出，相同情形下，小汽車的換道率遠高于重載貨車，這是由于重載貨車的換道條件更嚴格且重載貨車司機的換道欲望比小汽車司機更低。重載貨車與小汽車換道率的最大值隨堵塞持續(xù)時間的增加而增大，這是由于在車輛密度不大的情況下，隨著堵塞持續(xù)時間的增加會有更多的重載貨車與小汽車為了追求更快的速度從事故車道變道。在密度較高時，重載貨車與小汽車的換道率會隨著車輛密度增大而減小，這是因為，隨著車輛密度的增大，車輛之間的間隙變小，車輛之間沒有足夠的安全空間進行換道。

2.3 重載貨車比例對堵塞交通流的影響

分別令重載貨車比例k=0，k=0.1，k=0.2，k=0.3，k=0.4，k=0.5，k=0.6堵塞持續(xù)時間h=600 s，研究不同重載貨車比例下，交通流量、速度、換道率隨交通密度的變化規(guī)律。

2.3.1 重載貨車比例對交通流量的影響

圖6為隨著重載貨車比例逐漸增大，重載貨車比例與交通流量的流量—密度圖。由圖6可知，當車流中只有小汽車時平均車流量最大，其中最大流量約為1 772 veh/h，而重載貨車比例為0.5時平均車流量最小，其中最大流量約為1 379 veh/h。在堵塞持續(xù)時間不變的情形下，隨著重載貨車比例的增大，平均車流量逐漸減少。這是因為，與小汽車相比，重載貨車占用的元胞數(shù)量更多、道路空間更大，車輛行駛靈活性更低，隨著重載貨車比例的增大，重載貨車數(shù)量的增加降低了道路的通行能力，小汽車只能以較低的速度行駛。

2.3.2 重載貨車比例對車速的影響

圖7為隨著重載貨車比例逐漸增大，重載貨車比例與平均速度的速度—密度圖。由圖7可知，當ρ≤30 veh/km時，重載貨車比例對車流平均速度具有明顯影響，在堵塞持續(xù)時間不變的情形下，車流平均速度隨著重載貨車比例的增大而減小。重載貨車的加入限制了本車道后車的速度，使得本車道后車不得不采取減速或者換道，以保持兩車安全距離。當ρ≥30 veh/km時，重載貨車比例對車流平均速度的影響不再明顯，在堵塞持續(xù)時間不變的情形下，車流平均速度隨著重載貨車比例的增大而減小幅度不大。這是因為，隨著車流密度的增大，車輛之間的相互作用更加明顯，重載貨車對車流影響減小。

2.3.3 重載貨車比例對換道率的影響

圖8為隨著重載貨車比例逐漸增大，重載貨車及小汽車的換道率—密度圖。由圖8可知，當ρ≤30 veh/km時，在堵塞持續(xù)時間不變的情形下，重載貨車的最大換道率隨著重載貨車比例的增大而增大。這是由于，在低密度情況下，道路空間并不會隨著重載貨車數(shù)量的增加而產(chǎn)生很大的變化，當?shù)缆烦霈F(xiàn)交通事故時，重載貨車與其相鄰車輛之間仍有足夠的空間進行換道。在堵塞持續(xù)時間不變的情形下，當重載貨車比例k≤0.3時，小汽車的最大換道率是隨著重載貨車比例的增大而增大；當重載貨車比例k>0.3時，小汽車的最大換道率是隨著重載貨車比例的增大而減小。這是由于，在低密度情況下，小汽車的換道需求隨著重載貨車比例的增大而增大；但是當重載貨車比例超過臨界值時，重載貨車占用的空間較大，提升了小汽車的換道難度。當ρ>30 veh/km時，在堵塞持續(xù)時間不變的情形下，小汽車與重載貨車的換道率均隨重載貨車比例的增大而減小。這是因為，隨著重載貨車比例的增加，重載貨車占用的空間也越來越大，使得小汽車與重載貨車換道變得愈加困難。

3 結(jié)論

本文基于雙車道安全距離元胞自動機模型，引入了交通事故，模擬了高速公路異質(zhì)交通流發(fā)生事故時，堵塞持續(xù)時間以及重載貨車占比對交通流的影響。從模擬結(jié)果可以得出：在重載貨車占比不變的情形下，重載貨車與小汽車的換道頻率均隨著堵塞持續(xù)時間的增加而降低；當車輛密度ρ≤4veh/km或者ρ≥80veh/km時，交通堵塞持續(xù)時間對交通流的影響不大；當車輛密度4veh/km<ρ<80veh/km時，隨著堵塞持續(xù)時間的增加，車流量與速度均降低。在堵塞持續(xù)時間不變的情形下，隨著重載貨車占比增大，車流量與速度均逐漸減小；當車輛密度ρ≤30veh/km時，重載貨車與小汽車均有較高的換道頻率，重載貨車占比對車流影響也較大；隨著車輛密度的逐漸增大，重載貨車與小汽車的換道頻率均隨著重載貨車占比的增大而逐漸減小。

參考文獻：

［1］AGHABAYK K， SARVI M， YOUNG W. Attribute selection for modelling driver′s car-following behaviour in heterogeneous congested traffic conditions［J］. Transportmetrica A： Transport Science， 2014， 10（5）： 457-468.

［2］MESA-ARANGO R， UKKUSURI S V. Modeling the car-truck interaction in a system-optimal dynamic traffic assignment model［J］. Journal of Intelligent Transportation Systems， 2014， 18（4）： 327-338.

［3］PEETA S， ZHOU W， ZHANG P. Modeling and mitigation of car-truck interactions on freeways［J］. Transportation Research Record， 2004， 1899（1）： 117-126.

［4］NGODUY D. Effect of the car-following combinations on the instability of heterogeneous traffic flow［J］. Transportmetrica B： Transport Dynamics， 2015， 3（1）： 44-58.

［5］DAS A K， CHATTARAJ U. Heterogeneous traffic simulation for urban streets using cellular automata［J］. Arabian Journal for Science and Engineering， 2019， 44（10）： 8557-8571.

［6］王永明.基于元胞自動機的道路交通堵塞仿真研究［J］. 系統(tǒng)仿真學報，2010，22（9）：2149-2154.

WANG Y M. Study of traffic congestion′s simulation based on cellular automaton model［J］. Journal of System Simulation， 2010，22（9）：2149-2154.

［7］錢勇生，曾俊偉，杜加偉，等.考慮意外事件對交通流影響的元胞自動機交通流模型［J］.物理學報，2011，60（6）：103-112.

QIAN YS， ZENG J W， DU J W， et al. Cellular automaton traffic flow model considering influence of accidents［J］. Acta Physica Sinica， 2011，60（6）：103-112.

［8］劉霞，胡凱，龔鵬.基于元胞自動機的事故交通流模型仿真［J］.解放軍理工大學學報（自然科學版），2017，18（3）：285-288.

LIU X， HU K， GONG P. Simulation of accident traffic flow model based on cellular automata［J］. Journal of PLA University of Science and Technology（Natural Science Edition）， 2017，18（3）：285-288.

［9］姬浩，王永東，李佩，等.事故車輛影響下的城市三車道道路交通流仿真［J］.中國安全科學學報，2021，31（3）：112-120.

JI H， WANG Y D， LI P， et al. SU Bing. Traffic flow simulation of urban three-lane road considering influence of accident vehicle［J］. Journal of Safety Science and Technology， 2021， 31（3）：112-120.

［10］ 石京，羅靈，王恩達.基于CA模型的交通事故疊加風險量化分析方法［J］.交通信息與安全，2019，37（6）：56-62.

SHI J， LUO L， WANG E D. A quantitative analysis method of traffic crash risk based on cellular automata model［J］. Journal of Transport Information and Safety， 2019， 37（6）：56-62.

［11］ ME LARRAGA， ALVAREZ-ICAZA L. Cellular automaton model for traffic flow based on safe driving policies and human reactions［J］. Physica A： Statistical Mechanics and Its Applications， 2010， 389（23）：5425-5438.

［12］ LI X， XIAO Y， JIA B. Modeling mechanical restriction differences between car and heavy truck in two-lane cellular automata traffic flow model［J］. Physica A： Statistical Mechanics and Its Applications， 2016， 451（C）： 49-62.

［13］ 王凱. 基于元胞自動機的高速公路重載貨車移動瓶頸研究［D］. 武漢：華中科技大學，2017.

WANG K. Research on moving bottleneck of heavy trucks driving on freeway based on cellular automata［D］. Wuhan：Huazhong University of Science and Technology，2017.

［14］ NAGEL K， SCHRECKENBERG M. A cellular automaton model for freeway traffic［J］. Journal De Physique I， 1992， 2（12）：2221-2229.

（責任編輯 李 進）