基于出行特征的汽車客運(yùn)站最高聚集人數(shù)測算研究

唐秋生,許 浩,張永輝,熊 洋,孫齡波

(重慶交通大學(xué) 交通運(yùn)輸學(xué)院,重慶 400074)

0 引 言

汽車客運(yùn)站設(shè)計初期,因缺乏原始客流、數(shù)據(jù)獲取難度高、測算方法主觀性較強(qiáng),測算出的最高聚集人數(shù)存在偏差,導(dǎo)致在客運(yùn)站后續(xù)運(yùn)營中出現(xiàn)資源浪費(fèi)或緊缺的問題。為此需確定一套在最高聚集人數(shù)定義上,避免使用在設(shè)計初期無法獲取的數(shù)據(jù)進(jìn)行人數(shù)測算的方法。

汽車客運(yùn)站最高聚集人數(shù)是指按照設(shè)計年度平均日旅客發(fā)送量乘以固定百分比得出,何宇強(qiáng)等[1]在分析相關(guān)因素后提出概率法、最壞情景法、模擬法3種方法進(jìn)行人數(shù)預(yù)測,其中概率法是一種根據(jù)發(fā)車時刻表,假定到達(dá)旅客分布從而計算的客觀方法,也是后續(xù)相關(guān)問題研究中所使用的主流方法;葉玉玲等[2]對概率計算法進(jìn)行補(bǔ)充,綜合分析發(fā)車頻率與聚集人數(shù)的關(guān)系;張?zhí)靷サ萚3]指出規(guī)范中定義所存在的歧義,并在實地調(diào)研分析問卷結(jié)果后指出,旅客到站時間基本服從正態(tài)分布,在考慮列車晚點和旅客出行習(xí)慣的基礎(chǔ)上,采用概率法進(jìn)行人數(shù)預(yù)測;C.DING等[4]在此基礎(chǔ)上,繼續(xù)構(gòu)建客流預(yù)測模型,在考慮動態(tài)波動性后將其納入ARIMA模型進(jìn)行區(qū)間預(yù)測;王愛麗等[5]則繼續(xù)前人在行車時刻表上入手預(yù)測的基礎(chǔ)上認(rèn)為旅客到站時間服從一種復(fù)合負(fù)指數(shù)分布,根據(jù)北京南站的客流時空特征和行車時刻表完成仿真;陳喜春等[6]提出了一種基于改進(jìn)圖解法的方法,將傳統(tǒng)計算方法的圖表法結(jié)合旅客聚集過程完成測算;姚加林等[7]認(rèn)為同一車站因發(fā)車方向不同也將導(dǎo)致聚集人數(shù)產(chǎn)生差異。

從最高聚集人數(shù)研究中總結(jié)得出:概率法是一種比較科學(xué)的計算方法,但其中涉及旅客到達(dá)時間的概率分布函數(shù)以及參數(shù)等受其他因素影響。在此基礎(chǔ)上,李兆豐等[8]認(rèn)為日期和天氣特征對旅客到站時間存在影響,并以LSTM神經(jīng)網(wǎng)絡(luò)為框架進(jìn)行短時客流預(yù)測;狄智瑋等[9]認(rèn)為乘坐頻率、行李量、流程熟悉度、出行目的等特征對旅客出行時間選擇有較大影響;崔敘等[10]認(rèn)為區(qū)域路網(wǎng)密度、容積率、辦公設(shè)施、交通設(shè)施數(shù)量等環(huán)境因素與旅客出行時間存在非線性作用機(jī)制;J.LU等[11]認(rèn)為旅客在車站內(nèi)的聚集時間(到站時間)與站內(nèi)的設(shè)施服務(wù)水平有重要關(guān)系,并定義了一種基于AFC數(shù)據(jù)和關(guān)鍵設(shè)施旅客步行時間的SSDC概念,進(jìn)行旅客到站時間預(yù)測;T.STASKO等[12]認(rèn)為旅客的到站時間偏好受列車擁擠程度影響作用明顯,并嘗試采用FW方法將旅客按照最短路徑的思路分配至列車以此來計算旅客在站的逗留時間;V.CANTILLO等[13]認(rèn)為航空旅客主要考慮時間的主觀值、機(jī)票成本的彈性、旅行時間等因素進(jìn)行出發(fā)時間的選擇;LIU Yang等[14]在前人基礎(chǔ)上,綜合各方面因素,包括外部環(huán)境、時間相關(guān)性、空間特征和運(yùn)營特性構(gòu)建一種名為DEEPPF的深度學(xué)習(xí)模型進(jìn)行旅客流量預(yù)測。

