王 明 常云濤 邱紅桐 董可然

(公安部交通管理科學(xué)研究所1) 無錫 214151) (同濟大學(xué)道路與交通工程教育部重點實驗室2) 上海 201804)

0 引 言

行駛工況也稱運轉(zhuǎn)循環(huán)，是某一類型車輛(如乘用車、公交車、重型車輛等)在特定環(huán)境(如城區(qū)，郊區(qū)等)下車輛行駛的速度-時間變化規(guī)律[1].車輛在行駛過程中，影響車輛行駛工況的因素復(fù)雜，既包括駕駛員和車輛本身的內(nèi)在因素，也包括交通管理控制、交通流密度、道路物理幾何等外在因素.汽車行駛工況可以直接反映特定地區(qū)車輛在行駛過程中速度的變化情況，反映車輛行駛過程的平順性，也能夠評估車輛在行駛過程中的能源消耗量和污染物排放量，并助力新型車開發(fā)等工作[2-3].

國內(nèi)外學(xué)者對車輛行駛工況做了大量的研究.其中，對于車輛行駛工況構(gòu)建方法的研究，姜平等[4]采用離散小波變換的多分辨信號分解算法，以短行程為單元對汽車行駛工況進(jìn)行構(gòu)建.李寧[5]通過聚類分析和主成分分析對車輛運動學(xué)片段進(jìn)行劃分，并按照固定時間長度進(jìn)行車輛行駛工況構(gòu)建.李友文等[6]通過馬爾可夫過程來構(gòu)建不同狀態(tài)之間轉(zhuǎn)移矩陣，通過狀態(tài)轉(zhuǎn)移矩陣來構(gòu)建車輛行駛工況.另外，對于車輛行駛工況與道路環(huán)境關(guān)系的研究，毛林鋒等[7-9]研究城市道路交叉口間距對車輛行駛工況的影響，交叉口間距越小，車輛的行駛工況越不穩(wěn)定，道路上的交通事故越多，出行時間、能耗與排放越大.卓亞娟[10]研究城市道路交叉口信號配時方案的不同對車輛的行駛工況影響，針對不同的周期、綠信比和相位差，車輛行駛過程中的勻速、加速、減速、怠速狀態(tài)都會不同，由此帶來的車輛延誤、能耗和排放也不同.Greenshields[11]提出了流量-速度-密度之間的函數(shù)關(guān)系，體現(xiàn)了交通量對于車輛行駛工況影響.

然而，當(dāng)前普遍使用的車輛行駛工況是一種基于時間參數(shù)的車輛速度隨機變化序列，見圖 1a).這種行駛工況無法反映車輛的空間位置信息，難以反映行駛工況與道路不同位置的交通服務(wù)水平、交叉口間距、信號控制參數(shù)的關(guān)系.為了解決這一問題，本文將行駛工況重新定義為車速沿車輛行駛方向的空間變化過程，見圖 1b).這些交通行駛環(huán)境因素是與車輛所處的道路空間構(gòu)成直接關(guān)聯(lián)，車輛的行駛工況不再是一種隨時間變化的過程，而是被定義為沿道路行駛方向的空間變化過程.

圖1 行駛工況

行駛環(huán)境是指車輛所在道路的交通服務(wù)水平、距離交叉口的位置、信號控制參數(shù)等，這些環(huán)境參數(shù)與車輛在道路上所處的位置有關(guān).本文采集行駛工況數(shù)據(jù)的道路交叉口的信號控制參數(shù)無變化，并且選擇暢通狀態(tài)下行駛工況數(shù)據(jù)進(jìn)行分析，故而本文只探索車輛行駛工況隨著空間位置改變的變化規(guī)律，分析沿道路行駛方向車輛行駛狀態(tài)的變化規(guī)律以及車輛能耗和排放變化規(guī)律.

1 車輛行駛狀態(tài)劃分的基本原理

在聚類分析中，樣本之間的相似性通常采用樣本之間的距離來表示.數(shù)據(jù)之間的距離通常用明氏(Minkowski)距離表示，當(dāng)q=1時，稱為絕對距離；當(dāng)q=2時，稱為歐式距離；當(dāng)q=∞時，稱為切比雪夫距離.

(1)

車輛在行駛過程中，車輛行駛速度會隨著車輛位置的變化而變化.為了反映車輛在行駛過程中的車輛行駛速度的變化規(guī)律，按照車輛行駛過程中的速度和加速度值差異，可通過聚類分析將沿道路行駛方向車輛行駛狀態(tài)進(jìn)行聚類.

