全國大學(xué)生數(shù)學(xué)建模競賽A題車輛排隊長與事故持續(xù)時間道路實際通行能力路段上流流量間的關(guān)系

1、道路上不斷增加的交通流經(jīng)常導(dǎo)致?lián)頂D。 擁擠產(chǎn)生延誤、 降低流率、 帶來燃油損耗和負(fù)面的環(huán)境影響 。 為了提高道路系統(tǒng)的效率， 國內(nèi)外許多研究者一直致力于車流運行模型的研究。 Daganzo 1提出了一種和流體力學(xué)LWR 模型相一致的元胞傳輸模型 ， 這種模型能用來模擬和預(yù)測交通流的時空演化 ， 包括暫時的現(xiàn)象， 如排隊的形成、 傳播 、 和消散 。 Heydecker 和 Addison 2 通過研究車速和密度的因果關(guān)系分析和模擬了在變化的車速限制下的交通流 。 Jennifer 和 Sallissou 3 提出了一種混合宏觀模型有效地描述了路網(wǎng)的交通流。然而 ， 擁擠也會由交通異常事件引起。

2、 交通異常事件定義為影響道路通行能力的意外事件 4 ， 如交通事故、 車輛拋錨、 落物 、 短期施工等， 從廣義角度看 ， 還應(yīng)包括惡劣天氣與特殊勤務(wù)等。 異常事件往往造成局部車道阻塞或關(guān)閉， 形成交通瓶頸， 引起偶發(fā)性擁擠， 這已經(jīng)逐漸成為高速道路交通擁擠的主要原因 5 ， 越來越多地受到研究者們的重視。例如 M. Baykal-Gursoy 6 等人提出了成批服務(wù)受干擾下的穩(wěn)態(tài)M/M/c 排隊系統(tǒng)模擬了發(fā)生異常事件的道路路段的交通流。 Chung 7 依據(jù)韓國高速公路系統(tǒng)監(jiān)測的準(zhǔn)確記錄的大型交通事故數(shù)據(jù)庫提出了一種事故持續(xù)時間預(yù)測模型。 當(dāng)然 ， 這些研究最終都是為了幫助緩解異常事件引起的

3、交通擁擠。交通異常事件發(fā)生后， 事發(fā)地段通行能力減小 ， 當(dāng)交通需求大于事發(fā)段剩余通行能力時， 車輛排隊 ， 產(chǎn)生延誤， 行程時間增加 8 ， 交通流量發(fā)生變化 。 本文以高速公路基本路段發(fā)生交通事故為例，主要分析了交通事故發(fā)生后不同時間段內(nèi)事故點及其上游下游路段交通流量的變化， 用于以后進(jìn)一步的突發(fā)事件下交通流預(yù)測工作。1 交通事故影響時間分析由于從交通事故發(fā)生到檢測到事故、 接警 、 事故現(xiàn)場勘測、 處理 、 清理事故現(xiàn)場恢復(fù)交通， 以及恢復(fù)交通后車輛排隊不再增加都需要一定的時間。這部分時間主要由三部分構(gòu)成: 第一部分是事故發(fā)生到警察到達(dá)現(xiàn)場的時間T 1。 第二部分是交通事故現(xiàn)場處理時間T

4、2， 由現(xiàn)場勘測、 處理到事故族除、 恢復(fù)交通 。 第三部分是交通事故持續(xù)影響時間T3，這部分時間從恢復(fù)事故現(xiàn)場交通開始， 到事故上游車輛排隊不再增加， 即排隊開始減弱 9 。在 T1內(nèi) ， 事故現(xiàn)場保持原狀， 沒有進(jìn)行處理， 這里分兩種情況考慮: ( 1) 當(dāng)交通事故占用部分車道時，這時事故點的剩余通行能力QsW 0,交通事故越嚴(yán)重， 則相應(yīng)Qs越小。 若事故點上游的交通需求Q vQs， 則車輛以較低的速度通過事故點， 上游不會形成車輛擁擠排隊。若Q >Qs， 則交通流可按事故點的剩余斷面通行能力通過事故點， 超過該通行能力的車流在事故點上游排隊。 ( 2) 當(dāng)交通事故十分嚴(yán)重時 ，

