• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    考慮在能見度受限下行人跟隨行為特性的建模與模擬*

    2019-12-24 08:21:18楊燦陳群陳璐
    物理學(xué)報(bào) 2019年24期
    關(guān)鍵詞:元胞行人視野

    楊燦 陳群 陳璐

    (中南大學(xué)交通運(yùn)輸工程學(xué)院, 長(zhǎng)沙 410075)

    研究了行人在能見度受限情況下的疏散行為, 考慮行人對(duì)環(huán)境的熟悉程度, 將行人分為熟悉環(huán)境人群和不熟悉環(huán)境人群.對(duì)于房間內(nèi)熟悉環(huán)境的行人, 改進(jìn)勢(shì)函數(shù)元胞自動(dòng)機(jī)模型來(lái)模擬其疏散行為.對(duì)于不熟悉環(huán)境人群, 分析其在視野范圍內(nèi)的跟隨行為, 并制定了不同跟隨行為策略, 來(lái)研究其跟隨行為特性.仿真模擬了房間內(nèi)熟悉環(huán)境人群的人數(shù)占比、房間內(nèi)的視野半徑大小以及行人密度等參數(shù), 研究其對(duì)行人疏散的影響, 比較不熟悉環(huán)境人群采取的跟隨策略的優(yōu)劣.發(fā)現(xiàn)疏散時(shí)間的大小與房間內(nèi)視野半徑的大小和房間內(nèi)熟悉環(huán)境者密度的大小有關(guān).其次, 跟隨策略的有效性與視野半徑的大小和熟悉環(huán)境者密度有關(guān).而且在單一策略環(huán)境下, 有著同樣的規(guī)律.這些發(fā)現(xiàn)能對(duì)大型公共場(chǎng)所如超市、體育館的應(yīng)急疏散情況提供一些啟示, 有助于在視野受限情況下制定一些有效的指導(dǎo)策略.

    1 引 言

    隨著社會(huì)經(jīng)濟(jì)的不斷發(fā)展, 商場(chǎng)、劇院、體育館、博物館、車站、機(jī)場(chǎng)、綜合交通樞紐等大型公眾設(shè)施也在不斷地發(fā)展和完善, 這極大地豐富和方便了居民生活, 行人的出行頻率也大大提高.但是行人在建筑物內(nèi)易受突發(fā)事件(如擁擠、踩踏和火災(zāi)等情況) 的影響, 引起重大的群體事故.因此, 建筑物內(nèi)的行人疏散一直是交通研究者研究的熱點(diǎn)問(wèn)題之一.

    目前, 科學(xué)家們先后提出了一系列各具特色的計(jì)算機(jī)仿真模型來(lái)研究建筑物內(nèi)的行人疏散行為,越來(lái)越多的微觀仿真模型被應(yīng)用在建筑物疏散中.其中, 最為廣泛應(yīng)用的微觀模型包括:社會(huì)力模型[1-5]、格子氣模型[6-8]和元胞自動(dòng)機(jī)模型[9-17].Helbing等[1-4]提出的社會(huì)力模型是一個(gè)連續(xù)的仿真模型, 它研究了恐慌情形下行人的疏散行為, 發(fā)現(xiàn)了行人的自組織現(xiàn)象、出口處行人的聚集與成拱現(xiàn)象和快即是慢效應(yīng).Burstedde等[12]提出的元胞自動(dòng)機(jī)場(chǎng)域模型也能較好地模擬行人動(dòng)態(tài)的聚集效應(yīng)和自組織現(xiàn)象.目前, 大多數(shù)模型都是在此基礎(chǔ)上進(jìn)行改進(jìn)研究, 來(lái)模擬行人的疏散過(guò)程.

    目前, 有部分研究者開始研究能見度受限情況下行人的疏散行為, 即考慮房間內(nèi)部能見度的大小和行人對(duì)房間的熟悉程度.Hou等[18]和Ma等[19,20]建立社會(huì)力模型來(lái)研究有限能見度下行人的領(lǐng)導(dǎo)作用, 認(rèn)為有效的領(lǐng)導(dǎo)力與視野半徑范圍的大小息息相關(guān), 在有限能見度下行人的領(lǐng)導(dǎo)對(duì)疏散影響是雙重的.Guo等[21]通過(guò)真人疏散實(shí)驗(yàn)建立了零能見度下的元胞自動(dòng)機(jī)模型, 仿真表明在零能見度條件下的行人速度約為正常能見度情況下行人速度的一半, 并且認(rèn)為疏散演習(xí)有助于提高零能見度條件下行人疏散的效率.Frank和Dorso[22]以及Wang和Cao[23]都考慮了在視野受限的房間內(nèi), 行人疏散過(guò)程的不同行為(包括放牧、沿墻壁行走等).Li等[24]建立了考慮視野受限的元胞自動(dòng)機(jī)模型來(lái)分析視野半徑對(duì)疏散行人的影響, 認(rèn)為疏散時(shí)間依賴于房間內(nèi)的初始密度和視野半徑.Li和Han[25]研究疏散房間內(nèi)行人的環(huán)境熟悉程度和積極性行為來(lái)分析對(duì)疏散時(shí)間的影響, 認(rèn)為疏散時(shí)間隨著疏散人員對(duì)疏散路線熟悉程度的增加而呈指數(shù)下降.

    而研究者們對(duì)疏散過(guò)程中行人的跟隨行為研究較少, 大多數(shù)研究都是關(guān)于領(lǐng)導(dǎo)者和導(dǎo)游對(duì)行人的領(lǐng)導(dǎo)作用.在社會(huì)力模型中, Hou等[18]建立了具有領(lǐng)導(dǎo)效應(yīng)的社會(huì)力模型來(lái)研究領(lǐng)導(dǎo)者對(duì)普通行人的領(lǐng)導(dǎo)作用, 研究了疏散領(lǐng)導(dǎo)者的數(shù)量和位置對(duì)能見度有限的疏散動(dòng)力學(xué)的影響, 認(rèn)為更多的領(lǐng)導(dǎo)者促使單出口房間的疏散更快, 但是一個(gè)或兩個(gè)受過(guò)訓(xùn)練的領(lǐng)導(dǎo)者也足以加速動(dòng)態(tài)疏散.Ma等[19,20]改進(jìn)社會(huì)力模型來(lái)研究人群疏散過(guò)程中的領(lǐng)導(dǎo)行為, 認(rèn)為領(lǐng)導(dǎo)者和導(dǎo)游并不是越多越好, 一小部分領(lǐng)導(dǎo)者已經(jīng)足以有效地引導(dǎo)整個(gè)人群.而在元胞自動(dòng)機(jī)模型中, 典型的樓層元胞自動(dòng)機(jī)模型[12]的主要思想是基于仿生學(xué)來(lái)捕捉行人的放牧行為, 以及在疏散過(guò)程中尋找最短路徑的行為.在雙向行人流的研究過(guò)程中, Zhang[26]改進(jìn)勢(shì)函數(shù)元胞自動(dòng)機(jī)模型分析了雙向行人流的跟隨行為、推動(dòng)行為和避讓行為, 仿真結(jié)果顯示改進(jìn)模型更能表現(xiàn)行人合理性和多樣性, 而且更能反映集體行為特征.Guo等[27]建立了勢(shì)函數(shù)模型來(lái)分析行人跟隨行為和避碰行為的元胞自動(dòng)機(jī)模型, 對(duì)雙向行人流進(jìn)行了仿真研究, 認(rèn)為對(duì)于渠道的形成, 避碰行為更優(yōu)于跟隨行為.在房間內(nèi)行人疏散情況中, Yuan和 Tan[28]改進(jìn)了元胞自動(dòng)機(jī)模型, 模擬了行人“隨流流動(dòng)”的現(xiàn)象, 認(rèn)為引導(dǎo)者的能見度和數(shù)量能顯著影響疏散時(shí)間.而Lu等[29]改進(jìn)場(chǎng)域模型來(lái)研究行人小團(tuán)體的領(lǐng)導(dǎo)者-跟隨者行為規(guī)則, 能夠很好地捕獲疏散過(guò)程中群體移動(dòng)的特征, 并成功預(yù)測(cè)總疏散時(shí)間.Wang等[30]建立了擴(kuò)展的場(chǎng)域模型來(lái)研究疏散助手的位置和數(shù)量對(duì)行人疏散的影響, 建立了信息傳遞場(chǎng)表示對(duì)行人的領(lǐng)導(dǎo)作用, 認(rèn)為在疏散出口的疏散助手更能有效地領(lǐng)導(dǎo)行人的疏散, 而過(guò)多的疏散助手并不能提高疏散效率.Fang等[31]建立了基于代理的跟隨者模型, 提出了一個(gè)領(lǐng)導(dǎo)者-跟隨者模型來(lái)解釋社會(huì)集體行為, 并且建模了逆流、排隊(duì)行走和集體移動(dòng)三種場(chǎng)景.

    2012年, Zhang等[16]建立了勢(shì)函數(shù)場(chǎng)元胞自動(dòng)機(jī)模型, 該模型考慮了行人的行走時(shí)間和周圍行人引起的不適成本, 以便行人移動(dòng)到空的相鄰元胞.移動(dòng)產(chǎn)生的費(fèi)用將簡(jiǎn)單地測(cè)量當(dāng)前元胞和目的地之間的距離, 類似于場(chǎng)域模型中的靜態(tài)場(chǎng)[12].同時(shí)又考慮行人的行程中的負(fù)舒適度(周圍密度的增函數(shù)), 而負(fù)舒適度的作用類似于場(chǎng)域模型中的動(dòng)態(tài)場(chǎng)[12].2014年, Jian 等[17]提出了一種感知的勢(shì)函數(shù)元胞自動(dòng)機(jī)模型, 該模型考慮聚集力場(chǎng)和可視度的影響, 得到考慮聚集力場(chǎng)的預(yù)測(cè)費(fèi)用勢(shì)函數(shù)場(chǎng)模型, 并合理再現(xiàn)了有瓶頸或轉(zhuǎn)角區(qū)域內(nèi)的行人流現(xiàn)象.本文考慮行人對(duì)環(huán)境的熟悉程度, 將行人分為熟悉環(huán)境人群和不熟悉環(huán)境人群, 由于勢(shì)函數(shù)元胞自動(dòng)機(jī)模型[16]能模擬行人沿著疏散時(shí)間最短線路進(jìn)行疏散, 能更切合實(shí)際地模擬熟悉環(huán)境人群的疏散行為, 因此本文采用改進(jìn)勢(shì)函數(shù)元胞自動(dòng)機(jī)模型來(lái)模擬熟悉環(huán)境人群的疏散行為.而對(duì)于不熟悉環(huán)境人群來(lái)說(shuō), 考慮其在視野半徑內(nèi)的跟隨策略,元胞自動(dòng)機(jī)能較好地制定不熟悉環(huán)境者的跟隨規(guī)則, 這些在周金旺等[13]、永貴等[14]、岳昊等[15]的研究中也有所體現(xiàn).而且, 本文對(duì)行人高密度下的疏散進(jìn)行了仿真, 元胞自動(dòng)機(jī)模型能再現(xiàn)各類復(fù)雜的疏散現(xiàn)象的特點(diǎn).綜合考慮, 本文采取元胞自動(dòng)機(jī)模型進(jìn)行建模, 來(lái)模擬行人跟隨疏散行為.

    在大多數(shù)與上述研究相關(guān)的行人疏散模擬中,研究者們研究了房間內(nèi)能見度受限或視野較差的情況, 但忽視了行人對(duì)環(huán)境的熟悉程度這一重要的因素.因此本文考慮在房間內(nèi)能見度受限情況下,將行人分為熟悉環(huán)境人群和不熟悉環(huán)境人群綜合考慮其疏散行為, 這比以往的研究更貼近實(shí)際情況.而且, 上述文獻(xiàn)對(duì)疏散過(guò)程中行人的跟隨行為考慮較少, 大多數(shù)研究都是研究從眾行為、領(lǐng)導(dǎo)者和導(dǎo)游對(duì)行人的領(lǐng)導(dǎo)作用, 而本文考慮了在能見度受限情況下, 不熟悉環(huán)境人群的跟隨行為, 不熟悉環(huán)境人群采取不同跟隨策略, 依賴視野范圍內(nèi)的其他行人進(jìn)行移動(dòng), 從而達(dá)到安全疏散的目的.因此,本文相比其他文獻(xiàn), 分析了不熟悉環(huán)境人群的跟隨策略, 并充分考慮了其疏散過(guò)程中的跟隨行為, 更符合實(shí)際疏散情況.

    因此, 與以往的研究相比, 本文研究的行人疏散行為同時(shí)考慮了能見度受限和行人對(duì)環(huán)境熟悉度的不同情況.改進(jìn)費(fèi)用勢(shì)場(chǎng)元胞自動(dòng)機(jī)模型來(lái)模擬熟悉環(huán)境者人群的疏散, 同時(shí)考慮不熟悉環(huán)境者的跟隨行為, 制定其跟隨策略和規(guī)則, 對(duì)其跟隨行為進(jìn)行仿真研究, 并分析熟悉環(huán)境人群的比例數(shù)量、房間內(nèi)行人密度和房間內(nèi)不熟悉環(huán)境人群的可視半徑等參數(shù)對(duì)行人疏散時(shí)間以及行人疏散演化過(guò)程的影響.

    2 模 型

    疏散模型建立在房間內(nèi)部為 (L+ 2) × (L+ 2)的二維離散元胞網(wǎng)格系統(tǒng)內(nèi).黑色邊框表示障礙物墻壁, 把行人疏散空間的移動(dòng)區(qū)域分割為L(zhǎng)×L個(gè)大小相等的離散元胞空格.每個(gè)元胞網(wǎng)格可以是空的或被行人占用, 并且其大小可被近似視為0.4 m × 0.4 m[12], 而W表示房間出口寬度.每個(gè)行人也只能占據(jù)一個(gè)元胞空格, 行人的仿真過(guò)程也被離散化為相等的時(shí)間步長(zhǎng), 在每個(gè)離散的時(shí)間步長(zhǎng)內(nèi), 行人僅僅能移動(dòng)一個(gè)元胞的位置.每個(gè)時(shí)間步長(zhǎng)為0.3 s, 行人移動(dòng)速度為1時(shí)間步一個(gè)元胞,則行人移動(dòng)速度為1.33 m/s.在每個(gè)時(shí)間步長(zhǎng)內(nèi),行人或停止等待, 或以最大的速度移動(dòng), 且行人不能穿越和跨越圍墻.如圖1所示, 整個(gè)系統(tǒng)分為兩個(gè)區(qū)域:一個(gè)表示可見出口區(qū)域A、一個(gè)表示不可見出口區(qū)域B.白色區(qū)域?yàn)樾腥说淖咝袇^(qū)域, 并且假設(shè)每個(gè)行人的視野半徑相等, 為R, 如圖1所示,紅色實(shí)心圓表示位于視野半徑為R中心的行人.

    本文考慮行人對(duì)房間的熟悉程度, 將行人分為熟悉環(huán)境人群和不熟悉環(huán)境人群, 假定熟悉環(huán)境者能清楚地知道安全出口的位置, 即使是在視野較差的情況下, 也能知道最短路線和時(shí)間最短路線.因此本文采用勢(shì)函數(shù)元胞自動(dòng)機(jī)模型來(lái)進(jìn)行建模, 而不是考慮距離最短的場(chǎng)域模型進(jìn)行建模.

    2.1 改進(jìn)的勢(shì)函數(shù)場(chǎng)

    勢(shì)函數(shù)場(chǎng)元胞自動(dòng)機(jī)模型由Zhang等[16]定義, 表示為當(dāng)前小區(qū)移動(dòng)到目的地的最低成本.對(duì)于勢(shì)函數(shù)場(chǎng)元胞自動(dòng)機(jī)模型來(lái)說(shuō), 起關(guān)鍵性作用的因素是局部費(fèi)用分布函數(shù).引入局部費(fèi)用分布函數(shù), 使得勢(shì)函數(shù)場(chǎng)成為靜態(tài)場(chǎng)和動(dòng)態(tài)場(chǎng)的結(jié)合體,模型可以被應(yīng)用于模擬正常情況下行人運(yùn)動(dòng)達(dá)到系統(tǒng)費(fèi)用最優(yōu)時(shí)的運(yùn)動(dòng)狀態(tài).

    費(fèi)用分布函數(shù)表示:一個(gè)是由行人的步行速度和從當(dāng)前位置到目標(biāo)位置的距離確定的行進(jìn)時(shí)間的成本; 另一個(gè)是取決于周圍密度的不舒服的成本.對(duì)于費(fèi)用函數(shù)來(lái)說(shuō), 在固定時(shí)刻t, 在元胞相應(yīng)的位置(x,y)處的費(fèi)用分布函數(shù)為[16]

    其中等式右邊第一部分 1 /vmax代表行人的步行速度和從當(dāng)前位置到目標(biāo)出口的距離確定的移動(dòng)時(shí)間的成本, 等式右邊第二部分g(ρ(x,y,t)) 代表與密度有關(guān)的行人的不適成本.而假設(shè)[16]

    從當(dāng)前位置 (x,y) 走向安全出口位置 (x0,y0) 的總費(fèi)用為[16]

    假設(shè)積分與路徑無(wú)關(guān), 即[16]:

    上述的模型為Zhang等[16]建立的勢(shì)函數(shù)場(chǎng)域模型, 而本文認(rèn)為熟悉環(huán)境者在其視野范圍內(nèi), 在其移動(dòng)方向上存在和不熟悉環(huán)境者移動(dòng)方向的沖突, 因此考慮引入視野范圍內(nèi)不熟悉環(huán)境者的干擾函數(shù), 來(lái)進(jìn)一步表示影響熟悉環(huán)境者的不舒適度,而視野范圍為當(dāng)前行人行駛方向的半圓區(qū)域, 如圖2所示, 而干擾函數(shù)定義如下.

    圖2 其他行人對(duì)熟悉環(huán)境者行人影響圖Fig.2.Influence of the moving direction for surrounding pedestrian on the desired moving direction for the center pedestrian.

    考慮當(dāng)前行人行駛方向的左右90°扇形區(qū)間,在此視野范圍內(nèi), 影響該行人運(yùn)動(dòng)的干擾計(jì)算如下.

    對(duì)于左側(cè)90°區(qū)域, 計(jì)算此區(qū)域其他行人與當(dāng)前運(yùn)動(dòng)行人的沖突,θ為其他行人移動(dòng)方向和當(dāng)前行人移動(dòng)方向的逆時(shí)針夾角,

    如果θ為[180°, 360°), 則Nli=0 , 如果θ為(0, 180°),則Nli=1 .

    對(duì)于右側(cè)90°區(qū)域, 計(jì)算此區(qū)域其他行人與當(dāng)前運(yùn)動(dòng)行人的沖突,θ為其他行人移動(dòng)方向和當(dāng)前行人移動(dòng)方向的逆時(shí)針夾角,

    如果θ為(0, 180°)則Nri=0 , 如果θ為[180°, 360°),則Nri=1 .

    因此, 其干擾函數(shù)為

    同時(shí), 考慮原勢(shì)函數(shù)模型的不舒適成本, 也就是視野范圍內(nèi)行人密度引起的不舒適成本, 視野范圍內(nèi)行人密度p為視野范圍內(nèi)總?cè)藬?shù)Na與視野面積之比r2,

    參考文獻(xiàn) [16], 視野半徑為 1 m 時(shí),r= 3, 半徑大于等于 2 m 時(shí), 認(rèn)為r= 5.

    因此本文將不適的原始成本分為兩部分:一部分與視野范圍內(nèi)行人密度有關(guān), 另一個(gè)是由視野范圍內(nèi)其他行人運(yùn)動(dòng)變化引起的.因此, 將成本函數(shù)重新定義為

    參考文獻(xiàn)[16], 這里選擇相應(yīng)的系數(shù)為a=0.075,b= 0.025,α=β=2 .

    2.2 不熟悉環(huán)境者走行策略及規(guī)則

    在房間內(nèi)能見度受限時(shí), 不熟悉環(huán)境者的行人有多種疏散策略, 本文考慮以下五種跟隨策略作為不熟悉環(huán)境者的移動(dòng)策略, 不熟悉環(huán)境人群策略如圖3所示.

    策略S0如果不熟悉環(huán)境者處在可見安全出口區(qū)域, 意味著此類行人能夠清晰地知道出口位置, 此區(qū)域內(nèi)的行人遵守靜態(tài)場(chǎng)域模型[12]的行走規(guī)則.房間內(nèi)區(qū)域分布如圖4所示.

    圖3 不熟悉環(huán)境人群策略圖Fig.3.Diagram of uninformed pedestrian strategies.

    圖4 房間內(nèi)區(qū)域分布圖Fig.4.Illustration of different areas.

    在可見安全出口區(qū)域內(nèi)構(gòu)建靜態(tài)場(chǎng)模型, 模型采用元胞位置距安全出口的最短距離表示房間內(nèi)不同位置對(duì)行人的吸引程度.距安全出口內(nèi)最近元胞的距離值作為元胞位置距安全出口的最短距離.元胞位置距安全出口的最短距離由下式計(jì)算[12]:

    式中Sij為元胞 (x,y) 周邊鄰域距安全出口的最短距離,Sxy為元胞 (x,y) 距離安全出口的最短距離;(xi,yj)為元胞 (x,y) 在疏散系統(tǒng)中鄰域的坐標(biāo);(xn,yn)為出口內(nèi)第n個(gè)元胞在疏散系統(tǒng)中的坐標(biāo).

    如果不熟悉環(huán)境者處在不可見安全出口區(qū)域,采取四種跟隨策略進(jìn)行移動(dòng).

    策略S1視野范圍內(nèi)存在其他行人, 會(huì)等概率地選擇視野范圍內(nèi)一個(gè)行人進(jìn)行跟隨, 跟隨其方向進(jìn)行移動(dòng).在跟隨過(guò)程中, 若被跟隨者消失在自己視野范圍內(nèi), 則會(huì)重新選擇周邊行人繼續(xù)跟隨.如圖5所示, 紅圈代表跟隨的行人.其跟隨規(guī)則為

    其中 (x,y) 表示當(dāng)前元胞坐標(biāo), (x1,y1) 表示被跟隨元胞坐標(biāo), (xi,yj) 表示跟隨者周邊鄰域的八個(gè)位置的坐標(biāo),Sm表示鄰居元胞距被跟隨者的最短距離.

    策略S2行人考慮視野范圍內(nèi)的8個(gè)區(qū)域,向視野范圍內(nèi)行人最多的區(qū)域進(jìn)行移動(dòng), 表示從眾跟隨行為, 而跟隨方向?yàn)镸oore[12]鄰域的八個(gè)方向之一, 如圖6所示, 行人對(duì)八個(gè)方向區(qū)域的人數(shù)進(jìn)行判斷, 每?jī)蓧K扇形區(qū)域代表一個(gè)方向.其跟隨規(guī)則如下:

    其中h表示下一步的方向,Nq表示視野范圍內(nèi)區(qū)域的人數(shù).

    策略S3行人根據(jù)視野范圍內(nèi)大多數(shù)行人的移動(dòng)方向來(lái)確定自己的移動(dòng)方向, 行人跟隨多數(shù)人的移動(dòng)方向進(jìn)行移動(dòng).一個(gè)時(shí)間步計(jì)算一次, 即跟隨上一時(shí)間步視野范圍內(nèi)行人最多進(jìn)行移動(dòng)的方向, 采用 Moore[12]鄰域, 走行方向如圖6 所示, 表示跟隨眾多行人移動(dòng)方向行走.走行規(guī)則如下:

    其中h表示下一步的方向,Ne表示視野范圍內(nèi)八個(gè)移動(dòng)方向的人數(shù).

    策略S4行人堅(jiān)定自己的方向, 選擇不跟隨他人, 按照最初設(shè)定的方向一直移動(dòng), 到達(dá)可見墻壁區(qū)域后, 往墻壁方向移動(dòng), 隨機(jī)選擇順時(shí)針或逆時(shí)針方向, 跟隨墻壁行走, 直到找到出口.

    圖5 策略 S1 跟隨示意圖Fig.5.Following strategy S1.

    圖6 策略 S2, S3 跟隨示意圖Fig.6.Following strategy S2 and S3.

    對(duì)于選擇S1, S2和S3策略的行人, 若開始疏散時(shí)其視野范圍內(nèi)沒有其他行人, 則按照S4策略進(jìn)行移動(dòng), 直到視野范圍內(nèi)存在其他行人, 后采取跟隨策略.假定一旦選擇跟隨策略就不能改變, 一直以此策略運(yùn)動(dòng)下去, 直到走到出口或者停留在某個(gè)區(qū)域, 對(duì)于停留在某個(gè)區(qū)域的行人來(lái)說(shuō), 在停留t步之后, 按照S4策略所定規(guī)則進(jìn)行移動(dòng).

    2.3 移動(dòng)概率

    本文采用Moore[12]鄰域, 被占用的小區(qū)有8個(gè)相鄰小區(qū), 各個(gè)小區(qū)的移動(dòng)概率Pi,j計(jì)算過(guò)程如下, 行人鄰域及其移動(dòng)概率示意圖如圖7所示.

    圖7 行人鄰域及概率示意圖 (a)當(dāng)前元胞及8個(gè)鄰居元胞圖; (b)當(dāng)前元胞移動(dòng)概率圖Fig.7.Neighboring cell and moving probabilities of pedestrian:(a) An occupied cell and its eight neighboring cells;(b) the corresponding probabilities for occupation in the next update step.

    對(duì)于熟悉環(huán)境人群來(lái)說(shuō), 其移動(dòng)概率如下.

    對(duì)于元胞(i,j)∈S0=(i,j)|(i,j)∈S0,0為空元胞,計(jì)算差商其中,di,j為元 胞 (i,j) 和 (0,0) 之間的距離; 構(gòu)造集合Sm={(i,j)|δ?i,j(0)=δ?min(0)}?S0, 其 中 ,δ?min(0)=min(i,j)?S0δ?i,j(0); 因此, 定義可移動(dòng)概率為:

    否則,

    其中 |Sm| 代表集合Sm的元素個(gè)數(shù).

    對(duì)于不熟悉環(huán)境人群來(lái)說(shuō), 位于可見出口區(qū)域, 采取S0策略靜態(tài)場(chǎng)域進(jìn)行移動(dòng), 而對(duì)于靜態(tài)場(chǎng)移動(dòng)概率如下.

    對(duì){于靜態(tài)場(chǎng)移動(dòng)概率來(lái)說(shuō)},對(duì)于元胞(i,j)∈S0=(i,j)|(i,j)∈S0,0為空元胞, 計(jì)算差商?S=Sxy-Sij, 構(gòu)造集合Sn={(i,j)|max?S}?S0; 因此, 定義可移動(dòng)概率為:

    當(dāng)Sn?, 且max?S>0時(shí),

    否則,

    其中 |Sn| 代表集合Sn的元素個(gè)數(shù).

    而對(duì)于不熟悉環(huán)境人群采取其他四種跟隨策略的行人, 按照上節(jié)所述的規(guī)則進(jìn)行行走, 所以此處不再定義其移動(dòng)概率.

    2.4 模型仿真步驟

    仿真采用并行更新的規(guī)則來(lái)更新行人位置, 模型步驟描述如下(圖8所示).

    圖8 仿真步驟Fig.8.Simulation step diagram.

    步驟1在房間內(nèi)隨意分配密度為ρ的行人,其中密度ρ為總?cè)藬?shù)N除以內(nèi)部房間格子數(shù)L×L.

    步驟2設(shè)置熟悉環(huán)境者比例Ks, 即熟悉環(huán)境行人為N×Ks若干,N-N×Ks為不熟悉環(huán)境行人, 不熟悉環(huán)境行人采取四種跟隨策略進(jìn)行移動(dòng), 四種跟隨行人的數(shù)量為 (N-N×Ks)/4 .熟悉環(huán)境者計(jì)算勢(shì)函數(shù)場(chǎng)進(jìn)行移動(dòng), 其他行人按照疏散策略進(jìn)行移動(dòng).

    步驟3如果存在多個(gè)行人競(jìng)爭(zhēng)一個(gè)空元胞時(shí), 隨機(jī)選擇一個(gè)行人進(jìn)入該元胞, 其他行人退回自己原位置.

    步驟4當(dāng)疏散行人移動(dòng)到安全出口內(nèi)時(shí), 在下一時(shí)間步, 行人將移出系統(tǒng).當(dāng)所有疏散行人都移出系統(tǒng)后, 仿真過(guò)程結(jié)束.

    在并行更新機(jī)理的演化規(guī)則中, 模型在每個(gè)時(shí)間步長(zhǎng)內(nèi)需要并行地把行人的移動(dòng)規(guī)則應(yīng)用到每一個(gè)人, 并做相應(yīng)的數(shù)據(jù)統(tǒng)計(jì)工作.

    3 仿真及結(jié)果分析

    在仿真研究中, 疏散行人流的密度ρ定義為系統(tǒng)內(nèi)的行人數(shù)量N與元胞個(gè)數(shù)L×L之比,L=25個(gè)元胞.熟悉環(huán)境者密度定義為熟悉環(huán)境者的數(shù)量與初始時(shí)刻房間行人數(shù)量之比, 行人的疏散時(shí)間T定義為系統(tǒng)內(nèi)所有行人離開房間所需要的時(shí)間步數(shù), 為減少仿真初始狀態(tài)對(duì)各項(xiàng)統(tǒng)計(jì)指標(biāo)的影響, 每項(xiàng)統(tǒng)計(jì)指標(biāo)分別取20次運(yùn)行指標(biāo)的平均值.在仿真的初始時(shí)刻, 所有疏散行人隨機(jī)地分布在房間內(nèi), 安全出口內(nèi)沒有設(shè)置行人.待行人全部離開房間后, 仿真過(guò)程終止.對(duì)于停留在某個(gè)區(qū)域的不熟悉環(huán)境者, 認(rèn)為其停留t步之后, 按照策略S4進(jìn)行移動(dòng).由于停留步數(shù)t不能小于出口處正常排隊(duì)時(shí)行人停留的最大步數(shù), 通過(guò)仿真不同密度下疏散情況, 在密度 0.1 時(shí),t= 30, 密度 0.2 時(shí),t= 50,密度為 0.3 時(shí),t= 80, 密度為 0.4 時(shí),t= 100.仿真相關(guān)參數(shù)說(shuō)明如表1所列.

    3.1 視野半徑和熟悉環(huán)境者密度對(duì)疏散時(shí)間的影響

    疏散房間定義為 10 m × 10 m, 即為 25 ×25個(gè)元胞, 在本節(jié)仿真中, 考慮視野半徑R以及熟悉環(huán)境者密度Ks對(duì)疏散時(shí)間的影響.由于房間為10 m × 10 m, 考慮視野半徑為 1—5 m, 熟悉環(huán)境者密度Ks為0—0.9.如果超過(guò) 400個(gè)時(shí)間步, 房間內(nèi)行人還未能全部疏散, 則認(rèn)為最大疏散時(shí)間為400步.仿真結(jié)果如圖9所示.

    圖9(a)顯示了視野半徑為1 m時(shí), 在不同密度下, 疏散時(shí)間隨熟悉環(huán)境者密度的增加而減少.在人流密度較大時(shí), 疏散時(shí)間曲線先緩慢地下降,然后迅速地下降.而在人流密度為0.1時(shí), 疏散時(shí)間曲線一直平緩地下降.圖9(b)的變化趨勢(shì)和圖9(a)類似, 不同的是在熟悉環(huán)境者密度達(dá)到0.6時(shí), 疏散時(shí)間開始保持緩慢地減少.

    圖9(c)—(e)三幅圖的疏散時(shí)間變化趨勢(shì)相似, 隨著熟悉環(huán)境者密度的增加先緩慢地減少, 然后迅速地減少, 最后保持平緩.不同的是, 圖9(c)和圖9(d)中, 熟悉環(huán)境者密度達(dá)到0.4時(shí), 曲線保持平緩.而在圖9(e)中, 熟悉環(huán)境者密度達(dá)到0.2時(shí), 曲線開始保持平緩.

    因此, 視野半徑一定時(shí), 疏散時(shí)間隨著熟悉環(huán)境者密度的增加而持續(xù)減少.意味著房間內(nèi)熟悉環(huán)境者人數(shù)的增多, 更多的人準(zhǔn)確知道出口位置, 從而引導(dǎo)不熟悉環(huán)境者向出口方向移動(dòng), 這在不同視野半徑下有著相同的規(guī)律.而對(duì)于視野半徑R≥2 m時(shí), 疏散時(shí)間在熟悉環(huán)境者密度為0.4之后保持穩(wěn)定.因此, 考慮熟悉環(huán)境密度為 0—0.5 時(shí), 分析視野半徑對(duì)疏散時(shí)間的影響.

    從圖10(a)—(d)四幅圖可以看出, 相同人流密度下疏散時(shí)間隨視野半徑的變化保持相同的趨勢(shì), 隨著視野半徑的增加, 疏散時(shí)間保持緩慢的下降.而圖10(e), 在熟悉環(huán)境者密度為 0.4 時(shí), 當(dāng)視野半徑為1—3 m時(shí), 疏散時(shí)間下降的速率較大,視野半徑為4—5 m時(shí), 疏散時(shí)間基本保持不變.圖10(f)和圖10(e)的變化趨勢(shì)一致, 不同的是, 疏散時(shí)間曲線在視野半徑2 m處保持平緩.

    表1 參數(shù)說(shuō)明Table 1.Parameters description.

    圖9 環(huán)境熟悉度與疏散時(shí)間的關(guān)系圖 (a) R = 1 m; (b) R = 2 m; (c) R = 3 m; (d) R = 4 m; (e) R = 5 mFig.9.Relationship between environmental familiarity and evacuation time:(a) R = 1 m; (b) R = 2 m; (c) R = 3 m; (d) R = 4 m;(e) R = 5 m.

    因此, 熟悉環(huán)境者密度一定時(shí), 疏散時(shí)間隨著行人視野半徑的增加而持續(xù)減少.意味著房間內(nèi)能見度越好, 更多的行人能準(zhǔn)確知道出口位置, 向出口方向移動(dòng), 不熟悉環(huán)境者的跟隨策略也越有效,從而導(dǎo)致疏散時(shí)間減少.

    行人密度為0.3, 熟悉環(huán)境者密度為0.4, 視野半徑為3 m時(shí), 行人疏散流演化過(guò)程如圖11所示.

    3.2 視野半徑對(duì)四類跟隨策略疏散時(shí)間的影響

    由3.1節(jié)可知, 疏散時(shí)間在熟悉環(huán)境者密度Ks=0.4后開始趨于平緩, 因此本節(jié)考慮在熟悉環(huán)境者密度為0.4, 房間內(nèi)行人密度為0.1—0.4, 各疏散策略的疏散時(shí)間隨著視野半徑的變化情況.同時(shí), 將策略 S1, S2, S3 和 S4 的疏散時(shí)間分別用TS1,TS2,TS3和TS4表示.

    圖10 視野半徑與疏散時(shí)間的關(guān)系圖 (a) Ks = 0; (b) Ks = 0.1; (c) Ks = 0.2; (d) Ks = 0.3; (e) Ks = 0.4; (f) Ks = 0.5Fig.10.Relationship between different visibilities and evacuation time:(a) Ks = 0; (b) Ks = 0.1; (c) Ks = 0.2; (d) Ks = 0.3;(e) Ks = 0.4; (f) Ks = 0.5.

    圖11 行人疏散流演化過(guò)程 (ρ =0.3,Ks=0.4,R=3) (a) T = 0; (b) T = 20; (c) T = 40; (d) T = 60; (e) T = 90; (f) T = 110Fig.11.Evacuation evolution (ρ =0.3,Ks=0.4,R=3):(a) T = 0; (b) T = 20; (c) T = 40; (d) T = 60; (e) T = 90; (f) T = 110.

    圖12 不同密度下各策略視野半徑與疏散時(shí)間的關(guān)系 (a) ρ =0.1 ; (b) ρ =0.2 ; (c) ρ =0.3 ; (d)ρ=0.4Fig.12.Relationship between different visibilities and evacuation time under different density:(a) ρ =0.1 ; (b) ρ =0.2 ;(c) ρ =0.3 ; (d) ρ =0.4 .

    從圖12(a)可以看出, 房間內(nèi)行人密度為0.1,四種跟隨策略疏散的時(shí)間隨著視野半徑的增大而減少.在房間內(nèi)視野半徑小于7 m時(shí), 四種疏散時(shí)間大小關(guān)系為TS2>TS1>TS4>TS3.表明策略S2 (從眾跟隨)的疏散時(shí)間最大, 而策略 S1 (跟隨周圍行人)跟隨時(shí)間較大, 策略S4 (跟隨墻壁移動(dòng))疏散時(shí)間與策略S3 (跟隨較多行人方向)的時(shí)間較小, 這是因?yàn)椴扇4策略靠墻行走的行人在墻壁、拐角和出口處的行人存在沖突.而在視野半徑大于7 m時(shí), 四種策略的疏散時(shí)間相差不大, 這是由于此時(shí)行人都具有清晰的視野, 不論采取何種跟隨策略的行人都明確地知道疏散出口的位置, 從而采取相同的移動(dòng)規(guī)則直接向出口方向移動(dòng).

    圖12(b)顯示到, 房間內(nèi)行人密度為0.2, 四種跟隨策略的疏散時(shí)間隨著視野半徑的增大先迅速減少再緩慢地減少.在視野半徑小于7 m時(shí), 四種疏散時(shí)間大小關(guān)系為TS2>TS1>TS3≥TS4.而圖12(c)和圖12(d)變化趨勢(shì)和圖12(b)相似, 不同的是, 在視野半徑小于4 m時(shí), 四種疏散時(shí)間大小關(guān)系為TS2>TS1>TS3≥TS4, 在 4—7 m 時(shí),S1, S3和S4策略疏散時(shí)間基本一致.而且三幅圖相同的是, 在視野半徑大于7 m后, 四種策略疏散時(shí)間相差不大.

    因此, 對(duì)于跟隨策略S1和S2來(lái)說(shuō), 視野半徑越大則其疏散時(shí)間越小, 而且在R= 2, 3, 4 m 處疏散時(shí)間有顯著變化, 而在視野半徑大于5 m后,疏散時(shí)間緩慢減少而后趨于穩(wěn)定.對(duì)于跟隨策略S3來(lái)說(shuō), 視野半徑對(duì)其疏散時(shí)間影響較大, 密度大于0.1, 在視野半徑為1和2 m時(shí), 受視野的影響,疏散時(shí)間TS3較大.對(duì)于策略S4來(lái)說(shuō), 由于采取靠墻行走策略, 即使在視野較差時(shí), 也能找到出口.總體來(lái)說(shuō), 在密度為0.2—0.4時(shí), 各策略基本保持相同的變化趨勢(shì), 而S1, S2和S3受視野半徑影響較大, 在視野半徑小于 7 m 時(shí),TS2>TS1>TS3≥TS4, 而視野半徑大于 7 m 后,TS2=TS1=TS4=TS3, 各策略疏散時(shí)間無(wú)明顯差異.

    3.3 熟悉環(huán)境者密度對(duì)四類跟隨策略疏散時(shí)間的影響

    由 3.1節(jié)可知, 疏散時(shí)間在視野半徑R=3 m后保持緩慢地減少, 因此本節(jié)研究在房間內(nèi)行人密度為 0.1—0.4, 視野半徑為 3 m 時(shí), 四類疏散策略下的疏散時(shí)間TS1,TS2,TS3和TS4隨著熟悉環(huán)境者密度的變化情況, 來(lái)比較熟悉環(huán)境者密度對(duì)四種跟隨策略的影響.當(dāng)熟悉環(huán)境者密度較大時(shí),不熟悉環(huán)境者密度較小, 導(dǎo)致選擇跟隨策略人數(shù)較少, 為比較不熟悉環(huán)境者的四類跟隨策略的優(yōu)劣,本文考慮熟悉環(huán)境者比例為0—0.6.

    從圖13(a)可以看到, 在密度為0.1, 視野半徑為 3 m, 熟悉環(huán)境者密度小于 0.2時(shí), 四類跟隨疏散時(shí)間大小為TS2>TS1>TS3>TS4.而在熟悉環(huán)境者密度大于 0.2 時(shí),TS2>TS1>TS4>TS3.這是由于選擇策略S4的行人靠墻進(jìn)行運(yùn)動(dòng), 不受熟悉環(huán)境者密度大小的影響, 因此疏散時(shí)間保持平緩地減少.而對(duì)于選擇策略S3的行人來(lái)說(shuō), 熟悉環(huán)境者密度越大, 房間內(nèi)知道出口方向的行人越多, 因此策略 S3 就越有效.同樣, 圖13(b)和圖13(a)的變化趨勢(shì)一致.

    從圖13(c)和圖13(d)可以看出, 疏散時(shí)間隨熟悉環(huán)境者密度的變化曲線相似, 在熟悉環(huán)境者密度小于 0.5 時(shí),TS2>TS1>TS3>TS4.而在熟悉環(huán)境者密度大于0.5時(shí), 四種策略疏散時(shí)間相差不大, 這是因?yàn)榉块g內(nèi)熟悉環(huán)境者占了大多數(shù), 大多數(shù)人能準(zhǔn)確地知道出口方向, 對(duì)于 S1, S2, S3 這三種策略, 不論行人采取何種跟隨策略都能準(zhǔn)確地進(jìn)行跟隨, 從而找到正確的疏散出口, 因此這三種跟隨策略疏散時(shí)間相差不大.

    因此, 熟悉環(huán)境者密度的大小影響著S1, S2,S3這三種跟隨策略的疏散, 在熟悉環(huán)境者密度小于 0.5 時(shí),TS2>TS1>TS3, 這是由于熟悉環(huán)境者密度較小時(shí), 房間內(nèi)知道出口的行人較少, 而采取跟隨策略S1和S2的行人都無(wú)法進(jìn)行正確的跟隨,但隨著熟悉環(huán)境者密度的增加, 房間內(nèi)知道準(zhǔn)確出口的行人變多, 策略S1和S2也能準(zhǔn)確地進(jìn)行跟隨.而在熟悉環(huán)境者密度大于0.5時(shí), 三種策略都能取得較好的效果, 疏散時(shí)間相差不大, 因?yàn)椴扇∪N跟隨策略都能找到準(zhǔn)確的出口.對(duì)于疏散策略S4來(lái)說(shuō), 采取靠墻行走的方式, 幾乎不受熟悉環(huán)境者密度的影響.

    3.4 單一策略條件下視野半徑對(duì)疏散時(shí)間的影響

    圖13 不同密度下各策略環(huán)境熟悉度與疏散時(shí)間的關(guān)系 (a) ρ =0.1 ; (b) ρ =0.2 ; (c) ρ =0.3 ; (d)ρ=0.4Fig.13.Relationship between environmental familiarity and evacuation time under different density:(a) ρ =0.1 ; (b) ρ =0.2 ;(c) ρ =0.3 ; (d) ρ =0.4 .

    以上部分討論了在同一房間, 采取多種疏散策略的疏散情況.為了進(jìn)一步比較四類跟隨策略的優(yōu)劣, 本節(jié)考慮房間內(nèi)采取單一跟隨策略的疏散情況, 研究視野半徑的變化對(duì)各類跟隨策略疏散時(shí)間的影響, 從而進(jìn)一步比較四類跟隨策略的優(yōu)劣.

    從圖14(a)和圖14(b)可以看到, 策略S1和S2隨著視野半徑的增加, 疏散時(shí)間緩慢地減少, 直到視野半徑為8 m時(shí), 疏散時(shí)間基本保持不變.可以理解的是, 對(duì)于策略S1和S2, 當(dāng)視野半徑較小時(shí), 采取策略S1和S2的行人沒有具體的跟隨方向, 而隨著視野半徑的增加, 在出口附近的行人能看到疏散出口位置, 從而直接移動(dòng)到出口.

    從圖14(c)可以看到, 對(duì)于策略S3, 隨著視野半徑的增加, 疏散時(shí)間首先迅速減少, 然后緩慢減少, 最后疏散時(shí)間保持不變.而在視野半徑小于3 m時(shí), 疏散時(shí)間較大, 隨著視野半徑的增大疏散時(shí)間迅速減少, 這是由于策略S3對(duì)視野半徑比較敏感, 在可見出口區(qū)域有行人往出口移動(dòng), 會(huì)引導(dǎo)其他行人往出口方向移動(dòng), 從而疏散時(shí)間迅速減少.而在視野半徑R≥3 m 后, 疏散時(shí)間曲線保持緩慢減少.

    圖14(d)和圖14(c)的疏散時(shí)間變化趨勢(shì)類似, 不同的是選擇策略S4的行人在視野半徑較小時(shí), 也能沿著墻壁行走且順利找到出口, 所以疏散時(shí)間較S3來(lái)說(shuō)較小.同樣, 隨著房間內(nèi)視野條件變好, 可視出口的范圍越大, 在可視出口區(qū)域的行人可以直接向出口進(jìn)行移動(dòng), 從而導(dǎo)致疏散時(shí)間較少.

    圖15顯示了在視野半徑下四種疏散策略疏散時(shí)間的變化.可以發(fā)現(xiàn), 在不同密度下, 各策略疏散時(shí)間變化趨勢(shì)一致, 而且各策略的疏散時(shí)間關(guān)系與視野半徑的大小密切相關(guān), 如下所示:

    結(jié)果表明, 各類跟隨策略的疏散時(shí)間與視野半徑的大小有關(guān).在視野半徑較小時(shí), 如R= 1 m,跟隨策略S4效果最優(yōu); 而視野半徑較大時(shí), 視野半徑大于2 m小于5 m時(shí), 跟隨策略S3才是更好的疏散策略; 視野半徑更大時(shí), 如視野半徑大于5 m小于 8 m時(shí), 跟隨策略 S3和 S4的疏散時(shí)間相差不大, 效果都較好; 而視野半徑大于8 m時(shí),四種跟隨疏散策略沒有優(yōu)劣, 疏散時(shí)間基本一致.

    圖14 單一策略下不同密度的視野半徑與疏散時(shí)間的關(guān)系 (a) 策略 S1; (b)策略 S2; (c)策略 S3; (d)策略 S4Fig.14.Relationship between different visibilities and evacuation time at different densities under a single strategy:(a) Strategy S1;(b) strategy S2; (c) strategy S3; (d) strategy S4.

    圖15 不同密度下各策略的視野半徑與疏散時(shí)間的關(guān)系 (a) ρ =0.1 ; (b) ρ =0.2 ; (c) ρ =0.3 ; (d)ρ=0.4Fig.15.Relationship of different strategy between different visibilities and evacuation time at different densities:(a) ρ =0.1 ;(b) ρ =0.2 ; (c) ρ =0.3 ; (d) ρ =0.4 .

    4 結(jié)論和討論

    本文研究了疏散行人的跟隨行為特性, 建立了元胞自動(dòng)機(jī)模型來(lái)仿真能見度受限和熟悉環(huán)境者人數(shù)變化下的行人跟隨行為表現(xiàn), 改進(jìn)勢(shì)函數(shù)元胞自動(dòng)機(jī)模型來(lái)模擬熟悉環(huán)境人群的疏散行為, 對(duì)不熟悉環(huán)境人群制定跟隨策略來(lái)模擬其跟隨行為, 觀察模擬仿真情況后, 得出以下結(jié)論.

    1)房間內(nèi)視野條件越好或熟悉環(huán)境人群數(shù)量越多, 則行人疏散時(shí)間越少.而且當(dāng)熟悉環(huán)境者密度Ks= 0.4 或視野半徑R= 3 m 時(shí), 疏散時(shí)間開始變得較小, 意味著行人能開始得到快速的疏散.

    2)跟隨策略的有效性和視野半徑的大小有關(guān),在熟悉環(huán)境者密度一定, 視野半徑很小時(shí), 跟隨策略S4效果最好; 視野半徑較小時(shí), 跟隨策略S4和S3效果較好, 跟隨策略S1和S2效果較差, 此時(shí)疏散時(shí)間大小關(guān)系為TS2>TS1>TS3≥TS4.而在視野半徑繼續(xù)增大時(shí), 四種疏散策略時(shí)間沒有較大差異.

    3)當(dāng)視野半徑一定且熟悉環(huán)境者密度小于0.5 時(shí), 跟隨策略 S3 效果最好, 而跟隨策略 S1, S2,S4 效果較差, 此時(shí)疏散時(shí)間大小關(guān)系為TS2>TS1>TS4>TS3.在熟悉環(huán)境者密度大于 0.5 時(shí), 三種策略都能取得較好的效果, 疏散時(shí)間相差不大.而且策略S4受熟悉環(huán)境者密度影響不大.

    4)在單一策略下, 視野半徑對(duì)跟隨策略影響很大.在視野半徑較小時(shí), 疏散策略S4效果最優(yōu),此時(shí)跟隨墻壁行走為更好的策略; 在視野半徑較大時(shí), 跟隨策略S3效果最好, 這和Wang和Cao[23]的研究表現(xiàn)一致, 由于策略S4的行人在墻壁、拐角和出口處存在較多沖突, 而策略S3對(duì)視野半徑比較敏感, 在可見出口區(qū)域有行人往出口移動(dòng), 會(huì)引導(dǎo)S3策略行人往出口方向移動(dòng), 從而疏散時(shí)間較策略S4更小; 而當(dāng)房間內(nèi)視野半徑更好時(shí), 四種策略疏散時(shí)間沒有較大差異.

    本文的研究結(jié)果和其他研究者的結(jié)果有一定的區(qū)別和相似之處.Hou 等[18]、Ma 等[19,20]建立社會(huì)力模型來(lái)研究有限能見度下行人的領(lǐng)導(dǎo)作用, 結(jié)果表明有限能見度下領(lǐng)導(dǎo)作用對(duì)疏散影響是雙重的.Ma等[20]認(rèn)為在視野半徑較大時(shí), 未看到出口的行人都會(huì)向領(lǐng)導(dǎo)者靠攏, 而不同領(lǐng)導(dǎo)者所處房間位置不同, 而增加了疏散時(shí)間.而本文的研究結(jié)果為, 視野條件越好行人疏散時(shí)間越少, 因?yàn)橐曇皸l件越好, 知道出口的行人越多, 則跟隨效果越好,這與本文所考慮的跟隨行為有關(guān).Frank和Dorso[22]以及Wang和Cao[23]考慮了在視野受限的房間內(nèi),行人疏散過(guò)程的放牧、沿墻壁行走等行為, 他們的實(shí)驗(yàn)表明疏散時(shí)間隨著視野半徑的增大而變小, 這和本文有著相同的結(jié)論.Wang和Cao[23]也在考慮行人視野受限時(shí), 建立了社會(huì)力模型研究行人的疏散策略包括跟隨平均位置、平均移動(dòng)方向和沿墻行走等行為, 他們認(rèn)為單一策略下沿墻行走和跟隨平均移動(dòng)方向等策略更有效.這和本文研究結(jié)果在單一策略環(huán)境下有相似之處, 而與多種策略并存情況下有些區(qū)別, 由于Wang和Cao[23]只考慮了單一策略下疏散行為, 而本文考慮了熟悉環(huán)境人群疏散和多種跟隨策略下行人的疏散, 更綜合地分析了行人的疏散行為, 認(rèn)為在多種跟隨策略下, 能見度較差時(shí), 跟隨策略 S4 (沿墻行走)更有效, 而能見度較好時(shí), 跟隨策略S3 (跟隨最多移動(dòng)方向)更有效.

    本文研究的內(nèi)容更切合實(shí)際地分析了房間內(nèi)的疏散行為, 研究了不熟悉環(huán)境人群在能見度受限下采取不同跟隨策略的表現(xiàn), 此項(xiàng)工作能對(duì)大型公共場(chǎng)所如超市、體育館的應(yīng)急疏散情況和視野條件較差時(shí)提供一些啟示, 有助于制定一些有效的指導(dǎo)策略.但是本文研究的模擬環(huán)境較為基礎(chǔ), 在未來(lái)的研究過(guò)程中需要研究出口寬度、多出口條件下以及障礙物情況下的疏散情況, 以研究這些跟隨策略在未來(lái)的適用性.而且, 目前本文的研究?jī)H處于理論層面, 暫未進(jìn)行真人實(shí)驗(yàn)來(lái)標(biāo)定和驗(yàn)證, 在未來(lái)的研究中, 可結(jié)合實(shí)際條件對(duì)現(xiàn)實(shí)場(chǎng)景進(jìn)行真人實(shí)驗(yàn)以進(jìn)一步研究跟隨策略的有效性和適用性.

    猜你喜歡
    元胞行人視野
    毒舌出沒,行人避讓
    意林(2021年5期)2021-04-18 12:21:17
    居· 視野
    中華民居(2020年3期)2020-07-24 01:48:04
    路不為尋找者而設(shè)
    基于元胞自動(dòng)機(jī)下的交通事故路段仿真
    我是行人
    基于元胞數(shù)據(jù)的多維數(shù)據(jù)傳遞機(jī)制
    視野
    科學(xué)家(2015年2期)2015-04-09 02:46:46
    真相
    讀者(2014年18期)2014-05-14 11:40:56
    基于AIS的航道移動(dòng)瓶頸元胞自動(dòng)機(jī)模型
    零邊界條件下二維元胞自動(dòng)機(jī)矩陣可逆性分析
    最近最新中文字幕大全电影3 | 久久久久网色| 久久人人97超碰香蕉20202| 免费在线观看完整版高清| 久久精品国产亚洲av高清一级| 久久这里只有精品19| 脱女人内裤的视频| 亚洲熟女毛片儿| 欧美在线黄色| 欧美在线一区亚洲| 熟女少妇亚洲综合色aaa.| 香蕉久久夜色| 免费看a级黄色片| 亚洲,欧美精品.| 啪啪无遮挡十八禁网站| 欧美成人免费av一区二区三区 | 黄色视频在线播放观看不卡| 国产不卡一卡二| 中文字幕另类日韩欧美亚洲嫩草| 欧美黄色淫秽网站| 色播在线永久视频| 亚洲成人免费av在线播放| 久久精品人人爽人人爽视色| 久久这里只有精品19| 黄色毛片三级朝国网站| av天堂在线播放| 十八禁网站免费在线| 免费人妻精品一区二区三区视频| 久久人人爽av亚洲精品天堂| 亚洲精品中文字幕在线视频| 99在线人妻在线中文字幕 | 人人妻人人添人人爽欧美一区卜| 女人被躁到高潮嗷嗷叫费观| 国产免费现黄频在线看| 性少妇av在线| 中文字幕精品免费在线观看视频| 亚洲国产毛片av蜜桃av| 国产有黄有色有爽视频| 免费不卡黄色视频| 丰满少妇做爰视频| 搡老岳熟女国产| 俄罗斯特黄特色一大片| 国产成人欧美在线观看 | 天天躁狠狠躁夜夜躁狠狠躁| 久久99一区二区三区| 精品国产一区二区久久| 一进一出抽搐动态| 久久中文看片网| 在线天堂中文资源库| 国产精品久久久久久人妻精品电影 | 大香蕉久久成人网| 十分钟在线观看高清视频www| 欧美中文综合在线视频| 久久国产精品人妻蜜桃| 国产aⅴ精品一区二区三区波| 日日夜夜操网爽| 在线观看免费高清a一片| 国产亚洲精品一区二区www | 午夜两性在线视频| 女人被躁到高潮嗷嗷叫费观| 丝袜喷水一区| 午夜福利一区二区在线看| 久久av网站| 深夜精品福利| 成年人黄色毛片网站| 精品第一国产精品| 纯流量卡能插随身wifi吗| 黄色视频不卡| 日韩大片免费观看网站| 精品福利观看| 18禁国产床啪视频网站| 中文字幕人妻熟女乱码| 欧美日韩国产mv在线观看视频| 乱人伦中国视频| 一个人免费在线观看的高清视频| 精品午夜福利视频在线观看一区 | 老司机影院毛片| 老司机午夜十八禁免费视频| 色老头精品视频在线观看| h视频一区二区三区| 久久久久久久久久久久大奶| 欧美国产精品一级二级三级| 俄罗斯特黄特色一大片| 99久久精品国产亚洲精品| 精品久久久精品久久久| 国产日韩欧美在线精品| 色精品久久人妻99蜜桃| 精品少妇久久久久久888优播| 丝袜在线中文字幕| 90打野战视频偷拍视频| 亚洲综合色网址| 成年版毛片免费区| 欧美黑人欧美精品刺激| 午夜老司机福利片| 欧美性长视频在线观看| 丝袜喷水一区| 精品久久久久久久毛片微露脸| 亚洲熟女精品中文字幕| 久久ye,这里只有精品| 久久精品亚洲熟妇少妇任你| 亚洲精品中文字幕在线视频| 欧美变态另类bdsm刘玥| 90打野战视频偷拍视频| 亚洲精品自拍成人| 成人18禁高潮啪啪吃奶动态图| 久久久久久免费高清国产稀缺| 久久久久精品人妻al黑| 日韩一卡2卡3卡4卡2021年| 久久精品亚洲精品国产色婷小说| 一区福利在线观看| 男人操女人黄网站| 亚洲国产欧美网| 亚洲成人免费av在线播放| 国产成人精品无人区| 丰满人妻熟妇乱又伦精品不卡| 午夜两性在线视频| 99精国产麻豆久久婷婷| 亚洲欧洲精品一区二区精品久久久| 建设人人有责人人尽责人人享有的| 国产精品二区激情视频| 操出白浆在线播放| 五月天丁香电影| 国产精品久久久人人做人人爽| 久久精品亚洲av国产电影网| 99国产精品99久久久久| 99re在线观看精品视频| 欧美日韩一级在线毛片| 亚洲自偷自拍图片 自拍| 在线观看免费视频网站a站| 99riav亚洲国产免费| 色视频在线一区二区三区| 午夜精品国产一区二区电影| 国产亚洲精品久久久久5区| 欧美中文综合在线视频| 久久人人97超碰香蕉20202| 少妇裸体淫交视频免费看高清 | 亚洲精品中文字幕一二三四区 | 午夜福利免费观看在线| 亚洲avbb在线观看| 999久久久国产精品视频| 亚洲九九香蕉| 久久国产精品人妻蜜桃| 99久久人妻综合| 一夜夜www| 黑人猛操日本美女一级片| 久久婷婷成人综合色麻豆| 国产男女超爽视频在线观看| 老司机午夜福利在线观看视频 | 高清黄色对白视频在线免费看| www.熟女人妻精品国产| 天天躁夜夜躁狠狠躁躁| 国产精品 国内视频| 天堂8中文在线网| 黑人巨大精品欧美一区二区蜜桃| 午夜福利在线观看吧| 91老司机精品| 韩国精品一区二区三区| 这个男人来自地球电影免费观看| 免费一级毛片在线播放高清视频 | 香蕉丝袜av| 亚洲人成电影免费在线| 伦理电影免费视频| 纯流量卡能插随身wifi吗| 肉色欧美久久久久久久蜜桃| 黄色怎么调成土黄色| 亚洲伊人色综图| 精品福利永久在线观看| 在线观看免费高清a一片| 巨乳人妻的诱惑在线观看| 十八禁人妻一区二区| 国产1区2区3区精品| 国产一区有黄有色的免费视频| 麻豆av在线久日| 精品国产乱子伦一区二区三区| 人人妻,人人澡人人爽秒播| 电影成人av| 看免费av毛片| 搡老岳熟女国产| 亚洲五月婷婷丁香| 国产在视频线精品| 热99久久久久精品小说推荐| 国产成人欧美在线观看 | 高清视频免费观看一区二区| 婷婷丁香在线五月| 91字幕亚洲| 最近最新中文字幕大全免费视频| 日日爽夜夜爽网站| 久久毛片免费看一区二区三区| 欧美精品一区二区免费开放| 亚洲情色 制服丝袜| 欧美亚洲日本最大视频资源| 深夜精品福利| 久久午夜亚洲精品久久| 国产精品一区二区在线不卡| 在线观看一区二区三区激情| 99久久国产精品久久久| 亚洲欧美激情在线| 亚洲七黄色美女视频| 国产一区二区 视频在线| 国产深夜福利视频在线观看| 亚洲五月婷婷丁香| 岛国在线观看网站| 成人免费观看视频高清| 男女高潮啪啪啪动态图| 久久久国产欧美日韩av| 欧美成狂野欧美在线观看| 欧美乱码精品一区二区三区| 亚洲精品美女久久av网站| 精品欧美一区二区三区在线| www.自偷自拍.com| 国产精品成人在线| 热99国产精品久久久久久7| 亚洲av美国av| 十分钟在线观看高清视频www| 亚洲中文字幕日韩| 国内毛片毛片毛片毛片毛片| 久久 成人 亚洲| 亚洲国产看品久久| 菩萨蛮人人尽说江南好唐韦庄| 9热在线视频观看99| 久久久水蜜桃国产精品网| 婷婷丁香在线五月| 国产男靠女视频免费网站| 亚洲一码二码三码区别大吗| 亚洲av成人一区二区三| av欧美777| 国产黄色免费在线视频| 午夜免费鲁丝| 欧美亚洲日本最大视频资源| 亚洲av日韩在线播放| 青青草视频在线视频观看| 欧美精品高潮呻吟av久久| 亚洲熟女精品中文字幕| 水蜜桃什么品种好| 欧美激情高清一区二区三区| 国产成人免费无遮挡视频| 人妻 亚洲 视频| 极品人妻少妇av视频| 麻豆成人av在线观看| 91成年电影在线观看| av又黄又爽大尺度在线免费看| 在线 av 中文字幕| 亚洲成人免费av在线播放| 久久中文字幕一级| 黄色片一级片一级黄色片| 国产av国产精品国产| 国产精品久久久久久精品古装| 久久精品国产亚洲av香蕉五月 | 亚洲欧美一区二区三区黑人| 国产区一区二久久| 色播在线永久视频| 午夜精品国产一区二区电影| 人人妻人人爽人人添夜夜欢视频| 99国产精品一区二区蜜桃av | 少妇 在线观看| 国产av精品麻豆| 久久人妻av系列| 亚洲国产毛片av蜜桃av| 日本撒尿小便嘘嘘汇集6| 久久久精品区二区三区| 久久久久国内视频| 美女国产高潮福利片在线看| 精品免费久久久久久久清纯 | 老熟女久久久| 汤姆久久久久久久影院中文字幕| 国产精品99久久99久久久不卡| 午夜福利在线免费观看网站| 久久国产亚洲av麻豆专区| 18禁国产床啪视频网站| 精品午夜福利视频在线观看一区 | 亚洲色图 男人天堂 中文字幕| 欧美日韩av久久| 俄罗斯特黄特色一大片| 亚洲国产欧美网| kizo精华| 99热国产这里只有精品6| 一本大道久久a久久精品| 亚洲美女黄片视频| 久久久久久久精品吃奶| 狂野欧美激情性xxxx| 一本大道久久a久久精品| 午夜福利视频在线观看免费| 叶爱在线成人免费视频播放| 99久久精品国产亚洲精品| 9191精品国产免费久久| 国产高清视频在线播放一区| 午夜福利免费观看在线| 性少妇av在线| av电影中文网址| 丁香欧美五月| 最近最新中文字幕大全免费视频| 老司机福利观看| av免费在线观看网站| 亚洲av美国av| 国产区一区二久久| 国产精品一区二区精品视频观看| 国产精品久久久人人做人人爽| a级毛片在线看网站| 狂野欧美激情性xxxx| av免费在线观看网站| 他把我摸到了高潮在线观看 | 久久精品亚洲熟妇少妇任你| 久久狼人影院| 满18在线观看网站| 色尼玛亚洲综合影院| 久热爱精品视频在线9| 久久av网站| 汤姆久久久久久久影院中文字幕| 人人澡人人妻人| 天天躁夜夜躁狠狠躁躁| 久久国产精品大桥未久av| www.精华液| 午夜视频精品福利| 午夜精品久久久久久毛片777| 在线观看免费视频日本深夜| 又大又爽又粗| 亚洲欧美激情在线| a在线观看视频网站| 水蜜桃什么品种好| av超薄肉色丝袜交足视频| 亚洲va日本ⅴa欧美va伊人久久| 亚洲人成77777在线视频| 精品国内亚洲2022精品成人 | 久久国产精品大桥未久av| 国产熟女午夜一区二区三区| 免费少妇av软件| 亚洲五月婷婷丁香| 少妇猛男粗大的猛烈进出视频| 日韩成人在线观看一区二区三区| 国产av又大| 在线十欧美十亚洲十日本专区| 国产精品二区激情视频| 菩萨蛮人人尽说江南好唐韦庄| 在线十欧美十亚洲十日本专区| 日韩欧美三级三区| 亚洲天堂av无毛| 老司机深夜福利视频在线观看| 午夜福利乱码中文字幕| tocl精华| 99九九在线精品视频| 嫩草影视91久久| 亚洲国产精品一区二区三区在线| svipshipincom国产片| 老司机在亚洲福利影院| 欧美黄色片欧美黄色片| 色婷婷av一区二区三区视频| 精品福利观看| 捣出白浆h1v1| 国产97色在线日韩免费| 精品少妇一区二区三区视频日本电影| 男女下面插进去视频免费观看| 中文字幕精品免费在线观看视频| 亚洲七黄色美女视频| 欧美人与性动交α欧美软件| 18禁美女被吸乳视频| 色播在线永久视频| 国产有黄有色有爽视频| 热99国产精品久久久久久7| 黄频高清免费视频| 国产精品熟女久久久久浪| 免费女性裸体啪啪无遮挡网站| 久久精品亚洲精品国产色婷小说| 久久精品国产亚洲av高清一级| 亚洲自偷自拍图片 自拍| 肉色欧美久久久久久久蜜桃| 亚洲一区中文字幕在线| 亚洲国产欧美一区二区综合| 日韩一区二区三区影片| 正在播放国产对白刺激| 欧美激情 高清一区二区三区| www.熟女人妻精品国产| 黄片小视频在线播放| 伦理电影免费视频| 午夜福利,免费看| 亚洲情色 制服丝袜| 亚洲男人天堂网一区| 制服人妻中文乱码| 亚洲久久久国产精品| 我的亚洲天堂| 性高湖久久久久久久久免费观看| 欧美日韩国产mv在线观看视频| 极品教师在线免费播放| 免费日韩欧美在线观看| 欧美精品一区二区大全| 黄色成人免费大全| 一区二区三区精品91| 黄网站色视频无遮挡免费观看| 久久性视频一级片| 久久久久久久国产电影| 久9热在线精品视频| 欧美日韩黄片免| 久久香蕉激情| www.自偷自拍.com| 日本vs欧美在线观看视频| 在线观看人妻少妇| av网站在线播放免费| 久久亚洲精品不卡| 最近最新中文字幕大全电影3 | 国产精品麻豆人妻色哟哟久久| 夜夜夜夜夜久久久久| 新久久久久国产一级毛片| 一本色道久久久久久精品综合| 天天影视国产精品| 黄色丝袜av网址大全| 人妻一区二区av| 黄片播放在线免费| 亚洲中文av在线| 午夜91福利影院| 免费观看av网站的网址| 97人妻天天添夜夜摸| 99国产精品一区二区蜜桃av | 一进一出好大好爽视频| 国产日韩一区二区三区精品不卡| 国产欧美亚洲国产| 咕卡用的链子| 人成视频在线观看免费观看| 亚洲中文字幕日韩| 国产欧美日韩综合在线一区二区| 十八禁高潮呻吟视频| 国产单亲对白刺激| 少妇被粗大的猛进出69影院| 人人澡人人妻人| 国产在线免费精品| 一本大道久久a久久精品| 99国产精品一区二区蜜桃av | 久久av网站| 99精品在免费线老司机午夜| 久久中文字幕人妻熟女| 下体分泌物呈黄色| 久久九九热精品免费| 日韩免费av在线播放| 色婷婷av一区二区三区视频| 成年动漫av网址| 国产aⅴ精品一区二区三区波| 捣出白浆h1v1| 一区二区三区乱码不卡18| 高清黄色对白视频在线免费看| 国产精品1区2区在线观看. | 色婷婷av一区二区三区视频| 国产精品香港三级国产av潘金莲| 国产高清国产精品国产三级| 19禁男女啪啪无遮挡网站| 欧美日韩国产mv在线观看视频| 亚洲avbb在线观看| 国产91精品成人一区二区三区 | 肉色欧美久久久久久久蜜桃| 国产精品免费视频内射| 69精品国产乱码久久久| 精品亚洲成国产av| 黄色丝袜av网址大全| 精品人妻熟女毛片av久久网站| 色播在线永久视频| 交换朋友夫妻互换小说| 精品一区二区三区视频在线观看免费 | 亚洲专区中文字幕在线| 老司机午夜福利在线观看视频 | 亚洲国产欧美在线一区| kizo精华| 国产成人精品久久二区二区免费| 69av精品久久久久久 | 极品人妻少妇av视频| 啦啦啦免费观看视频1| 国产精品一区二区在线不卡| 久久久国产精品麻豆| 国产精品熟女久久久久浪| 久久久久国内视频| 大型黄色视频在线免费观看| 精品午夜福利视频在线观看一区 | 国产男女内射视频| 亚洲精品美女久久久久99蜜臀| 黑人操中国人逼视频| 亚洲黑人精品在线| 飞空精品影院首页| 老司机午夜福利在线观看视频 | 欧美日韩亚洲国产一区二区在线观看 | 在线观看一区二区三区激情| 欧美大码av| 国产精品九九99| 美国免费a级毛片| 日韩人妻精品一区2区三区| 国产成人欧美| 18禁裸乳无遮挡动漫免费视频| 国产主播在线观看一区二区| 久久99热这里只频精品6学生| 免费高清在线观看日韩| 黄色成人免费大全| 国产极品粉嫩免费观看在线| 国产淫语在线视频| 91老司机精品| 日韩欧美免费精品| 成人特级黄色片久久久久久久 | 中文字幕高清在线视频| 露出奶头的视频| 天天躁日日躁夜夜躁夜夜| 美女高潮喷水抽搐中文字幕| 国产亚洲精品久久久久5区| 亚洲精品av麻豆狂野| 99国产精品免费福利视频| √禁漫天堂资源中文www| 国产淫语在线视频| av网站在线播放免费| 国产xxxxx性猛交| 两个人看的免费小视频| 免费高清在线观看日韩| 久久国产精品影院| 人人妻人人澡人人看| 成年动漫av网址| 啦啦啦中文免费视频观看日本| 免费女性裸体啪啪无遮挡网站| 亚洲三区欧美一区| 国产97色在线日韩免费| 免费在线观看完整版高清| 免费一级毛片在线播放高清视频 | 一区二区av电影网| 巨乳人妻的诱惑在线观看| 男女高潮啪啪啪动态图| 亚洲精品美女久久av网站| 欧美日韩国产mv在线观看视频| cao死你这个sao货| 精品一区二区三区视频在线观看免费 | 亚洲欧美激情在线| 亚洲成a人片在线一区二区| 久久天堂一区二区三区四区| 精品卡一卡二卡四卡免费| 日韩人妻精品一区2区三区| 亚洲专区中文字幕在线| 多毛熟女@视频| 久久久国产精品麻豆| 午夜激情av网站| 最黄视频免费看| 久久精品熟女亚洲av麻豆精品| 国产av国产精品国产| 亚洲欧美一区二区三区黑人| 亚洲三区欧美一区| 亚洲国产av影院在线观看| 久久狼人影院| 亚洲欧美日韩高清在线视频 | 在线 av 中文字幕| 国产色视频综合| av天堂在线播放| 国产成人啪精品午夜网站| 蜜桃在线观看..| 王馨瑶露胸无遮挡在线观看| 交换朋友夫妻互换小说| 91字幕亚洲| 高清在线国产一区| 大陆偷拍与自拍| 亚洲三区欧美一区| 中文字幕色久视频| 99在线人妻在线中文字幕 | 久久精品国产亚洲av香蕉五月 | 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲国产毛片av蜜桃av| 免费在线观看日本一区| 国产有黄有色有爽视频| 亚洲精品成人av观看孕妇| 99精品久久久久人妻精品| 丁香欧美五月| 淫妇啪啪啪对白视频| 一个人免费在线观看的高清视频| 精品福利观看| 菩萨蛮人人尽说江南好唐韦庄| 国产91精品成人一区二区三区 | 极品教师在线免费播放| 少妇 在线观看| 人人妻人人澡人人看| 久久精品熟女亚洲av麻豆精品| 精品国产乱子伦一区二区三区| 中文字幕人妻丝袜制服| 亚洲五月色婷婷综合| 热99re8久久精品国产| av欧美777| 国产精品电影一区二区三区 | 桃红色精品国产亚洲av| 精品亚洲成a人片在线观看| 两个人看的免费小视频| 最近最新中文字幕大全免费视频| 18禁观看日本| 欧美精品av麻豆av| 欧美激情高清一区二区三区| 少妇粗大呻吟视频| 国产精品欧美亚洲77777| 精品久久久久久久毛片微露脸| 成人特级黄色片久久久久久久 | 女人高潮潮喷娇喘18禁视频| 国产日韩一区二区三区精品不卡| 亚洲av第一区精品v没综合| 在线亚洲精品国产二区图片欧美| 国产成人啪精品午夜网站| 成在线人永久免费视频| 999久久久精品免费观看国产| 婷婷丁香在线五月| 久久精品人人爽人人爽视色| 黄网站色视频无遮挡免费观看| 桃花免费在线播放| 1024香蕉在线观看| 一区二区三区乱码不卡18| 欧美精品一区二区免费开放| 精品一品国产午夜福利视频| 99久久99久久久精品蜜桃| 欧美日韩成人在线一区二区| 日本欧美视频一区| 日日爽夜夜爽网站| 日韩大码丰满熟妇| 欧美在线黄色| 9191精品国产免费久久| 99riav亚洲国产免费| 新久久久久国产一级毛片| 久久影院123| 91av网站免费观看| 久久久久久免费高清国产稀缺| 可以免费在线观看a视频的电影网站| 91精品三级在线观看| 久久久久网色| 国产精品美女特级片免费视频播放器 |