2 空間序列數(shù)據(jù)的獲取方法

2.1 工況數(shù)據(jù)采集范圍

本文車輛行駛工況數(shù)據(jù)來源于美國佛羅里達(dá)州NW 13th street道路(雙向4車道、雙黃線隔離)，實驗車(小汽車、功率147 kW、排量2.0 L手動)共有效采集了單程97組車輛行駛工況數(shù)據(jù)，根據(jù)各組行駛工況標(biāo)準(zhǔn)化后的平均速度和速度變異系數(shù)差異，經(jīng)狀態(tài)聚類，獲得暢通工況67組，協(xié)同流工況25組，擁堵工況5組，暢通工況可作為理想的分析對象，因此本文只是對暢通狀態(tài)下的車輛行駛工況數(shù)據(jù)進(jìn)行分析.車輛行駛工況數(shù)據(jù)采集范圍從距SW 10th street交叉口停車線后100 m到SW 16th street交叉口停車線后140 m，整個行駛區(qū)域共計長650 m，見圖 2.

圖2 工況數(shù)據(jù)采集范圍

2.2 基于空間位置的工況數(shù)據(jù)獲取

實驗車輛配備有OBD設(shè)備和GPS設(shè)備，在車輛行駛過程中，OBD設(shè)備主要采集車輛行駛工況數(shù)據(jù)，包括速度、加速度、能耗和排放等參數(shù)，GPS設(shè)備主要記錄車輛的位置和行駛路徑信息，包括時間、經(jīng)度、緯度和高程等參數(shù).設(shè)置OBD設(shè)備和GPS設(shè)備同時按秒進(jìn)行數(shù)據(jù)采集，統(tǒng)一使用GPS時間，因此兩個設(shè)備采集的數(shù)據(jù)包含相同的時間參數(shù)，通過時間參數(shù)來實現(xiàn)車輛行駛工況數(shù)據(jù)與具體空間位置的結(jié)合.

2.3 空間序列的行駛工況數(shù)據(jù)獲取

通過時間參數(shù)相關(guān)聯(lián)的OBD設(shè)備和GPS設(shè)備采集的數(shù)據(jù)是按照每秒統(tǒng)計的，雖然具有了空間位置信息，但是依然屬于時間序列數(shù)據(jù).將每秒輸出的車輛速度、加速度進(jìn)行變換為沿道路行駛方向每相等間隔的車輛行駛速度和加速度數(shù)據(jù)，加速度通過相鄰速度差獲得.

車輛速度的變換公式

(2)

圖3為每10 m間隔的空間序列工況數(shù)據(jù)，可以發(fā)現(xiàn)，路段范圍車輛速度較為平穩(wěn)，當(dāng)車輛到達(dá)交叉口范圍時通常會出現(xiàn)明顯的減速停車過程，接著啟動加速駛離交叉口，基于空間位置的車輛行駛工況能夠清楚反映車輛行駛速度隨著道路空間位置轉(zhuǎn)移的變化規(guī)律.

圖3 67條行駛工況每10 m間隔的速度值

3 沿道路行駛方向車輛行駛工況及能耗和排放分析

3.1 車輛行駛狀態(tài)的劃分

當(dāng)機動車通過信號控制交叉口時，由于受到信號控制的影響，車輛會頻繁出現(xiàn)減速、怠速、加速等狀態(tài)，相對于路段勻速行駛過程，車輛能耗和尾氣排放因此會有較大增加.為了研究沿道路行駛方向車輛行駛狀態(tài)的變化規(guī)律以及能耗和排放的變化規(guī)律，將沿道路行駛方向每10 m間隔的車輛速度和加速度值進(jìn)行標(biāo)準(zhǔn)化，經(jīng)聚類分析，共獲得7類車輛行駛狀態(tài)，也就是7個模式事件，其中：狀態(tài)1代表車輛處于高速勻速行駛狀態(tài)、狀態(tài)2代表車輛處于高速減速行駛狀態(tài)、狀態(tài)3代表車輛處于中速加速行駛狀態(tài)、狀態(tài)4代表車輛處于中速勻速行駛狀態(tài)、狀態(tài)5代表車輛處于中速減速行駛狀態(tài)、狀態(tài)6代表車輛處于低速加速行駛狀態(tài)，狀態(tài)7代表車輛處于低速勻速行駛狀態(tài)(包含怠速狀態(tài)).各狀態(tài)聚類中心見表1，聚類效果見圖4.

表1 7類行駛工況狀態(tài)的聚類中心

圖4 沿道路行駛方向車輛行駛狀態(tài)劃分