5、事故點的剩余通行能力Qs= 0， 造成事發(fā)路段斷流 ， 事故點上游車輛排隊， 發(fā)生交通擁擠堵塞， 進(jìn)而排隊一直向上游延伸。在 T2內(nèi) ， 確認(rèn)交通事故發(fā)生后， 相關(guān)部門到現(xiàn)場處理異常事件， 在此過程中， 事故點交通可能會受到進(jìn)一步影響， 事故斷面通行能力也隨之發(fā)生變化 5 ， 一般會變小， 甚至變?yōu)?( 全封閉處理) ， 視事件處理具體情況而定， 事發(fā)點上游交通處于嚴(yán)重?fù)頂D狀態(tài) ， 車輛排隊增加。由于在交通事故接警時間T 1和處理時間T2階段事故點上游交通車輛產(chǎn)生排隊， 若沒有車輛排隊 ， 則 T3= 0。 若有車輛排隊， 則當(dāng)事故處理完畢、 道路恢復(fù)交通時， 排隊車輛開始消散。 交通事故持續(xù)

6、影響時間T3是事故處理完畢、 道路恢復(fù)交通至車輛排隊不再增加這段時間， 即交通流消散波從車輛排隊隊列的頭部傳到尾部這段時間 9 。2 事故路段車輛排隊長度分析如圖 1 所示 ， 設(shè)某高速公路基本路段長度為L( m) ， 單方向車道數(shù)為n， 單方向車道寬度為D( m) ，在道路上t = 0 時刻發(fā)生了一起交通事故， 事故車輛占用道路寬度為b( m) ， 長度為 a( m) ， 事故點上游路段長度為L' 。 假設(shè)車輛的到達(dá)率為Q， 在同級服務(wù)水平上事故發(fā)生斷面通行能力為Qs， 道路在正常條件下的單方向的通行能力為Qi。圖 1 發(fā)生交通事故的高速公路基本路段本文暫只考慮如圖1所示的基本路段內(nèi)

7、的車輛排隊長度，這里不同于以往文獻(xiàn)的“排隊長度”，以往文獻(xiàn)中的“排隊長度”沒有區(qū)分不同的“阻塞行車 道寬度”。這里的“阻塞行車道寬度”不只是事故車 輛實際占用寬度，還包括虛擬占用寬度，比如事故 發(fā)生位置橫跨在兩車道之間，導(dǎo)致事故點只能通行 一個車道寬度的車流，那么此時“阻塞行車道寬度” 為兩個車道的寬度。設(shè)Q >Qs,m( T i)為時間內(nèi)事故點阻塞行車道寬度(本文把單個車道寬度和車輛寬度看作同寬)，Lm(t)為t時刻事故點上游路段L'內(nèi)車流以阻塞行車道寬度m的排隊長度，且w( T i)=u f 1 kii+ k i ()2k j為Ti時間內(nèi)新產(chǎn)生的交通波10的速度，其中Uf為該

8、事故路段的自由流速度，即該路段的設(shè)計車速，可以通過城市地理信息平臺GIS得到道路基本數(shù)據(jù)。kii、ki2分別為Ti內(nèi)事故點上游、事故點瓶頸段的交通密度，可以由交通 檢測系統(tǒng)監(jiān)測得到。kj為該路段的交通堵塞密度，由道路的基本數(shù)據(jù)可以計算得到9。2. 1 0 <t<Ti(1)若 L'w( T 1),則t v L'w( T 1)時，L m( T) ( t) = w( T 1)t v L'。t>L'w( T 1)時，L m( T) ( t) = L'。L m( T) m( T)( t) = 0 。3 / 12(2)若 TiWL'w( T

9、 1)，則L m( T) ( t) = w( T 1) t， Lm( T) m( T) ( t) = 0 ，3. 2 T i <t< Ti+ T2在 T2時間內(nèi) ， 事發(fā)點斷面通行能力一般會變化 ， 設(shè)變?yōu)?Q'S， 則 m( Ti) 也會相應(yīng)發(fā)生變化。 這里 ，還需要考慮一個時間， 就是交通波w ( T2) 趕上w( T 1) 的時間 ( 設(shè)為 T' 2) ， 趕上之后車流以w( T 2) 的速度 、 m( T i) 的寬度繼續(xù)排隊。 本文由于只考慮L'內(nèi)的排隊長度， 所以考慮在T2時間內(nèi)且在L' 段內(nèi)交通波 w( T 2) 是否趕上w( T 1)

