王宗直，路 堃，陳 濤，姚娟娟

(1.清華大學 工程物理系 公共安全研究院，北京 100084；2.北京市地鐵運營有限公司地鐵運營技術(shù)研發(fā)中心 地鐵運營安全保障技術(shù)北京市重點實驗室，北京 102208)

0 引言

在現(xiàn)代都市中，地鐵作為一種高效、便捷、大運力的公共交通工具，承擔著城市中極其重要的交通運輸任務(wù)，很大程度上緩解了大型城市地面交通擁堵的情況，為城市居民日常生活提供了很大的便利。截至2015年底，全國共有37個直轄市或地級市已經(jīng)運營或正在鋪設(shè)地鐵線路，隨著中國城鎮(zhèn)化進程的不斷發(fā)展，地鐵的使用將越來越廣泛，每天搭乘地鐵出行的人數(shù)將不斷上升。以北京地鐵為例，北京地鐵是全世界日均運送旅客數(shù)量最多的城市軌道交通系統(tǒng)，日均流量約為一千萬人次，擔負著北京全市大部分公共運輸任務(wù)。

對于城市軌道交通運輸系統(tǒng)來說，最重要的一項任務(wù)就是保護旅客生命安全，雖然自倫敦地鐵投入使用以來，全世界范圍內(nèi)地鐵事故發(fā)生概率都極小，但是依然存在一定的安全風險。除了惡意破壞和設(shè)備故障兩大事故誘因外，人流量過大導(dǎo)致的人群密集、擁塞對軌道交通公共安全也帶來了極大的挑戰(zhàn)。尤其在中國都市中，人口密度極大，導(dǎo)致地鐵站內(nèi)擁擠的情況也較為嚴重，存在一定的安全隱患。

針對這一問題，本文提出了基于格子氣模型的地鐵站內(nèi)人流移動模型，通過數(shù)值模擬的方式，預(yù)算站內(nèi)人流擁堵情況，通過可視化的輸出結(jié)果顯示站內(nèi)人群的在未來某一時刻內(nèi)的分布情況以及關(guān)鍵位置發(fā)生擁堵的情況，實現(xiàn)預(yù)測預(yù)警功能。

1 使用模型

本文在格子氣模型的基礎(chǔ)上，根據(jù)地鐵站臺實際情況，以閘機作為分界點，搭建站內(nèi)人流移動模型?；舅枷霝槔妹商乜_方法，結(jié)合人群移動特性與統(tǒng)計學方法，利用實證數(shù)據(jù)修正模型參數(shù)，得到精確性符合要求的人群移動模型。

1.1 模型原理

在本文工作中使用的模型以格子氣模型為基礎(chǔ)，這一模型是由Muramatsu和Nagatani等[1-2]將擴散限制凝聚(DLA，Diffusion Limited Aggregation)模擬中的模型應(yīng)用到行人流中而搭建起來的。在這一模型中，首先要將研究區(qū)域進行離散化處理，分成W×L的矩陣，將行人抽象化為服從一定游走規(guī)律的隨機游走粒子進行處理，其滿足的移動規(guī)則如圖1所示。圖中，圓點表示即將移動的行人，圖中箭頭以及其上標注的概率表示可能移動方向和移動概率，被交叉覆蓋的格點表示被其他行人或障礙物占據(jù)，不能向其移動。這一模型滿足周期性邊界條件，從一側(cè)離開矩陣的行人將會從另一側(cè)重新進入，專門用于分析通道通行效率。

圖1 格子氣模型移動規(guī)則范例Fig.1 Example of moving regulation of lattice gas model