根據(jù)上述行駛狀態(tài)的聚類結(jié)果，選取5條車輛行駛工況繪制沿道路行駛方向車輛行駛狀態(tài)變化圖，見圖5.

圖5 沿道路行駛方向車輛行駛狀態(tài)變化圖

3.2 各車輛行駛狀態(tài)下車輛能耗和排放率

為了對沿道路行駛方向每10 m間隔車輛能耗和排放進(jìn)行分析，并得到車輛行駛過程中各類行駛狀態(tài)對應(yīng)的能耗和排放數(shù)據(jù)，將按每秒間隔采集的車輛能耗和排放數(shù)據(jù)變換為沿道路行駛方向每10 m間隔車輛的能耗和排放數(shù)據(jù).車輛能耗的變換公式為

(3)

式中：C為10 m細(xì)分段的能耗值，L；ci為10 m細(xì)分段內(nèi)第i個非怠速的能耗值，L/h；cj為10 m細(xì)分段內(nèi)第j個怠速的能耗值，L/h；vi為10 m細(xì)分段內(nèi)第i個非怠速的速度值，km/h；k1為10 m細(xì)分段內(nèi)速度非怠速的個數(shù)；k2為10 m細(xì)分段內(nèi)怠速的個數(shù).

車輛排放的變換公式為

(4)

根據(jù)上述能耗和排放計算方法，獲得車輛行駛過程中每10 m間隔車輛的能耗和排放數(shù)據(jù)，圖 6為7類狀態(tài)的能耗和排放分布，可以發(fā)現(xiàn)，每類狀態(tài)車輛能耗和排放分布區(qū)間差異明顯，但都比較集中，每類狀態(tài)能耗和排放的均值可用于代表該類狀態(tài)的能耗和排放.

圖6 狀態(tài)的能耗和排放分布圖

根據(jù)7類行駛狀態(tài)下車輛行駛能耗和排放值以及分布頻率，計算7類車輛行駛狀態(tài)下的車輛平均能耗和排放值，見表 2.

表2 各車輛行駛狀態(tài)平均能耗和排放

圖7為7類車輛行駛狀態(tài)下車輛平均能耗和排放的關(guān)系圖，其中，無填充柱型代表每個狀態(tài)的平均能耗，斜線填充柱型代表每個狀態(tài)的平均排放.可以發(fā)現(xiàn)：①車輛能耗和排放呈現(xiàn)正相關(guān)性；②狀態(tài)1和狀態(tài)4的車輛能耗和排放基本相同，表明車輛在保持中高速勻速行駛狀態(tài)行駛時，車輛能耗和排放處于穩(wěn)定狀態(tài)；③狀態(tài)2和狀態(tài)5車輛能耗和排放都很低，大小基本相同，表明減速狀態(tài)車輛能耗和排放與減速這一過程有關(guān)，而與減速過程中的速度大小關(guān)系很??；④狀態(tài)3和狀態(tài)6車輛能耗和排放都明顯較大，表明加速過程會增加車輛能耗和排放，并且起步加速階段能耗和排放更大；⑤狀態(tài)7車輛能耗和排放最大，因為這個狀態(tài)包含了車輛怠速狀態(tài)，車輛在怠速狀態(tài)車輛發(fā)動機空轉(zhuǎn)，車輛在交叉口停車等待時間越長，這一過程能耗和排放越大.

圖7 各車輛行駛狀態(tài)下車輛平均能耗和排放

3.3 車輛行駛狀態(tài)占比的空間分布規(guī)律

實驗數(shù)據(jù)共包含了67組行駛工況數(shù)據(jù)，按照10 m間隔進(jìn)行劃分，單條650 m車輛行駛工況包含了65個車輛行駛狀態(tài)，統(tǒng)計沿道路行駛方向每10 m間隔每類狀態(tài)所占比例，其結(jié)果見圖8，其中，第510 m處為交叉口停車線位置.

1) 狀態(tài)1為高速勻速行駛，前部路段區(qū)域，狀態(tài)1所占比例比較穩(wěn)定，隨著離交叉口停車線越來越近，狀態(tài)1所占比例逐漸減少，隨著車輛遠(yuǎn)離交叉口，狀態(tài)1所占比例逐漸增加.這是因為車輛在通過交叉口時受到信號控制影響，車輛逐漸減速由狀態(tài)1轉(zhuǎn)為狀態(tài)2，通過交叉口后車輛加速又逐漸變?yōu)闋顟B(tài)1.

