鄭波濤，毛曉汶，楊萍茹，許 石，郝彧露

(1.重慶城市綜合交通樞紐(集團)有限公司，重慶 400060;2.西南交通大學地球科學與環(huán)境工程學院，四川 成都 611756)

隨著世界經濟的不斷發(fā)展，在高層、超高層建筑出現的同時，多層深埋地下建筑的數量也在逐漸增加，而人員通過樓梯上行過程中需要消耗大量的體力，會嚴重影響災害時人員的疏散效率[1]。已有研究表明，電梯協助疏散有利于縮短人員疏散時間[2-3]。目前，很多學者對災害時人員電梯疏散的疏散時間進行了相關研究，如在合理調度電梯以縮短疏散時間方面，Ma等[4]通過建立電梯疏散網絡模型，總結了電梯輔助撤離的原則；Chen等[5]利用Pathfinder制定出不同疏散情景，研究了電梯疏散的策略；王晶等[6]將緊急疏散電梯調度問題描述為整數規(guī)劃問題，并給出了一種啟發(fā)式緊急疏散電梯調度的算法。此外，李釗等[7]研究了不同電梯性能、樓層高度、人員分布等因素對人員電梯疏散時間的影響。目前，在計算人員電梯疏散時間方面，一些主要電梯公司和制造商開發(fā)了一些計算人員疏散時間的仿真模型如BTS[8]、ELEVATE[9]和ELVAC[10]等模型，但對于管理者或非專業(yè)人員來說，使用這些軟件進行計算時具有一定的難度。

排隊論(Queuing Theory)是運籌學中的一個重要分支，是在各種排隊系統(tǒng)概率規(guī)律的基礎上，探究擁擠現象(排隊、等待)的科學[11]。排隊論理論已廣泛應用到生活中的各個領域，如體檢、交通等[12-15]。在災害時使用電梯和樓梯協同疏散的過程中，也存在人員需要排隊等待的問題，因此本文根據排隊論相關理論推導出電梯總疏散時間TRSET的計算公式，以為管理人員或科研工作者計算人員疏散時間提供參考。

1 電梯疏散排隊模型的建立

排隊模型是20世紀50年代初由Kendall引入的，它大大簡化了排隊系統(tǒng)的描述。具體記法如下[10]：

X/Y/Z/N

其中:X表示顧客相繼到達的時間間隔的分布；Y表示服務時間的分布；Z表示服務臺的個數；N表示等待空間數。若N為無窮大時，即為等待制時，省去N而只用X/Y/Z記一個排隊系統(tǒng)。這種分類方法是以排隊系統(tǒng)中的主要特征，即顧客相繼到達的時間間隔的分布、服務時間的分布和服務臺的個數對排隊模型進行分類的。

在排隊論中，一般令M表示負指數分布，此時隨機排隊模型稱為馬爾可夫型；D為確定型分布；Ek為k階愛爾蘭(Erlang)分布；G為一般分布；Gl為一般獨立的分布。

電梯疏散的問題，可以利用排隊論理論進行研究。電梯疏散中待疏散的人員為乘客，電梯為服務設施，Barney[16]的統(tǒng)計分析證明，乘客到達電梯服從泊松分布，即乘客到達的時間間隔X服從負指數分布(M)；電梯的服務時間指電梯運送乘客后往返一次的時間，一般假設電梯服務時間Y也服從負指數分布(M)[10]；同時假設有n部電梯，電梯排隊規(guī)則為先到先服務、后到后服務。因此，電梯疏散排隊模型可描述如下：

M/M/n

2 基于排隊論的電梯疏散時間計算公式推導

2. 1 電梯往返時間的計算

任一部電梯的往返時間tr指電梯從大廳層出發(fā)接送乘客后再次返回大廳層往返一次的運行時間。根據相關文獻[9],電梯往返時間tr包括三個部分，其表達式為

tr=tm+ts+tp

(1)

式中：tm為電梯空載情況下往返一次的運行時間(s)；ts為電梯往返一次?？繒r開關門和預啟動的時間(s)；tp為乘客進出電梯的時間(s)。

電梯往返一次可能在多個樓層停靠去接送乘客，設電梯往返一次到達的最高樓層為第h層，為了方便計算，本文將地面大廳層認定為第0層，則空載情況下電梯往返一次的運行時間tm可近似用下式表示：