格子氣模型在應(yīng)用過程中也由研究者進行了一系列改進，為了解決格子氣模型移動速度離散化的問題，Jiang[3]另外引入一組概率以描述行人在一個時間周期內(nèi)跨越2個格點的情形，通過調(diào)整概率大小可以實現(xiàn)移動速度的變化，使用改進后的模型對通道通行效率進行了分析；隨著格子氣模型的發(fā)展，格點轉(zhuǎn)移也逐漸復(fù)雜，力求能更好地模擬行人運動，引入Moore型和擴展的Moore型轉(zhuǎn)移方式，但這2種轉(zhuǎn)移方式都會出現(xiàn)實際直線距離較短、程序模擬路徑較長的情況； Guo[4]提出了圓形概率轉(zhuǎn)移模型，解決了這一問題；在此基礎(chǔ)上，Hartmann[5]提出了六邊形格子氣模型，兼顧了行人路徑選擇、路徑距離的一致性和編程上的便利性，在近幾年逐漸被人重視并采用；Yu[6]，Zheng[7]等將格子氣模型與實際結(jié)合，利用該模型研究了房間內(nèi)的人群疏散問題；在此基礎(chǔ)上，學者們針對人群信息傳遞[8]、路徑選擇[9-10]、擁堵[11-12]以及在國內(nèi)場景中的實際應(yīng)用[13-15]進行了研究。本文利用動態(tài)參數(shù)的格子氣模型進行建模計算，因其需要設(shè)置每一處概略分布矩陣的值，設(shè)置參數(shù)過程較繁瑣，故編程工作量大，作為離散性的模型，格子氣模型在移動概率矩陣變化處的連續(xù)性較差，但是，此算法邏輯較連續(xù)性模型簡單，能夠保證一定的精確度和運算速度，比較適合用于解決本文中提出的快速短時人流模擬問題。

1.2 建模過程

本文以北京地鐵五號線與十號線的換乘站——惠新西街南口站作為研究對象，該站位于北京東北部，處在北三環(huán)與北四環(huán)之間，是重要的地鐵交通樞紐，也是北京地鐵最繁忙的站臺之一，建站以來人流壓力很大，一直通過更改部分扶梯上下行的方向、高峰期采取限流措施等方式對站內(nèi)人流進行疏導(dǎo)。2014年11月07日，在該站由于過度擁擠發(fā)生了1起車門事故，造成了1名乘客死亡。因此，對該站臺進行研究具有很強的必要性與代表性。

惠新西街南口站站內(nèi)設(shè)施較少，通過閘機進站后僅有2層，分別為位于地下一層的地鐵十號線站臺和位于地下二層的地鐵五號線站臺。閘機口均設(shè)置在地下一層，是北京地鐵站臺設(shè)計中較少的進站后不需上下樓即可直接乘車的地鐵站之一。五號線站臺設(shè)計上和十號線站臺垂直，位于十號線站臺的正下方，通過一段樓梯或扶梯就可以實現(xiàn)換乘，換乘時間大約為30 s，換乘時間很短。

本文在惠新西街南口站施工平面圖的基礎(chǔ)上，以閘機口為邊界，將站內(nèi)行人可以到達的區(qū)域進行網(wǎng)格化處理，共建模區(qū)域面積4 040 m2?？紤]到站臺整體情況，負一層為十號線側(cè)式站臺，將其劃分為南北2個模塊，負二層為五號線島式站臺，將其整體劃為1個模塊，因此，將整個站臺劃分為3個模塊進行處理。模塊內(nèi)采取網(wǎng)格化的模式建模處理，模塊之間利用通道相連，通道內(nèi)情況較單一，通過測算移動速度與通道長度，給出行人在通道中的轉(zhuǎn)移時間。

1.3 模型假設(shè)

為編寫算法的需要，除格子氣模型基本假設(shè)外，結(jié)合站臺實際情況，對人流特性額外做出了一系列假設(shè)：

1)模型對分析區(qū)域進行了網(wǎng)格化劃分，每一個網(wǎng)格的大小為0.5 m×0.5 m，每1個行人占據(jù)1個網(wǎng)格；

2)人群移動速度約為1.5 m/s，相對應(yīng)地，模型3個時間步長對應(yīng)實際運動時長為1 s；

3)在沒有阻礙的情況下，每個行人優(yōu)先選擇的路線一定是符合其效益最大化原則的，并且每個行人清楚地知道自己的目的地以及前往該目的地最短的路線；

4)在行人進入模型時，使用蒙特卡羅方法確定其目的地，并認為在應(yīng)用模型的過程中，行人出發(fā)-到達選擇(O-D對)概率分布在短時間內(nèi)不變；