10、 ， 即 T' 2若同時滿足以下兩個條件才需被考慮:L m( T) ( T1)< L' , Lm( T) ( Ti) + w( T 1)T' 2V L'。T'2<T2o28 期陳誠 ， 等 : 交通事故影響下事發(fā)路段交通流量變化分析69052. 2. 1 Lm( T) ( T1)< L'若 Lm( T) ( T1) + w( T 1)T' 2 W L'且T'2 <T2，則:當(dāng) T2<( L' Lm( T) ( Ti) w( T 1) T' 2) /w( T 2)+ T'

11、 2時，I m( T) ( t) = L m( T) ( T1)+ | w( T i)1( t -Ti),當(dāng) t TiWT'2時II m( T) ( t) = L m( T) ( T1)+ | w( T 1)| T'2 + | w( T 2)I(t T i - T'2),當(dāng)t Ti >T'2IIIII時oL m( T) m( T) (t) = w( T 2)(t -Ti),當(dāng)t TiWL' / w( T 2)時Lm( T) m( T)( t) = L',當(dāng) t T > L' / W( T 2) 時o當(dāng) T2K L' L

12、m( T) ( Ti)w( T i)T'2)/ w( T 2)+T'2時，L m( T) ( t) = L m( T) ( T1)+ | w( T 1)I ( t -T1),當(dāng) t TiWT'2時L m( T) ( t) = L m( T) ( T1)+ | w( T 1)| T'2+ | w( T 2)1(t T i T' 2),當(dāng)(L' Lm( T)( T1)+| w( T1)I T2)/I w( T2)I + T' 2 >t T1 >t'2時Lm( T) ( t) = L',當(dāng)(L' Lm( T)

13、 ( Tl)+ | w( T i)I T2)/I w( T2)I + T'2Vt -TKt2min min 時oL m( T) m( T) (t) = w( T 2)(t -Ti),當(dāng)t TiWL' / w( T 2)時Lm( T) m( T)( t) = L',當(dāng) t T > L' / W( T 2)時 。若 Lm( T) ( Ti)+ w( T i)T' 2 > L',則L m( T) ( t) =L m( T) ( T1)+ w( T i)(t -Ti5 / 12),當(dāng) t TiWL' Lm( T)( T i )/ w(