綜上所述,針對最高聚集人數(shù)以及旅客出發(fā)時間(到站時間)的研究已存在較多成果,但在人數(shù)測算過程中可進(jìn)一步綜合考慮旅客到站時間對最高聚集人數(shù)產(chǎn)生的微觀影響。為此,筆者提出一種結(jié)合旅客出行特征的最高聚集人數(shù)計算方法,以重慶市陳家坪汽車站為例進(jìn)行計算。

1 最高聚集人數(shù)影響因素

采用概率法進(jìn)行最高聚集人數(shù)測算時,主要考慮旅客在候車廳內(nèi)的聚集和消散過程:其中旅客聚集過程主要與到站時間、概率分布、發(fā)車時刻相關(guān);消散過程主要與站內(nèi)設(shè)置的閘機(jī)數(shù)量、預(yù)留隊列空間等客觀因素有關(guān)。

1.1 集散規(guī)律

1.1.1 旅客聚集概率分布

旅客聚集概率分布是指旅客在何時到達(dá)客運(yùn)站候車廳并開始候車直到進(jìn)入發(fā)車程序開始登車的概率分布。為準(zhǔn)確獲得其分布密度函數(shù),由陳家坪汽車站信息辦公室通過調(diào)取站內(nèi)閘機(jī)刷卡數(shù)據(jù)獲得,包含乘客到達(dá)時間和所乘坐班車信息。筆者通過MySQL軟件調(diào)取所選班次的刷卡數(shù)據(jù)信息,將乘客到達(dá)時間和所乘坐班車發(fā)車時間相減換算得出乘客候車時間,終而得到120名乘客的到達(dá)時間統(tǒng)計數(shù)據(jù)。同時為保證數(shù)據(jù)時效性,所選班次盡量緊湊,故選擇3趟2021年8月15日由重慶始發(fā)的3趟客運(yùn)距離大于200 km,定員40人以上,發(fā)車時間分別為12:00、13:00、14:00的班次作為乘客到達(dá)時間數(shù)據(jù)收集范圍,如圖1。

圖1 旅客到站時間統(tǒng)計Fig. 1 Statistics of passenger arrival time

據(jù)觀察,原分布逼近正態(tài)分布,針對統(tǒng)計數(shù)據(jù)以標(biāo)準(zhǔn)正態(tài)擬合如圖2,擬合優(yōu)度R2為0.837。

圖2 旅客到站時間擬合Fig. 2 Fitting of passenger arrival time

擬合得到的正態(tài)分布概率密度函數(shù)為:

(1)

參數(shù)為:μ=42.96,σ=10.32。

1.1.2 旅客消散概率分布

旅客消散概率分布是指旅客在發(fā)車時間前,開始檢票后由候車廳向班車內(nèi)消散過程的概率分布,同時還存在部分陪站人員離站過程的概率分布。由于汽車客運(yùn)班次單趟運(yùn)輸量較小,站內(nèi)設(shè)置的檢票閘機(jī)基本充足,近似認(rèn)為是均勻分布。在汽車站工作人員幫助下對2021年8月15日3趟由重慶始發(fā)的班車旅客離站時間進(jìn)行統(tǒng)計分析得到,旅客消散過程自檢票時間開始起至發(fā)車時間,基本勻速進(jìn)行。故認(rèn)為,旅客消散函數(shù)為:

(2)

參數(shù)為:ε=3

1.2 發(fā)車時刻及上車人數(shù)

發(fā)車時刻表或行車計劃表是客運(yùn)站內(nèi)反映行車調(diào)度計劃的重要表單,其中可以獲得發(fā)車時間、發(fā)車間隔、行車?yán)锍?車輛信息的重要數(shù)據(jù),可根據(jù)發(fā)車時間表確定概率法中所需要的時間節(jié)點以此生成對應(yīng)分布進(jìn)行聚集人數(shù)的仿真計算;根據(jù)發(fā)車間隔可以粗略判斷站內(nèi)同時間聚集人數(shù)情況;針對車輛定員信息,可以確定到站人數(shù);針對行車?yán)锍炭梢宰龀鼋羞h(yuǎn)途班車的基本判斷,為后續(xù)確定出行特征提供數(shù)據(jù)基礎(chǔ)。