5)不考慮列車滿員無法上車的情況，并簡化了乘客對屏蔽門的選擇過程，只考慮下車人流從屏蔽門進入站臺，沒有設(shè)計上車乘客的上車排隊過程；

6)從地下二層前往地下一層2個模塊的人群全部通過扶梯進行模塊之間的轉(zhuǎn)移，沒有行人會選擇樓梯進行上行的換乘。

1.4 算法邏輯

圖2 算法流程Fig.2 Algorithm flow chart

在本程序中所采用的算法流程圖如圖2所示，程序開始運行后，先完成矩陣和變量的初始化，在這一階段預(yù)設(shè)各個閘機口和列車屏蔽門進入站臺的人流情況、設(shè)置程序運行的時間周期，隨后進入時間周期循環(huán)。在每一個循環(huán)內(nèi)，對網(wǎng)格化的站內(nèi)各個格點進行分析，調(diào)取概率轉(zhuǎn)移矩陣，利用蒙特卡羅方法模擬實際的移動情況，將1個時間周期內(nèi)每個格點都計算完成后進入下一個格點的分析。

2 初步模擬與實證分析

2.1 輸出范例

程序輸出界面如圖3所示，其中每一個不同顏色的格點表示目的地不同的行人，利用圖形界面，直觀地輸出程序運行結(jié)果。圖中黑色部分為站臺及站廳中乘客不可到達的區(qū)域或通往其他模塊的通道區(qū)域，深色部分為建模區(qū)域，乘客可以在其中按移動規(guī)則進行移動，淺色方塊代表行人。

圖3 程序輸出結(jié)果界面Fig.3 The output interface of the algorithm

2.2 實證研究

為了確定模型中的各項參數(shù)，在模型基本框架搭建完成后，制定了實驗方案，組織實驗人員赴惠新西街南口站進行1次實地實驗。進行此次實驗的目的主要有以下4點：

1)驗證模型中各項移動參數(shù)的準確性；

2)對站臺內(nèi)部行人O-D對進行記錄；

3)測量閘機、扶梯、樓梯的最大人流量；

4)拍攝視頻資料用于模型精確度分析。

通過實地考察，得到各個通道的最大通行效率、通道通過時間、下班高峰期行人O-D對數(shù)據(jù)、現(xiàn)場視頻資料等重要數(shù)據(jù)。

2.3 模擬結(jié)果與實證對比

利用在現(xiàn)場采集得到的同步視頻，對視頻進行后處理以進行模型精確度分析，具體操作方式為將視頻中錄制得到的各個通道出入口的人流情況輸入對應(yīng)描述通道人流的矩陣中，運行程序，可以得到程序在各個關(guān)鍵通道、節(jié)點人流擁塞的情況，記錄模型中人流開始擁塞和擁塞結(jié)束的時刻，和視頻中實際產(chǎn)生擁塞的時間對比，間接確定模型的精確性程度,如表1所示。

表1 模擬結(jié)果與實際結(jié)果對比Table 1 The comparison of simulation result and actual result

運行結(jié)果表明，在地下一層南側(cè)和地下一層北側(cè)模塊中，模擬結(jié)果吻合較好，在地下二層的模擬與實際結(jié)果差別較大。精確度不足的原因主要是地下二層的地形相對較復(fù)雜，視野容易受阻，下車人流需要經(jīng)過一段時間確定路線，導(dǎo)致了在實際過程中人流移動更加平滑，不容易產(chǎn)生擁堵，而程序中假設(shè)了行人清楚自己的路線選擇，下車的一瞬間就開始朝一個目的地移動，導(dǎo)致了擁堵發(fā)生時間較早；在地下一層中，下車屏蔽門正對換乘通道，行人很容易確定其應(yīng)該前進的方向，所以程序的這一假設(shè)符合的情況較好，開始擁堵的時間估計較為準確。擁堵結(jié)束的時間比較準確，是由于擁堵結(jié)束時，擁堵在通道入口處的人流基本上都通過通道轉(zhuǎn)移，其通過時間和通道最大通行效率密切相關(guān)，因此，擁堵結(jié)束時間的準確性結(jié)果也證明了實驗測量所得通道最大通行效率數(shù)據(jù)的準確性。對于這一特定情形，綜合模型平均誤差控制在5 s以內(nèi)。