tm=2ht1

(2)

式中：h為電梯往返一次運行的樓層數(層)；t1為電梯勻速運行一層樓層距離所需的時間(s)。

由于電梯往返一次的過程中必定會在大廳層?？恳淮?，設s為除大廳層外電梯往返一次停靠的次數，則電梯停靠的總次數為(s+1)次，t2為電梯每次?？繒r開關門的時間和電梯加減速延誤的時間，則電梯往返一次?？繒r開關門和預啟動的時間ts可以描述為

ts=(s+1)t2

(3)

式中：(s+1)為電梯?？康目偞螖?次)；t2為電梯每次?？繒r開關門的時間和電梯加減速延誤的時間(s)。

乘客進出電梯的時間tp可以表示為

tp=2Nit3

(4)

式中：Ni為該次往返電梯運送的人數(人)；t3為單個乘客進入電梯或走出電梯轎廂的時間(s)，假設乘客走進或走出電梯轎廂的時間相同。

因此，電梯往返時間tr的計算公式為：

tr=2ht1+(s+1)t2+2Nit3

(5)

2.2 基于乘客到達率的電梯往返一次運行時間期望值計算

設Tr、H、S、N分別為tr、h、s、Ni的期望值，為了得到Tr與乘客到達率λ之間的關系，需將變量H、S和N用乘客到達率λ表示。那么，電梯往返一次運行時間的期望值Tr為

Tr=2Ht1+(S+1)t2+2Nt3

(6)

上式中，t1、t2和t3均為常數，變量為H、S和N。由于乘客到達電梯服從泊松分布，通過極大似然估計可以得出乘客到達率λ的極大似然估計量為乘客的平均到達率。因此可以將變量H、S和N用乘客到達率λ表示，從而得到電梯往返一次運行時間的期望值Tr與乘客到達率λ之間的關系。

2.2.1 電梯往返一次運行的樓層數期望值H的計算

定義電梯相繼到達大廳層的時間間隔t為電梯往返時間T除以建筑內的電梯數量L，即t=T/L，設單部電梯各個樓層的乘客到達率為λ，在t這段時間內有N(t)=k個乘客來到某個樓層某部電梯前的概率Pk(N(t))服從參數為λ的泊松分布，即：

Pk(N(t)=k)=(λt)ke-λt/k! (k=0,1,2,…)

(7)

假設電梯在某個樓層不停留的概率為q，則停留的概率為p=1-q，當k=0時說明這個樓層沒有乘客，則電梯在該樓層不停留的概率q=e-λt，那么停留的概率p=1-q=1-e-λt。

考慮地下疏散電梯需向上運行的情形，設電梯服務的地下樓層有n層，電梯運行1層相當于電梯最低下到-1層，沒有繼續(xù)往下運行，即-1層以下樓層乘客數k均為0，則電梯到達-1層的概率為P1{H=1}=pqn-1。電梯運行樓層數為2層相當于電梯最低下到-2層，沒有繼續(xù)往下運行，則電梯到達-2層的概率為P2{H=2}=pqn-2。以此類推，則得到電梯到達的最低樓層為m的概率為Pm{H=m}=pqn-m(m=1，2,…,n)。那么電梯運行的樓層數期望值H可表示為

(8)

式中：n為樓層數；m為電梯到達的最低樓層；Pm為電梯到達的最低樓層為m的概率。

同理，可以推導出電梯在地面上運行時，電梯運行的樓層數期望值H的計算公式與上式一致。本文將地面大廳層假定為第0層，因此在描述地面上的樓層時，第m層實際為日常樓層中的第m+1層。

2.2.2 除大廳層外電梯往返一次?？看螖灯谕礢的計算

電梯往返一次運行時間內?？康拇螖悼梢杂脴菍訑党艘噪娞菰谀骋粚油Ａ舻母怕蔖t來表示，即有：

S=nPt=n(1-e-λt)

(9)

式中：n為樓層數；Pt為電梯在某一層停留的概率。

2.2.3 電梯往返一次運送的人數期望值N的計算

假設到達電梯的乘客較少，電梯往返一次運送的人數期望值N為時間間隔t內所有到達電梯的乘客，由于每層的乘客到達率為λ，那么n層的乘客到達率為nλ，即有：

N=E(n)t=nλt