2) 狀態(tài)2為高速減速行駛，隨著車輛靠近交叉口，狀態(tài)2所占的比例先逐漸增加然后逐漸減少.前期逐漸增加是因為高速行駛的車輛遇到交叉口開始減速從而由狀態(tài)1轉(zhuǎn)化為狀態(tài)2，后期逐漸減少是因為車輛速度降到一定程度后保持中速行駛或者繼續(xù)減速行駛或者開始加速行駛，即狀態(tài)2又轉(zhuǎn)為狀態(tài)4或狀態(tài)3或狀態(tài)5.

3) 狀態(tài)3為中速加速行駛，在車輛通過交叉口后開始逐漸增加，后期逐漸減少.這是因為狀態(tài)3為中速加速過程，車輛通過交叉口后逐漸加速，當(dāng)車輛達(dá)到中速后繼續(xù)加速就變?yōu)橹兴偌铀龠M(jìn)入狀態(tài)3.隨著車輛速度進(jìn)一步加大，又逐漸轉(zhuǎn)為高速行駛進(jìn)入狀態(tài)1.

4) 狀態(tài)4為中速勻速行駛，在距離交叉口停車線前約200 m位置開始出現(xiàn)，往后小幅度增加，說明在交叉口范圍保持中速勻速行駛的車輛不多.這一狀態(tài)是由車輛高速行駛減速到中速并勻速行駛，即狀態(tài)2轉(zhuǎn)化為狀態(tài)4的過程.

5) 狀態(tài)5為中速減速行駛，在距離交叉口停車線前約200 m位置開始出現(xiàn)，在70 m處達(dá)到最大值，隨后逐漸下降直到停車線位置.這一狀態(tài)是由狀態(tài)2(高速減速狀態(tài))或者狀態(tài)4(中速勻速狀態(tài))轉(zhuǎn)變而來，由于交叉口停車線后車輛開始加速駛離交叉口，所以停車線后就不存在狀態(tài)5.

6) 狀態(tài)6為低速加速行駛，在距離交叉口停車線前約100 m逐漸增加，并且在停車線附近達(dá)到最大值，從停車線往后約100 m狀態(tài)6逐漸減少直到消失.這一狀態(tài)由低速狀態(tài)轉(zhuǎn)變而來，即由狀態(tài)7轉(zhuǎn)為狀態(tài)6.由于車輛一般是在交叉口上游集聚，等待綠燈啟亮?xí)r開始加速駛離交叉口，因此這一狀態(tài)占比會在停車線附近達(dá)到最大值.

7) 狀態(tài)7為低速勻速行駛，主要在交叉口停車線前100 m范圍內(nèi)存在，包含了車輛在交叉口停車怠速過程以及低速通過交叉口過程.

圖8 沿道路行駛方向各狀態(tài)占比變化規(guī)律

3.4 沿道路行駛方向車輛能耗和排放

根據(jù)以上求得7類車輛行駛狀態(tài)下車輛平均能耗和排放值，結(jié)合沿道路行駛方向車輛行駛狀態(tài)的變化，獲取沿道路行駛方向每10 m間隔內(nèi)車輛平均能耗和排放的變化值.各細(xì)分段平均能耗計算公式為

(5)

式中：Ci為第i個細(xì)分段車輛能耗值(i=1,2,…,65)，L/100 km；pk為第i個細(xì)分段內(nèi)狀態(tài)k所占的比例(k=1,2,…,7)；ck為車輛行駛狀態(tài)k的平均能耗值，L/100 km.

各細(xì)分段平均能耗計算公式為

(6)

式中：Di為第i個細(xì)分段車輛排放值(i=1,2，…，65)，g/km；pk為第i個細(xì)分段內(nèi)狀態(tài)k所占的比例(k=1,2，…，7)；dk為車輛行駛狀態(tài)k的平均排放值，g/km.

根據(jù)上述計算方法，生成沿道路行駛方向每10 m間隔內(nèi)車輛平均能耗和排放(CO2)變化情況，見圖9.

注：虛線代表停車線位置圖9 沿道路行駛方向車輛平均能耗和排放

由圖9可知：①路段范圍車輛行駛過程中車輛能耗和排放較為穩(wěn)定，能耗值在5～8 L/100 km，排放量在130～200 g/km；②車輛能耗和排放在交叉口停車線前150 m左右到停車線后150 m左右先逐漸上升后逐漸下降，并且在交叉口停車線附近達(dá)到最大值.停車線前能耗和排放逐漸增大是由于車輛在停車線前低速行駛和怠速狀態(tài)產(chǎn)生的過高能耗和排放；停車線后車輛能耗和排放逐漸下降說明前期車輛加速過程能耗和排放較大，隨著車速達(dá)到一定程度后車輛處于勻速行駛狀態(tài)，車輛能耗和排放相對較小，并且趨于穩(wěn)定.