3 人流移動模擬

3.1 擁堵情況分析

圖4 站臺關(guān)鍵位置擁堵時間與人流量的關(guān)系Fig.4 Relation of congestion and pedestrian volume at key location

利用模型精確度較好的模塊對站內(nèi)關(guān)鍵位置進行擁堵分析，結(jié)果如圖4所示。一次列車進站人流量超過70 人就可能發(fā)生擁堵，尤其在換乘通道入口處，隨著人流量增大，擁堵時間逐漸增大。主要原因是換乘通道的通過能力有限，不能完成短時間大規(guī)模人流的移動，在通道入口處存在一定的安全隱患。

針對目前站臺換乘通道通行效率不足的問題，探討如果設(shè)計伊始將通道適當加寬對站內(nèi)通行能力的提升作用。對模型中與通道寬度變動相關(guān)的部分的轉(zhuǎn)移概率矩陣進行修改，將站臺布局進行一定變動，重復(fù)進行列車進站人流量對擁堵時間的影響分析。

考慮將換乘通道增寬2 m對車站擁堵情況的影響，控制進站人流情況相同，得到擁堵時間變化情況如圖5所示。從圖中曲線變化情況可以看出，將站臺拓寬后，不僅擁堵承載能力得到提升，在擁堵發(fā)生后，站臺擁堵持續(xù)時間縮短，站臺整體通行效率得到很大的提升，擁堵時間-人流曲線整體下移趨勢明顯。

圖5 通道擴寬與現(xiàn)有設(shè)計通行能力對比Fig.5 Comparison of traffic capacity between widened and actual passageway

將通道擴寬后，在人流量較少時的緩解擁堵的變化十分明顯，能夠?qū)⒄九_保持順暢通行的人數(shù)上限從70人/車次提升到110人/車次。對比用于調(diào)試模型在下班高峰期拍攝到的視頻資料，車站在工作日下班高峰期間某次列車到站下車人次約為120人/車次，所以將通道擴寬2 m就可以基本滿足車站換乘需求。

3.2 安全提升建議

對模型運行結(jié)果進行分析，惠新西街南口站在設(shè)計伊始沒有充分考慮站內(nèi)人流情況，對換乘人流量估計不足，導(dǎo)致設(shè)計換乘通道過于狹窄，不能滿足站內(nèi)人流進行快速換乘的需求，建議在未來規(guī)劃相似換乘布局時，適當拓寬站臺換乘通道的寬度,并將較低一層的島式站臺的寬度也適當加寬,從源頭緩解站臺擁堵情況。

對于包括惠新西街南口站以內(nèi)的北京地鐵人流密集站臺，由于設(shè)計上已經(jīng)存在一定的缺陷，所以會在人流高峰時期通過臨時設(shè)置的隔離欄改變站內(nèi)地形，起到控制人流移動模式的作用，保障乘客安全；或是在特殊時期(如公眾假期或特殊活動時期)關(guān)閉某些站臺，采取地鐵列車過站不停車的方式緩解站臺壓力。這些策略有一定的效果，但是對地鐵的運營效率和便利性造成了很大的影響，所以，應(yīng)當盡可能地減短限流措施或封閉站臺措施的時間，提高地鐵系統(tǒng)的運行能力和旅客的整體乘車體驗，在站臺建設(shè)已經(jīng)無法更改的情況下，利用本模型提前分析站臺的擁堵情況，在從10號線換乘5號線單次列車人流量超過某個閾值(如90～110人)后，決定采取相關(guān)安全措施，盡可能在保證安全的前提下減少車站限流措施。

4 結(jié)論

1)北京地鐵惠新西街南口站站臺設(shè)計中存在部分不足與缺陷，主要不足為十號線與五號線之間換乘通道設(shè)計通行能力較低，導(dǎo)致?lián)矶掳l(fā)生在換乘通道口，存在一定的安全隱患。

2)如果在前期設(shè)計中將換乘通道拓寬2 m，能夠使車站通行能力大大提高，基本滿足目前該站在工作日下班高峰期的換乘需求。