(10)

式中：N為電梯往返一次運送的人數期望值，且1≤N≤R(R為電梯的荷載，N大于或等于1說明電梯運行時不能空載，N小于或等于R說明電梯不能過載運行)；E(n)為乘客在單位時間內到達電梯的人數n的期望值。

將公式(8)、(9)和(10)代入公式(6)中，可得到電梯往返一次的平均運行時間期望值Tr，即

[n(1-e-λt)+1]t2+2nλtt3

(11)

2. 3 電梯總疏散時間的計算

電梯總疏散時間為

TRSET=TA+TR+Te

式中TA為探測報警時間，保守設為60 s；TR為人員響應準備時間，保守設為120 s。

則電梯總疏散時間為

n(e-λt)mt1+[n(1-e-λt)+1]t2+2nλtt3}

(12)

2. 4 電梯疏散時間計算公式的簡化和求解

(13)

由公式(13)求出乘客平均到達率λ后，再根據不同建筑的具體情況對t1、t2、t3、n和L進行賦值，則可由公式(11)求出電梯往返時間Tr。但由于公式(11)存在指數項，直接求得相應的解具有一定的難度，因此可根據牛頓迭代法，將其轉化為求解下列線性方程的根，并通過Matlab編程求解。其線性方程為

[n(1-e-λt)+1]t2+2nλtt3

(14)

式中：t、n、λ、L為變量，m的數值隨n的變化而變化;其余為常數。

為了得到t與λ的關系式，需將n和L的值固定，則可以得到特定條件下t與λ的對應關系式。

為了得到某建筑物電梯疏散時間的通用計算公式，利用Matlab編程獲得了一組{λi，ti|i=1,2,…,n}的測量數據，將這組數據進行多項式擬合，可得到簡化后的λ與t的擬合關系曲線和關系式，進而得到電梯往返時間簡化后的表達式。在不同的建筑中，只要改變源程序中的t1、t2、t3、n和L的值，并對得到的{λi,ti}測量數據進行多項式擬合，即可得到λ與t的擬合關系曲線和關系式，進而得到適用于該建筑電梯總疏散時間TRSET的計算公式。

3 實例應用與分析

本文以重慶沙坪壩地下大型綜合交通樞紐為實例，利用上述得到的電梯疏散時間計算公式對該建筑電梯總疏散時間進行了計算與分析，以驗證該計算公式的有效性。

該地下建筑共負七層(n=7)，共有34個載客電梯(L=34)；每層高5 m，電梯運行額定速度v取2.5 m/s，在相鄰的兩個樓層運行所需的時間t1取3 s；電梯為中分自動門，每次開關門的時間為4 s，考慮到電梯每次?？繒r加減速造成的時間延誤，電梯?？恳淮蔚臅r間t2取5 s；經統(tǒng)計得出單個乘客進入或走出電梯轎廂的平均時間t3為2 s，具體參數詳見表3。

表1 重慶沙坪壩地下大型綜合交通樞紐的參數設置

本文利用Matlab編程，將建筑物參數代入到Matlab程序中，可得到簡化后的λ與t的擬合關系方程為

t=0.000 415λ3-0.040 176λ2+1.215 564λ-9.637 067

根據公式(12)可計算得到，沙坪壩地下綜合交通樞紐電梯總疏散時間TRSET為

(15)

本文選取某種工況來進行模擬驗證，該工況為：負六到負一層的人員走樓梯，負七層的人員坐電梯。由模擬結果可知，該建筑電梯總疏散時間為1 740.5 s，約為29 min。負七層中乘坐電梯的總人數為3 659人，計算得乘客平均到達率λ為126 人/min，將λ的數值代入公式(15)求解，可計算得到電梯總疏散時間為27.43 min。該計算結果與模擬結果29 min相差 1.57 min，相對誤差為5.4%，證明電梯疏散時間的計算公式是有效的。

4 結 論

(2) 以重慶沙坪壩地下大型綜合交通樞紐為例，通過模擬得到的電梯疏散時間與通過本文公式計算得到的電梯疏散時間之間的相對誤差為5.4%，誤差較小， 證明了本文提出的電梯疏散時間計算公式的有效性。利用本文的計算公式可以在較短時間內計算出建筑物電梯疏散時間，可為緊急疏散情況下管理人員的決策提供幫助。