1.3 旅客到站時間

筆者針對旅客聚集規(guī)律進(jìn)行了一次實地調(diào)研,對3班列車的旅客候車時間進(jìn)行擬合得出近似標(biāo)準(zhǔn)正態(tài)分布的概率分布密度函數(shù)。但是分別對3趟列車候車數(shù)據(jù)進(jìn)行擬合發(fā)現(xiàn),其μ取值存在差異,如圖3,表明在不同客觀條件下,旅客到站時間存在明顯差異,該差異將直接導(dǎo)致站內(nèi)聚集人數(shù)變化。

圖3 分班次旅客到站時間擬合Fig. 3 Fitting of passenger arrival time in different shifts

結(jié)合相關(guān)文獻(xiàn),進(jìn)一步探討旅客到站時間與出行特征的影響機(jī)制,同時保證所選擇出行特征在客運(yùn)站設(shè)計初期的可獲取性[10-12]。

1.3.1 日內(nèi)時段影響

單日內(nèi)的發(fā)車時間可以被劃分為0:00—08:00區(qū)間的早間時段車輛,08:00—12:00區(qū)間的上午時段車輛,12:00—18:00的下午時段車輛,18:00—0:00的晚間時段車輛。根據(jù)旅客作為自然人的生理特性,在面對不同時段發(fā)車的班車時,會根據(jù)發(fā)車時間所在的具體自然時段調(diào)整到站時間。在面對早間時段發(fā)車的班車,因為城市公共交通還未開行,只能選擇網(wǎng)約車或自駕車前往場站,并且早間時段道路交通狀況良好,時間耗費(fèi)容易掌握,該時段班次的旅客更傾向靠近發(fā)車時間到達(dá)場站。

1.3.2 周內(nèi)日期影響

周內(nèi)日期主要劃分為工作日和非工作日。工作日開行的班車除開在通勤時間的早高峰和晚高峰時段,大部分時間交通狀況較好,旅客傾向于靠近發(fā)車時間到站;非工作日開行的班車因節(jié)假日自駕出行比例較高、交通狀況不易估計等原因,更傾向于提前發(fā)車時間到站。并且在工作日與非工作日的同一時段上,旅客到站時間也存在明顯差異。

1.3.3 客運(yùn)距離影響

客運(yùn)距離可作為出行目的判斷依據(jù)之一,同樣也可作為隨身行李量多少的判斷指標(biāo)?？瓦\(yùn)距離較短的班車主要服務(wù)于城市近郊,多為本地乘客周末出行度假、商務(wù)目的,其服務(wù)對象均對乘車流程以及時間掌握度較高,傾向于靠近發(fā)車時刻到站;客運(yùn)距離較長的班車多為外地游客返鄉(xiāng),本地游客出遠(yuǎn)門目的,其服務(wù)對象從行李量、乘車時間等方面來說均長于普通旅客,更傾向于提前發(fā)車時間到站。

1.3.4 季節(jié)氣候影響

季節(jié)氣候主要與當(dāng)?shù)丶竟?jié)溫度相關(guān),以重慶亞熱帶季風(fēng)性濕潤氣候為例,全年平均氣溫較高,單日溫差變化大。當(dāng)旅客在氣溫較高或較低的時刻出發(fā)時,更傾向于靠近發(fā)車時刻到達(dá)場站,以減少候車時刻。同理,在不同季節(jié)時,因溫度原因?qū)е侣每妥龀龅某霭l(fā)選擇不盡相同。

綜上所述,筆者選取日內(nèi)時段,周內(nèi)日期,客運(yùn)距離,季節(jié),溫度5個出行特征作為旅客到達(dá)時間預(yù)測依據(jù),如表1。

表1 特征選取Table 1 Feature selection

為進(jìn)一步探求該5個特征與旅客到站時間的影響機(jī)制,考慮到預(yù)測模型需要,于重慶市陳家坪汽車客運(yùn)站開展4次實地調(diào)研,具體調(diào)研內(nèi)容為旅客到站時間與所乘坐班次信息,在獲取每一位受訪旅客的車次信息后,通過到站時間與發(fā)車時間換算求得提前到站時間,同時記錄每位旅客受訪時當(dāng)前時段、日期、季節(jié)以及當(dāng)前氣溫,以完成樣本數(shù)據(jù)收集。實際獲得318份問卷,篩選部分缺失信息的問卷得到280份有效問卷。將各項特征數(shù)據(jù)按照表1進(jìn)行嚴(yán)格數(shù)值轉(zhuǎn)換,得到最終調(diào)研數(shù)據(jù)如表2。

