楊潘云,傅麗碧,方 捷,劉宇星
(福州大學(xué)土木工程學(xué)院,福州,350116)
隨著經(jīng)濟(jì)增長,各類娛樂場(chǎng)所數(shù)量逐年攀升,大型綜合性商場(chǎng)琳瑯滿目,吸引大量行人聚集。行人聚集使安全問題突出,行人動(dòng)力學(xué)的研究可為火災(zāi)等緊急情況下的安全疏散研究提供理論支撐。群組作為行人流中重要的組成部分[1,2],在行人動(dòng)力學(xué)研究中不可忽略。與獨(dú)立行人相比群組行人的運(yùn)動(dòng)特征更為復(fù)雜,群組行人不僅需考慮與周圍環(huán)境的相互作用,還需考慮同組成員之間的相互牽絆。屬于同一組的成員總是尋求以相同的速度和偏移角度移動(dòng)[3,4],以確保行走的舒適性和連貫性,其存在對(duì)速度、流量有明顯的負(fù)面影響[3,5,6]。
不同規(guī)模群組在行走過程中呈現(xiàn)不同的構(gòu)型。對(duì)于雙人組而言,行走構(gòu)型主要為肩并肩并列行走,該構(gòu)型在單向及雙向流的可控實(shí)驗(yàn)中進(jìn)一步得到證實(shí)[7,8],在雙向流中雙人組行人為減少與對(duì)向行人的摩擦?xí)霈F(xiàn)拉鏈效應(yīng),呈現(xiàn)錯(cuò)位行走的構(gòu)型[7],成員之間的距離則隨著密度的增加逐漸減小。三人組的構(gòu)型主要為“V”字形[2,9],該構(gòu)型及前后行人之間的距離不會(huì)隨密度的變化而發(fā)生改變[9]。四人組在行走時(shí)呈現(xiàn)“U”字構(gòu)型[2]。在高密度情況下,為減小與對(duì)向行人的沖突,三個(gè)或三個(gè)以上的群組會(huì)分裂單獨(dú)的個(gè)體或雙人組等小群體[2]。
已有學(xué)者將群組的運(yùn)動(dòng)特征考慮進(jìn)行人模型中,構(gòu)建了介觀模型、微觀模型[6,10-13]等。從這些模型中學(xué)者總結(jié)了一些現(xiàn)象,如合適的領(lǐng)導(dǎo)和穩(wěn)定的團(tuán)隊(duì)組成有助于減少疏散時(shí)間,減少與辦公室內(nèi)障礙物的碰撞[6],疏散過程中的行人橫向移動(dòng)可能會(huì)被垂直分組阻塞[13]等。在考慮群組運(yùn)動(dòng)特征的微觀模型中,Wang等[14]考慮了行人在運(yùn)動(dòng)過程中的轉(zhuǎn)向、避障、慣性等因素來構(gòu)建群組場(chǎng)模型,以再現(xiàn)群組運(yùn)動(dòng)特征。部分學(xué)者將群組成員劃分為領(lǐng)導(dǎo)者及跟隨者[6,10],領(lǐng)導(dǎo)者有絕對(duì)的領(lǐng)導(dǎo)權(quán),僅考慮周圍環(huán)境及行人的影響,跟隨者服從于領(lǐng)導(dǎo)者的決策。但對(duì)于大群組而言,與群組中某一位成員維持較近距離較為困難。領(lǐng)導(dǎo)者的概念在模擬雙人組行走特性時(shí)效果較好,但在模擬三人組、四人組等較大群組時(shí)存在一定的局限性。由此可見,涉及群組行人的動(dòng)力學(xué)研究具有現(xiàn)實(shí)意義,仍需要深入探究。本文基于前人對(duì)群組行人運(yùn)動(dòng)特征研究的基礎(chǔ)上,對(duì)群組場(chǎng)進(jìn)行拓展,并引入預(yù)測(cè)場(chǎng)進(jìn)行雙向行人流研究。
場(chǎng)域模型[15]由蟻群行為特征啟發(fā)而來,模型中介紹了靜態(tài)場(chǎng)、動(dòng)態(tài)場(chǎng)兩種基本勢(shì)場(chǎng)。靜態(tài)場(chǎng)與周圍環(huán)境相關(guān),反應(yīng)的是行人尋找最短路徑到達(dá)目的地的期望,在模擬過程中始終不變。動(dòng)態(tài)場(chǎng)與周圍行人運(yùn)動(dòng)軌跡相關(guān),反應(yīng)的是行人的跟隨作用,會(huì)隨著時(shí)間的推移進(jìn)行擴(kuò)散和衰退。除此之外為避免與周圍行人或障礙物發(fā)生碰撞,該模型還考慮了格子是否被行人占據(jù)及是否存在障礙物的問題。Kirchner和Schadschneider[15]依據(jù)以上假設(shè)制定的行人轉(zhuǎn)移概率公式見公式(1):
pij=Nexp(kDDij)exp(kSSij)(1-nij)ζij
(1)
式中:N:歸一化參數(shù),N=[∑exp(kDDij)·exp(kSSij)(1-nij)ζij)]-1;kD:動(dòng)態(tài)場(chǎng)敏感性參數(shù);kS:靜態(tài)場(chǎng)敏感性參數(shù);Sij:與距離相關(guān)的靜態(tài)場(chǎng);Dij:與周圍行人運(yùn)動(dòng)軌跡相關(guān)的動(dòng)態(tài)場(chǎng);ζij:障礙物參數(shù),若為障礙物則取0,否則為1;nij:行人占據(jù)參數(shù),若該位置被行人占據(jù)則為1,否則為0。
靜態(tài)場(chǎng)可以使用歐式距離計(jì)算:
(2)
其中(iTs,jTs)為目的地元胞坐標(biāo),(il,jl)為測(cè)量區(qū)域內(nèi)所有的元胞坐標(biāo),(i,j)為計(jì)算時(shí)的元胞坐標(biāo)。
動(dòng)態(tài)場(chǎng)Dij是一個(gè)隨時(shí)間變化的參數(shù),如螞蟻留下的化學(xué)軌跡那般會(huì)隨時(shí)間的推移而向周圍擴(kuò)散及自身產(chǎn)生衰退,Nishinari等[16]將公式進(jìn)行量化:
(3)
1.2.1 考慮群組凝聚力
場(chǎng)域模型中一個(gè)分支即為考慮群組運(yùn)動(dòng)特征的研究[10,11,17,18],如Müller等[18]將群組行人分為領(lǐng)導(dǎo)者及跟隨者,領(lǐng)導(dǎo)者遺留下來的波爾粒子m>1,作用于同組行人,領(lǐng)導(dǎo)者不受其他領(lǐng)導(dǎo)者的影響,而跟隨者受到領(lǐng)導(dǎo)者的影響。Qiu和Hu[19]在建立群組行人模型時(shí),將群組凝聚方式分為兩種,分別模擬了行人并列行走及跟隨領(lǐng)導(dǎo)者行走的構(gòu)型。魏曉鴿[20]定義群組成員傾向于向群組重心靠攏或跟隨群組中走在最前方的行人,依此提出了群組中心生成法、群組領(lǐng)導(dǎo)者生成法。由此可見,在群組重心或領(lǐng)導(dǎo)者位置處存在的凝聚力點(diǎn),能使群組成員凝聚在一起。
(4)
(5)
式中xmn和ymn分別代表格點(diǎn)(m,n)的X坐標(biāo)和Y坐標(biāo)。xi-sourcecell、yi-sourcecell為i群組重心的X坐標(biāo)和Y坐標(biāo)。
使用上述方法的模擬場(chǎng)景中,無法保障群組在運(yùn)動(dòng)過程中構(gòu)型的穩(wěn)定性,特別是對(duì)于四人組來說,受同組行人影響而處于不穩(wěn)定狀態(tài)情形較多。例如,圖1(a)中為三人組,該構(gòu)型無法滿足1號(hào)行人往正前方行走的空間需求,雖然是與群組保持較近的距離,但不是穩(wěn)定的構(gòu)型。圖1(b)、圖1(c)中為四人組,該構(gòu)型能夠滿足成員維持聯(lián)系及移動(dòng)所需空間,往重心方向靠近的行為反而會(huì)增加同組成員之間的沖突,因此4號(hào)、8號(hào)行人會(huì)選擇往正前方行走而不是為了要走在群組在重心位置而選擇往重心方向靠近。
按照日常生活經(jīng)驗(yàn)及學(xué)者的理論研究[2,20,21],在行人運(yùn)動(dòng)過程中,群組構(gòu)型會(huì)有一定的震蕩現(xiàn)象,而非嚴(yán)格的肩并肩行走,當(dāng)該形態(tài)足以維持群組交流等需求時(shí),會(huì)以該形態(tài)行走,圖2~圖4為結(jié)合前人研究結(jié)果[2,20,21]總結(jié)出的群組部分穩(wěn)定構(gòu)型示意圖。
(a)圖為不穩(wěn)定構(gòu)型,(b)、(c)為穩(wěn)定構(gòu)型
圖2 雙人組部分穩(wěn)定構(gòu)型
圖3 三人組部分穩(wěn)定構(gòu)型
圖4 四人組部分穩(wěn)定構(gòu)型
由此可見,群組成員不會(huì)無限地向群組重心靠近,因此本文引入預(yù)判距離閾值,改進(jìn)成員無限靠近群組重心的機(jī)械性問題,使行人能以維持較近距離的相對(duì)穩(wěn)定構(gòu)型行走。我們給定預(yù)判距離閾值,當(dāng)距離小于該閾值時(shí)則不再考慮群組場(chǎng)。
大群組在行走過程中比較容易分成小群組[17],因此本文將四人組的群組重心進(jìn)行偏移,在行走方向上向左向右移動(dòng)0.5格,劃分為兩個(gè)小重心,靠近原群組重心右方的行人使用右方偏移的重心,靠近左方的行人使用左方偏移的重心。以圖5(a)為例,上方兩個(gè)行人在原群組重心位置的上方則使用向上偏移的群組重心,下方兩個(gè)行人則使用向下偏移的群組重心。
圖5 重心偏移示意圖
1.2.2 考慮群組避讓行為
(6)
式中:Np:預(yù)測(cè)區(qū)域內(nèi)反向行人總?cè)藬?shù);Npg:預(yù)測(cè)區(qū)域內(nèi)反向行人中不同規(guī)模群組總?cè)藬?shù),下標(biāo)g表示群組的規(guī)模,取1時(shí)為單人組,2為雙人組以此類推;ε:群組規(guī)模影響權(quán)重值(見表1);Sp:距離阻力影響權(quán)重值(見圖6、圖7)。
群組規(guī)模的大小影響著行人的避讓決策,在有足夠空間避讓行人時(shí),獨(dú)立行人傾向于主動(dòng)避讓群組行人。群組規(guī)模越大,較小的群組避讓大群組的意識(shí)也越強(qiáng)?;谝陨霞僭O(shè)本文依據(jù)群組的規(guī)模賦予群組不同的權(quán)重值。
表1 群組規(guī)模影響權(quán)重值
存在沖突的行人之間的距離是行人避讓決策的考慮因素,發(fā)生沖突的行人之間距離越近壓迫感就越明顯,本文依據(jù)對(duì)向行人到預(yù)判位置的距離來確定距離阻力影響權(quán)重值。圖6、圖7為非后退及后退情況下距離阻力影響權(quán)重值的取值情況示意圖,行人運(yùn)動(dòng)方向?yàn)閺淖笸?,其中虛線圓為t時(shí)刻行人位置,實(shí)線圓為t+1時(shí)刻行人預(yù)判位置,淺灰色的區(qū)域?yàn)轭A(yù)判區(qū)域,預(yù)判區(qū)域?yàn)?×5元胞。
圖6 非后退情況下距離阻力影響權(quán)重值
圖7 后退情況下距離阻力影響權(quán)重值
1.2.3 拓展的群組模型
本文使用有后退的摩爾鄰居,行人可往8個(gè)方向運(yùn)動(dòng)或維持原位不動(dòng),具體見圖8。行人在行走過程中一般不會(huì)選擇往斜對(duì)角行走及后退,因此本文將5、6、7、8方向的靜態(tài)場(chǎng)進(jìn)行削弱處理,并采用不會(huì)連續(xù)后退兩步的后退模式。經(jīng)過多次試驗(yàn),當(dāng)5、6、7、8方向的靜態(tài)場(chǎng)值為1、2、3、4方向的0.5倍時(shí),能減少同組成員之間的沖突,維持穩(wěn)定的行走方向,因此本文5、6、7、8方向的靜態(tài)場(chǎng)為1、2、3、4方向的0.5倍。
將預(yù)測(cè)場(chǎng)與群組場(chǎng)考慮進(jìn)原始場(chǎng)域模型中構(gòu)建概率轉(zhuǎn)移公式(7):
qij=Nexp(kSSij+kDDij+
(7)
其余參數(shù)含義與基本場(chǎng)域模型相同,具體見公式1。
圖8 對(duì)行人周圍9個(gè)位置進(jìn)行編號(hào),0為行人占據(jù)位置
模型的框架如圖9所示,行人每一時(shí)間步的轉(zhuǎn)移概率均由靜態(tài)場(chǎng)、動(dòng)態(tài)場(chǎng)、預(yù)測(cè)場(chǎng)、群組場(chǎng)共同決定。由轉(zhuǎn)移概率預(yù)先判斷移動(dòng)的方向。在無沖突或競(jìng)爭(zhēng)成功,且不是連續(xù)后退兩步時(shí)往預(yù)判方向移動(dòng),否則保持原位不動(dòng)。依此循環(huán)至模擬結(jié)束。
為尋找符合行人運(yùn)動(dòng)特性的模型參數(shù),對(duì)上述模型進(jìn)行雙向流通道場(chǎng)景下的情景模擬。通道長為L,寬為W,每個(gè)格子采用0.4 m×0.4 m的大小。密度為通道內(nèi)行人總數(shù)除以通道面積。四人組、三人組、雙人組、單人組是行人流中的主要組成部分[1,2],因此本文僅研究這幾類群組,如圖10所示,初始時(shí)行人在通道兩端隨機(jī)產(chǎn)生。
圖11為在靜態(tài)場(chǎng)及群組的協(xié)調(diào)作用下模擬得到的平均速度。模擬場(chǎng)景為6×24個(gè)格子的循環(huán)通道,通道中設(shè)置有雙向行人,兩端行人均為等比例的單人組、雙人組、三人組、四人組,總共有80人,每個(gè)場(chǎng)景模擬十次,每次實(shí)驗(yàn)?zāi)M5 000個(gè)時(shí)間步。模擬發(fā)現(xiàn)靜態(tài)場(chǎng)對(duì)速度具有積極的影響,而群組場(chǎng)則為消極影響。當(dāng)速度在0.8以上時(shí)群組構(gòu)型穩(wěn)定性不佳,速度小于0.3時(shí)往目的地行走意識(shí)弱。因此本文模擬參數(shù)ks取值為3.5,kg取3。
Kirchner和Schadschneider[15]提到:在kS取1.0、10的情況下,kD均可以在取1時(shí)達(dá)到最小值。α、β取不同值的情況下kD亦可在取1時(shí)達(dá)到最小值,因此本文將動(dòng)態(tài)場(chǎng)敏感性kD參數(shù)取值為1,α、β均取0.3。
圖9 拓展場(chǎng)域模型框架
圖10 通道場(chǎng)景示意圖
圖11 kS及kg協(xié)調(diào)作用下的速度分析
本文依據(jù)行人行走過程中需要的前行空間,即移動(dòng)空間需求確定預(yù)判距離閾值。研究發(fā)現(xiàn)男、女性平均步幅分別為66.6 cm、60.6 cm[22],因此本文的預(yù)判距離閾值取值在1~2個(gè)元胞之間,又由于行人往正前方移動(dòng)時(shí),發(fā)現(xiàn)成員落后的反應(yīng)時(shí)間比較遲鈍,因此正前方的預(yù)測(cè)距離閾值較其他方向大。在使用預(yù)判距離閾值的模擬中,以從左往右行人為例,正前方即2方向(從右往左行人則為4方向)的預(yù)判距離閾值四人組、三人組為2元胞,單人組、雙人組為0元胞,其余方向的預(yù)判距離閾值四人組、三人組、雙人組均為1個(gè)元胞。
圖12為無干擾情況下的群組運(yùn)動(dòng)軌跡圖,在僅考慮靜態(tài)場(chǎng)、動(dòng)態(tài)場(chǎng)、群組場(chǎng)對(duì)群組的影響場(chǎng)景下能夠很好地模擬行人在無干擾(無對(duì)向行人沖突)下的群組并肩行走的運(yùn)動(dòng)特征。按照日常生活經(jīng)驗(yàn)及學(xué)者的理論研究[2,20,21],在行人運(yùn)動(dòng)過程中,群組構(gòu)型會(huì)有一定的震蕩現(xiàn)象,模擬中雙人組出現(xiàn)的輕微震蕩現(xiàn)象符合群組運(yùn)動(dòng)特征。
圖12 無干擾情況下群組運(yùn)動(dòng)軌跡圖
參考Suma等[23]的分析方法,本文將不往2方向(以從左往右行人為例,從右往左行人則為4方向)行走的時(shí)間步認(rèn)為是時(shí)間損失,將該時(shí)間損失進(jìn)行細(xì)化分類,圖13為時(shí)間損失類別示意圖,其中:
圖13 時(shí)間損失類別示意圖
(a)Lv為空元胞的時(shí)間損失,即2方向?yàn)榭赵麉s不往正前方行走的時(shí)間損失;
(b)Lc為反向行人的時(shí)間損失,即2方向被反向行人占據(jù)無法往正前方行走的時(shí)間損失;
(c)Lf為同向行人的時(shí)間損失,即2方向被同向行人占據(jù)無法往正前方行走的時(shí)間損失。
圖14為在6 m×8 m的通道中模擬收集得到的時(shí)間損失情況。模擬場(chǎng)景為包含混合群組的雙向行人流,兩端各30人,平均分布有單人組、雙人組、三人組、四人組。每一場(chǎng)景模擬1 000次。圖14(a)中Lc隨著ka的增大而降低,各規(guī)模群組均在ka取0.8前有大幅度的下降,而后下降幅度變緩,最終趨于穩(wěn)定。從群組規(guī)模來看,Lc隨著群組的增大逐漸減小。圖14(b)中,從群組規(guī)模來看,Lf隨規(guī)模的增大而增大,ka<2時(shí)不同規(guī)模群組的差異小于ka>2。圖14(c)中,不同規(guī)模群組的Lv值基本趨于一致。Lv、Lf、T的曲線波動(dòng)軌跡幾乎相同。ka增大時(shí),行人提前改變方向避讓對(duì)向行人的概率增大導(dǎo)致Lv增大。避讓對(duì)向行人時(shí)改變方向等行為,影響到同向行人的正常通行,與同向行人沖突增多導(dǎo)致Lf增大。同向行人正常通行受到了干擾,會(huì)增加通過的時(shí)間,最終導(dǎo)致了總的通過時(shí)間T增加。Lv、Lf、T在ka取值大于5時(shí)的震蕩現(xiàn)象是由于預(yù)測(cè)場(chǎng)所產(chǎn)生的避讓意識(shí)過強(qiáng),在沖突區(qū)域存在疏散失敗現(xiàn)象導(dǎo)致的。
圖14 ka變化時(shí)對(duì)時(shí)間損失及總通過時(shí)間的影響
我們對(duì)比使用預(yù)判距離閾值及不使用預(yù)判距離閾值時(shí)不同規(guī)模群組的構(gòu)型差異。兩個(gè)場(chǎng)景的預(yù)測(cè)場(chǎng)敏感性參數(shù)ka均為0.8。使用以上參數(shù)模擬在6 m×16 m通道中的雙向行人運(yùn)動(dòng),初始化時(shí)通道兩端各10人,即單人組、雙人組、三人組、四人組每種類型群組各一組,共重復(fù)進(jìn)行5 000次實(shí)驗(yàn)。將空間等分為0.05 m×0.05 m的小區(qū)域,記錄所有群組成員每一時(shí)間步以群組重心作為坐標(biāo)原點(diǎn)的相對(duì)位置,每次實(shí)驗(yàn)的平均相對(duì)位置作為群組成員的占據(jù)位置。以此方法統(tǒng)計(jì)出5 000次重復(fù)實(shí)驗(yàn)中以群組重心為坐標(biāo)的周圍區(qū)域被占據(jù)的次數(shù),如圖15所示。
(a、c、e中不考慮預(yù)判距離閾值,b、d、f中考慮預(yù)判距離閾值)
圖15(a)、圖15(c)、圖15(e)為不考慮預(yù)判距離閾值時(shí)繪制的相對(duì)位置占有次數(shù)圖,圖15(b)、圖15(d)、圖15(f)為考慮預(yù)判距離閾值時(shí)繪制的相對(duì)位置占有次數(shù)圖。由圖15可見,不同規(guī)模群組均為并肩行走出現(xiàn)的次數(shù)最高,兩種場(chǎng)景下均能模擬出群組內(nèi)部成員之間維持較近距離的運(yùn)動(dòng)特征,考慮預(yù)判距離閾值可以減少同組成員之間的沖突。圖15(a)中大量四人組成員在重心半徑1元胞的范圍內(nèi)凝聚,群組成員間沖突較為嚴(yán)重,而圖15(b)中當(dāng)考慮預(yù)判距離閾值時(shí),四人組成員在重心半徑1元胞的范圍內(nèi)凝聚較為分散,即成員間的沖突較少,因此構(gòu)型較為穩(wěn)定。圖15(c)、圖15(d)中三人組的構(gòu)型主要為肩并肩行走,“∧”、“∨”的構(gòu)型也可被觀察到。對(duì)于雙人組來說均能較好得維持肩并肩行走。
我們對(duì)比分析行人流中含有一類行人及含有多類行人的速度差異。圖16中的單純?nèi)航M指通道中僅存在一類群組。圖17中的混合群組指通道中同時(shí)存在單人組、雙人組、三人組、四人組,群組組數(shù)比例為1∶1∶1∶1。在長為24 m的循環(huán)通道中進(jìn)行模擬,每一種密度下的實(shí)驗(yàn)均重復(fù)20次,每次運(yùn)行5 000時(shí)間步。模擬中的期望速度為1 m/s,即每一個(gè)時(shí)間步為0.4 s,圖16、圖17中的速度為往正前方即2方向行走(以從左往右行人為例,從右往左行人則為4方向)的距離除以時(shí)間。
圖16 不同密度下單純?nèi)航M的速度
圖17 不同密度下混合群組的速度
從圖16、圖17中看出,兩種場(chǎng)景下的速度均隨著密度的增加逐漸減低,最終到達(dá)阻塞狀態(tài)速度為0。單純?nèi)航M中雙人組、三人組、四人組與單人組平均速度差值均大于0.510 m/s,而混合群組中的平均速度差值均在0.08 m/s以下,單純?nèi)航M的速度差異性大于混合群組,混合群組各密度下的速度差異性較微小。密度在1.70人/m2之前,兩種場(chǎng)景下不同規(guī)模群組的速度均為單人組>雙人組>三人組>四人組。密度在1.70人/m2之后,單純?nèi)航M場(chǎng)景中速度最快交替出現(xiàn)在三人組與四人組之間,單人組的速度小于雙人組。這一現(xiàn)象可能是由于,低密度時(shí)群組為維持構(gòu)型需要花費(fèi)一定的時(shí)間,群組規(guī)模越大花費(fèi)時(shí)間越多,則速度越慢。而在高密度情況下,沒有足夠的行走空間,行走速度緩慢,群組之間自成的小群體抱團(tuán)行走自動(dòng)分層,減小了行人之間的摩擦力,因此大群組的速度反而比小群組快,即群組的存在能夠提高通過效率,這與Zhang等[6]的結(jié)果相識(shí)。
本文在場(chǎng)域模型中引入了群組場(chǎng)及預(yù)測(cè)場(chǎng),在含有群組的雙向流中進(jìn)行模擬。結(jié)果表明增加預(yù)測(cè)場(chǎng)能夠模擬群組行人提前避讓的行為,預(yù)測(cè)場(chǎng)敏感性參數(shù)取值過大會(huì)影響通過效率。增加預(yù)測(cè)距離閾值與未考慮預(yù)測(cè)距離閾值相比可以減小群組成員之間的沖突,使行走構(gòu)型穩(wěn)定。在含有群組的通道中,行人的速度隨著密度增加而減小。在僅考慮一種類型群組時(shí),在密度較低情況下,群組規(guī)模越小速度越快,而在高密度情況下,最快速度在三人組、四人組之間交替出現(xiàn),即高密度場(chǎng)景中大群組的存在能夠提高通過速度?;旌先航M場(chǎng)景與單純?nèi)航M場(chǎng)景的速度有一定的差異性,群組存在對(duì)通過效率的影響值得深入研究。