14、 T i)時Lm( T) ( t) =L',當(dāng) t T1 >L' Lm( T) ( Ti () J/ w( T 1)時 。L m( T) m( T) (t) = w( T 2)(t -Ti),當(dāng)t TiWL' / w( T 2)時Lm( T) m( T)( t) = L',當(dāng) t Ti > L' / W( T 2)時 。若T' 2 >T2,則L m( T) ( t) = L m( T) ( T1)+ w( T i)(t -Ti),當(dāng) t TiWL' L m( t) ( T i )/ w( T i)時Lm( T) ( t)

15、 =L',當(dāng)t Ti >L' -Lm( T) ( Ti () y w( T i)時 。L m( T) m( T) (t) = w( T 2)(t -Ti),當(dāng)t TiWL' / w( T 2)時Lm( T) m( T)( t) = L',當(dāng) t Ti > L' / W( T 22. 2. 2 L m( T) ( Ti)=L'Lm( T) ( t) = L',L m( T) m( T) (t) = w( T 2)(t Ti),當(dāng)t TiWL' / w( T 2)時L m( T) m( T)( t) = L',當(dāng)

16、t - Ti > L' / W( T 2)時 。2. 2. 3 T i+ T2 <t<Ti+ T2 + T3這里同樣要考慮交通波 w( T 3)趕上w( T 2)、w( Ti)的時間 T'3、T'3。w( T2)追趕w( Ti),趕上之后以 w( T2)的速度向上延伸排隊，w( T3)追趕前面兩者，先趕上w( T 2)后趕上w( T 1)(假如w( T2)還沒趕± w( T 1),趕上之后排隊不再增加，考慮在T3內(nèi)L段內(nèi)w( T 3)是否趕上w( T 2)、w( T i),而w( T3)是消散波何，t時刻排隊長度為w( Ti)、w( T 2)

17、到t時刻為止產(chǎn)生的排隊長度減去消散波w( T 3)向上游傳播延伸的長度L( t),相關(guān)分析及計算式類似上述，不 再贅述。3不同時間段內(nèi)不同路段的交通流量變化分析3. 1 QW Qs如果上游交通量需求低于剩余可通行車道的正常通行能力，那么，盡管事發(fā)路段存在通行能力瓶頸，但不會導(dǎo)致交通擁擠，當(dāng)上游的流量到達(dá) 瓶頸處，密度增大，車速降低，車輛以較低的速度通 過瓶頸點，不會形成排隊 ，該事故路段的交通流 量(veh /h)受影響很小，近似為Q ( veh /h)。在T2 時間內(nèi)，由于事故點斷面通行能力可能進(jìn)一步變?yōu)镼's需要重新判斷Q與Q's的大小：若QWQ's,則交通流影響不

18、大，該基本路段流量(veh /h)仍與Q相當(dāng)。若 Q >Q's,上游車輛排隊，過程類似于Q >Qs時Ti內(nèi)車流變化，見下面分析，不再贅述。3. 2 Q >Qs此時事故點上游路段車流排隊，如上文所述，6906科學(xué)技術(shù)與工程11卷由于在事故發(fā)生后事故影響時間的不同時間段內(nèi)車輛排隊的不同變化，上游路段不同位置不同時刻 的斷面短時交通流量會有所不同：設(shè)Lp為上游斷面P距離事故點的長度，Qp(t)為斷面P在t( s)時刻的短時交通流量，則Qp(t)=Q,當(dāng) L P > L m( T) ( t) , L P > L m( T) m( T) (t)時Qs,當(dāng) L P w

19、 L m( T) ( t) , L P > L m( T) m( T) ( t)時 Q's,當(dāng) LPWLm( T) ( t) , L P < L m( T) m( T)(t)時(0 <t<T1+ T2)。Qp(t) =Qi(Lp< L( t) , T1+ T2 <t< T1+ T2+ T 3)(即疏散波到達(dá)的斷面排隊消散，車流以最大的通行量前進(jìn))。Ti+ T2 + T3后整個路段恢復(fù)到事故發(fā)生之前的交通流量。事故點和事故點下游路段各斷面短時交通流量在事故發(fā)生后的不同時間段內(nèi)變化基本一致：0 Vt<Ti 內(nèi) ， 事故點和事故點下游的各斷面短

20、時交通流量相當(dāng)于事故點的斷面通行能力QS。Ti <t<Ti+ T2內(nèi) ， 事故點和事故點下游的各斷面短時交通流量相當(dāng)于事故點的斷面通行能力Q'S。Ti+ T2 vtw+T2 + T3內(nèi) ， 事故點和事故點下游的各斷面短時交通流量相當(dāng)于道路通行能力Qi。 T1+ T2 + T3后道路暢通 ， 恢復(fù)到事故發(fā)生前的交通流量。4 仿真驗證本文采用德國PTV 公司的交通微觀仿真軟件VISSIM 11 對上述分析進(jìn)行仿真驗證。 仿真中建立如圖 1 所示的基本路段， 路段長度L = 500 m ， 由單方向三車道組成。 事故車輛?？康奈恢迷O(shè)在右車道 L' = 346 m 處 ，

21、在上游路段均勻設(shè)置四個數(shù)據(jù)采集點 ， 事故點斷面設(shè)置一個采集點， 下游路段均勻設(shè)置三個采集點。 事故發(fā)生的時間設(shè)置在320 s 左右 ， 即上述的t = 0 時刻 ， T1= 280 s， T2 = 330 s。 T2內(nèi)是事故現(xiàn)場處理， 仿真中采用信號燈控制的方式 ， 信號燈設(shè)置在靠近事故點前方， 考慮到事故處理時占用更多的車道， 這里設(shè)置為控制右車道和中間車道， 即阻斷兩車道進(jìn)行事故處理。 在 600 s 時紅燈亮， 同時阻斷右車道和中間車道進(jìn)行事故處理 ， 930 s 時綠燈亮， 事故處理完畢， 兩閉塞車道暢通 ， 道路恢復(fù)， 仿真時間設(shè)為3 600 s。 仿真中 ， 分別以 6 000

22、veh /h 和 800 veh /h 的流量進(jìn)行道路車輛到達(dá)率輸入（輸入車輛類型比例為重型貨車：小汽車=0. 02 : 0. 98）:圖 2 事故點上游不同斷面短時交通流量變化圖 3 事故點斷面短時交通流量變化4 . 1 車輛輸入為6 000 veh /h 時各數(shù)據(jù)采集點每隔30 s 采集一次斷面車輛通過量數(shù)據(jù)， 圖 2 圖 6 為不同斷面每隔30 s 交通量變化圖 ， 橫坐標(biāo)表示的是180 s 1 200 s 共 34 個 305 的時間間隔， 縱坐標(biāo)表示在相應(yīng)的30 s 內(nèi)通過該11 / 12斷面的交通量(veh - ( 30 s) i) , ser01ser08表示28 期陳誠 ， 等

23、 : 交通事故影響下事發(fā)路段交通流量變化分析6907圖1 所示的8個數(shù)據(jù)采集點斷面 ，觀察各圖可以看出這樣幾點:事故點斷面和下游各斷面流量變化較一致 ， 曲線較重合， 在事故發(fā)生后， 下游斷面流量和事故點相同。事故點上游各斷面流量變化趨勢大體一致 ， 越接近事故點的斷面流量變化趨勢越提前。各斷面流量變化都經(jīng)歷這樣一個過程：Q- Qsf Q'sfQe Q。圖 4 事故點下游不同斷面短時交通流量變化圖 5 事故點及其下游不同斷面短時交通流量變化4 . 2 車輛輸入為800 veh /h ( 非隨機(jī) ) 時各數(shù)據(jù)采集點每隔3 600 s 采集一次斷面車輛通過量數(shù)據(jù)， VISSIM 中道路輸入

24、流量是以veh /h 來計的 ， 輸入車輛類型為大貨車、 小汽車 ， 而且車輛進(jìn)入路段服從泊松分布， 所以在車輛輸入量較小的情況下 ， 數(shù)據(jù)波動較大， 采取短時間( 30 s) 內(nèi)采集各斷面數(shù)據(jù)是不準(zhǔn)確的。 這里采取3 600 s( 一小時 ) 的間隔從 200 s 時開始采集小時交通量來考察。 仿真得到從上游到下游各斷面到3 800 s 時流量分別為796 veh /h 、 795 veh /h 、 794 veh /h 、 794 veh /h 、 794 veh /h 、 793 veh /h 、 793 veh /h 、 793 veh /h ， 說明車輛到達(dá)率較小時， 交通事故對該路

25、段的交通流量影響不大 ， 可認(rèn)為交通流量不變。5 結(jié)論本文主要考慮了不同時間段內(nèi)不同阻塞行車道寬度的車輛排隊長度， 對交通事故影響下不同時間段內(nèi)交通流量變化進(jìn)行了分析， 為以后的異常事件下的交通流預(yù)測奠定基礎(chǔ)， 下一步的工作是考慮非參數(shù)回歸法在交通事故等異常事件影響下的交通流量預(yù)測算法。圖 6 事故點及其上游不同斷面短時交通流量變化參考文獻(xiàn)1 Daganzo C F. The cell transmission model: A dynamic representationof highway traffic consistent with the hydrodynamic theory. T

26、ransportation Research Part B ， 1994。28B( 4) : 269 2876908 科學(xué)技術(shù)與工程11 卷2 Heydecker B G ， Addison J D. Analysis and modeling of traffic flow under variable speed limits. Transportation Research Part C ， 2011。19: 206 2173 Mccrea J ， Moutari S. A hybrid macroscopic-based model for traffic flow in road n

27、etworks. European Journal of Operational Research ， 2010。 207: 676 6844 Skabardonis A ， Petty K ， Varaiya P ， et al. Evaluation of the freewayservice patrol in Los Angeles UCB-ITS-PRR-98-31 ， CaliforniaPATH Research Report ， Institute of Transportation Studies ， Universityof California ， Berkeley ，

28、19985 張海軍 ， 張玨 ， 楊曉光 . 異常事件下高速道路交通狀態(tài)的分析與仿真 . 交通運輸工程學(xué)報， 2008。8( 2) : 116 1216 Baykal-gursoy M ， Xiao W， Ozbay K. Modeling traffic flow interrupted by incidents. European Journal of Operational Research ， 2009。195: 127 1387 Chung Y. Development of an accident duration prediction model on the Korean Freeway Systems. Accident Analysis and Prevention ，2010。