4 結(jié) 論

1) 沿道路行駛方向的車輛行駛狀態(tài)能夠進(jìn)行有效聚類.包括：高速勻速行駛狀態(tài)、高速減速行駛狀態(tài)、中速加速行駛狀態(tài)、中速勻速行駛狀態(tài)、中速減速行駛狀態(tài)、低速加速行駛狀態(tài)和低速勻速行駛狀態(tài)(包含怠速狀態(tài))，并且分析發(fā)現(xiàn)，每類狀態(tài)內(nèi)的能耗和排放分布都比較集中，每類狀態(tài)的平均值能夠有效代表該狀態(tài)的能耗和排放.

2) 車輛行駛狀態(tài)受道路行駛環(huán)境影響變化顯著.根據(jù)沿道路行駛方向車輛行駛狀態(tài)變化規(guī)律，路段范圍車輛主要包含狀態(tài)1和狀態(tài)2，并且狀態(tài)比較穩(wěn)定；交叉口范圍車輛行駛狀態(tài)從狀態(tài)1到狀態(tài)7都有分布，并且狀態(tài)之間轉(zhuǎn)換頻繁，這與實際車輛行駛過程相吻合.

3) 受交叉口信號控制影響，沿道路行駛方向車輛能耗和排放變化規(guī)律明顯.①車輛能耗和排放基本成正相關(guān)關(guān)系；②狀態(tài)1和狀態(tài)4車輛能耗和排放基本相同；③狀態(tài)2和狀態(tài)5車輛能耗和排放都很低，大小基本相同；④狀態(tài)3和狀態(tài)6車輛能耗和排放明顯較大；⑤狀態(tài)7車輛能耗和排放最大.⑥路段范圍車輛能耗和排放較為穩(wěn)定，能耗值在5～8 L/100 km，排放量在130～200 g/km；⑦車輛能耗和排放在交叉口停車線前150 m左右到停車線后150 m左右先逐漸增加后逐漸減少，并且在交叉口停車線附近達(dá)到最大值.

下一步將在本文研究基礎(chǔ)上，深入研究車輛行駛工況與道路行駛環(huán)境的內(nèi)在聯(lián)系，探索從比較容易大規(guī)模獲取的幾何、交通控制參數(shù)出發(fā)，生成行駛工況數(shù)據(jù)，這些工況數(shù)據(jù)可以作為車輛能耗和排放分析研究的輸入，應(yīng)用于能耗和排放模型之中，對車輛能耗和排放進(jìn)行分析和計算.

參考文獻(xiàn)

[1] 余志生．汽車?yán)碚揫M]．第五版，北京：機械工業(yè)出版社，2011.

[2] 馬冬，丁焰，劉志華，等.輕型汽車實際行駛工況的排放研究[J].安全與環(huán)境學(xué)，2008,8(5):66-68.

[3] 鄭殿宇,吳曉剛,陳漢等.哈爾濱城區(qū)乘用車行駛工況的構(gòu)建[J].公路交通科技,2017,34(4):101-107.

[4] 姜平,石琴,陳無畏，等.基于小波分析的城市道路行駛工況構(gòu)建的研究[J].汽車工程,2011,33(1):70-73.

[5] 李寧.城市道路車輛行駛工況的構(gòu)建與研究[D].保定：河北農(nóng)業(yè)大學(xué),2013.

[6] 李友文,石琴,姜平，等.基于馬爾科夫過程的行駛工況構(gòu)建中數(shù)據(jù)處理與分析[J].合肥工業(yè)大學(xué)學(xué)報(自然科學(xué)版),2010,33(4):491-494.

[7] 毛林鋒.公路平面交叉口安全間距研究[D].南京：東南大學(xué),2007.

[8] 張?zhí)m芳，方守恩.公路平面交叉口合理間距研究[J].重慶交通大學(xué)學(xué)報(自然科學(xué)版).2008(2):276-278.

[9] 龍雪琴，關(guān)宏志.基于仿真的城市道路交叉口最小安全間距研究[J].武漢理工大學(xué)學(xué)報(交通科學(xué)與工程版)，2012,36(3):510-514.

[10] 卓亞娟.城市干道單行系統(tǒng)信號配時方法研究[D].西安：長安大學(xué),2013.

[11] GREENSHIELDS B D. A study in highway capacity[J]. Highway Research Board Proceedings,1935(1):448-477.