3)在未來地鐵換乘站類似設(shè)計方案確定前，可以應(yīng)用本模型的參數(shù)，搭建規(guī)劃站臺人流移動模型，合理估計站臺人流量，模擬測試站臺通行能力。

4)未來模型有望整合北京市軌道交通大數(shù)據(jù)平臺，獲取地鐵閘機數(shù)據(jù)、攝像頭影像數(shù)據(jù)、列車運行情況數(shù)據(jù)等，不斷完善模型精確度，增強模型功能。目前模型已經(jīng)搭好框架與算法部分，只需要獲得更加精確的輸入數(shù)據(jù)，或者結(jié)合大數(shù)據(jù)預(yù)測算法獲得的其他數(shù)據(jù)，就可以對站內(nèi)情況進行更廣泛準確的研究。

[1] MURAMATSU M, IRIE T, NAGATANI T. Jamming transition in pedestrian counter flow[J]. Physica A Statistical Mechanics & Its Applications, 1999, 267(S 3-4):487-498.

[2] MURAMATSU M, NAGATANI T. Jamming transition in two-dimensional pedestrian traffic[J]. Physica A Statistical Mechanics & Its Applications, 2000, 275(1-2):281-291.

[3] JIANG R, WU Q S. Pedestrian behaviors in a lattice gas model with large maximum velocity[J]. Physica A Statistical Mechanics & Its Applications, 2007, 373(36):683-693.

[4] GUO R Y, HUANG H J. A mobile lattice gas model for simulating pedestrian evacuation[J]. Physica A Statistical Mechanics & Its Applications, 2008, 387(2-3):580-586.

[5] HARTMANN D. Adaptive pedestrian dynamics based on geodesics[J]. New Journal of Physics, 2010, 12(4):1155-1208.

[6] YU X, CHANG R, ZHANG C. Evacuation of pedestrians using lattice gas model and floor field model[C] //International Conference on Audio, Language and Image Processing， IEEE, 2015.

[7] ZHENG Y, XI X, ZHUANG Y, et al. Dynamic parameters cellular automaton model for passengers in subway[J]. Tsinghua Science and Technology, 2015, 20(6): 594-601.

[8] HAN Y, LIU H. Modified social force model based on information transmission toward crowd evacuation simulation[J]. Physica A: Statistical Mechanics and its Applications, 2017(469): 499-509.

[9] PEREIRA L A, BURGARELLI D, DUCZMAL L H, et al. Emergency evacuation models based on cellular automata with route changes and group fields[J]. Physica A: Statistical Mechanics and its Applications, 2017(473): 97-110.

[10] ANTONINI G, BIERLAIRE M, WEBER M. Discrete choice models of pedestrian walking behavior[J]. Transportation Research Part B: Methodological, 2006, 40(8): 667-687.

[11] PEREIRA L A, DUEZMAL L H, CRUZ F R B. Congested emergency evacuation of a population using a finite automata approach[J]. Safety science, 2013, 51(1): 267-272.

[12] YANG X, DONG H, YAO X. Passenger distribution modelling at the subway platform based on ant colony optimization algorithm[J]. Simulation Modelling Practice & Theory, 2017(77):228-244.

[13] 史聰靈, 鐘茂華, 張嵐,等. 地鐵換乘車站客流疏運模擬及風險分析(2)——單通道換乘車站[J]. 中國安全生產(chǎn)科學技術(shù), 2011, 7(8):21-28.

SHI Congling, ZHONG Maohua, ZHANG Lan, et al. Simulation and risk Analysis on passenger transportation of metro transfer station(2)——Single passage transferring Station[J]. Journal of Safety Science & Technology, 2011, 7(8):21-28.

[14] 陳海濤, 劉占, 靳紅雨,等. 雙出口疏散吸引區(qū)域模型建立與模擬分析[J]. 中國安全生產(chǎn)科學技術(shù), 2017, 13(7):55-61.

CHEN Haitao, LIU Zhan, JIN Hongyu, et al. Modeling and simulation analysis on evacuation attraction region model for double exits[J]. Journal of Safety Science & Technology, 2017, 13(7):55-61.

[15] YE J, CHEN X, JIAN N. Impact analysis of human factors on pedestrian traffic characteristics[J]. Fire Safety Journal, 2012(52):46-54.