表2 樣本節(jié)選Table 2 Sample excerpt

根據(jù)有效樣本進(jìn)行Pearson相關(guān)性檢驗。由表3可得,日內(nèi)時段對旅客到站時間影響最大,其次為客運(yùn)距離、季節(jié)、溫度等,周內(nèi)日期相關(guān)程度最小但相關(guān)系數(shù)也達(dá)到0.882 5。結(jié)果表明筆者所選擇的5個出行特征與旅客到站時間存在影響機(jī)制,為后文預(yù)測模型的建立提供理論基礎(chǔ)。

表3 Pearson相關(guān)性分析Table 3 Pearson correlation analysis

2 考慮出行特征的最高聚集人數(shù)測算

2.1 符號定義

根據(jù)模型中出現(xiàn)的符號,現(xiàn)定義如表4。

表4 符號定義Table 4 Symbol definition

2.2 旅客到站時間預(yù)測模型

針對前文所述5個出行特征,其與旅客到站時間的具體影響機(jī)制難以直接描述,筆者選用BP神經(jīng)網(wǎng)絡(luò)預(yù)測模型進(jìn)行建模。BP神經(jīng)網(wǎng)絡(luò)模型具有良好的分類識別能力,在足夠的樣本量支撐下也具有一定自適應(yīng)能力,訓(xùn)練完成后的神經(jīng)網(wǎng)絡(luò)可認(rèn)為是特征與結(jié)論的一個函數(shù)映射關(guān)系。因此,將日內(nèi)時段、周內(nèi)時期等5個特征作為輸入,旅客到站時間作為輸出,進(jìn)行神經(jīng)網(wǎng)絡(luò)訓(xùn)練。

輸入層:

(3)

隱含層:

(4)

其中,f(x)是激發(fā)函數(shù),一般采用Sigmoid函數(shù)為:

f(x)=1/(1+e-x)

(5)

輸出層:

(6)

2.3 聚集旅客量

根據(jù)旅客集散規(guī)律和行車時刻表,旅客到達(dá)分布為:

(7)

參數(shù)為:μ=Ti,σ=Ti/3

式中:σ=Ti/3是為了控制大部分旅客在此時間段內(nèi)到達(dá)候車廳。因正態(tài)分布密度函數(shù)在無窮大時間軸上積分時總和為1的特性,故無法保證所有旅客按規(guī)定時間范圍內(nèi)到達(dá)。引入修正因子將原概率密度函數(shù)修正至指定時間范圍總和為1,有修正過程為:

(8)

(9)

經(jīng)修正后概率密度函數(shù)則有k時刻內(nèi)進(jìn)站人數(shù)為:

(10)

2.4 消散旅客量

按照旅客消散規(guī)律,乘坐i趟列車的旅客以及陪站人員消散函數(shù)為:

(11)

可分別獲得旅客以及陪站人員在k時刻乘坐i趟列車消散人數(shù)為:

(12)

(13)

2.5 聚集人數(shù)測算

根據(jù)k時刻聚集旅客量和消散旅客量,可得:

(14)

H=max(Hk)

(15)

3 模型求解及算例

3.1 模型求解

Step1樣本數(shù)據(jù)處理。在進(jìn)行BP神經(jīng)網(wǎng)絡(luò)模型訓(xùn)練之前,需指定模型輸出變量以及輸入變量,本研究中該模型考慮出行特征對乘客到站時間的具體作用機(jī)制,故將樣本數(shù)據(jù)的第1列定義為因變量,第2列～6列定義為自變量。

Step2訓(xùn)練網(wǎng)絡(luò)。使用測試集數(shù)據(jù)訓(xùn)練神經(jīng)網(wǎng)絡(luò),設(shè)置參數(shù)如下:訓(xùn)練方法為“Trainlm”,神經(jīng)元個數(shù)3,訓(xùn)練集70%,驗證和測試集各15%,驗證方法選擇均方差法,輸出訓(xùn)練完成的網(wǎng)絡(luò)以及訓(xùn)練情況。檢驗集擬合優(yōu)度為0.930 6,訓(xùn)練結(jié)果較好,保留該網(wǎng)絡(luò)。

