鄭 翔 劉 龍 鄒曉磊
(1.廣州地鐵設(shè)計(jì)研究院股份有限公司,510010,廣州;2.同濟(jì)大學(xué)道路與交通工程教育部重點(diǎn)實(shí)驗(yàn)室,201804,上海∥第一作者,高級(jí)工程師)
城市軌道交通車站的運(yùn)營(yíng)效率和安全性已成為運(yùn)營(yíng)方和專家學(xué)者廣泛關(guān)注的問題。目前我國(guó)多個(gè)城市的軌道交通線網(wǎng)客流進(jìn)入了快速增長(zhǎng)期,部分車站的客流壓力激增。針對(duì)車站的常態(tài)大客流,學(xué)者們?cè)诜治隹土鳉v史數(shù)據(jù)的基礎(chǔ)上,提出了客流組織的優(yōu)化方案[1]、減少進(jìn)站量的限流方法[2]、提升高峰時(shí)段車站客流疏解效率的列車運(yùn)行調(diào)整方法[3];針對(duì)車站的突發(fā)大客流,學(xué)者們分別提出了運(yùn)用視頻[4]、Wi-Fi(無線保真)[5]、手機(jī)信令[6]等技術(shù)手段及時(shí)獲取客流數(shù)據(jù),實(shí)時(shí)感知車站客流擁擠狀態(tài),并提出了基于三級(jí)大客流響應(yīng)標(biāo)準(zhǔn)采取客流管控措施[7]。然而,既有的大客流管控大多基于“見堵治疏”原則,常態(tài)大客流的管控措施常導(dǎo)致限流時(shí)段的延長(zhǎng),突發(fā)大客流應(yīng)急管控措施的實(shí)施時(shí)間也常滯后于需求,進(jìn)而造成一定的資源浪費(fèi),乘客的乘坐體驗(yàn)、車站的運(yùn)營(yíng)效率及安全性也受到了一定的影響。
為此,運(yùn)營(yíng)方希望能夠?qū)崿F(xiàn)大客流管控的“預(yù)堵治疏”,及時(shí)獲取站內(nèi)客流的聚集狀態(tài)及其10~15 min內(nèi)的客流發(fā)展趨勢(shì),以便及早采取管控措施。有學(xué)者基于短期客流預(yù)測(cè)方法建立了城市軌道交通車站客流限流動(dòng)態(tài)預(yù)警的數(shù)學(xué)模型,實(shí)現(xiàn)了高峰時(shí)段車站進(jìn)出站客流量及站臺(tái)客流密度的預(yù)測(cè)[8]。行人仿真作為重要的研究手段被廣泛應(yīng)用于城市軌道交通車站客流組織研究,用以還原不同場(chǎng)景下客流在站內(nèi)接受服務(wù)和集疏的過程[9-10]。但是,既有研究和應(yīng)用場(chǎng)景大多為基于車站設(shè)計(jì)和客流組織方案優(yōu)化的離線客流分析,可用于車站客流動(dòng)態(tài)組織和管控的在線實(shí)時(shí)客流態(tài)勢(shì)仿真和預(yù)警的研究成果還相對(duì)缺乏。
本文提出基于客流態(tài)勢(shì)的快速仿真推演方法,用以對(duì)車站的客流組織進(jìn)行調(diào)整,輔助決策者發(fā)現(xiàn)客流組織瓶頸和風(fēng)險(xiǎn)并及時(shí)進(jìn)行調(diào)整,提高車站的運(yùn)營(yíng)效率,確保車站管理安全。
客流態(tài)勢(shì)快速推演方法基于車站實(shí)時(shí)客流數(shù)據(jù)和歷史同期客流數(shù)據(jù),對(duì)車站未來10~15 min短時(shí)客流(以下簡(jiǎn)稱“短時(shí)客流”)的變化趨勢(shì)及變化過程進(jìn)行仿真推演,以用于車站客流聚集異常情況的研究判斷和及時(shí)干預(yù)。其系統(tǒng)架構(gòu)如圖1所示,主要包含3個(gè)模塊:實(shí)時(shí)數(shù)據(jù)和歷史數(shù)據(jù)接入模塊、面向快速推演的行人微觀仿真模塊、態(tài)勢(shì)推演結(jié)果復(fù)現(xiàn)及調(diào)整迭代模塊。這3個(gè)模塊分別發(fā)揮了數(shù)據(jù)處理、快速推演及決策支持的作用。
圖1 基于客流態(tài)勢(shì)快速推演方法的系統(tǒng)架構(gòu)
車站客流態(tài)勢(shì)快速推演的目的是獲取目標(biāo)時(shí)間段車站各區(qū)域的客流態(tài)勢(shì)變化過程,及時(shí)發(fā)現(xiàn)短時(shí)客流組織及客流分布的瓶頸及風(fēng)險(xiǎn),以便車站提出客流組織調(diào)整及管控措施,并在驗(yàn)證措施的合理性后迅速予以實(shí)施。基于客流態(tài)勢(shì)快速推演方法的技術(shù)流程如圖2所示。
圖2 基于客流態(tài)勢(shì)快速推演方法的技術(shù)流程
車站客流態(tài)勢(shì)推演的時(shí)機(jī)與耗時(shí)應(yīng)滿足城市軌道交通車站客流組織風(fēng)險(xiǎn)即時(shí)研判和應(yīng)急調(diào)整的要求,因此,將推演的目標(biāo)時(shí)間段Ta設(shè)為10~15 min,以保證短期客流態(tài)勢(shì)推演的精確度及客流管控措施調(diào)整在時(shí)間上的可行性。對(duì)某城市軌道交通車站在t0(推演目標(biāo)時(shí)間段開始時(shí)刻)至t1(推演目標(biāo)時(shí)間段結(jié)束時(shí)刻)時(shí)段內(nèi)的客流變化態(tài)勢(shì)進(jìn)行快速推演,其時(shí)間邏輯如圖3所示。
圖3 車站客流態(tài)勢(shì)快速推演方法的時(shí)間邏輯示意圖
圖3中:Tp、Tr和Ta在推演中可視為常量;Td為變量,受仿真場(chǎng)景復(fù)雜度和軟硬件運(yùn)算效率影響;Td1與Td數(shù)值接近;Tp1與Tp數(shù)值接近。要獲取Ta內(nèi)完整的客流態(tài)勢(shì)變化情況,需要在tr時(shí)刻提前導(dǎo)入對(duì)應(yīng)的實(shí)時(shí)客流數(shù)據(jù)和歷史客流數(shù)據(jù),并推演tr至t1時(shí)刻的客流變化過程。為確??土鲬B(tài)勢(shì)推演對(duì)客流管控措施調(diào)整的有效性,在推演總用時(shí)中預(yù)留了Td1,供車站決策者進(jìn)行二次仿真推演,以驗(yàn)證采取客流管控措施后站內(nèi)客流態(tài)勢(shì)變化情況。一般情況下Td1=Td,Tp1=Tp。
1.3.1 客流態(tài)勢(shì)快速推演方法的全流程效率
定義η為客流態(tài)勢(shì)快速推演方法的全流程效率,其計(jì)算式為:
(1)
由式(1)可知:要實(shí)現(xiàn)客流態(tài)勢(shì)快速推演,需使η∈(1, +∞)。當(dāng)η趨向于正無窮大時(shí),意味著客流態(tài)勢(shì)快速推演方法的效率趨于最大,車站初始狀態(tài)時(shí)刻tr趨向于tf1;當(dāng)η趨向于1時(shí),意味著tf1趨向于t1,即無法通過提前開始推演方式獲取Ta的客流變化態(tài)勢(shì)。η≤1時(shí),無法實(shí)現(xiàn)在線推演。
1.3.2 快速推演的仿真效率
定義λ為快速推演(包括一次推演和二次推演)的仿真效率,其計(jì)算式為:
(2)
由式(2)可知:λ受仿真場(chǎng)景復(fù)雜度及軟硬件運(yùn)算效率影響,主要由Td決定。η也可采用以下的計(jì)算式進(jìn)行計(jì)算:
(3)
由式(3)可知:因Tp、Tr和Ta可視為常量,則η由λ確定。在Ta、η和Tr均已確定的情況下,可得tr的計(jì)算式如下:
(4)
由上述分析可知,優(yōu)化快速推演的軟硬件條件,可縮短Td和Ts,提高λ和η,有利于提高仿真推演結(jié)果的精確度和快速推演的時(shí)效性。
綜合考慮快速推演過程的數(shù)據(jù)來源、預(yù)警方式、客流調(diào)整途徑及推演效率,本文認(rèn)為客流態(tài)勢(shì)快速推演方法的智能化水平可劃分5個(gè)等級(jí),如表1所示。
表1 車站客流態(tài)勢(shì)快速推演方法的智能化水平等級(jí)
表1中:在線、離線兩種推演方式的主要區(qū)別在于能否在tr提前開始推演,并獲取Ta的客流變化態(tài)勢(shì),即是否有足夠高的λ和η,可完成短時(shí)客流推演并立即用于車站客流組織調(diào)整?;谀壳暗募夹g(shù)水平,本文主要聚焦Level-3對(duì)應(yīng)的實(shí)時(shí)客流數(shù)據(jù)應(yīng)用和快速推演技術(shù)等方面的研究。
車站客流推演場(chǎng)景主要包括仿真推演目標(biāo)時(shí)段內(nèi)的乘客運(yùn)動(dòng)特征、車站設(shè)施設(shè)備布局、列車運(yùn)行方案、客流組織方案、客流初始狀態(tài)及客流到達(dá)分布等方面的特征。其中:乘客運(yùn)動(dòng)特征、車站設(shè)施設(shè)備布局、列車運(yùn)行方案、客流組織方案等場(chǎng)景特征應(yīng)在推演需求提出前預(yù)先處理;tr的初始客流狀態(tài)及tr—t1時(shí)段內(nèi)的客流到達(dá)分布則需在ts前快速生成。這些場(chǎng)景的生成主要運(yùn)用以下技術(shù):
1) 基于實(shí)時(shí)監(jiān)測(cè)客流數(shù)據(jù)和歷史客流數(shù)據(jù)生成tr的客流初始狀態(tài)。運(yùn)用自動(dòng)檢票閘機(jī)、視頻監(jiān)控、Wi-Fi探針等實(shí)時(shí)監(jiān)測(cè)手段獲取車站的實(shí)時(shí)客流數(shù)據(jù)。使用基于實(shí)時(shí)和歷史數(shù)據(jù)的站內(nèi)乘客個(gè)體出行選擇還原方法,處理得到乘客的位置、速度及流程環(huán)節(jié)等狀態(tài)參數(shù),進(jìn)而生成站內(nèi)客流群體在tr的空間、流程分布狀態(tài)和運(yùn)動(dòng)狀態(tài),并將這些數(shù)據(jù)快速導(dǎo)入至推演場(chǎng)景車站空間。
2) 基于歷史客流數(shù)據(jù)生成Ta內(nèi)的客流到達(dá)分布。結(jié)合歷史同期數(shù)據(jù)與仿真時(shí)段客流初始狀態(tài),運(yùn)用時(shí)間序列法和貝葉斯模型,分析客流發(fā)展規(guī)律并進(jìn)行預(yù)測(cè),生成Ta內(nèi)的車站客運(yùn)量及到達(dá)分布情況。
本文采用同濟(jì)大學(xué)自主研發(fā)的基于A*尋路算法和社會(huì)力模型的客流微觀仿真原型系統(tǒng),引入JPS(跳點(diǎn)尋路)算法和線程池技術(shù)(見圖4),在保證仿真精度的基礎(chǔ)上,提高尋路效率,改進(jìn)仿真架構(gòu),使η和λ滿足要求。JPS算法在A*尋路算法的基礎(chǔ)上,通過找尋障礙物附近可能改變乘客個(gè)體移動(dòng)方向的特定跳點(diǎn),使得尋點(diǎn)數(shù)量減少,進(jìn)而使尋路效率在不同場(chǎng)景下提高35%~100%。此外,基于線程池技術(shù)的仿真系統(tǒng)可實(shí)現(xiàn)多線程并行計(jì)算乘客個(gè)體位置變化,使仿真效率進(jìn)一步提高30%~80%。
圖4 線程池技術(shù)示意圖
基于上述技術(shù),客流態(tài)勢(shì)快速推演方法的η在常見場(chǎng)景下可達(dá)到4.5左右,足以滿足在線實(shí)時(shí)客流態(tài)勢(shì)仿真的時(shí)效性需求。
仿真推演完成后,需對(duì)Ta內(nèi)客流整體運(yùn)動(dòng)過程及主要特征進(jìn)行直觀的快速?gòu)?fù)現(xiàn),以便及時(shí)發(fā)現(xiàn)客流組織的瓶頸區(qū)域和擁擠風(fēng)險(xiǎn),研判車站短時(shí)客流的發(fā)展趨勢(shì)??土髅芏鹊臅r(shí)空分布變化情況最能直觀反映站內(nèi)的擁擠風(fēng)險(xiǎn)和發(fā)展趨勢(shì),本文優(yōu)化了客流狀態(tài)數(shù)據(jù)存儲(chǔ)及展示模塊,實(shí)現(xiàn)了仿真推演過程中客流空間密度數(shù)據(jù)的快速記錄、存儲(chǔ)、讀取,以及客流密度時(shí)空分布的動(dòng)態(tài)回放復(fù)現(xiàn)。
如果在推演結(jié)果中發(fā)現(xiàn)Ta內(nèi)存在較嚴(yán)重的客流疏散瓶頸,站內(nèi)擁擠風(fēng)險(xiǎn)有增加趨勢(shì)時(shí),車站決策者可以提出客流組織調(diào)整措施,并建立客流組織調(diào)整的推演場(chǎng)景。基于推演場(chǎng)景對(duì)客流組織調(diào)整方案下的客流態(tài)勢(shì)變化進(jìn)行快速仿真推演及數(shù)據(jù)再處理,以驗(yàn)證客流調(diào)整措施的實(shí)施效果。
本文以上海軌道交通8號(hào)線楊思站為例,選取工作日早高峰時(shí)段的8:30—8:40作為車站客流態(tài)勢(shì)推演的Ta,測(cè)試客流態(tài)勢(shì)快速推演方法的執(zhí)行效率和實(shí)施效果。
本案例取Ta=600 s,Tp1=Tp=10 s,Tr=60 s。經(jīng)測(cè)試可得Td≈68 s。根據(jù)式(2)求得λ=12,根據(jù)式(3)求得η≈5.23。根據(jù)式(4),取tr為08:25,即將Ta提前5 min作為初始時(shí)刻。根據(jù)實(shí)時(shí)客流數(shù)據(jù)生成車站的客流分布狀態(tài),并參考?xì)v史同期客流數(shù)據(jù)預(yù)測(cè)08:25—08:40該站的客流到達(dá)分布。
獲取tr的有效AFC(自動(dòng)售檢票)客流數(shù)據(jù),并根據(jù)事先已仿真獲取的乘客出行軌跡經(jīng)驗(yàn)數(shù)據(jù)集,生成該站的乘客空間分布及運(yùn)動(dòng)狀態(tài)。將這些數(shù)據(jù)導(dǎo)入仿真系統(tǒng),可得到tr的實(shí)時(shí)客流分布,如圖5所示。
圖5 楊思站客流實(shí)時(shí)分布截圖
考慮乘客動(dòng)態(tài)行走與靜態(tài)等候的空間服務(wù)水平評(píng)價(jià)差異,本文對(duì)原來僅適用于動(dòng)態(tài)行走場(chǎng)景的Fruin行人空間服務(wù)水平評(píng)價(jià)指標(biāo)加以改進(jìn),以粒度為50 cm(長(zhǎng))×50 cm(寬)的底層地圖柵格進(jìn)行車站客流密度統(tǒng)計(jì),并評(píng)估站內(nèi)空間的服務(wù)水平等級(jí)。在車站平面圖中,用熱力圖動(dòng)態(tài)顯示站內(nèi)客流密度的變化情況。
通過系統(tǒng)內(nèi)置的瞬時(shí)密度動(dòng)態(tài)復(fù)現(xiàn)技術(shù),可查看車站各時(shí)刻行人位置及密度分布,以輔助車站管理人員快速了解客流態(tài)勢(shì)變化情況,識(shí)別可能出現(xiàn)的客流組織瓶頸及擁擠風(fēng)險(xiǎn),研判客流的變化趨勢(shì)。針對(duì)可能發(fā)生客流擁擠的位置和時(shí)段,車站人員應(yīng)根據(jù)車站應(yīng)急處置預(yù)案提前進(jìn)行客流管控和客流組織調(diào)整。
楊思站站廳層一次推演后客流態(tài)勢(shì)變化展示界面如圖6所示。由圖6可知:站廳層左側(cè)樓扶梯區(qū)域、1號(hào)口走行通道及左側(cè)安檢通道區(qū)域易出現(xiàn)高密度行人聚集,可采取通道分流、安檢前繞行及樓扶梯前繞行等措施,以降低擁擠風(fēng)險(xiǎn)。
圖6 楊思站站廳層一次推演后客流態(tài)勢(shì)變化展示界面截圖
對(duì)采取客流組織調(diào)整措施后的楊思站開展全站客流態(tài)勢(shì)的二次快速推演,得到站廳層客流態(tài)勢(shì)二次推演后展示界面截圖,如圖7所示。由圖7可知:實(shí)施客流組織調(diào)整措施后,上述3個(gè)易擁擠區(qū)域的客流聚集情況得到了有效改善。
圖7 楊思站站廳層客流態(tài)勢(shì)二次推演后展示界面截圖
推演各階段的用時(shí)如表2所示。由表2可知,二次推演用時(shí)為67 s。完整的客流態(tài)勢(shì)推演從tr提出需求至tf1二次推演完畢,共計(jì)用時(shí)155 s,可滿足車站客流態(tài)勢(shì)在線推演和運(yùn)營(yíng)優(yōu)化的時(shí)效性要求。
表2 推演各階段用時(shí)表
本文針對(duì)城市軌道交通車站大客流管控問題,提出了基于多源客流數(shù)據(jù)和微觀仿真技術(shù)的車站客流態(tài)勢(shì)快速推演方法。該方法構(gòu)建了客流快速推演的系統(tǒng)架構(gòu),梳理了客流態(tài)勢(shì)推演和客流組織調(diào)整的完整流程,解決了客流快速路徑搜索和多線程仿真推演等關(guān)鍵問題。以上海軌道交通8號(hào)線楊思站工作日早高峰客流為案例,測(cè)試了客流態(tài)勢(shì)快速推演和客流組織措施調(diào)整的可行性和有效性。本文的研究成果在廣州地鐵22號(hào)線廣州南站站的互聯(lián)互通換乘服務(wù)客流態(tài)勢(shì)推演及優(yōu)化項(xiàng)目中得到了成功應(yīng)用。該方法可較高效地完成車站10~15 min短時(shí)客流態(tài)勢(shì)推演,并為車站客流組織方案的調(diào)整提供輔助決策。
經(jīng)驗(yàn)證,該方法總體達(dá)到了本文所提Level-3級(jí)的車站客流態(tài)勢(shì)推演技術(shù)水平。下一階段可基于本文所提的客流態(tài)勢(shì)快速推演系統(tǒng)架構(gòu),進(jìn)一步對(duì)智能化客流擁擠風(fēng)險(xiǎn)識(shí)別、發(fā)展趨勢(shì)研判和動(dòng)態(tài)預(yù)警方法等進(jìn)行深入研究,使該方法達(dá)到智能化和自動(dòng)化程度更高的車站客流動(dòng)態(tài)推演層級(jí)水平,實(shí)現(xiàn)車站客運(yùn)組織方案的自適應(yīng)調(diào)整和車站大客流的智能化管控。