任競穎
(商丘師范學(xué)院 計算機科學(xué)系,河南 商丘 476000)
近年來,在各種公共場所的人群聚集活動越來越多,各種事故和緊急事件頻繁發(fā)生,直接威脅著人民的生命財產(chǎn)安全并帶來不同程度的經(jīng)濟損失。故對于大量人群聚集的場所,應(yīng)預(yù)先對人員的疏散方案進行充分地準(zhǔn)備,做好各種應(yīng)急的疏散預(yù)案。地理信息系統(tǒng)(GIS)[1]在信息管理、空間分析和可視化表現(xiàn)等方面具有強大的優(yōu)勢,可以為人群的疏散模型研究提供全新的技術(shù)支持。
人員疏散量化的基礎(chǔ)是對自然現(xiàn)象的觀測,疏散仿真系統(tǒng)也正是以幾條對疏散人群的觀測結(jié)果為建模的基礎(chǔ)的。從大量的分析結(jié)果可以看出,其中的一些是疏散人群行為的根本特征,而另一些則是由人的一些特征所導(dǎo)致的結(jié)果。為了盡量真實地模擬疏散人群的動態(tài),必須對模型進行一些假設(shè),從而體現(xiàn)出那些根本的、起決定性作用的因素的影響。
在疏散模型中,主要將體現(xiàn)以下幾條人群疏散時的特征[2]:
1)在疏散時,個人總是想選擇到最近的出口的最短路線;
2)疏散人員是清醒的,在疏散開始的時刻,同時同刻井然有序地進行疏散,且在疏散過程中,不會中途返回選擇其他疏散路徑;
3)個人十分不愿意偏離其意愿的行進方向,即使改方向已經(jīng)比較擁擠了,個人也不愿意繞路;
4)在一定的情況下,個人十分容易受周圍人的行為的影響,從而表現(xiàn)出一種盲從行為;
5)火災(zāi)的發(fā)展趨勢、指示信息和命令對疏散的人有很大的影響。
在模型中,建筑平面被劃分成大小相等的正方形單元格,格子邊長為0.4 M[3]。待仿真的區(qū)域為被墻和出口包圍的單元格集合,而待疏散的每個人都用一個單獨的Agent來表示(在仿真術(shù)語中,作為研究對象的個人稱為Agent)[4]。
在空間中,文中主要考慮5種不同的元素:墻(包括障礙物),人,門,走道,火災(zāi)發(fā)生區(qū)域。
因此,加上單元格不被占據(jù)的情況,單元格共有5種狀態(tài):被墻(或被障礙物)占據(jù),被人占據(jù),被門占據(jù),被火災(zāi)區(qū)域占據(jù),單元格為空。5種情況必居其一,且只能居其一。一個人也只能占據(jù)一個單元格。因此,人的位置不僅可以由其坐標(biāo)來顯示,也可以由其所處的單元格編號來表示。
對于空間及空間中的元素,采用圖層(FeatureLayer)來表示,主要包括格網(wǎng)圖層(線層,主要用來作為空間編輯圖形的參考)、墻圖層(面層)、火災(zāi)發(fā)生區(qū)域圖層(面)和Agent圖層(點層)等。
在疏散過程中,個人的行為要受諸多因素的影響。大體上說,這些因素可以分為兩大類。第一類,靜態(tài)的物理場景(即建筑物幾何結(jié)構(gòu))對人行為的影響[5];第二類,動態(tài)的周圍環(huán)境對人的影響。
首先,建筑物的空間結(jié)構(gòu)是人做出行為判斷的最根本依據(jù),其中出口的位置又是最重要的。疏散時人總是要選擇最近的出口,在模型中用一個參量—“距離場”來描述。對在仿真區(qū)域中任意的不為障礙物所占據(jù)的單元格,都有一個與之對應(yīng)的靜態(tài)場強值,它的大小為該單元格到達最近的出口的最短距離。顯然,對于表示出口的單元格,表示個人的Agent只要始終沿著靜態(tài)場的梯度方向,從場強高的地方移向低的地方,它的行動路徑就是到達最近出口的最短路徑。在本模型中的靜態(tài)場中,一般給墻和被障礙物占據(jù)格設(shè)置比較大的場強值。這種用一個場來描述建筑空間信息的方法已經(jīng)為很多疏散仿真模型所采用,幾乎所有的離散微觀仿真模型和一部分連續(xù)模型都采用了類似的方法。 對于一個房間多個出口的空間結(jié)構(gòu),模型中我們?nèi)匀灰噪x房間門的距離確定其靜態(tài)場,通過比較到多個出口的距離,按照一定的權(quán),給門分配一定的單元格,在此單元格內(nèi)的人,則從該門出去。
其次,疏散中人的行為還要受一些別的因素的影響。比如說,某個地方著火了或者有清晰的指令要求人們到達某處,這些都會對人產(chǎn)生很大的影響。由于這些因素會隨著時間的流逝而變化,有可能哪個地方不能去,也有可能那個地方是人們要去的,在模型中統(tǒng)一以另一個參量—“動態(tài)場”來描述。與靜態(tài)場類似,每一個不為障礙物所占據(jù)的單元格,都有一個與之對應(yīng)的動態(tài)場強值。某單元格對個人的吸引力越高,則該格的場值越小(與距離場統(tǒng)一)。比如說,在疏散過程中有指令要求人員疏散至某處,則該地點的動態(tài)場強值就小。在仿真時,就會吸引更多的Agent到達該處[5]??傊谀P椭?,所有的對疏散的影響因素都被用“距離場”和“動態(tài)場”描述。
當(dāng)房間中的疏散人員離開房間后,連接房間和樓層出口之間的部分都屬于走道部分。對于復(fù)雜的建筑物來說,走道可能是多個入口、一個出口的空間,也可能是多個入口、多個出口的空間。入口是指與上一疏散通道相連的各房間單元,出口是指樓層的出口。
當(dāng)走道是多個入口、一個出口時,所有從房間等通道單元出來的人,都進入了走道部分,他們的移動目標(biāo)都是樓道出口。人們只需要按照離樓道出口最短的路線移動,就可以完成疏散動態(tài)的模擬。但是,當(dāng)走道是多個入口、多個出口時,即樓道出口有多個,那么所有人并不一定從一個出口出去,問題就復(fù)雜很多,國內(nèi)在這方面的研究成果非常少。此時,需要考慮哪些房間的哪部分人會從哪個出口出去,即需要給各個樓道出口分配一定的通過人數(shù)。文中根據(jù)房間的出口處的單元格到各個樓層出口的通道長度來定權(quán)分配,并確定房間內(nèi)疏散人員離開房間后會以哪個樓層出口為移動的目的地。
人員的行進速度不僅與火災(zāi)事故現(xiàn)場的空間結(jié)構(gòu)狀態(tài)和事故狀態(tài)有關(guān),而且受疏散人員年齡、性別、生理、心理及工作狀態(tài)等因素的影響,與人員的著裝及荷載情況、疏散人員的分布特點以及人員相互之間的社會關(guān)系有關(guān)[3]。由于時間和能力有限,將人群僅分為男人、女人、小孩和老人這4大類,且未對人群種類的速度與密度之間的關(guān)系以及各種影響人員速度的因素做深入的研究。假定人群的速度為常量值,如男人和女人的步行速度為1.2 m/s,用半徑為0.2 m的小圓表示。小孩和老人的步行速度為0.8 m/s,用半徑為0.12 m的小圓表示。通過對各種類型的人的屬性參數(shù)的計算,顯示其在建筑空間中的位置和動態(tài)疏散的情況。人的屬性參數(shù)包括人所在單元格的狀態(tài)、單元格的靜態(tài)場值、疏散人員的速度等多種屬性。
國外學(xué)者Fruin建議一個人占據(jù)0.45 m×0.61 m的矩形面積[6]。根據(jù)我國人體的體厚與肩寬測量數(shù)據(jù):按橢圓計算時人體的平均投影面積為0.164 m2;按矩形計算時人體的平均投影面積為0.197 m2。因此,從以上人體投影面積的相關(guān)數(shù)據(jù),取一個中間量,選擇0.4 m×0.4 m的正方形方格為疏散人員的投影面積[7]。
微觀模型的重點在于單個人的建模,這也真是本模型的核心。本模型對于單個人,采用一個點圖層來表示,用人的初始位置的單元格編號來標(biāo)識人。同時,考慮到人移動的速度,定義人的4方向為鄰居,如圖1所示,在仿真模型中,任意人在下一個時間段內(nèi)可能的位置有5個,即原來的位置加上相鄰的4個單元格。人依照一定的規(guī)則或算法在這5個格子中選取一個作為下一步所處的位置。
圖1 人的四方向行走圖Fig.1 The four direction map
對于某個人而言,它依照自己所處的單元格子加上相鄰4格共5格的狀態(tài),從中選擇一個單元格作為自己下一時間段所處的位置。假設(shè)它在t時刻的位置的單元格編號為(i,j),則其他相鄰的 4個格子的單元格編號分別為 (i,j-1)、(i-1,j)、(i,j+1)和(i+1,j)。 若 5 個格子中某個單元格被別的人或者障礙物占據(jù),稱之為不可達格,賦予其對應(yīng)的屬性值;否則稱為可達格,賦予其對應(yīng)的屬性值。
首先所有個體做出決策,對應(yīng)的單元格編號和將移動到的下一個單元格的編號存放到一個數(shù)組中,當(dāng)兩個個體同時走向同一單元格時,規(guī)定只有一個個體不移動,而另外一個個體不移動,保持原來的位置。
模型是以元胞自動機更新規(guī)則為基礎(chǔ),結(jié)合了多自主體、人工智能等研究思路建立起來的。模型的基本框架是將建筑物平面進行均勻的網(wǎng)格劃分,每個網(wǎng)格或被墻壁或其他障礙物占據(jù)、或被人員占據(jù)、或者為空。每個網(wǎng)格必須是這3種狀態(tài)中的一種,而且只能有一種狀態(tài)。模型中人員的運動將遵循一定的局部規(guī)則并具有一定的智能性,其行為會隨周圍環(huán)境的變化做出相應(yīng)的調(diào)整。
模型中每個元胞對應(yīng)0.4 m×0.4 m的空間,規(guī)定每個時間步每個人員都只能動一格。模型局部的運動規(guī)則在每個時間步有兩個基本問題需要解決,一是路線的選擇問題,二是如何解決當(dāng)多于一人同時競爭一個空格點時產(chǎn)生的沖突。模型中,賦予每個人員所在網(wǎng)格一個靜態(tài)場值,場值越小,表明該元胞(即網(wǎng)格)對人員的吸引力的越大,各出口點的場強值為0,對人們的吸引力最大。所有人員將根據(jù)其周圍鄰域內(nèi)的所有元胞的狀態(tài)和靜態(tài)場值來選擇一個相鄰的元胞作為下一步的移動目標(biāo),當(dāng)多個人員競爭同一個網(wǎng)格時,只能有一個留下,其他的回到原來的網(wǎng)格。如圖2所示。
圖2 四方向行走及其行走概率Fig.2 Four direction and walking probability
用元胞自動機模擬一個物理過程,省去了用微分方程作為過渡而直接通過指定簡單的規(guī)則來模擬非線性物理現(xiàn)象[3]。也就是說,元胞自動機具有利用簡單的、局部規(guī)則的和離散的方法描述復(fù)雜的、全局的、連續(xù)系統(tǒng)的能力。
疏散模型在應(yīng)用前,需要驗證以確定其預(yù)測功能。火災(zāi)中的人員疏散是一項十分復(fù)雜的過程,即使試驗條件控制得非常好,也需要進行重復(fù)試驗。因為,即使是重復(fù)進行由多人參加的真實疏散演習(xí),得出的結(jié)果也不會相同。由于篇幅有限,僅對單一房間內(nèi)人員的疏散進行描述。
單一房間是指疏散空間只有一個出口的單元建筑物空間。下面選用一個長6 m、寬4 m的單一房間,出口寬度為1.2 m,隨機分布了10個人的狀態(tài)。其在幾種時刻下的疏散狀態(tài)如圖3所示。
從圖3中,可以看出,對于人數(shù)來說,24 m2的空間相對很大,人口密度僅為0.42人/m2,總共疏散時間為6 s,人群沒有發(fā)生過度擁擠的情況。同時,可以發(fā)現(xiàn)模型可以很好的重現(xiàn)實際觀察到的“同向人同道”現(xiàn)象。因此,該模型可以在行人運動和人員疏散研究工作中起到相應(yīng)的作用。
圖3 單一房間人員疏散圖Fig.3 Single room evacuation map
由于實際工作的需要和其理論上的價值,對人群運動規(guī)律,特別是疏散中人群的運動規(guī)律的研究,近來逐步成為公共安全研究的一個熱點。最主要的研究方法是建立疏散仿真系統(tǒng),對各種情況下的運動做出模擬。文中在現(xiàn)有的各種仿真模型的基礎(chǔ)上,建立了一個基于GIS的微觀離散疏散仿真模型,并在此模型基礎(chǔ)上開發(fā)了相應(yīng)的仿真系統(tǒng)。通過仿真系統(tǒng)的試驗,可以看出通過仿真來研究人群的疏散有很大的優(yōu)勢。人員疏散本身是比較復(fù)雜的,涉及到人的心理素質(zhì)、文化背景、教育、生活習(xí)慣等難以量化的因素,這些影響因素很難準(zhǔn)確地用數(shù)學(xué)模型來進行描述,必然造成模擬計算結(jié)果的偏差,仿真軟件不可能完全模擬人員在火災(zāi)過程中的行為,模擬得到的疏散時間不可能做到十分精確,強調(diào)某一個模型的計算精度在幾秒鐘之內(nèi)是沒有必要的,值得重視的是應(yīng)該對不同疏散場景的比較,以找出可能出現(xiàn)的最不利位置,從而優(yōu)化疏散路線的設(shè)計。
[1]龔建雅.地理信息系統(tǒng)基礎(chǔ)[M].北京:科學(xué)出版社,2001.
[2]Owen M.The EXODUS evacuation model applied to building evacuation scenarios[J].Fire Engineers Journal,2010,20(7):27-28.
[3]陳靜.基于ArcGIS的動態(tài)顯示模型開發(fā)與應(yīng)用[D].上海:華東師范大學(xué),2009:34-36.
[4]蔣波濤.ArcObjects開發(fā)基礎(chǔ)與技巧—基于VB.NET[M].武漢:武漢大學(xué)出版社,2006.
[5]Lo S M,F(xiàn)ang Z,Lin P,et al.An evacuation model:the SGEM package[J].Fire Safety Journal, 2009(39):16-17.
[6]胡廣霞.火災(zāi)?;▽Υ罂臻g建筑人員疏散安全評估[J].沈陽航空工業(yè)學(xué)院學(xué)報,2010,2(3):27-29.HU Guang-xia.The fire modeling method for large space building evacuation safety assessment[J].Journal of Shenyang Aerospace University, 2010,2(3):27-29.
[7]何大治.基于幾何方法的人員疏散最佳路徑[J].工程圖學(xué)學(xué)報,2010,3(5):56-57.HE Da-zhi.Based on geometric method an optimal path[J].Journal of Engineering Graphics, 2010, 3 ( 5):56-57.