Step3旅客到站時間預(yù)測。根據(jù)訓(xùn)練好的神經(jīng)網(wǎng)絡(luò)進(jìn)行旅客到站時間預(yù)測,根據(jù)單日內(nèi)行車時刻表,確定每班列車的5個特征取值,輸入網(wǎng)絡(luò)進(jìn)行預(yù)測,得到每趟班車旅客傾向的到站時間。

Step4建立仿真模型。根據(jù)列車發(fā)車時刻以及旅客到站時間,生成若干符合發(fā)車條件的聚集正態(tài)分布和消散分布,以此模擬旅客到站的數(shù)量變化。

Step5計算聚集人數(shù)。導(dǎo)入班車定員、陪站人員系數(shù)、滿載率等系數(shù),最終計算得出日內(nèi)人數(shù)變化。

3.2 算 例

結(jié)合筆者提出的聚集人數(shù)計算模型,針對陳家坪汽車客運(yùn)站現(xiàn)有運(yùn)行情況,選取年內(nèi)發(fā)車班次較多的單日行車計劃表進(jìn)行仿真計算,共計發(fā)車215輛均由站內(nèi)始發(fā)。根據(jù)車型統(tǒng)計和經(jīng)驗公式,滿載率在運(yùn)營高峰時取1,陪站系數(shù)取1.3,車次定員取值如表5。

表5 車型定員Table 5 Vehicle model seating capacity

根據(jù)模型計算得出12:57時刻出現(xiàn)旅客聚集最大值,旅客數(shù)量為867人,按定義計算得最高聚集人數(shù)788人,陳家坪汽車客運(yùn)站單位時間旅客人數(shù)變化如圖4,一日內(nèi)旅客聚集數(shù)量變化如圖5。

圖4 旅客人數(shù)變化時刻統(tǒng)計Fig. 4 Statistics of changes in passenger numbers at different times

圖5 聚集人數(shù)變化Fig. 5 Change in the number of the gathered people

圖5中,08:00至11:00出現(xiàn)第1波高峰聚集,主因該時段有較多班次發(fā)車,于12:00至16:00出現(xiàn)日內(nèi)第2波高峰,該時段發(fā)車班列車型多為大高一級和大高二級,單車輛乘坐旅客數(shù)較多且大部分為長途出行,導(dǎo)致站內(nèi)人數(shù)保持在較高聚集水平。

如圖6為根據(jù)站內(nèi)刷卡數(shù)據(jù)繪制而成的人數(shù)變化圖。實際情況中,最高聚集人數(shù)出現(xiàn)在14:00為766人,與測算結(jié)果相差約1 h;出現(xiàn)該情況主要考慮因測算模型中針對旅客的消散過程設(shè)定為勻速消散,而實際情況中旅客的消散多為非線性消散,由此而產(chǎn)生滯后作用的干擾;在整體趨勢上,將原始數(shù)據(jù)列與測算數(shù)據(jù)列進(jìn)行相關(guān)性分析,求得復(fù)相關(guān)系數(shù)為0.761 3,說明測算模型結(jié)果較準(zhǔn)確,適用于汽車客運(yùn)站最高聚集人數(shù)測算。

圖6 實測人數(shù)變化Fig. 6 Changes in the measured number of people

4 結(jié) 論

1) 陳家坪汽車站旅客平均到站時間為29.7 min,提前15～20 min的旅客占21.07%,21～30 min的旅客占33.93%,31～40 min的旅客占30.71%,40 min以上的旅客占14.29%。

2) 汽車客運(yùn)站的旅客到站時間與出發(fā)環(huán)境的不同存在一定影響機(jī)制。根據(jù)出行特征的差異從而確定旅客到站概率密度分布,構(gòu)建適用于汽車客運(yùn)站的最高聚集人數(shù)測算模型,以此來避免設(shè)計前期基礎(chǔ)客流資料主觀性較強(qiáng)等問題。在新建客運(yùn)站無法確定行車時刻表時,同樣也可利用規(guī)劃發(fā)車線路、間隔等資料進(jìn)行聚集人數(shù)測算。

3) 以陳家坪汽車客運(yùn)站作為算例分析。陳家坪汽車站內(nèi)聚集人數(shù)最大值在12:57出現(xiàn),人數(shù)為867人,最高聚集人數(shù)為788人,考慮遠(yuǎn)景汽車客運(yùn)站轉(zhuǎn)型升級為旅游集散中心或是商務(wù)中心等因素帶來的客流量上浮,可能出現(xiàn)的最高聚集人數(shù)為1 024人。